A Finite Element Method for Angular Discretization of the Radiation Transport Equation on Spherical Geodesic Grids
arXiv:2212.01409 · doi:10.1016/j.jcp.2023.112365
Abstract
Discrete ordinate () and filtered spherical harmonics () based schemes have been proven to be robust and accurate in solving the Boltzmann transport equation but they have their own strengths and weaknesses in different physical scenarios. We present a new method based on a finite element approach in angle that combines the strengths of both methods and mitigates their disadvantages. The angular variables are specified on a spherical geodesic grid with functions on the sphere being represented using a finite element basis. A positivity-preserving limiting strategy is employed to prevent non-physical values from appearing in the solutions. The resulting method is then compared with both and schemes using four test problems and is found to perform well when one of the other methods fail.
24 pages, 13 figures
References in corpus (7)
- Neutrino-driven winds from neutron star merger remnants
- Neutrino signatures and the neutrino-driven wind in Binary Neutron Star Mergers
- AT2017gfo: an anisotropic and three-component kilonova counterpart of GW170817
- Post-merger evolution of a neutron star-black hole binary with neutrino transport
- Impact of an improved neutrino energy estimate on outflows in neutron star merger simulations
- A New Spherical Harmonics Scheme for Multi-Dimensional Radiation Transport I: Static Matter Configurations
- A DG-IMEX method for two-moment neutrino transport: Nonlinear solvers for neutrino-matter coupling
Cited by in corpus (6)
- The Guided Moments formalism: a new efficient full-neutrino treatment for astrophysical simulations
- Performance-Portable Binary Neutron Star Mergers with AthenaK
- Neutrino Flavor Transformation in Neutron Star Mergers
- Impact of Neutrino Flavor Conversions on Neutron Star Merger Dynamics, Ejecta, Nucleosynthesis, and Multi-Messenger Signals
- Hierarchical Search of Tree Tensor Networks for High-Dimensional Data
- Non-Newtonian corrections to radiative viscosity: Israel-Stewart theory as a viscosity limiter