Well-balanced finite volume schemes for nearly steady adiabatic flows
arXiv:2006.05960 · doi:10.1016/j.jcp.2020.109805
Abstract
We present well-balanced finite volume schemes designed to approximate the Euler equations with gravitation. They are based on a novel local steady state reconstruction. The schemes preserve a discrete equivalent of steady adiabatic flow, which includes non-hydrostatic equilibria. The proposed method works in Cartesian, cylindrical and spherical coordinates. The scheme is not tied to any specific numerical flux and can be combined with any consistent numerical flux for the Euler equations, which provides great flexibility and simplifies the integration into any standard finite volume algorithm. Furthermore, the schemes can cope with general convex equations of state, which is particularly important in astrophysical applications. Both first- and second-order accurate versions of the schemes and their extension to several space dimensions are presented. The superior performance of the well-balanced schemes compared to standard schemes is demonstrated in a variety of numerical experiments. The chosen numerical experiments include simple one-dimensional problems in both Cartesian and spherical geometry, as well as two-dimensional simulations of stellar accretion in cylindrical geometry with a complex multi-physics equation of state.
References in corpus (13)
- Stability of Standing Accretion Shocks, With an Eye Toward Core Collapse Supernovae
- The explosion mechanism of core-collapse supernovae: progress in supernova theory and experiments
- High-order conservative reconstruction schemes for finite volume methods in cylindrical and spherical coordinates
- CWENO: uniformly accurate reconstructions for balance laws
- Well-Balanced Schemes for the Euler Equations with Gravitation: Conservative Formulation Using Global Fluxes
- Effect of Rotation on the Stability of a Stalled Cylindrical Shock and its Consequences for Core-Collapse Supernovae
- An advanced leakage scheme for neutrino treatment in astrophysical simulations
- Arbitrary order finite volume well-balanced schemes for the Euler equations with gravity
- Well balanced Arbitrary-Lagrangian-Eulerian finite volume schemes on moving nonconforming meshes for the Euler equations of gasdynamics with gravity
- New insights on the spin-up of a neutron star during core-collapse
- Second order finite volume scheme for Euler equations with gravity which is well-balanced for general equations of state and grid systems
- A well-balanced scheme for the simulation tool-kit A-MaZe: implementation, tests, and first applications to stellar structure
- Universality in the run-up of shock waves to the surface of a star