A comparative study of limiting strategies in discontinuous Galerkin schemes for the M1 model of radiation transport

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:342

Article

Chidyagwai, Prince¹  Frank, Martin²  Schneider, Florian³  Seibold, Benjamin⁴ 

[1] Loyola Univ Maryland, Dept Math & Stat, 4501 N Charles St, Baltimore, MD 21210 USA

[2] Rhein Westfal TH Aachen, Dept Math, Schinkelstr 2, D-52062 Aachen, Germany

[3] TU Kaiserslautern, Fachbereich Math, Erwin Schrodinger Str, D-67663 Kaiserslautern, Germany

[4] Temple Univ, Dept Math, 1805 N Broad St, Philadelphia, PA 19122 USA

关键词: Discontinuous Galerkin;    Moment models;    Minimum entropy;    Realizability limiter;   

DOI  :  10.1016/j.cam.2018.04.017

来源: Elsevier

【 摘 要 】

The M-1 minimum entropy moment system is a system of hyperbolic balance laws that approximates the radiation transport equation, and has many desirable properties. Among them are symmetric hyperbolicity, entropy decay, moment realizability, and correct behavior in the diffusion and free-streaming limits. However, numerical difficulties arise when approximating the solution of the M1 model by high order numerical schemes; namely maintaining the realizability of the numerical solution and controlling spurious oscillations. In this paper, we extend a previously constructed one-dimensional realizability limiting strategy to 2D. In addition, we perform a numerical study of various combinations of the realizability limiter and the TVBM local slope limiter on a third order Discontinuous Galerkin (DG) scheme on both triangular and rectangular meshes. In several test cases, we demonstrate that in general, a combination of the realizability limiter and a TVBM limiter is necessary to obtain a robust and accurate numerical scheme. Our code is published so that all results can be reproduced by the reader. (C) 2018 Elsevier B.V. All rights reserved.

【 授权许可】

Free