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 low-rank method for time-dependent transport calculations

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

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

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

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

A Deterministic Dynamical Low-rank Approach for Charged Particle Transport

Asymptotic-Preserving and Energy Stable Dynamical Low-Rank Approximation

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

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

A Reduced Basis Method for Radiative Transfer Equation

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

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

High-order Adaptive Rank Integrators for Multi-scale Linear Kinetic Transport Equations in the Hierarchical Tucker Format

Solving optimal control problems governed by random Navier-Stokes equations using low-rank methods

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

Dynamical low-rank approximation of the Vlasov-Poisson equation with piecewise linear spatial boundary

Low-rank solutions to the stochastic Helmholtz equation

Low-rank approximation in the numerical modeling of the Farley-Buneman instability in ionospheric plasma

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

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

Low rank approximation method for perturbed linear systems with applications to elliptic type stochastic PDEs