Unconditionally positive and conservative third order modified Patankar-Runge-Kutta discretizations of production-destruction systems

DOI10.1007/S10543-018-0705-1zbMATH Open1397.65102arXiv1703.05052OpenAlexW2603908784MaRDI QIDQ1783380FDOQ1783380

Authors: Stefan Kopecz, Andreas Meister 

Publication date: 20 September 2018

Published in: BIT (Search for Journal in Brave)

Abstract: Modified Patankar-Runge-Kutta (MPRK) schemes are numerical methods for the solution of positive and conservative production-destruction systems. They adapt explicit Runge-Kutta schemes to ensure positivity and conservation irrespective of the time step size. The first two members of this class were introduced in [Burchard et. al.: A high-order conservative Patankar-type discretisation for stiff systems of production-destruction equations. Appl. Numer. Math., 47(1):1-30, 2003] and have been successfully applied in a large number of applications. Recently, we introduced a general definition of MPRK schemes and presented a thorough investigation of first and second order MPRK schemes in [Kopecz. S, Meister, A.: On order conditions for modified patankar-runge-kutta schemes. arXiv:1702.04589 [math.NA], 2017.]. A potentially third order Patankar-type method was introduced in [Formaggia L., Scotti, A.: Positivity and conservation properties of some integration schemes for mass action kinetics. SIAM J. Numer. Anal., 49(3):1267-1288, 2011.]. This method uses the MPRK22(1) scheme of [Burchard et. al.: A high-order conservative Patankar-type discretisation for stiff systems of production-destruction equations. Appl. Numer. Math., 47(1):1-30, 2003] as a predictor and a modification of the BDF(3) multistep method as a corrector. This method is at most third order and at least second order accurate. In this paper we continue the work of [Kopecz. S, Meister, A.: On order conditions for modified patankar-runge-kutta schemes. arXiv:1702.04589 [math.NA], 2017.] and present necessary and sufficient conditions for third order MPRK schemes. For the first time, we introduce MPRK schemes, which are third order accurate independent of the specific positive and conservative system under consideration. The theoretical results are confirmed by numerical experiments.

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

Recommendations

• On the existence of three-stage third-order modified Patankar-Runge-Kutta schemes
• On order conditions for modified Patankar-Runge-Kutta schemes
• On Lyapunov stability of positive and conservative time integrators and application to second order modified Patankar-Runge-Kutta schemes
• Extension of modified Patankar-Runge-Kutta schemes to nonautonomous production-destruction systems based on Oliver's approach
• A third-order unconditionally positivity-preserving scheme for production-destruction equations with applications to non-equilibrium flows

zbMATH Keywords

conservationunconditional positivityPatankar-type schemesmodified Patankar-Runge-Kutta 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

• Efficient implementation of essentially nonoscillatory shock-capturing schemes
• Title not available (Why is that?)
• Solving Ordinary Differential Equations I
• Title not available (Why is that?)
• A second-order, unconditionally positive, mass-conserving integration scheme for biochemical systems
• Total variation diminishing Runge-Kutta schemes
• A first course in numerical analysis.
• Positivity of Runge-Kutta and diagonally split Runge-Kutta methods
• 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
• Chemical instabilities and sustained oscillations
• Title not available (Why is that?)
• Numerical methods for ordinary differential equations. Nonstiff, stiff, and differential-algebraic equations
• An improved and generalized second order, unconditionally positive, mass conserving integration scheme for biochemical systems
• Unconditionally strong stability preserving extensions of the TR-BDF2 method
• Phosphorus cycles in lakes and rivers: modeling, analysis, and simulation
• On unconditionally positive implicit time integration for the DG scheme applied to shallow water flows
• 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

Cited In (27)

• Positivity-preserving and elementary stable nonstandard method for a COVID-19 SIR model
• The effect of splitting strategy on qualitative property preservation
• 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
• Third-order conservative sign-preserving and steady-state-preserving time integrations and applications in stiff multispecies and multireaction detonations
• On the existence of three-stage third-order modified Patankar-Runge-Kutta schemes
• Extension of modified Patankar-Runge-Kutta schemes to nonautonomous production-destruction systems based on Oliver's approach
• High-order bound-preserving finite difference methods for multispecies and multireaction detonations
• Second-order unconditional positive preserving schemes for non-equilibrium reactive flows with mass and mole balance
• Efficient inequality-preserving integrators for differential equations satisfying forward Euler conditions
• 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
• On the stability of strong-stability-preserving modified Patankar–Runge–Kutta schemes
• Positivity-preserving methods for ordinary differential equations
• On order conditions for modified Patankar-Runge-Kutta schemes
• On the Stability of Unconditionally Positive and Linear Invariants Preserving Time Integration Schemes
• Qualitative property preservation of high-order operator splitting for the SIR model
• Two Nonlinear Positivity-Preserving Finite Volume Schemes for Three-Dimensional Heat Conduction Equations on General Polyhedral Meshes
• Positivity-preserving adaptive Runge-Kutta methods
• A third-order unconditionally positivity-preserving scheme for production-destruction equations with applications to non-equilibrium flows
• On Lyapunov stability of positive and conservative time integrators and application to second order modified Patankar-Runge-Kutta schemes
• An arbitrary high order and positivity preserving method for the shallow water equations
• Issues with positivity-preserving Patankar-type schemes
• A high-order conservative Patankar-type discretisation for stiff systems of production--destruction equations
• A comprehensive theory on generalized BBKS schemes
• High-order bound-preserving discontinuous Galerkin methods for stiff multispecies detonation

This page was built for publication: Unconditionally positive and conservative third order modified Patankar-Runge-Kutta discretizations of production-destruction systems

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