Bound-Preserving Discontinuous Galerkin Methods for Conservative Phase Space Advection in Curvilinear Coordinates

Abstract

We extend the positivity-preserving method of Zhang & Shu (2010, JCP, 229, 3091-3120) to simulate the advection of neutral particles in phase space using curvilinear coordinates. The ability to utilize these coordinates is important for non-equilibrium transport problems in general relativity and also in science and engineering applications with specific geometries. The method achieves high-order accuracy using Discontinuous Galerkin (DG) discretization of phase space and strong stability-preserving, Runge-Kutta (SSP-RK) time integration. Special care in taken to ensure that the method preserves strict bounds for the phase space distribution function $f$; i.e., $f\in [0,1]$. The combination of suitable CFL conditions and the use of the high-order limiter proposed in Zhang & Shu (2010) is sufficient to ensure positivity of the distribution function. However, to ensure that the distribution function satisfies the upper bound, the discretization must, in addition, preserve the divergence-free property of the phase space flow. Proofs that highlight the necessary conditions are presented for general curvilinear coordinates, and the details of these conditions are worked out for some commonly used coordinate systems (i.e., spherical polar spatial coordinates in spherical symmetry and cylindrical spatial coordinates in axial symmetry, both with spherical momentum coordinates). Results from numerical experiments --- including one example in spherical symmetry adopting the Schwarzschild metric --- demonstrate that the method achieves high-order accuracy and that the distribution function satisfies the maximum principle.

Submitted to JCP

References in corpus (14)

• Explosion Mechanisms of Core-Collapse Supernovae
• Perspectives on Core-Collapse Supernova Theory
• Nodal Discontinuous Galerkin Methods on Graphics Processors
• Explosion Mechanism, Neutrino Burst, and Gravitational Wave in Core-Collapse Supernovae
• A Finite Difference Representation of Neutrino Radiation Hydrodynamics in Spherically Symmetric General Relativistic Space-Time
• 3+1 Formalism and Bases of Numerical Relativity
• A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov-Poisson equations
• Interplay of Neutrino Opacities in Core-collapse Supernova Simulations
• Three-dimensional Boltzmann-Hydro code for core-collapse in massive stars I. special relativistic treatments
• A New Spherical Harmonics Scheme for Multi-Dimensional Radiation Transport I: Static Matter Configurations
• Third order Maximum-Principle-Satisfying Direct discontinuous Galerkin methods for time dependent convection diffusion equations on unstructured triangle mesh
• Conservative 3+1 General Relativistic Boltzmann Equation
• Hyperbolic character of the angular moment equations of radiative transfer and numerical methods
• Conservative formulations of general relativistic kinetic theory

Cited by in corpus (10)

• SpECTRE: A Task-based Discontinuous Galerkin Code for Relativistic Astrophysics
• Design of Provably Physical-Constraint-Preserving Methods for General Relativistic Hydrodynamics
• Formulation of discontinuous Galerkin methods for relativistic astrophysics
• Realizability-Preserving DG-IMEX Method for the Two-Moment Model of Fermion Transport
• A DG-IMEX method for two-moment neutrino transport: Nonlinear solvers for neutrino-matter coupling
• Simulating magnetized neutron stars with discontinuous Galerkin methods
• thornado-hydro: a discontinuous Galerkin method for supernova hydrodynamics with nuclear equations of state
• A hp-adaptive discontinuous Galerkin solver for elliptic equations in numerical relativity
• Positivity preserving DG schemes for a Boltzmann - Poisson model of electrons in semiconductors in curvilinear momentum coordinates
• Entropy-stable positivity-preserving DG schemes for Boltzmann-Poisson models of collisional electronic transport along energy bands