A realizability-preserving discontinuous Galerkin scheme for entropy-based moment closures for linear kinetic equations in one space dimension

From MaRDI portal
Publication:350100

DOI10.1016/J.JCP.2015.04.034zbMATH Open1349.82067arXiv1501.03660OpenAlexW2079085406MaRDI QIDQ350100FDOQ350100

Florian Schneider, Graham Alldredge

Publication date: 5 December 2016

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

Abstract: We implement a high-order numerical scheme for the entropy-based moment closure, the so-called MN model, for linear kinetic equations in slab geometry. A discontinuous Galerkin (DG) scheme in space along with a strong-stability preserving Runge-Kutta time integrator is a natural choice to achieve a third-order scheme, but so far, the challenge for such a scheme in this context is the implementation of a linear scaling limiter when the numerical solution leaves the set of realizable moments (that is, those moments associated with a positive underlying distribution). The difficulty for such a limiter lies in the computation of the intersection of a ray with the set of realizable moments. We avoid this computation by using quadrature to generate a convex polytope which approximates this set. The half-space representation of this polytope is used to compute an approximation of the required intersection straightforwardly, and with this limiter in hand, the rest of the DG scheme is constructed using standard techniques. We consider the resulting numerical scheme on a new manufactured solution and standard benchmark problems for both traditional MN models and the so-called mixed-moment models. The manufactured solution allows us to observe the expected convergence rates and explore the effects of the regularization in the optimization.


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





Cites Work


Cited In (20)

Uses Software






This page was built for publication: A realizability-preserving discontinuous Galerkin scheme for entropy-based moment closures for linear kinetic equations in one space dimension

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