Well-balanced Arbitrary-Lagrangian-Eulerian finite volume schemes on moving nonconforming meshes for the Euler equations of gas dynamics with gravity

In this work we present a novel second order accurate well balanced Arbitrary-Lagrangian-Eulerian (ALE) finite volume scheme on moving nonconforming meshes for the Euler equations of compressible gasdynamics with gravity in cylindrical coordinates. The main feature of the proposed algorithm is the capability of preserving many of the physical properties of the system exactly also on the discrete level: besides being conservative for mass, momentum and total energy, also any known steady equilibrium between pressure gradient, centrifugal force and gravity force can be exactly maintained up to machine precision. Perturbations around such equilibrium solutions are resolved with high accuracy and with minimal dissipation on moving contact discontinuities even for very long computational times. This is achieved by the novel combination of well balanced path-conservative finite volume schemes, that are expressly designed to deal with source terms written via nonconservative products, with ALE schemes on moving grids, which exhibit only very little numerical dissipation on moving contact waves. In particular, we have formulated a new HLL-type and a novel Osher-type flux that are both able to guarantee the well balancing in a gas cloud rotating around a central object. Moreover, to maintain a high level of quality of the moving mesh, we have adopted a nonconforming treatment of the sliding interfaces that appear due to the differential rotation. A large set of numerical tests has been carried out in order to check the accuracy of the method close and far away from the equilibrium, both, in one and two space dimensions.

[1]  Mikhail J. Shashkov,et al.  One-step hybrid remapping algorithm for multi-material arbitrary Lagrangian-Eulerian methods , 2012, J. Comput. Phys..

[2]  Bruno Després,et al.  A one-mesh method for the cell-centered discretization of sliding , 2014 .

[3]  Jérôme Breil,et al.  Two-step hybrid conservative remapping for multimaterial arbitrary Lagrangian-Eulerian methods , 2011, J. Comput. Phys..

[4]  Bruno Després,et al.  Lagrangian Gas Dynamics in Two Dimensions and Lagrangian systems , 2005 .

[5]  Michael Dumbser,et al.  High order ADER schemes for a unified first order hyperbolic formulation of Newtonian continuum mechanics coupled with electro-dynamics , 2016, J. Comput. Phys..

[6]  Manuel J. Castro,et al.  WELL-BALANCED NUMERICAL SCHEMES BASED ON A GENERALIZED HYDROSTATIC RECONSTRUCTION TECHNIQUE , 2007 .

[7]  Pierre-Henri Maire,et al.  Multi-scale Godunov-type method for cell-centered discrete Lagrangian hydrodynamics , 2009, J. Comput. Phys..

[8]  P. Murdin.  MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY , 2005 .

[9]  Christian Klingenberg,et al.  A Well-Balanced Scheme for the Euler Equation with a Gravitational Potential , 2014 .

[10]  Manuel Jesús Castro Díaz,et al.  High order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products. Applications to shallow-water systems , 2006, Math. Comput..

[11]  A. Ferrari,et al.  PLUTO: A Numerical Code for Computational Astrophysics , 2007, astro-ph/0701854.

[12]  F. Bouchut.  Nonlinear Stability of Finite Volume Methods for Hyperbolic Conservation Laws: and Well-Balanced Schemes for Sources , 2005 .

[13]  Carlos Parés,et al.  On the well-balance property of Roe?s method for nonconservative hyperbolic systems , 2004 .

[14]  J. Greenberg,et al.  A well-balanced scheme for the numerical processing of source terms in hyperbolic equations , 1996 .

[15]  C. Parés.  Numerical methods for nonconservative hyperbolic systems: a theoretical framework. , 2006 .

[16]  Alfredo Bermúdez,et al.  Upwind methods for hyperbolic conservation laws with source terms , 1994 .

[17]  Carlos Parés,et al.  A Q-SCHEME FOR A CLASS OF SYSTEMS OF COUPLED CONSERVATION LAWS WITH SOURCE TERM. APPLICATION TO A TWO-LAYER 1-D SHALLOW WATER SYSTEM , 2001 .

[18]  Chiara Simeoni,et al.  Convergence of the Upwind Interface Source Method for Hyperbolic Conservation Laws , 2003 .

[19]  Michael Dumbser,et al.  A direct Arbitrary-Lagrangian-Eulerian ADER-WENO finite volume scheme on unstructured tetrahedral meshes for conservative and non-conservative hyperbolic systems in 3D , 2014, J. Comput. Phys..

[20]  J. Greenberg,et al.  Analysis and Approximation of Conservation Laws with Source Terms , 1997 .

[21]  Pavel B. Bochev,et al.  Fast optimization-based conservative remap of scalar fields through aggregate mass transfer , 2013, J. Comput. Phys..

[22]  Enrique D. Fernández Nieto,et al.  A family of stable numerical solvers for the shallow water equations with source terms , 2003 .

[23]  Bruno Després,et al.  A cell-centered Lagrangian hydrodynamics scheme on general unstructured meshes in arbitrary dimension , 2009, J. Comput. Phys..

[24]  Eleuterio F. Toro,et al.  On some fast well-balanced first order solvers for nonconservative systems , 2009, Math. Comput..

[25]  Pavel Váchal,et al.  Synchronized flux corrected remapping for ALE methods , 2011 .

[26]  Alfredo Bermúdez,et al.  Numerical solution of non-isothermal non-adiabatic flow of real gases in pipelines , 2016, J. Comput. Phys..

[27]  Emmanuel Audusse,et al.  A Fast and Stable Well-Balanced Scheme with Hydrostatic Reconstruction for Shallow Water Flows , 2004, SIAM J. Sci. Comput..

[28]  Michael Dumbser,et al.  Reprint of: Direct Arbitrary-Lagrangian-Eulerian finite volume schemes on moving nonconforming unstructured meshes , 2016, Computers & Fluids.

[29]  Raphaël Loubère,et al.  Enhancement of Lagrangian slide lines as a combined force and velocity boundary condition , 2013 .

[30]  Michael Dumbser,et al.  Arbitrary-Lagrangian-Eulerian One-Step WENO Finite Volume Schemes on Unstructured Triangular Meshes , 2013, 1302.3076.

[31]  S. Mishra,et al.  Well-balanced schemes for the Euler equations with gravitation , 2014, J. Comput. Phys..

[32]  Siddhartha Mishra,et al.  A well-balanced finite volume scheme for the Euler equations with gravitation - The exact preservation of hydrostatic equilibrium with arbitrary entropy stratification , 2016 .

[33]  V. Springel.  E pur si muove: Galilean-invariant cosmological hydrodynamical simulations on a moving mesh , 2009, 0901.4107.

[34]  Chi-Wang Shu,et al.  High order conservative Lagrangian schemes with Lax-Wendroff type time discretization for the compressible Euler equations , 2009, J. Comput. Phys..

[35]  Chi-Wang Shu,et al.  A high order ENO conservative Lagrangian type scheme for the compressible Euler equations , 2007, J. Comput. Phys..

[36]  L. Gosse.  A well-balanced flux-vector splitting scheme designed for hyperbolic systems of conservation laws with source terms☆ , 2000 .

[37]  Michael Dumbser,et al.  Arbitrary-Lagrangian-Eulerian Discontinuous Galerkin schemes with a posteriori subcell finite volume limiting on moving unstructured meshes , 2016, J. Comput. Phys..

[38]  M. Dumbser,et al.  A well-balanced path conservative SPH scheme for nonconservative hyperbolic systems with applications to shallow water and multi-phase flows , 2017 .

[39]  Chaowei Hu,et al.  No . 98-32 Weighted Essentially Non-Oscillatory Schemes on Triangular Meshes , 1998 .

[40]  Kun Xu,et al.  A gas-kinetic scheme for shallow-water equations with source terms , 2004 .

[41]  E. J. Caramana,et al.  The implementation of slide lines as a combined force and velocity boundary condition , 2009, J. Comput. Phys..

[42]  Michael Dumbser,et al.  A new efficient formulation of the HLLEM Riemann solver for general conservative and non-conservative hyperbolic systems , 2016, J. Comput. Phys..

[43]  P. Roe,et al.  On Godunov-type methods near low densities , 1991 .

[44]  P. Colella,et al.  THE PLUTO CODE FOR ADAPTIVE MESH COMPUTATIONS IN ASTROPHYSICAL FLUID DYNAMICS , 2011, 1110.0740.

[45]  Randall J. LeVeque,et al.  Balancing Source Terms and Flux Gradients in High-Resolution Godunov Methods , 1998 .

[46]  Lucas O. Müller,et al.  Well-balanced high-order numerical schemes for one-dimensional blood flow in vessels with varying mechanical properties , 2013, J. Comput. Phys..

[47]  Christian Klingenberg,et al.  A well‐balanced scheme to capture non‐explicit steady states in the Euler equations with gravity , 2016 .

[48]  Manuel Jesús Castro Díaz,et al.  Approximate Osher-Solomon schemes for hyperbolic systems , 2016, Appl. Math. Comput..

[49]  Michael Dumbser,et al.  High order ADER schemes for a unified first order hyperbolic formulation of continuum mechanics: Viscous heat-conducting fluids and elastic solids , 2015, J. Comput. Phys..

[50]  Laurent Gosse,et al.  A Well-Balanced Scheme Using Non-Conservative Products Designed for Hyperbolic Systems of Conservati , 2001 .

[51]  Silvia Bertoluzza,et al.  A conservative slide line method for cell-centered semi-Lagrangian and ALE schemes in 2D , 2016 .

[52]  Michael Dumbser,et al.  A Simple Extension of the Osher Riemann Solver to Non-conservative Hyperbolic Systems , 2011, J. Sci. Comput..

[53]  Manuel Jesús Castro Díaz,et al.  A Class of Computationally Fast First Order Finite Volume Solvers: PVM Methods , 2012, SIAM J. Sci. Comput..

[54]  Manuel Jesús Castro Díaz,et al.  On Well-Balanced Finite Volume Methods for Nonconservative Nonhomogeneous Hyperbolic Systems , 2007, SIAM J. Sci. Comput..

[55]  Ilya Peshkov,et al.  On a pure hyperbolic alternative to the Navier-Stokes equations , 2014 .

[56]  B. Perthame,et al.  A kinetic scheme for the Saint-Venant system¶with a source term , 2001 .

[57]  Michael Dumbser,et al.  High order cell-centered Lagrangian-type finite volume schemes with time-accurate local time stepping on unstructured triangular meshes , 2014, J. Comput. Phys..

[58]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[59]  Enrique D. Fernández Nieto,et al.  Asymptotically balanced schemes for non-homogeneous hyperbolic systems – application to the Shallow Water equations , 2004 .

[60]  Christian Klingenberg,et al.  Astrophysical hydrodynamics with a high-order discontinuous Galerkin scheme and adaptive mesh refinement , 2015, 1506.06140.

[61]  E. Toro.  Shock-Capturing Methods for Free-Surface Shallow Flows , 2001 .

[62]  Rupert Klein,et al.  Well balanced finite volume methods for nearly hydrostatic flows , 2004 .

[63]  Bruno Després,et al.  A new method to introduce constraints in cell-centered Lagrangian schemes , 2013 .

[64]  Guglielmo Scovazzi,et al.  A geometrically-conservative, synchronized, flux-corrected remap for arbitrary Lagrangian-Eulerian computations with nodal finite elements , 2011, J. Comput. Phys..

[65]  Michael Dumbser,et al.  On Universal Osher-Type Schemes for General Nonlinear Hyperbolic Conservation Laws , 2011 .

[66]  Rémi Abgrall,et al.  A Cell-Centered Lagrangian Scheme for Two-Dimensional Compressible Flow Problems , 2007, SIAM J. Sci. Comput..

[67]  M. Dumbser,et al.  High-Order Unstructured Lagrangian One-Step WENO Finite Volume Schemes for Non-Conservative Hyperbolic Systems: Applications to Compressible Multi-Phase Flows , 2013, 1304.4816.