[论文解读] Logarithmic stability of parabolic Cauchy problems

Mourad Choulli, Masahiro Yamamoto|arXiv (Cornell University)|Feb 21, 2017

Numerical methods in inverse problems参考文献 17被引用 10

一句话总结

本文在利普希茨区域中建立了抛物型柯西问题的全局对数稳定性估计，证明了在 $ Q = \Omega \times (t_0,t_1) $ 上满足 $ Lu = 0 $ 的解，其唯一性由子边界 $ \Gamma_0 \times (t_0,t_1) $ 上的柯西数据确定，稳定性以三重对数衰减率 $ \Psi_{\varrho_0,\mu}(\varrho) = (\ln \ln \ln |\ln \varrho|)^{-\mu} $ 量化。该结果适用于 $ \mathscr{Z}(Q) $ 中的光滑解，且在最小几何假设下将唯一性结果推广至全域区域控制。

ABSTRACT

The uniqueness of parabolic Cauchy problems is nowadays a classical problem and since Hadamard \cite{Ha} these kind of problems are known to be ill-posed and even severely ill-posed. Until now there are only few partial results concerning the quantification of the stability for parabolic Cauchy problems. In the present article, we bring the complete answer to this issue, provided that the space domain has finite diameter with respect to the geodesic distance and assuming that solutions are sufficiently smooth.

研究动机与目标

• 量化抛物型柯西问题在整个区域 $ Q = \Omega \times (t_0,t_1) $ 上解的唯一性，而非仅局部唯一性。
• 为 $ Q $ 上的 $ Lu = 0 $ 建立全局稳定性估计，其中 $ L $ 为具有可测系数的二阶抛物型算子。
• 推导出在对数衰减意义下最优的稳定性速率，即使数据仅在小的子边界 $ \Gamma_0 \subset \partial\Omega $ 上给出亦成立。
• 通过考虑 $ \mathscr{Z}^j(Q) $ 中的解，将稳定性结果推广至高阶时间导数情形，确保在正则性假设下的鲁棒性。

提出的方法

• 作者采用专为具有 $ W^{1,\infty} $ 系数的抛物型算子设计的卡莱曼不等式，以边界数据控制解。
• 推导出三圆柱插值不等式，以连接空间与时间正则性，从而实现对 $ \Omega \times (t_0,t_1) $ 上 $ H^1 $-范数的控制。
• 证明依赖于对哈代型不等式的创新应用，以及基于区域利普希茨正则性的几何论证，确保测地距离的统一有界性。
• 通过数据范数的双重对数变换构造稳定性估计，从而导出三重对数稳定性函数 $ \Psi_{\varrho_0,\mu}(\varrho) $。
• 通过归纳法将方法扩展至高阶时间导数，利用 $ Lu = 0 $ 蕴含 $ L\partial_t^j u = 0 $ 的事实，实现对 $ H^{j+3} $-范数的递归控制。
• 关键的技术工具是使用截断函数与变量变换，以实现问题的局部化，并在有界区域中应用卡莱曼估计。

实验结果

研究问题

• RQ1 抛物型柯西问题解的唯一性是否可在全球区域 $ Q $ 上被量化，而非仅在局部？
• RQ2 当数据仅在小的子边界 $ \Gamma_0 \subset \partial\Omega $ 上给出时，抛物型柯西问题的最优稳定性速率是什么？
• RQ3 稳定性估计如何依赖于区域的正则性及抛物型算子系数的性质？
• RQ4 稳定性结果是否可推广至具有更高阶时间正则性的解，如 $ u \in \mathscr{Z}^j(Q) $？
• RQ5 是否可能在不依赖 $ \Gamma_0 $ 具有特殊几何构型（如等值面或光滑空洞边界的边界）的条件下，实现全局稳定性估计？

主要发现

• 本文建立了全局对数稳定性估计：$ C\|u\|_{L^2((t_0,t_1),H^1(\Omega))} \leq \|u\|_{\mathscr{Z}(Q)} \Psi_{\varrho_0,\mu}\left(\frac{\mathcal{C}(u,\Gamma_0)}{\|u\|_{\mathscr{Z}(Q)}}\right) $，其中 $ \Psi_{\varrho_0,\mu}(\varrho) = (\ln \ln \ln |\ln \varrho|)^{-\mu} $ 对于 $ \varrho \leq \varrho_0 $ 成立。
• 该稳定性速率在对数衰减意义下是最优的，与已知不适定问题的对数衰减一致，且在给定假设下无法获得更优速率。
• 常数 $ \mu = \min\{\alpha, s\}/4 $ 依赖于霍尔德正则性参数 $ \alpha $ 与时间局部化参数 $ s \in (0,1/2) $，反映了正则性与衰减速率之间的权衡。
• 该结果对任意利普希茨区域 $ \Omega $ 成立，无需 $ \Gamma_0 $ 为等值面或空洞边界的特殊几何结构，这使其区别于以往工作。
• 稳定性估计被推广至高阶时间导数情形：对于 $ u \in \mathscr{Z}^j(Q) $，有界为 $ C\|u\|_{H^j((t_0,t_1),H^1(\Omega))} \leq (j+1)\|u\|_{\mathscr{Z}^j(Q)} \Psi_{\varrho_0,\mu}\left(\frac{\mathcal{C}_j(u,\Gamma_0)}{\|u\|_{\mathscr{Z}^j(Q)}}\right) $。
• 证明依赖于一项新颖的几何论证，表明利普希茨区域中的内在测地距离被欧氏距离一致有界，这对卡莱曼估计的全局成立至关重要。