本专题前2期内容如下所示,建议从第1期开始阅读:
关于水力设计及网格划分部分,可以查看上述2篇专题,今天来看看如何在Fluent中采用MRF及SMM模型进行离心泵性能模拟,得到基本的扬程、效率、压力场、速度场等参数。
文中涉及到很多内容之前都已经写过了,内嵌了相关的文章链接,可以点击查看详细说明!
首先将划分好并装配好的网格导入Fluent中,最好是能够装配完成再导入,能省去不少报错和尺寸不匹配的麻烦!
一、稳态求解
- scale网格尺寸
导入的模型很有可能尺寸不是很合适,可能需要缩放网格,所以导入网格后先scale,将显示单位修改为mm后看尺寸是否和自己的模型一致,如果不在一个数量级则需要缩放,在后边的“Mesh Was Created In”中选择对应要转换的目标单位,点击下方Scale即可,直到确认尺寸合适!
- 勾选“Gravity”
介质为水,需要考虑重力,勾选“Gravity”选项并设置合适的方位及重力加速度的值!此处为-x方向,所以x方向重力加速度设置为-9.81
- 选择合适的湍流模型
关于湍流模型的选择,可以查看下胡老师总结的选择方法
选用原则:
- 通常情况下建议选择Realizable k-epsilon模型及SST k-omega模型
- 当需要高精度求解边界层,如涉及流动分离或传热细节时,建议使用SST komega模型
- 当仅需要粗略估计湍流时,可以使用Standard k-epsilon模型
此处我们选择标准k-epsilon模型(Standard k-epsilon)
- 增加介质“water_liquid”
添加完成后会出现在介质中
- 定义计算域
计算域的定义包括计算域介质的选择以及旋转属性的定义,可以查看下方文档:
MRF(多参考系)旋转模型
- 定义边界条件
边界包括进口,出口,旋转壁面,静止壁面,关于进出口应该怎么选择边界,可以查看下方的文档:
边界条件
对于进口INLET,已知流量和进口段横截面积,可以计算出流速,所以采用速度进口(Velocity inlet)
对于出口Outlet,已知设计扬程,可以计算出出口总压,所以可以定义为压力出口(Pressure outlet)
对于旋转壁面,定义如下:
- 定义进口压力、出口压力、扬程等参数的监控
这个之前已经写过详细的教程了,可以查看下方文档:
- 初始化
初始化采用标准初始化,从进口INLET初始化即可
- 求解计算
稳态求解,给定总步数为1200步,计算完成!
通过查看监测的扬程,泵的扬程在42m左右,通过查看监测的详细数据,H=42.5m,这比我们的额定扬程高出了5m,通常模拟出来的结果都是要比实际高一点的。
- 后处理
采用CFD-POST来看一下结果
叶片压力
叶轮压力
速度云图
速度矢量图
流线图
- 稳态模拟结束
稳态求解是与时间没有关系的,如果要看瞬态叶轮的旋转情况,咱们可以采用SMM模型也就是滑移网格来在稳态的基础上进行瞬态的求解。
二、瞬态求解
- 瞬态滑移网格
关于SMM模型的使用方法及设置要点可以查看下方的文档:
SMM(滑移网格)旋转模型
- 修改求解方式为Transient
在稳态计算结果的基础上直接修改,无需重新初始化,稳态的结果将自动作为瞬态的初始值
- 将叶轮域的MFR模型copy到SMM模型
其他部分无需改动!
- 设置监控界面
比如我要直观的感受叶轮的旋转,并创建动画,那么就需要用到这个功能了,之前有写过详细的教程,可查看下文:
- 设置步长及步数
这个也写过很多很多遍了,如何确定步长时间以及步数
在计算瞬态的时候需要设计时间步长和步数。在有关书籍中的介绍是叶轮转过6-8个周期就可以判断流场已经趋于稳定了,根据我们设定的转速为1450r/min,也就是24.167r/s,这样算下来叶轮转一圈需要的时间是1/24.167=0.04137s,一圈是360°,那么如果我们给定叶轮每转动1次转过的角度是2°,可以根据0.04137/360*2算出时间步长为2.3e-4s,步数为0.04137*8/2.3e-4s=1440步!
所以根据以上计算结果,设定时间步为2.3e-04。设定步数为1440
但是真实的步数应该是1440*20=28000步,大家如果仔细看Fluent的求解设置界面的话,会看到这样一个参数“Max Iterations/Time Step”,默认数值是20,意思是在一个时间步长内还会迭代计算20次,也就是实际计算20步才是一个完整的步长。所以实际计算步数应该是1440*20=28800步!
根据我电脑1秒1步的速度,算完大概需要8小时,豁!
通过瞬态计算能够得到更精确的扬程的数值,瞬态求解也更接近物理实际情况,但是计算量也要比MRF模型大很多,所以群里老哥讲一般只跑稳态也是ok的。这玩意现在分析已经到了人人都会的阶段,真正应该学习的是叶轮的优化设计,学会优化设计收入蹭蹭蹭!哈哈
- 开始计算
瞬态求解的残差曲线和稳态是截然大不相同的,看看胡老师针对这俩的表述吧!
文章链接:点击直达
文中不仅对残差这个概念进行了解释,还讲了如何判断是否收敛,收敛的标准以及残差曲线的类型,看完这篇文档秒懂!yyds
-
- 叶轮动画展示
- 泵后处理
这个在年初做轴流泵案例的时候详细写过文章,包括如何获取扭矩、压力、计算扬程、功率、效率等等,可查看下方文章:
具体的数值我就不算了,感兴趣的同学可以根据上述文档自己算一下,看看效率和设计效率差多少!
- 结束
三、结语
这一年来Fluent用的最多,忽略了CFX的学习,下一期专题更新CFX求解的教程,今天上手发现很多东西都忘了。
免责声明
泵小丫收集资源均来自网络
仅用于内部交流之用,如有需要请支持正版
非常好的资源?,推荐给其他泵开发爱好者
我看行!?
能问一下fluent里在哪儿设置的叶轮每次转过的角度是2°呢,可以改成1°或者更高吗,求管理员 解惑
你好,模型导进去这样是怎么了
这不是很正常么,显示一下网格display
网格显示还是这样
你可以改一下网格显示颜色的
这个在哪里改啊
Display旁边的Colors
导入mesh文件进行check以后报错,怎么解决呢
先创建interface再check
感谢?
又出现了报错
从入门到放弃
从一个坑再到另一个坑
[…] 离心泵模拟专题3:Fluent离心泵性能模拟 […]
请问,旋转壁面是哪三块呀
前盖板shroud、后盖板hub以及叶片blade
为啥不给下载呀
哪个不能下载啊
监控扬程998.2/9.81还是998.2*9.81
H=(P2-P1)/ρ/g
老哥,我下载的你的网格,出现下面的问题,不知怎么解决?interface zones overlap for mesh interface
网格导入fluent后需要重新对边界属性进行设置
对着的对着的
如何获取监测点的压力脉动值呀
网站有相应的教程
在初始化的时候,出现More than one cell is equidistant from the speciied reference pressure location 是为啥呀?是因为网格的原因嘛?求丫总解答。
只要能计算,这些都不是问题
这个是在哪打开
初始化的时候这些选项是灰的没法编辑 是哪里的问题
Failed to correct face handedness of 2 out of 10 left-handed faces on the sliding interface zone 8. 请问博主:在初始化的时候出现这个报错,是什么原因呢,如何去解决这个问题呢?
这个问题有解决方案,我找找
好勒,期待老哥的回复,感谢!
好大哥,找到解决方案了吗?还是卡在这里了,麻烦啦。
解决办法: 1、mesh/repair-improve/repair-face-handedness; disable repair face handeness ? [Yes] 2、Domain/Check/performcheck mesh; 3、mesh/repair-improve/repair; 4、检查网格质量,左手面网格消除; 5、“mesh/repair-improve/allow-repair-boundaries”重复以上操作; 再次进行“Check‘”,问题解决;
哥,根据您的过程操作一遍,问题还是存在
首先感谢老哥的回复。
老哥问题解决了,用的不是您的方法,但是还是谢谢老哥的指导意见,我把所有网格文件全部重新画了,发现问题是在icem设置网格大小的选项中,勾选了scale factor和min size limit 这样导致网格边界并不顺滑,具体原因我就不清楚了,取消勾选之后,问题就解决了。简单反馈一下~
恭喜哇
老哥您好,我遇到和您一样的问题。请问您用的是哪个版本的icem,我用的是19.2,上述scale factor和min size limit的设置是在全局网格尺寸那里吗?第一scale factor一打开就是亮的哈,怎么做到不勾选,下面的那个min size limit我也没选过,还烦请老哥不吝赐教~!
我遇到了这个问题,是因为在交界面法向上的加密导致的,取消加密后就消除了左手面
能请问一下具体怎么操作吗,我的网格导入后也出现了左手面问题,弄了好久😭
如图所示😎
请问您解决这个问题了吗,我也遇到了
能问一下fluent里在哪儿设置的叶轮每次转过的角度是2°呢,可以改成1°或者更高吗,求解惑
可以
为啥我的计算结果感觉失真严重阿
看数据明显不合理嘛
请问胡老师对湍流模型的选择总结有参考文献吗
在一本书上《ANSYS CFD疑难问题实例详解》
丫总,请问出口边界条件为啥是500000pa,设计扬程不是37.3m吗
👍
👍
请问老师,我在计算这个模型的空化时遇到了回流现象,请问是什么原因造成的?该怎么解决?非常感谢