Numerical implementations of dynamical probe method for non-stationary heat equation
Introduction
In the recent 15 years, several reconstruction schemes for identifying an unknown inclusion for steady state heat conductor have been proposed. They are the linear sampling method, probe method, singular sources method, no response test, range test, factorization method, enclosure method etc, see [2], [6], [7], [9], [11], [12], [13] and the references therein. None of these schemes is data fitting scheme using Tikhonov’s regularization. Each of them gave an exact scheme theoretically to recover the unknown inclusion directly from given measurement data.
On the other hand, for non steady state conductors, Nakamura and his collaborators have recently developed a probe typed reconstruction scheme for the one space dimensional case [3] for heat conduction process, which aims to detect the unknown inclusion of heat conductive material from the boundary measurement data. We refer to this scheme by dynamical probe method. Compared with the standard optimization method which considers the solution to this inverse problem as an minimizer of some cost functional, the advantage of this new scheme is that it is a theoretically exact method and no iteration process is required. The basic idea of this scheme is to construct an indicator function which can detect the boundary of the inclusion by using the boundary measurement data of heat field. For this function, the variable is the parameter . When approaches to the boundary of inclusion from outside the inclusion, this indicator blows up. The indicator function is defined via the pre-indicator function by taking some limiting procedure. For the numerical realization of this scheme, it is more convenient to use the pre-indicator function than the indicator function in order to avoid the limiting procedure.
Due to the diffusion of heat conduction with respect to time t, the information about the inclusion contained in the measurement data is relatively weak. Therefore, an efficient numerical realization scheme for this new dynamical probe method should be studied, which is the purpose of this paper.
On the other hand, it also should be noticed that the dynamical probe method is not the only one method to detect the inclusion in heat conductive media. A single one sided measurement version of the enclosure method has been applied in one dimensional case by Ikehata recently, see [8].
The inverse problem for identifying an unknown inclusion D inside a non steady state heat conductor in 1-dimensional case is formulated as follows. Let be an open connected interval in and be a fixed given time. For our forward problem, we consider the heat conduction problemwhere denotes the boundary of . We assume that the heat conductivity has the following formfor some constant k satisfying , where is some unknown open, connected interval and is the characteristic function of D. In other words, the function has jump discontinuity across , the boundary of D. For given , there exists a unique weak solutionto the forward problem (1.1) in the sensefor all with , see [14]. Here, and are the dual spaces of and , respectively. We also denote the solution u of (1.1) by . Based on this weak solution, we define the Neumann to Dirichlet mapwhich maps from to . The inverse problem considered in this paper is to identify from numerically.
This paper is organized as follows. The indicator function (i.e., a mathematical testing machine to identify the unknown inclusion) and pre-indicator function are defined in Section 2. We also give the behavior of the pre-indicator function in the 1-dimensional case using the behavior of the so called reflected solution. The key to define the pre-indicator function is the Runge approximation. In Sections 3 Realization of Runge approximation by regularization scheme, 4 Denseness of, the Runge approximation is given in a constructive way by using a regularizing method. That is in Section 3, we propose a regularizing scheme to determine the density function of the integral operator , and in Section 4, we prove the denseness of the range of operator which guarantees the existence of this density function. We also give in Section 4 the discrete argument for the numerical realization of the dynamical probe method. Finally in Section 5, some examples are presented to show the validity of the dynamical probe method.
Section snippets
Pre-indicator function and indicator function of dynamical probe method
Let with and be the closure of open domain . Also, let with a small be a domain
- (H1)
,
- (H2)
is Lipschitz,
- (H3)
is an empty set, or connected to for each , which means that each point in can be connected to a point in by a continuous curve in .
In the rest of this paper, we call such a
Realization of Runge approximation by regularization scheme
Let be intervals , with the two sets of ends and . We assume . In this case, for any given or and with , we can construct in a very simple way, that is, the configuration of is independent of , which will be denoted by in the sequel. Moreover, when we approximate by with being the right-hand side boundary of , the left-hand side keeps unchanged, see
Denseness of and discrete scheme
In this section, we prove the denseness of the range of map and the discrete scheme, that is the detailed discretization for computing and the simulation process for generating the boundary measurement data. As stated previously, we begin with ((3.7), (3.8), (3.9), (3.10), (3.11), (3.12), (3.13)). It is easy to see thatdefined by (3.6) has the adjoint operatorTherefore, the regularizing equation in (3.8) for given
Numerical results
From the simulation procedure for and the construction of given in Sections 3 Realization of Runge approximation by regularization scheme, 4 Denseness of, we can finally compute the indicator function and check the numerical performance of dynamical probe method in 1-dimensional case.
Based on the above scheme, we use the configurationsin the following Example 1, Example 2. Example 1 For a given inclusion , we
Acknowledgements
The authors would like to thank one of the referee for pointing out the mistakes in the original proof of Theorem 2.1 and providing the references [4], [5], which leads to an improved version of this paper. The third author is supported by Natural Science Foundation of China, No.10771033. The other authors are supported by Grant-in-Aid for Scientific Research (B) (No. 19340028) of Japan Society for Promotion of Science. The second and third authors thank Faculty of Science of Hokkaido
References (15)
Stability estimates and reconstructions in inverse scattering using singular sources
J. Comput. Appl. Math.
(2000)The One-Dimensional Heat Equation
(1984)- et al.
A simple method for solving inverse scattering problems in the resonance region
Inverse Probl.
(1996) - et al.
A probe method for the inverse boundary value problem of non-stationary heat equations
Inverse Probl.
(2007) A Runge theorem for solutions of the heat equation
Proc. Amer. Math. Soc.
(1980)An approximation theorem of Runge type for the heat equation
Proc. Amer. Math. Soc.
(1975)Reconstruction of the shape of the inclusion by boundary measurements
Comm. PDE
(1998)
Cited by (10)
The framework of the enclosure method with dynamical data and its applications
2011, Inverse ProblemsImproved asymptotic analysis for dynamical probe method
2016, Journal of Inverse and Ill-Posed ProblemsReconstruction algorithm for unknown cavities via Feynman–Kac type formula
2015, Computational Optimization and Applications