On robustly convergent and efficient iterative methods for anisotropic radiative transfer

DOI10.1007/S10915-021-01757-9zbMATH Open1487.65036arXiv2102.09038OpenAlexW3129746008MaRDI QIDQ2113633FDOQ2113633

Jürgen Dölz, Olena Palii, Matthias Schlottbom

Publication date: 14 March 2022

Published in: Journal of Scientific Computing (Search for Journal in Brave)

Abstract: This paper considers the iterative solution of linear systems arising from discretization of the anisotropic radiative transfer equation with discontinuous elements on the sphere. In order to achieve robust convergence behavior in the discretization parameters and the physical parameters we develop preconditioned Richardson iterations in Hilbert spaces. We prove convergence of the resulting scheme. The preconditioner is constructed in two steps. The first step borrows ideas from matrix splittings and ensures mesh independence. The second step uses a subspace correction technique to reduce the influence of the optical parameters. The correction spaces are build from low-order spherical harmonics approximations generalizing well-known diffusion approximations. We discuss in detail the efficient implementation and application of the discrete operators. In particular, for the considered discontinuous spherical elements, the scattering operator becomes dense and we show that $mathcalH$- or $mathcal{H}^2$-matrix compression can be applied in a black-box fashion to obtain almost linear or linear complexity when applying the corresponding approximations. The effectiveness of the proposed method is shown in numerical examples.

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

zbMATH Keywords

preconditioningcompressioniterative solutionanisotropic radiative transfer

Mathematics Subject Classification ID

Preconditioners for iterative methods (65F08) Iterative numerical methods for linear systems (65F10) Numerical solution of discretized equations for boundary value problems involving PDEs (65N22) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Method of contraction of the boundary for boundary value problems involving PDEs (65N45)

Cites Work

• Title not available (Why is that?)
• Title not available (Why is that?)
• Title not available (Why is that?)
• Efficient numerical methods for non-local operators. \(\mathcal H^2\)-matrix compression, algorithms and analysis.
• A fast algorithm for particle simulations
• Optical tomography: forward and inverse problems
• The method of alternating projections and the method of subspace corrections in Hilbert space
• A mixed variational framework for the radiative transfer equation
• Hierarchical Matrices: Algorithms and Analysis
• Discrete-Ordinate Discontinuous Galerkin Methods for Solving the Radiative Transfer Equation
• Compression Techniques for Boundary Integral Equations---Asymptotically Optimal Complexity Estimates
• Preconditioning of complex symmetric linear systems with applications in optical tomography
• Multi-group full-spectrum \(k\)-distribution database for water vapor mixtures in radiative transfer calculations.
• A Robust Multigrid Preconditioner for $S_N$DG Approximation of Monochromatic, Isotropic Radiation Transport Problems
• The black-box fast multipole method
• A Hybrid Finite-Volume, Discontinuous Galerkin Discretization for the Radiative Transport Equation
• A Discrete-Ordinate Discontinuous-Streamline Diffusion Method for the Radiative Transfer Equation
• Discontinuous Galerkin for the Radiative Transport Equation
• Comparison of fast boundary element methods on parametric surfaces
• A two-mesh adaptive mesh refinement technique for \(S_N\) neutral-particle transport using a higher-order DGFEM
• A Boundary Functional for the Least-Squares Finite- Element Solution of Neutron Transport Problems
• Stationary radiative transfer with vanishing absorption
• Preconditioning
• A space-angle DGFEM approach for the Boltzmann radiation transport equation with local angular refinement
• A Perfectly Matched Layer Approach for $P_N$-Approximations in Radiative Transfer
• Sweep-Net: an artificial neural network for radiation transport solves
• A Multigrid Framework for $S_n$ Discretizations of the Boltzmann Transport Equation
• Preconditioned Krylov subspace methods for solving radiative transfer problems with scattering and reflection
• An adaptive nested source term iteration for radiative transfer equations
• Diffusion Synthetic Acceleration Preconditioning for Discontinuous Galerkin Discretizations of $S_N$ Transport on High-Order Curved Meshes
• Improvements to a class of hybrid methods for radiation transport: Nyström reconstruction and defect correction methods
• On a convergent DSA preconditioned source iteration for a DGFEM method for radiative transfer
• A cascadic multigrid asymptotic-preserving discrete ordinate discontinuous streamline diffusion method for radiative transfer equations with diffusive scalings
• Certified reduced-order methods for optimal treatment planning
• Low-Memory, Discrete Ordinates, Discontinuous Galerkin Methods for Radiative Transport

Cited In (2)

• Fast, low-memory numerical methods for radiative transfer via \(hp\)-adaptive mesh refinement
• Numerical radiative transfer with state-of-the-art iterative methods made easy

Uses Software

• hlib
• Bembel
• H2Lib

This page was built for publication: On robustly convergent and efficient iterative methods for anisotropic radiative transfer

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