28、扩散方程离散化与求解方法详解

扩散方程离散化与求解方法详解

1. 离散模型稳定性分析

在对扩散方程进行离散化处理时，我们遇到了一些有趣的现象。对于某些显式格式，得出了 (D > 0 \Rightarrow|\lambda ± | < 1, \forall\varphi) 的结果，这意味着该显式格式具有无条件稳定性。然而，这一结果令人怀疑，因为它明显违反了 Dahlquist 的第二稳定性屏障。

为了深入分析，我们需要研究格式的收敛性，进而考虑其一致性。通过推导，得到了等价的微分方程：
(\frac{\partial\psi}{\partial t} + \kappa(\frac{\Delta t}{\Delta x})^2\frac{\partial^2\psi}{\partial t^2} = \kappa\frac{\partial^2\psi}{\partial x^2})
这表明我们实际上求解的是一个修正后的微分方程，而非原始方程，明显违反了一致性条件。不过，在实际应用中，若取扩散数 (D = \kappa\frac{\Delta t}{\Delta x^2}) 为常数，DuFort–Frankel 格式几乎是一致的，因为 (\Delta t) 比 (\Delta x) 更快趋近于 0。

2. 扩散算子的特征谱

接下来，我们研究对应不同边界条件下矩阵 (A = \frac{\kappa}{\Delta x^2}(E - 2I + E^{-1})) 的特征谱。
- Dirichlet 边界条件 ：设 (\psi_0(0) = a)，(\psi_N = b)，端点方程为 (\frac{d\psi_1}{