WENO 5th in Generalized Coordinate (Curvilinear Coordinate)

 Register Blogs Members List Search Today's Posts Mark Forums Read

 January 1, 2021, 06:19 WENO 5th in Generalized Coordinate (Curvilinear Coordinate) #1 New Member   Join Date: Jan 2021 Posts: 3 Rep Power: 3 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

 January 1, 2021, 12:10 #2 Senior Member     - Join Date: Jul 2012 Location: Germany Posts: 158 Rep Power: 12 First step would be to simplify your test problem/mesh. __________________ Check out my side project: A multiphysics discontinuous Galerkin framework: Youtube, Gitlab.

 January 1, 2021, 19:54 #3 New Member   Join Date: Jan 2021 Posts: 3 Rep Power: 3 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?

 January 1, 2021, 20:41 #4 Senior Member     - Join Date: Jul 2012 Location: Germany Posts: 158 Rep Power: 12 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 __________________ Check out my side project: A multiphysics discontinuous Galerkin framework: Youtube, Gitlab.

 January 1, 2021, 20:58 #5 New Member   Join Date: Jan 2021 Posts: 3 Rep Power: 3 Thank you so much! I will try some simulations along with your suggestions.