4、CFD基础理论（二）：边界层理论、地形与粗糙度对风廓线的影响、大气稳定度
各位好，咱们接着聊CFD基础理论。上一节我们把流体力学的基本方程过了一遍，今天要聊的这三个概念——边界层、地形粗糙度、大气稳定度——说实话，是风电场微观选址里最容易出问题的地方。
我个人习惯是，拿到一个新场址，先不看风玫瑰，先看这三个东西。为什么？因为它们决定了风到底怎么吹到你的风机面前。你想想看，同样的风速，在平原和在山脊上，能量差了一倍都不止。
4.1 大气边界层：风从哪里来，到哪里去
大气边界层，说白了就是地球表面往上大概几百米到一公里那层空气。这层空气直接受地面摩擦影响，风速从地面为零，慢慢增加到边界层顶部的“地转风”速度。
我记得刚入行那会儿，有个老工程师跟我说：“小张，别把大气边界层当成均匀的，它分三层呢。”确实如此：
- 近地层（0-100m）：风切变最剧烈，湍流最强。咱们风机正好在这层里转。
- 埃克曼层（100m-1km）：风向开始偏转，科里奥利力起作用了。
- 自由大气（1km以上）：地面摩擦基本消失，风就是地转风。
做CFD仿真时，边界层厚度怎么取？我一般用这个经验公式：
h = 0.1 * u* / f其中u*是摩擦速度，f是科里奥利参数。算出来大概300-600米。但要注意，山区可能只有100米，平原能到800米。别死套公式。
核心要点：边界层厚度决定了你的计算域高度。设得太矮，风被“压扁”了；设得太高，浪费网格。我一般取风机轮毂高度的10-15倍。
4.2 风廓线：指数律 vs 对数律
风廓线就是风速随高度变化的曲线。做微观选址，你总得知道轮毂高度处的风速吧？但测风塔往往只有10m、50m、80m几个高度，怎么推？
两种主流方法：
| 方法 | 公式 | 适用场景 |
|---|---|---|
| 指数律 | V(z) = V_ref * (z/z_ref)^α | 工程快速估算，α取0.1-0.4 |
| 对数律 | V(z) = (u*/κ) * ln(z/z0) | 科研、CFD入口条件，精度更高 |
指数律简单，但有个坑：α不是常数。我在内蒙古一个项目里，白天α=0.12，晚上能到0.35。为什么？因为夜间稳定层结，风切变更大。你如果只用白天数据推全年发电量，误差能到15%。
对数律更物理，但需要知道粗糙度长度z0。z0怎么取？看下表：
| 地表类型 | z0 (m) |
|---|---|
| 水面/冰面 | 0.0001 - 0.001 |
| 草地 | 0.01 - 0.05 |
| 农田 | 0.05 - 0.15 |
| 森林 | 0.5 - 1.0 |
| 城市 | 1.0 - 3.0 |
我的经验：做CFD时，入口风廓线最好用对数律。指数律虽然简单，但会导致近地面风速梯度失真，影响尾流计算精度。我曾经对比过，用指数律算出的尾流恢复长度比实测短了20%。
4.3 地形与粗糙度：CFD里怎么处理
地形对风的影响，说白了就是“加速”和“减速”。山脊顶部风速能比平地高30%，山谷里可能只有一半。
CFD里处理地形，我建议用贴体网格。什么意思？就是网格跟着地形走，而不是用阶梯近似。阶梯近似在陡峭地形上误差很大，尤其是计算分离流时。
粗糙度怎么设？两种思路：
- 均匀粗糙度：整个计算域用一个z0值。适合平坦均匀地形。
- 非均匀粗糙度：每个网格单元单独设z0。适合复杂地形，比如一半农田一半森林。
我个人强烈建议用非均匀粗糙度。为什么？因为粗糙度突变会导致内边界层发展，风速和湍流都会剧烈变化。你想想看，风从草地吹到森林，风速骤降，湍流增强，这对风机载荷影响很大。
避坑指南：我曾经在一个项目里，把森林区域的z0设成了0.1（应该是0.5），结果算出来的发电量比实际高了8%。后来才发现是粗糙度设错了。记住：粗糙度每差一个数量级，轮毂高度风速能差0.5m/s。
4.4 大气稳定度：被忽视的关键因素
大气稳定度，说白了就是大气“想不想”上下混合。它分三种：
- 不稳定：地面热，空气往上跑。湍流强，风切变小。
- 中性：温度随高度递减刚好等于干绝热递减率（约1℃/100m）。湍流中等。
- 稳定：地面冷，空气“赖”在地上不动。湍流弱，风切变大。
为什么这个对风电场重要？因为稳定层结下，尾流恢复特别慢。我做过一个对比：中性条件下，尾流在下游10D处恢复80%；稳定条件下，20D才恢复80%。这意味着什么？意味着稳定层结多的地区，风机间距要拉大。
CFD里怎么考虑稳定度？用Monin-Obukhov相似理论。核心参数是莫宁-奥布霍夫长度L：
L = - (u*^3 * θ_v) / (κ * g * (w'θ_v'))其中θ_v是虚位温，w'θ_v'是湍流热通量。L的绝对值越小，稳定度影响越大。
实际工程中，我一般这么处理：
- 如果L < -100m：不稳定，用对流边界层模型
- 如果|L| > 500m：中性，用标准对数律
- 如果L > 100m：稳定，用稳定边界层模型
重要提醒：很多CFD软件默认用中性层结。但实际风电场，尤其是内陆地区，夜间稳定层结占比可能超过50%。忽略稳定度，你的尾流分析会严重偏乐观。
4.5 知识体系总览
下面这张图，是我自己总结的边界层与风廓线核心逻辑。你看一遍，应该能串起来：
嗯，这张图把今天讲的内容串起来了。边界层是基础，风廓线是工具，地形、粗糙度、稳定度是修正项。做CFD时，这三块缺一不可。
最后说一句：理论归理论，实际做项目时，我建议先用测风塔数据校准你的风廓线模型。别信任何默认参数，每个场址都有自己的脾气。你花一天时间调参数，可能省下后面一个月的返工。
公众号:蓝海资料掘金营，微信deep3321