用ANSYS CFX对离心泵内部流动做瞬态仿真的方法与结果分析

使用 ANSYS CFX 对离心泵内的流动进行瞬态仿真。湍流模型采用 SST。同时包含 CFX 定义文件。Centrifugal-Pump-Simulation.def

RMS residual is above 1e-5 after 100 iterations 弹出来的时候，心里第一反应是稳态算得挺好的怎么瞬态就跑不动了。稳态残差降到1e-5作为瞬态初始场，CFX里瞬态计算的单步收敛标准可以放到1e-4，先把流动建立起来再收紧。

CFX里离心泵瞬态计算分两个域，叶轮用旋转域，蜗壳用静止域。CFX-Pre 里右键 Insert > Domain，叶轮域 Domain Motion 选 Rotating，转速按实际给，单位rad/s或者rpm看清楚。蜗壳域保持 Stationary。两个域之间的 Domain Interface 是关键设置，稳态阶段用 Frozen Rotor 先跑一个初始流场出来，残差降到1e-5之后切到瞬态。瞬态计算必须用 Transient Rotor-Stator，Frozen Rotor 在瞬态里只做参考系转换不做时间推进，跨交界面传递的流动信息是冻结的，叶轮转动过程中动静干涉的真实脉动会被抹掉。

时间步长按叶轮旋转角度算。转速2700r/min的时候，每步转3°对应的时间步是0.006s，转一圈120步。这个步长下CFL数在1附近，物理时间尺度和湍流时间尺度能匹配上。转速高的话步长还得缩，6000r/min的泵用每步2°算过，时间步0.00023s，采样频率4348Hz，奈奎斯特频率2174Hz，第十阶叶频谐波差不多1450Hz，频率混叠够不着。

CEL表达式写转速变化。启动过程用 DEFINE_TRANSIENT_PROFILE 或者直接在表达式里用 step() 函数写分段。线性加速的表达式大概是 Nt = (n/Ttol)*t*step(Ttol-t)+n*(1-step(Ttol-t))，Ttol 是从零到额定转速的时间，n 是额定转速。三种启动方式里，恒定加速度的扬程上升最快，软启动的扬程先降后升，启动瞬间扬程大于设计值，多出来的那部分是瞬态附加项。

压力脉动监测点在 CFX-Pre 的 Output Control 里建 Monitor Point，叶轮进口、叶轮出口、蜗壳隔舌附近各放一个。瞬态结果里看主频，设计工况下蜗壳监测点主频等于叶片通过频率，叶轮和蜗壳交界面附近脉动强度最大，离叶轮出口越远幅值越小。监测到的数据用 Chart 画时域图，FFT在CFX-Post里做 Tools > Macro Calculator 或者导到MATLAB里处理。前几个旋转周期不算，等相邻周期的转矩和出口压力波动小于1%之后再取数据，初始瞬态还没耗散掉。

路径汇总：域设置 Insert > Domain，叶轮旋转蜗壳静止。交界面 Frozen Rotor 跑稳态，瞬态切 Transient Rotor-Stator。时间步长按旋转角度换算，每步2°到3°。监测点在 Output Control > Monitor Point 建。瞬态结果读 CFX-Post，时域波形加FFT看主频，叶频是主要关注频率。

去看看

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删