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.

Asymptotic-Preserving and Energy Stable Dynamical Low-Rank Approximation

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

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

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

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

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

A new approach to almost-sure asymptotic stability of stochastic systems of higher dimension

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

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

Dynamical Low-Rank Approximation for Stochastic Differential Equations

Randomised methods for dynamical low-rank approximation

A Deterministic Dynamical Low-rank Approach for Charged Particle Transport

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

Projected exponential methods for stiff dynamical low-rank approximation problems

Robust and conservative dynamical low-rank methods for the Vlasov equation via a novel macro-micro decomposition

A Reduced Basis Method for Radiative Transfer Equation

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

Second-order robust parallel integrators for dynamical low-rank approximation

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

Low rank approximation for the numerical simulation of high dimensional Lindblad and Riccati equations

An efficient scheme for approximating long-time dynamics of a class of non-linear models