- Journal Home
- Volume 38 - 2025
- Volume 37 - 2025
- Volume 36 - 2024
- Volume 35 - 2024
- Volume 34 - 2023
- Volume 33 - 2023
- Volume 32 - 2022
- Volume 31 - 2022
- Volume 30 - 2021
- Volume 29 - 2021
- Volume 28 - 2020
- Volume 27 - 2020
- Volume 26 - 2019
- Volume 25 - 2019
- Volume 24 - 2018
- Volume 23 - 2018
- Volume 22 - 2017
- Volume 21 - 2017
- Volume 20 - 2016
- Volume 19 - 2016
- Volume 18 - 2015
- Volume 17 - 2015
- Volume 16 - 2014
- Volume 15 - 2014
- Volume 14 - 2013
- Volume 13 - 2013
- Volume 12 - 2012
- Volume 11 - 2012
- Volume 10 - 2011
- Volume 9 - 2011
- Volume 8 - 2010
- Volume 7 - 2010
- Volume 6 - 2009
- Volume 5 - 2009
- Volume 4 - 2008
- Volume 3 - 2008
- Volume 2 - 2007
- Volume 1 - 2006
Commun. Comput. Phys., 38 (2025), pp. 1130-1172.
Published online: 2025-09
Cited by
- BibTex
- RIS
- TXT
In this paper, we propose a high order semi-implicit well-balanced finite difference scheme for all-Mach Euler equations with a gravitational source. We start with the conservative form of full compressible Euler equations including conservation of total energy, for shock capturing in the high Mach regime. For asymptotic preserving in the low Mach regime and to address the difficulty of strong coupling between the stiff gravitational source and conservative variables, we add the evolution equation of the perturbation of potential temperature, which corresponds to weak potential temperature stratification under a hydrostatic background potential temperature. The resulting system is then split into a (non-stiff) nonlinear low dynamic material wave to be treated explicitly, and (stiff) fast acoustic and gravity waves to be treated implicitly. With the aid of explicit time evolution for the perturbation of potential temperature, we design a novel well-balanced finite difference weighted essentially non-oscillatory (WENO) scheme, which can be proven to be both asymptotic preserving and asymptotically accurate in the incompressible limit. Numerical experiments are provided to validate these properties.
}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.OA-2023-0118}, url = {http://global-sci.org/intro/article_detail/cicp/24355.html} }In this paper, we propose a high order semi-implicit well-balanced finite difference scheme for all-Mach Euler equations with a gravitational source. We start with the conservative form of full compressible Euler equations including conservation of total energy, for shock capturing in the high Mach regime. For asymptotic preserving in the low Mach regime and to address the difficulty of strong coupling between the stiff gravitational source and conservative variables, we add the evolution equation of the perturbation of potential temperature, which corresponds to weak potential temperature stratification under a hydrostatic background potential temperature. The resulting system is then split into a (non-stiff) nonlinear low dynamic material wave to be treated explicitly, and (stiff) fast acoustic and gravity waves to be treated implicitly. With the aid of explicit time evolution for the perturbation of potential temperature, we design a novel well-balanced finite difference weighted essentially non-oscillatory (WENO) scheme, which can be proven to be both asymptotic preserving and asymptotically accurate in the incompressible limit. Numerical experiments are provided to validate these properties.