On Lyapunov stability of positive and conservative time integrators and application to second order modified Patankar-Runge-Kutta schemes

DOI10.1051/M2AN/2022031zbMATH Open1492.65191arXiv2202.01099OpenAlexW4220769286MaRDI QIDQ5086979FDOQ5086979

Authors: Thomas Izgin, Stefan Kopecz, Andreas Meister 

Publication date: 8 July 2022

Published in: ESAIM: Mathematical Modelling and Numerical Analysis (Search for Journal in Brave)

Abstract: 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($alpha$) and MPRK22ncs($alpha$) schemes. We prove that MPRK22($alpha$) schemes are unconditionally stable and derive the stability regions of MPRK22ncs($alpha$) schemes. Finally, numerical experiments are presented, which confirm the theoretical results.

Full work available at URL: https://arxiv.org/abs/2202.01099

Recommendations

• On the Stability of Unconditionally Positive and Linear Invariants Preserving Time Integration Schemes
• Unconditionally positive and conservative third order modified Patankar-Runge-Kutta discretizations of production-destruction systems
• On the stability of strong-stability-preserving modified Patankar–Runge–Kutta schemes
• On order conditions for modified Patankar-Runge-Kutta schemes
• On the existence of three-stage third-order modified Patankar-Runge-Kutta schemes

zbMATH Keywords

Lyapunov stability analysismodified Patankar-Runge-Kutta schemescenter manifold theorem for mapsproduction-destruction systemsunconditionally positive and conservative schemes

Mathematics Subject Classification ID

Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations (65L06) Stability and convergence of numerical methods for ordinary differential equations (65L20) Numerical methods for initial value problems involving ordinary differential equations (65L05)

Cites Work

• The Hopf bifurcation and its applications. With contributions by P. Chernoff, G. Childs, S. Chow, J. R. Dorroh, J. Guckenheimer, L. Howard, N. Kopell, O. Lanford, J. Mallet-Paret, G. Oster, O. Ruiz, S. Schecter, D. Schmidt, and S. Smale
• On the Inverse of the Sum of Matrices
• Applications of centre manifold theory
• Title not available (Why is that?)
• Title not available (Why is that?)
• Title not available (Why is that?)
• A high-order conservative Patankar-type discretisation for stiff systems of production--destruction equations
• Positivity and conservation properties of some integration schemes for mass action kinetics
• Unconditionally strong stability preserving extensions of the TR-BDF2 method
• Phosphorus cycles in lakes and rivers: modeling, analysis, and simulation
• A conservative, positivity preserving scheme for advection-diffusion-reaction equations in biochemical applications
• On order conditions for modified Patankar-Runge-Kutta schemes
• A detailed view of filaments and sheets in the warm-hot intergalactic medium. I: Pancake formation
• Unconditionally positive and conservative third order modified Patankar-Runge-Kutta discretizations of production-destruction systems
• Arbitrary high-order, conservative and positivity preserving Patankar-type deferred correction schemes
• Positivity-preserving time discretizations for production-destruction equations with applications to non-equilibrium flows
• On the existence of three-stage third-order modified Patankar-Runge-Kutta schemes
• \textit{GeCo}: geometric conservative nonstandard schemes for biochemical systems
• Qualitative property preservation of high-order operator splitting for the SIR model
• A third-order unconditionally positivity-preserving scheme for production-destruction equations with applications to non-equilibrium flows
• Positivity-preserving adaptive Runge-Kutta methods
• A comprehensive theory on generalized BBKS schemes
• Extension of modified Patankar-Runge-Kutta schemes to nonautonomous production-destruction systems based on Oliver's approach

Cited In (10)

• Extension of modified Patankar-Runge-Kutta schemes to nonautonomous production-destruction systems based on Oliver's approach
• A study of the local dynamics of modified Patankar DeC and higher order modified Patankar–RK methods
• On the dynamics of first and second order GeCo and gBBKS schemes
• A necessary condition for non-oscillatory and positivity preserving time-integration schemes
• Unconditionally positive and conservative third order modified Patankar-Runge-Kutta discretizations of production-destruction systems
• On the stability of strong-stability-preserving modified Patankar–Runge–Kutta schemes
• On order conditions for modified Patankar-Runge-Kutta schemes
• On the Stability of Unconditionally Positive and Linear Invariants Preserving Time Integration Schemes
• An arbitrary high order and positivity preserving method for the shallow water equations
• Issues with positivity-preserving Patankar-type schemes

This page was built for publication: On Lyapunov stability of positive and conservative time integrators and application to second order modified Patankar-Runge-Kutta schemes

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5086979)