K-e湍流模型(2)

湍流

合求解粘性流场），发现湍动能产生项（雷诺应力和一个速度张量相乘组成的项）出现负

值，请问是不是一种错误现象？

如果是错误现象一般怎样避免。另外处理湍动能产生项采

用什么样的差分格式最好。而且因为源项的影响，使得程序总是不稳定，造成k,e值出现负

值，请问有什么办法克服这种现象。

你可以试试这里计算的时候加一个判断，出现负值的时候强制为一个很小的正值。

这可能是因为你采用的数值格式的问题，一般计算程序对k方程都要做一定处理，

以保证k的正定。

比如，强制规定源项与0的关系，以使数值计算稳定。

就ke模型而言。

它是problem dependent.对简单的无弯曲无旋转无...的湍流问题，它能算而且能给出好的结果，但对复杂的流动问题，它就不能使用了。

出现负的ke不仅仅是计算格式的问题，

更重要的是模型问题，没有谁能证明ke模型在任何流动问题中都能保证ke是正的。

有这么一些办法避免ke出现负值

1。对K＝ln(k)和E=ln(e)求解，问题：壁面ke=0难处理，

2。先用层流计算500步，然后再用ke算

3。各种强制限制办法

4。源项局部线性化

5。算到一定程度，如果k值趋势对了，就干脆不求ke方程