超音速绕流发散

breathless

最近在做一个超音速绕流问题
大致形状是一个长圆柱（流速方向与母线相同），但是圆柱前段有5mm长陷进去的凹口。
边界距离圆柱大约10倍与圆柱长度，故选用远场压力条件，MA 1.2，根据公式算出湍流强度为3%，湍流尺寸为0.0004mm
控制方程为realized k-e
凹口处有patch为速度0m/s
密度基，steady.。
现在的问题是总是在凹口的zone有time step reduced due to excessive temperature change 然后再算几步就temperature limted to ........再算就发散了。
修改过库郎数至0.05，湍动能0.1，发散率0.1 但是只是治标不治本，还是发散，只不过迭代步数增加而已
网格质量check良好，网格数280w
怀疑是物理条件不合适，但无奈流体力学无能。
求高手答疑解惑。