Research reports

Well-balanced finite volume schemes for nearly steady adiabatic flows

by L. Grosheintz-Laval and R. Käppeli

(Report number 2020-37)

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.

Keywords: Numerical methods Hydrodynamics Source terms Well-balanced scheme

BibTeX
@Techreport{GK20_910,
  author = {L. Grosheintz-Laval and R. K\"appeli},
  title = {Well-balanced finite volume schemes for nearly steady adiabatic flows},
  institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
  number = {2020-37},
  address = {Switzerland},
  url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2020/2020-37.pdf },
  year = {2020}
}

Disclaimer
© Copyright for documents on this server remains with the authors. Copies of these documents made by electronic or mechanical means including information storage and retrieval systems, may only be employed for personal use. The administrators respectfully request that authors inform them when any paper is published to avoid copyright infringement. Note that unauthorised copying of copyright material is illegal and may lead to prosecution. Neither the administrators nor the Seminar for Applied Mathematics (SAM) accept any liability in this respect. The most recent version of a SAM report may differ in formatting and style from published journal version. Do reference the published version if possible (see SAM Publications).

JavaScript has been disabled in your browser