Volume 56, Number 3, May-June 2022
|Page(s)||1053 - 1080|
|Published online||13 May 2022|
On Lyapunov stability of positive and conservative time integrators and application to second order modified Patankar–Runge–Kutta schemes
Department of Mathematics, University of Kassel, Kassel, Germany
* Corresponding author: email@example.com
Accepted: 22 March 2022
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.
Mathematics Subject Classification: 65L05 / 65L06 / 65L20
Key words: Modified Patankar–Runge–Kutta schemes / production-destruction systems / unconditionally positive and conservative schemes / Lyapunov stability analysis / center manifold theorem for maps
© The authors. Published by EDP Sciences, SMAI 2022
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.