Source author record

Stefan Kopecz

Stefan Kopecz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

3works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

3 published item(s)

preprint2022arXiv

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

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.

preprint2022arXiv

On the stability of strong-stability-preserving modified Patankar Runge-Kutta schemes

In this paper, we perform stability analysis for a class of second and third order accurate strong-stability-preserving modified Patankar Runge-Kutta (SSPMPRK) schemes, which were introduced in [4,5] and can be used to solve convection equations with stiff source terms, such as reactive Euler equations, with guaranteed positivity under the standard CFL condition due to the convection terms only. The analysis allows us to identify the range of free parameters in these SSPMPRK schemes in order to ensure stability. Numerical experiments are provided to demonstrate the validity of the analysis.

preprint2022arXiv

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

Higher-order time integration methods that unconditionally preserve the positivity and linear invariants of the underlying differential equation system cannot belong to the class of general linear methods. This poses a major challenge for the stability analysis of such methods since the new iterate depends nonlinearly on the current iterate. Moreover, for linear systems, the existence of linear invariants is always associated with zero eigenvalues, so that steady states of the continuous problem become non-hyperbolic fixed points of the numerical time integration scheme. Altogether, the stability analysis of such methods requires the investigation of non-hyperbolic fixed points for general nonlinear iterations. Based on the center manifold theory for maps we present a theorem for the analysis of the stability of non-hyperbolic fixed points of time integration schemes applied to problems whose steady states form a subspace. This theorem provides sufficient conditions for both the stability of the method and the local convergence of the iterates to the steady state of the underlying initial value problem. This theorem is then used to prove the unconditional stability of the MPRK22($α$)-family of modified Patankar-Runge-Kutta schemes when applied to arbitrary positive and conservative linear systems of differential equations. The theoretical results are confirmed by numerical experiments.