Abstract:Linear kinetic transport equations play a critical role in optical tomography, radiative transfer and neutron transport. The fundamental difficulty hampering their efficient and accurate numerical resolution lies in the high dimensionality of the physical and velocity/angular variables and the fact that the problem is multiscale in nature. Leveraging the existence of a hidden low-rank structure hinted by the diffusive limit, in this work, we design and test the angular-space reduced order model for the linear radiative transfer equation, the first such effort based on the celebrated reduced basis method (RBM). Our method is built upon a high-fidelity solver employing the discrete ordinates method in the angular space, an asymptotic preserving upwind discontinuous Galerkin method for the physical space, and an efficient synthetic accelerated source iteration for the resulting linear system. Addressing the challenge of the parameter values (or angular directions) being coupled through an integration operator, the first novel ingredient of our method is an iterative procedure where the macroscopic density is constructed from the RBM snapshots, treated explicitly and allowing a transport sweep, and then updated afterwards. A greedy algorithm can then proceed to adaptively select the representative samples in the angular space and form a surrogate solution space. The second novelty is a least squares density reconstruction strategy, at each of the relevant physical locations, enabling the robust and accurate integration over an arbitrarily unstructured set of angular samples toward the macroscopic density. Numerical experiments indicate that our method is effective for computational cost reduction in a variety of regimes.

Low-rank tensor product Richardson iteration for radiative transfer in plane-parallel geometry

Iterative space-angle discontinuous Galerkin method for radiative transfer equation

Asymptotic-preserving and energy stable dynamical low-rank approximation for thermal radiative transfer equations

A Reduced Basis Method for Radiative Transfer Equation

A class of Galerkin schemes for time-dependent radiative transfer

Preconditioned low-rank Riemannian optimization for linear systems with tensor product structure

An Adaptive Angular Domain Compression Scheme for Solving Multiscale Radiative Transfer Equation

An Adaptive Angular Aomain Compression Scheme For Solving Multiscale Radiative Transfer Equation

A Nonlinear Projection-Based Iteration Scheme with Cycles over Multiple Time Steps for Solving Thermal Radiative Transfer Problems

Discrete-Ordinate Discontinuous Galerkin Methods for Solving the Radiative Transfer Equation

A Theoretical Study for RTE-based Parameter Identification Problems

Exact Analytical Solution of the One-Dimensional Time-Dependent Radiative Transfer Equation with Linear Scattering

Simulation of the 3D Radiative Transfer with Anisotropic Scattering for Convective Trails

A low-rank method for two-dimensional time-dependent radiation transport calculations

On accelerated iterative schemes for anisotropic radiative transfer using residual minimization

A Uniformly Convergent Scheme for Radiative Transfer Equation in the Diffusion Limit Up to the Boundary and Interface Layers

An Effective Approach to Solve Radiative Transfer Equations in Strongly Forward Scattering Environments and Its Applications in Vegetation Scattering

Reflective Conditions for Radiative Transfer in Integral Form with H-Matrices

Multiscale Solutions of Radiative Heat Transfer by the Discrete Unified Gas Kinetic Scheme

Three-dimensional Radiative Transfer with Multilevel Atoms

Radiative Heat Transfer and Effective Transport Coefficients