[论文解读] Determining the potential in a wave equation without a geometric condition. Extension to the heat equation

Kaïs Ammari, Mourad Choulli+1 more|arXiv (Cornell University)|Sep 30, 2015

Numerical methods in inverse problems参考文献 23被引用 12

一句话总结

本文在不施加测量子边界几何条件的前提下，建立了从边界测量反演波动方程中势函数的对数稳定性估计。通过利用谱分解和罗比安诺（Robbiano）提出的插值不等式，而非经典可观测性不等式，作者获得了如下形式的稳定性：$ \|q_0 - q\|_{L^2} \leq \Psi(\|\widetilde{\mathscr{C}}_q^\tau - \widetilde{\mathscr{C}}_{q_0}^\tau\|) $，其中 $ \Psi(\gamma) = |\ln \gamma|^{-1/(8+2n)} + \gamma $。该方法被进一步推广至热方程，得到了类似的对数稳定性估计。

ABSTRACT

We prove a logarithmic stability estimate for the inverse problem of determining the potential in a wave equation from boundary measurements obtained by varying the first component of the initial condition. The novelty of the present work is that no geometric condition is imposed to the sub-boundary where the measurements are made. Our results improve those obtained by the first and second authors in [2]. We also show how the analysis for the wave equation can be adapted to an inverse coefficient problem for the heat equation

研究动机与目标

• 建立从边界测量反演波动方程中势函数的对数稳定性估计。
• 消除对测量子边界上几何条件（如非捕获或非凹性）的依赖。
• 将分析从波动方程推广至热方程的反系数问题。
• 提供一种仅依赖于 Dirichlet-to-Neumann 算子差值的范数、而不依赖于区域几何假设的稳定性结果。

提出的方法

• 利用波动算子的谱分解，将解表示为特征函数与特征值的组合。
• 在缺乏几何约束的条件下，用罗比安诺（2009）提出的插值不等式替代经典可观测性不等式。
• 以插值不等式的形式应用一种卡莱曼型估计，以控制低频分量。
• 通过在 Sobolev 空间 $ H^1_0(\Omega) $ 和 $ \mathcal{H} = H^1_0(\Omega) \cap H_\Delta(\Omega) $ 中的能量估计与对偶论证，推导出稳定性估计。
• 通过分析半群解及其在边界上的法向导数，将该方法适配至热方程。
• 使用算子 $ \mathcal{N}_q $ 将初值映射到解在边界上的法向导数，并通过估计其范数差来恢复势函数。

实验结果

研究问题

• RQ1 在测量子边界 $ \Upsilon \subset \Gamma $ 上不施加几何条件时，是否能为波动方程中的反势问题建立对数稳定性？
• RQ2 在缺乏几何假设的条件下，如何用其他方法替代经典可观测性不等式以推导稳定性估计？
• RQ3 用于波动方程的方法在多大程度上可被推广至热方程的反系数问题？
• RQ4 在缺乏几何控制的条件下，可达到的最优稳定性速率是什么？其与测量算子范数差的关系如何？

主要发现

• 在测量子边界 $ \Upsilon \subset \Gamma $ 上不施加任何几何条件的前提下，建立了波动方程中反势问题的对数稳定性估计，其稳定性速率形式为 $ \|q_0 - q\|_{L^2} \leq \Psi(\|\widetilde{\mathscr{C}}_q^\tau - \widetilde{\mathscr{C}}_{q_0}^\tau\|) $，其中 $ \Psi(\gamma) = |\ln \gamma|^{-1/(8+2n)} + \gamma $。
• 该稳定性估计通过谱分解与罗比安诺（2009）提出的插值不等式推导得出，从而避免了传统方法中为保证可观测性而必须施加的几何条件。
• 该方法被推广至热方程，得到了类似的对数稳定性估计 $ \|q - q_0\|_{L^2} \leq \Theta(\|\mathcal{N}_q - \mathcal{N}_{q_0}\|) $，其中 $ \Theta(\gamma) = |\ln \gamma|^{-1/(1+4n)} + \gamma $。
• 稳定性估计在 $ L^\infty(\Omega) $ 中有界势函数集合内一致成立，且界仅依赖于区域、测量子边界及势函数的 $ L^\infty $-范数上界。
• 分析依赖于 Sobolev 空间 $ \mathcal{H} $ 和 $ \mathcal{H}_0 $ 中的估计，确保解的法向导数属于 $ H^1((0,\tau); L^2(\Upsilon)) $，从而可应用迹估计与能量估计。
• 关键技术突破在于：在缺乏几何控制的条件下，利用插值不等式控制低频模态，从而在无需非捕获或非凹性假设的前提下实现稳定性推导。