realizable k-e模型初始化

问题简介:400m*200m的计算域，轴对称，底部0.14m的圆孔入射高温309℃的流体，0.2米每秒。
模型中密度等属性随温度变化，温度变化范围2-309摄氏度，密度变化范围700-1000摄氏度。采用relizable k-e模型，在初始化时k和e初始化为零时，没算几步湍流粘度出现非正常的极值点，并很快发散？如果把密度换成常数，则比较正常。请问问题出在哪里？