高超声速平板或尖锥怎么避免等温壁面下经常发散的情况

[Briareos]

本人目前在算10马赫左右的高超声速平板和尖锥的问题，雷诺数为10e4量级，流动为层流，壁面网格雷诺数为0.1的量级，壁面条件是等温边界条件，为了提高稳定性采用的是马赫数逐步提高的方式，但还是经常在头部区域会出现发散或者物理量不正常的情况，虚拟网格温度采用线性外插，压力梯度为零，密度根据理想气体方程算出来。不知道有没有大神做过或了解过什么方法能算出来正确的物理量并且不发散。图附的是在平板头部出现的奇怪的东西。。。