[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

[Eifoehn4]

First step would be to simplify your test problem/mesh.

[MFurufuru]

Thank you so much for your reply.
I thought this problem setting was the most simple case in 2D.

Do you have some suggestions for the problem setting?
What problem and mesh should I use?

[Eifoehn4]

I suppose you have implemented the Euler or Navier-Stokes equations:

• Check if your code fulfills the free stream preservation property.
• Check for conservativity e.g. mass/energy/momentum.
• Check for symmetries, e.g. simulate a moving contact discontinuity from left to right or from down to top or vice versa.
• Implement the manufacture solution procedure and check your convergence order.
• Modify your mesh only in one direction.
• Start with a more simple equation system e.g. advection/diffusion or burgers equations
• ...

There are so many things you can try just proceed according to the exclusion procedure.

Regards

[MFurufuru]

Thank you so much!
I will try some simulations along with your suggestions.