2. 流体力学基础回顾：连续性方程、动量方程（N-S方程）、湍流模型简介

各位同学，欢迎来到第二讲。做风电CFD仿真，说白了就是跟风打交道。但风这东西，看不见摸不着，想算准它，得先搞明白它背后的物理规律。

这一章，咱们不搞复杂的数学推导，我带你从工程应用的角度，把流体力学里最核心的三个家伙捋一遍：连续性方程、N-S方程，还有那几个常用的湍流模型。嗯，这些都是后续仿真的“内功心法”。

2.1 连续性方程：质量守恒的“记账本”

先问个问题：风从风机前面吹过去，质量会变少吗？

当然不会。这就是连续性方程的核心——质量守恒。你想想看，空气流过一段管道，入口进了多少质量，出口就得出去多少，除非中间有地方漏了或者存起来了。

对于咱们风电仿真里最常见的不可压缩流动（马赫数小于0.3，风速也就十几米每秒，完全满足），方程可以简化成这样：

∂u/∂x + ∂v/∂y + ∂w/∂z = 0

这个公式看着简单，但意义重大。它告诉我们：流入一个微元体的流量，必须等于流出的流量。我在做第一个风机尾流项目时，网格质量没把控好，结果算出来的流量不平衡，后处理一看，速度场全是乱的。后来排查了半天，就是连续性方程没满足，网格扭曲太严重了。

2.2 动量方程（N-S方程）：流体的“牛顿第二定律”

连续性方程只管“量”，不管“力”。那风为什么转？为什么有升力？这就得请出N-S方程了。

N-S方程，说白了就是F=ma在流体上的应用。它描述了流体微团的加速度与所受压力、粘性力、体积力之间的关系。对于不可压缩流体，它的矢量形式长这样：

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

别被这个公式吓到，咱们拆开看：

• 左边： 惯性项，代表流体的加速度。其中V·∇V是对流项，说白了就是“风吹着风跑”，这是N-S方程非线性的根源，也是最难算的部分。
• 右边第一项： 压力梯度力。风从高压区往低压区吹，就是这个力在驱动。
• 右边第二项： 粘性力。空气有粘性，会“拖拽”相邻的流体，造成能量耗散。
• 右边第三项： 体积力，比如重力。在风电里，除非研究大气边界层，一般可以忽略。

我个人习惯，在调试模型时，会先关掉湍流模型，用层流N-S方程算一个简单算例。如果连层流都算不对，那加上湍流只会更糟。这就像学走路，先站稳了再跑。

2.3 湍流模型简介：工程上的“妥协艺术”

现实中的风，几乎全是湍流。直接求解N-S方程（DNS）当然最准，但计算量大到无法用于工程。一台风机，雷诺数动辄上千万，DNS算一个工况，超级计算机都得跑几个月。所以，咱们得用湍流模型来“近似”。

湍流模型的核心思想，就是把N-S方程做时间平均（RANS），然后引入额外的方程来模拟湍流带来的“额外应力”——雷诺应力。这里我重点讲两个最常用的：

2.4.1 k-epsilon 模型：工业界的“老黄牛”

这个模型求解两个输运方程：湍动能k和湍流耗散率ε。它假设湍流是各向同性的，说白了就是认为湍流在各个方向上“乱”的程度差不多。

• 优点： 稳定、收敛快、计算量小。对于远离壁面的自由剪切流（比如风机尾流远场），效果不错。
• 缺点： 对近壁面流动、强分离流、旋转流预测不准。你想想，风机叶片附近流动分离很严重，用标准k-epsilon算，结果往往偏乐观。

我记得刚入行时，用k-epsilon算一个失速工况，结果升力系数一直偏高。后来前辈告诉我：“这模型在分离区就是个瞎子。”从那以后，我算叶片表面流动，再也没用过标准k-epsilon。

2.4.2 SST k-omega 模型：风电仿真的“主力选手”

SST k-omega是Menter提出的一个混合模型。它在近壁面用k-omega模型（对粘性底层解析好），在远场用k-epsilon模型（对自由流稳定），中间通过一个混合函数平滑过渡。

• 优点： 对逆压梯度、流动分离的预测能力远超k-epsilon。在风电领域，无论是叶片翼型分析，还是整机气动计算，SST k-omega都是首选。
• 缺点： 对入口边界条件中的omega值比较敏感，设置不当容易发散。

2.5 知识体系总览

说了这么多，咱们用一张图把这一章的核心逻辑串起来。你想想看，从物理规律到工程应用，路径其实很清晰：

这张图把咱们这一章的内容串起来了。从最底层的物理规律（连续性方程、N-S方程），到工程上不得不做的近似（湍流模型），再到具体选型（k-epsilon vs SST k-omega）。你以后做仿真，遇到问题就往这张图上靠，思路会清晰很多。