On a convergent DSA preconditioned source iteration for a DGFEM method for radiative transfer

Olena Palii*, Matthias Schlottbom*

*Corresponding author for this work

Research output: Contribution to journalArticleAcademicpeer-review

4 Citations (Scopus)
99 Downloads (Pure)


We consider the numerical approximation of the radiative transfer equation using discontinuous angular and continuous spatial approximations for the even parts of the solution. The even-parity equations are solved using a diffusion synthetic accelerated source iteration. We provide a convergence analysis for the infinite-dimensional iteration as well as for its discretized counterpart. The diffusion correction is computed by a subspace correction, which leads to a convergence behavior that is robust with respect to the discretization. The proven theoretical contraction rate deteriorates for scattering dominated problems. We show numerically that the preconditioned iteration is in practice robust in the diffusion limit. Moreover, computations for the lattice problem indicate that the presented discretization does not suffer from the ray effect. The theoretical methodology is presented for plane-parallel geometries with isotropic scattering, but the approach and proofs generalize to multi-dimensional problems and more general scattering operators verbatim.

Original languageEnglish
Pages (from-to)3366-3377
Number of pages12
JournalComputers & mathematics with applications
Issue number12
Early online date22 Feb 2020
Publication statusPublished - 15 Jun 2020


  • Convergence rates
  • Diffusion synthetic acceleration
  • Discontinuous angular approximation
  • Discrete ordinates method
  • Radiative transfer


Dive into the research topics of 'On a convergent DSA preconditioned source iteration for a DGFEM method for radiative transfer'. Together they form a unique fingerprint.

Cite this