Problems about the normal shock capturing in PyFR

Dear all,
I recently initiated a normal shock in the flow field and found that it is unstable with the shock capturing methods including roem or roe. I think this phenomenon may be caused by the zero velocity in y direction. Any solutions to this problem?

Yours
Kiny

Hi Kiny,

For high-order methods it is necessary to also have a dedicated shock
capturing scheme inside of the volume of each element. In PyFR this is
configured through the shock-capturing parameter. We have some
documentation in the user guide on how to configure this and others on
the mailing list may be able to share some example cases with you.

Regards, Freddie.

Thank you. It is difficult to keep a normal shock wave stable with roem, and maybe other shock capturing schemes are reqiured.