滑移网格法
滑移网格法的理论基础
概述
在非定常涡轮机械解析中经常听到滑移网格法，与MRF有什么区别？
MRF将转子位置固定以求伪稳定解，而滑移网格法则在每个时间步长中物理上旋转旋转域的网格。通过非共形网格界面（Non-Conformal Interface）进行信息插值。
界面处的网格滑动，所以叫滑移网格法。
正是如此。由此可以物理上捕捉叶轮-叶轮干涉（尾流通过、势流干涉）。压力脉动和非定常叶面荷载预测是不可或缺的方法。
界面处理
界面的插值是如何进行的？
在每个时间步长中计算旋转侧和静止侧界面网格的重叠，并以保守方式插值通量。CFX的非定常转子-定子采用基于GGI的加权插值。Fluent的滑移网格采用基于面-面交叉判断的插值。
适用范围
在什么情况下需要使用滑移网格法？
| 分析目的 | 滑移网格法必要程度 |
|---|---|
| 设计点效率预测 | 不必要（MRF/混合平面就够了） |
| 叶片通过引起的压力脉动 | 必需 |
| 尾流-叶轮干涉 | 必需 |
| 非定常叶面荷载（振动评估） | 必需 |
| 噪声预测（FW-H输入） | 必需 |
| 浪涌/失速 | 必需（全周） |
滑移网格法的诞生——非定常转子-定子干涉CFD的黎明
滑移网格（Sliding Mesh）法的工业推广始于1990年代后期。在此之前，涡轮机械CFD仅限于定常冻结转子或混合平面法，无法捕捉转子与定子的非定常干涉（尾流干涉、势流干涉）。转折点是ANSYS Fluent在5.0版本（1998年）首次将"滑动界面"功能作为商用CFD软件正式实现。这使得转子到定子的涡通输，叶片通过频率（BPF）的脉动力计算和声学压力变化的预测成为可能。在现代涡轮机械开发中（喷气发动机压缩机失速解析、泵的水力脉动降低），非定常滑移网格已成为必不可少的工具，自1998年实现以来25年间改变了CFD的常规认识。
滑移网格法的数值计算方法
时间步的确定
滑移网格的时间刻度如何确定？
每个叶片通过时一般取20～50个时间步长。
例：3000 rpm、叶片12片、每个叶片通过30步 → $\Delta t = 60/(3000 \times 12 \times 30) = 55.6 \mu s$
需要更细的时间刻度时如何处理？
在噪声预测或DES/LES中需要解析BPF的10阶谐波，需要100～200步/叶片通过。
周期定常的判定
需要计算多少转数才够？
通过监测压力脉动来判断。
1. 在合适的点（如喘振附近）放置压力监测点
2. 比较连续两转的压力波形
3. 振幅变化在2%以内时认为达到周期定常
通常5～15转达到稳定。使用MRF或冻结转子稳定解作为初值可以缩短至3～5转。
节距比的处理
转子和定子的叶片数不是整数比时怎么办？
理想情况是建立全周模型，但计算成本巨大。可以按叶片数的最大公约数建立扇形模型。例如转子7片、定子12片时最大公约数为1，需要全周；转子6片、定子12片时最大公约数为6，可以建立1/6扇形模型。
CFX的时间变换法或FINE/Turbo的NLH法可以在非整数节距比条件下用单叶片计算近似非定常干涉。
滑移网格法的数值实现——AMI与通量修正的作用
滑移网格（Sliding Mesh）法通过"任意网格界面"（AMI: Arbitrary Mesh Interface）连接旋转的转子区和静止的定子区，每个时间步更新界面网格的相对位置来计算非定常转子-定子干涉。AMI处的变量插值采用双线性插值或加权最小二乘法（WLS），当界面单元非共形时需要"通量修正"以保证通量守恒。OpenFOAM中对AMI面的通量修正已自动化，但已知在两个转子转数后会出现质量守恒误差累积问题，建议定期进行额外的压力修正。
滑移网格法的实务应用
压力脉动的提取
如何从滑移网格结果中评估压力脉动？
按照以下步骤进行。
1. 在关注点（蜗壳壁面、管道连接部等）放置压力监测点
2. 采集达到周期定常后2～5转的数据
3. 使用FFT（快速傅里叶变换）进行频谱分析
4. 评估BPF及其谐波的峰值振幅
BPF振幅多大属于正常？
取决于机型，对于离心泵，BPF压力振幅通常为平均扬程的1～5%。超过这个值会增加管道振动和结构共振风险。
叶面非定常荷载
能否评估叶面非定常力？
可以。输出每个叶面受到的力的x、y、z分量时间历程，进行FFT频谱分析。当叶片固有频率与BPF谐波一致时存在共振（颤振）的风险。
后处理注意事项
滑移网格后处理有什么需要注意的？
洗衣机滚筒的CFD——用滑移网格法优化衣物的水分配
家电产品的意外CFD应用是洗衣机滚筒的水流解析。采用VOF法＋滑移网格法来解析洗涤液（含表面活性剂的水）的运动，以提高与衣物的接触效率。滚筒内设置的提升式挡板（baffle）的形状、数量、配置决定了洗涤液的掬起和衣物全体的均匀分配。松下的案例中，采用CFD（滑移网格+VOF+粒子追踪）与衣物污垢分散模型进行耦合解析，通过识别洗涤效率高的新型挡板形状，实现了较传统产品节水30%的洗衣机开发。这是制造业CFD应用作为高端产品差别化技术发挥作用的好例子。
滑移网格法的软件比较
求解器间比较
滑移网格的实现在求解器间是否有差异？
| 功能 | Ansys CFX | Ansys Fluent | STAR-CCM+ | OpenFOAM |
|---|---|---|---|---|
| 名称 | 非定常转子-定子 | 滑移网格 | 刚体运动 | cyclicAMI |
| 界面插值 | GGI（通量守恒） | 面交叉判断 | 内插/外插 | AMI（任意网格界面） |
| 节距比 | 剖面变换 | 整数比推荐 | 剖面缩放 | 整数比推荐 |
| 并行效率 | 优良 | 优良 | 优良 | 略低(AMI通信) |
| GPU支持 | 无 | 2024R1版起支持 | 部分支持 | 部分(AmgX) |
GPU支持进展如何？
Fluent 2024R1发布了GPU原生求解器，滑移网格非定常计算也能进行GPU加速。报告显示比CPU快5～10倍。这意味着原本需要1周的计算现在可以在1天内完成。
计算成本估计
滑移网格的计算成本大概是多少？
| 情景 | 网格数 | 核心数 | 转数 | 所需时间 |
|---|---|---|---|---|
| 单级泵(1节距) | 200万 | 32 | 10转 | 12～24小时 |
| 单级风机(全周) | 1000万 | 128 | 10转 | 2～4天 |
| 多级涡轮(1节距) | 500万 | 64 | 15转 | 3～5天 |
| 全周压缩机 | 5000万 | 512 | 20转 | 1～2周 |
全周计算确实需要很大的计算资源。
利用HPC集群或云计算（AWS、Azure等）是现实的选择。可以用Ansys Cloud或Simcenter Cloud Engineering在需要时申请大规模计算资源。
滑移网格CFD工具比较——ANSYS Fluent和OpenFOAM的AMI设置差异
对比ANSYS Fluent和OpenFOAM的滑移网格实现。在Fluent中，只需在Dynamic Mesh面板激活"Sliding Mesh"功能，并指定转子-定子间的Interface Zones，AMI会自动设置。转数的输入也可在Cell Zone Conditions的Frame Motion中直接设置，操作相当简洁。OpenFOAM需要在dynamicMeshDict和fvConstraints中手动设置AMI边界条件，在0目录的patchField中指定cyclicAMI。需要在多个设置文件间保证一致性。不过OpenFOAM在并行网格分割和AMI界面一致性方面曾有报告bug，v2112版起基本已解决。最新版本的确认和上下游patch对的设置一致性验证是稳定运用的关键。
滑移网格法的前沿研究
滑移网格法 + DES
什么时候将滑移网格法和DES结合使用？
需要同时解析非定常涡结构和叶轮-叶轮干涉时。典型应用包括风机噪声预测和涡轮冷却孔-主流干涉解析。
网格要求与URANS不同吗？
要求严格得多。DES/SDES在离主体较远区域会切换到LES模式，因此叶栅流道的主流区也需要LES适配的等向网格。
| 项目 | URANS | SDES/SAS |
|---|---|---|
| 叶面y+ | < 2 | < 1 |
| 叶栅网格尺寸 | 弦长的2～5% | 弦长的0.5～1% |
| 每个节距的网格数 | 50～100万 | 300～1000万 |
| 时间步 | 20～50/叶片通过 | 100～500/叶片通过 |
声学解析的连接
用滑移网格 + DES + FW-H进行噪声预测的流程是？
1. 用滑移网格+SDES/LES获取叶面和透过面（permeable surface）的压力时间历程
2. 累积最少10个BPF周期的数据
3. 输入至FW-H求解器，计算远场音压谱
4. 用A加权获评全面噪声等级（dBA）
最近的发展动向
滑移网格法最前沿的话题是什么？
滑移网格法的最前沿——数亿网格的大规模非定常LES解析
随着超级计算机的普及，滑移网格法与LES相结合进行大规模涡轮机械非定常解析正逐渐成为现实的研究手段。特别是在涡轮泵管系诱发振动（FIV）和直升机转子的声学预测中，RANS无法捕捉的非定常涡结构详细解析成为必需。在ECN的PROEL3D（欧洲水力涡轮研究）中进行了2亿网格级Kaplan水轮机的LES解析，成功将喷水管中非定常压力脉动（Vortex Rope）的RMS值预测值和实验值控制在±10%以内。计算需要3000核超级计算机100小时，但对叶片疲劳寿命评估的实用价值远大于计算成本。
滑移网格法的故障排除
典型故障
滑移网格计算常见的问题有哪些？
1. 初期发散
直接从滑移网格开始非定常计算容易发散，因为初期流场还未充分发展。用MRF或冻结转子的稳定解作为初值可以规避这一问题。
2. 时间步过大
库朗数大会导致发散吗？
CFX的耦合求解器隐式求解CFL限制宽松，但滑移网格界面的插值精度依赖时间步。界面处库朗数超过10时插值精度急剧下降。每个叶片通过角度应该在5度以内。
3. 界面网格密度不足
转子侧和定子侧界面附近网格密度若差异大，插值误差会在非定常计算中累积，导致振荡或漂移。应将界面前后的网格尺寸控制在1:2以内。
结果验证清单
如何验证滑移网格的结果？
如果时间平均效率与MRF相差较大怎么办？
差异超过2个百分点时，可能是时间步不足、网格质量不良或初期过渡切断不充分。应将时间平均的起始点足够延后（从第5转开始平均等）后重新评估。
滑移网格CFD中非定常力振荡——时间步与AMI插值的关系
"转子上的力随时间步大幅振荡，无法进行收敛判定"是滑移网格CFD中的常见问题，通常由时间步与AMI插值的组合引起。当转子在一个时间步内旋转角度过大时，AMI界面单元相对位置变化大，插值误差增加。指导准则是：确保叶片通过角度内的时间步数至少为20步（从叶片数Z和转数N计算出Dt_max = 1/(20*N*Z)）。另一对策是将转子-定子力按"叶片通过周期（BPP: Blade Passing Period）"进行平均，可以消除短周期数值振荡。进一步优化AMI插值的阶数（一阶 vs 二阶）也能降低振荡，但会增加计算成本——应在确认收敛性后再选用。
相关话题
更多
错误