4. Fluent旋转机械仿真工作流程

旋转机械仿真与普通内流仿真最大的区别在于：叶轮在旋转，静止部件不动，两者之间存在相对运动。因此整个工作流程必须围绕“如何描述旋转”这条主线展开。本节以一台单级离心泵为例（叶轮转速 2900 r/min，设计流量 50 m³/h，介质为常温水），把从几何到性能曲线的完整链路串起来，并给出每一步的关键输入、操作要点和结果判读依据。

一、几何处理：把“能画出来”变成“能算出来”

旋转机械的几何通常来自 CAD 或叶轮设计软件，但直接导入 Fluent 前必须做三件事。

第一是抽取流体域。叶轮本体是固体，仿真关心的是它“扫过”的水体。做法是用叶轮外轮廓包络出旋转域，再减去叶片实体，得到叶轮通道流体域。进口段、蜗壳、出口段则分别抽取为静止域。最终得到三个计算域：进口域（静止）、叶轮域（旋转）、蜗壳域（静止）。

第二是确定交界面。叶轮域与蜗壳域之间、叶轮域与进口域之间各有一个圆柱面，这就是后面要设为 interface 的位置。交界面的直径应等于叶轮外径，位置尽量避开叶片尾缘，避免交界面直接切过叶片。

第三是几何清理。叶片前缘、尾缘常有微小倒角或尖角，网格质量会急剧恶化。通常把尾缘厚度控制在 0.5～1 mm 以上，前缘做圆角过渡。若几何存在干涉、缝隙或重复面，需在 SpaceClaim 或 DesignModeler 中用“填充”“修复”工具处理，否则网格划分阶段会出现负体积。

二、网格划分：旋转域的质量决定成败

旋转机械推荐使用 Fluent Meshing 的 Watertight Geometry 工作流，或 ANSYS Meshing 中的扫掠/多面体方案。核心原则是：叶轮域用六面体或六面体核心多面体，静止域可适当放宽。

以叶轮域为例，叶片表面设置边界层，第一层高度按 y+ 目标估算。若采用 SST k-ω 模型，目标 y+≈1，第一层高度约 0.02 mm，增长比 1.2，层数 10～15。叶片表面网格尺寸 1.5 mm，轮毂和轮盖面 2 mm，交界面 2 mm。这样单叶轮域网格量约 120 万，三个域合计约 300 万，属于工程可接受范围。

网格质量检查重点看三项：最小正交质量（Orthogonal Quality）应大于 0.15，最好大于 0.2；最大长宽比控制在 100 以内；交界面两侧网格尺寸尽量接近，避免插值误差过大。若交界面两侧尺寸差异超过 3 倍，建议在交界面附近做局部加密。

三、求解器设置：把“旋转”翻译成 Fluent 能懂的语言

打开 Fluent，选择 3D、压力基、瞬态或稳态。对于泵、风机这类定常工况，先用 MRF（多重参考系） 快速得到初场，再视需要切换滑移网格。MRF 的做法是：在 Cell Zone Conditions 中把叶轮域设为 Moving Reference Frame，转速 2900 r/min，旋转轴为 Z 轴；进口域和蜗壳域保持静止。

交界面处理：Mesh Interfaces 中新建一对 interface，分别选择叶轮域侧和蜗壳域侧的面。MRF 下选择“Mixing Plane”或“Frozen Rotor”均可，前者做周向平均，后者保持相对位置冻结。若关注叶片通过频率的非定常脉动，则改用滑移网格，时间步长按每步转过 1°～2° 设置，2900 r/min 对应每步约 5.7×10⁻⁵ s。

湍流模型选 SST k-ω，它兼顾近壁面与主流区，适合有逆压梯度的叶轮通道。压力-速度耦合用 Coupled，压力离散用 Second Order，动量、湍动能、比耗散率均用 Second Order Upwind。进口给质量流量 50 m³/h（约 13.89 kg/s），出口给静压 0 Pa（表压），壁面绝热无滑移。

收敛判断不能只看残差。除连续性、动量残差降到 1×10⁻⁴ 以下外，更应监控出口总压、叶轮扭矩、扬程三条物理量，当它们随迭代步数趋于平稳（波动小于 0.5%）时才算收敛。若残差下降但扬程持续漂移，说明尚未稳定，需继续迭代或减小松弛因子。

四、后处理与性能曲线提取

收敛后先看整体流场：在 CFD-Post 或 Fluent 后处理中创建叶轮中间截面的压力云图，正常应看到从进口到出口压力逐级升高，叶片工作面压力高于背面。再看速度流线，蜗壳内应形成稳定的周向流动，隔舌附近不应出现大面积回流。

性能参数按以下方式提取：

• 扬程 H = (p_out − p_in)/(ρg) + (v_out² − v_in²)/(2g) + Δz。示例中进出口同高，Δz=0，若算得 p_out−p_in = 3.2×10⁵ Pa，则 H ≈ 32.6 m。
• 轴功率 P = M·ω，M 为叶轮所有叶片扭矩之和，ω = 2πn/60 = 303.7 rad/s。若 M = 5.8 N·m，则 P ≈ 1.76 kW。
• 效率 η = ρgQH/P。Q = 50/3600 = 0.01389 m³/s，代入得 η ≈ 0.01389×9.81×32.6/1.76 ≈ 0.63，即 63%。

要得到完整性能曲线，需改变进口流量重复计算。建议取 0.6Q、0.8Q、1.0Q、1.2Q 四个工况点，每个工况重新初始化或从上一点结果续算。把四个点的 (Q, H)、(Q, η)、(Q, P) 分别描点连线，即得扬程-流量曲线、效率-流量曲线和功率-流量曲线。示例数据如下表，均为教学假设值，实际应以计算收敛结果为准。

五、流程闭环与常见返工点

整个流程并非单向。若后处理发现扬程明显低于设计值，常见原因有三类：一是交界面设置错误，叶轮域未真正旋转；二是网格在叶片前缘过粗，导致流动分离被抹平；三是出口边界离蜗壳太近，产生人为反射。此时应回到对应阶段修正，而不是在求解器里反复调参。

建议在项目开始时建立一份检查清单：几何是否封闭、交界面是否成对、叶轮域转速方向是否符合右手定则、边界条件量纲是否一致、收敛监控量是否平稳。把这五项确认到位，旋转机械仿真的一次成功率会显著提高。