Modified Collocation Trefftz Method for the Geometry Boundary Identification Problem of Heat Conduction

Abstract

In this article, a meshless numerical algorithm is proposed for the boundary identification problem of heat conduction, one kind of inverse problem. In the geometry boundary identification problem, the Cauchy data is given for part of the boundary. The Neumann boundary condition is given for the other portion of the boundary, whose spatial position is unknown. In order to stably solve the inverse problem, the modified collocation Trefftz method, a promising boundary-type meshless method, is adopted for discretizing this problem. Since the spatial position for part of the boundary is unknown, the numerical discretization results in a system of nonlinear algebraic equations (NAEs). Then, the exponentially convergent scalar homotopy algorithm (ECSHA) is used to efficiently obtain the convergent solution of the system of NAEs. The ECSHA is insensitive to the initial guess of the evolutionary process. In addition, the efficiency of the computation is greatly improved, since calculation of the inverse of the Jacobian matrix can be avoided. Four numerical examples are provided to validate the proposed meshless scheme. In addition, some factors that might influence the performance of the proposed scheme are examined through a series of numerical experiments. The stability of the proposed scheme can be proven by adding some noise to the boundary conditions.

Notes

n _b2 = 20, m = 10, R ₀ = 3.0, icc = 0.2, icb = 1.2.

n _b1 = 30, n _b2 = 20, R ₀ = 3.0, icc = 0.2, icb = 1.2.

n _b1 = 30, n _b2 = 20, m = 10, R ₀ = 3.0, icc = 0.2, icb = 1.2.

n _b1 = 40, n _b2 = 20, m = 14, R ₀ = 3.0, icc = 2, icb = 1.4.

n _b1 = 50, n _b2 = 30, m = 20, R ₀ = 2.0, R ₁ = 0.2, icc = 0, icb = 0.6.