# How PyFR handle shock capture problem?

**URL:** <https://pyfr.discourse.group/t/how-pyfr-handle-shock-capture-problem/664>\
**Category:** General\
**Created:** [4 August 2022 07:06 UTC](https://pyfr.discourse.group/t/how-pyfr-handle-shock-capture-problem/664 "2022-08-04T07:06:48Z")\
**Posts on this page:** 3\
**Page:** 1

<div class="post-metadata">

**Author:** ![luli](https://avatars.discourse-cdn.com/v4/letter/l/5f9b8f/32.png) [@luli](https://pyfr.discourse.group/u/luli)\
**Post date:** [4 August 2022 07:06 UTC](https://pyfr.discourse.group/t/how-pyfr-handle-shock-capture-problem/664/1 "2022-08-04T07:06:48Z")

</div>

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 ?

---

<div class="post-metadata">

**Author:** ![WillT](https://avatars.discourse-cdn.com/v4/letter/w/bc79bd/32.png) [@WillT](https://pyfr.discourse.group/u/WillT)\
**Post date:** [4 August 2022 10:07 UTC](https://pyfr.discourse.group/t/how-pyfr-handle-shock-capture-problem/664/2 "2022-08-04T10:07:44Z")

</div>

At the moment the only option available is the artifical viscousity method similar to that of [Persson and Peraire](https://doi.org/10.2514/6.2006-112).

In the future I think theb plan is to implement the method @tdzanic developed which is shown [here](https://doi.org/10.1016/j.jcp.2022.111501). 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.

---

<div class="post-metadata">

**Author:** ![nnunn](https://avatars.discourse-cdn.com/v4/letter/n/c89c15/32.png) [@nnunn](https://pyfr.discourse.group/u/nnunn)\
**Post date:** [5 August 2022 10:28 UTC](https://pyfr.discourse.group/t/how-pyfr-handle-shock-capture-problem/664/3 "2022-08-05T10:28:33Z")

</div>

Just wanted to add that Tarik’s work on this is VERY exciting!