Thornado-transport: Anderson- And GPU-accelerated nonlinear solvers for neutrino-matter coupling

Research output: Contribution to journalConference articlepeer-review

3 Scopus citations

Abstract

Algorithms for neutrino-matter coupling in core-collapse supernovae (CCSNe) are investigated in the context of a spectral two-moment model, which is discretized in space with the discontinuous Galerkin method, integrated in time with implicit-explicit (IMEX) methods, and implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado). The model considers electron neutrinos and antineutrinos and tabulated opacities from Bruenn (1985), which includes neutrino-electron scattering and pair processes. The nonlinear system arising from implicit time discretization of the equations governing neutrino-matter coupling is iterated to convergence using Anderson-accelerated fixed-point methods, which avoid formation of Jacobians and inversion of dense linear systems. Numerical experiments show that, for a given tolerance, a nested iteration scheme which aims to reduce opacity evaluations can lower the computational cost. Our initial port to GPUs, using both OpenMP and OpenACC, shows an overall speedup of up to ∼ 100 when compared to results using a single CPU core. These results indicate that the algorithms implemented in thornado are well-suited to GPU acceleration.

Original languageEnglish
Article number012013
JournalJournal of Physics: Conference Series
Volume1623
Issue number1
DOIs
StatePublished - Sep 24 2020
Event14th International Conference on Numerical Modeling of Space Plasma Flows, ASTRONUM 2019 - Paris, France
Duration: Jul 1 2019Jul 5 2019

Bibliographical note

Publisher Copyright:
© Published under licence by IOP Publishing Ltd.

Fingerprint

Dive into the research topics of 'Thornado-transport: Anderson- And GPU-accelerated nonlinear solvers for neutrino-matter coupling'. Together they form a unique fingerprint.

Cite this