On the non-global linear stability and spurious fixed points of MPRK schemes with negative RK parameters

dc.date.accessioned2024-08-02T11:21:38Z
dc.date.available2024-08-02T11:21:38Z
dc.date.issued2024-02-13
dc.description.sponsorshipGefördert im Rahmen des Projekts DEAL
dc.identifierdoi:10.17170/kobra-2024080210628
dc.identifier.urihttp://hdl.handle.net/123456789/15954
dc.language.isoengeng
dc.relation.doidoi:10.1007/s11075-024-01770-7
dc.rightsNamensnennung 4.0 International*
dc.rights.urihttp://creativecommons.org/licenses/by/4.0/*
dc.subjectConservative schemeeng
dc.subjectUnconditional positivityeng
dc.subjectModified Patankar–Runge–Kutta methodseng
dc.subjectLinear stabilityeng
dc.subject.ddc510
dc.subject.swdRunge-Kuta-Verfahrenger
dc.subject.swdStabilitätger
dc.subject.swdLineares Modellger
dc.titleOn the non-global linear stability and spurious fixed points of MPRK schemes with negative RK parameterseng
dc.typeAufsatz
dc.type.versionpublishedVersion
dcterms.abstractRecently, 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.eng
dcterms.accessRightsopen access
dcterms.creatorIzgin, Thomas
dcterms.creatorKopecz, Stefan
dcterms.creatorMeister, Andreas
dcterms.creatorSchilling, Amandine
dcterms.source.identifiereissn:1572-9265
dcterms.source.issueIssue 3
dcterms.source.journalNumerical Algorithmseng
dcterms.source.pageinfo1221 - 1242
dcterms.source.volumeVolume 96
kup.iskupfalse

Files

Original bundle

Now showing 1 - 1 of 1
Thumbnail Image
Name:
s11075_024_01770_7.pdf
Size:
623 KB
Format:
Adobe Portable Document Format
Description:

License bundle

Now showing 1 - 1 of 1
No Thumbnail Available
Name:
license.txt
Size:
3.03 KB
Format:
Item-specific license agreed upon to submission
Description:

Collections