Realizability-preserving DG-IMEX method for the two-moment model of fermion transport

DOI10.1016/J.JCP.2019.03.037MaRDI QIDQ2218196zbMATH OpenOpenAlexDBLPWikidataFDO

Authors Ran Chu, Eirik Endeve, Cory D. Hauck, Anthony Mezzacappa

Publication date 15 January 2021

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1809.06949

zbMATH Keywords

Boltzmann equationdiscontinuous Galerkinhyperbolic conservation lawsimplicit-explicitmoment realizabilityradiation transport

Mathematics Subject Classification ID

Hyperbolic conservation laws (35L65) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Applications to the sciences (65Z05) Fermionic systems in quantum theory (81V74) Computational methods for problems pertaining to astronomy and astrophysics (85-08) Radiative transfer in astronomy and astrophysics (85A25)

Abstract: Building on the framework of Zhang & Shu cite{zhangShu_2010a,zhangShu_2010b}, we develop a realizability-preserving method to simulate the transport of particles (fermions) through a background material using a two-moment model that evolves the angular moments of a phase space distribution function $f$. The two-moment model is closed using algebraic moment closures; e.g., as proposed by Cernohorsky & Bludman cite{cernohorskyBludman_1994} and Banach & Larecki cite{banachLarecki_2017a}. Variations of this model have recently been used to simulate neutrino transport in nuclear astrophysics applications, including core-collapse supernovae and compact binary mergers. We employ the discontinuous Galerkin (DG) method for spatial discretization (in part to capture the asymptotic diffusion limit of the model) combined with implicit-explicit (IMEX) time integration to stably bypass short timescales induced by frequent interactions between particles and the background. Appropriate care is taken to ensure the method preserves strict algebraic bounds on the evolved moments (particle density and flux) as dictated by Pauli's exclusion principle, which demands a bounded distribution function (i.e., $fin[0,1]$). This realizability-preserving scheme combines a suitable CFL condition, a realizability-enforcing limiter, a closure procedure based on Fermi-Dirac statistics, and an IMEX scheme whose stages can be written as a convex combination of forward Euler steps combined with a backward Euler step. Numerical results demonstrate the realizability-preserving properties of the scheme. We also demonstrate that the use of algebraic moment closures not based on Fermi-Dirac statistics can lead to unphysical moments in the context of fermion transport.

Recommendations

• Kershaw closures for linear transport equations in slab geometry. II: High-order realizability-preserving discontinuous-Galerkin schemes
• A realizability-preserving discontinuous Galerkin method for the \(M_1\) model of radiative transfer
• A realizability-preserving discontinuous Galerkin scheme for entropy-based moment closures for linear kinetic equations in one space dimension
• A comparative study of limiting strategies in discontinuous Galerkin schemes for the M₁ model of radiation transport
• Kershaw-type transport equations for fermionic radiation

Cites work

• A comparison of moment closures for linear kinetic transport equations: the line source benchmark
• A new spherical harmonics scheme for multi-dimensional radiation transport. I: Static matter configurations
• A realizability-preserving discontinuous Galerkin method for the \(M_1\) model of radiative transfer
• Asymptotic analysis of upwind discontinuous Galerkin approximation of the radiative transport equation in the diffusive limit
• Asymptotic preserving implicit-explicit Runge--Kutta methods for nonlinear kinetic equations
• Asymptotic solutions of numerical transport problems in optically thick, diffusive regimes. II
• Asymptotic-preserving and positivity-preserving implicit-explicit schemes for the stiff BGK equation
• Bound-preserving discontinuous Galerkin methods for conservative phase space advection in curvilinear coordinates
• Convex Duality and Entropy-Based Moment Closures: Characterizing Degenerate Densities
• Domain of definition of Levermore's five-moment system
• Efficient implementation of essentially nonoscillatory shock-capturing schemes
• Entropy-based mixed three-moment description of fermionic radiation transport in slab and spherical geometries
• Formulation of discontinuous Galerkin methods for relativistic astrophysics
• High order WENO and DG methods for time-dependent convection-dominated PDEs: A brief survey of several recent developments
• High-order accurate physical-constraints-preserving finite difference WENO schemes for special relativistic hydrodynamics
• scientific article; zbMATH DE number 914383 (Why is no real title available?)
• scientific article; zbMATH DE number 5179816 (Why is no real title available?)
• scientific article; zbMATH DE number 3107192 (Why is no real title available?)
• Implicit filtered \(P_N\) for high-energy density thermal radiation transport using discontinuous Galerkin finite elements
• Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations
• Implicit-explicit Runge-Kutta schemes and applications to hyperbolic systems with relaxation
• Kershaw-type transport equations for fermionic radiation
• Low-storage implicit/explicit Runge-Kutta schemes for the simulation of stiff high-dimensional ODE systems
• Maximum-principle-satisfying second order discontinuous Galerkin schemes for convection-diffusion equations on triangular meshes
• Moment closure hierarchies for kinetic theories.
• Nodal discontinuous Galerkin methods on graphics processors
• Nonoscillatory High Order Accurate Self-Similar Maximum Principle Satisfying Shock Capturing Schemes I
• Numerical schemes for hyperbolic conservation laws with stiff relaxation terms
• On maximum-principle-satisfying high order schemes for scalar conservation laws
• On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes
• Perturbed, entropy-based closure for radiative transfer
• Positivity-preserving DG and central DG methods for ideal MHD equations
• Positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations with source terms
• Relativistic radiative transfer: moment formalisms
• Relativistic transport theory
• Robust and accurate filtered spherical harmonics expansions for radiative transfer
• Runge--Kutta discontinuous Galerkin methods for convection-dominated problems
• Semi-implicit time integration for \(P_N\) thermal radiative transfer
• Spectral maximum entropy hydrodynamics of fermionic radiation: a three-moment system for one-dimensional flows
• Steady state and sign preserving semi-implicit Runge-Kutta methods for ODEs with stiff damping term
• Strong stability-preserving high-order time discretization methods
• Truncated Moment Formalism for Radiation Hydrodynamics in Numerical Relativity
• Two-dimensional time dependent Riemann solvers for neutron transport

Cited in

(4)

• Conservative DG method for the micro-macro decomposition of the Vlasov-Poisson-Lenard-Bernstein model
• DG-IMEX method for a two-moment model for radiation transport in the \(\mathcal{O}(v/c)\) limit
• Sparse-grid discontinuous Galerkin methods for the Vlasov-Poisson-Lenard-Bernstein model
• Adaptive Sparse-grid Discontinuous Galerkin Approximations the Bhatnagar--Gross--Krook Model

This page was built for publication: Realizability-preserving DG-IMEX method for the two-moment model of fermion transport

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2218196)