When applying the higher-order scheme, the capture of shock is a difficult problem. How does pyfr deal with it ? I saw “artificial viscosity” option in user guide , is there any limiter can I choose ?
At the moment the only option available is the artifical viscousity method similar to that of Persson and Peraire.
In the future I think theb plan is to implement the method @tdzanic developed which is shown here. This is currently work in progress as the method works well but some modifications are needed to get the code working at the computational efficiency of main line PyFR.
Just wanted to add that Tarik’s work on this is VERY exciting!