bounding k,bounding epsilon,均超限,连续性方程不收敛



  • 流场湍流模型中k,epsilon值、边界条件以及差分格式等方面如何选取,才能尽量避免出现bouding k,bounding epsilon值超限而导致的流场发散的错误?大家有什么好建议么?


  • OpenFOAM教授

    1. 网格质量如何?如果网格太差(checkMesh有严重报错,最大非正交大于85,skewness 超过20,等),建议回炉重造。
    2. k 和 epsilon 的对流项离散格式用 linearUpwind,grad(k) 和 grad(epsilon) 用 cellimited,如果还不行,直接上 upwind试试。也可以先用upwind算一段时间再换成 linearUpwind。
    3. 边界条件的设置取决于你使用的湍流模型,对于高雷诺数模型,入口 k 和 epsilon 给固定值,根据湍流度来估算;壁面上 k 用 kqrWallFunction,epsilon 用 epsilonWallFunction,nut 用 nutUSpaldingWallFunction。高雷诺数k-epsilon湍流模型还要注意壁面 y+,一般认为需要 y+>30,如果y+太小,结果反而不好,也可能会引起计算不稳定。


  • 感谢!



  • @xpqiu 老师您好,我在求解器中添加瞬态温度方程,浮升力加到动量方程中,使用Komega模型算矩形管中的磁流体传热,也是一直bounding k,bounding omega,K,omega的初始条件怎么设都发散,是K,Omega初始条件太敏感了嘛,不断在试错当中,感觉拼运气,请问有什么好的建议嘛:chouchou:


  • OpenFOAM副教授

    @一颗鸭蛋 是只有kOmega才这样还是别的模型也这样?RANS的话入口条件都挺标准的,好像没有什么可以大改的地方,参考这里

    RANS对壁面要求可能稍微严格一点,倒不是说需要很细密,而是不同的模型和壁面处理适用于不同的网格,FOAM里面壁面函数可以参考这里。不过听起来你的不像是算不对,而是直接发散。建议先换个模型试试,确定到底是不是kOmega的问题。



  • @cccrrryyy 我把模型改成了KomegaSST了 然后直接可以计算了,还挺快的。所以我现在就觉得KOmega模型对入口条件太敏感了,其次可能是我求解器里添加了浮力项后更不稳定了。我再按您的推荐测试一下Komega。


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2