PDE 离散的一致性、稳定性与收敛性

形式陈述 ​

离散接口与一致性 ​

一个 PDE 数值问题必须连同区域、初边值条件和解空间一起离散。把连续问题抽象写成

$$\mathcal{L}u=f,\mathcal{B}u=g,u(0)=u_0$$

时，$\mathcal{B}$ 与初值项是否出现取决于问题类型。给定空间尺度 $h$ 和时间步长 $\tau$，有限差分格式公理库偏微分方程有限差分法Finite-difference method for PDE · Finite-difference PDE discretization · Difference stencil在整张网格上用差分 stencil 替换偏微分算子，并连同边界闭合组装为稀疏离散系统。在离散空间 $X_h$ 中产生网格解 $U_{h,\tau }$。限制算子 $R_h$ 负责采样连续函数，重构算子 $P_h$ 把网格数据解释回连续空间；节点值、单元平均与插值系数不能在没有这两个接口时直接同 $u$ 比较。

一致性检查精确解是否近似满足离散方程。若 ${\mathcal{L}}_{h,\tau }$ 是包含内部 stencil、边界闭合和初值离散的完整算子，$Q_{h,\tau }$ 离散连续数据，则缺陷为

$$d_{h,\tau }(u)={\mathcal{L}}_{h,\tau }(R_{h}u)-Q_{h,\tau }(\mathcal{L}u,\mathcal{B}u,u_0).$$

对规定的光滑解类，若在选定离散范数中

$$\Vert d_{h,\tau }(u){\Vert }_h⟶0(h,\tau \to 0),$$

称该网格族上的格式一致。若把演化格式的单步残差除以 $\tau$ 再定义缺陷，其阶数会与未缩放版本相差一；任何 $O(h^p+{\tau }^q)$ 声明都要说明采用哪种规范。

局部内部 stencil 一致还不等于整个初边值离散一致。边界公式可能低阶，初值投影可能改变约束，曲线边界的几何近似也可能产生独立缺陷。只有把这些部分放进同一个 ${\mathcal{L}}_{h,\tau }$，局部 Taylor 展开才真正对应所求的离散问题。

稳定性 ​

稳定性控制离散问题对数据和扰动的放大，并且必须相对于指定范数公理库赋范向量空间Normed vector space带满足正定、齐次与三角不等式范数的向量空间。、网格族和有限时间区间陈述。对线性齐次一步演化

$$U^{n+1}=C_{h,\tau }{U}^n,$$

一种常用定义是：对每个固定 $T>0$，存在与足够小的 $h,\tau$ 以及 $n$ 无关的常数 $K_T$，使

$$\Vert C_{h,\tau }^n{\Vert }_{X_h\to X_h}\le K_T,0\le n\tau \le T.$$

若还有离散源项 $q^n$，相应估计应控制

$$\Vert U^n{\Vert }_h\le K_T(\Vert U^0{\Vert }_h+\tau \sum _{m=0}^{n-1}\Vert q^m{\Vert }_h),$$

并把边界扰动放入正确的数据范数。这里允许 $K_T$ 随终止时间增长，但不允许它随网格加密发散；只观察一次计算没有爆炸，不能代替这项一致界。

收敛与三者关系 ​

收敛比较离散解与连续真解。对固定有限区间 $[0,T]$，可在网格空间中要求

$$\underset{n\tau \le T}{max}\Vert U^n-R_{h}u(t_n){\Vert }_h⟶0,$$

也可在连续空间中要求

$$\underset{n\tau \le T}{max}\Vert P_h{U}^n-u(t_n){\Vert }_X⟶0.$$

两种定义是否等价取决于限制、重构和范数之间的一致控制。收敛是关于一族网格解的序列收敛公理库序列收敛Convergence of a sequence序列项最终任意接近某个极限值。，不是单张细网格图像与预期曲线相似。

在线性问题中，误差常满足“稳定传播旧误差，再加入本步缺陷”的递推。稳定估计若能统一控制传播算子，逐步缺陷便可累积成全局误差界；反之，未受控的寄生模态会把趋零缺陷放大。这个推理是分析框架，不是对所有 PDE、所有离散和所有范数都成立的无条件定理；精确的等价结论由Lax–Richtmyer 定理公理库Lax–Richtmyer 等价定理Lax-Richtmyer equivalence theorem · Lax equivalence theorem在良定线性演化初值问题、共同 Banach 空间和有限时间网格族中，刻画一致差分格式的稳定性与收敛性等价。在良定线性演化初值问题中给出。

相邻概念的边界 ​

本页的“格式稳定性”在精确算术下研究离散解算子随网格加密的扰动放大。ODE 绝对稳定公理库ODE 方法的绝对稳定域Absolute stability region · Region of absolute stability · A-stability以标量衰减测试方程的放大因子定义一步法稳定域，并区分 A-stability、L-stability、精度与其他稳定概念。研究时间积分公式作用于 $y^{\prime }=\lambda y$ 时的标量放大因子；它可以参与方法线分析，却没有包含空间边界、离散范数或非正规传播。浮点数值稳定性公理库数值稳定性Numerical stability · Backward stability以允许的小输入扰动刻画算法的有限精度行为，并与问题条件性及其他稳定性概念分开。则问实际算法是否等价于一个小扰动问题，舍入模型与 backward error 属于另一层。

连续 PDE 的良定性也要单独说明。它要求连续解对初边值数据存在、唯一且连续依赖；格式稳定性要求离散算子族拥有与网格无关的界。一个连续问题良定，并不保证任意离散格式稳定；连续问题若不良定，离散计算看似平滑也可能只是额外正则化的结果。

直觉

一致性像检查每块局部拼片的形状：把真解放上 stencil，它留下的缝隙应随网格缩小。稳定性检查拼片之间的连接会不会把微小缝隙越传越大。收敛则观察整幅离散图景在重构后是否靠近连续解；三者关注的对象分别是方程缺陷、传播增益与最终轨道。

范数决定什么叫“放大”。最大范数关注单个网格点的最坏偏差，离散 $L^2$ 范数按网格权重衡量整体能量，能量范数还会把离散导数计入。有限维中每张固定网格上的范数虽然等价，等价常数可能随 $h$ 发散，所以不能在网格族之间随意换范数后保留同一个稳定结论。

例子与边界

一致但不稳定的平流格式 ​

考虑周期区间上的常系数平流方程

$$u_t+a{u}_x=0$$

及前向时间、中心空间格式

$$U_j^{n+1}=U_j^n-\frac{\nu }{2}(U_{j+1}^n-U_{j-1}^n),\nu =\frac{a\tau }{h}.$$

将光滑精确解代入并用 PDE 消去 $u_t+a{u}_x$，归一化单步缺陷为 $O(\tau +h^2)$。它在时间上一阶、空间上二阶，因而确实一致；这个局部结论没有检查误差怎样经过许多时间步传播。

取网格振荡

$$V_j^n=G(\xi {)}^n{e}^{ij\xi }.$$

代入格式得到

$$G(\xi )=1-i\nu \sin \xi ,|G(\xi ){|}^2=1+{\nu }^2{\sin }^2\xi .$$

沿固定非零 Courant 数 $\nu$ 的周期加密网格子族，可取网格点数为 $4$ 的倍数，使 $\xi =\pi /2$ 是可达频率；一般网格也可取最接近它的离散频率。连续频率上的最坏单步增益为

$$\underset{\xi }{max}|G(\xi )|=\sqrt{1+{\nu }^2}.$$

当 $n\tau \le T$ 时，用 $\log (1+x)\le x$ 可得

$$|G{|}^n\le \exp \left(\frac{a^{2}T}{2}\frac{\tau }{h^2}\right).$$

固定非零 $\nu =a\tau /h$ 时，指数随 $h^{-1}$ 发散；局部一致性因此不足以推出收敛。一般 Fourier 工具和其他 CFL 推导由von Neumann 分析公理库von Neumann 稳定性分析von Neumann stability analysis · Fourier stability analysis · Amplification factor analysis用离散 Fourier 模态计算常系数规则网格格式的放大符号，并由全部频率的幂界判断有限时间稳定性。集中说明。

上式也显示结论怎样随网格族改变：若 $\tau /h^2$ 一致有界，增长至多为与网格无关的 $e^{C_T}$，单步增益相应降到 $1+O(\tau )$。稳定性因此不能只问“某个固定网格的 $|G|$ 是否略大于一”，而要说明 $h$ 与 $\tau$ 怎样共同趋零。

标量谱图之外 ​

非正规更新矩阵还会破坏“特征值都在单位圆内就稳定”的推断。即使每个特征值模不超过一，特征向量病态或单位圆上的 Jordan 块仍可能令 $\Vert C_{h,\tau }^n\Vert$ 随 $n$ 大幅增长。稳定性定义使用算子幂的统一范数界，正是为了保留这些特征值集合看不到的方向耦合。

非线性格式中，两条离散轨道之差不再由同一个固定线性算子推进。稳定分析可能改用离散能量、不变量区域、单调性、熵不等式或局部 Lipschitz 估计；非光滑解会使 Taylor 缺陷阶失效，激波问题还需要弱收敛与熵条件。把线性光滑理论的符号替换成非线性函数，并不会自动得到同一结论。

非一致网格需要局部尺度 $h_K$、带网格权重的范数和形状正则性。若最小单元退化，固定网格上的范数等价常数、逆不等式或离散算子界可能随网格族恶化。此时“最大网格宽度趋零”只是几何条件的一部分。

推论与应用

分析一个新格式时，可先用离散缺陷公理库离散化误差与截断误差Discretization error · Truncation error区分连续问题的离散缺陷与离散解的全局误差，并说明一致性为何不能单独保证收敛。确定一致阶，再寻找离散能量估计、最大值原则、Fourier 符号或算子范数界，最后把两者组合成重构误差估计。若数值实验没有同时报告范数、边界处理、$h$—$\tau$ 关系和固定终止时间，观测到的斜率无法与理论命题一一对应。

一般框架还负责标记理论没有覆盖的误差层。迭代线性求解未收敛会改变离散方程，浮点舍入会持续注入扰动，模型数据误差则先于离散化存在。它们可以通过稳定界传播，却不能未经说明就合并进“截断误差”或“格式不稳定”。

参考资料

• Robert D. Richtmyer and K. W. Morton, Difference Methods for Initial-Value Problems, 2nd ed., Interscience, 1967.
• John C. Strikwerda, Finite Difference Schemes and Partial Differential Equations, 2nd ed., SIAM, 2004.
• Lloyd N. Trefethen, Finite Difference and Spectral Methods for Ordinary and Partial Differential Equations, unpublished text, 1996, Ch. 4.