Fluent的计算结果与理论公式值不匹配
-
求助各位大佬,我正在模拟光滑通道的流动以复现一篇文献,但是计算结果比理论公式值低了约15%
文献中使用相同设置能得出正确结果,网格无关性已验证,增加到2000W网格也无法得出正确结果。我想知道问题出在哪里?以下是我的模拟细节:
那篇文献的网址:https://doi.org/10.1016/j.ijheatmasstransfer.2021.121573 -
@李东岳 补充下具体细节:
网格:X方向流向,总长5m,网格长约6mm,均匀分布
Y方向展向,总长0.5m,网格长约3mm,均匀分布
Z方向高度方向,总长0.125m,底面和顶面第一层网格高度0.3mm,增长率1.1,网格长度最大值3mm,根据算出的结果完全发展后的y+值约为0.8
方法:SIMPLEC,全二阶迎风格式
物性参数:
320K空气的参数,文献里没给物性常数,但是可以根据h+的公式算出的√ρ/μ的值,进而推断出文献里貌似用的是默认参数,改成默认参数之后压降还是不对而且应力更小了
收敛判定:残差能量1e-9其余1e-5
-
@李东岳 感谢大佬。网格实际上是因为想获得完全发展湍流,把通道弄很长所以才会有这么多。一个光滑通道的网格也不怕泄密啥的。
网格链接(两个网址是同一个文件):
https://pan.baidu.com/s/1qTxoeUEgPGzIB2kSDf0YVA?pwd=1ruw
https://www.aliyundrive.com/s/XHWpd4c8FfU 提取码: 2n2l网格是用icem画的,导出的solver里没找到openfoam所以用的是fluent
-
@李东岳 project1.blk project1.tin
emmmmmmm我不太了解没有涡指的是什么,这一步也只是为了下一步大涡模拟获得完全发展的速度入口,如果没有影响的话倒也无所谓。
另外,刚刚发现我模拟出的结果好像是层流(也可能是顶端速度变化很小的湍流?),Re=20000的流动应该是湍流,这也许就是问题所在?
-
@李东岳 研一,这几个月才开始接触CFD,目前属于赶鸭子上架,有些概念确实来不及深入了解 可能之前我话没说清楚有点歧义,我知道光滑通道的RANS模拟不会有涡,但是我在帖子里也没提到涡所以我不太清楚大佬你提涡是为了什么。RANS模拟不会有湍流的那种状态我也清楚,我主要想表达的是速度分布的状态是官方教程里左边图的样子,但是从雷诺数来看应该是右边的样子所以会不会是这里的问题,关于这方面的知识我确实需要深入了解一下。
一维网格的概念没接触过,最开始倒是有考虑过二维网格,不过因为涉及到雷诺数参考长度的选取以及profile文件导入三维网格可能出现的问题就转做三维了。模型搞这么长主要还是为了保险,毕竟官方教程里的湍流公式按水力直径算得4.6m才能完全发展。目前对我来说网格多倒不是问题,反正还得上课,课上完了差不多也算完了。
现在我最想深入了解的还是条件一样为啥算出来的东西不一样而且结果还是错的的问题