吾生有涯 学海无涯
析模有界 知识无界

MRF|滑移网格|重叠网格

上一篇文章(SRF、壁面运动及动网格)中探讨了SRF与动网格的区别。本文在前一篇文章的基础上进一步对比并探讨MRF、SMM以及重叠网格的差异。

还是上篇文章中的几何,如下图所示。密闭容器中包含一个以角速度200 RPM的十字桨叶。本次比较采用MRF、滑移网格以及重叠网格时计算结果的差异。

注意,重叠网格采用的方案与MRF及滑移网格不同,其需要准备两套计算网格。

1 计算网格

MRF及滑移网格计算采用的网格如下图所示。

需要注意的是,滑移网格必须使用交界面,否则会提示出现负体积。MRF模型可以使用内部面也可以使用交界面(推荐使用内部面)。为了便于比较,全部使用interface。交界面的位置取距离桨叶最远端附近(其实很容易理解,交界面取得更远等同于使用SRF)。

注:不过交界面会造成较大误差,也不宜放置在桨叶前端附近,这里只是测试分界面距离运动边界过近会有什么后果。通常情况下将分界面放置在桨叶前端与壁面的中间位置,尽量避免分界面对桨叶前端及壁面造成扰动,但是也要尽可能减少将静止区域人为指定为动区域。

为确保产生interface,在SCDM中不要对几何进行共享拓扑。将内部区域命名为inner,外部区域命名为outer。为了与前文动网格结果进行比较,这里也采用三角形网格划分,并采用相同的网格尺寸。生成的计算网格如下图所示。

2 MRF计算

  • 稳态数据
  • 采用RNG K-Epsilon湍流模型
  • 采用MRF模型,指定内部区域角速度为300 RPM
  • 创建分界面
  • 设置桨叶壁面速度为相对区域速度为0
  • 计算得到的速度分布(残差收敛到1e-7)

3 滑移网格

滑移网格需要计算瞬态,这里可以在MRF基础上进行修改计算。

  • 选择Transient进行瞬态计算
  • 指定inner区域运动为Mesh Motion,并指定其旋转速度为300 RPM
  • 速度动画如下图所示

4 重叠网格

重叠网格的几何与前面MRF及滑移网格不同,其需要两个重叠在一起的独立几何,并分别划分计算网格。

为了便于比较,这里采用了三角形网格,网格尺寸与前面保持一致。生成计算网格如下图所示。

  • 激活Transient采用瞬态计算
  • 指定前景网格区域inner旋转角速度为300 rpm
  • 定义重叠网格区域
  • 计算结果如下图所示

注:重叠网格开启了并行计算就出错,用1个CPU就没问题。现在不清楚是软件安装有问题还是版本有问题。

5 总结

下面是相同的工况条件下分别使用SRF、MRF、SMM、动网格以及滑移网格得到的结果。放在一起便于比较。

  • SRF计算结果
  • MRF计算结果
  • 滑移网格计算结果

  • 动网格计算结果

  • 重叠网格计算结果

注:上面的动图中桨叶旋转快慢不同是由于截图采用的帧率不同所造成。

以动网格结果作为目标(毕竟与真实工况最吻合),从上面的计算结果来看:

  • 对于文中所示的旋转运动,上面所有的方法其实都是可以使用的
  • 稳态计算时MRF要好于SRF,原因在于SRF将更多原本为静止的区域处理成了运动区域,而MRF可以缩小这些区域的范围
  • 瞬态计算时优先选用滑移网格,能最大程度保证网格质量。不过滑移网格需要使用interface,在计算时依然存在插值
  • 重叠网格计算精度并不高,因为重叠域内需要不断插值,不过好在重叠网格不需要担心网格运动带来的负体积问题
  • 动网格是与实际工况最为符合的,不过动网格方法在运动时网格质量会变得很差

相关文件:

提取码:8960

本篇文章来源于微信公众号: CFD之道

赞(4) 打赏
版权声明:未经允许,请勿随意用于商业用途。
文章名称:《MRF|滑移网格|重叠网格》
文章链接:https://www.topcfd.cn/17517/
本站资源仅供个人学习交流,请于下载后24小时内删除,不允许用于商业用途,否则法律问题自行承担。
分享到

说两句 抢沙发

评论前必须登录!

 

觉得文章有用就打赏一下文章作者吧

非常感谢你的打赏,我们将继续给力更多优质内容,让我们一起创建更加美好的网络世界!

支付宝扫一扫

微信扫一扫

登录

找回密码

注册