有限体积法

有限体积法（finite volume method, FVM）把守恒律直接在控制体上积分，让"物理上守恒"变成"代数上守恒"：相邻单元共享同一界面通量，进出相消，激波以正确速度传播。它是可压缩流与工业 CFD 的绝对主力——Fluent、STAR-CCM+、SU2 到航空外流求解器几乎清一色基于 FVM。

守恒离散：从积分形式出发

可压缩流写成守恒律形式（见 navier-stokes）：

$$\frac{\partial \mathbf{U}}{\partial t}+\nabla \cdot \mathbf{F}(\mathbf{U})=0,$$

$\mathbf{U}=(\rho ,\rho \mathbf{u},\rho E{)}^{\mathsf{T}}$ 是守恒变量向量。对单元 $V_p$ 积分并用高斯定理：

$$\frac{\mathrm{d}}{\mathrm{d}t}{\int }_{V_p}\mathbf{U}\mathrm{d}V=-{\oint }_{\partial V_p}\mathbf{F}\cdot \mathbf{n}\mathrm{d}S\approx -\sum _{e\in \partial V_p}{\mathbf{F}}_e\cdot {\mathbf{n}}_e{A}_e,$$

即单元平均值 ${\bar{\mathbf{U}}}_p$ 的演化由界面通量 ${\mathbf{F}}_e$ 决定。对两个相邻单元，同一界面取同一个 ${\mathbf{F}}_e$——通量守恒在离散层面精确成立，总守恒量按构造守恒（求和时内界面通量两两相消）。这是 FVM 与有限差分的本质区别：后者只在网格节点近似微分算子，守恒是误差阶的结果而非恒等式。

通量重构与 Riemann 求解器

界面通量怎么算？两侧单元平均值给出左、右状态 $({\mathbf{U}}_L,{\mathbf{U}}_R)$。Godunov 的洞见：把它们看成局部Riemann 问题（左右两常态、界面间断的精确解），用其界面处自相似解的通量作为 ${\mathbf{F}}_e$：

$${\mathbf{F}}_e=\mathbf{F}({\mathbf{U}}_{\mathrm{RIEMANN}}({\mathbf{U}}_L,{\mathbf{U}}_R;x/t=0)).$$

精确 Riemann 求解器对 Euler 方程需要迭代，实用中用近似格式：Roe（线性化 $\partial \mathbf{F}/\partial \mathbf{U}$ 在平均态取值、解线性化 Riemann 问题）、HLL/HLLC（两/三波速框架估计通量，HLLC 恢复接触间断）、Lax-Friedrichs/Rusanov（最稳健、最耗散）。高阶精度通过空间重构（MUSCL：由单元平均重构界面两侧的线性分布）与时间推进（TVD Runge-Kutta）获得；为抑制激波附近的振荡，重构处加限制器（minmod、MC 等）把解限制在单调区间。

对不可压 FVM，压力-速度耦合仍需特殊处理（Rhie-Chow 动量插值防棋盘振荡、SIMPLE/PISO 压力修正族，见 projection-method）。

为什么工业选择 FVM

• 守恒性对激波、燃烧、多相流是硬需求；
• 天然适配非结构网格：任意多面体单元，复杂几何（全机、发动机）无需坐标变换；
• 粘性/无粘、层流/湍流（RANS 见 turbulence-simulation）在同一框架下装配；
• 代价：非结构稀疏矩阵结构不规则，线性求解器与并行通信模式比结构网格复杂（见 multigrid）。

与量子算法的接口

FVM 离散后同样是 $A\mathbf{u}=\mathbf{b}$，量子提案可以直接挂上。代表性工作：本源量子（Origin Quantum）的混合量子-经典 CFD 框架用有限体积离散、以变分量子线性求解器（VQLS）解线性系统，在 Poiseuille 与方腔等基准上完成数值验证Ye et al. 2024 (Phys. Fluids)——工程化讨论见 hybrid-hpc-integration。FVM 的非结构矩阵对块编码（见 block-encoding）不如有限差分的规则模板友好，这是量子 FVM 需要面对的结构性差异；各家离散化选择对照见 quantum-navier-stokes。

延伸阅读

• finite-difference-method
• projection-method
• vqls
• quantum-navier-stokes
• hybrid-hpc-integration