An adaptive nested source term iteration for radiative transfer equations

DOI10.1090/MCOM/3505zbMATH Open1437.65183arXiv1810.07035OpenAlexW2990455384WikidataQ126651684 ScholiaQ126651684MaRDI QIDQ4960070FDOQ4960070

O. Mula, F. Gruber, Wolfgang Dahmen

Publication date: 8 April 2020

Published in: Mathematics of Computation (Search for Journal in Brave)

Abstract: 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.

Full work available at URL: https://arxiv.org/abs/1810.07035

Recommendations

• Sparse adaptive finite elements for radiative transfer
• On robustly convergent and efficient iterative methods for anisotropic radiative transfer
• On a convergent DSA preconditioned source iteration for a DGFEM method for radiative transfer
• Reliable anisotropic-adaptive discontinuous Galerkin method for simplified \(\mathbf{P}_{\mathbf{N}}\) approximations of radiative transfer
• A posteriori error analysis in radiative transfer

zbMATH Keywords

radiative transfermatrix compressionkinetic problemslinear Boltzmanna posteriori boundsHilbert-Schmidt decompositionDPG transport solverfast application of scattering operatoriteration in function space

Mathematics Subject Classification ID

Scattering theory for PDEs (35P25) Boltzmann equations (35Q20) Stability and convergence of numerical methods for boundary value problems involving PDEs (65N12) Error bounds for boundary value problems involving PDEs (65N15) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Mesh generation, refinement, and adaptive methods for boundary value problems involving PDEs (65N50)

Cites Work

• Inverse transport theory and applications
• Title not available (Why is that?)
• Title not available (Why is that?)
• AnL^ptheory for stationary radiative transfer
• Adaptive wavelet methods. II: Beyond the elliptic case
• A mixed variational framework for the radiative transfer equation
• Adaptive wavelet methods for elliptic operator equations: Convergence rates
• Compression Techniques for Boundary Integral Equations---Asymptotically Optimal Complexity Estimates
• Adaptive Petrov--Galerkin Methods for First Order Transport Equations
• Sparse discrete ordinates method in radiative transfer
• Asymptotic Analysis of Upwind Discontinuous Galerkin Approximation of the Radiative Transport Equation in the Diffusive Limit
• Convergence of a Fully Discrete Scheme for Two-Dimensional Neutron Transport
• A robust \(S_{N}\)-DG-approximation for radiation transport in optically thick and diffusive regimes
• A Class of Bases in $L^2$ for the Sparse Representation of Integral Operators
• Solution of Radiative Transfer Problems with Finite Elements
• Spatial approximation of the radiation transport equation using a subgrid-scale finite element method
• L₂-error estimates for the discrete ordinates method for three-dimensional neutron transport
• A Finite Element Method for the Neutron Transport Equation in an Infinite Cylindrical Domain
• Adaptive strategies for transport equations
• On the stability of DPG formulations of transport equations

Cited In (8)

• On robustly convergent and efficient iterative methods for anisotropic radiative transfer
• Localization of discontinuity surfaces of the scattering coefficient according to the time-angular distribution of the radiation flux density
• A phase-space discontinuous Galerkin scheme for the radiative transfer equation in slab geometry
• Efficient approximation of solutions of parametric linear transport equations by ReLU DNNs
• On the problem of determining the scattering coefficient in frequency modulated sounding of a medium
• Sparse grid time-discontinuous Galerkin method with streamline diffusion for transport equations
• Fast, low-memory numerical methods for radiative transfer via \(hp\)-adaptive mesh refinement
• Comparative analysis of the error of the single scattering approximation when solving one inverse problem in two-dimensional and three- dimensional cases

Uses Software

• DUNE

This page was built for publication: An adaptive nested source term iteration for radiative transfer equations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4960070)