多相流超声速计算出现负温度
-
-
用CompressibleIsoInterFoam算超声速射流有两个本质上无法解决的问题：
- 压力基Rhie-Chow插值属于中心插值，可压缩出现shock时会振荡，如果开启了choi修正加入了耗散会把激波抹平（当然看你关注什么，如果只是关注推力之类的统计量可能没啥影响，但是研究流场结构就够呛了）
- 如果你认真推一遍公式，会发现它的能量（温度）方程压根就是错的，算出负温度不足为奇。另外，我记得thermo里的mixture函数在每个时间步计算时都会initial一下，温度会重新算，可以按“Negative initial”关键字搜索定位一下代码。
总之吧，这个求解器勉强能用，如果你换套网格、换个格式（比如snGrad用limited）或许到这个时刻就能算通了，然而隐患一直都在，不知啥时候又会出现。
如果真想把浸没射流算好，建议用多相的density-based solver，我们算的效果长这样：
-
