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

[Alvin]

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

[xpqiu]

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+太小，结果反而不好，也可能会引起计算不稳定。

[heike256]

感谢！

[一颗鸭蛋]

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

[cccrrryyy]

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

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

[一颗鸭蛋]

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

[bit_hypersonic]

@xpqiu 老师好，我按照您方法的第二条调试成功了，我比较好奇的是，这样的原理是什么？ cellimited是什么含义呢？

[xpqiu]

@bit_hypersonic
cellLimited 这个东西是一种 gradient limiter，你可以搜一下这个关键字了解一下其物理含义。也可以参考如下论文看看：Michalak, K., & Ollivier-Gooch, C. (2008). Limiters for unstructured higher-order accurate solutions of the euler equations. 46th AIAA Aerospace Sciences Meeting and Exhibit, January. https://doi.org/10.2514/6.2008-776

[bit_hypersonic]

@xpqiu 好的，谢谢老师