Premium
A reconstructed discontinuous Galerkin method for multi‐material hydrodynamics with sharp interfaces
Author(s) -
Pandare Aditya K.,
Waltz Jacob,
Bakosi Jozsef
Publication year - 2020
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.4810
Subject(s) - discretization , discontinuous galerkin method , total variation diminishing , mathematics , work (physics) , flow (mathematics) , finite element method , order of accuracy , finite volume method , galerkin method , mathematical optimization , numerical analysis , computer science , mathematical analysis , mechanics , physics , geometry , numerical stability , thermodynamics
Summary Discontinuous Galerkin (DG) methods have been well established for single‐material hydrodynamics. However, consistent DG discretizations for non‐equilibrium multi‐material (more than two materials) hydrodynamics have not been extensively studied. In this work, a novel reconstructed DG (rDG) method for the single‐velocity multi‐material system is presented. The multi‐material system being considered assumes stiff velocity relaxation, but does not assume pressure and temperature equilibrium between the multiple materials. A second‐order DG(P 1 ) method and a third‐order least‐squares based rDG(P 1 P 2 ) are used to discretize this system in space, and a third‐order total variation diminishing (TVD) Runge‐Kutta method is used to integrate in time. A well‐balanced DG discretization of the non‐conservative system is presented and is verified by numerical test problems. Furthermore, a consistent interface treatment is implemented, which ensures strict conservation of material masses and total energy. Numerical tests indicate that the DG and rDG methods are, indeed, the second‐ and third‐order accurate. Comparisons with the second‐order finite volume method show that the DG and rDG methods are able to capture the interfaces more sharply. The DG and rDG methods are also more accurate in the single‐material regions of the flow. This work focuses on the general multidimensional rDG formulation of the non‐equilibrium multi‐material system and a study of properties of the method via one‐dimensional numerical experiments. The results from this research will be the foundation for a multidimensional high‐order rDG method for multi‐material hydrodynamics.