[论文解读] Inverse Problem of Finding the Time-dependent Coefficient of Heat Equation from Integral Overdetermination Condition Data

Mansur I. Ismailov, Fatma Kanca|arXiv (Cornell University)|Apr 30, 2010

Numerical methods in inverse problems参考文献 7被引用 14

一句话总结

本文研究了一维热方程中具有非局部边界条件和积分过定条件的逆问题，旨在识别时变热扩散系数 $ a(t) $ 和温度分布 $ u(x,t) $。通过傅里叶方法与不动点原理，作者在特定光滑性与相容性条件下证明了经典解的存在性与唯一性，并提出了一种基于Crank-Nicolson格式与预测-校正迭代的数值方法，通过一个数值例子验证了该方法的高精度。

ABSTRACT

In this paper we consider the problem of simultaneously determining the time-dependent thermal diffusivity and the temperature distribution in one-dimensional heat equation in the case of nonlocal boundary and integral overdetermination conditions. We establish conditions for the existence and uniqueness of a classical solution of the problem under considerations. We present some results on the numerical solution with an example.

研究动机与目标

• 建立在具有非局部边界条件和积分过定条件的一维热方程中，识别 $ a(t) $ 和 $ u(x,t) $ 的逆问题存在经典解的条件。
• 分析解对问题数据的连续依赖性。
• 基于Crank-Nicolson格式与迭代预测-校正方法，开发求解逆问题的数值算法。
• 通过一个具有已知精确解的数值例子验证该方法，展示其收敛性与精度。

提出的方法

• 利用傅里叶方法将问题转化为不动点问题，借助空间域内的特征函数展开。
• 分析与非局部边界条件相关的谱问题，构造特征函数与相关函数的双正交系统。
• 将解表示为特征函数的级数展开形式，并通过积分过定条件 $ \int_0^1 u(x,t)dx = E(t) $ 推导时变系数 $ a(t) $。
• 设计预测-校正迭代格式：从 $ E(t) $ 的时间导数与 $ u $ 的空间梯度预测系数 $ a^{j+1} $，再通过Crank-Nicolson格式进行校正。
• 应用具有空间与时间二阶精度的Crank-Nicolson格式，采用高斯消去法迭代求解系统。
• 数值实现采用有限差分半径离散化，其中 $ h = 0.005 $，$ \tau = h/4 $，并通过连续迭代之间的差值监控收敛性。

实验结果

研究问题

• RQ1 在具有非局部边界条件与积分过定条件的一维热方程中，识别时变热扩散系数 $ a(t) $ 与温度 $ u(x,t) $ 的逆问题在何种条件下存在经典解？
• RQ2 为使解的存在性与唯一性成立，初始函数、源项与过定函数需满足哪些必要的光滑性与相容性条件？
• RQ3 如何从积分条件 $ \int_0^1 u(x,t)dx = E(t) $ 重构时变系数 $ a(t) $，且该重构过程的稳定性如何？
• RQ4 何种数值格式能确保在求解具有非局部条件的逆问题时具备二阶精度与稳定性？
• RQ5 与精确解相比，所提出的迭代数值方法在恢复 $ a(t) $ 与 $ u(x,t) $ 时的精度如何？

主要发现

• 在假设 (A1)–(A3) 下，逆问题在 $ Q_T $ 中存在经典解，这些假设包括对数据的光滑性、相容性与正性条件。
• 解在 $ Q_{T_0} $ 中唯一确定，其中 $ T_0 < T $，且 $ T_0 $ 依赖于问题数据。
• 基于Crank-Nicolson与迭代预测-校正格式的数值方法，能以高精度成功恢复 $ a(t) $ 与 $ u(x,t) $。
• 对于测试例子 $ T = 1/4 $，$ h = 0.005 $，$ \tau = h/4 $，该方法与精确解 $ a(t) = \frac{1}{(2\pi)^2} + \exp(4t) $ 和 $ u(x,t) = (1-x)\sin(2\pi x)\exp(-t) $ 高度吻合。
• 迭代格式收敛，且在每个时间步，当连续迭代之间的差值低于预设容差时，即接受计算得到的 $ a(t) $ 与 $ u(x,t) $ 值。
• 数值结果（表1与表2）证实该方法具有稳定性与二阶精度。