A study of the local dynamics of modified Patankar DeC and higher order modified Patankar-RK methods. (English) Zbl 1525.65067
Summary: Patankar schemes have attracted increasing interest in recent years because they preserve the positivity of the analytical solution of a production-destruction system (PDS) irrespective of the chosen time step size. Although they are now of great interest, for a long time it was not clear what stability properties such schemes have. Recently a new stability approach based on Lyapunov stability with an extension of the center manifold theorem has been proposed to study the stability properties of positivity-preserving time integrators. In this work, we study the stability properties of the classical modified Patankar-Runge-Kutta schemes (MPRK) and the modified Patankar Deferred Correction (MPDeC) approaches. We prove that most of the considered MPRK schemes are stable for any time step size and compute the stability function of MPDeC. We investigate its properties numerically revealing that also most MPDeC are stable irrespective of the chosen time step size. Finally, we verify our theoretical results with numerical simulations.
MSC:
| 65L06 | Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations |
| 65L05 | Numerical methods for initial value problems involving ordinary differential equations |
| 65L20 | Stability and convergence of numerical methods for ordinary differential equations |
Keywords:
Lyapunov stability; conservative; positivity-preserving; modified Patankar schemes; Runge-Kutta; deferred correctionSoftware:
GitHubReferences:
| [1] | R. Abgrall, P. Öffner and H. Ranocha, Reinterpretation and extension of entropy correction terms for residual distribution and discontinuous Galerkin schemes: application to structure preserving discretization. J. Comput. Phys. 453 (2022) 24. · Zbl 07517718 |
| [2] | R. Abgrall, E.L. Mélédo, P. Öffner and D. Torlo, Relaxation deferred correction methods and their applications to residual distribution schemes. SMAI J. Comput. Math. 8 (2022) 125-160. · Zbl 1501.65061 · doi:10.5802/smai-jcm.82 |
| [3] | L. Benvenuti and L. Farina, Eigenvalue regions for positive systems. Syst. Control Lett. 51 (2004) 325-330. · Zbl 1157.93401 · doi:10.1016/j.sysconle.2003.09.009 |
| [4] | H. Burchard, E. Deleersnijder and A. Meister, A high-order conservative Patankar-type discretisation for stiff systems of production-destruction equations. Appl. Numer. Math. 47 (2003) 1-30. · Zbl 1028.80008 |
| [5] | A. Chertock, S. Cui, A. Kurganov and T. Wu, Steady state and sign preserving semi-implicit Runge-Kutta methods for ODEs with stiff damping term. SIAM J. Numer. Anal. 53 (2015) 2008-2029. · Zbl 1327.65128 · doi:10.1137/151005798 |
| [6] | M. Ciallella, L. Micalizzi, P. Öffner and D. Torlo, An arbitrary high order and positivity preserving method for the shallow water equations. Comput. Fluids 247 (2022) 21. · Zbl 1521.76408 |
| [7] | P. Deuflhard and F. Bornemann, Scientific computing with ordinary differential equations, in Texts Appl. Math.. Vol. 42. Springer, New York, NY (2002). · Zbl 1001.65071 |
| [8] | A. Dutt, L. Greengard and V. Rokhlin, Spectral deferred correction methods for ordinary differential equations. BIT 40 (2000) 241-266. · Zbl 0959.65084 · doi:10.1023/A:1022338906936 |
| [9] | M. Han Veiga, P. Öffner and D. Torlo, DeC and ADER: similarities, differences and a unified framework. J. Sci. Comput. 87 (2021) 35. · Zbl 1459.76086 · doi:10.1007/s10915-020-01397-5 |
| [10] | J. Huang and C.-W. Shu, Positivity-preserving time discretizations for production-destruction equations with applications to non-equilibrium flows. J. Sci. Comput. 78 (2019) 1811-1839. · Zbl 1420.35190 · doi:10.1007/s10915-018-0852-1 |
| [11] | J. Huang, W. Zhao and C.-W. Shu, A third-order unconditionally positivity-preserving scheme for production-destruction equations with applications to non-equilibrium flows. J. Sci. Comput. 79 (2019) 1015-1056. · Zbl 1444.35125 · doi:10.1007/s10915-018-0881-9 |
| [12] | J. Huang, T. Izgin, S. Kopecz, A. Meister and C.-W. Shu, On the stability of strong-stability-preserving modified Patankar-Runge-Kutta schemes. Preprint Preprint arxiv:2205.01488 (2022). · Zbl 1515.65192 |
| [13] | T. Izgin and P. Öffner, A study of the local dynamics of MPDeC and higher order MPRK methods (code). https://github.com/IzginThomas/MPRK43II.git (2023). |
| [14] | T. Izgin, S. Kopecz and A. Meister, Recent developments in the field of modified Patankar-Runge-Kutta-methods. PAMM 21 (2021) e202100027. · Zbl 1492.65191 · doi:10.1002/pamm.202100027 |
| [15] | T. Izgin, S. Kopecz and A. Meister, On lyapunov stability of positive and conservative time integrators and application to second order modified Patankar-Runge-Kutta schemes. ESAIM: M2AN 56 (2022) 1053-1080. · Zbl 1492.65191 · doi:10.1051/m2an/2022031 |
| [16] | T. Izgin, S. Kopecz and A. Meister, On the stability of unconditionally positive and linear invariants preserving time integration schemes. SIAM J. Numer. Anal. 60 (2022) 3029-3051. · Zbl 1506.65098 · doi:10.1137/22M1480318 |
| [17] | T. Izgin, P. Öffner and D. Torlo, Modified Patankar: Oscillations and Lyapunov Stability (code). https://github.com/accdavlo/Modified-Patankar-Oscillations-and-Lyapunov-Stability (2022). |
| [18] | S. Kopecz and A. Meister, Unconditionally positive and conservative third order modified Patankar-Runge-Kutta discretizations of production-destruction systems. BIT 58 (2018) 691-728. · Zbl 1397.65102 · doi:10.1007/s10543-018-0705-1 |
| [19] | S. Kopecz and A. Meister, On the existence of three-stage third-order modified Patankar-Runge-Kutta schemes. Numer. Algorithms 81 (2019) 1473-1484. · Zbl 1416.65205 |
| [20] | D.G. Luenberger, Introduction to Dynamic Systems. Theory, Models, and Applications. Vol. XIV. John Wiley & Sons, New York (1979) 446. · Zbl 0458.93001 |
| [21] | A. Meister and S. Ortleb, On unconditionally positive implicit time integration for the DG scheme applied to shallow water flows. Int. J. Numer. Methods Fluids 76 (2014) 69-94. · Zbl 1455.65173 · doi:10.1002/fld.3921 |
| [22] | P. Öffner and D. Torlo, Arbitrary high-order, conservative and positivity preserving Patankar-type deferred correction schemes. Appl. Numer. Math. 153 (2020) 15-34. · Zbl 1437.65073 |
| [23] | H. Ranocha, L. Lóczi and D.I. Ketcheson, General relaxation methods for initial-value problems with application to multistep schemes. Numer. Math. 146 (2020) 875-906. · Zbl 1456.65051 · doi:10.1007/s00211-020-01158-4 |
| [24] | E.M. Stein and R. Shakarchi, Complex analysis, in Princeton Lectures in Analysis. Vol. 2. Princeton University Press, Princeton, NJ (2003). · Zbl 1020.30001 |
| [25] | A. Stuart and A.R. Humphries, Dynamical Systems and Numerical Analysis. Vol. 2. Cambridge University Press (1998). · Zbl 0913.65068 |
| [26] | E.C. Titchmarsh, The Theory of Functions, second ed., Oxford University Press, Oxford (1939). · JFM 65.0302.01 |
| [27] | D. Torlo, P. Öffner and H. Ranocha, Issues with positivity-preserving Patankar-type schemes. Appl. Numer. Math. 182 (2022) 117-147. · Zbl 1506.65106 · doi:10.1016/j.apnum.2022.07.014 |
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. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
