第2章流体力学基础

2、流体力学基础：流体属性、控制方程（连续性、动量、能量）、湍流模型简介

各位好，欢迎来到流体力学基础这一节。说实话，很多做结构出身的朋友一听到流体就头大。我当年刚接触流固耦合时也一样，觉得流体力学公式又多又抽象。但后来我发现，搞懂几个核心概念，剩下的就是查表和调参数了。咱们今天就把这些基础捋一遍。

做仿真前，第一件事就是搞清楚你面对的是什么流体。空气和水，虽然都是流体，但行为差远了。我个人习惯先把这几个属性列出来：

我的小经验： 在设置材料属性时，别直接用默认值。我遇到过好几次因为忘记修改温度相关的粘度参数，导致仿真结果偏差很大的情况。尤其是高风速工况，空气温度变化会显著影响粘度。

流体运动遵循三大守恒定律：质量守恒、动量守恒、能量守恒。用数学语言表达出来，就是下面这三个方程。你不需要背下来，但得理解它们各自在说什么。

说白了就是：流入控制体的质量，等于流出控制体的质量，加上内部质量的变化。对于不可压缩流体（比如低速空气），方程可以简化为：

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

如果密度是常数（不可压缩），就变成：

∇·u = 0

嗯，这里要注意：风机叶片周围的空气，在低速时通常当作不可压缩处理。但如果你仿真的是高速叶片（比如大型海上风机），就得考虑可压缩性了。

这个方程描述的是：流体的加速度等于它受到的各种力（压力、粘性力、体积力）之和。写成张量形式就是：

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

左边是惯性项，右边依次是压力项、粘性项、体积力项。这个方程是非线性的，很难直接求解。我刚开始做仿真时，总想着能不能手算验证一下，后来发现根本不可能——只能靠数值方法。

  关键点： N-S方程是流体仿真的核心。所有湍流模型、边界层处理，本质上都是在想办法近似求解这个方程。

如果你关心温度场（比如叶片结冰分析、热应力分析），就需要能量方程。它描述的是：流体能量的变化率，等于热传导、粘性耗散、对外做功的总和。

ρc_p(∂T/∂t + u·∇T) = k∇²T + Φ

其中Φ是耗散函数，代表粘性摩擦产生的热量。在大多数风机叶片流固耦合中，如果不涉及温度变化，能量方程可以暂时忽略。但我建议你至少知道它的存在，万一哪天需要分析热效应呢？

你想想看，真实的风机叶片周围，气流是乱七八糟的——有漩涡、有分离、有脉动。直接求解N-S方程（DNS）需要网格小到微米级，计算量大到离谱。所以我们需要湍流模型，用平均量来近似描述湍流效应。

常用的湍流模型有这几类：

避坑指南： 我曾经在一个叶片失速工况仿真中，用了标准的k-ε模型，结果分离点预测完全不对。后来换成k-ω SST模型，结果才合理。所以我的建议是：对于风机叶片这种有分离、有逆压梯度的流动，优先选k-ω SST。

下面这张图，是我自己总结的流体力学基础在风机叶片仿真中的位置。你可以把它当作一个导航图：

从这张图可以看出来，流体属性、控制方程、湍流模型是三大基础。它们决定了边界条件怎么设、网格怎么画、求解器怎么调。最终所有工作都指向流固耦合仿真这个目标。

我的建议： 刚开始学的时候，别急着把所有模型都搞懂。先把不可压缩流动、k-ω SST模型、稳态求解这三个组合练熟。我当年就是这么过来的，等基础打牢了，再慢慢扩展。

好了，这一节的内容就到这里。流体力学基础虽然看起来理论性强，但只要你动手做几个案例，很快就能找到感觉。记住：仿真不是数学考试，理解物理意义比背公式重要得多。

公众号:蓝海资料掘金营，微信deep3321