loganoz / horses3d

HORSES3D: A high-order discontinuous Galerkin solver for flow simulations and multi-physics applications
https://loganoz.github.io/horses3d/
MIT License
112 stars 24 forks source link

Add more explicit schemes and stage limiters #116

Closed Andres-MG closed 1 year ago

Andres-MG commented 1 year ago

I will add more explicit schemes to the code since the ones we have are pretty basic and not specific to hyperbolic equations or the NS equations. So far I have implemented two Strong Stability Preserving RK schemes. Both are 3rd order; however, the SSPKR43 has 4 stages and its SSP constant is 2.

There are more ideas here and here. I would also like to add some stage limiters to avoid crashes when the density and/or the pressure become negative.

loganoz commented 1 year ago

Great addition @Andres-MG. Don't forget to add the keywords to the documentation 👍

loganoz commented 1 year ago

@Andres-MG , go ahead and merge in master. Don't forget to add a test case and include in the actions.

Andres-MG commented 1 year ago

After having been using it for some time now, I would like to modify the limiter a bit. Anyway, it seems to help a lot when small shocks and discontinuities appear (M~1 and initialization).