buoyantSimpleFoam求解器使用求助

[文刀白勺]

本人之前使用simpleFoam求解器计算钢液流场时可以很好收敛，但现在采用buoyantSimpleFoam求解器计算考虑浮升力（Boussinesq）的钢液流场时却发散，发散提示是出现负温度。

--> 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)

[李东岳]

这个算例是参考的某sci么？为什么128m这么大啊

[文刀白勺]

@李东岳 老师好，这个算例是针对工程实际的，不是参考文献里面的

[李东岳]

我下来看看

[文刀白勺]

谢谢老师。
我现在将k和omega的初始场的值设为0.1和10，计算就没有发散了，虽然计算的收敛速度较慢，但把松弛因子调小后也还是在不断减小。
fluent里面初始化k和omega场时也是按入口计算的值，为什么到OpenFoam里面这样初始化会导致收敛性变差呢？