[MFurufuru]

Hi.

I am developing WENO 5th finite difference method Fortran code because I want to simulate hypersonic flow with bodies.

I developed 2D-rectanglar uniform grid solver as shown in the figure.(Oblique shock wave, AOA=10 deg, Mach number is 2.)
I think this result is reasonable.

So, I decided to add metrics and Jacobian and simulate same problem in Curvilinear Grid.
However, the result is weird.

Did you have same problem like these figure?
It seems my code is correct. (metrics and jacobian are reasonably added, I think...)

I'll appreciate your comments!
Thank you.

WENO 5th:
Flux : Lax-Friedrichs flux splitting (Global)
parameter epsilon=1e-6

Time Discretization:
3rd order TVD Runge-Kutta