An angular multigrid preconditioner for the radiation transport equation with Fokker–Planck scattering

Research output: Contribution to journalArticleScientificpeer-review

4 Downloads (Pure)

Abstract

In a previous paper (Hennink and Lathouwers, 2017) we developed a finite element discretization for the Boltzmann transport equation with forward peaked scatter modeled by the Fokker–Planck approximation. The discretization was based on the discontinuous Galerkin method in both space and angle. It was expected and found that the regular source iteration algorithm for the Boltzmann equation is not effective in solving the discretized system and becomes excessively expensive for problems with many angular degrees of freedom. The purpose of this paper is to develop a multigrid scheme as preconditioner for the above mentioned discretization. The method exploits the nested nature of the meshes and the natural prolongation/restriction between meshes by Galerkin projection. A set of test problems ranging from pure spherical diffusion to the complete Boltzmann transport problem in 3D are presented to illustrate that the method is very effective, resulting in iteration counts nearly independent of problem size even for highly non-isotropically refined angular meshes.

Original languageEnglish
Pages (from-to)165-177
JournalJournal of Computational and Applied Mathematics
Volume350
DOIs
Publication statusPublished - 2019

Keywords

  • Discontinuous Galerkin
  • Fokker–Planck
  • Interior penalty
  • Multigrid
  • Particle transport
  • Radiation transport

Fingerprint Dive into the research topics of 'An angular multigrid preconditioner for the radiation transport equation with Fokker–Planck scattering'. Together they form a unique fingerprint.

  • Cite this