Abstract:This paper introduces a novel computational approach termed the Reduced Augmentation Implicit Low-rank (RAIL) method by investigating two predominant research directions in low-rank solutions to time-dependent partial differential equations (PDEs): dynamical low-rank (DLR), and step and truncation (SAT) tensor methods. The RAIL method, along with the development of the SAT approach, is designed to enhance the efficiency of traditional full-rank implicit solvers from method-of-lines discretizations of time-dependent PDEs, while maintaining accuracy and stability. We consider spectral methods for spatial discretization, and diagonally implicit Runge-Kutta (DIRK) and implicit-explicit (IMEX) RK methods for time discretization. The efficiency gain is achieved by investigating low-rank structures within solutions at each RK stage using a singular value decomposition (SVD). In particular, we develop a reduced augmentation procedure to predict the basis functions to construct projection subspaces. This procedure balances algorithm accuracy and efficiency by incorporating as many bases as possible from previous RK stages and predictions, and by optimizing the basis representation through SVD truncation. As such, one can form implicit schemes for updating basis functions in a dimension-by-dimension manner, similar in spirit to the K-L step in the DLR framework. We also apply a globally mass conservative post-processing step at the end of each RK stage. We validate the RAIL method through numerical simulations of advection-diffusion problems and a Fokker-Planck model, showcasing its ability to efficiently handle time-dependent PDEs while maintaining global mass conservation. Our approach generalizes and bridges the DLR and SAT approaches, offering a comprehensive framework for efficiently and accurately solving time-dependent PDEs with implicit treatment.

Sylvester-Preconditioned Adaptive-Rank Implicit Time Integrators for Advection-Diffusion Equations with Inhomogeneous Coefficients

Krylov-based Adaptive-Rank Implicit Time Integrators for Stiff Problems with Application to Nonlinear Fokker-Planck Kinetic Models

Implicit low-rank Riemannian schemes for the time integration of stiff partial differential equations

High order implicit time integration schemes on multiresolution adaptive grids for stiff PDEs

CUR for Implicit Time Integration of Random Partial Differential Equations on Low-Rank Matrix Manifolds

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

Efficient ADI Schemes and Preconditioning for a Class of High-Dimensional Spatial Fractional Diffusion Equations with Variable Diffusion Coefficients

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

Adaptively Implicit Advection for Atmospheric Flows

Preconditioning Low Rank Generalized Minimal Residual Method (GMRES) for Implicit Discretizations of Matrix Differential Equations

A rank-adaptive robust integrator for dynamical low-rank approximation

Adaptive rational Krylov methods for exponential Runge--Kutta integrators

Adaptive Rational Krylov Methods for Exponential Runge–Kutta Integrators

Reduced Augmentation Implicit Low-rank (RAIL) integrators for advection-diffusion and Fokker-Planck models

A parallel rank-adaptive integrator for dynamical low-rank approximation

A Semi-Lagrangian Adaptive-Rank (SLAR) Method for Linear Advection and Nonlinear Vlasov-Poisson System

Coordinate-Adaptive Integration of PDEs on Tensor Manifolds

An ADI Iteration Method for Solving Discretized Two-Dimensional Space-Fractional Diffusion Equations

High order multiscale methods for advection-diffusion equation in highly oscillatory regimes: application to surfactant diffusion and generalization to arbitrary domains

Multi-Level Solution Strategies for Implicit Time Discontinuous Galerkin Discretizations

A space-time adaptive low-rank method for high-dimensional parabolic partial differential equations