[论文解读] Reconstruction of a source domain from the Cauchy data: II. Three dimensional case

Masaru Ikehata|arXiv (Cornell University)|Jul 12, 2021

Numerical methods in inverse problems参考文献 15被引用 3

一句话总结

本文将包围法扩展至三维由固定波数下Helmholtz方程控制的反源与障碍物问题。通过使用复几何光学解和振荡积分，建立了一个框架，以重构具有锥形奇点的源域支撑，并证明了四面体和圆锥形奇点下一个关键复系数的非零性。主要贡献是基于单组柯西数据，给出了源域几何结构与强度的显式重构公式，该公式在波数和源正则性满足一定条件时成立。

ABSTRACT

This paper is concerned with reconstruction issue of some typical inverse problems and consists of three parts. First a framework of the enclosure method for an inverse source problem governed by the Helmholtz equation at a fixed wave number in three dimensions is introduced. It is based on the nonvanishing of the coefficient of the leading profile of an oscillatory integral over a domain having a conical singularity. Second an explicit formula of the coefficient for a domain having a circular cone singularity and its implication under the framework are given. Third, an application under the framework to an inverse obstacle problem governed by an inhomogeneous Helmholtz equation at a fixed wave number in three dimensions is given.

研究动机与目标

• 使用单组柯西数据解决三维Helmholtz方程的反源与障碍物问题。
• 基于包围法构建框架，以提取未知源域的几何与物理信息。
• 证明振荡积分在锥形区域上产生的复系数的非零性，这是重构的关键。
• 在几何与正则性假设下，建立源域支撑函数与源强度的显式重构公式。
• 将先前的二维结果扩展至三维情形，特别针对具有四面体或圆锥形奇点的区域。

提出的方法

• 使用形如 $ v \sim e^{x \cdot z} $ 的复几何光学（CGO）解，其中 $ z = \tau(\omega + i\vartheta) $，且 $ \tau \to \infty $，以探测源域。
• 分析振荡积分 $ \int_D \rho(x) e^{\tau x \cdot (\omega + i\vartheta)} dx $ 在 $ \tau \to \infty $ 时的渐近行为，识别主阶系数。
• 引入一个复常数 $ C(p,\omega)(\delta, Q, \vartheta) $，用于表征以 $ Q $ 为底面、顶点为 $ p $ 的锥体上振荡积分的主阶轮廓。
• 证明在三维四面体锥形区域上，$ C(p,\omega)(\delta, Q, \vartheta) $ 非零，从而确保该系数可从数据中识别。
• 通过指标函数及其导数，推导出源域支撑函数 $ h_D(\omega) $ 与乘积 $ \tilde{\rho}(p(\omega)) u(p(\omega)) V(\theta) $ 的显式重构公式。
• 将该框架应用于非齐次Helmholtz方程，实现从边界柯西数据对源域及其强度的重构。

实验结果

研究问题

• RQ1 如何将包围法从二维扩展至由Helmholtz方程控制的三维反源问题？
• RQ2 在三维空间中，何种条件可确保振荡积分在锥形区域上的渐近展开主系数非零？
• RQ3 如何从单组柯西数据重构具有圆锥形奇点的源域支撑函数？
• RQ4 可导出哪些关于源强度与几何结构的显式公式，利用指标函数及其导数？
• RQ5 波数 $ k $ 如何影响源域及其性质重构的可行性？

主要发现

• 对于任意三维四面体锥形区域，控制振荡积分主阶轮廓的复系数 $ C(p,\omega)(\delta, Q, \vartheta) $ 非零。
• 通过极限 $ \lim_{\tau \to ∞} \frac{I'_{\omega,\vartheta_j}(\tau)}{I_{\omega,\vartheta_j}(\tau)} $ 推导出支撑函数 $ h_D(\omega) = p(\omega) \cdot \omega $ 的显式重构公式，其中 $ I_{\omega,\vartheta} $ 为指标函数。
• 乘积 $ \tilde{\rho}(p(\omega)) u(p(\omega)) V(\theta) $（结合源强度、解值与立体角）可通过公式 (4.7)、(4.8) 和 (4.9) 从柯西数据中重构。
• 在 $ \omega \approx n $（顶点处的法向量）的条件下，当且仅当 $ \omega = n $ 时，函数 $ I(\omega, \vartheta) $ 为常数，从而可识别法向方向。
• 在小 $ k $ 条件下，推导出源域内 $ |u(x)| $ 的显式下界，确保重构所需的条件 $ |u(p(\omega))| > 0 $ 成立。
• 在适当假设下，若 $ \tilde{\rho}(p(\omega)) $ 为实数，则可模 $ 2\pi $ 恢复复值解 $ u(p(\omega)) $ 的相位，从而实现对源强度的完整重构。