各位老师下午好，我用compressibleInterFoam计算横射流时总会出现Negative initial temperature的报错，这个问题困扰我很久啦，我也看了论坛内和cfd-online上各位老师的回答，但还是没有解决。
主要的尝试有：

采用可变的时间步长，对Co数进行限制（主要参考文献中的数值）
在主程序里面对温度进行限制，但这会导致下图所示非常异常的结果
在constant/fvOptions中添加limitTemperature限制
修改边界条件（目前p出口采用waveTransmissive边界条件，温度和速度出口采用inletOutlet边界条件）

此外，我对全场的温度最大值和最小值进行了监视，发现温度并不是逐步下降而是在某一迭代过程中骤降，感觉像是没有求解好一样

Time = 0.041271304
MAX::after TEqn_correct.................=max(T) [0 0 0 1 0 0 0] 1317.1702
   MIN::after TEqn_correct.................=min(T) [0 0 0 1 0 0 0] 299.99994
Time = 0.041271308
MAX::after TEqn_correct.................=max(T) [0 0 0 1 0 0 0] 2007.9887
   MIN::after TEqn_correct.................=min(T) [0 0 0 1 0 0 0] -4652.3157

并且在发散最后输出的时刻，就已经出现了异常的压力和速度（每个发散的case都是这样子，温度全场没有负温）：
case1：

case2：

能不能请各位老师给我一点建议，非常感谢！先谢过各位老师了！

CFD中文网

Rachel0096