Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers

Edelmann, P. F.; Horst, L.; Berberich, J. P.; Andrassy, R.; Higl, J.; Leidi, G.; Klingenberg, C.; Röpke, F. K.

doi:10.1051/0004-6361/202140653

Title: Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers

Journal Article · 11 August 2021 · Astronomy and Astrophysics

DOI: https://doi.org/10.1051/0004-6361/202140653 · OSTI ID:1873334

^[1];

^[2];

^[3]; Andrassy, R. ^[2];

^[2];

^[4];

^[3];

^[4]

Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Heidelberger Institut für Theoretische Studien, Heidelberg (Germany); Wurzburg University (Germany)
Heidelberger Institut für Theoretische Studien, Heidelberg (Germany)
Wurzburg University (Germany)
Heidelberger Institut für Theoretische Studien, Heidelberg (Germany); Heidelberg Univ. (Germany)

We report accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are merely tiny perturbations on top of a close balance between gravity and the pressure gradient, such simulations place heavy demands on numerical hydrodynamics schemes. We demonstrate how discretization errors on grids of reasonable size can lead to spurious flows orders of magnitude faster than the physical flow. Well-balanced numerical schemes can deal with this problem. Three such schemes were applied in the implicit, finite-volume SEVEN-LEAGUE HYDRO code in combination with a low-Mach-number numerical flux function. We compare how the schemes perform in four numerical experiments addressing some of the challenges imposed by typical problems in stellar hydrodynamics. We find that the α-β and deviation well-balancing methods can accurately maintain hydrostatic solutions provided that gravitational potential energy is included in the total energy balance. They accurately conserve minuscule entropy fluctuations advected in an isentropic stratification, which enables the methods to reproduce the expected scaling of convective flow speed with the heating rate. The deviation method also substantially increases accuracy of maintaining stationary orbital motions in a Keplerian disk on long timescales. The Cargo–LeRoux method fares substantially worse in our tests, although its simplicity may still offer some merits in certain situations. Overall, we find the well-balanced treatment of gravity in combination with low Mach number flux functions essential to reproducing correct physical solutions to challenging stellar slow-flow problems on affordable collocated grids.

View Accepted Manuscript (DOE)

Cite

Export

Save

Research Organization:: Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)

Sponsoring Organization:: German Research Foundation (DFG); Klaus Tschira Foundation; USDOE National Nuclear Security Administration (NNSA)

Grant/Contract Number:: 89233218CNA000001

OSTI ID:: 1873334

Report Number(s):: LA-UR-21-21056

Journal Information:: Astronomy and Astrophysics, Journal Name: Astronomy and Astrophysics Vol. 652; ISSN 0004-6361

Publisher:: EDP SciencesCopyright Statement

Country of Publication:: United States

Language:: English

References (70)

A numerical model for the flooding and drying of irregular domains: ID AND 2D MODELLING OF SHALLOW WATER FLOW Brufau, P.; Vázquez-Cendón, M. E.; García-Navarro, P. International Journal for Numerical Methods in Fluids, Vol. 39, Issue 3 https://doi.org/10.1002/fld.285	journal	May 2002
L ² Roe: a low dissipation version of Roe's approximate Riemann solver for low Mach numbers : L Oßwald, K.; Siegmund, A.; Birken, P. International Journal for Numerical Methods in Fluids, Vol. 81, Issue 2 https://doi.org/10.1002/fld.4175	journal	September 2015
A well-balanced scheme to capture non-explicit steady states in the Euler equations with gravity: WELL-BALANCED SCHEME FOR THE EULER EQUATIONS WITH GRAVITY Desveaux, Vivien; Zenk, Markus; Berthon, Christophe International Journal for Numerical Methods in Fluids, Vol. 81, Issue 2 https://doi.org/10.1002/fld.4177	journal	October 2015
Third‐ and fourth‐order well‐balanced schemes for the shallow water equations based on the CWENO reconstruction Castro, Manuel J.; Semplice, Matteo International Journal for Numerical Methods in Fluids, Vol. 89, Issue 8 https://doi.org/10.1002/fld.4700	journal	December 2018
A New Flux Splitting Scheme Liou, Meng-Sing; Steffen, Christopher J. Journal of Computational Physics, Vol. 107, Issue 1 https://doi.org/10.1006/jcph.1993.1122	journal	July 1993
A Sequel to AUSM: AUSM+ Liou, Meng-Sing Journal of Computational Physics, Vol. 129, Issue 2 https://doi.org/10.1006/jcph.1996.0256	journal	December 1996
Balancing Source Terms and Flux Gradients in High-Resolution Godunov Methods: The Quasi-Steady Wave-Propagation Algorithm LeVeque, Randall J. Journal of Computational Physics, Vol. 146, Issue 1 https://doi.org/10.1006/jcph.1998.6058	journal	October 1998
Transparent Boundary Conditions for MHD Simulations in Stratified Atmospheres Dedner, A.; Kröner, D.; Sofronov, I. L. Journal of Computational Physics, Vol. 171, Issue 2 https://doi.org/10.1006/jcph.2001.6779	journal	August 2001
Physics, Formation and Evolution of Rotating Stars Maeder, André Astronomy and Astrophysics Library https://doi.org/10.1007/978-3-540-76949-1	book	December 2008
Stellar Structure and Evolution Kippenhahn, Rudolf; Weigert, Alfred; Weiss, Achim Astronomy and Astrophysics Library https://doi.org/10.1007/978-3-642-30304-3	book	January 2012
�ber die partiellen Differenzengleichungen der mathematischen Physik Courant, R.; Friedrichs, K.; Lewy, H. Mathematische Annalen, Vol. 100, Issue 1 https://doi.org/10.1007/BF01448839	journal	December 1928
Riemann Solvers and Numerical Methods for Fluid Dynamics Toro, Eleuterio F. https://doi.org/10.1007/b79761	book	January 2009
A Numerical Scheme for the Compressible Low-Mach Number Regime of Ideal Fluid Dynamics Barsukow, Wasilij; Edelmann, Philipp V. F.; Klingenberg, Christian Journal of Scientific Computing, Vol. 72, Issue 2 https://doi.org/10.1007/s10915-017-0372-4	journal	January 2017
Approximate Riemann solvers, parameter vectors, and difference schemes Roe, P. L. Journal of Computational Physics, Vol. 43, Issue 2 https://doi.org/10.1016/0021-9991(81)90128-5	journal	October 1981
Preconditioned methods for solving the incompressible and low speed compressible equations Turkel, Eli Journal of Computational Physics, Vol. 72, Issue 2 https://doi.org/10.1016/0021-9991(87)90084-2	journal	October 1987
Efficient implementation of essentially non-oscillatory shock-capturing schemes Shu, Chi-Wang; Osher, Stanley Journal of Computational Physics, Vol. 77, Issue 2 https://doi.org/10.1016/0021-9991(88)90177-5	journal	August 1988
Upwind methods for hyperbolic conservation laws with source terms Bermudez, Alfredo; Vazquez, Ma Elena Computers & Fluids, Vol. 23, Issue 8 https://doi.org/10.1016/0045-7930(94)90004-3	journal	November 1994
Analysis and implementation of TR-BDF2 Hosea, M. E.; Shampine, L. F. Applied Numerical Mathematics, Vol. 20, Issue 1-2 https://doi.org/10.1016/0168-9274(95)00115-8	journal	February 1996
On the behaviour of upwind schemes in the low Mach number limit Guillard, Hervé; Viozat, Cécile Computers & Fluids, Vol. 28, Issue 1 https://doi.org/10.1016/S0045-7930(98)00017-6	journal	January 1999
Well-balanced central finite volume methods for the Ripa system Touma, R.; Klingenberg, C. Applied Numerical Mathematics, Vol. 97 https://doi.org/10.1016/j.apnum.2015.07.001	journal	November 2015
On the behavior of upwind schemes in the low Mach number limit: II. Godunov type schemes Guillard, Hervé; Murrone, Angelo Computers & Fluids, Vol. 33, Issue 4 https://doi.org/10.1016/j.compfluid.2003.07.001	journal	May 2004
High order well-balanced finite volume methods for multi-dimensional systems of hyperbolic balance laws Berberich, Jonas P.; Chandrashekar, Praveen; Klingenberg, Christian Computers & Fluids, Vol. 219 https://doi.org/10.1016/j.compfluid.2021.104858	journal	April 2021
A sequel to AUSM, Part II: AUSM+-up for all speeds Liou, Meng-Sing Journal of Computational Physics, Vol. 214, Issue 1 https://doi.org/10.1016/j.jcp.2005.09.020	journal	May 2006
An All-Speed Roe-type scheme and its asymptotic analysis of low Mach number behaviour Li, Xue-song; Gu, Chun-wei Journal of Computational Physics, Vol. 227, Issue 10 https://doi.org/10.1016/j.jcp.2008.01.037	journal	May 2008
A low-Mach number fix for Roe’s approximate Riemann solver Rieper, Felix Journal of Computational Physics, Vol. 230, Issue 13 https://doi.org/10.1016/j.jcp.2011.03.025	journal	June 2011
Well-balanced schemes for the Euler equations with gravitation Käppeli, R.; Mishra, S. Journal of Computational Physics, Vol. 259 https://doi.org/10.1016/j.jcp.2013.11.028	journal	February 2014
High-order well-balanced finite volume schemes for the Euler equations with gravitation Grosheintz-Laval, L.; Käppeli, R. Journal of Computational Physics, Vol. 378 https://doi.org/10.1016/j.jcp.2018.11.018	journal	February 2019
Internal Gravity Waves Sutherland, Bruce R. https://doi.org/10.1017/CBO9780511780318	book	October 2010
On multigrid solution of the implicit equations of hydrodynamics: Experiments for the compressible Euler equations in general coordinates Kifonidis, K.; Müller, E. Astronomy & Astrophysics, Vol. 544 https://doi.org/10.1051/0004-6361/201116979	journal	July 2012
New numerical solver for flows at various Mach numbers Miczek, F.; Röpke, F. K.; Edelmann, P. V. F. Astronomy & Astrophysics, Vol. 576 https://doi.org/10.1051/0004-6361/201425059	journal	March 2015
A well-balanced finite volume scheme for the Euler equations with gravitation: The exact preservation of hydrostatic equilibrium with arbitrary entropy stratification Käppeli, R.; Mishra, S. Astronomy & Astrophysics, Vol. 587 https://doi.org/10.1051/0004-6361/201527815	journal	February 2016
Spherical-shell boundaries for two-dimensional compressible convection in a star Pratt, J.; Baraffe, I.; Goffrey, T. Astronomy & Astrophysics, Vol. 593 https://doi.org/10.1051/0004-6361/201628296	journal	September 2016
Benchmarking the Multidimensional Stellar Implicit Code MUSIC Goffrey, T.; Pratt, J.; Viallet, M. Astronomy & Astrophysics, Vol. 600 https://doi.org/10.1051/0004-6361/201628960	journal	March 2017
Testing a one-dimensional prescription of dynamical shear mixing with a two-dimensional hydrodynamic simulation Edelmann, P. V. F.; Röpke, F. K.; Hirschi, R. Astronomy & Astrophysics, Vol. 604 https://doi.org/10.1051/0004-6361/201629873	journal	July 2017
A well-balanced scheme for the simulation tool-kit A-MaZe: implementation, tests, and first applications to stellar structure Popov, M. V.; Walder, R.; Folini, D. Astronomy & Astrophysics, Vol. 630 https://doi.org/10.1051/0004-6361/201834180	journal	October 2019
Comparison of 2D and 3D compressible convection in a pre-main sequence star Pratt, J.; Baraffe, I.; Goffrey, T. Astronomy & Astrophysics, Vol. 638 https://doi.org/10.1051/0004-6361/201834736	journal	June 2020
Fully compressible simulations of waves and core convection in main-sequence stars Horst, L.; Edelmann, P. V. F.; Andrássy, R. Astronomy & Astrophysics, Vol. 641 https://doi.org/10.1051/0004-6361/202037531	journal	September 2020
Calibrating core overshooting parameters with two-dimensional hydrodynamical simulations Higl, J.; Müller, E.; Weiss, A. Astronomy & Astrophysics, Vol. 646 https://doi.org/10.1051/0004-6361/202039532	journal	February 2021
Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star Horst, L.; Hirschi, R.; Edelmann, P. V. F. Astronomy & Astrophysics, Vol. 653 https://doi.org/10.1051/0004-6361/202140825	journal	September 2021
A low-Mach Roe-type solver for the Euler equations allowing for gravity source terms Barsukow, Wasilij; Edelmann, Philipp V. F.; Klingenberg, Christian ESAIM: Proceedings and Surveys, Vol. 58 https://doi.org/10.1051/proc/201758027	journal	January 2017
Riemann Solvers for some Hyperbolic Problems with a Source Term Leroux, Alain Yves ESAIM: Proceedings, Vol. 6 https://doi.org/10.1051/proc:1999047	journal	January 1999
The Accuracy, Consistency, and Speed of an Electron‐Positron Equation of State Based on Table Interpolation of the Helmholtz Free Energy Timmes, F. X.; Swesty, F. Douglas The Astrophysical Journal Supplement Series, Vol. 126, Issue 2 https://doi.org/10.1086/313304	journal	February 2000
Mapping Initial Hydrostatic Models in Godunov Codes Zingale, M.; Dursi, L. J.; ZuHone, J. The Astrophysical Journal Supplement Series, Vol. 143, Issue 2 https://doi.org/10.1086/342754	journal	December 2002
Simulations of Core Convection in Rotating A‐Type Stars: Differential Rotation and Overshooting Browning, Matthew K.; Brun, Allan Sacha; Toomre, Juri The Astrophysical Journal, Vol. 601, Issue 1 https://doi.org/10.1086/380198	journal	January 2004
Low Mach Number Modeling of Type Ia Supernovae. I. Hydrodynamics Almgren, A. S.; Bell, J. B.; Rendleman, C. A. The Astrophysical Journal, Vol. 637, Issue 2 https://doi.org/10.1086/498426	journal	February 2006
Active Carbon and Oxygen Shell Burning Hydrodynamics Meakin, Casey A.; Arnett, David The Astrophysical Journal, Vol. 637, Issue 1 https://doi.org/10.1086/500544	journal	January 2006
Turbulent Convection in Stellar Interiors. I. Hydrodynamic Simulation Meakin, Casey A.; Arnett, David The Astrophysical Journal, Vol. 667, Issue 1 https://doi.org/10.1086/520318	journal	September 2007
Turbulent Convection in Stellar Interiors. iii. Mean-Field Analysis and Stratification Effects Viallet, Maxime; Meakin, Casey; Arnett, David The Astrophysical Journal, Vol. 769, Issue 1 https://doi.org/10.1088/0004-637X/769/1/1	journal	April 2013
Internal Gravity Waves in Massive Stars: Angular Momentum Transport Rogers, T. M.; Lin, D. N. C.; McElwaine, J. N. The Astrophysical Journal, Vol. 772, Issue 1 https://doi.org/10.1088/0004-637X/772/1/21	journal	July 2013
HYDRODYNAMIC SIMULATIONS OF H ENTRAINMENT AT THE TOP OF He-SHELL FLASH CONVECTION Woodward, Paul R.; Herwig, Falk; Lin, Pei-Hung The Astrophysical Journal, Vol. 798, Issue 1 https://doi.org/10.1088/0004-637X/798/1/49	journal	December 2014
Fish: a Three-Dimensional Parallel Magnetohydrodynamics code for Astrophysical Applications Käppeli, R.; Whitehouse, S. C.; Scheidegger, S. The Astrophysical Journal Supplement Series, Vol. 195, Issue 2 https://doi.org/10.1088/0067-0049/195/2/20	journal	August 2011
Well-balanced schemes to capture non-explicit steady states: Ripa model Desveaux, Vivien; Zenk, Markus; Berthon, Christophe Mathematics of Computation, Vol. 85, Issue 300 https://doi.org/10.1090/mcom/3069	journal	January 2016
Idealized hydrodynamic simulations of turbulent oxygen-burning shell convection in 4π geometry Jones, S.; Andrassy, R.; Sandalski, S. Monthly Notices of the Royal Astronomical Society, Vol. 465, Issue 3 https://doi.org/10.1093/mnras/stw2783	journal	October 2016
3D hydrodynamic simulations of carbon burning in massive stars Cristini, A.; Meakin, C.; Hirschi, R. Monthly Notices of the Royal Astronomical Society, Vol. 471, Issue 1 https://doi.org/10.1093/mnras/stx1535	journal	June 2017
Well-balanced Arbitrary-Lagrangian-Eulerian finite volume schemes on moving nonconforming meshes for the Euler equations of gas dynamics with gravity Gaburro, Elena; Castro, Manuel J.; Dumbser, Michael Monthly Notices of the Royal Astronomical Society, Vol. 477, Issue 2 https://doi.org/10.1093/mnras/sty542	journal	March 2018
3D hydrodynamic simulations of C ingestion into a convective O shell Andrassy, R.; Herwig, F.; Woodward, P. Monthly Notices of the Royal Astronomical Society, Vol. 491, Issue 1 https://doi.org/10.1093/mnras/stz2952	journal	October 2019
Dependence of convective boundary mixing on boundary properties and turbulence strength Cristini, A.; Hirschi, R.; Meakin, C. Monthly Notices of the Royal Astronomical Society, Vol. 484, Issue 4 https://doi.org/10.1093/mnras/stz312	journal	February 2019
Logically Rectangular Grids and Finite Volume Methods for PDEs in Circular and Spherical Domains Calhoun, Donna A.; Helzel, Christiane; LeVeque, Randall J. SIAM Review, Vol. 50, Issue 4 https://doi.org/10.1137/060664094	journal	January 2008
A Second Order Well-Balanced Finite Volume Scheme for Euler Equations with Gravity Chandrashekar, Praveen; Klingenberg, Christian SIAM Journal on Scientific Computing, Vol. 37, Issue 3 https://doi.org/10.1137/140984373	journal	January 2015
Well-Balanced Unstaggered Central Schemes for the Euler Equations with Gravitation Touma, R.; Koley, U.; Klingenberg, C. SIAM Journal on Scientific Computing, Vol. 38, Issue 5 https://doi.org/10.1137/140992667	journal	January 2016
A Fast and Stable Well-Balanced Scheme with Hydrostatic Reconstruction for Shallow Water Flows Audusse, Emmanuel; Bouchut, François; Bristeau, Marie-Odile SIAM Journal on Scientific Computing, Vol. 25, Issue 6 https://doi.org/10.1137/S1064827503431090	journal	January 2004
The Evolution of Rotating Stars Maeder, André; Meynet, Georges Annual Review of Astronomy and Astrophysics, Vol. 38, Issue 1 https://doi.org/10.1146/annurev.astro.38.1.143	journal	September 2000
Low-Diffusion Flux-Splitting Methods for Flows at All Speeds Edwards, Jack R.; Liou, Meng-Sing AIAA Journal, Vol. 36, Issue 9 https://doi.org/10.2514/2.587	journal	September 1998
An Advanced Leakage Scheme for Neutrino Treatment in Astrophysical Simulations Perego, A.; Cabezón, R. M.; Käppeli, R. The Astrophysical Journal Supplement Series, Vol. 223, Issue 2 https://doi.org/10.3847/0067-0049/223/2/22	journal	April 2016
The last Minutes of Oxygen Shell Burning in a Massive star Müller, Bernhard; Viallet, Maxime; Heger, Alexander The Astrophysical Journal, Vol. 833, Issue 1 https://doi.org/10.3847/1538-4357/833/1/124	journal	December 2016
A High-performance and Portable All-Mach Regime Flow Solver Code with Well-balanced Gravity. Application to Compressible Convection Padioleau, Thomas; Tremblin, Pascal; Audit, Edouard The Astrophysical Journal, Vol. 875, Issue 2 https://doi.org/10.3847/1538-4357/ab0f2c	journal	April 2019
Three-dimensional Simulations of Massive Stars. I. Wave Generation and Propagation Edelmann, P. V. F.; Ratnasingam, R. P.; Pedersen, M. G. The Astrophysical Journal, Vol. 876, Issue 1 https://doi.org/10.3847/1538-4357/ab12df	journal	April 2019
An Extension of the Athena++ Framework for Fully Conservative Self-gravitating Hydrodynamics Mullen, P. D.; Hanawa, Tomoyuki; Gammie, C. F. The Astrophysical Journal Supplement Series, Vol. 252, Issue 2 https://doi.org/10.3847/1538-4365/abcfbd	journal	February 2021
Capturing Near-Equilibrium Solutions: A Comparison between High-Order Discontinuous Galerkin Methods and Well-Balanced Schemes Veiga, Maria Han Communications in Computational Physics, Vol. 26, Issue 1 https://doi.org/10.4208/cicp.OA-2018-0071	journal	February 2019
Second Order Finite Volume Scheme for Euler Equations with Gravity which is Well-Balanced for General Equations of State and Grid Systems Berberich, Jonas P. Communications in Computational Physics, Vol. 26, Issue 2 https://doi.org/10.4208/cicp.OA-2018-0152	journal	February 2019

Similar Records

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

Journal Article · 2024 · Astronomy and Astrophysics · OSTI ID:2367416

Leidi, G.; Andrassy, R.; Barsukow, W.; +3 more

Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star

Journal Article · 2021 · Astronomy and Astrophysics · OSTI ID:1833258

Horst, Leonard; Hirschi, Raphael; Edelmann, Philipp Valentin Ferdinand; +2 more

A finite-volume scheme for modeling compressible magnetohydrodynamic flows at low Mach numbers in stellar interiors

Journal Article · 2022 · Astronomy and Astrophysics · OSTI ID:1922024

Leidi, Giovanni; Birke, Claudius; Andrassy, Robert; +5 more

Related Subjects

79 ASTRONOMY AND ASTROPHYSICS
convection
hydrodynamics
methods: numerical

Title: Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers

Citation Formats

References (70)

Similar Records

Related Subjects