×

zbMATH — the first resource for mathematics

Arbitrary high-order, conservative and positivity preserving Patankar-type deferred correction schemes. (English) Zbl 1437.65073
Summary: Production-destruction systems (PDS) of ordinary differential equations (ODEs) are used to describe physical and biological reactions in nature. The considered quantities are subject to natural laws. Therefore, they preserve positivity and conservation of mass at the analytical level. In order to maintain these properties at the discrete level, the so-called modified Patankar-Runge-Kutta (MPRK) schemes are often used in this context. However, up to our knowledge, the family of MPRK has been only developed up to third order of accuracy. In this work, we propose a method to solve PDS problems, but using the Deferred Correction (DeC) process as a time integration method. Applying the modified Patankar approach to the DeC scheme results in provable conservative and positivity preserving methods. Furthermore, we demonstrate that these modified Patankar DeC schemes can be constructed up to arbitrarily high order. Finally, we validate our theoretical analysis through numerical simulations.

MSC:
65L06 Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations
37N25 Dynamical systems in biology
92D25 Population dynamics (general)
PDF BibTeX XML Cite
Full Text: DOI
References:
[1] Abgrall, R., High order schemes for hyperbolic problems using globally continuous approximation and avoiding mass matrices, J. Sci. Comput., 73, 461-494 (2017) · Zbl 1398.65242
[2] Abgrall, R.; Torlo, D., Asymptotic preserving deferred correction residual distribution schemes (2018), arXiv preprint
[3] Abgrall, R.; Meledo, E.l.; Öffner, P., On the connection between residual distribution schemes and flux reconstruction (2018), arXiv preprint
[4] Burchard, H.; Deleersnijder, E.; Meister, A., A high-order conservative Patankar-type discretisation for stiff systems of production-destruction equations, Appl. Numer. Math., 47, 1-30 (2003) · Zbl 1028.80008
[5] Burchard, H.; Deleersnijder, E.; Meister, A., Application of modified Patankar schemes to stiff biogeochemical models for the water column, Ocean Dyn., 55, 326-337 (2005)
[6] Christlieb, A.; Ong, B.; Qiu, J.-M., Integral deferred correction methods constructed with high order Runge-Kutta integrators, Math. Comput., 79, 761-783 (2010) · Zbl 1209.65073
[7] Dutt, A.; Greengard, L.; Rokhlin, V., Spectral deferred correction methods for ordinary differential equations, BIT Numer. Math., 40, 241-266 (2000) · Zbl 0959.65084
[8] Gottlieb, S.; Ketcheson, D. I.; Shu, C.-W., Strong Stability Preserving Runge-Kutta and Multistep Time Discretizations (2011), World Scientific · Zbl 1241.65064
[9] Hairer, E.; Nørsett, S. P.; Wanner, G., Solving Ordinary Differential Equations. 1, Nonstiff Problems (1991), Springer-Vlg
[10] Hense, I.; Beckmann, A., The representation of cyanobacteria life cycle processes in aquatic ecosystem models, Ecol. Model., 221, 2330-2338 (2010)
[11] Huang, J.; Shu, C.-W., Positivity-preserving time discretizations for production-destruction equations with applications to non-equilibrium flows, J. Sci. Comput., 78, 1811-1839 (2019) · Zbl 1420.35190
[12] Huang, J.; Zhao, W.; Shu, C.-W., A third-order unconditionally positivity-preserving scheme for production-destruction equations with applications to non-equilibrium flows, J. Sci. Comput., 1-42 (2018)
[13] Ketcheson, D. I., Highly efficient strong stability-preserving Runge-Kutta methods with low-storage implementations, SIAM J. Sci. Comput., 30, 2113-2136 (2008) · Zbl 1168.65382
[14] Kopecz, S.; Meister, A., On order conditions for modified Patankar-Runge-Kutta schemes, Appl. Numer. Math., 123, 159-179 (2018) · Zbl 1377.65089
[15] Kopecz, S.; Meister, A., Unconditionally positive and conservative third order modified Patankar-Runge-Kutta discretizations of production-destruction systems, BIT Numer. Math., 1-38 (2018) · Zbl 1397.65102
[16] Kopecz, S.; Meister, A., On the existence of three-stage third-order modified Patankar-Runge-Kutta schemes, Numer. Algorithms, 1-12 (2019) · Zbl 1416.65205
[17] Liu, Y.; Shu, C.-W.; Zhang, M., Strong stability preserving property of the deferred correction time discretization, J. Comput. Math., 633-656 (2008) · Zbl 1174.65036
[18] Meister, A.; Ortleb, S., On unconditionally positive implicit time integration for the DG scheme applied to shallow water flows, Int. J. Numer. Methods Fluids, 76, 69-94 (2014)
[19] Minion, M. L., Semi-implicit spectral deferred correction methods for ordinary differential equations, Commun. Math. Sci., 1, 471-500 (2003) · Zbl 1088.65556
[20] Patankar, S., Numerical Heat Transfer and Fluid Flow (1980), CRC Press · Zbl 0521.76003
[21] Rackauckas, C.; Nie, Q., Differentialequations.jl-a performant and feature-rich ecosystem for solving differential equations in Julia, J. Open Res. Softw., 5 (2017)
[22] Ranocha, H.; Öffner, P., \( L_2\) stability of explicit Runge-Kutta schemes, J. Sci. Comput., 75, 1040-1056 (2018) · Zbl 1398.65188
[23] Ranocha, H.; Öffner, P.; Sonar, T., Summation-by-parts operators for correction procedure via reconstruction, J. Comput. Phys., 311, 299-328 (2016) · Zbl 1349.65524
[24] Wanner, G.; Hairer, E., Solving Ordinary Differential Equations II (1996), Springer Berlin Heidelberg · Zbl 0859.65067
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.