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.

A sweep-based low-rank method for the discrete ordinate transport equation

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

A low-rank method for time-dependent transport calculations

A robust collision source method for rank adaptive dynamical low-rank approximation in radiation therapy

Asymptotic-Preserving and Energy Stable Dynamical Low-Rank Approximation

A Reduced Basis Method for Radiative Transfer Equation

A Deterministic Dynamical Low-rank Approach for Charged Particle Transport

A semi-implicit dynamical low-rank discontinuous Galerkin method for space homogeneous kinetic equations. Part I: emission and absorption

Energy stable and conservative dynamical low-rank approximation for the Su-Olson problem

Interpolatory dynamical low-rank approximation for the 3+3d Boltzmann-BGK equation

A stable multiplicative dynamical low-rank discretization for the linear Boltzmann-BGK equation

A Sparse Grid Discrete Ordinate Discontinuous Galerkin Method for the Radiative Transfer Equation

High order conservative Lagrangian schemes for two-dimensional radiation hydrodynamics equations in the equilibrium-diffusion limit

An Implicit Finite Volume Scheme to Solve the Time Dependent Radiation Transport Equation Based on Discrete Ordinates

Lagrange Discrete Ordinates: a new angular discretization for the three dimensional linear Boltzmann equation

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

Least-squares Finite-Element SN Method for Solving Three-Dimensional Transport Equation

A local character based method for solving linear systems of radiation diffusion problems

A Discontinuous Galerkin Finite Element Method Based Axial SN for the 2D/1D Transport Method

The Extended Diamond Difference - Constant Nodal Method with Decoupled Cell Iteration Scheme in Two-Dimensional Discrete Ordinate Transport Problems

A review of low-rank methods for time-dependent kinetic simulations