Richtmyer Meshkov - Split form discontinuous Galerkin - FEM - Euler Equations - Python/Matlab/C++
Richtmyer Meshkov Instability with Split-DG, supersonic-outflow BC and improved Shockcapturing
Density, Entropy, Velocity field
Python/Matlab/C++ Code:
Euler Equations
Cartesian Mesh with Nx x Ny = 128 x 384
Split-form discontinuous Galerkin Method (SDGSEM)
Rusanov flux on the faces and Chandrashekar flux in the volume
Lagrange Polynomials
Gauss-Lobatto Distribution
Polynomial Degree=2
Order of Convergence=3
Shockindicator: Persson Indicator
Shockcapturing: Finite Volume Subcell
Simulation area:
dx=13.333...
dy=40.0
State Top:
r=0.25
p=1
vx=0
vy=0
State Middle: (y≥18+2*cos(2*pi*3/40*x)
r=1
p=1
vx=0
vy=0
State Bottom: (2≤y≤6)
r=4.22
p=4.9
vx=0
vy=0