i give the value of k and e like this :114.84375,1619396.7 ,which is caculated in line with the formulation of the NUMECA menu.i am not sure where the problem lies in ,mesh or the boundary conditions.when it runs to the fine mesh ,it just appears the information that reads pressure negative and density negative .

I think it is the problem of mesh,maybe you should check the quality of the mash such as orthogonality,Aspect ratio, expansion ratio and so on,I have met such problem,good lucky.

thank you a lot !i then try it .i appreciate your kindness.