OpenFOAM里kOmega模型的问题
-
@东岳 老师我把壁面改成omegaWallFunction又跑了一遍,结果改善的不是很明显,我猜是壁面的计算对喷流流场计算影响有限。
然后我去对比了Fluent中和 OpenFOAM中kOmega方程的区别,其实之前也比较过,发现其实大部分系数都一致,但有两个和流动相关的系数在OpenFOAM里做了简化。具体是这样的:
标准kOmega方程:
其他项和系数都一样,但是OpenFOAM里这两项
的系数和Fluent中存在差别。OpenFOAM把
分别定义为beta和Cmu,取值分别为0.072和0.09.也就是说忽略了 的影响,将其默认为1.
fluent中的这两项有更细节的计算方法:
由于我也找不到什么别的不同了,所以我就在OpenFOAM的kOmega模型里加上了这两项的计算,又重新跑了一遍case。这回跑出来的结果倒是和Fluent里的计算结果比较接近了,几次对比的结果如下:
不过我也不敢确定我这样就是对的,或许这两项系数的修改对于我的算例影响很大,对于其他类型的流动影响就不是很大。所以也想看看大家有没有用OpenFOAM里的kOmega模型算过其他的算例,或许对于有些算例使用OpenFOAM默认的kOmega模型就可以算的很好,不需要考虑那两项系数的变化,对于有些算例需要考虑系数的变化,希望能确定一下OpenFOAM里标准kOmega模型的适用范围。
我把我改的模型放上来,如果老师同学有兴趣可以用这个模型算一下试试。不过改的时候代码编得十分粗糙,就简单实现了功能,运算优化之类的都没有考虑。如果有人愿意帮我检查检查或者看看能不能优化就更好了,谢谢大家~mykOmega.zip -
omega模型适合计算高速流的问题吗?
-
@huoyanxiadebing 哈哈,其实也能看出来用kOmega模型算我这个算例结果和实验差别挺大的,算的不准,或许是不适合这种高速流动。
做这个对比其实是想说openfoam里的kOmega模型和fluent里的标准kOmega模型不太一样,和我在书上看到的也不太一样,主要是有的系数做了简化。
这种变化对我的算例的计算结果有较大影响,但对于其他的算例或者流动情况这种简化的影响可能就不会很明显。而且这种简化对计算的影响是积极的还是消极的也不好判断,对于有些流动可能简化了算的更准,对于有些流动可能不简化算的更准,这可以提醒我们以后用OpenFoam的kOmega模型算东西可能要考虑一下这种影响。万一也有朋友想对比一下fluent和openfoam的计算结果,也用了komega模型,也发现差别很大,那我们做的东西还能相互映证一下,没准还更有底气,还能少走个弯路。
然后kOmega模型的使用范围,我在网上查到的是低雷诺数预测效果较好,但好像高雷诺数也还行。对于自由剪切流比较适合。不过这些概括都比较笼统,我也就只算过一个喷流,没啥工程经验,如果有新算例我也就只能按着网上的建议去试,然后合不合适,有的时候算了才知道,还得具体问题具体分析。我觉得湍流模型这块其实挺有意思的,以我的算例为例,为了算喷流好多人拿不同的湍流模型在那测试,各有各的优劣。不过我觉得我还得多学习学习,吸取前人的经验。 -
@shampoo 在 OpenFOAM里kOmega模型的问题 中说:
明显
从理论上来说,komega模型不适合高速流流动,你如果单纯的研究边界层或剪切流中这种低速流动,可以用komega模型进行计算,但是你计算的喷管流域中既包括高速流(湍流核心区)也包括低速流(边界层或剪切层),低速区和高速区的湍动粘度系数求解是不相同的,理应使用SST komega。另外,试湍流模型的方法非常不可取,在你计算之前你就应该知道所使用的湍流模型是否适合你的计算域,是直接求解边界层中的流动还是使用壁面函数的经验模型对这种分离流很大的数值结果肯定有影响。
-
@huoyanxiadebing 哇您解释的好清楚,学到了很多,感谢感谢!我积累的实在是太少了,还得加油呢