有限差分法(finite difference method, FDM)用泰勒展开把导数替换成相邻节点值的加权差,是最直观、最古老的 PDE 数值方法。规则网格上它给出干净的稀疏模板——正是这个"规则稀疏"让有限差分离散的算子成为量子算法(HHL、哈密顿量模拟)最容易接入的对象。
方腔问题的有限差分网格示例
图:方腔问题的有限差分网格示例——贴壁加密的非均匀网格可在边界层用更少节点达到相同精度(来源:QC076)

泰勒展开与差分公式

对光滑函数 ,在节点 处展开:

两式相减(消除偶数阶项)得中心差分:

相加得二阶导数:

一阶迎风()只有一阶精度但带数值耗散,能稳定对流占优问题——精度与稳定性的经典权衡。差分近似可以装配成矩阵:节点值向量 、导数矩阵 ,整个离散就是把 PDE 变成 或 。

显式、隐式与 CFL 条件

以一维平流方程 ()为例。显式格式(时间前向、空间迎风)

做 von Neumann 分析:设 ,放大因子

这就是 CFL 条件(Courant-Friedrichs-Lewy 1928)的实例:信息传播的数值依赖域必须包含物理依赖域,;多维多波速时 (声速限制见 dimensionless-numbers)。隐式格式(如 Crank-Nicolson、后向欧拉)无条件稳定,但每步要解线性系统——把成本从”每步便宜但步数多”搬到”每步贵但步数少”,求解器见 multigrid。

泊松方程与七点模板

三维泊松方程 在均匀网格上差分离散为七点模板:

装配出的矩阵 对称、负定、极稀疏(每行 7 个非零元)。它有显式的正弦特征分解(离散正弦变换对角化):

条件数 ——直接迭代收敛慢,引出多重网格(见 multigrid)与预条件技术。

与量子算法的联系

七点模板是量子泊松求解器的标准离散对象:Cao 等的电路设计即针对此类 stencil 矩阵做相位估计Cao et al. 2013 (New J. Phys.),HHL 原文的复杂度分析也以稀疏、良好条件化的 stencil 为前提Harrow et al. 2009 (PRL)(见 quantum-poisson-solver)。更一般地,FDM 的 型半离散可直接映射为稀疏哈密顿量模拟(见 hamiltonian-simulation);差分算子的平移不变结构还允许块编码高效实现(见 block-encoding)。顶盖方腔的流函数-涡量 FDM 求解(Ghia 基准复现,见 classical-benchmarks)是这条路线的经典验证场Mao 2024 (arXiv)。

延伸阅读