3. 旋转机械性能参数体系
在 Fluent 中做旋转机械仿真,真正决定“算得对不对、有没有工程意义”的,不是网格有多密,而是你是否清楚每个性能参数的定义、参照系和相互换算关系。本节把泵、风机、压缩机、水轮机等叶轮机械共用的参数体系梳理成一条可计算、可验证的链条,并说明它们在 CFD 后处理中如何提取。
一、基本性能参数:从流量到效率
体积流量 Q 是单位时间通过叶轮机械的流体体积,单位 m³/s 或 m³/h。在 Fluent 中通常用进口面或出口面的面积分获得:Report → Fluxes → Mass Flow Rate 得到质量流量 ṁ,再除以当地密度得到 Q。对可压缩介质(如压缩机),进出口密度不同,必须明确 Q 指的是进口体积流量还是出口体积流量,工程上泵与风机习惯用进口体积流量。
扬程 H 是泵最核心的参数,定义为单位重量流体获得的能量,单位 m。对不可压缩流体:
H = (p_out - p_in)/(ρg) + (v_out² - v_in²)/(2g) + (z_out - z_in)
三项分别是静压头、动压头和位势头。Fluent 中可在进出口面用面积加权平均取总压 Total Pressure,则 H = (P_t,out − P_t,in)/(ρg),位势头在进出口高度差明显时需单独补入。风机常用全压 P_t(Pa)代替扬程,二者关系为 P_t = ρgH。
轴功率 P_shaft 是原动机传给叶轮的机械功率,P_shaft = M·ω,其中 M 为扭矩、ω 为角速度(rad/s)。有效功率 P_e 是流体实际获得的功率,P_e = ρgQH(泵)或 P_e = Q·P_t(风机)。效率 η = P_e / P_shaft,它把水力损失、容积损失、机械损失全部打包成一个可对比的指标。在 CFD 中,P_shaft 由叶轮壁面扭矩积分得到,P_e 由进出口总压差与流量得到,因此 η 是后处理中最容易“算错”的参数——务必确认扭矩积分方向与旋转方向一致。
二、扭矩、推力与滑移系数
扭矩 M 在 Fluent 中通过叶轮所有叶片与轮毂壁面的压力与黏性力对旋转轴的力矩求和获得,可用 Report → Forces → Moment,指定旋转轴原点和方向矢量。注意:力矩参考点必须落在旋转轴线上,否则会引入虚假的力臂误差。
轴向推力 F_axial 是流体对叶轮产生的沿轴向的合力,直接决定轴承选型与平衡孔设计。离心泵中叶轮前后盖板压力分布不均,会产生指向进口的轴向力;Fluent 中同样用 Forces 报告,方向取轴向单位矢量。多级泵或高压比压缩机中,推力可达数万牛,是必须校核的量。
滑移系数 σ 描述有限叶片数导致的出口气流偏离理想无限叶片假设的程度:
σ = c_u2_actual / c_u2_ideal = 1 - π·sinβ₂ / Z (Stodola 近似)
其中 Z 为叶片数,β₂ 为叶片出口安装角。σ 越小,说明滑移越大、实际扬程越低。CFD 的价值正在于:它不需要任何滑移经验公式,直接解出真实出口速度三角形,从而反算出 σ,用于校核一维设计结果。后处理时在叶轮出口取周向平均的绝对速度切向分量 c_u2,与欧拉方程给出的理想值比较即可。
三、无量纲系数与相似定律
把性能参数无量纲化,才能在不同尺寸、不同转速的机器之间比较。常用系数包括:
- 流量系数 φ = Q /(u₂·D₂²·π/4),反映通流能力;
- 扬程系数 ψ = 2gH / u₂²,反映做功能力;
- 功率系数 λ = P /(ρ·u₂³·D₂²);
- 比转速 n_s = n·√Q / H^(3/4),是决定叶轮几何形状的相似准则。
其中 u₂ = πD₂n/60 为叶轮出口圆周速度。相似定律则给出同一台机器在变转速下的换算关系:Q ∝ n,H ∝ n²,P ∝ n³。这三条关系是泵与风机变频调速、性能曲线外推的理论基础,也是 CFD 变转速算例的验证依据——如果你在 Fluent 中把转速从 1450 r/min 提到 2900 r/min,算出的扬程应约为原来的 4 倍,偏差超过 5% 就说明网格或湍流模型有问题。
四、后处理实操:一次完整的参数提取
以一台离心泵为例,设计工况 Q = 50 m³/h、n = 2900 r/min、D₂ = 0.22 m、叶片数 Z = 6、β₂ = 30°。在 Fluent 中完成 MRF 稳态计算后,按以下顺序提取参数:
- 进口面报告质量流量,除以密度得 Q,与设计值对比,偏差应小于 2%;
- 进出口面报告面积加权平均总压,代入 H 公式,得到扬程;
- 叶轮所有壁面报告绕 Z 轴的力矩 M,乘以角速度 ω = 2πn/60 = 303.7 rad/s,得轴功率;
- 计算 η = ρgQH / (Mω),与设计效率对比;
- 叶轮出口取周向平均 c_u2,与理想值 u₂ − c_m2·cotβ₂ 比较,得滑移系数 σ;
- 叶轮壁面报告轴向力,校核轴承载荷。