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?
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.