A discontinuous Galerkin method for the mono-energetic Fokker-Planck equation based on a spherical interior penalty formulation

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:330

Article

Hennink, Aldo¹  Lathouwers, Danny¹ 

[1] Delft Univ Technol, Dept Radiat Sci & Technol, Mekelweg 15, NL-2629 JB Delft, Netherlands

关键词: Discontinuous Galerkin;    Fokker-Planck;    Particle transport;    Radiation transport;    Upwinding;    Interior penalty;   

DOI  :  10.1016/j.cam.2017.08.019

来源: Elsevier

【 摘 要 】

We present a new discretization of the mono-energetic Fokker-Planck equation. We build on previous work (Kophazi and Lathouwers, 2015) where we devised an angular discretization for the Boltzmann equation, allowing for both heterogeneous and anisotropic angular refinement. The angular discretization is based on a discontinuous finite element method on the unit sphere. Here we extend the methodology to include the effect of the Fokker-Planck scatter operator describing small angle particle scatter. We describe the construction of an interior penalty method on the sphere surface. Results are provided for a variety of test cases, ranging from purely angular to fully three-dimensional. The results show that the scheme can resolve highly forward-peaked flux distributions with forward peaked scatter. (C) 2017 Elsevier B.V. All rights reserved.

【 授权许可】

Free