Aufsatz
On the non-global linear stability and spurious fixed points of MPRK schemes with negative RK parameters
Abstract
Recently, a stability theory has been developed to study the linear stability of modified Patankar–Runge–Kutta (MPRK) schemes. This stability theory provides sufficient conditions for a fixed point of an MPRK scheme to be stable as well as for the convergence of an MPRK scheme towards the steady state of the corresponding initial value problem, whereas the main assumption is that the initial value is sufficiently close to the steady state. Initially, numerical experiments in several publications indicated that these linear stability properties are not only local but even global, as is the case for general linear methods. Recently, however, it was discovered that the linear stability of theMPDeC(8) scheme is indeed only local in nature. Our conjecture is that this is a result of negative Runge–Kutta (RK) parameters ofMPDeC(8) and that linear stability is indeed global if the RK parameters are nonnegative. To support this conjecture, we examine the family of MPRK22(α) methods with negative RK parameters and show that even among these methods there are methods for which the stability properties are only local. However, this local linear stability is not observed for MPRK22(α) schemes with nonnegative Runge–Kutta parameters. In particular, it is shown that MPRK22(α) schemes with 0 < α < 0.5 or −0.5 < α < 0 are only stable if the time step size is sufficiently small. But schemes with α ≤ −0.5 are stable and converge towards the steady state of the initial value problems for all time step sizes, at least for the test problem under consideration. Furthermore, it is shown that for some of the latter systems, the initial values must actually be close enough to the steady state to guarantee this result.
Sponsorship
Gefördert im Rahmen des Projekts DEALCitation
@article{doi:10.17170/kobra-2024080210628,
author={Izgin, Thomas and Kopecz, Stefan and Meister, Andreas and Schilling, Amandine},
title={On the non-global linear stability and spurious fixed points of MPRK schemes with negative RK parameters},
journal={Numerical Algorithms},
year={2024}
}
0500 Oax
0501 Text $btxt$2rdacontent
0502 Computermedien $bc$2rdacarrier
1100 2024$n2024
1500 1/eng
2050 ##0##http://hdl.handle.net/123456789/15954
3000 Izgin, Thomas
3010 Kopecz, Stefan
3010 Meister, Andreas
3010 Schilling, Amandine
4000 On the non-global linear stability and spurious fixed points of MPRK schemes with negative RK parameters / Izgin, Thomas
4030
4060 Online-Ressource
4085 ##0##=u http://nbn-resolving.de/http://hdl.handle.net/123456789/15954=x R
4204 \$dAufsatz
4170
5550 {{Runge-Kuta-Verfahren}}
5550 {{Stabilität}}
5550 {{Lineares Modell}}
7136 ##0##http://hdl.handle.net/123456789/15954
<resource xsi:schemaLocation="http://datacite.org/schema/kernel-2.2 http://schema.datacite.org/meta/kernel-2.2/metadata.xsd"> 2024-08-02T11:21:38Z 2024-08-02T11:21:38Z 2024-02-13 doi:10.17170/kobra-2024080210628 http://hdl.handle.net/123456789/15954 Gefördert im Rahmen des Projekts DEAL eng Namensnennung 4.0 International http://creativecommons.org/licenses/by/4.0/ Conservative scheme Unconditional positivity Modified Patankar–Runge–Kutta methods Linear stability 510 On the non-global linear stability and spurious fixed points of MPRK schemes with negative RK parameters Aufsatz Recently, a stability theory has been developed to study the linear stability of modified Patankar–Runge–Kutta (MPRK) schemes. This stability theory provides sufficient conditions for a fixed point of an MPRK scheme to be stable as well as for the convergence of an MPRK scheme towards the steady state of the corresponding initial value problem, whereas the main assumption is that the initial value is sufficiently close to the steady state. Initially, numerical experiments in several publications indicated that these linear stability properties are not only local but even global, as is the case for general linear methods. Recently, however, it was discovered that the linear stability of theMPDeC(8) scheme is indeed only local in nature. Our conjecture is that this is a result of negative Runge–Kutta (RK) parameters ofMPDeC(8) and that linear stability is indeed global if the RK parameters are nonnegative. To support this conjecture, we examine the family of MPRK22(α) methods with negative RK parameters and show that even among these methods there are methods for which the stability properties are only local. However, this local linear stability is not observed for MPRK22(α) schemes with nonnegative Runge–Kutta parameters. In particular, it is shown that MPRK22(α) schemes with 0 < α < 0.5 or −0.5 < α < 0 are only stable if the time step size is sufficiently small. But schemes with α ≤ −0.5 are stable and converge towards the steady state of the initial value problems for all time step sizes, at least for the test problem under consideration. Furthermore, it is shown that for some of the latter systems, the initial values must actually be close enough to the steady state to guarantee this result. open access Izgin, Thomas Kopecz, Stefan Meister, Andreas Schilling, Amandine doi:10.1007/s11075-024-01770-7 Runge-Kuta-Verfahren Stabilität Lineares Modell publishedVersion eissn:1572-9265 Issue 3 Numerical Algorithms 1221 - 1242 Volume 96 false </resource>
The following license files are associated with this item:
Except where otherwise noted, this item's license is described as Namensnennung 4.0 International