Fast algorithms for problems on thermal tomography

Abstract

In this paper, we study an ill-posed, nonlinear inverse problem in heat conduction and hydrology applications. In [2], the problem is linearized to give a linear integral equation, which is then solved by the Tikhonov method with the identity as the regularization operator. We prove in this paper that the resulting equation is well-condition and has clustered spectrum. Hence if the conjugate gradient method is used to solve the equation, we expect superlinear convergence. However, we note that the identity operator does not give good solution to the original equation in general. Therefore in this paper, we use the Laplacian operator as the regularization operator instead. With the Laplacian operator, the regularized equation is ill-conditioned and hence a preconditioner is required to speed up the convergence rate if the equation is solved by the conjugate gradient method. We here propose to use the Laplacian operator itself as preconditioner. This preconditioner can be inverted easily by fast sine-transforms and we prove that the resulting preconditioned system is well-conditioned and has clustered spectrum too. Hence the conjugate gradient method converges superlinearly for the preconditioned system. Numerical results are given to illustrate the fast convergence.

Research supported by CUHK Direct Grant 220600680.

Research supported by CUHK Grants for Overseas Academic Activities.