[Victor Mourão]

Hi! I am new in the forum, I work with cosmology (a RHD approach of QCD phase transition with a realist EoS) and am developing a numerical code using WENO and/or Monotonicity Preserving schemes (Suresh 1997).

To start, I am implementing a HD using a algorithm described in Balsara et al. 2000, Mignone et al. 2010. and Radice et al. 2012

Basically, I use a Lax-Friedrichs flux split and in a characteristic space apply a reconstruction method (WENO3, WENO5, WENO7 and MP5).

But my validation tests don't give a satisfactory results in a strong shock test and in a Lax problem. At the top of the shocks are creating oscillations which should not exist. Someone already had these problems?

I am adding some pictures of the tests that I did.

Best Reagards

http://imageshack.us/a/img201/2378/strong.png

http://imageshack.us/a/img211/6888/sodx.png

http://imageshack.us/a/img845/6573/lax.png

http://imageshack.us/a/img685/9017/blastwave.png

http://imageshack.us/a/img163/994/123ckk.png

[reza_65]

Hi,

Did you implement this high order method in OpenFOAM? And can I ask you whether there are available somewhere or not?

Regards,
Reza