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.

An extension of the low-rank Lyapunov ADI to non-zero initial values and its applications

Low-rank generalized alternating direction implicit iteration method for solving matrix equations

Low-Rank Methods for Solving Discrete-Time Projected Lyapunov Equations

A class of ADI method for the nonsymmetric coupled algebraic Riccati equations

Inexact linear solves in the low-rank ADI iteration for large Sylvester equations

Low-rank alternating direction doubling algorithm for solving large-scale continuous time algebraic Riccati equations

A general alternating direction implicit iteration method for solving complex symmetric linear systems

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

Alternately linearized implicit iteration methods for the minimal nonnegative solutions of the nonsymmetric algebraic Riccati equations

Linearized Alternating Direction Method with Adaptive Penalty for Low-Rank Representation.

Linearized Alternating Direction Method with Adaptive Penalty for Low-Rank Representation

Splitting ADI scheme for fractional Laplacian wave equations

Using $LDL^{T}$ factorizations in Newton's method for solving general large-scale algebraic Riccati equations

A Low-Rank ADMM Splitting Approach for Semidefinite Programming

An increasing rank Riemannian method for generalized Lyapunov equations

Riemannian Linearized Proximal Algorithms for Nonnegative Inverse Eigenvalue Problem.

An RADI-type method for stochastic continuous-time algebraic Riccati equations

Preconditioned Low-Rank Riemannian Optimization for Symmetric Positive Definite Linear Matrix Equations

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

Arbitrary high order ADER-DG method with local DG predictor for solutions of initial value problems for systems of first-order ordinary differential equations