Regresar

Extension of modified Patankar–Runge–Kutta schemes to nonautonomous production–destruction systems based on Oliver's approach

Abstract:

The mathematical modeling of various real life applications leads to systems of ordinary differential equations which include crucial properties like the positivity of the solution as well as the conservation of mass or energy. Based on the fundamental work of Burchard et al. (2003), unconditionally positive and conservative modified Patankar–Runge–Kutta schemes (MPRK) are available. These methods are highly stable and often outperform standard Runge–Kutta schemes. In this article, we extend MPRK methods, named MPRKO methods, using Oliver's approach (Oliver, 1975) to improve the accuracy of these schemes in the field of nonautonomous systems. The approach does not demand Ae=c in the Butcher tableau (A,b,c), where e=(1,…,1)T. Following the general analysis of MPRK schemes described in Kopecz and Meister (2018), positivity and mass conservation fundamental properties are proven and even conditions concerning the Patankar weights are given to get second order accuracy of the MPRKO methods. Finally, we consider different linear models and a non-linear epidemiological SEIR problem to confirm the theoretical results and to give reliable statements about the accuracy of the novel class of MPRKO methods.