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

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, the first order MPE scheme and the second order MPRK22(1) scheme, have been successfully applied in a large number of applications. Recently, a general definition of MPRK schemes was introduced and necessary as well as sufficient conditions to obtain first and second order MPRK schemes were presented. In this paper we derive necessary and sufficient conditions for third order MPRK schemes and introduce the first family of such schemes. The theoretical results are confirmed by numerical experiments considering linear and nonlinear as well as nonstiff and stiff systems of differential equations.

[1]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[2]  Eric Deleersnijder,et al.  A high-order conservative Patankar-type discretisation for stiff systems of production-destruction equations , 2003 .

[3]  Inga Hense,et al.  The representation of cyanobacteria life cycle processes in aquatic ecosystem models , 2010 .

[4]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[5]  E. Hairer,et al.  Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

[6]  Anna Scotti,et al.  Positivity and Conservation Properties of Some Integration Schemes for Mass Action Kinetics , 2011, SIAM J. Numer. Anal..

[7]  Andreas Meister,et al.  An improved and generalized second order, unconditionally positive, mass conserving integration scheme for biochemical systems , 2008 .

[8]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[9]  Luca Bonaventura,et al.  Unconditionally Strong Stability Preserving Extensions of the TR-BDF2 Method , 2016, Journal of Scientific Computing.

[10]  Chia-Jung Hsu.  Numerical Heat Transfer and Fluid Flow , 1981 .

[11]  Helmut Podhaisky,et al.  Numerik gewöhnlicher Differentialgleichungen , 2012 .

[12]  Patankar-Type Runge-Kutta Schemes for Linear PDEs , 2016, 1610.02715.

[13]  Andreas Meister,et al.  On unconditionally positive implicit time integration for the DG scheme applied to shallow water flows , 2014 .

[14]  Inga Hense,et al.  Modelling cyanobacteria in shallow coastal seas. , 2010 .

[15]  A positive and multi-element conserving time stepping scheme for biogeochemical processes in marine ecosystem models , 2015 .

[16]  Bob W. Kooi,et al.  A second-order, unconditionally positive, mass-conserving integration scheme for biochemical systems , 2006 .

[17]  Andreas Meister,et al.  Phosphorus Cycles in Lakes and Rivers: Modeling, Analysis, and Simulation , 2010 .

[18]  Inga Hense,et al.  Modelling the life cycle of dinoflagellates: a case study with Biecheleria baltica , 2013 .

[19]  Z. Horváth,et al.  Positivity of Runge-Kutta and diagonally split Runge-Kutta methods , 1998 .

[20]  Oliver Gressel.  Toward realistic simulations of magneto-thermal winds from weakly-ionized protoplanetary disks , 2017 .

[21]  Eric Deleersnijder,et al.  Application of modified Patankar schemes to stiff biogeochemical models for the water column , 2005 .

[22]  S. Kopecz,et al.  On Order Conditions for modified Patankar-Runge-Kutta schemes , 2017, 1702.04589.

[23]  David I. Ketcheson,et al.  Strong stability preserving runge-kutta and multistep time discretizations , 2011 .

[24]  J. Klar,et al.  A detailed view of filaments and sheets in the warm-hot intergalactic medium. I. Pancake formation , 2010, 1008.2311.

[25]  G. Nicolis,et al.  Chemical instabilities and sustained oscillations. , 1971, Journal of theoretical biology.

[26]  Rosenbrock methods in biogeochemical modelling – A comparison to Runge–Kutta methods and modified Patankar schemes , 2011 .

[27]  Andreas Meister,et al.  Description of a flexible and extendable physical–biogeochemical model system for the water column , 2006 .

[28]  Ashu Dastoor,et al.  Development of a global ocean mercury model with a methylation cycle: Outstanding issues , 2017 .

[29]  Angelika Fruehauf,et al.  A First Course In Numerical Analysis , 2016 .

[30]  Kellen Petersen August.  Real Analysis , 2009 .

[31]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..