On Lyapunov Stability of Positive and Conservative Time Integrators and Application to Second Order Modified Patankar–Runge–Kutta Schemes

by   Thomas Izgin, et al.

Since almost twenty years, modified Patankar–Runge–Kutta (MPRK) methods have proven to be efficient and robust numerical schemes that preserve positivity and conservativity of the production-destruction system irrespectively of the time step size chosen. Due to these advantageous properties they are used for a wide variety of applications. Nevertheless, until now, an analytic investigation of the stability of MPRK schemes is still missing, since the usual approach by means of Dahlquist's equation is not feasible. Therefore, we consider a positive and conservative 2D test problem and provide statements usable for a stability analysis of general positive and conservative time integrator schemes based on the center manifold theory. We use this approach to investigate the Lyapunov stability of the second order MPRK22(α) and MPRK22ncs(α) schemes. We prove that MPRK22(α) schemes are unconditionally stable and derive the stability regions of MPRK22ncs(α) schemes. Finally, numerical experiments are presented, which confirm the theoretical results.


On the Stability of Modified Patankar Methods

Patankar schemes have attracted more and more interests as a time-integr...

On the Stability of Unconditionally Positive and Linear Invariants Preserving Time Integration Schemes

Higher-order time integration methods that unconditionally preserve the ...

Qualitative properties of space-dependent SIR models with constant delay and their numerical solutions

In this article a space-dependent epidemic model equipped with a constan...

gPAV-Based Unconditionally Energy-Stable Schemes for the Cahn-Hilliard Equation: Stability and Error Analysis

We present several first-order and second-order numerical schemes for th...

Stable and accurate numerical methods for generalized Kirchhoff-Love plates

Efficient and accurate numerical algorithms are developed to solve a gen...

Local Multiple Traces Formulation for Electromagnetics: Stability and Preconditioning for Smooth Geometries

We consider the time-harmonic electromagnetic transmission problem for t...