[论文解读] Identification of a boundary influx condition in a one-phase Stefan problem

Chifaa Ghanmi, Saloua Mani Aouadi+1 more|arXiv (Cornell University)|Feb 21, 2020

Numerical methods in inverse problems参考文献 32被引用 11

一句话总结

本文通过从测量的移动界面位置反演边界热流入条件，解决了单相Stefan反问题。利用积分方程与解析函数延拓方法，建立了对数稳定性估计，并提出了一种Tikhonov正则化数值算法，在噪声数据下相对误差低于1.3%，验证了理论上的不适定性以及在t=0附近的数值鲁棒性。

ABSTRACT

We consider a one-dimensional one-phase inverse Stefan problem for the heat equation. It consists in recovering a boundary influx condition from the knowledge of the position of the moving front, and the initial state. We derived a logarithmic stability estimate that shows that the inversion is ill-posed. The proof is based on integral equations and unique continuation for holomorphic functions. We also proposed a direct algorithm with a regularization term to solve the nonlinear inverse problem. Several numerical tests using noisy data are provided with relative errors analysis.

研究动机与目标

• 从移动界面位置和初始温度信息中，恢复单相Stefan问题的边界热流入条件。
• 为反问题建立理论稳定性估计，证明其不适定性。
• 基于正则化线性积分方程，开发一种求解非线性反问题的数值算法。
• 通过带噪声数据的数值实验验证该方法，并分析相对误差。

提出的方法

• 将反Stefan问题表述为具有自由边界的热方程非标准柯西问题。
• 利用积分方程与解析函数的唯一延拓性质，推导出对数稳定性估计。
• 将反问题转化为第一类正则化线性积分方程。
• 应用Tikhonov正则化，引入参数λ以稳定不适定问题的求解。
• 采用中点法则与高斯-勒让德求积法对积分方程进行离散化，实现数值求解。
• 针对三个示例进行数值测试，包括解析与非解析的边界函数，在不同高斯噪声水平下进行评估。

实验结果

研究问题

• RQ1 能否为从移动界面位置反演边界热流入的单相Stefan反问题，推导出对数稳定性估计？
• RQ2 反问题的不适定性在稳定性方面如何表现，是否可量化？
• RQ3 基于正则化积分方程的直接数值算法能否在噪声数据下有效重建边界热流入？
• RQ4 重建精度如何随噪声水平和t=0附近区域的变化而变化？
• RQ5 正则化参数λ在数值解中对相对误差的影响程度如何？

主要发现

• 反问题表现出对数稳定性，证实其不适定性，源于对自由边界数据微小扰动的高度敏感性。
• Tikhonov正则化数值算法在所有测试噪声水平下均成功重建边界热流入，相对误差低于1.3%。
• 对于解析的边界热流入函数，重建精度在t=0附近较低，且随着噪声减小，稳定性提高。
• 在5%高斯噪声下，方法表现出鲁棒性，示例3中相对误差从0.0040增至0.0124。
• 数值结果验证了理论稳定性估计，尤其在λ ≈ 10⁻²时表现更优，具有更好的性能与稳定性。
• 由于正则化与离散化的影响，在t=0附近观察到弱数值奇异性，导致该区域相对误差略有增加。