rhoCentralFoam发散，出现Negative Initial Temperature

[RebelYoung]

@东岳 我其他的边界条件pT都是一样的，控制雷诺数，计算了不同的速度。因为是层流，所以上传的是最小速度1m/s的。

[李东岳]

不能用rhoCentralFoam算低音速流动，这个速度太低了

[RebelYoung]

@东岳 那么高速导致雷诺数量级太大，层流模型还能算么？

[RebelYoung]

https://www.cfd-china.com/topic/1737/可压缩求解器为什么密度残差一直为0-迭代步也一直是0

这里面回答了rho，rhoE，rhoU残差为0的原因。

[RebelYoung]

@东岳 重新试了一下250m/s，相当于0.83Ma，算了20步就出现负初始温度的问题

[李东岳]

@RebelYoung 网格是否可以使用六面体网格？

[RebelYoung]

@东岳 我是用pointwise画了导出的，这个球比较难画 ，
我的另一个简单算例cube (提取码: ywvc)是有附面层的全结构网格，上传的设置是初始条件是类似rhoS的0012算例的freestream type。
来流速度30m/s，就是上面贴图的龟派气功的结果，在几百步后也会发散，出现负温度。

[RebelYoung]

250m/s的cube，60步出现negative initial temperature。

# Field minima and maxima
# Time        	field         	min           	location(min) 	max           	location(max) 
3.57143e-05   	T             	1.454278e+02	(1.646393e-02 1.393888e-19 -1.709886e-18)	3.108947e+02	(-1.646393e-02 1.393888e-19 -1.725321e-18)
3.57143e-05   	e             	1.043157e+05	(1.646393e-02 1.393888e-19 -1.709886e-18)	2.230054e+05	(-1.646393e-02 1.393888e-19 -1.725321e-18)
3.64286e-05   	T             	1.335562e+02	(5.517241e-03 1.105814e-02 1.105814e-02)	3.098790e+02	(-1.646393e-02 1.393888e-19 -1.725321e-18)
3.64286e-05   	e             	9.580017e+04	(5.517241e-03 1.105814e-02 1.105814e-02)	2.222768e+05	(-1.646393e-02 1.393888e-19 -1.725321e-18)
3.71429e-05   	T             	2.917792e+01	(6.206897e-03 1.105814e-02 1.105814e-02)	3.078215e+02	(-1.646393e-02 1.393888e-19 -1.725321e-18)
3.71429e-05   	e             	2.092939e+04	(6.206897e-03 1.105814e-02 1.105814e-02)	2.208010e+05	(-1.646393e-02 1.393888e-19 -1.725321e-18)
3.58333e-05   	T             	1.400088e+02	(5.517241e-03 1.105814e-02 1.179091e-02)	3.108729e+02	(-1.646393e-02 1.393888e-19 -1.725321e-18)
3.58333e-05   	e             	1.004286e+05	(5.517241e-03 1.105814e-02 1.179091e-02)	2.229898e+05	(-1.646393e-02 1.393888e-19 -1.725321e-18)
3.65278e-05   	T             	6.218740e+01	(5.517241e-03 1.105814e-02 1.105814e-02)	3.097428e+02	(-1.646393e-02 1.393888e-19 -1.725321e-18)
3.65278e-05   	e             	4.460715e+04	(5.517241e-03 1.105814e-02 1.105814e-02)	2.221792e+05	(-1.646393e-02 1.393888e-19 -1.725321e-18)

最后几步的温度最值，最后两步$T_{min}$下降一个量级。

[RebelYoung]

沉了么。。？有大佬试试我的算例了么？

[李东岳]

周末我下载你结构网格的试试

[Tens]

不知道题主有没有解决，最近也碰到Negative Initial Temperature T0的问题

[李东岳]

@tens 这个是一种发散的表象。但是不是特别好处理，原因有很多，也不太好找是哪个。openfoam里面用了一种limit来强制处理，但有的时候也不能保证没有负温度出现。只能继续深入debug

[酱爆鸡丁]

lz解决了吗？我也遇到同样的问题了，高压气流冲击真空环境，使用rhoCentralFoam，算了一段时间之后就出现negative initial temperature t0，改了好久，可以推迟出现的时间，但总会出现，不知道怎么弄了