rhoCentralFoam发散，出现Negative Initial Temperature
-
算例文件在文中已经贴出
-
尝试画粗网格计算正方体外流，结果确实可以计算，deltat=10^-7算到0.05s，未出现以上问题。
对应精细网格，计算到300个时间步，又出现了温度陡增的发散问题。
细网格内能残差
粗网格内能残差两个算例的迭代中能量残差都是initial residual< final residual，不过粗网格没有发散。
在这些rhoC算例中计算发现，温度压强都有“龟派气功”的特效。。。前缘驻点附近形成的高压会脱体。。。后缘驻点类似。。。
到最后消失了。。。 -
@东岳 是的哇
-
@东岳 我其他的边界条件pT都是一样的，控制雷诺数，计算了不同的速度。因为是层流，所以上传的是最小速度1m/s的。
-
@东岳 那么高速导致雷诺数量级太大，层流模型还能算么？
-
https://www.cfd-china.com/topic/1737/可压缩求解器为什么密度残差一直为0-迭代步也一直是0
这里面回答了rho，rhoE，rhoU残差为0的原因。
-
@东岳 重新试了一下250m/s，相当于0.83Ma，算了20步就出现负初始温度的问题
-
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}$下降一个量级。
-
沉了么。。？有大佬试试我的算例了么？
-
，