2. 流体力学控制方程

外流场CFD仿真的本质，是把空气绕过机翼、车身、翼型等物体时的流动，翻译成一组偏微分方程，再交给计算机离散求解。这组方程就是流体力学控制方程。理解它们的物理含义、守恒形式与非守恒形式的差别，以及外流场中哪些项可以合理简化，直接决定了你后面设置求解器、选择湍流模型和判断结果可信度的方式。

一、控制方程的物理来源：三个守恒律

控制方程不是凭空写出的数学公式，而是三条守恒律在连续介质上的表达：质量守恒、动量守恒、能量守恒。它们共同描述一个流体微团在运动过程中，质量、动量和能量如何随空间与时间变化。

连续性方程表达质量守恒。对可压缩非定常流动，其微分形式为：

∂ρ/∂t + ∇·(ρV) = 0

其中 ρ 为密度，V 为速度矢量。物理含义是：控制体内密度的局部变化率，等于流出控制体的净质量通量。对于不可压缩流动，ρ 为常数，方程退化为 ∇·V = 0，即速度场散度为零。外流场中，马赫数低于约 0.3 时通常按不可压缩处理，此时连续性方程只约束速度场，不再直接求解密度。

动量方程表达牛顿第二定律。以守恒形式写为：

∂(ρV)/∂t + ∇·(ρVV) = -∇p + ∇·τ + ρf

左端两项分别是动量密度的局部变化和对流通量；右端依次为压力梯度力、黏性应力散度和体积力。外流场中体积力通常只有重力，且常被忽略，因此核心是压力梯度与黏性应力的平衡。压力梯度驱动流动，黏性应力则通过边界层把物面影响传入流场。

能量方程表达热力学第一定律，守恒形式为：

∂(ρE)/∂t + ∇·(ρEV) = -∇·(pV) + ∇·(τ·V) + ∇·(k∇T) + S

其中 E 为单位质量总能量，k 为导热系数，T 为温度。右端分别代表压力做功、黏性做功、热传导和源项。外流场中，低速不可压缩流动的动能变化和黏性耗散对温度影响很小，能量方程常与动量方程解耦，甚至可先不求解；但在高速可压缩流动、气动加热或热防护问题中，能量方程必须与连续性、动量方程耦合求解。

二、Navier-Stokes方程：动量方程的完整形态

把牛顿黏性本构关系代入动量方程，就得到Navier-Stokes方程。对常黏度不可压缩流动，其向量形式为：

ρ(∂V/∂t + V·∇V) = -∇p + μ∇²V + ρf

这里 μ 为动力黏度。N-S方程是CFD求解的核心。它的非线性来自对流项 V·∇V，正是这一项使流动在雷诺数增大时出现分离、涡脱落和湍流。外流场仿真中，直接数值求解N-S方程（DNS）计算量极大，工程上通常采用雷诺平均（RANS）或大涡模拟（LES）对N-S方程做统计或滤波处理，再配合湍流模型封闭方程。

三、守恒形式与非守恒形式

同一个物理定律可以写成守恒形式和非守恒形式。守恒形式把时间导数写在通量散度之外，例如 ∂(ρV)/∂t + ∇·(ρVV)；非守恒形式则把密度移入对流项，例如 ρ(∂V/∂t + V·∇V)。两者在数学上可通过连续性方程互相转换，但在数值离散中差别明显。

守恒形式在有限体积法中天然对应“控制体通量平衡”，能保证质量、动量和能量在离散层面严格守恒，适合存在激波、强压缩或相变的问题。非守恒形式形式简洁，便于理论推导和量纲分析，但在跨激波等间断处可能给出错误结果。因此，外流场CFD求解器普遍采用守恒形式的积分方程，再对控制体表面通量进行离散。

四、外流场中的简化假设

完整N-S方程求解代价高，外流场工程仿真通常根据流动特征做分层简化。以下假设需要结合具体对象判断是否成立。

• 不可压缩假设：当马赫数 Ma < 0.3 时，密度变化小于约 5%，可令 ρ 为常数。此时连续性方程变为 ∇·V = 0，能量方程与动量方程解耦。例如汽车外流场、低速无人机机翼绕流，通常按不可压缩处理。
• 定常假设：若来流条件和物体姿态不随时间变化，且流动无大尺度非定常分离，可忽略 ∂/∂t 项，用定常RANS求解。但大迎角机翼、钝体尾流存在周期性涡脱落时，必须用非定常求解。
• 无黏与边界层近似：在高雷诺数外流场中，物面附近薄边界层内黏性主导，边界层外可视为无黏。工程上常用势流加边界层修正，或直接在全场求解RANS，由湍流模型提供等效黏性。
• 忽略体积力：除浮力驱动流动外，外流场中重力对气动力影响极小，动量方程中的 ρf 项通常略去。
• 量热完全气体假设：空气在常温至中温范围内，比热容近似为常数，状态方程用 p = ρRT 封闭。高温高速流动需考虑真实气体效应。

五、从方程到求解设置：一个外流场示例

以某低速翼型绕流为例，来流速度 30 m/s，弦长 1 m，空气密度 1.225 kg/m³，动力黏度 1.789×10⁻⁵ Pa·s。先估算：

Ma = V/a ≈ 30/340 ≈ 0.088  → 不可压
Re = ρVL/μ = 1.225×30×1 / 1.789e-5 ≈ 2.05×10⁶  → 湍流

据此可确定：采用不可压定常RANS，连续性方程用 ∇·V = 0，动量方程保留对流项和湍流等效黏性，能量方程暂不求解。湍流模型可选SST k-ω，以兼顾边界层内和自由剪切层。若后续把来流提高到 250 m/s，Ma ≈ 0.74，出现跨声速激波，则必须切换为可压缩求解，保留能量方程，并采用守恒形式的通量分裂格式。

实际设置时，重点检查三项：一是参考系与来流方向是否与方程中速度分量一致；二是压力-速度耦合算法是否与不可压假设匹配；三是离散格式是否对对流项采用足够低耗散的迎风格式，避免非物理振荡。求解收敛后，用残差下降、升阻力系数稳定性和网格无关性三重判据确认结果可信。

掌握控制方程的关键，不是背下所有公式，而是能根据流动的马赫数、雷诺数和分离特征，判断哪些项必须保留、哪些项可以简化，并让求解器设置与这些判断保持一致。这是外流场CFD仿真从“能算”走向“算得准”的第一步。