3、CFD基础理论（一）：流体力学控制方程（N-S方程）、湍流模型简介（k-epsilon, k-omega SST）
各位同学，欢迎来到CFD基础理论的第一讲。
说实话，很多做风资源评估的工程师，一听到「N-S方程」就头大。我当年刚入行时也一样，觉得这东西太数学了，跟实际工作没啥关系。直到有一次，我在一个复杂山地项目里，用工程模型算出来的结果跟实测差了20%以上……嗯，从那以后，我老老实实把CFD基础补了一遍。
今天我们就聊聊，风电场微观选址里最核心的流体力学控制方程，以及常用的湍流模型。我不打算堆公式，而是想让你明白：这些方程到底在算什么？我们为什么要选这个模型而不是那个？
核心观点：CFD仿真的本质，就是在计算机里求解流体运动的控制方程。你选的方程和模型，直接决定了仿真结果的精度和计算成本。
3.1 流体力学控制方程——N-S方程
先说说最基础的。流体运动遵循三个物理定律：质量守恒、动量守恒、能量守恒。在CFD里，我们把这三大定律写成数学形式，就是控制方程。
连续性方程（质量守恒）
说白了就是：流入控制体的质量，等于流出控制体的质量，加上内部质量的变化。对于不可压缩流体（比如低速空气），方程可以简化为：
∂u/∂x + ∂v/∂y + ∂w/∂z = 0这个公式的意思是：速度在三个方向上的变化率加起来等于零。我习惯把它理解成「流进来的和流出去的一样多」。在风电场仿真中，我们通常假设空气是不可压缩的，因为风速远低于音速。
动量方程（N-S方程）
这才是重头戏。N-S方程描述的是流体微团的动量变化，它考虑了压力、粘性力和惯性力的平衡。写成张量形式是这样的：
ρ(∂u_i/∂t + u_j ∂u_i/∂x_j) = -∂p/∂x_i + μ ∂²u_i/∂x_j∂x_j + f_i别被符号吓到。我来拆解一下：
- 左边第一项：时间变化率（非定常项）
- 左边第二项：对流项（惯性力）—— 这是最难算的部分
- 右边第一项：压力梯度
- 右边第二项：粘性扩散项
- 右边第三项：体积力（比如重力）
我在项目中遇到过一个问题：用稳态求解器算一个尾流场，结果怎么都不收敛。后来发现，是因为尾流区存在明显的周期性脱落涡，必须用非稳态（瞬态）求解才能捕捉到。所以，对流项的处理，是N-S方程求解的难点所在。
我的经验：对于风电场微观选址，90%的情况可以用稳态RANS（雷诺平均N-S方程）搞定。但如果你要研究尾流中的涡结构或动态载荷，就得考虑非稳态方法了。
3.2 湍流模型简介
直接求解N-S方程（DNS）需要极其精细的网格和巨大的计算资源，目前工程上根本做不到。所以我们用湍流模型来近似模拟湍流效应。
你想想看，风电场里的气流，经过风机叶片、地形起伏，到处都是湍流。如果不处理湍流，仿真结果基本没法看。
为什么需要湍流模型？
因为N-S方程本身是封闭的，但雷诺平均之后，多出了雷诺应力项，方程不封闭了。湍流模型的作用，就是提供额外的方程来封闭这个系统。
3.3 常用湍流模型：k-epsilon 与 k-omega SST
在风工程领域，最常用的两个模型就是k-epsilon和k-omega SST。我分别说说它们的特点和适用场景。
| 模型 | 核心思想 | 优点 | 缺点 | 适用场景 |
|---|---|---|---|---|
| 标准 k-epsilon | 求解湍动能k和耗散率ε | 计算稳定，收敛性好，应用广泛 | 近壁面处理差，对逆压梯度不敏感 | 自由剪切流、远离壁面的区域 |
| k-omega SST | 结合k-epsilon和k-omega的优点 | 近壁面精度高，能捕捉分离流 | 对入口边界条件敏感，计算量稍大 | 尾流区、复杂地形、分离流 |
标准 k-epsilon 模型
这个模型在工业界用得最多。它假设湍流是各向同性的，通过两个输运方程来求解湍动能k和耗散率ε。优点是鲁棒性好，收敛快。但有一个致命问题：在近壁面区域，它需要壁面函数来修正，否则误差很大。
我曾经在一个平坦地形项目中，用k-epsilon算出来的尾流恢复速度偏快，跟实测对不上。后来换成k-omega SST，结果就好多了。
k-omega SST 模型
这个模型是Menter在1994年提出的，它巧妙地结合了k-omega在近壁面的优势和k-epsilon在远场的优势。具体来说：
- 近壁面使用k-omega，能准确捕捉边界层
- 远场切换到k-epsilon，避免k-omega对自由流条件的敏感性
对于风电场尾流分析，我个人强烈推荐k-omega SST。为什么呢？因为尾流区存在明显的速度亏损和湍流混合，k-omega SST能更好地捕捉这些特征。
避坑指南：我曾经在一个山地项目中，直接用默认的k-epsilon模型，结果在背风坡出现了严重的分离流，仿真直接发散。后来改用k-omega SST，并调整了入口湍流强度，才得到合理的结果。记住：复杂地形一定要用k-omega SST。
3.4 知识体系框架
为了让你更直观地理解本章的知识结构，我画了一张图：
从这张图可以看得很清楚：我们从N-S方程出发，经过雷诺平均后，需要引入湍流模型来封闭方程组。而k-epsilon和k-omega SST，就是最常用的两个选择。
3.5 小结与建议
好了，这一讲的内容就到这里。总结几个要点：
- N-S方程是CFD的基石，理解它的物理意义比记住公式更重要
- 湍流模型是工程妥协的产物，没有完美的模型，只有合适的模型
- k-omega SST是风电场尾流分析的首选，尤其对于复杂地形和近尾流区
我的建议：刚开始做CFD仿真时，先用k-epsilon跑通流程，感受一下收敛性和计算时间。等你对结果有了一定判断力，再切换到k-omega SST做精细分析。别一上来就追求「最精确」的模型，先把基础打牢。
下一讲，我们会深入讨论网格生成和边界条件设置。到时候我会分享一些我在网格划分上踩过的坑，保证让你少走弯路。
公众号:蓝海资料掘金营，微信deep3321