A New Class of Simple, General and Efficient Finite Volume Schemes for Overdetermined Thermodynamically Compatible Hyperbolic Systems

Abstract

In this paper, a new efficient, and at the same time, very simple and general class of thermodynamically compatible finite volume schemes is introduced for the discretization of nonlinear, overdetermined, and thermodynamically compatible first-order hyperbolic systems. By construction, the proposed semi-discrete method satisfies an entropy inequality and is nonlinearly stable in the energy norm. A very peculiar feature of our approach is that entropy is discretized directly, while total energy conservation is achieved as a mere consequence of the thermodynamically compatible discretization. The new schemes can be applied to a very general class of nonlinear systems of hyperbolic PDEs, including both, conservative and non-conservative products, as well as potentially stiff algebraic relaxation source terms, provided that the underlying system is overdetermined and therefore satisfies an additional extra conservation law, such as the conservation of total energy density. The proposed family of finite volume schemes is based on the seminal work of Abgrall [1], where for the first time a completely general methodology for the design of thermodynamically compatible numerical methods for overdetermined hyperbolic PDE was presented. We apply our new approach to three particular thermodynamically compatible systems: the equations of ideal magnetohydrodynamics (MHD) with thermodynamically compatible generalized Lagrangian multiplier (GLM) divergence cleaning, the unified first-order hyperbolic model of continuum mechanics proposed by Godunov, Peshkov, and Romenski (GPR model) and the first-order hyperbolic model for turbulent shallow water flows of Gavrilyuk et al. In addition to formal mathematical proofs of the properties of our new finite volume schemes, we also present a large set of numerical results in order to show their potential, efficiency, and practical applicability.

Keywords

Overdetermined thermodynamically compatible hyperbolic systems / Hyperbolic and thermodynamically compatible (HTC) finite volume schemes / Abgrall framework / Discrete entropy inequality / Nonlinear stability in the energy norm / Applications to ideal magnetohydrodynamics (MHD), Godounov-Peshkov-Romenski (GPR) / Turbulent shallow water (TSW) flows

References

[1.] Abgrall R. A general framework to construct schemes satisfying additional conservation relations. Application to entropy conservative and entropy dissipative schemes. J. Comput. Phys., 2018, 372: 640-666

[2.] Abgrall R, Bacigaluppi P, Tokareva S. A high-order nonconservative approach for hyperbolic equations in fluid dynamics. Comput. Fluids, 2018, 169: 10-22

[3.] Abgrall R, Busto S, Dumbser M. A simple and general framework for the construction of thermodynamically compatible schemes for computational fluid and solid mechanics. Appl. Math. Comput., 2023, 440

[4.] Abgrall, R., Öffner, P., Ranocha, H.: Reinterpretation and extension of entropy correction terms for residual distribution and discontinuous Galerkin schemes: application to structure preserving discretization. J. Comput. Phys. 453, 110955 (2022)

[5.] Balsara D. Second-order accurate schemes for magnetohydrodynamics with divergence-free reconstruction. Astrophys. J. Suppl. Ser., 2004, 151: 149-184

[6.] Balsara D. Multidimensional HLLE Riemann solver: application to Euler and magnetohydrodynamic flows. J. Comput. Phys., 2010, 229: 1970-1993

[7.] Balsara D, Dumbser M. Divergence-free MHD on unstructured meshes using high order finite volume schemes based on multidimensional Riemann solvers. J. Comput. Phys., 2015, 299: 687-715

[8.] Balsara D, Spicer D. A staggered mesh algorithm using high order Godunov fluxes to ensure solenoidal magnetic fields in magnetohydrodynamic simulations. J. Comput. Phys., 1999, 149: 270-292

[9.] Becker R. . Stosswelle und Detonation. Physik, 1923, 8: 321

[10.] Bell JB, Coletta P, Glaz HM. A second-order projection method for the incompressible Navier-Stokes equations. J. Comput. Phys., 1989, 85: 257-283

[11.] Bhole A, Nkonga B, Gavrilyuk S, Ivanova K. Fluctuation splitting Riemann solver for a non-conservative modeling of shear shallow water flow. J. Comput. Phys., 2019, 392: 205-226

[12.] Bonnet A, Luneau J. . Aérodynamique, 1989 Toulouse Théories de la dynamique des fluides. Cepadues Editions

[13.] Brock R. Development of roll-wave trains in open channels. J. Hydraul. Div., 1969, 95: 1401-1428

[14.] Brock R. Periodic permanent roll waves. J. Hydraul. Div., 1970, 96: 2565-2580

[15.] Busto, S., Dumbser, M.: A new thermodynamically compatible finite volume scheme for magnetohydrodynamics. SIAM J. Numer. Anal. (2023). In press

[16.] Busto S, Dumbser M, Gavrilyuk S, Ivanova K. On thermodynamically compatible finite volume methods and path-conservative ADER discontinuous Galerkin schemes for turbulent shallow water flows. J. Sci. Comput., 2021, 88: 28

[17.] Busto S, Dumbser M, Peshkov I, Romenski E. On thermodynamically compatible finite volume schemes for continuum mechanics. SIAM J. Sci. Comput., 2022, 44: A1723-A1751

[18.] Busto S, Río-Martín L, Vázquez-Cendón ME, Dumbser M. A semi-implicit hybrid finite volume/finite element scheme for all Mach number flows on staggered unstructured meshes. Appl. Math. Comput., 2021, 402

[19.] Castro M, Gallardo J, Parés C. High-order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products. Applications to shallow-water systems. Math. Comput., 2006, 75: 1103-1134

[20.] Castro MJ, Fjordholm US, Mishra S, Parés C. Entropy conservative and entropy stable schemes for nonconservative hyperbolic systems. SIAM J. Numer. Anal., 2013, 51(3): 1371-1391

[21.] Chan J, Lin Y, Warburton T. Entropy stable modal discontinuous Galerkin schemes and wall boundary conditions for the compressible Navier-Stokes equations. J. Comput. Phys., 2022, 448

[22.] Chan J, Taylor CG. Efficient computation of Jacobian matrices for entropy stable summation-by-parts schemes. J. Comput. Phys., 2022, 448

[23.] Chandrashekar P, Klingenberg C. Entropy stable finite volume scheme for ideal compressible MHD on 2-D Cartesian meshes. SIAM J. Numer. Anal., 2016, 54(2): 1313-1340

[24.] Chandrashekar P, Nkonga B, Meena AM, Bhole A. A path conservative finite volume method for a shear shallow water model. J. Comput. Phys., 2020, 413

[25.] Cheng T, Shu C-W. Entropy stable high order discontinuous Galerkin methods with suitable quadrature rules for hyperbolic conservation laws. J. Comput. Phys., 2017, 345: 427-461

[26.] Dedner A, Kemm F, Kröner D, Munz CD, Schnitzer T, Wesenberg M. Hyperbolic divergence cleaning for the MHD equations. J. Comput. Phys., 2002, 175: 645-673

[27.] Derigs D, Winters AR, Gassner G, Walch S, Bohm M. Ideal GLM-MHD: about the entropy consistent nine-wave magnetic field divergence diminishing ideal magnetohydrodynamics equations. J. Comput. Phys., 2018, 364: 420-467

[28.] Dhaouadi F, Dumbser M. A first order hyperbolic reformulation of the Navier-Stokes-Korteweg system based on the GPR model and an augmented Lagrangian approach. J. Comput. Phys., 2022, 470

[29.] Dhaouadi F, Favrie N, Gavrilyuk S. Extended Lagrangian approach for the defocusing nonlinear Schrödinger equation. Stud. Appl. Math., 2019, 142: 336-358

[30.] Dumbser M, Peshkov I, Romenski E, Zanotti O. High order ADER schemes for a unified first order hyperbolic formulation of continuum mechanics: viscous heat-conducting fluids and elastic solids. J. Comput. Phys., 2016, 314: 824-862

[31.] Falle SAEG. Rarefaction shocks, shock errors and low order of accuracy in ZEUS. Astrophys. J., 2002, 577: L123-L126

[32.] Falle, S.A.E.G., Komissarov, S.: On the inadmissibility of non-evolutionary shocks. J. Plasma Phys. 65, 29–58 (2001)

[33.] Favrie N, Gavrilyuk S. A rapid numerical method for solving Serre-Green-Naghdi equations describing long free surface gravity waves. Nonlinearity, 2017, 30: 2718-2736

[34.] Fjordholm US, Mishra S. Accurate numerical discretizations of non-conservative hyperbolic systems. ESAIM Math. Model. Numer. Anal., 2012, 46(1): 187-206

[35.] Fjordholm US, Mishra S, Tadmor E. Arbitrarily high-order accurate entropy stable essentially nonoscillatory schemes for systems of conservation laws. SIAM J. Numer. Anal., 2012, 50(2): 544-573

[36.] Freistühler H. Relativistic barotropic fluids: a Godunov-Boillat formulation for their dynamics and a discussion of two special classes. Arch. Ration. Mech. Anal., 2019, 232: 473-488

[37.] Friedrichs K. Symmetric positive linear differential equations. Comm. Pure Appl. Math., 1958, 11: 333-418

[38.] Friedrichs K, Lax P. Systems of conservation equations with a convex extension. Proc. Nat. Acad. Sci. USA, 1971, 68: 1686-1688

[39.] Gaburro E, Öffner P, Ricchiuto M, Torlo D. High order entropy preserving ADER-DG schemes. Appl. Math. Comput., 2022, 440

[40.] Gassner G, Winters A, Kopriva D. A well balanced and entropy conservative discontinuous Galerkin spectral element method for the shallow water equations. Appl. Math. Comput., 2016, 272: 291-308

[41.] Gavrilyuk S, Ivanova K, Favrie N. Multi-dimensional shear shallow water flows: problems and solutions. J. Comput. Phys., 2018, 366: 252-280

[42.] Ghia U, Ghia KN, Shin CT. High-Re solutions for incompressible flow using Navier-Stokes equations and multigrid method. J. Comput. Phys., 1982, 48: 387-411

[43.] Godunov SK. An interesting class of quasilinear systems. Dokl. Akad. Nauk, 1961, 139(3): 521-523

[44.] Godunov SK. Symmetric form of the equations of magnetohydrodynamics. Numer. Methods Mech. Contin. Media., 1972, 3(1): 26-31

[45.] Godunov SK. Thermodynamic formalization of the fluid dynamics equations for a charged dielectric in an electromagnetic field. Comput. Math. Math. Phys., 2012, 52: 787-799

[46.] Godunov SK, Romenski E. Nonstationary equations of the nonlinear theory of elasticity in Euler coordinates. J. Appl. Mech. Tech. Phys., 1972, 13: 868-885

[47.] Godunov SK, Romenski EI. . Elements of Continuum Mechanics and Conservation Laws, 2003 New York Kluwer Academic/Plenum Publishers

[48.] Guermond J, Popov P. Viscous regularization of the Euler equations and entropy principles. SIAM J. Appl. Math., 2014, 74: 284-305

[49.] Hennemann S, Rueda-Ramírez AM, Hindenlang FJ, Gassner GJ. A provably entropy stable subcell shock capturing approach for high order split form DG for the compressible Euler equations. J. Comput. Phys., 2021, 426

[50.] Hiltebrand A, Mishra S. Entropy stable shock capturing space-time discontinuous Galerkin schemes for systems of conservation laws. Numer. Math., 2014, 126(1): 103-151

[51.] Hu C, Shu C-W. Weighted essentially non-oscillatory schemes on triangular meshes. J. Comput. Phys., 1999, 150: 97-127

[52.] Ivanova K, Gavrilyuk S. Structure of the hydraulic jump in convergent radial flows. J. Fluid Mech., 2019, 860: 441-464

[53.] Jiang G, Wu C. A high-order WENO finite difference scheme for the equations of ideal magnetohydrodynamics. J. Comput. Phys., 1999, 150: 561-594

[54.] Liu Y, Shu C-W, Zhang M. Entropy stable high order discontinuous Galerkin methods for ideal compressible MHD on structured meshes. J. Comput. Phys., 2018, 354: 163-178

[55.] Lukácová-Medvidóvá, M., Puppo, G., Thomann, A.: An all Mach number finite volume method for isentropic two-phase flow. J. Numer. Math. (2022)

[56.] Munz C, Omnes P, Schneider R, Sonnendrücker E, Voss U. Divergence correction techniques for maxwell solvers based on a hyperbolic model. J. Comput. Phys., 2000, 161: 484-511

[57.] Nkonga, B., Chandrashekar, P.: Exact solution for Riemann problems of the shear shallow water model. ESAIM Math. Modell. Numer. Anal. 56(4), 1115–1150 (2022)

[58.] Parés C. Numerical methods for nonconservative hyperbolic systems: a theoretical framework. SIAM J. Numer. Anal., 2006, 44: 300-321

[59.] Peshkov I, Pavelka M, Romenski E, Grmela M. Continuum mechanics and thermodynamics in the Hamilton and the Godunov-type formulations. Contin. Mech. Thermodyn., 2018, 30(6): 1343-1378

[60.] Peshkov I, Romenski E. A hyperbolic model for viscous Newtonian flows. Contin. Mech. Thermodyn., 2016, 28: 85-104

[61.] Peshkov I, Romenski E, Dumbser M. Continuum mechanics with torsion. Contin. Mech. Thermodyn., 2019, 31: 1517-1541

[62.] Ray D, Chandrashekar P. An entropy stable finite volume scheme for the two dimensional Navier-Stokes equations on triangular grids. Appl. Math. Comput., 2017, 314: 257-286

[63.] Ray D, Chandrashekar P, Fjordholm US, Mishra S. Entropy stable scheme on two-dimensional unstructured grids for Euler equations. Commun. Comput. Phys., 2016, 19(5): 1111-1140

[64.] Romenski E. Hyperbolic systems of thermodynamically compatible conservation laws in continuum mechanics. Math. Comput. Modell., 1998, 28(10): 115-130

[65.] Romenski E, Belozerov A, Peshkov I. Conservative formulation for compressible multiphase flows. Q. Appl. Math., 2016, 74: 113-136

[66.] Romenski E, Drikakis D, Toro E. Conservative models and numerical methods for compressible two-phase flow. J. Sci. Comput., 2010, 42: 68-95

[67.] Romenski E, Peshkov I, Dumbser M, Fambri F. A new continuum model for general relativistic viscous heat-conducting media. Philos. Trans. R. Soc. A, 2020, 378: 20190175

[68.] Romenski E, Resnyansky A, Toro E. Conservative hyperbolic formulation for compressible two-phase flow with different phase pressures and temperatures. Q. Appl. Math., 2007, 65: 259-279

[69.] Rueda-Ramírez, A.M., Hennemann, S., Hindenlang, F.J., Winters, A.R., Gassner, G.J.: An entropy stable nodal discontinuous Galerkin method for the resistive MHD equations. Part II: subcell finite volume shock capturing. J. Comput. Phys. 444 (2021)

[70.] Ruggeri, T., Strumia, A.: Main field and convex covariant density for quasilinear hyperbolic systems Relativistic fluid dynamics. Ann. Inst. Henri Poincaré 34 65–84 (1981)

[71.] Schnücke, G., Krais, N., Bolemann, T., Gassner, G.J.: Entropy stable discontinuous Galerkin schemes on moving meshes for hyperbolic conservation laws. J. Sci. Comput. 82, 69 (2020)

[72.] Tadmor E. The numerical viscosity of entropy stable schemes for systems of conservation laws I. Math. Comput., 1987, 49: 91-103

[73.] Tavelli M, Dumbser M. A pressure-based semi-implicit space-time discontinuous Galerkin method on staggered unstructured meshes for the solution of the compressible Navier-Stokes equations at all Mach numbers. J. Comput. Phys., 2017, 341: 341-376

[74.] Thein F, Romenski E, Dumbser M. Exact and numerical solutions of the Riemann problem for a conservative model of compressible two-phase flows. J. Sci. Comput., 2022, 93: 83

[75.] Toro E. . Riemann Solvers and Numerical Methods for Fluid Dynamics, 2009 Cham Springer

Funding

Ministerio de Ciencia, Innovación y Universidades(PID2021-122625OB-I00); Ministero dell’Istruzione, dell’Università e della Ricerca(PRIN 2017); European Commission(NextGenerationEU (PNRR, Spoke 7 CN HPC))