Problems about the normal shock capturing in PyFR

[Kiny Wan]

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

[Freddie Witherden]

Hi Kiny,

On 23/11/2019 20:25, Kiny Wan wrote:
>         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.

Regards, Freddie.

[Kiny Wan]

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