A high-order / low-order (HOLO) algorithm for preserving conservation in time-dependent low-rank transport calculations

Abstract

Dynamical low-rank (DLR) approximation methods have previously been developed for time-dependent radiation transport problems. One crucial drawback of DLR is that it does not conserve important quantities of the calculation, which limits the applicability of the method. Here we address this conservation issue by solving a low-order equation with closure terms computed via a high-order solution calculated with DLR. We observe that the high-order solution well approximates the closure term, and the low-order solution can be used to correct the conservation bias in the DLR evolution. We also apply the linear discontinuous Galerkin method for the spatial discretization to obtain the asymptotic limit. We then demonstrate with the numerical results that this so-called high-order / low-order (HOLO) algorithm is conservative without sacrificing computational efficiency and accuracy.

References in corpus (9)

• Conservative model reduction for finite-volume models
• A New Spherical Harmonics Scheme for Multi-Dimensional Radiation Transport I: Static Matter Configurations
• A low-rank method for two-dimensional time-dependent radiation transport calculations
• A mass, momentum, and energy conservative dynamical low-rank scheme for the Vlasov equation
• A low-rank projector-splitting integrator for the Vlasov--Maxwell equations with divergence correction
• Implicit Filtered PN for High-Energy Density Thermal Radiation Transport using Discontinuous Galerkin Finite Elements
• Scalable angular adaptivity for Boltzmann transport
• Moment Closures Based on Minimizing the Residual of the P${}_N$ Angular Expansion in Radiation Transport
• Low-memory, discrete ordinates, discontinuous Galerkin methods for radiative transport

Cited by in corpus (5)

• A mass, momentum, and energy conservative dynamical low-rank scheme for the Vlasov equation
• A low-rank power iteration scheme for neutron transport criticality problems
• Efficient 6D Vlasov simulation using the dynamical low-rank framework Ensign
• Conservative DG Method for the Micro-Macro Decomposition of the Vlasov-Poisson-Lenard-Bernstein Model
• Multi-Frequency Implicit Semi-analog Monte-Carlo (ISMC) Radiative Transfer Solver in Two-Dimensions (without Teleportation)