Abstract:We develop two new sets of stable, rank-adaptive Dynamically Orthogonal Runge-Kutta (DORK) schemes that capture the high-order curvature of the nonlinear low-rank manifold. The DORK schemes asymptotically approximate the truncated singular value decomposition at a greatly reduced cost while preserving mode continuity using newly derived retractions. We show that arbitrarily high-order optimal perturbative retractions can be obtained, and we prove that these new retractions are stable. In addition, we demonstrate that repeatedly applying retractions yields a gradient-descent algorithm on the low-rank manifold that converges superlinearly when approximating a low-rank matrix. When approximating a higher-rank matrix, iterations converge linearly to the best low-rank approximation. We then develop a rank-adaptive retraction that is robust to overapproximation. Building off of these retractions, we derive two rank-adaptive integration schemes that dynamically update the subspace upon which the system dynamics are projected within each time step: the stable, optimal Dynamically Orthogonal Runge-Kutta (so-DORK) and gradient-descent Dynamically Orthogonal Runge-Kutta (gd-DORK) schemes. These integration schemes are numerically evaluated and compared on an ill-conditioned matrix differential equation, an advection-diffusion partial differential equation, and a nonlinear, stochastic reaction-diffusion partial differential equation. Results show a reduced error accumulation rate with the new stable, optimal and gradient-descent integrators. In addition, we find that rank adaptation allows for highly accurate solutions while preserving computational efficiency.

Diagonally drift-implicit Runge-Kutta methods of weak order one and two for Itô SDEs and stability analysis

Construction and mean-square stability analysis of a new family of stochastic Runge-Kutta methods

A brief review on stability investigations of numerical methods for systems of stochastic differential equations

An Exponential Stochastic Runge-Kutta Type Method of Order up to 1.5 for SPDEs of Nemytskii-type

Algebraically Stable SDIRK Methods with Controllable Numerical Dissipation for First/second-Order Time-Dependent Problems

Convergence Aspects of Step-Parallel Iteration of Runge-Kutta Methods for Delay Differential Equations

Continuous Stage Stochastic Runge–Kutta Methods

A study on generalized balanced split drift stochastic Runge- Kutta methods for stochastic differential equations

Convergence of the tamed-Euler-Maruyama method for SDEs with discontinuous and polynomially growing drift

Structure-preserving Stochastic Runge–Kutta–Nyström Methods for Nonlinear Second-Order Stochastic Differential Equations with Multiplicative Noise

The Euler-Maruyama Scheme for SDEs with Irregular Drift: Convergence Rates via Reduction to a Quadrature Problem

Analysis for Implicit and Implicit-Explicit ADER and DeC Methods for Ordinary Differential Equations, Advection-Diffusion and Advection-Dispersion Equations

Stochastic C-stability and B-consistency of explicit and implicit Euler-type schemes

Stability-Optimized High Order Methods and Stiffness Detection for Pathwise Stiff Stochastic Differential Equations

Runge-Kutta methods are stable

Three-stage Semi-implicit Stochastic Runge-Kutta Methods for Stochastic Differential Equations

A strong order $3/4$ method for SDEs with discontinuous drift coefficient

Stable rank-adaptive Dynamically Orthogonal Runge-Kutta schemes

Stability Analysis of a Class of Parallel Runge-Kutta Methods for Delay Differential Equations

Necessary and sufficient conditions for strong stability of explicit Runge-Kutta methods

Strong stability preserving two-step Runge-Kutta methods