冻结转子法
冻结转子法的理论基础
概要
冻结转子法这个名字很酷,但它具体是做什么的呢?
旋转体与静止体的相对位置固定进行定常解析的手法。动叶的坐标置于旋转系中,静叶的坐标置于静止系中,但在界面处不进行周向平均,而是"直接"传递信息。
与混合平面法的区别是什么?
| 手法 | 界面处理 | 叶片数约束 | 计算成本 | 精度 |
|---|---|---|---|---|
| 冻结转子 | 位置固定、直接插值 | 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 Fluent | MRF + 界面 | 无(前提为1:1) | 通过单元区域动作设置旋转 |
| STAR-CCM+ | 冻结转子界面 | 轮廓缩放 | 刚体运动指定 |
| OpenFOAM | MRF (定常) | 无 | 通过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有规范吗?
水轮机的希尔图绘制
什么是希尔图?
以单位流量 $Q_{11}$ vs 单位转速 $N_{11}$ 的平面上,用等高线显示效率的图表。
通过改变导叶开度、转速、流量,用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时是插值误差的主因,必须实施品质检查。
相关专题
价值
更详细
误