1.2 CFD 核心理论基石
1.2 CFD 核心理论基石
工具地图铺好后，本节补上发动机原理：Fluent 求解器昼夜不停在解的到底是什么方程、方程里每一项是什么物理含义。这是全书的理论地基——第二章的有限体积法、第四章的湍流模型、第九章的收敛判读，全部从这节的三条守恒律与一条通用输运方程出发。读不懂这里，后面的操作就只能靠死记；读懂了，后面的每个设置项都能对号入座。
三条守恒律：流动世界的宪法
流体运动受三条守恒律约束。质量守恒：流进控制体的净质量流量等于控制体内质量的增长速率，数学上写成连续性方程。动量守恒：牛顿第二定律在流体微团上的表达，就是纳维-斯托克斯方程——压力梯度、黏性力与体积力共同决定流体的加速度。能量守恒：能量随对流输入、热传导扩散、黏性耗散与体积热源之间保持账面平衡。这三条律不看任何软件脸色，Fluent 只是它们的数值搬运工。
把动量方程写成工程手册里常见的形式（以 x 方向为例）：
∂(ρu)/∂t + ∂(ρuu)/∂x + ∂(ρvu)/∂y + ∂(ρwu)/∂z = -∂p/∂x + ∂τxx/∂x + ∂τxy/∂y + ∂τxz/∂z + ρfx
左边是对流输运（流体带着动量跑），右边依次是驱动项（压力）、耗散项（黏性应力）与外力项。真正的数值误差、发散、伪扩散，全部藏在这几项的离散方式里——这也是为什么第二章要花整章篇幅讲离散。
通用输运方程：一张表管住所有方程
CFD 教材最优雅的一步，是把质量、动量、能量乃至组分、湍动能的方程统一写成通用输运方程：瞬态项加对流项等于扩散项加固源项。四个槽位换上不同的物理量，就得到不同的方程。理解了这个结构，你在 Fluent 里勾选任何模型，本质上都是在往方程组里添加新的"φ 变量"。
图 1-2：通用输运方程四项解剖
四个槽位与常见方程的对应关系，值得动手抄一遍：
| 方程 | φ 的含义 | 扩散系数 Γ | 源项 S 的典型来源 |
|---|---|---|---|
| 连续性方程 | 1（退化为特殊形式） | 0 | 无 |
| 动量方程 | 速度分量 | 分子黏度 μ | 压力梯度、重力、旋转力 |
| 能量方程 | 温度 T 或焓 h | 导热系数除以比热 | 黏性耗散、辐射、反应热 |
| 组分方程 | 组分质量分数 | 组分扩散系数 | 化学反应生成率 |
| 湍动能方程 | 湍动能 k | 有效黏度除以系数 | 生成项与耗散项之差 |
连续介质假设与物性：方程成立的前提
纳维-斯托克斯方程把流体看成连续介质，这个假设用克努森数衡量——分子平均自由程与特征尺度的比值。常规工程流动里该比值远小于阈值，连续假设成立；但高真空设备或微尺度流动里它可能失效，那时就该换稀薄气体方法而不是硬调 Fluent。物性方面，密度如何随压力温度变化决定了可压缩性处理方式：马赫数低于约零点三时按不可压处理误差可忽略；接近或超过音速则必须启用可压缩流动与相应的耦合算法，这在第二章压力-速度耦合一节会展开。
一个估算案例：雷诺数决定建模路线
背景。新员工拿到任务：模拟水在内径五十毫米圆管内的流动，流量每秒十升，判断该用层流还是湍流模型。这是每个项目开工前的必答题。
操作。查物性：二十摄氏度水的密度约九百九十八千克每立方米，动力黏度约千分之一帕秒。算平均流速：流量除以截面积，约每秒五点一米。算雷诺数：密度乘速度乘管径再除以黏度，代入数字约二十五万。查判据：管内流动临界雷诺数约两千三百，当前值高出两个数量级，充分发展后为旺盛湍流。
结果与解读。建模路线定为湍流模型（第四章会具体到选哪一族），入口需要足够长的发展段或直接给定充分发展的速度型。假如算出来的雷诺数只有几百，路线则完全相反——层流求解，关闭湍流模型，网格近壁面要求也随之放松。同一段管子，数字一行，路线全变。
变式。非圆截面通道用当量直径（水力直径四倍截面积除以湿周）代替管径；变截面流动取最不利截面复核；非牛顿流体还要先确认黏度模型。估算这一步永远不会白做——它同时是验收单上"模型选择栏"的第一条证据。
常见判据速查与误区
开工前的三个量纲判断，值得贴在工位上。雷诺数：惯性对黏性之比，管流临界约两千三百，绕流临界约数十万，直接决定层流或湍流路线。马赫数：速度对当地音速之比，低于零点三按不可压处理，跨音速必须可压缩加激波捕捉。克努森数：分子自由程对特征尺度之比，超过约百分之一就要警惕连续假设。三个数算完，建模路线的骨架就定了，剩下的都是填肉。
| 常见误区 | 后果 | 正确姿势 |
|---|---|---|
| 不算雷诺数直接选默认湍流模型 | 层流被硬算成湍流，阻力传热全错 | 先量纲后选型，数字入档 |
| 高速流动仍按不可压设置 | 压力场与马赫波完全失真 | 马赫数超阈即转可压缩配置 |
| 物性取常数用到宽温域 | 温度结论漂移且难察觉 | 按温度范围决定分段或查表物性 |
| 把源项当开关随手加 | 守恒账本被无据源项污染 | 每个源项写明物理出处与量级 |
三个高频疑问。方程会解错吗：会，离散误差、迭代误差与模型误差三层叠加，第九章的三重签字与 GCI 就是为拆这层账准备的。需要背下 N-S 方程吗：需要的是结构而不是默写——对流、扩散、源项的三段式账本结构记牢，软件里每个设置都能归位。湍流为什么不是方程的问题而是封闭问题：方程仍然成立，只是平均处理后多出不认识的雷诺应力项，模型的全部性格都来自对这项的封堵假设，第四章由此展开。
从方程到软件：翻译链条的全景
把这一节的知识串成一条链：守恒律给出方程组，通用输运方程给出统一结构，无量纲判据给出建模路线，物性给出介质身份——Fluent 里每个下拉选项都是这条链上的一环。看懂链条的好处是排错有方向：结果可疑时，先问是方程层级的问题（模型选错）、离散层级的问题（格式与网格不够格）、还是输入层级的问题（边界与物性设错），三层的检查动作完全不同。第一章之后全书就沿着这个层级展开：第二章与第三章治离散层，第四章到第八章治方程层，第九章把三层证据合账。你现在读的是地图的总纲，之后每一章都会回到这条链条上加盖一层理解。
本节验收点
- 能默写通用输运方程的四个槽位，并给任一方程对号入座。
- 开工前已算过雷诺数（或马赫数、克努森数等无量纲判据），结论写进建模说明。
- 连续介质假设在当前工况下成立，有数字佐证。
- 物性来源与温度范围已记录，物性取常数还是随温度变化已明确。