旋转机械的几种流体仿真方式

发表评论

引言：当”转起来”遇上仿真

离心泵叶轮的高速旋转、涡轮发动机压气机级间的复杂流动、风力发电机叶片切割气流……旋转机械几乎无处不在，支撑着现代工业的每一个角落。然而，当工程师们试图用数值仿真去”看见”这些旋转部件内部的流动时，却面临一个根本性难题：

传统 CFD 方法假设计算域静止，而旋转机械的核心就是”动”。

如何在计算机中忠实还原叶片的旋转、叶道中的分离、动静叶之间的相互干扰？这一问题催生了一系列专门针对旋转流动的仿真方法。不同方法在精度、计算成本、适用边界上各有取舍，选择合适的方法往往决定了仿真项目的成败。

本文将系统梳理旋转机械 CFD 中最主流的几种仿真方式，深入讲解其物理机制、适用场景和工程要点，并附上选型对比，希望对从事泵、风机、压气机、涡轮等旋转机械仿真的工程师有所帮助。

一、问题的本质：旋转域与静止域的耦合

在进入具体方法之前，先理解核心矛盾。

旋转机械通常由转子（旋转域）和定子/外壳（静止域）两部分构成。转子（叶轮、叶片、转鼓）以角速度 ω 旋转，而进出口管道、蜗壳等结构静止不动。两个域之间存在动静交界面，流体需要在这里完成参考系的切换。

在数学上，描述旋转流动有两个参考系选择：

  • 惯性参考系（Lab Frame）：坐标系固定，旋转边界随时间变化，描述绝对运动
  • 旋转参考系（Rotating Frame）：坐标系随转子一起转，转子在此系中静止，但 N-S 方程中引入科里奥利力和离心力等附加体积力项

∂u′∂t+(u′⋅∇)u′=−1ρ∇p+ν∇2u′−2ω×u′⏟科氏力−ω×(ω×r)⏟离心力\frac{\partial \mathbf{u}’}{\partial t} + (\mathbf{u}’ \cdot \nabla)\mathbf{u}’ = -\frac{1}{\rho}\nabla p + \nu \nabla^2 \mathbf{u}’ – \underbrace{2\boldsymbol{\omega} \times \mathbf{u}’}_{\text{科氏力}} – \underbrace{\boldsymbol{\omega} \times (\boldsymbol{\omega} \times \mathbf{r})}_{\text{离心力}}∂t∂u′​+(u′⋅∇)u′=−ρ1​∇p+ν∇2u′−科氏力2ω×u′​​−离心力ω×(ω×r)​​

不同仿真方法的本质区别，正在于如何处理这两个参考系的切换，以及是否在时间上推进。

二、多参考系法（MRF）

2.1 基本原理

多参考系法（Multiple Reference Frame，MRF）是旋转机械 CFD 中最常用、最简单的稳态仿真方法。

其核心思想是：将计算域划分为旋转子域和静止子域，在各自的参考系下分别求解 N-S 方程。旋转子域在旋转参考系中求解（叶轮在此系中静止，但附加体积力存在），静止子域在惯性参考系中求解。两个子域通过内部交界面进行流量和压力的插值交换，整个计算在稳态条件下进行。

关键点：MRF 假设转子与定子之间不存在时间相关的相互作用，交界面上的速度和压力通过坐标变换直接匹配，而不随时间推进。这意味着转子被”冻结”在某个位置，动静叶的相对位置固定不变。

2.2 适用场景

MRF 最适合以下情况：

  • 工况预测与效率曲线绘制：叶轮泵、离心式压气机的性能曲线（Q-H 曲线、效率曲线）
  • 动静干扰相对较弱的机型：径流式叶轮、轴流风机（非对转）
  • 设计优化阶段：需要大量工况点快速评估，计算资源有限
  • 无需时域信息的场景：只关心平均性能，不关心叶片通过频率（BPF）等非定常效应

典型应用设备：离心泵、离心风机、离心压气机、螺旋桨（远场近似）。

2.3 优缺点

优点：

  • 计算量小，收敛速度快（通常只需稳态迭代）
  • 网格无需随时间变化，前处理简单
  • 各大商业软件（Fluent、CFX、StarCCM+、OpenFOAM）均有成熟实现

缺点：

  • 冻结转子假设：动静叶相对位置固定，无法捕捉叶片通过频率等非定常效应
  • 动静干扰强时误差明显（如导叶与叶轮间隙小、通流面积突变大）
  • 不适合有强时变效应的工况（如失速、喘振）
  • 交界面的插值方式对结果有影响，需谨慎设置

2.4 软件实现要点

  • ANSYS Fluent：在 Cell Zone Conditions 中将旋转子域设置为 MRF，指定旋转轴和角速度；交界面设置为 Interface，激活 Mesh Interface 功能
  • ANSYS CFX：使用 Frozen Rotor 接口类型（即 MRF 的 CFX 叫法），在 Domain Interface 中选择 Frozen Rotor
  • OpenFOAM：使用 MRFZones 字典在 constant/ 目录下定义旋转区域，配合 steadyMRFFoam 或 simpleFoam（带MRF选项）求解器
  • StarCCM+：在 Physics Model 中启用 Moving Reference Frame，为旋转子域指定 Rotating Frame

网格建议：交界面两侧网格尺寸应尽量接近，避免插值误差；旋转子域网格节点无需匹配（非共形网格可用），但面积保持量应满足一致性。

三、冻结转子法与混合面法

3.1 冻结转子法（Frozen Rotor）

冻结转子法本质上是 MRF 的另一种说法，在 ANSYS CFX 中被明确称为 Frozen Rotor 接口。其与 MRF 的区别在于：

  • MRF 常用于整机稳态，计算域为全周或扇区
  • Frozen Rotor 强调动静交界面处的坐标系变换处理方式——转子被”冻结”在特定角度位置，交界面两侧流场直接做参考系变换

当关心某个特定转子位置下的载荷分布时（如叶片在最不利位置的受力），Frozen Rotor 是快速评估工具。

3.2 混合面法（Mixing Plane）

混合面法（Mixing Plane Method）是一种在动静交界面上进行周向平均的稳态处理方式，常用于多级轴流/离心压缩机仿真。

原理：将动静交界面上的流量、压力、速度等参数沿周向方向进行平均，消去周向不均匀性，然后将平均值作为边界条件传递给对侧域。这样，转子”看到”的是定子出口的平均流场，定子”看到”的是转子出口的平均流场。

优势：

  • 比 MRF 更物理合理——动静两侧交换的是”周向平均”信息，消除了由于某个特定转子位置引起的非物理扰动
  • 特别适合多级轴流压气机/涡轮，各级之间的干扰通过混合面传递

局限性：

  • 周向平均本质上是一种稳态近似，无法反映叶片通过频率（BPF）等真正的非定常效应
  • 对于动静干扰强烈（如导叶与叶片通道数接近，强势流干扰）的机型，误差较大

软件实现：CFX 中的 Stage 接口类型即混合面法；Fluent 通过 Mixing Plane Model 实现；OpenFOAM 的 cyclicGgi 可实现类似功能。

四、滑动网格法（Sliding Mesh）

4.1 基本原理

滑动网格法（Sliding Mesh Method）是目前精度最高、最常用的非定常旋转机械仿真方法。

与 MRF”冻结”转子不同，滑动网格法让旋转子域的网格真实地随时间旋转。在每个时间步，旋转子域的网格相对于静止子域旋转一个角度 Δθ = ω·Δt，动静交界面（滑移面）两侧的节点位置不断变化，通过非共形界面插值（如加权插值、区域耦合）实时交换数据。

整个仿真在瞬态条件下推进，能够真实捕捉：

  • 叶片通过频率（BPF）引起的压力脉动
  • 动静叶之间的势流干扰和尾迹干扰
  • 叶轮内部流动的时变特征（如旋转失速的初生与发展）

4.2 适用场景

  • 需要非定常信息的所有旋转机械场景
  • 叶片通过频率分析、噪声预测（气动声学）
  • 旋转失速、喘振、流动不稳定性研究
  • 叶轮与蜗壳/扩压器之间干扰强烈的情况（如叶轮出口与导叶进口间隙小）
  • 多级压缩机/涡轮的非定常级间干扰研究

典型应用：轴流风机（气动噪声分析）、涡轮机组（叶片载荷时变）、水力透平（压力脉动）、混流泵（驼峰区失稳研究）。

4.3 优缺点

优点：

  • 物理保真度高：完整捕捉旋转引起的所有非定常效应
  • 可获得时域压力脉动、速度脉动信号，可进行 FFT 分析
  • 无”冻结”假设，动静干扰真实反映
  • 对强动静干扰问题（如低比转速泵、密集导叶涡轮）尤其准确

缺点：

  • 计算量大：需要推进足够的旋转圈数直至周期收敛（通常 5-20 圈），每圈若干时间步
  • 时间步需满足 CFL 稳定性条件（ΔθCFL ≈ 0.5°~2° per step 为常见经验值）
  • 前处理更复杂：需要明确设置滑移面、时间步、数据采集窗口
  • 存储需求大：瞬态场文件数量多

4.4 软件实现要点

  • ANSYS Fluent：在 Cell Zone Conditions 中设置 Mesh Motion（而非 MRF），激活 Dynamic Mesh 中的 Rigid Body Motion；交界面设置为 Mesh Interface，勾选 Periodic Repeats（如扇区模型）
  • ANSYS CFX：动静交界面选择 Transient Rotor Stator 类型，在 Analysis Type 中设置瞬态求解
  • OpenFOAM：使用 pimpleFoam 或 pimpleDyMFoam 配合 dynamicMeshDict 中的 solidBodyMotionFvMesh 类型，设置旋转轴和角速度；交界面使用 cyclicAMI（Arbitrary Mesh Interface）处理非共形滑移面
  • StarCCM+：在 Physics Model 中启用 Rigid Body Motion，在 Interface 中设置 Sliding 类型的 In-place Interface

时间步选择建议：

  • 每个叶片通过一个网格单元的时间 ≈ 合理时间步上限
  • 经验值：每转360步（即1°/step）用于大多数工况；关注BPF时可细化至0.5°/step
  • 通常需要推进 3-5 个旋转周期达到统计收敛，再取后续周期作为有效数据

五、动网格法（Dynamic Mesh）

5.1 基本原理

动网格法（Dynamic Mesh Method）是比滑动网格更通用的瞬态网格运动技术。在滑动网格中，旋转域的网格拓扑保持不变，仅做整体刚体旋转；而动网格允许网格本身发生变形、重划分，能够处理更复杂的运动形式。

动网格技术通常包含三种网格处理策略：

  1. 弹簧光顺（Spring Smoothing）：将网格边视作弹簧，节点随边界运动而”弹性变形”，适合小变形问题
  2. 局部重网格（Local Remeshing）：当网格质量劣化（高偏斜率）时，自动删除并重新生成局部网格单元，适合大变形
  3. 层铺法（Layering）：在运动边界附近周期性地增加或删除网格层，常用于活塞式运动（往复压缩机）

5.2 适用场景

动网格法在旋转机械中常用于以下特殊场景：

  • 往复式压缩机/内燃机：活塞往复运动，气缸容积周期性变化
  • 叶片可调机构：变桨距螺旋桨、可调导叶（IGV），叶片角度随工况变化
  • 旋转机械与外部流场的强耦合：如风力机叶片大变形时的气弹分析
  • 非规则旋转：进动、章动等非恒定转速运动

对于常规恒速旋转（泵、压气机、涡轮），通常优先选用滑动网格（更简单、更稳定），动网格更多用于以上特殊情形。

5.3 优缺点

优点：

  • 通用性极强，可处理任意形式的运动（旋转、平动、变形的组合）
  • 能模拟叶片变形、间隙变化等几何非线性效应

缺点：

  • 前处理复杂：需要精细配置运动函数（UDF/表达式），网格质量控制难
  • 重网格引入数值误差：局部重网格后网格拓扑改变，保守量插值带来误差
  • 计算成本高：每个时间步可能需要重划分，速度比滑动网格慢
  • 调试和收敛控制难度大，工程应用门槛高

5.4 软件实现要点

  • ANSYS Fluent：Dynamic Mesh 功能最完善，支持 UDF 自定义运动、Profile 文件驱动运动；六自由度（6-DOF）求解器可与流体计算双向耦合
  • OpenFOAM：dynamicFvMesh 库提供多种动网格选项，sixDoFRigidBodyMotion 适合流固耦合场景
  • StarCCM+：通过 DFBI（Dynamic Fluid-Body Interaction）实现六自由度运动，图形化设置较友好

六、周期性边界条件与扇区模型

6.1 基本原理

对于叶轮机械，转子和定子通常在周向上是周期对称的——相邻两个叶片通道之间的流场完全相同。利用这一对称性，可以只计算一个叶片通道（扇区），大幅降低计算量。

旋转周期性边界条件将扇区一侧边界上的流场变量（速度、压力、湍流量）与另一侧边界关联：ϕ(r,θ+Δθ,z)=ϕ(r,θ,z)\phi(r, \theta + \Delta\theta, z) = \phi(r, \theta, z)ϕ(r,θ+Δθ,z)=ϕ(r,θ,z)

其中 Δθ = 2π/N（N 为叶片数），这一关系在稳态情况下精确成立。

6.2 与 MRF/滑动网格的组合使用

周期性边界条件本身不是独立的旋转处理方法，而是配合 MRF 或滑动网格使用的网格规模缩减技术：

  • 扇区 MRF：单叶片通道 + MRF，稳态计算，计算量极小，是设计优化的首选
  • 扇区滑动网格：单叶片通道瞬态计算，要求转子与定子叶片数之比为简单有理数（如 2:3），通过”时间移位”（Time-Shift / Phase-Lag）技术处理相位差

**相位滞后法（Phase-Lag Method / Time-Transformation）**是扇区非定常仿真的重要突破：即使叶片数之比不是简单分数，通过将时间延迟引入周期边界，也能用单通道计算模拟全周非定常流动，大幅节省计算资源。

  • ANSYS CFX：内置 Profile Transformation 和 Time Transformation 两种相位滞后方法
  • Fluent：SCORE（Single Channel Oscillating Reference Expander）方法支持非整数叶片比的扇区计算

七、格子玻尔兹曼方法（LBM）在旋转机械中的应用

7.1 基本原理

格子玻尔兹曼方法（Lattice Boltzmann Method，LBM）是近年来兴起的 CFD 方法，基于介观尺度的粒子分布函数演化方程，而非直接求解 N-S 方程。LBM 天然并行、边界处理简便，在工程 CFD 中逐渐崭露头角。

在旋转机械仿真中，LBM 主要通过浸入边界法（Immersed Boundary Method, IBM）或非惯性系格子处理旋转叶片，避免了传统 FVM 方法中复杂的动网格或 MRF 设置。

7.2 工程应用现状

LBM 在旋转机械中目前主要应用于：

  • 气动噪声预测（风机叶片噪声、汽车冷却风扇噪声）：LBM 在低马赫数湍流和声场耦合方面有天然优势
  • 汽车散热风扇系统：PowerFLOW（现为 Simulia 产品）在汽车行业有成熟商业应用
  • 电机冷却流道：复杂几何下的低速湍流

代表软件：Simulia PowerFLOW（基于 LBM）；开源方面 OpenLB、Palabos 等。

LBM 在高马赫数、强旋转、多相流等场景还有一定局限，目前尚未成为旋转机械仿真的主流，但其在噪声和低速湍流领域的精度优势值得关注。

八、方法对比与选型建议

8.1 综合对比表

方法求解类型计算量精度非定常信息适用场景主要软件支持
MRF（多参考系）稳态★☆☆中等✗性能曲线、设计优化Fluent, CFX, StarCCM+, OpenFOAM
冻结转子（Frozen Rotor）稳态★☆☆中等✗特定位置载荷评估CFX, Fluent
混合面（Mixing Plane）稳态★★☆较高✗多级压气机/涡轮CFX（Stage）, Fluent
滑动网格（Sliding Mesh）瞬态★★★高✓非定常分析、噪声、失速Fluent, CFX, StarCCM+, OpenFOAM
动网格（Dynamic Mesh）瞬态★★★高✓往复机械、变桨距、气弹Fluent, StarCCM+, OpenFOAM
相位滞后扇区法瞬态（扇区）★★☆高✓多级机（计算资源受限）CFX, Fluent（SCORE）
LBM瞬态★★★高（低马赫）✓气动噪声、低速流PowerFLOW, OpenLB

8.2 工程选型决策树

你的仿真目标是什么？
│
├── 只需要平均性能（效率、扬程、流量）？
│   └── → MRF + 扇区模型（最快）
│       → 多级机器 → 混合面法
│
├── 需要非定常信息（压力脉动、BPF、噪声）？
│   └── 计算资源充足？
│       ├── 是 → 滑动网格（全圆或扇区）
│       └── 受限 → 相位滞后扇区法
│
├── 往复运动或变几何（活塞、变桨距）？
│   └── → 动网格法
│
└── 气动噪声为主要目标？
    └── → 考虑 LBM（PowerFLOW）或 Fluent DES + FW-H 方程

8.3 几点实战经验

1. MRF 结果作为滑动网格初始场

对于滑动网格（非定常），先跑一个 MRF 稳态解，将结果作为瞬态初始场，可以显著加速非定常迭代的收敛，减少”暖机”圈数。

2. 交界面设置是关键

动静交界面的设置（共形 vs 非共形、插值方式）对结果影响明显。推荐在交界面两侧保持相近的网格尺寸；对于 MRF，优先选择周向面积一致的圆柱交界面位置，避免放在叶片尾迹区域。

3. 时间步长与叶片数的关系

滑动网格的时间步建议满足：每个叶片通过一个相邻叶片通道的时间至少分成 20 步以上（通常取叶片通过时间的 1/20~1/50 作为时间步）。叶片数多时可适当放宽。

4. 验证与确认（V&V）

旋转机械 CFD 结果应与试验台数据（扬程、效率特性曲线）进行系统验证。常见验证指标：额定工况效率误差 < 2%，压力脉动主频幅值误差 < 10%。

九、未来趋势展望

9.1 GPU 加速与大规模并行

随着 NVIDIA GPU 和 AMD GPU 在 HPC 集群中的普及，传统 CPU 集群运行的 MPI 并行 CFD 正在向 GPU 异构计算迁移。ANSYS Fluent 2024 版本已实现完整的 GPU Native 求解，滑动网格的单工况计算时间有望从数天压缩至数小时，这将推动非定常仿真（滑动网格）进一步普及。

9.2 基于机器学习的代理模型

深度学习正在进入旋转机械设计流程。通过对大量 MRF 仿真数据进行训练，可以构建叶轮性能预测的神经网络代理模型（Surrogate Model），实现毫秒级的工况预测，支持多目标优化算法（如 NSGA-II）的大规模搜索。目前已有团队在离心泵和轴流风机叶型优化中取得突破。

9.3 LES/DNS 在旋转机械中的应用扩展

随着算力提升，基于大涡模拟（LES）乃至直接数值模拟（DNS）的旋转机械仿真正在从学术走向工程。尤其是在旋转失速、喘振起始等强非定常、强分离场景，URANS 湍流模型已逐渐暴露出局限性，Hybrid RANS-LES（如 DES、SAS、WMLES）方法正在成为高精度工业仿真的新选择。

9.4 数字孪生与在线仿真

将快速 CFD 模型（简化的 MRF 或代理模型）与实时传感器数据融合，构建旋转机械的流体动力学数字孪生（CFD Digital Twin），是工业互联网时代的重要方向。这要求仿真模型不仅精度足够，还需具备近实时的响应速度。

结语

旋转机械的流体仿真，本质上是在精度与效率之间找到最优平衡。MRF 以极小的计算代价提供设计阶段所需的性能评估；滑动网格以时域真实性满足非定常分析的严苛需求；动网格和相位滞后法则覆盖了更复杂的工程边界；而 LBM 等新兴方法正在特定赛道上展现独特价值。

没有”最好”的方法，只有最适合当前工程问题的选择。理解每种方法背后的物理假设和数学本质，才能在面对新问题时做出有据可依的决策——这也是每一位 CFD 工程师从”会用软件”到”真正懂仿真”的必经之路。

如果本文对你有帮助，欢迎点赞收藏，也欢迎在评论区分享你在旋转机械仿真中的实战经验！

本文配图均为示意性图解，实际仿真界面因软件版本不同略有差异。文中涉及的软件功能描述基于 2024 年主流商业版本。

参考资料方向（供扩展阅读）：

  1. ANSYS Fluent Theory Guide – Moving Reference Frame Formulation
  2. ANSYS CFX Solver Theory Guide – Rotating Frame of Reference
  3. Lakshminarayana, B. Fluid Dynamics and Heat Transfer of Turbomachinery, Wiley
  4. Casey, M. & Winterbone, D. The Thermodynamics and Gas Dynamics of Internal Combustion Engines
  5. OpenFOAM Documentation – AMI (Arbitrary Mesh Interface) for Rotating Machinery

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注