Table of Contents
Fetching ...

An Iterative Direct Sampling Method for Reconstructing Moving Inhomogeneities in Parabolic Problems

Bangti Jin, Fengru Wang, Jun Zou

TL;DR

This work introduces the Iterative Direct Sampling Method (IDSM) for reconstructing moving inhomogeneities in parabolic inverse problems from limited boundary data. By linking boundary measurements to the inclusion through a semi-linear forward map and constructing a probing kernel via a low-rank-augmented resolvent, IDSM directly images inclusions and handles nonlinear and mixed-type anomalies. The method operates on time segments with a fixed-point iteration to address nonlinearity, employs an $L^{2}$ inner product for stability, and updates a low-rank kernel on the fly, achieving robust reconstructions with only a handful of PDE solves per data set. Numerical experiments across merging/splitting, mixed-type, nonlinear, and fading/diminishing inclusions demonstrate accurate tracking under noise and limited data, underscoring the approach’s practical potential for time-dependent parabolic imaging.

Abstract

We propose in this work a novel iterative direct sampling method for imaging moving inhomogeneities in parabolic problems using boundary measurements. It can efficiently identify the locations and shapes of moving inhomogeneities when very limited data are available, even with only one pair of lateral Cauchy data, and enjoys remarkable numerical stability for noisy data and over an extended time horizon. The method is formulated in an abstract framework, and is applicable to linear and nonlinear parabolic problems, including linear, nonlinear, and mixed-type inhomogeneities. Numerical experiments across diverse scenarios show its effectiveness and robustness against the data noise.

An Iterative Direct Sampling Method for Reconstructing Moving Inhomogeneities in Parabolic Problems

TL;DR

This work introduces the Iterative Direct Sampling Method (IDSM) for reconstructing moving inhomogeneities in parabolic inverse problems from limited boundary data. By linking boundary measurements to the inclusion through a semi-linear forward map and constructing a probing kernel via a low-rank-augmented resolvent, IDSM directly images inclusions and handles nonlinear and mixed-type anomalies. The method operates on time segments with a fixed-point iteration to address nonlinearity, employs an inner product for stability, and updates a low-rank kernel on the fly, achieving robust reconstructions with only a handful of PDE solves per data set. Numerical experiments across merging/splitting, mixed-type, nonlinear, and fading/diminishing inclusions demonstrate accurate tracking under noise and limited data, underscoring the approach’s practical potential for time-dependent parabolic imaging.

Abstract

We propose in this work a novel iterative direct sampling method for imaging moving inhomogeneities in parabolic problems using boundary measurements. It can efficiently identify the locations and shapes of moving inhomogeneities when very limited data are available, even with only one pair of lateral Cauchy data, and enjoys remarkable numerical stability for noisy data and over an extended time horizon. The method is formulated in an abstract framework, and is applicable to linear and nonlinear parabolic problems, including linear, nonlinear, and mixed-type inhomogeneities. Numerical experiments across diverse scenarios show its effectiveness and robustness against the data noise.

Paper Structure

This paper contains 20 sections, 1 theorem, 65 equations, 5 figures, 1 table, 1 algorithm.

Key Result

Lemma 1

\newlabellemma10 For $y \in L^{2}(0, \delta t; L^{2}(\Gamma))$, let $z\in L^{2}(0,\delta t;H^{1}(\Omega))$ solve the backward heat equation: Then, $z(x, t) = ( F(\cdot; x, t), y )_{L^{2}(0, \delta t; L^{2}(\Gamma))} = (\mathcal{T}\mathcal{F})^{*} y$.

Figures (5)

  • Figure 1: Numerical results for Example \ref{['exam:1']}, where two moving inclusions merge into one and then split. Rows 1 and 2 are for $\varepsilon = 5\%$, while rows 3 and 4 are for $\varepsilon = 10\%$. \newlabelfig10
  • Figure 2: Numerical results for Example \ref{['exam:2']}, which involves mixed-type inhomogeneities. The top and bottom two rows are for the conductivity and potential inclusions, respectively. \newlabelfig20
  • Figure 3: Numerical results for Example \ref{['exam:3']}, which involves nonlinear inhomogeneities. \newlabelfig30
  • Figure 4: Numerical results for Example \ref{['exam:4']}, which involves fading inhomogeneities. \newlabelfig40
  • Figure 5: Numerical results for Example \ref{['exam:5']}, which involves diminishing inhomogeneities. \newlabelfig50

Theorems & Definitions (3)

  • Lemma 1
  • Proof 1
  • Example 4.1