鼓风机模流仿真：进风口形状改1mm，叶轮出口速度均匀性差8%

风机设计有个特点：实验测试只能告诉你”行不行”——效率85%还是82%、噪声72dB还是68dB。但如果不行，实验不会告诉你是因为叶轮的入口角偏了2°，还是蜗壳的截面面积变化率太陡。

鼓风机模流仿真恰好能告诉你。数值风洞拆开看——每一个流道的静压恢复、每一圈叶片的做功分布、每一段蜗壳的动能-压力能转化率——这些量在台架上根本测不到。

一个做工业离心鼓风机的项目带着两张设计图来：一张是原始方案（叶轮外径320mm，12片后弯叶片，蜗壳矩形截面），一张是凭借经验改了几轮的新方案。台架测试显示新方案的效率从81%提到了84%，但全压从3200Pa掉到了2980Pa。

“压力掉下来了，看看流道里哪里堵了。”

旋转域和静止域的交接处理

鼓风机模流仿真第一个容易出错的点是旋转域（叶轮区）和静止域（进口段+蜗壳）之间的界面。Fluent和CFX都支持两种转子-定子界面类型：Frozen Rotor（冻结转子）和Transient Rotor-Stator（瞬态转子-定子）。Frozen Rotor假设界面两侧的流场在相对参考系中是”冻结”的，不做周向平均，计算量小，但只适用于界面附近流动在周向均匀的情况。

这个鼓风机的叶轮出口离蜗舌（cut-off）只有18mm——在叶轮旋转一圈的过程中，每个叶片通过蜗舌时都会产生一个压力脉冲。蜗舌区域的流动在周向高度非均匀，用Frozen Rotor会抹掉这个脉冲，把蜗舌处的静压抬高约5%。

改用了Transient Rotor-Stator，每个时间步内做界面上的相位对齐。叶轮每转2°（对应时间步长），界面网格做滑移。跑了完整5圈（1800个时间步），蜗舌处的静压波动从Frozen Rotor的”平坦直线”变成了±120Pa的正弦脉动——与台架实测蜗舌处壁面静压的脉动幅值（±135Pa）对上了。

进风口改的那1mm

新方案效率84%比原始方案81%高了3个百分点，但全压掉了220Pa。逐段拆开看：

叶轮入口到出口的静压升是2870Pa——比原始方案（3010Pa）低了140Pa。这个降幅对应流量从设计点（6200 m³/h）偏到了6880 m³/h——叶片出口角的后弯度（β₂=35°）让叶轮对背压敏感，压力一掉流量就上冲。

追到进风口——新方案为了降低进口损失，把进风口圆锥段缩短了15mm，进口直径从280mm扩大到290mm。扩大1厘米意味着进口面积增加了7%，但圆锥段的引导长度变短了，进口气流在进入叶轮前没有得到充分发展。叶轮入口截面的速度均匀性指数（Uniformity Index）从原方案的0.91降到了0.84——有接近10%的转速能量被入口分离区消耗了，没转成压力。

仿真在进风口做了三个几何迭代：原方案（D=280mm, L=120mm）、中间方案（D=285mm, L=115mm）、现方案（D=290mm, L=105mm）。D=285mm+L=115mm的进口速度均匀性恢复到0.89，全压回到3120Pa，只比设计值低2.5%——近权衡点。

蜗壳损失的两面性

蜗壳把叶轮出口的动能转化为静压恢复。转化率取决于蜗壳螺旋角的渐变设计——如果截面积随周角的增量过快，气流动能来不及转成压力能就被排出去了。

这个项目的蜗壳截面从90°到360°的扩张比是1:3.2——偏激进的设计，牺牲了一些静压恢复来换紧凑的几何尺寸。仿真能直接算出蜗壳各截面上的总压损失系数：从90°到180°是0.08（低损），从180°到360°跳到0.22——后半个蜗壳损失大了近三倍。问题在出口锥管的扩散角设了12°，过陡，附面层分离了。

调整扩散角到8°之后（保持出口法兰尺寸不变），全压在仿真里恢复了约110Pa。台架复测全压3125Pa——仿真的几何迭代真的把压力追回来了。

鼓风机模流仿真做对了，不是在Post里出一张漂亮的流线图，而是把整条压力-流量曲线拆成叶轮做功和蜗壳恢复两段来查账。