High-order enforcement of jumps conditions between compressible viscous phases: An extended interior penalty discontinuous Galerkin method for sharp interface simulation

• Published in: Computer methods in applied mechanics and engineering 2023-10, Vol.415, p.116215, Article 116215
• Main Authors: Henneaux, David, Schrooyen, Pierre, Chatelain, Philippe, Magin, Thierry
• Format: Article
• Language: English
• Subjects: Extended; Interior penalty discontinuous Galerkin; Sharp interface; Two-phase flow; Unfitted; Viscous compressible
• ISSN: 0045-7825; 1879-2138
• DOI: 10.1016/j.cma.2023.116215

In this paper, we develop a high-order method for the steady two-phase compressible Navier–Stokes equations closed by generic equations of state, adapted to gas–liquid flows. Our framework relies on the extended discontinuous Galerkin method that is tailored here to a general viscous compressible two-phase configuration. While the imposition of interface jumps conditions is commonly recognized as one of the main difficulties to overcome in the context of geometrically unfitted methods, the setting considered in this work makes this aspect even more challenging. To enforce the non-linear coupling conditions at the interface in a sharp manner, we devise a novel strategy combining multiphase Riemann solvers to handle the convective fluxes across the interface and a weighted stabilized Nitsche’s method for the viscous jumps. The weak enforcement of the jumps conditions along with implicit steady-state iterations and cell-agglomeration procedure make the approach robust to treat arbitrary large contrast phase interface problems governed by stiff models, irrespective of the interface location on the mesh. The study is here restricted to static interfaces. Based on the method of manufactured solution, reference solutions are designed. This allows to conduct a detailed numerical study investigating the influence of several numerical parameters. It is shown that to reach optimal error convergence, the use of pressure-velocity-temperature variables set, appropriate all-speed bulk Riemann solver and Symmetric Interior Penalty method combined with tensorial penalty is necessary. [Display omitted] •Derivation of a geometrically unfitted extended discontinuous Galerkin discretization.•Development of a novel strategy to weakly enforce the non-linear interface jumps.•Application to a steady two-phase compressible Navier–Stokes system.•Use of multiphase Riemann solvers and a revisited interior penalty method.•Study about error convergence, parameters sensitivity and geometrical robustness.•Design of complex manufactured reference solutions.