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 nonrelativistic 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 degrees 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. * This manuscript has been authored in part by UT-Battelle, LLC, under contract DE-AC05-00OR22725 with the US Department of Energy (DOE). The US government retains and the publisher, by accepting the article for publication, acknowledges that the US government retains a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for US government purposes. The DOE will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan).
- Publication:
-
The Astrophysical Journal Supplement Series
- Pub Date:
- April 2021
- DOI:
- arXiv:
- arXiv:2102.02186
- Bibcode:
- 2021ApJS..253...52L
- Keywords:
-
- Computational methods;
- Core-collapse supernovae;
- Radiative transfer simulations;
- Supernova neutrinos;
- 1965;
- 304;
- 1967;
- 1666;
- Astrophysics - High Energy Astrophysical Phenomena
- E-Print:
- This paper has been accepted for publication in ApJS