Implicit shock fitting for multimaterial shock dynamics using a high-order space-time discontinuous finite-element method

A new numerical method for high-order discretization of fluid flows with discontinuities is introduced. This kind of fluid flow is likely to happen in reactor cores in case of accident. The method combines a space-time finite-element discretization with shock-fitting methodology, implicitly solving for both the finite-element degrees of freedom of conservation laws and the mesh geometry, incorporating interfacial jump conditions into the Newton-based non-linear solver, forcing mesh alignment with all solution and solution derivative discontinuities. The method feature unique capabilities to truly converge with high-order, even in the presence of shocks