Communications on Applied Mathematics and Computation ›› 2024, Vol. 6 ›› Issue (4): 2119-2154.doi: 10.1007/s42967-023-00313-6

• ORIGINAL PAPERS • Previous Articles     Next Articles

High-Order ADER Discontinuous Galerkin Schemes for a Symmetric Hyperbolic Model of Compressible Barotropic Two-Fluid Flows

Laura Río-Martín1,2, Michael Dumbser1   

  1. 1 Department of Civil, Environmental and Mechanical Engineering, University of Trento, Via Mesiano 77, Trento 38123, Italy;
    2 Department of Information Engineering and Computer Science, University of Trento, via Sommarive 9, Povo, Trento 38123, Italy
  • Received:2023-03-28 Revised:2023-07-31 Accepted:2023-08-24 Published:2024-12-20
  • Contact: Michael Dumbser,E-mail:michael.dumbser@unitn.it;Laura Río-Martín,E-mail:laura.delrio@unitn.it E-mail:michael.dumbser@unitn.it;laura.delrio@unitn.it
  • Supported by:
    Research was supported in part by the Office of Naval Research (ONR) N00014-13-1-0346, ONR N00014-17-1-2174, ARL AHPCRC W911NF-07-0027, and generous gifts from Amazon and Toyota.

Abstract: This paper presents a high-order discontinuous Galerkin (DG) finite-element method to solve the barotropic version of the conservative symmetric hyperbolic and thermodynamically compatible (SHTC) model of compressible two-phase flow, introduced by Romenski et al. in[59, 62], in multiple space dimensions. In the absence of algebraic source terms, the model is endowed with a curl constraint on the relative velocity field. In this paper, the hyperbolicity of the system is studied for the first time in the multidimensional case, showing that the original model is only weakly hyperbolic in multiple space dimensions. To restore the strong hyperbolicity, two different methodologies are used:(i) the explicit symmetrization of the system, which can be achieved by adding terms that contain linear combinations of the curl involution, similar to the Godunov-Powell terms in the MHD equations; (ii) the use of the hyperbolic generalized Lagrangian multiplier (GLM) curl-cleaning approach forwarded. The PDE system is solved using a high-order ADER-DG method with a posteriori subcell finite-volume limiter to deal with shock waves and the steep gradients in the volume fraction commonly appearing in the solutions of this type of model. To illustrate the performance of the method, several different test cases and benchmark problems have been run, showing the high order of the scheme and the good agreement when compared to reference solutions computed with other well-known methods.

Key words: Compressible two-fluid flows, Symmetric hyperbolic and thermodynamically compatible(SHTC)systems, Hyperbolic systems with curl involutions, High-order ADER discontinuous Galerkin(DG)schemes with subcell finite-volume limiter, Conservative form of hyperbolic models