有限元法(finite element method, FEM)把偏微分方程改写成"弱形式"(分布意义下的积分恒等式),再用分片多项式基函数逼近。复杂几何、混合边界条件与自适应加密都在它的能力范围,使它成为结构力学的事实标准与流体力学的重要选项——对量子计算,它的"装配再求解"两段式结构直接给出了量子有限元的切入方式。

弱形式:从逐点到积分

以 Poisson 方程为例(散度型算子的原型):

强形式要求方程在每点成立;弱形式只要求对所有试验函数 (满足齐次 Dirichlet 条件)成立:两边乘 、分部积分、

分部积分把二阶导数降为一阶——解只需 正则性(弱解),Neumann 边界自动进入右端。这与物理上的虚功原理同构,是 FEM 的数学与物理双重根基。

离散化:基函数与刚度矩阵

把区域剖分为单元(三角形/四面体/六面体),在每个单元上用分片低阶多项式插值:全局基函数 取值为 1 于节点 、0 于其它节点。逼近 ,代入弱形式并取 :

即刚度矩阵:对称(自伴问题)、正定( coercivity)、稀疏(基函数支集局部, 仅当节点相邻)。装配按单元循环:每个单元算局部刚度矩阵 (可变换到参考单元上的数值积分),散射叠加进全局矩阵。误差理论由 Céa 引理支撑:——离散误差由最佳逼近控制, 阶单元给出 收敛,光滑解配合升阶(-refinement)可达指数收敛(谱元即其极限,见 spectral-methods)。

在流体力学中的位置

不可压 NS 的 FEM 需要处理压力-速度 inf-sup(LBB)稳定性——速度空间要足够丰富,否则压力棋盘振荡;标准解法是 Taylor-Hood 元(-)或稳定化格式。工程上 FEM 在固体、结构、传热与生物流体(低 Re)占优;高 Re 气动外流因激波捕捉与耗散控制难,仍是 FVM 的地盘(见 finite-volume-method、turbulence-simulation)。

量子接口:量子有限元

Montanaro 与 Pallister 指出 FEM 的两段式结构对量子算法是天然切口Montanaro & Pallister 2016 (PRA):

  1. 装配: 的分量是积分,用量子求和(振幅估计的迭代加法,见 quantum-amplitude-estimation)在对数时间内估计;
  2. 求解: 交给 HHL 型线性系统算法(见 hhl),复杂度对网格规模呈 polylog,代价在条件数 与输出精度。

完整分析与资源账见 quantum-fem。要点是:FEM 的矩阵不规则稀疏(非结构网格),块编码访问开销高于规则 stencil(对照 finite-difference-method);且 Hamiltonian 模拟要求厄米化处理(见 hamiltonian-simulation)。

延伸阅读