CFD中文网

@shifc 不知道你现在解决了没有，我最近做格栅湍流也碰到了同样的问题，我的解决方法是谱空间的复速度对三个坐标轴必须都是共轭对称的。

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