冻结转子法

分类:流体分析(CFD) | 统一版 2026-04-06
CAE visualization for frozen rotor theory - technical simulation diagram
冻结转子法 — 理论与应用范围

冻结转子法的理论基础

概要

🧑‍🎓

冻结转子法这个名字很酷,但它具体是做什么的呢?


🎓

旋转体与静止体的相对位置固定进行定常解析的手法。动叶的坐标置于旋转系中,静叶的坐标置于静止系中,但在界面处不进行周向平均,而是"直接"传递信息。


🧑‍🎓

与混合平面法的区别是什么?


🎓
手法界面处理叶片数约束计算成本精度
冻结转子位置固定、直接插值1个节距不同则需要节距比修正存在位置依赖性
混合平面周向质量平均任意节距比平滑化周向变动
滑动网格时间推进中旋转整数比较优最精确

何时使用

🧑‍🎓

冻结转子法适用于什么场景?


🎓

以下情况适用。


  • 初期设计阶段的高速筛选:比较多个形状方案时比混合平面更快
  • 带蜗壳的离心机:蜗壳非轴对称,混合平面难以应用。用冻结转子近似叶片-蜗壳干涉
  • 水力涡轮机:与尾水管的干涉评估

🧑‍🎓

像蜗壳这样混合平面无法使用的情况下很有用啊。


🎓

正是如此。但结果取决于叶片与蜗壳的相对位置。为了精确评估,应在多个相对位置进行计算并求平均,或最终过渡到滑动网格。


CFX中的冻结转子设置

🧑‍🎓

请教一下CFX中如何设置冻结转子。


🎓

将域界面的类型设定为"冻结转子"即可。当节距比不是1:1时,"节距变化"选项会自动进行缩放插值。注意是GGI面两侧的网格节距应大致相同。

咖啡小休 杂谈

涡轮机械CFD的历史——从1970年代的细缝槽法到3D非粘性解析

涡轮机械内部流动的CFD解析开始于1970年代。由于当时计算机能力限制,主要采用2D叶栅解析,无法考虑3D叶片形状。转折点出现在1970年代后期,伴随ANSYS Fluent前身代码的开发,以及Denton(1982)用IBM大型主机实现的3D非粘性(欧拉方程)涡轮机械解析。随后1990年代,Harvey & Denton、Arnone等实现了包含转子-定子非定常干涉的粘性解析,现代涡轮CFD(RANS+滑动网格)的原型得以确立。冻结转子法在这段历史中作为"定常近似的合理首选方案"而应运而生,半个世纪后的当今仍作为设计探索阶段的标准手法继续发挥作用。

冻结转子法的数值计算手法

位置依赖性问题

🧑‍🎓

冻结转子法的结果会因叶片位置而改变,这是真的吗?


🎓

是的。比如在离心泵中,叶片位于蜗壳舌部正面的情况与偏离的情况相比,扬程可能会变化5~10%。这是冻结转子法本质上的局限性。


🧑‍🎓

那就不可信了?


🎓

把单一位置的结果直接作为性能值报告是危险的。推荐做法是在叶片节距的1/3~1/2间隔处计算3~5个相位位置,然后求平均。


节距比修正

🧑‍🎓

什么是节距比?


🎓

是转子和定子的1个节距的角度比。比如转子7片(节距51.4度),定子12片(节距30度),节距比就是51.4/30=1.71。在冻结转子法中,界面两侧必须以某种方式处理这个节距差。


CFX中设定界面的节距比后,会进行周向缩放插值。但当节距比超过2时,精度会大幅下降。


混合平面法与冻结转子法的使用区分

🧑‍🎓

实务中怎样区分使用?


🎓
情况推荐手法
轴对称扩压管/无蜗壳混合平面
带蜗壳的离心泵冻结转子(多相位)→ 滑动网格
多级轴流混合平面
预研设计的参数化冻结转子(快速)
压力脉动·噪声评估滑动网格(必须)
🧑‍🎓

最终还是要用滑动网格验证一下比较放心。


🎓

正是这样。冻结转子是"快速近似",混合平面是"稳定的定常近似",滑动网格是"物理上正确的非定常解析"。这是三者的定位。

咖啡小休 杂谈

冻结转子法的数值设定——界面处理与收敛的诀窍

冻结转子(Frozen Rotor)法是在MRF固定条件下进行定常解析的方法,在转子-定子间周向均匀的假设下成立。实现中"转子-定子界面"的处理决定了精度。此界面处理中进行坐标系变换(从旋转坐标系到静止坐标系),但当流动周向不均匀时(尾迹干涉强烈),假设条件失效,在不同周向位置评估解时会产生"周向依赖性"。对策是在多个周向位置执行冻结转子计算,取周向平均(节距平均)。此外,界面流场变量的插值精度影响整体精度,因此在重要部位进行滑动网格(精度更高)的对比验证是推荐做法。

冻结转子法的实务应用

离心泵的模型构成

🧑‍🎓

离心泵的整体模型怎样组建?


🎓

典型的构成如下。


  • 吸入管道:静止域,无旋转
  • 叶轮:旋转域,MRF或滑动网格
  • 蜗壳:静止域,非轴对称
  • 叶轮-蜗壳界面:冻结转子(定常)或滑动网格(非定常)

由于蜗壳是非轴对称的,无法应用混合平面。这就是冻结转子大有用武之地的原因。


🧑‍🎓

蜗壳的网格怎样生成?


🎓

蜗壳的截面形状沿螺旋方向变化,TurboGrid无法生成。用Ansys Meshing或Fluent Meshing的非结构网格,或STAR-CCM+的自动网格。将截面形状沿扫掠方向排列,形成六面体为主的网格会比较有利。


蜗壳舌部的处理

🧑‍🎓

舌部(截断口)附近有什么难度?


🎓

舌部是叶轮出口流与再循环流相碰撞的区域,压力梯度陡峭。网格必须特别细密。另外,冻结转子中叶片与舌部的相对位置决定了流场,所以压力脉动评估必须用滑动网格。


性能特性曲线的创建

🧑‍🎓

用冻结转子法能生成离心泵的H-Q曲线吗?


🎓

可以,但推荐每个流量点都进行多相位平均。步骤如下。


1. 设计流量处冻结转子计算收敛

2. 叶片每隔10度旋转,进行3~5个相位的重新计算

3. 各相位的扬程进行算术平均

4. 改变流量,重复步骤2-3


🧑‍🎓

工作量相当大啊…


🎓

所以实务中常见的做法是:用冻结转子快速评估设计点附近的1~2个点,最终性能特性曲线通过滑动网格非定常计算来确定。

咖啡小休 杂谈

水力涡轮机(卡普兰型)的冻结转子解析——导叶角最优化

拥有可动导叶(进口导向羽片)的卡普兰水轮机在CFD设计中,冻结转子法用于多个导叶开度(通常10个阶段)与转速的组合,可快速生成性能特性曲线。在全流量·全水头运行条件范围内计算η-Q曲线,为最高效率点(BEP: Best Efficiency Point)附近的设计指导原则奠定基础。每个导叶角度的冻结转子计算(定常RANS)每个工况仅需数小时完成,10个工况的性能特性曲线可在1~2天内生成——而滑动网格同样的计算需要数周。瑞典水力发电设备制造商公开的案例显示,冻结转子+优化可将设计周期从3个月缩短到3周。

冻结转子法的软件对比

求解器的冻结转子实现

🧑‍🎓

不同求解器的冻结转子实现方式不同吗?


🎓

名称与设置方法有差异。


求解器名称/设置节距比修正备注
Ansys CFX域界面 → 冻结转子节距变化自动最常用的实现
Ansys FluentMRF + 界面无(前提为1:1)通过单元区域动作设置旋转
STAR-CCM+冻结转子界面轮廓缩放刚体运动指定
OpenFOAMMRF (定常)通过cyclicAMI进行界面处理
🧑‍🎓

Fluent没有节距比修正功能吗?


🎓

Fluent的MRF中界面两侧必须覆盖相同的角度范围。节距比不同的情况需要构建全周模型或整数比的扇形模型。这点CFX更灵活。


需要注意的限制条件

🧑‍🎓

冻结转子法有什么风险请教一下。


🎓

1. 尾迹冻结:叶片后流被固定在某个位置,对下游的影响可能过大或过小

2. 非轴对称进口条件:蜗壳下游的扩压管会产生非轴对称进口流,冻结转子中位置依赖

3. 转矩变化被忽视:实机中叶片与舌部的相对位置会导致转矩变化,冻结转子中被固定为单一值

4. 全周模型必要性:非轴对称系统中往往无法使用扇形模型


🧑‍🎓

所以冻结转子是"定常近似"这个事实必须时刻铭记。


🎓

正是这样。直接把冻结转子结果作为最终设计依据是危险的。必须用滑动网格进行验证,或至少通过多相位平均来评估不确定性。

咖啡小休 杂谈

涡轮机械专用CFD工具——NUMECA FINE/Turbo vs ANSYS CFX vs OpenFOAM的设置难度

在涡轮机械定常解析(冻结转子·MRF)中,主要3款工具按设置难度对比。NUMECA FINE/Turbo具有涡轮专用的GUI,通过向导式界面完成旋转方向·转速·转子-定子界面设置,初期设置最简便。网格品质指标和错误检查充实,能自动诊断初始收敛错误。ANSYS CFX的CFX-Pre向导具有几乎相同的便利性,在拥有Ansys许可证的企业中普及率高。OpenFOAM需要手动设置rotatingWallVelocity边界条件和MRFfvPatchField,需要对设置文件有理解。但一旦脚本化,多个参数变更也能自动化,适合优化循环。

冻结转子法的前沿研究

水力涡轮机中的冻结转子

🧑‍🎓

水力涡轮机的CFD中也用冻结转子法吗?


🎓

弗朗西斯水车和卡普兰水车的初期设计中经常使用。由导流叶→导叶→转轮→尾水管的4个域构成,转轮前后设置冻结转子或混合平面。


🧑‍🎓

听说水力涡轮机中尾水管会成为问题。


🎓

尾水管内会产生转轮出口余旋流形成的涡绳(涡核)。这是非定常现象,冻结转子无法捕捉。部分负荷下的涡绳引起的压力脉动评估必须用滑动网格。


IEC 62006指南

🧑‍🎓

水力涡轮机的CFD有规范吗?


🎓

IEC 62006是关于水力涡轮机CFD应用的指南。对网格品质湍流模型、边界条件的推荐事项有所规定。冻结转子进行的定常计算广泛用于生成效率希尔图。


水轮机的希尔图绘制

🧑‍🎓

什么是希尔图?


🎓

以单位流量 $Q_{11}$ vs 单位转速 $N_{11}$ 的平面上,用等高线显示效率的图表。


$$ Q_{11} = \frac{Q}{D^2 \sqrt{H}}, \quad N_{11} = \frac{ND}{\sqrt{H}} $$

通过改变导叶开度、转速、流量,用CFD计算100多个运行点,生成效率特性曲线。冻结转子的高速性在这里发挥了优势。


🧑‍🎓

100多个点! 没有冻结转子真的做不了。


🎓

是的。每个点用冻结转子仅需数小时,并行计算数天即可完成。但滑动网格每点需1天,100点就不现实了。

咖啡小休 杂谈

冻结转子法的精度限界——非线性干涉强的情况识别法

冻结转子(Frozen Rotor)法在定常MRF解析中计算很快,但有适用物理条件的限制。转子-定子间干涉强的条件——①转子-定子间距短(s/R<0.05)、②叶片负荷高产生强尾迹(Zweiffel数>0.8)、③再循环区跨界面的低流量·高负荷运行——会导致周向依赖性增大,解对"界面相对位置"变得敏感。定量判断标准是"在周向3~5个位置进行冻结转子计算,性能(效率·全压)的变动≥±2%则需要非定常滑动网格"这样的实务指南。先用冻结转子进行敏感性检查,必要时再切换到非定常——这个两阶段方法是最合理的选择。

冻结转子法的故障排除

冻结转子界面的收敛问题

🧑‍🎓

冻结转子计算中界面附近的残差下不来。什么原因?


🎓

有多个原因。


1. 节距比过大

🎓

节距比超过2时,插值导致的数值过冲容易发生。考虑用最大公约数构建扇形模型。


2. 界面两侧速度差大

🧑‍🎓

旋转域和静止域的速度差很大会有问题吗?


🎓

当然。从旋转坐标系相对速度向静止坐标系绝对速度变换时,大的不连续会恶化收敛。初期条件用前面运行点的解,或采用转速斜坡(逐步提升转速)有效。


3. GGI面网格不匹配

🎓

界面两侧网格大小相差过大会降低插值精度。旋转侧与静止侧界面附近的单元大小应大致相同(比例1:2以内)。


转速斜坡

🧑‍🎓

转速斜坡具体怎么做?


🎓

CFX中通过Expert Parameter设置,逐步增加转速。比如针对目标3000rpm,先用300rpm运行100迭代,再1000rpm,最后3000rpm。Fluent中在UDF内将区域动作转速定义为时间函数。


冻结转子→滑动网格的过渡

🧑‍🎓

冻结转子的结果能用作滑动网格的初值吗?


🎓

可以。CFX可以将冻结转子结果文件(.res)作为滑动网格计算的初始条件重启动。这样能大幅缩短非定常计算的启动时间。


🧑‍🎓

很方便呢。


🎓

设计流程来说,冻结转子快速概算→滑动网格精细化这个两阶段方法最有效率。

咖啡小休 杂谈

冻结转子CFD中压力比不稳定——界面插值失配问题

"转子-定子界面跨越处压力出现急跳"的问题源于界面流量守恒精度不足。特别在非一致界面(Non-Conformal Interface)中,转子侧与定子侧单元数·大小不同,面积加权平均进行的通量交换可能不守恒。诊断步骤:①界面两侧积分质量流量对比,差异≥1%则设置有误。②检查界面设置是否混合了混合平面与冻结转子的同时指定。③转子侧与定子侧界面附近网格大小不同时需重新细化。界面近物网格宽高比>10时是插值误差的主因,必须实施品质检查。

相关仿真工具

通过本领域的交互式仿真器体验理论

仿真工具列表

相关领域

热分析V&V·品质保证结构分析
本文评价
感谢您的反馈!
有参考
价值
需要
更详细
有错
有参考价值
0
需要更详细
0
有错误
0
由NovaSolver贡献者撰写
匿名工程师 & AI — 网站地图
查看个人资料