有限体积法：守恒离散的核心思路
理解控制体积分、面通量和离散系数。
入门
本文要点
建立物理量、守恒形式和离散方程之间的联系；识别模型假设；用量纲与极限工况检查推导结果。
1. 问题与物理背景
计算流体力学的核心不是把连续方程直接交给计算机，而是把一个物理问题转化为有限、守恒且可以验证的代数系统。任何可靠的求解都应明确研究对象、控制体、时间尺度与边界条件。
对不可压缩牛顿流体，质量守恒要求速度场满足零散度条件。动量守恒则刻画惯性、压力、黏性和体力之间的平衡：
∂u/∂t + (u · ∇)u = −(1/ρ)∇p + ν∇²u + f(1)
2. 守恒形式与控制体
积分守恒形式是有限体积离散的自然起点。通过高斯定理，体积分中的散度项被转换为控制体各表面的通量之和，从而使相邻控制体共享同一面通量。
工程提示
守恒不等于精确。严格守恒的离散格式仍可能有较大的耗散或色散误差，因此必须同时检查精度、稳定性和网格收敛性。
| 方程项 | 物理含义 | 数值关注点 |
|---|---|---|
| ∂u/∂t | 局部加速度 | 时间步长与阶数 |
| u · ∇u | 对流惯性 | 有界性与数值耗散 |
| ∇p | 压力驱动 | 压力—速度耦合 |
| ν∇²u | 黏性扩散 | 网格正交性 |
3. 从连续方程到离散系统
空间离散决定各控制体之间如何交换质量、动量与能量。时间离散决定不同时间层之间的信息传播。离散后通常得到稀疏线性或非线性代数系统。
while residual > tolerance:\n solve_momentum()\n solve_pressure_correction()\n correct_fluxes()4. 验证与可信度
计算结束不意味着结论成立。至少需要检查残差、全局守恒、网格敏感性以及与解析解、实验数据或公开基准的偏差。
- 先确认输入量纲与边界条件。
- 观察残差是否稳定下降，避免只看单一方程。
- 核对质量或能量的全局不平衡。
- 使用至少三组网格评估离散误差。
参考与修订
[1] Ferziger, Perić & Street, Computational Methods for Fluid Dynamics.
[2] Versteeg & Malalasekera, An Introduction to Computational Fluid Dynamics.