Shock-density wave interaction

This test case is concerned with assessing the accuracy in predicting the interaction of a density wave with a shock front. For this, the Shu–Osher problem is extended by considering a binary mixture of Helium and Nitrogen which is initially separated by a shock. The shock Mach number is Ma = 3, and the computational domain is x ∈ [−5, 5]. Boundary conditions of supersonic outflow are imposed at both ends of the computational domain, and the initial conditions are given as:

The grid is uniformly discretised into 0.01 m elements size (mesh.geo file), time step is 1e-6 s.

The non-diagonal compressible tau of 2001 and the 2006 discontinuity capturing operator are adopted (setup.txt file).

Initial and boundary conditions are set in the fluid_ic_bc.hpp file, while the two material properties are specified in the props.txt file.

Results are compared at t = 1.8 s.

Gales computations well agree with the reference solution.

References

Yu Lv , Matthias Ihme “Discontinuous Galerkin method for multicomponent chemically reacting flows and combustion” (2014) J. Comp. Phy.

Amareshwara Sainadh Chamarthi “Consistent Interface Capturing Adaptive Reconstruction Approach for Viscous Compressible Multicomponent Flows”