Shock/Mach diamonds - Discontinuous Galerkin - FEM - Euler Equations - Python/Matlab/C++
Shock/Mach diamonds
Density (Schlieren like)
Python/Matlab/C++ Code:
Euler Equations
Stretched mesh with Nx x Ny = 350 x 182
Standard discontinuous Galerkin Method (DGSEM)
Rusanov flux
Lagrange Polynomials
Gauss-Legrende Distribution
Polynomial Degree=2
Order of Convergence=3
Shockindicator: Persson Indicator
Shockcapturing: Finite Volume Subcell
Simulation area:
D=0.001
dx=0.015
dy=0.03
Chamber State:
r=1.16102755
p=100000
vx=0
vy=0
Inlet State:
r=0.49062949
p=64300
vx=1045.15934123
vy=0