偏微分方程数值解：有限差分法的稳定性与收敛性分析

引言

偏微分方程（PDE）是描述物理场量时空演化的核心数学语言。大多数工程PDE无解析解，需依赖数值方法。本文系统分析有限差分法（FDM）的精度、稳定性和收敛性，重点讨论Lax等价性定理的应用。

理论基础

差分格式的构造

考虑一维热传导方程：

$$\frac{\partial u}{\partial t} = \alpha \frac{\partial^2 u}{\partial x^2}, \quad x \in (0, L), t > 0$$

FTCS（Forward-Time Central-Space）离散格式：

$$\frac{u_j^{n+1} - u_j^n}{\Delta t} = \alpha \frac{u_{j+1}^n - 2u_j^n + u_{j-1}^n}{\Delta x^2}$$

截断误差与精度

将精确解 $u(x_j, t_n)$ 代入差分格式，Taylor 展开得截断误差（truncation error）：

$$\tau = \frac{\Delta t}{2} \frac{\partial^2 u}{\partial t^2} - \alpha \frac{\Delta x^2}{12} \frac{\partial^4 u}{\partial x^4} + \mathcal{O}(\Delta t^2, \Delta x^4)$$

FTCS 格式时间一阶精度、空间二阶精度。

Von Neumann 稳定性分析

设解为 $u_j^n = \xi^n e^{ikj\Delta x}$，代入 FTCS 格式得放大因子：

$$g(k) = 1 - 2r(1 - \cos k\Delta x), \quad r = \frac{\alpha \Delta t}{\Delta x^2}$$

稳定条件 $|g| \le 1$ 导出 FTCS 的 CFL 条件：

$$r = \frac{\alpha \Delta t}{\Delta x^2} \le \frac{1}{2}$$

Lax 等价定理

对于适定的线性初值问题，相容 + 稳定 ⟺ 收敛。此定理保证：只要差分格式正确逼近原方程（相容性）且误差不随时间放大（稳定性），数值解必然收敛到精确解。

数值实现

import numpy as np

def heat_fdm_ftcs(u0, dx, dt, alpha, Nt):
    """FTCS求解1D热传导，stability check included"""
    r = alpha * dt / dx**2
    if r > 0.5:
        print(f"Warning: r={r:.3f} > 0.5, may be unstable")
    u = u0.copy()
    Nx = len(u)
    for t in range(Nt):
        u_new = u.copy()
        for j in range(1, Nx-1):
            u_new[j] = u[j] + r * (u[j+1] - 2*u[j] + u[j-1])
        u = u_new
    return u

Nx = 50; L = 1.0; dx = L/(Nx-1); dt = 0.0004; alpha = 1.0
u0 = np.sin(np.pi * np.linspace(0, L, Nx))
u_final = heat_fdm_ftcs(u0, dx, dt, alpha, Nt=1000)

常见差分格式对比

格式	时间精度	空间精度	稳定性	类型
FTCS	一阶	二阶	条件稳定（$r \le 1/2$）	显式
BTCS（向后 Euler）	一阶	二阶	无条件稳定	隐式
Crank-Nicolson	二阶	二阶	无条件稳定	隐式
Du Fort-Frankel	二阶	二阶	无条件稳定（有条件相容）	显式

结论

FDM 是 PDE 数值解法中最直观的方法。其缺点在于对复杂几何的适应力弱（需结构化网格），这恰是 FEM 的优势所在。Crank-Nicolson 隐式格式（无条件稳定、二阶精度）是热传导问题的工程标准。