Triple-Point Shock Interaction - Discontinuous Galerkin - FEM - Euler-Equations - Matlab/Python/C++
Thanks to Tarik Dzanic for the nice example!
Logarithm of Density
Code Runtime: 1 day on only 4 cores
Python/C++ Code:
Positivity Preserving (Some random own stuff)
Mesh Nx x Ny: 862x370
Nelems=317709, Ndofs = 5083344
Compressible Euler Equations
Polynomial Degree=3
Order of Convergence=4
Standard discontinuous Galerkin Method (DGSEM)
Modified HLLC Riemann solver
Lagrange Polynomials
Gauss-Legrende Distribution
Shockindicator: Persson Indicator
Shockcapturing: Finite Volume Subcell