算例发散问题

[先先贝]

以前台阶层流算例为例，采用hePsiThermo<multiComponentMixture<const<hConst<perfectGas<specie>>,sensibleEnthalpy>>>模型，入口边界条件设定为1Pa、1K、3m/s时，能稳定求解运行；当采用thermoType hePsiThermo<multiComponentMixture<sutherland<janaf<perfectGas<specie>>,sensibleEnthalpy>>>模型，入口边界条件调到普通超声速入口边界条件数值时运行就出现以下报错。求解大佬们指点，该怎么解决呢! [替代文字]

[wangfei9088]

@先先贝 看上去是#4 sqrt 开方函数有问题，是不是出现负数了？

[先先贝]

这个地方计算了声速，但是输出了一下值，也没有负数欸：

[先先贝]

[wangfei9088]

@先先贝 那瞧不出来啥问题，得慢慢debug了。

[先先贝]

@wangfei9088 谢谢！只能慢慢试了