[论文解读] Lipschitz stability in an inverse problem for the wave equation

Lucie Baudouin|arXiv (Cornell University)|Jun 8, 2011

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

一句话总结

本文建立了涉及波动方程的逆问题的Lipschitz稳定性，证明了在 $ \partial_{tt}u - \Delta u + p(x)u = 0 $ 中，势函数 $ p(x) $ 可以从 $ \Gamma_0 \subset \partial\Omega $ 上的法向导数 $ \partial_\nu u $ 测量值唯一且稳定地重构。关键结果通过全局Carleman估计得出，该估计在 $ \Gamma_0 $ 和观测时间 $ T $ 满足几何条件时，可证明稳定性，并显式依赖于数据正则性和初始数据的下界。

ABSTRACT

We are interested in the inverse problem of the determination of the potential $p(x), x\inΩ\subset\mathbb{R}^n$ from the measurement of the normal derivative $\partial_νu$ on a suitable part $Γ_0$ of the boundary of $Ω$, where $u$ is the solution of the wave equation $\partial_{tt}u(x,t)-Δu(x,t)+p(x)u(x,t)=0$ set in $Ω imes(0,T)$ and given Dirichlet boundary data. More precisely, we will prove local uniqueness and stability for this inverse problem and the main tool will be a global Carleman estimate, result also interesting by itself.

研究动机与目标

• 建立波动方程中通过边界法向导数 $ \partial_\nu u $ 测量值反演势函数 $ p(x) $ 的局部Lipschitz稳定性。
• 在观测区域 $ \Gamma_0 $ 和观测时间 $ T $ 满足几何条件时，证明唯一性和稳定性结果，确保重构所需的信息充分。
• 通过全局Carleman估计提供稳定性的直接证明，避免先前工作中使用的紧致性-唯一性论证。
• 弱化对解正则性的假设，并证明仅在边界子集 $ \Gamma_0 $ 上的测量已足够，而非整个边界。

提出的方法

• 为波动算子 $ \partial_{tt} - \Delta + q(x) $ 推导全局Carleman估计，使用依赖于 $ |x - x_0|^2 $ 和时间 $ t $ 的权重函数 $ \varphi $，其中 $ x_0 \notin \Omega $。
• 将Carleman估计应用于逆问题的线性化版本，其中解 $ u $ 满足齐次初值和边界条件，并由源项 $ f(x)R(x,t) $ 驱动。
• 使用截断函数 $ \chi $ 局域化解在空间和时间上的行为，从而控制Carleman不等式中的能量和通量项。
• 通过估计源项 $ f $ 的加权 $ L^2 $-范数与 $ \Gamma_0 \times (0,T) $ 上的边界通量 $ \partial_\nu u $ 之间的关系，利用Carleman估计建立类似可观测性的不等式。
• 结合能量估计和权重函数的性质，吸收低阶项，推导出 $ \|f\|_{L^2} $ 与边界测量 $ \|\partial_\nu u\|_{L^2(\Gamma_0 \times (0,T))} $ 之间的最终不等式。
• 最终得出Lipschitz稳定性估计：$ \|f\|_{L^2(\Omega)} \leq C \|\partial_\nu u\|_{L^2(\Gamma_0 \times (0,T))} $，其中常数 $ C $ 依赖于 $ m, T, K, r $。

实验结果

研究问题

• RQ1 在足够观测时间 $ T $ 下，是否能从 $ \Gamma_0 \subset \partial\Omega $ 上法向导数 $ \partial_\nu u $ 的测量值唯一确定波动方程中的势函数 $ p(x) $？
• RQ2 是否可以为该逆问题建立Lipschitz稳定性估计，将 $ q - p $ 的 $ L^2 $-范数以边界通量差 $ \partial_\nu y[q] - \partial_\nu y[p] $ 的 $ L^2 $-范数为界？
• RQ3 全局Carleman估计是否允许在弱正则性假设和仅使用子集 $ \Gamma_0 $ 测量（而非整个边界）下获得稳定性结果？
• RQ4 在 $ \Gamma_0 $ 和 $ T $ 的何种几何与时间条件下，波动方程的逆问题具有局部Lipschitz稳定性？
• RQ5 稳定性常数是否可显式地用数据正则性和初始数据下界表示？

主要发现

• 本文证明了逆问题的局部Lipschitz稳定性：对于 $ q, p \in L^\infty_{\leq m}(\Omega) $，差值 $ \|q - p\|_{L^2(\Omega)} $ 被常数倍的 $ \|\partial_\nu y[q] - \partial_\nu y[p]\|_{H^1(0,T;L^2(\Gamma_0))} $ 控制，其中常数 $ C = C(m, T, K, r) $ 依赖于 $ L^\infty $-界 $ m $、时间 $ T $、解正则性 $ K $ 和初始数据下界 $ r $。
• 稳定性结果在如下几何条件下成立：存在 $ x_0 \notin \Omega $，使得 $ \Gamma_0 \supset \{x \in \partial\Omega : (x - x_0) \cdot \nu(x) \geq 0\} $，且 $ T > \sup_{x \in \Omega} |x - x_0| $。
• 该证明依赖于波动算子的全局Carleman估计，使得可使用边界子集 $ \Gamma_0 $ 上的测量，而非整个边界。
• 该方法避免了先前工作中使用的紧致性-唯一性论证，通过Carleman估计直接证明稳定性。
• Carleman估计应用于初值和边界条件为零的线性化问题，其中源项 $ f(x)R(x,t) $ 由边界通量 $ \partial_\nu u $ 恢复。
• 最终的稳定性估计为 $ \|f\|_{L^2(\Omega)} \leq C \|\partial_\nu u\|_{L^2(\Gamma_0 \times (0,T))} $，其中 $ C $ 依赖于 $ m, T, K, r $，且该结果通过线性化方法可推广至非线性逆问题。