--> FOAM FATAL ERROR:
Negative initial temperature T0: -283.834

我尝试把温度边界全都设置为同一个值或zeroGradient（我理解等同于绝热边界），这时又会出现Tnew过高（达到4.6e+15）的报错，如下图(iter 1以后的数值与iter 1基本一样)。

iter    Test         e/h     Cv/p    Tnew
0    1880.54     1.12e+6     712    4.61e+15
1    4.61e+15    3.28e+18    712    4.61e+15

我发现计算时连续性残差也一路飙升，觉得可能是流动与传热耦合出问题，流动和传热主要通过密度rho耦合。于是我把密度设为定值rhoConst，试图将传热和流动解耦，但计算仍会发散，提示负温度，且连续性残差仍一路飙升。

--> FOAM FATAL ERROR:
Negative initial temperature T0: -905.973

现在就很纳闷，既然之前能使用simpleFoam求解器求解流场，那应该流动的边界条件设置没问题，温度的边界条件又是绝热边界，也没有热源，温度场应该就是设定的均匀的初始温度，可为什么计算却发散了
[案例文件](链接：https://pan.baidu.com/s/13hw5KDph_17AmRPrmVlIhQ
提取码：OFv8)