Improved Coupling of Hydrodynamics and Nuclear Reactions via Spectral Deferred Corrections

Abstract

Simulations in stellar astrophysics involve the coupling of hydrodynamics and nuclear reactions under a wide variety of conditions, from simmering convective flows to explosive nucleosynthesis. Numerical techniques such as operator splitting (most notably Strang splitting) are usually employed to couple the physical processes, but this can affect the accuracy of the simulation, particularly when the burning is vigorous. Furthermore, Strang splitting does not have a straightforward extension to higher-order integration in time. We present a new temporal integration strategy based on spectral deferred corrections and describe the second- and fourth-order implementations in the open-source, finite-volume, compressible hydrodynamics code Castro. One notable advantage to these schemes is that they combine standard low-order discretizations for individual physical processes in a way that achieves an arbitrarily high order of accuracy. We demonstrate the improved accuracy of the new methods on several test problems of increasing complexity.

arXiv commentaccepted to ApJ; supplementary Jupyter/SymPy notebook with some derivations here: https://github.com/AMReX-Astro/supplementary_material/tree/master/castro_sdc

References in corpus (17)

• The NumPy array: a structure for efficient numerical computation
• A Multi-Code Analysis Toolkit for Astrophysical Simulation Data
• Turbulent Convection in Stellar Interiors. I. Hydrodynamic Simulation
• CASTRO: A New Compressible Astrophysical Solver. I. Hydrodynamics and Self-Gravity
• Mapping Initial Hydrostatic Models in Godunov Codes
• Beyond second-order convergence in simulations of magnetised binary neutron stars with realistic microphysics
• MAESTRO: An Adaptive Low Mach Number Hydrodynamics Algorithm for Stellar Flows
• A fourth-order accurate finite volume method for ideal MHD via upwind constrained transport
• High-Order Algorithms for Compressible Reacting Flow with Complex Chemistry
• White Dwarf Mergers on Adaptive Meshes I. Methodology and Code Verification
• A high-order spectral deferred correction strategy for low Mach number flow with complex chemistry
• A Fourth-Order Adaptive Mesh Refinement Algorithm for the Multicomponent, Reacting Compressible Navier-Stokes Equations
• On the Piecewise Parabolic Method for Compressible Flow with Stellar Equations of State
• Numerical Stability of Detonations in White Dwarf Simulations
• APSARA: A multi-dimensional unsplit fourth-order explicit Eulerian hydrodynamics code for arbitrary curvilinear grids
• pyro: a framework for hydrodynamics explorations and prototyping
• Toward Resolved Simulations of Burning Fronts in Thermonuclear X-ray Bursts

Cited by in corpus (10)

• Dynamics of Laterally Propagating Flames in X-ray Bursts. I. Burning Front Structure
• An Improved Method for Coupling Hydrodynamics with Astrophysical Reaction Networks
• Deep Neural Networks for Modeling Astrophysical Nuclear Reacting Flows
• Comparing Early Evolution of Flames in X-ray Bursts in Two and Three Dimensions
• SUNDIALS Time Integrators for Exascale Applications with Many Independent ODE Systems
• The Castro AMR Simulation Code: Current and Future Developments
• Performance of explicit and IMEX MRI multirate methods on complex reactive flow problems within modern parallel adaptive structured grid frameworks
• Toward First-Principles Multi-Messenger Predictions: Coupling Nuclear Networks with GR Radiation-MHD in {\tt Gmunu}
• Quantifying Advantages of a Moving Mesh in Nuclear Hydrodynamics
• Modelling low Mach number stellar hydrodynamics with MAESTROeX