台风来了别只刷预警，用 Fluent 把风场自己算一遍
台风"白海豚"那几天，朋友圈一半人在转气象台路径预报，一半人在吐槽"又跑偏了"。
我这个做 CFD 的，看这事角度不太一样：你们看到的是一张预报图，我看到的是——这玩意儿本质就是一次大气尺度的流体仿真。
气象台背后跑的 WRF、GRAPES 那些模式，跟咱们工程上用的 Fluent，底层是同一套方程：Navier-Stokes。区别只在尺度，人家算几千公里，咱们算几米到几公里。
所以今天聊点实在的：手上要有 Fluent，能不能自己把一场台风过境时的城市风场、地形绕流算出来？
能。而且步骤比你想的规矩。下面这套 Fluent 的边界条件、湍流模型、UDF 怎么设，我的 ima 知识库《流体仿真》里都有整章原文，本文就是照着那些资料捋出来的实操。
图 1 · 用 Fluent 模拟气象的完整工作流（六步）
一、先搭个"风洞"——计算域和几何
气象模拟的计算域简单粗暴：一个大方块。要算城市风场，就把楼群的几何体塞进去；要算地形绕流（比如台风翻山），就把山体 STL 放进去。
知识库里《赵玉新 FLUENT 教程》第一章就讲：先想清楚要计算的物理域，再谈网格。域太小说不清边界影响，太大又浪费算力。常规做法是上游留 5–10 倍特征长度，下游留 10–15 倍，顶部伸到边界层高度以上。
图 2 · 计算域与边界条件布置（入口挂风廓线，地面为粗糙壁面）
二、网格：近地面一定加密
大气边界层（ABL）的麻烦在于——风速在离地几十米内变化最猛（对数律），所以贴地那层网格必须密。
知识库《FLUENT 边界条件》章讲壁面时专门提了：模拟大气边界层、机翼、船体这些近壁湍流，壁面函数对第一层网格高度 y⁺ 很敏感。地面那层用边界层网格（inflation），增长率别超过 1.2。
这一步直接决定你算出来的风剖面准不准。
图 3 · 近地边界层网格加密（第一层高度 y⁺ 敏感）
三、边界条件：这是气象模拟的灵魂
Fluent 边界条件类型一堆，做风场真正用得上的就四种（知识库《边界条件》章原文就是这么分类的）：
入口：速度入口（velocity-inlet）。气象里风速不是常数，是随高度变化的廓线，所以入口要挂一个风廓线——要么用 profile 文件，要么写个 UDF。指数律 v(z)=v_ref·(z/z_ref)^α，或者更准的对数律 v(z)=(u*/κ)·ln(z/z₀)。
出口：压力出口（pressure-outlet）。出现回流别慌，给个合理回流条件就行。
顶部、远场：对称边界（symmetry）或压力远场（pressure-far-field，仅可压流用）。稳态风场常用对称。
地面：壁面（wall），而且一定要设壁面粗糙度 K_s。知识库讲得明白：大气边界层就是典型的粗糙壁面——"流过机翼、船体、涡轮机、换热器以及管系统的流动，还有具有各种粗糙度的地面上的大气边界层"。光滑壁面算出来的风剖面是错的。
图 4 · 大气边界层风速廓线：指数律 vs 对数律
四、物理模型：开能量方程，别忽略浮力
风场模拟开不开能量方程？看你关不关心温度。台风这种大尺度，温差驱动的密度变化（浮力）是主因。
知识库《基本物理模型》章"浮力驱动流动和自然对流"一节写得透：密度随温度变 → 重力下产生浮力 → 自然对流。可以用 Boussinesq 近似（小温差、收敛快），也可以直接让密度等于温度函数（理想气体）。
湍流模型选 Realizable k-ε 或 RNG k-ε，配壁面函数，基本够用。入口湍流参数怎么给？知识库给了现成公式：完全发展管流核心湍流强度 I≈0.16·Re_Dh^(−1/8)；壁面限制流动用"湍流强度+长度尺度"，l=0.07L（或 l=0.4δ_99）。外部流自由流湍流强度能低到 0.05%，城市冠层里能到 10% 以上——这个量级千万别给错。
五、求解：SIMPLE 走起
压力-速度耦合用 SIMPLE 或 SIMPLEC，二阶迎风格式，亚松弛因子先跑默认，不收敛再调。
想看台风眼的涡旋结构、楼间角隅涡，就把工况设成非定常（transient），时间步按涡脱落周期取 20 步左右。知识库"时间相关模拟"章讲得很细。
六、后处理：把风场"看"出来
算完不是看一堆数字。打开云图：
· 风速云图——看哪里风速被放大、哪里是死区；
· 流线 / 矢量图——看涡在哪儿、回流在哪儿；
· 要算污染物扩散，就开组分输运，看浓度云图。
这些在后处理面板里点几下就出来。
图 5 · 结果云图：风速分布与流线
六步跑下来，你会发现"用 Fluent 算天气"不是玄学，就是一套规范流程：几何 → 网格 → 边界条件 → 模型 → 求解 → 后处理。跟算一个汽车外流、一个换热器，骨架一模一样，区别只在边界条件和物理模型怎么贴你的气象场景。
