A Finite Element Method for Angular Discretization of the Radiation Transport Equation on Spherical Geodesic Grids

Abstract

Discrete ordinate ($S_N$) and filtered spherical harmonics ($F{P}_N$) 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 $S_N$ and $F{P}_N$ schemes using four test problems and is found to perform well when one of the other methods fail.

arXiv comment24 pages, 13 figures

References in corpus (29)

• Neutrino-driven winds from neutron star merger remnants
• Dynamical Mass Ejection from Binary Neutron Star Mergers
• Neutrino signatures and the neutrino-driven wind in Binary Neutron Star Mergers
• Dynamical mass ejection from the merger of asymmetric binary neutron stars: Radiation-hydrodynamics study in general relativity
• Gravitational waves and neutrino emission from the merger of binary neutron stars
• AT2017gfo: an anisotropic and three-component kilonova counterpart of GW170817
• An Open-Source Neutrino Radiation Hydrodynamics Code for Core-Collapse Supernovae
• Low mass binary neutron star mergers : gravitational waves and neutrino emission
• 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
• Fornax: a Flexible Code for Multiphysics Astrophysical Simulations
• General Relativistic Three-Dimensional Multi-Group Neutrino Radiation-Hydrodynamics Simulations of Core-Collapse Supernovae
• Properties of Neutrino-driven Ejecta from the Remnant of Binary Neutron Star Merger : Purely Radiation Hydrodynamics Case
• 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
• Two Dimensional Core-Collapse Supernova Explosions Aided by General Relativity with Multidimensional Neutrino Transport
• Two-dimensional, Time-dependent, Multi-group, Multi-angle Radiation Hydrodynamics Test Simulation in the Core-Collapse Supernova Context
• A New Monte Carlo Method for Time-Dependent Neutrino Radiation Transport
• Three-Dimensional Core-Collapse Supernova Simulations with Multi-Dimensional Neutrino Transport Compared to the Ray-by-Ray-plus Approximation
• A New Spherical Harmonics Scheme for Multi-Dimensional Radiation Transport I: Static Matter Configurations
• Monte Carlo Neutrino Transport Through Remnant Disks from Neutron Star Mergers
• Neutrino pair annihilation above merger remnants: implications of a long-lived massive neutron star
• A Novel Multi-Dimensional Boltzmann Neutrino Transport Scheme for Core-Collapse Supernovae
• Rank-3 Moment Closures in General Relativistic Neutrino Transport
• NADA-FLD: A General Relativistic, Multi-dimensional Neutrino-hydrodynamics Code Employing Flux-limited Diffusion
• General relativistic neutrino transport using spectral methods
• A DG-IMEX method for two-moment neutrino transport: Nonlinear solvers for neutrino-matter coupling
• A New Discrete Implicit Monte Carlo Scheme for Simulating Radiative Transfer Problems
• A Parallel Adaptive Event Horizon Finder for Numerical Relativity

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