×

zbMATH — the first resource for mathematics

Collocation techniques for structured populations modeled by delay equations. (English) Zbl 1447.37077
Aguiar, Maira (ed.) et al., Current trends in dynamical systems in biology and natural sciences. Selected contributions presented at the ninth international workshop of dynamical systems applied to biology and natural sciences, DSABNS, Turin, Italy, February 7–9, 2018. Cham: Springer. SEMA SIMAI Springer Ser. 21, 43-62 (2020).
The authors in this paper study coupled equations such as the Daphnia model [O. Diekmann et al., J. Math. Biol. 61, No. 2, 277–318 (2010; Zbl 1208.92082)], which is an age-structured population model competing for an unstructured resource. This model, named after the organisms in the genus Daphnia, can be written as a coupled system of integral equations whose function terms are themselves solutions to an initial value problem ordinary differential equations.
The perspective is through continuation or homotopy methods. The primary problem stated by the authors is that in real-world models, solving the ordinary differential equations defining the terms in the coupled system is often a computationally inefficient procedure, since the ordinary differential equations are numerically resolved over each new iteration of the continuation parameter.
The authors combine this with polynomial collocation to quickly find solutions. Collocation is a method for solving ODEs which looks at a restricted space of functions such as polynomial (in polynomial collocation) that satisfy the ODE system at certain collocation points.
The main thrust of this paper is a series of demonstrations, starting with a simplified Daphnia model, on how one would implement polynomial collocation and continuation. They explore a variety of problems, including computing periodic solutions of models like the Daphnia model. In the penultimate section, explicit examples are demonstrated using the delayed logistic equation and a renewal equation for which an exact period solution is known.
In this work, the authors have not yet applied their method to the full Daphnia model but state that they plan to do this in a future work. They say they also plan to explore the error and convergence analysis when it comes to computing periodic solutions. In this realm, and if their method can be implemented in various popular software packages, it may provide a welcomed improvement to numerical simulation of differential equations. Personally, I would also like to see some example comparisons in terms of actual time taken per computation with various methods, to gain an intuitive understanding of the efficiency gain with a practical implementation.
For the entire collection see [Zbl 1445.92001].
MSC:
37N25 Dynamical systems in biology
92D25 Population dynamics (general)
92B05 General biology and biomathematics
34K13 Periodic solutions to functional-differential equations
65L60 Finite element, Rayleigh-Ritz, Galerkin and collocation methods for ordinary differential equations
PDF BibTeX XML Cite
Full Text: DOI
References:
[1] Allgower, E.L., Georg, K.: Introduction to numerical continuation methods. Classics in Applied Mathematics, vol. 45. SIAM, Philadelphia (2003) · Zbl 1036.65047
[2] Ascher, U., Mattheij, R.M.M., Russell, R.D.: Numerical Solution of Boundary Value Problems for Ordinary Differential Equations. Prentice-Hall, New Jersey (1988) · Zbl 0671.65063
[3] AUTO. http://indy.cs.concordia.ca/auto/ · Zbl 0786.68067
[4] Bader, G.: Numerische behandlung von randwertproblemen für funktionaldifferentialgleichungen. Technical Report 227, Universität Heidelberg (1983)
[5] Bader, G.: Solving boundary value problems for functional-differential equations by collocation. In: Ascher, U.M., Russell, R.D. (eds.) Numerical Boundary Value ODEs (Vancouver, B.C., 1984). Progress in Scientific Computing, vol. 5, pp. 227-243. Birkhäuser, Basel (1985)
[6] Bellen, A.: One-step collocation for delay differential equations. J. Comput. Appl. Math. 10, 275-283 (1984) · Zbl 0538.65047
[7] Bellen, A., Zennaro, M.: Numerical methods for delay differential equations. Numerical Mathematics and Scientific Computing series. Oxford University Press, Oxford (2003) · Zbl 1038.65058
[8] Breda, D., Liessi, D.: Approximation of eigenvalues of evolution operators for linear renewal equations. SIAM J. Numer. Anal. 56(3), 1456-1481 (2018) · Zbl 1404.65319
[9] Breda, D., Diekmann, O., de Graaf, W., Pugliese, A., Vermiglio, R.: On the formulation of epidemic models (an appraisal of Kermack and McKendrick). J. Biol. Dyn. 6(2), 103-117 (2012)
[10] Breda, D., Maset, S., Vermiglio, R.: Approximation of eigenvalues of evolution operators for linear retarded functional differential equations. SIAM J. Numer. Anal. 50(3), 1456-1483 (2012) · Zbl 1252.34089
[11] Breda, D., Diekmann, O., Maset, S., Vermiglio, R.: A numerical approach for investigating the stability of equilibria for structured population models. J. Biol. Dyn. 7(1), 4-20 (2013)
[12] Breda, D., Getto, P., Sánchez Sanz, J., Vermiglio, R.: Computing the eigenvalues of realistic Daphnia models by pseudospectral methods. SIAM J. Sci. Comput. 37(6), 2607-2629 (2015) · Zbl 1335.92072
[13] Breda, D., Diekmann, O., Gyllenberg, M., Scarabel, F., Vermiglio, R.: Pseudospectral discretization of nonlinear delay equations: new prospects for numerical bifurcation analysis. SIAM J. Appl. Dyn. Syst. 15(1), 1-23 (2016) · Zbl 1352.34101
[14] Breda, D., Diekmann, O., Liessi, D., Scarabel, F.: Numerical bifurcation analysis of a class of nonlinear renewal equations. Electron. J. Qual. Theory Differ. Equ. 65, 1-24 (2016) · Zbl 1413.65298
[15] Brunner, H.: Collocation methods for Volterra integral and related functional differential equations. No. 15 in Cambridge Monographs on Applied and Computational Mathematics. Cambridge University Press, Cambridge (2004)
[16] Canuto, C., Hussaini, M.Y., Quarteroni, A., Zang, T.: Spectral methods. Evolution to Complex Geometries and Applications to Fluid Dynamics. Scientific Computation Series. Springer, Berlin (2007) · Zbl 1121.76001
[17] de Roos, A.M., Diekmann, O., Getto, P., Kirkilionis, M.A.: Numerical equilibrium analysis for structured consumer resource models. Bull. Math. Biol. 72, 259-297 (2010) · Zbl 1185.92088
[18] Dhooge, A., Govaerts, W.J.F., Kuznetsov, Y.A., Meijer, H.G.E., Sautois, B.: New features of the software MatCont for bifurcation analysis of dynamical systems. Math. Comput. Model. Dyn. Syst. 14, 147-175 (2008) · Zbl 1158.34302
[19] Diekmann, O., Gyllenberg, M., Metz, J.A.J., Thieme, H.R.: On the formulation and analysis of general deterministic structured population models. I. linear theory. J. Math. Biol. 36, 349-388 (1998) · Zbl 0909.92023
[20] Diekmann, O., Gyllenberg, M., Huang, H., Kirkilionis, M., Metz, J.A.J., Thieme, H.R.: On the formulation and analysis of general deterministic structured population models. II. nonlinear theory. J. Math. Biol. 43, 157-189 (2001) · Zbl 1028.92019
[21] Diekmann, O., Gyllenberg, M., Metz, J.A.J., Nakaoka, S., de Roos, A.M.: Daphnia revisited: local stability and bifurcation theory for physiologically structured population models explained by way of an example. J. Math. Biol. 61(2), 277-318 (2010) · Zbl 1208.92082
[22] Diekmann, O., Getto, P., Nakata, Y.: On the characteristic equation λ = α_1 + (α_2 + α_3λ)e^−λ and its use in the context of a cell population model. J. Math. Biol. 72, 877-908 (2016) · Zbl 1344.34077
[23] Doedel, E.: Lecture notes on numerical analysis of nonlinear equations. In: Osinga, H.M., Krauskopf, B., Galán-Vioque, J. (eds.) Numerical Continuation Methods for Dynamical Systems, Understanding Complex Systems, pp. 1-49. Springer, Berlin (2007) · Zbl 1130.65119
[24] Engelborghs, K., Luzyanina, T., in ’t Hout, K.J., Roose, D.: Collocation methods for the computation of periodic solutions of delay differential equations. SIAM J. Sci. Comput. 22(5), 1593-1609 (2001) · Zbl 0981.65082
[25] Ermentrout, B.: Simulating, Analyzing, and Animating Dynamical Systems - A Guide to XPPAUT for Researchers and Students. Software - Environment - Tools series. SIAM, Philadelphia (2002) · Zbl 1003.68738
[26] Govaerts, W.J.F.: Numerical Methods for Bifurcations of Dynamical Equilibria. SIAM, Philadelphia (2000) · Zbl 0935.37054
[27] Hairer, E., Nörsett, S.P., Wanner, G.: Solving ordinary differential equations I: nonstiff problems. No. 8 in Computational Mathematics. Springer, Berlin (1987) · Zbl 0638.65058
[28] Inaba, H.: Age-Structured Population Dynamics in Demography and Epidemiology. Springer, New York (2017) · Zbl 1370.92010
[29] Jagers, P.: The deterministic evolution of general branching populations. In: de Gunst, M., Klaassen, C., van der Vaart, A. (eds.) State of the Art in Probability and Statistics. Lecture Notes-Monograph Series, vol. 36. Institute of Mathematical Statistics, Beachwood (2001) · Zbl 1371.60149
[30] Liessi, D.: Pseudospectral methods for stability of periodic solutions of delay models. Ph.D. Thesis, PhD in Computer Science, Mathematics and Physics, Università di Udine (2018) · Zbl 1404.65319
[31] Maset, S.: The collocation method in the numerical solution of boundary value problems for neutral functional differential equations. Part I: Convergence results. SIAM J. Numer. Anal. 53(6), 2771-2793 (2015) · Zbl 1330.65122
[32] Maset, S.: The collocation method in the numerical solution of boundary value problems for neutral functional differential equations. Part II: Differential equations with deviating arguments. SIAM J. Numer. Anal. 53(6), 2794-2821 (2015) · Zbl 1330.65123
[33] Maset, S.: An abstract framework in the numerical solution of boundary value problems for neutral functional differential equations. Numer. Math. 133(3), 525-555 (2016) · Zbl 1348.65104
[34] MatCont. https://sourceforge.net/projects/matcont/
[35] Metz, H., Diekmann, O.: The dynamics of physiologically structured populations. No. 68 in Lecture Notes in Biomathematics. Springer, New York (1986)
[36] Sánchez Sanz, J., Getto, P.: Numerical bifurcation analysis of physiologically structured populations: Consumer-resource, cannibalistic and trophic models. Bull. Math. Biol. 78(7), 1546-84 (2016) · Zbl 1367.92108
[37] Shurenkov, V.M.: On the theory of Markov renewal. Theory Probab. Appl. 29, 247-265 (1984) · Zbl 0557.60078
[38] Trefethen, L.N.: Spectral methods in MATLAB. Software - Environment - Tools series. SIAM, Philadelphia (2000)
[39] Weideman, J.A., Reddy, S.C.: A MATLAB differentiation matrix suite. ACM Trans. Math. Softw. 26(4), 465-519 (2000)
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.