A Jacobian-free Newton-Krylov method for time-implicit multidimensional hydrodynamics

Abstract

This work is a continuation of our efforts to develop an efficient implicit solver for multidimensional hydrodynamics for the purpose of studying important physical processes in stellar interiors, such as turbulent convection and overshooting. We present an implicit solver that results from the combination of a Jacobian-Free Newton-Krylov method and a preconditioning technique tailored to the inviscid, compressible equations of stellar hydrodynamics. We assess the accuracy and performance of the solver for both 2D and 3D problems for Mach numbers down to $10^{-6}$. Although our applications concern flows in stellar interiors, the method can be applied to general advection and/or diffusion-dominated flows. The method presented in this paper opens up new avenues in 3D modeling of realistic stellar interiors allowing the study of important problems in stellar structure and evolution.

Accepted for publication in A&A

References in corpus (5)

• Turbulent Convection in Stellar Interiors. I. Hydrodynamic Simulation
• Global non-spherical oscillations in 3-D $4\pi$ simulations of the H-ingestion flash
• Towards a new generation of multi-dimensional stellar evolution models: development of an implicit hydrodynamic code
• Multidimensional hydrodynamic simulations of the hydrogen injection flash
• Comparison of different nonlinear solvers for 2D time-implicit stellar hydrodynamics

Cited by in corpus (28)

• The Last Minutes of Oxygen Shell Burning in a Massive Star
• Lithium depletion in solar-like stars: effect of overshooting based on realistic multi-dimensional simulations
• Extreme value statistics for two-dimensional convective penetration in a pre-Main Sequence star
• Angular momentum fluctuations in the convective helium shell of massive stars
• A study of convective core overshooting as a function of stellar mass based on two-dimensional hydrodynamical simulations
• Benchmarking the Multi-dimensional Stellar Implicit Code MUSIC
• MAESTROeX: A Massively Parallel Low Mach Number Astrophysical Solver
• Spherical-shell boundaries for two-dimensional compressible convection in a star
• Two-dimensional simulations of solar-like models with artificially enhanced luminosity -- I. Impact on convective penetration
• Multi-dimensional structure of accreting young stars
• Comparison of two- and three-dimensional compressible convection in a pre-main sequence star
• Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes
• The KOSMA-$\tau$ PDR Model -- I. Recent updates to the numerical model of photo-dissociated regions
• Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star
• Two-dimensional simulations of solar-like models with artificially enhanced luminosity. II. Impact on internal gravity waves
• Multidimensional Simulations of Core Convection
• Two-dimensional simulations of internal gravity waves in a 5 $M_{\odot }$ Zero-Age-Main-Sequence model
• Accelerated evolution of convective simulations
• A well-balanced scheme for the simulation tool-kit A-MaZe: implementation, tests, and first applications to stellar structure
• Local heating due to convective overshooting and the solar modelling problem
• Impact of radial truncation on global 2D hydrodynamic simulations for a Sun-like model
• Convective shells in the interior of Cepheid variable stars: overshooting models based on hydrodynamic simulations
• Mixing by Internal Gravity Waves in Stars: Assessing Numerical Simulations Against Theory
• The shape of convection in 2D and 3D global simulations of stellar interiors
• Effects of stratification on overshooting and waves atop the convective core of $5M_{\odot }$ main-sequence stars
• Suppression of lithium depletion in young low-mass stars from fast rotation
• Modelling low Mach number stellar hydrodynamics with MAESTROeX
• Effects of radiative losses on the relativistic jets of high-mass microquasars