Hey All, I am running a case in hy2Foam I did implement a Riemann solver, AUSMPlus. I was trying a high speed scenario over a axisymmetric double cone. It seems that the temperature diverges early on and an aberration forms on the area with some cell skewness (the mean Co is small while the MaxCo is over 1). I changed the Co calculation from cell centers to face values which helps a bit. It seems to work with the Kurganov method and was wondering if (a) just improve the mesh (b) something I was missing in how I handle flux calculations. Thanks in advance. I can post the relevant code(s) as well.