Abstract:High-order accurate and stable explicit methods are powerful in solving differential equations efficiently. In this work, we propose a systematic framework to trade off accuracy for stability, especially the unconditional preservation of strong stability, positivity, range boundedness and contractivity. The whole algorithm consists of three steps: (1) Introducing a stabilizing term in the continuous system?(2) Integrating the system using an explicit exponential method? (3) Substituting the exponential functions with suitable approximations. We first show that a class of first-and second-order exponential time difference Runge-Kutta schemes are capable to preserve structures unconditionally when suitable stabilization parameter is chosen. Then by adopting the integrating factor approach with high-order Runge-Kutta and multi-step schemes as underlying schemes, three different approximation techniques are developed to make high-order schemes unconditionally structure-preserving, i.e., (1) a Taylor polynomial approximation? (2) a recursive approximation? (3) an approximation using combinations of exponential and linear functions. The proposed parametric schemes can be deployed to stiff problems straightforwardly by treating the stiff linear term as an integrating factor. The resulting time integration methods retain the explicitness and convergence orders of underlying time-marching schemes, yet with unconditional preservation of structures. The proposed framework using the second and third approximations has relatively mild requirement on underlying schemes, i.e., all coefficients are non-negative. Thus the parametric Runge-Kutta schemes can reach up to the fourth-order, and there is no order barrier in parametric multi-step schemes. The only free parameter-the stablization parameter in the framework can be determined a priori based on the forward Euler conditions. Unlike implicit methods, the parametric methodology allows for solving nonlinear problems stably and explicitly. As an alternative to conditionally structure-preserving methods, the proposed schemes are promising for the efficient computation of stiff and nonlinear problems. Numerical tests on benchmark problems with different stiffness are carried out to assess the performance of parametric methods.

Unconditionally positivity-preserving approximations of the Ait-Sahalia type model: Explicit Milstein-type schemes

Unconditionally positivity-preserving explicit Euler-type schemes for a generalized Ait-Sahalia model

An Explicit Positivity-Preserving Scheme for the Heston 3/2-Model with Order-One Strong Convergence

Mean-square convergence rates of implicit Milstein type methods for SDEs with non-Lipschitz coefficients

Antithetic multilevel Monte Carlo method for approximations of SDEs with non-globally Lipschitz continuous coefficients

Positivity-preserving truncated Euler and Milstein methods for financial SDEs with super-linear coefficients

Multilevel Monte Carlo methods for positivity-preserving approximations of the Heston 3/2-model

A General Framework to Construct High-order Unconditionally Structure-preserving Parametric Methods

Weak approximation schemes for SDEs with super-linearly growing coefficients

Positivity preserving logarithmic Euler-Maruyama type scheme for stochastic differential equations

An Advanced Numerical Scheme for Multi-Dimensional Stochastic Kolmogorov Equations with Superlinear Coefficients.

An extended Milstein scheme for effective weak approximation of diffusions

First order strong approximation of Ait-Sahalia-type interest rate model with Poisson jumps

Strongly asymptotically optimal schemes for the strong approximation of stochastic differential equations with respect to the supremum error

An Antithetic Multilevel Monte Carlo-Milstein Scheme for Stochastic Partial Differential Equations

Expectation-based Positivity-Preserving Noncompact Numerical Schemes

Approximation of the invariant measure for stable SDE by the Euler-Maruyama scheme with decreasing step-sizes

Weak error analysis for strong approximation schemes of SDEs with super-linear coefficients II: finite moments and higher-order schemes

Mittag-Leffler stability of complete monotonicity-preserving schemes for time-dependent coefficients sub-diffusion equations

Higher order time discretization method for a class of semilinear stochastic partial differential equations with multiplicative noise