Fornax: a Flexible Code for Multiphysics Astrophysical Simulations

Abstract

This paper describes the design and implementation of our new multi-group, multi-dimensional radiation hydrodynamics (RHD) code Fornax and provides a suite of code tests to validate its application in a wide range of physical regimes. Instead of focusing exclusively on tests of neutrino radiation hydrodynamics relevant to the core-collapse supernova problem for which Fornax is primarily intended, we present here classical and rigorous demonstrations of code performance relevant to a broad range of multi-dimensional hydrodynamic and multi-group radiation hydrodynamic problems. Our code solves the comoving-frame radiation moment equations using the M1 closure, utilizes conservative high-order reconstruction, employs semi-explicit matter and radiation transport via a high-order time stepping scheme, and is suitable for application to a wide range of astrophysical problems. To this end, we first describe the philosophy, algorithms, and methodologies of Fornax and then perform numerous stringent code tests, that collectively and vigorously exercise the code, demonstrate the excellent numerical fidelity with which it captures the many physical effects of radiation hydrodynamics, and show excellent strong scaling well above 100k MPI tasks.

arXiv commentAccepted to the Astrophysical Journal Supplement Series; A few more textual and reference updates; As before, one additional code test included

References in corpus (43)

• Athena: A New Code for Astrophysical MHD
• Composition and thermodynamics of nuclear matter with light clusters
• Core-collapse supernova equations of state based on neutron star observations
• Nucleosynthesis and Remnants in Massive Stars of Solar Metallicity
• Statistical Model for a Complete Supernova Equation of State
• Perspectives on Core-Collapse Supernova Theory
• Delayed neutrino-driven supernova explosions aided by the standing accretion-shock instability
• Radiation hydrodynamics with neutrinos: Variable Eddington factor method for core-collapse supernova simulations
• Shock Breakout in Core-Collapse Supernovae and its Neutrino Signature
• A Finite Difference Representation of Neutrino Radiation Hydrodynamics in Spherically Symmetric General Relativistic Space-Time
• Exploring the relativistic regime with Newtonian hydrodynamics: An improved effective gravitational potential for supernova simulations
• Truncated Moment Formalism for Radiation Hydrodynamics in Numerical Relativity
• Criteria for Core-Collapse Supernova Explosions by the Neutrino Mechanism
• A new multi-dimensional general relativistic neutrino hydrodynamics code for core-collapse supernovae. I. Method and code tests in spherical symmetry
• General-Relativistic Simulations of Three-Dimensional Core-Collapse Supernovae
• The Development of Explosions in Axisymmetric Ab Initio Core-Collapse Supernova Simulations of 12-25 $M_{\odot }$ Stars
• Equations and Algorithms for Mixed Frame Flux-Limited Diffusion Radiation Hydrodynamics
• A new multidimensional, energy-dependent two-moment transport code for neutrino-hydrodynamics
• General Relativistic Three-Dimensional Multi-Group Neutrino Radiation-Hydrodynamics Simulations of Core-Collapse Supernovae
• High-order conservative reconstruction schemes for finite volume methods in cylindrical and spherical coordinates
• High-Resolution Three-Dimensional Simulations of Core-Collapse Supernovae in Multiple Progenitors
• Crucial Physical Dependencies of the Core-Collapse Supernova Mechanism
• Fully General Relativistic Simulations of Core-Collapse Supernovae with An Approximate Neutrino Transport
• Two Dimensional Core-Collapse Supernova Explosions Aided by General Relativity with Multidimensional Neutrino Transport
• Simulating Self-gravitating Hydrodynamic Flows
• On the Impact of Three Dimensions in Simulations of Neutrino-Driven Core-Collapse Supernova Explosions
• Simulations of core-collapse supernovae in spatial axisymmetry with full Boltzmann neutrino transport
• A Validated Nonlinear Kelvin-Helmholtz Benchmark for Numerical Hydrodynamics
• Correlated Gravitational Wave and Neutrino Signals from General-Relativistic Rapidly Rotating Iron Core Collapse
• Revival of the Fittest: Exploding Core-Collapse Supernovae from 12 to 25 M${}_{\odot }$
• Numerical Relativity in Spherical Polar Coordinates: Evolution Calculations with the BSSN Formulation
• CASTRO: A New Compressible Astrophysical Solver. II. Gray Radiation Hydrodynamics
• A numerical model for multigroup radiation hydrodynamics
• CASTRO: A New Compressible Astrophysical Solver. III. Multigroup Radiation Hydrodynamics
• Should One Use the Ray-by-Ray Approximation in Core-Collapse Supernova Simulations?
• Turbulence in Core-Collapse Supernovae
• Radiation hydrodynamics integrated in the code PLUTO
• A New Algorithm for 2-D Transport for Astrophysical Simulations: I. General Formulation and Tests for the 1-D Spherical Case
• 2D and 3D Core-Collapse Supernovae Simulation Results Obtained with the CHIMERA Code
• Neutrino Signals of Core-Collapse Supernovae in Underground Detectors
• BETHE-Hydro: An Arbitrary Lagrangian-Eulerian Multi-dimensional Hydrodynamics Code for Astrophysical Simulations
• Detecting the Supernova Breakout Burst in Terrestrial Neutrino Detectors
• A Magellan-IMACS-IFU Search for Dynamical Drivers of Nuclear Activity. I. Reduction Pipeline and Galaxy Catalog

Cited by in corpus (89)

• Core-Collapse Supernova Explosion Theory
• The Overarching Framework of Core-Collapse Supernova Explosions as Revealed by 3D Fornax Simulations
• Characterizing the Gravitational Wave Signal from Core-Collapse Supernovae
• Hydrodynamics of core-collapse supernovae and their progenitors
• Three-Dimensional Supernova Explosion Simulations of 9-, 10-, 11-, 12-, and 13-M${}_{\odot }$ Stars
• Three-dimensional Models of Core-collapse Supernovae From Low-mass Progenitors With Implications for Crab
• Neutrino oscillations in supernovae: angular moments and fast instabilities
• Temporal and Angular Variations of 3D Core-Collapse Supernova Emissions and their Physical Correlations
• Simulating Turbulence-aided Neutrino-driven Core-collapse Supernova Explosions in One Dimension
• H-AMR: A New GPU-accelerated GRMHD Code for Exascale Computing With 3D Adaptive Mesh Refinement and Local Adaptive Time-stepping
• A systematic study of proto-neutron star convection in three-dimensional core-collapse supernova simulations
• Full Transport General Relativistic Radiation Magnetohydrodynamics for Nucleosynthesis in Collapsars
• Supernova Model Discrimination with Hyper-Kamiokande
• Core-collapse supernova neutrino emission and detection informed by state-of-the-art three-dimensional numerical models
• Neutrino transport in general relativistic neutron star merger simulations
• Neutrino Emission as Diagnostics of Core-Collapse Supernovae
• Neutrino emission characteristics of black hole formation in three-dimensional simulations of stellar collapse
• An Implicit Finite Volume Scheme to Solve the Time Dependent Radiation Transport Equation Based on Discrete Ordinates
• The Collapse and Three-Dimensional Explosion of Three-Dimensional, vis à vis One-Dimensional, Massive-star Supernova Progenitor Models
• Gravitational Waves from Neutrino Asymmetries in Core-Collapse Supernovae
• Towards an Understanding of the Resolution Dependence of Core-Collapse Supernova Simulations
• Chimera: A massively parallel code for core-collapse supernova simulation
• The Essential Character of the Neutrino Mechanism of Core-Collapse Supernova Explosions
• Binary-Stripped Stars as Core-Collapse Supernovae Progenitors
• Kicks and Induced Spins of Neutron Stars at Birth
• Neutrinos and nucleosynthesis of elements
• Supernova neutrino signals based on long-term axisymmetric simulations
• Parthenon -- a performance portable block-structured adaptive mesh refinement framework
• $\nu \text{bhlight}$: Radiation GRMHD for Neutrino-Driven Accretion Flows
• Numerical relativity in spherical coordinates: A new dynamical spacetime and general relativistic MHD evolution framework for the Einstein Toolkit
• Gmunu: Paralleled, grid-adaptive, general-relativistic magnetohydrodynamics in curvilinear geometries in dynamical spacetimes
• Beyond moments: relativistic Lattice-Boltzmann methods for radiative transport in computational astrophysics
• Resolution Study for Three-dimensional Supernova Simulations with the Prometheus-Vertex Code
• New method for detecting fast neutrino flavor conversions in core-collapse supernova models with two-moment neutrino transport
• Constructing angular distributions of neutrinos in core collapse supernova from zero-th and first moments calibrated by full Boltzmann neutrino transport
• Applications of Machine Learning to Predicting Core-collapse Supernova Explosion Outcomes
• LSTM and CNN application for core-collapse supernova search in gravitational wave real data
• Quokka: A code for two-moment AMR radiation hydrodynamics on GPUs
• Multi-Frequency General Relativistic Radiation-Hydrodynamics with ${\mathbf{M}}_1$ Closure
• Retrieval of energy spectra for all flavor of neutrinos from core-collapse supernova with multiple detectors
• Two-Moment Neutrino Flavor Transformation with applications to the Fast Flavor Instability in Neutron Star Mergers
• Convection-Aided Explosions in One-Dimensional Core-Collapse Supernova Simulations I: Technique and Validation
• A 3D Simulation of a Type II-P Supernova: from Core Bounce to Beyond Shock Breakout
• A Novel Multi-Dimensional Boltzmann Neutrino Transport Scheme for Core-Collapse Supernovae
• VETTAM: A scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method
• Three approaches for the classification of protoneutron star oscillation modes
• General-relativistic radiation transport scheme in $\text{Gmunu}$ I: Implementation of two-moment based multi-frequency radiative transfer and code tests
• Addition of tabulated equation of state and neutrino leakage support to IllinoisGRMHD
• MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics
• Efficient estimation method for time evolution of proto-neutron star mass and radius from supernova neutrino signal
• NADA-FLD: A General Relativistic, Multi-dimensional Neutrino-hydrodynamics Code Employing Flux-limited Diffusion
• Nucleosynthesis in Outflows from Black Hole-Neutron Star Merger Disks With Full GR$\nu$RMHD
• A DG-IMEX method for two-moment neutrino transport: Nonlinear solvers for neutrino-matter coupling
• Measuring the properties of $f-$mode oscillations of a protoneutron star by third generation gravitational-wave detectors
• The Guided Moments formalism: a new efficient full-neutrino treatment for astrophysical simulations
• Radiation transport methods in star formation simulations
• Smoothed Particle Radiation Hydrodynamics: Two-Moment method with Local Eddington Tensor Closure
• A Finite Element Method for Angular Discretization of the Radiation Transport Equation on Spherical Geodesic Grids
• Consistent perturbative modeling of pseudo-Newtonian core-collapse supernova simulations
• General Relativistic Hydrodynamics on a Moving-mesh I: Static Spacetimes
• thornado-hydro: a discontinuous Galerkin method for supernova hydrodynamics with nuclear equations of state
• CAFE-R a code that solves the special relativistic radiation hydrodynamics equations
• Calibration of the Advanced Spectral Leakage scheme for neutron star merger simulations, and extension to smoothed-particle hydrodynamics
• A Force Explosion Condition for Spherically Symmetric Core-collapse Supernovae
• Visualizing convolutional neural network for classifying gravitational waves from core-collapse supernovae
• An Asymptotically-Correct Implicit-Explicit Time Integration Scheme for Finite Volume Radiation-Hydrodynamics
• Angle-dependent in-situ fast flavor transformations in post-neutron star merger disks
• Relativistic radiation hydrodynamics in a reference-metric formulation
• Magnetorotational neutron star kicks
• Effects of Different Closure Choices in Core-Collapse Supernova Simulations
• Principal-Axis Analysis of the Eddington Tensor for the Early Post-Bounce Phase of Rotational Core-Collapse Supernovae
• Non-LTE radiation hydrodynamics in PLUTO
• Global high-order numerical schemes for the time evolution of the general relativistic radiation magneto-hydrodynamics equations
• Energy conserving and well-balanced discontinuous Galerkin methods for the Euler-Poisson equations in spherical symmetry
• Explosion Energies for Core-collapse Supernovae I: Analytic, Spherically Symmetric Solutions
• Dynamics of small grains in transitional discs
• Magnetic Field Strength Effects on Nucleosynthesis from Neutron Star Merger Outflows
• A New Residual Distribution Hydrodynamics Solver for Astrophysical Simulations
• Uncovering the neutrino mass ordering with the next galactic core-collapse supernova neutrino burst using water Cherenkov detectors
• A novel numerical method for mixed-frame multigroup radiation-hydrodynamics with GPU acceleration implemented in the QUOKKA code
• A Parametric Study of the SASI Comparing General Relativistic and Nonrelativistic Treatments
• Individual Neutrino Masses From a Supernova
• Where, when and why: occurrence of fast-pairwise collective neutrino oscillation in three-dimensional core-collapse supernova models
• A New Moment-Based General-Relativistic Neutrino-Radiation Transport Code: Methods and First Applications to Neutron Star Mergers
• On the Origin of Pulsar and Magnetar Magnetic Fields
• A Discontinuous Galerkin Method for General Relativistic Hydrodynamics in thornado
• GenASiS: General Astrophysical Simulation System. II. Self-gravitating Baryonic Matter
• Stochastic Gravitational-Wave Background from Stellar Core-Collapse Events
• A Black Hole is Born: 3D GRMHD Simulation of Black Hole Formation from Core-Collapse