- Journal Article
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. Show more
Journal / seriesJournal of Computational Physics
Pages / Article No.
SubjectNumerical methods; Hydrodynamics; Source terms; Well-balanced schemes
Organisational unit03851 - Mishra, Siddhartha / Mishra, Siddhartha
169631 - Structure preserving numerical methods for astrophysics (SNF)
MoreShow all metadata