2.1 有限体积法深度解析

2.1 有限体积法深度解析

理论基石就位，本节开始拆数值引擎的第一根柱子：Fluent 把连续的流场切成什么样、在每个小块上做什么样的物理账本。它是理解后面两节耦合与格式问题的前提——所有"结果不守恒"的悬案，最终都要回到控制体的通量平衡上来查账。

把流域切成控制体：守恒从哪里来

有限体积法（FVM）把计算域分割成互不重叠的控制体，对每个控制体积分守恒方程：流入减流出等于积累，再加源项。它的杀手锏是面通量的配对消去——相邻两个控制体共享同一个面，你从这个面流出的量就是我流入的量，全域求和时内部面的通量两两抵消，只剩边界通量。这意味着只要每个面的通量计算一致，离散格式在任意网格上都严格守恒，哪怕网格歪歪扭扭。这正是 CFD 选择 FVM 的根本原因：有限元在贴体坐标下守恒要额外设计，有限差分在非结构网格上几乎无法保证通量配对，而工程流体恰恰最在乎"质量能量账本对不对得上"。

Fluent 里每个网格单元（cell）就是一个控制体，网格面（face）是相邻控制体的公共边界，边界条件施加在边界面（boundary face）上。这三个词是第三章网格验收的主角，此处先立此存照。

图 2-1：控制体通量平衡示意

从积分到代数方程：面值重构与通量计算

对控制体积分后，时间项用单元平均值表示，麻烦全部集中在面上：方程里的对流项与扩散项都需要"面上的值"，而求解器手里只有单元中心的值。于是产生两个动作。其一是面值重构——用单元中心值插值出面值，插值方式决定精度与稳定性，这是 2.3 节的主角。其二是通量计算——对流通量等于面值乘质量流量，扩散通量等于面法向梯度乘扩散系数，梯度本身还要用高斯定理从单元值反推。

把这些代数关系拼起来，每个控制体得到一个只含邻居单元值的方程，全域构成大型稀疏线性方程组。Fluent 用点隐式高斯-塞德尔加代数多重网格（AMG）求解这个方程组：先在细网格上松弛，残差转移到粗网格快速消解，再插值回细网格。你不需要实现它，但需要知道两件事：矩阵求解的代价与网格量近线性，这就是第十一章并行扩展的物理基础；多重网格只加速矩阵求解，不改变物理模型。

手算一维对流扩散：控制体方法的第一课

背景。一根一维管道，网格只有三个单元，解一个既有对流又有扩散的标量（比如温度），边界左端固定值为百、右端固定值为零。目标：亲手组装 FVM 方程，体会"面值怎么取"这个细小决定的分量。

操作。取均匀网格、间距为一。对中间单元写通量平衡：左面流入减右面流出加源为零。对流项的面值直接取上游单元值（一阶迎风），扩散项的面值取两侧单元值的线性插值。整理得到中间单元的方程：左系数、右系数与中心系数构成一行三对角矩阵。三个单元依次写出，就得到一个三行方程组，手算即可求解。

结果。解出的温度分布明显偏向高值边界——把浓度剖面拉得比精确解更"钝"。这不是算错，恰恰是一阶迎风引入的数值扩散在起作用：它等效于给扩散系数偷偷加了料，加料多少与网格佩克莱数成正比。

解读。这个小实验复现了工程 CFD 最常见的慢性病：结果看着合理、剖面系统性偏肥、加密网格才逐渐好转。验收单上"网格无关性核对"一栏，防的就是它。判断标准也在这里埋下伏笔：网格佩克莱数大于二时一阶迎风的数值扩散不可忽略，工程上要么加密网格要么换高阶格式。

变式。把面值改成一阶迎风与中心差分的加权（Fluent 里的二阶迎风其实是沿特征线方向的二阶重构），数值扩散立刻显著减轻，但高佩克莱数下面值可能越过两侧单元值的合理区间，产生新振荡——精度与稳定的交易在 2.3 节正式展开。

Fluent 中的对应设置

有限体积法的理论在软件里浓缩为少量选项：Methods 面板的梯度方案（Least Squares Cell Based 默认、Green Gauss 节点型或单元型）决定梯度重构方式；方程离散与压力插值选项决定面值行为；Residual 定义里的归一化方式决定你看到的残差数字尺度。初学阶段保持默认梯度方案即可，但对强各向异性网格（比如拉伸严重的边界层网格），Least Squares 通常更稳。这些开关何时动、动了怎么验证，按验收单纪律，每动一项都要在报告里留一行说明。

梯度方案与网格量的实用对照

有限体积法在软件里的第一排开关是梯度重构，它与网格类型有稳定的搭配关系。Least Squares Cell Based 默认且在多边形网格上表现均衡；Green Gauss Node Based 在三角形为主的网格上梯度更准（曾经是燃烧类推荐的组合）；Green Gauss Cell Based 省内存，精度要求不苛刻的初算可用。判断原则：网格越歪、各向异性越强，越值得把梯度方案当回事；结构化六面体加温和拉伸时，各方案差距通常淹没在模型误差里。

三个高频疑问。守恒格式还会算出质量不守恒吗：会，网格质量、线性求解不完全收敛、源项离散都可能破坏账面，所以第九章的通量平衡检查永远要做。单元中心格式与节点格式怎么选：常规任务单元中心足够，需要精细梯度（后处理应力、精细传热）时节点格式值得多付的内存。手算案例的结论能外推到三维吗：结论的方向（一阶格式数值扩散）完全外推，量级随维度与网格变化——所以工程上靠 GCI 而不是靠手算定误差带。

还有一个值得记牢的工程视角：有限体积法的守恒性是把双刃剑。它保证了离散账面的自洽，但账面守恒不等于解正确——一个数值扩散严重的解同样可以完美守恒。所以工程报告里"质量守恒"与"结果正确"是两个独立命题，前者靠通量平衡签字，后者靠网格无关性与对表签字。把这对概念分开的人，才不会被"账平了所以没问题"的说辞带偏。守恒是必要条件，这句话值得写在每个模板的开头。

本节验收点

• 能向同事解释守恒性来自面通量配对，与网格质量解耦。
• 知道面值重构与通量计算是数值误差的两大源头。
• 手算过一维对流扩散，亲眼见过一阶迎风的数值扩散。
• 网格佩克莱数概念已建立，能在报告里用它为格式选择辩护。