2. 旋转机械流体力学基础
旋转机械流场仿真的底层依据是三大守恒方程，以及它们在旋转坐标系下的改写形式。理解这些方程的关键，不是背公式，而是弄清每个方程在 Fluent 中对应哪些设置、哪些项会因为旋转而多出来、以及相对速度与绝对速度在什么位置切换。本节先建立方程与仿真设置之间的对应关系，再重点讲清旋转坐标系带来的附加项和速度转换。
一、三大守恒方程的物理含义与 Fluent 对应
连续性方程表达质量守恒：控制体内密度的变化率加上通过表面的质量净流出为零。对不可压缩流动，密度为常数，方程退化为速度散度为零，即流入等于流出。在 Fluent 中，这一方程由求解器自动满足，用户需要保证的是边界条件质量流量闭合——例如一个离心泵叶轮，进口给定质量流量，出口若给定压力，则计算域内质量守恒由压力场自动调整；若进出口都给定流量，则必须保证两者相等，否则残差无法收敛。
动量方程表达牛顿第二定律：控制体内动量的变化率等于对流出的动量通量、表面压力、黏性应力和体积力之和。对旋转机械，体积力中最重要的是离心力和科氏力，它们只在旋转坐标系下出现。动量方程在 Fluent 中对应动量离散格式、压力-速度耦合算法（SIMPLE、SIMPLEC 或 Coupled）以及湍流模型引入的雷诺应力项。选择二阶迎风格式通常能降低数值扩散，但会牺牲稳定性，旋转机械中常先用一阶迎风启动，再切换到二阶迎风。
能量方程表达热力学第一定律：内能与动能的变化率等于导热、做功和外部热源之和。对大多数水泵、风机、压缩机，若马赫数低于 0.3 且温升可忽略，可以关闭能量方程以节省计算量；但对高速离心压缩机、涡轮，必须开启能量方程，并考虑黏性耗散产生的温升。Fluent 中能量方程的开关位于模型设置面板，开启后还需指定流体材料随温度变化的比热和导热系数。
二、旋转坐标系下的控制方程
旋转机械仿真通常不采用绝对坐标系直接模拟叶轮转动，而是把计算域固定在旋转参考系中，让叶轮静止、流体相对叶轮流动。这样做的代价是动量方程中多出两项：离心力和科氏力。设旋转角速度为 ω，相对速度为 w，则动量方程中附加的体积力为：
F_add = -ρ [ 2ω × w + ω × (ω × r) ]
其中：2ω × w 为科氏力，ω × (ω × r) 为离心力
ρ 为密度，r 为位置矢量离心力沿半径向外，对叶轮做功，是泵和压缩机提升压力的主要来源；科氏力垂直于相对速度，不直接做功，但会改变流道内的速度分布，使叶轮出口速度出现偏转。在 Fluent 中，旋转参考系通过 Cell Zone Conditions 设置，选择 Frame Motion 并指定旋转轴原点、方向和角速度。若采用多参考系模型（MRF），则叶轮区域设为旋转参考系，静止区域设为绝对参考系，交界面处做速度转换；若采用滑移网格，则叶轮区域网格随时间转动，交界面处做通量插值。
三、相对速度与绝对速度的转换
旋转坐标系下求解得到的是相对速度 w，而工程上关心的扬程、流量、力矩通常基于绝对速度 v。两者关系为：
v = w + ω × r
其中：v 为绝对速度，w 为相对速度，ω × r 为牵连速度在 Fluent 后处理中，Velocity 默认显示的是绝对速度还是相对速度，取决于参考系设置。查看叶轮内部流线时，应使用相对速度，才能看到流体相对叶片的真实流动；查看蜗壳和进出口时，应使用绝对速度，才能与实验测量的速度场对比。判断方法很简单：若流线在叶轮内呈明显的旋转跟随形态，说明用的是相对速度；若流线在叶轮内几乎不随叶片弯曲，说明用的是绝对速度。
四、方程与 Fluent 设置的对应关系
把方程落到操作上，可以形成一张清晰的对应表。连续性方程对应边界条件的质量流量闭合；动量方程对应压力-速度耦合算法、离散格式和旋转参考系设置；能量方程对应能量开关和材料热物性；旋转附加项对应 Cell Zone Conditions 中的 Frame Motion；速度转换对应后处理中的参考系选择。理解这层对应，就能在残差曲线异常或结果不合理时，快速定位是哪个方程对应的设置出了问题。
| 方程/概念 | Fluent 中的对应设置 | 常见异常 |
|---|---|---|
| 连续性方程 | 进出口边界质量流量、压力出口 | 进出口流量不闭合，连续性残差居高不下 |
| 动量方程 | 压力-速度耦合、离散格式、湍流模型 | 一阶迎风收敛但结果偏平，二阶迎风震荡 |
| 能量方程 | Energy 开关、材料比热与导热系数 | 高速压缩时温升被忽略，密度计算失真 |
| 旋转附加项 | Frame Motion 的旋转轴与角速度 | 旋转轴方向反，压力不升反降 |
| 速度转换 | 后处理参考系选择 | 叶轮流线不跟随叶片，误判为绝对速度 |
五、一个可复现的检查流程
以一台离心泵叶轮为例，假设转速 2900 r/min，叶轮外径 0.2 m，进口质量流量 2 kg/s，介质为水。按以下顺序检查方程与设置是否一致：
- 确认连续性：进口质量流量 2 kg/s，出口若为压力出口，则计算完成后检查出口质量流量是否约为 -2 kg/s（Fluent 中流出为负），偏差应小于 0.5%。
- 确认动量方程：先以一阶迎风、SIMPLE 算法迭代 500 步，观察残差是否降至 1e-3 以下；再切换二阶迎风继续迭代至 1e-5。
- 确认能量方程：若泵送温升小于 1 K，可关闭能量方程；若介质为高温油或高速压缩机，必须开启并检查出口温度。
- 确认旋转设置：在 Cell Zone Conditions 中检查旋转轴原点是否在叶轮中心、方向矢量是否与泵轴一致、角速度是否为 2900 r/min 换算后的 303.7 rad/s。
- 确认速度转换：在 CFD-Post 或 Fluent 后处理中，分别用相对速度和绝对速度显示叶轮中截面流线，对比流线是否跟随叶片弯曲。