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

MP Laiu, E Endeve, R Chu, JA Harris… - The Astrophysical …, 2021 - iopscience.iop.org
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 …

A collision-based hybrid method for the BGK equation

M Shin, CD Hauck, RG McClarren - Journal of Computational Physics, 2024 - Elsevier
We apply the collision-based hybrid method introduced by Hauck and McClarren [1] to the
Boltzmann equation with the BGK operator and a hyperbolic scaling. An implicit treatment of …

Anderson acceleration with approximate calculations: applications to scientific computing

M Lupo Pasini, MP Laiu - Numerical Linear Algebra with …, 2024 - Wiley Online Library
We provide rigorous theoretical bounds for Anderson acceleration (AA) that allow for
approximate calculations when applied to solve linear problems. We show that, when the …

An Asymptotic Preserving Discontinuous Galerkin Method for a Linear Boltzmann Semiconductor Model

VP DeCaria, CD Hauck, SR Schnake - SIAM Journal on Numerical Analysis, 2024 - SIAM
A key property of the linear Boltzmann semiconductor model is that as the collision
frequency tends to infinity, the phase space density converges to an isotropic function, called …

Analysis of a new implicit solver for a semiconductor model

VP DeCaria, CD Hauck, MTP Laiu - SIAM Journal on Scientific Computing, 2021 - SIAM
We present and analyze a new iterative solver for implicit discretizations of a simplified
Boltzmann--Poisson system. The algorithm builds on recent work that incorporated a …

High Order IMEX Stochastic Galerkin Schemes for Linear Transport Equation with Random Inputs and Diffusive Scalings

Z Chen, L Mu - Communications on Applied Mathematics and …, 2024 - Springer
In this paper, we consider the high order method for solving the linear transport equations
under diffusive scaling and with random inputs. To tackle the randomness in the problem …