A DG-IMEX method for two-moment neutrino transport: Nonlinear solvers for neutrino-matter coupling

Abstract

Neutrino-matter interactions play an important role in core-collapse supernova (CCSN) explosions as they contribute to both lepton number and/or four-momentum exchange between neutrinos and matter, and thus act as the agent for neutrino-driven explosions. Due to the multiscale nature of neutrino transport in CCSN simulations, an implicit treatment of neutrino-matter interactions is desired, which requires solutions of coupled nonlinear systems in each step of the time integration scheme. In this paper we design and compare nonlinear iterative solvers for implicit systems with energy coupling neutrino-matter interactions commonly used in CCSN simulations. Specifically, we consider electron neutrinos and antineutrinos, which interact with static matter configurations through the Bruenn~85 opacity set. The implicit systems arise from the discretization of a non-relativistic two-moment model for neutrino transport, which employs the discontinuous Galerkin (DG) method for phase-space discretization and an implicit-explicit (IMEX) time integration scheme. In the context of this DG-IMEX scheme, we propose two approaches to formulate the nonlinear systems -- a coupled approach and a nested approach. For each approach, the resulting systems are solved with Anderson-accelerated fixed-point iteration and Newton's method. The performance of these four iterative solvers has been compared on relaxation problems with various degree of collisionality, as well as proto-neutron star deleptonization problems with several matter profiles adopted from spherically symmetric CCSN simulations. Numerical results suggest that the nested Anderson-accelerated fixed-point solver is more efficient than other tested solvers for solving implicit nonlinear systems with energy coupling neutrino-matter interactions.

arXiv commentThis paper has been accepted for publication in ApJS

References in corpus (45)

• Explosion Mechanisms of Core-Collapse Supernovae
• Core-collapse supernova equations of state based on neutron star observations
• Perspectives on Core-Collapse Supernova Theory
• Radiation hydrodynamics with neutrinos: Variable Eddington factor method for core-collapse supernova simulations
• Nodal Discontinuous Galerkin Methods on Graphics Processors
• Shock Breakout in Core-Collapse Supernovae and its Neutrino Signature
• Neutrino Interactions in Hot and Dense Matter
• Supernova Simulations with Boltzmann Neutrino Transport: A Comparison of Methods
• A Finite Difference Representation of Neutrino Radiation Hydrodynamics in Spherically Symmetric General Relativistic Space-Time
• Three-dimensional core-collapse supernova simulated using a 15 $M_{\odot }$ progenitor
• The Overarching Framework of Core-Collapse Supernova Explosions as Revealed by 3D Fornax Simulations
• Neutrino Opacities in Nuclear Matter
• A New Multi-Dimensional General Relativistic Neutrino Hydrodynamics Code for Core-Collapse Supernovae II. Relativistic Explosion Models of Core-Collapse Supernovae
• Truncated Moment Formalism for Radiation Hydrodynamics in Numerical Relativity
• The consequences of nuclear electron capture in core collapse supernovae
• Supernova Neutrino Opacity from Nucleon-Nucleon Bremsstrahlung and Related Processes
• Three-Dimensional Simulations of Neutrino-Driven Core-Collapse Supernovae from Low-Mass Single and Binary Star Progenitors
• An Open-Source Neutrino Radiation Hydrodynamics Code for Core-Collapse Supernovae
• The Status of Multi-Dimensional Core-Collapse Supernova Models
• Neutrino-driven supernova of a low-mass iron-core progenitor boosted by three-dimensional turbulent convection
• Hydrodynamics of core-collapse supernovae and their progenitors
• Fornax: a Flexible Code for Multiphysics Astrophysical Simulations
• A new multidimensional, energy-dependent two-moment transport code for neutrino-hydrodynamics
• Global Comparison of Core-Collapse Supernova Simulations in Spherical Symmetry
• General Relativistic Three-Dimensional Multi-Group Neutrino Radiation-Hydrodynamics Simulations of Core-Collapse Supernovae
• Discontinuous Galerkin algorithms for fully kinetic plasmas
• Crucial Physical Dependencies of the Core-Collapse Supernova Mechanism
• Neutrino Transfer in Three Dimensions for Core-Collapse Supernovae. I. Static Configurations
• A New Multi-Energy Neutrino Radiation-Hydrodynamics Code in Full General Relativity and Its Application to Gravitational Collapse of Massive Stars
• Interplay of Neutrino Opacities in Core-collapse Supernova Simulations
• General Relativistic Effects in the Core Collapse Supernova Mechanism
• Core-collapse supernova simulations in one and two dimensions: comparison of codes and approximations
• A New Spherical Harmonics Scheme for Multi-Dimensional Radiation Transport I: Static Matter Configurations
• On the Requirements for Realistic Modeling of Neutrino Transport in Simulations of Core-Collapse Supernovae
• Conservative 3+1 General Relativistic Variable Eddington Tensor Radiation Transport Equations
• Mu and Tau Neutrino Thermalization and Production in Supernovae: Processes and Timescales
• Chimera: A massively parallel code for core-collapse supernova simulation
• Conservative 3+1 General Relativistic Boltzmann Equation
• 2D and 3D Core-Collapse Supernovae Simulation Results Obtained with the CHIMERA Code
• A Class of Multirate Infinitesimal GARK Methods
• Formulation of discontinuous Galerkin methods for relativistic astrophysics
• High Order Asymptotic Preserving Nodal Discontinuous Galerkin IMEX Schemes for the BGK Equation
• Realizability-Preserving DG-IMEX Method for the Two-Moment Model of Fermion Transport
• Bound-Preserving Discontinuous Galerkin Methods for Conservative Phase Space Advection in Curvilinear Coordinates
• A fast implicit solver for semiconductor models in one space dimension

Cited by in corpus (9)

• Flash-X, a multiphysics simulation software instrument
• 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
• 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
• A Finite Element Method for Angular Discretization of the Radiation Transport Equation on Spherical Geodesic Grids
• Global high-order numerical schemes for the time evolution of the general relativistic radiation magneto-hydrodynamics equations
• Time integration for neutrino radiation transport using minimally implicit Runge-Kutta methods
• A New Moment-Based General-Relativistic Neutrino-Radiation Transport Code: Methods and First Applications to Neutron Star Mergers