4. CFD微尺度模型原理：RANS与LES对比，地形网格划分策略，边界层参数化

各位同行，咱们今天聊点硬核的。CFD微尺度模型，说白了就是给风电场做「CT扫描」。你想想看，宏观的天气图只能告诉你「这片山风大」，但具体到每个机位点，风是怎么绕山的、湍流有多强，就得靠CFD来精细刻画了。

我个人习惯把CFD模型分成两类：RANS和LES。这俩就像「快照」和「录像」的区别。咱们一个个说。

4.1 RANS与LES：两种核心思路

RANS（雷诺平均Navier-Stokes），我刚开始做风资源评估时，用的就是它。它的核心思想很简单——把湍流运动分解成「平均流」和「脉动」，然后只求解平均流场，脉动部分用模型来参数化。

我在云南一个项目里遇到过，用RANS算出来的尾流区风速，跟实测差了将近15%。后来换成LES才把问题找出来——原来那个山脊后面有个巨大的回流区，RANS根本没捕捉到。

LES（大涡模拟），说白了就是「只过滤小涡，直接算大涡」。它把湍流分成「可解尺度」和「亚格子尺度」，大尺度的涡直接求解，小尺度的用亚格子模型近似。

嗯，这里要注意。LES不是万能的。我曾经在内蒙古一个平坦地形项目里，用LES和RANS对比，结果两者差异不到3%。所以我的建议是：复杂山地、强湍流区域，优先考虑LES；平坦地形、简单流场，RANS完全够用。

对比项	RANS	LES
计算成本	低（小时级）	高（天级）
湍流解析度	完全参数化	大涡直接求解
适用场景	平坦/缓坡地形	复杂山地/强分离流
工程常用软件	WindSim, Meteodyn WT	OpenFOAM, SOWFA

4.2 地形网格划分策略

网格划分，这是CFD里最「玄学」的部分。我见过太多人，模型参数设得再准，网格一塌糊涂，结果全白搭。

核心原则：网格要「疏密有致」。地形变化剧烈的地方加密，平坦区域可以粗一些。

我曾经在贵州一个项目里，网格划分太粗，结果算出来的风速比实测低了20%。后来把机位点附近的网格从50m加密到15m，误差直接降到5%以内。你想想看，网格这东西，真不能省。

另外，网格质量检查也很关键。我建议用两个指标：

• 正交性：网格夹角尽量接近90°，小于30°的网格要避免
• 长宽比：一般控制在5:1以内，边界层内可以放宽到10:1

4.3 边界层参数化

边界层参数化，说白了就是「近地面风怎么算」。大气边界层的高度一般在几百米到一两公里，而我们的风机轮毂高度通常在80-150m，正好落在边界层底层。

常用的边界层模型：

• 对数律：最简单，适用于中性层结、平坦地形。公式是 u(z) = (u*/κ) * ln(z/z0)
• 指数律：工程上常用，但精度一般。指数α通常取0.1-0.3
• Monin-Obukhov相似理论：考虑了大气稳定度，精度更高，但需要输入感热通量等参数

我个人习惯用Monin-Obukhov理论，尤其是在复杂山地项目里。为什么？因为山地的热力效应很明显，白天山坡受热上升气流强，晚上冷空气下沉，这些都会影响风切变。对数律在这种场景下，误差能到10%以上。

嗯，这里有个坑。很多人在设置入口边界时，直接用对数律算一个风廓线就完事了。但实际项目中，上游地形的影响会让风廓线严重偏离理论值。我建议有条件的话，用测风塔数据或激光雷达实测廓线来标定入口条件。

4.4 知识体系总览

下面这张图，是我自己总结的CFD微尺度模型核心逻辑。你把它理清楚了，做微观选址心里就有底了。

这张图你看懂了吗？从上到下，就是咱们做CFD微观选址的标准流程。先根据地形复杂度选模型，再精细划分网格，最后用合理的边界层参数收尾。每一步都不能马虎。