An Adaptive Nested Source Term Iteration for Radiative Transfer Equations
Dahmen, Wolfgang; Gruber, Felix; Mula, Olga (2020), An Adaptive Nested Source Term Iteration for Radiative Transfer Equations, Mathematics of Computation, 89, p. 1605-1646. 10.1090/mcom/3505
Type
Article accepté pour publication ou publiéExternal document link
https://hal.archives-ouvertes.fr/hal-01899058Date
2020Journal name
Mathematics of ComputationVolume
89Publisher
American Mathematical Society
Pages
1605-1646
Publication identifier
Metadata
Show full item recordAuthor(s)
Dahmen, WolfgangUniversity of South Carolina [Columbia]
Gruber, Felix
Rheinisch-Westfälische Technische Hochschule Aachen University [RWTH]
Mula, Olga

CEntre de REcherches en MAthématiques de la DEcision [CEREMADE]
Abstract (EN)
We propose a new approach to the numerical solution of radiative transfer equations with certified a posteriori error bounds. A key role is played by stable Petrov--Galerkin type variational formulations of parametric transport equations and corresponding radiative transfer equations. This allows us to formulate an iteration in a suitable, infinite dimensional function space that is guaranteed to converge with a fixed error reduction per step. The numerical scheme is then based on approximately realizing this iteration within dynamically updated accuracy tolerances that still ensure convergence to the exact solution. To advance this iteration two operations need to be performed within suitably tightened accuracy tolerances. First, the global scattering operator needs to be approximately applied to the current iterate within a tolerance comparable to the current accuracy level. Second, parameter dependent linear transport equations need to be solved, again at the required accuracy of the iteration. To ensure that the stage dependent error tolerances are met, one has to employ rigorous a posteriori error bounds which, in our case, rest on a Discontinuous Petrov--Galerkin (DPG) scheme. These a posteriori bounds are not only crucial for guaranteeing the convergence of the perturbed iteration but are also used to generate adapted parameter dependent spatial meshes. This turns out to significantly reduce overall computational complexity. Since the global operator is only applied, we avoid the need to solve linear systems with densely populated matrices. Moreover, the approximate application of the global scatterer accelerated through low-rank approximation and matrix compression techniques. The theoretical findings are illustrated and complemented by numerical experiments with non-trivial scattering kernels.Subjects / Keywords
A posteriori bounds; Matrix compression; Hilbert–Schmidt decomposition; Fast application of scattering operator; Iteration in function space; DPG transport solver; Kinetic problems; Linear Boltzmann; Radiative transferRelated items
Showing items related by title and author.
-
Cohen, Albert; Dahmen, Wolfgang; DeVore, Ron; Fadili, Jalal M.; Mula, Olga; Nichols, James (2020) Article accepté pour publication ou publié
-
Cohen, Albert; Dahmen, Wolfgang; Mula, Olga; Nichols, James (2022) Article accepté pour publication ou publié
-
Gruber, Felix; Klewinghaus, Angela; Mula, Olga (2017) Article accepté pour publication ou publié
-
Maday, Yvon; Mula, Olga (2020) Article accepté pour publication ou publié
-
Aghili, Joubine; Mula, Olga (2020) Document de travail / Working paper