×

A numerical framework for computing steady states of structured population models and their stability. (English) Zbl 1364.35018

Summary: Structured population models are a class of general evolution equations which are widely used in the study of biological systems. Many theoretical methods are available for establishing existence and stability of steady states of general evolution equations. However, except for very special cases, finding an analytical form of stationary solutions for evolution equations is a challenging task. In the present paper, we develop a numerical framework for computing approximations to stationary solutions of general evolution equations, which can also be used to produce approximate existence and stability regions for steady states. In particular, we use the Trotter-Kato Theorem to approximate the infinitesimal generator of an evolution equation on a finite dimensional space, which in turn reduces the evolution equation into a system of ordinary differential equations. Consequently, we approximate and study the asymptotic behavior of stationary solutions. We illustrate the convergence of our numerical framework by applying it to a linear Sinko-Streifer structured population model for which the exact form of the steady state is known. To further illustrate the utility of our approach, we apply our framework to nonlinear population balance equation, which is an extension of well-known Smoluchowski coagulation-fragmentation model to biological populations. We also demonstrate that our numerical framework can be used to gain insight about the theoretical stability of the stationary solutions of the evolution equations. Furthermore, the open source Python program that we have developed for our numerical simulations is freely available from our GitHub repository (github.com/MathBioCU).

MSC:

35A35 Theoretical approximation in context of PDEs
35L50 Initial-boundary value problems for first-order hyperbolic systems
35B40 Asymptotic behavior of solutions to PDEs
92B05 General biology and biomathematics
65N40 Method of lines for boundary value problems involving PDEs
PDFBibTeX XMLCite
Full Text: DOI arXiv

References:

[1] A. S. Ackleh, Parameter estimation in a structured algal coagulation-fragmentation model,, Nonlinear Anal. Theory Methods Appl., 28, 837 (1997) · Zbl 0869.35025 · doi:10.1016/0362-546X(95)00195-2
[2] A. S. Ackleh, Modeling aggregation and growth processes in an algal population model: analysis and computations,, J. Math. Biol., 35, 480 (1997) · Zbl 0867.92024
[3] V. I. Arnold, <em>Ordinary Differential Equations</em>,, Second printing of the 1992 edition. Universitext. Springer-Verlag (1992)
[4] J. Banasiak, Blow-up of solutions to some coagulation and fragmentation equations with growth,, Dyn. Syst., 1, 126 (2011) · Zbl 1306.35008
[5] J. Banasiak, Coagulation, fragmentation and growth processes in a size structured population,, Discrete Contin. Dyn. Syst. - Ser. B, 11, 563 (2009) · Zbl 1178.35136
[6] H. T. Banks, Transformation semigroups and L1-approximation for size structured population models,, Semigroup Forum, 38, 141 (1989) · Zbl 0687.47033
[7] C. Biggs, Modelling activated sludge flocculation using population balances,, Powder Technol., 124, 201 (2002)
[8] D. M. Bortz, Chapter 17: Modeling and simulation for nanomaterials in fluids: Nanoparticle self-assembly., In Tewary, 419 (2015)
[9] D. M. Bortz, Klebsiella pneumoniae Flocculation Dynamics,, Bull. Math. Biol., 70, 745 (2008) · Zbl 1144.92022
[10] D. Breda, Solution operator approximations for characteristic roots of delay differential equations,, Applied Numerical Mathematics, 56, 305 (2006) · Zbl 1095.65072 · doi:10.1016/j.apnum.2005.04.010
[11] D. Breda, A numerical approach for investigating the stability of equilibria for structured population models,, J. Biol. Dyn., 7, 4 (2013) · Zbl 1448.92178
[12] D. Breda, Computing the characteristic roots for delay differential equations,, IMA J Numer Anal, 24, 1 (2004) · Zbl 1054.65079
[13] D. Breda, Pseudospectral Differencing Methods for Characteristic Roots of Delay Differential Equations,, SIAM J. Sci. Comput., 27, 482 (2005) · Zbl 1092.65054
[14] D. Breda, Pseudospectral approximation of eigenvalues of derivative operators with non-local boundary conditions,, Applied Numerical Mathematics, 56, 318 (2006) · Zbl 1099.65064
[15] D. Breda, Numerical approximation of characteristic values of partial retarded functional differential equations,, Numer. Math., 113, 181 (2009) · Zbl 1175.65115
[16] E. Byrne, Postfragmentation density function for bacterial aggregates in laminar flow,, Phys. Rev. E, 83 (2011)
[17] V. Calvez, Self-similarity in a general aggregation-fragmentation problem, Application to fitness analysis,, Journal de Mathématiques Pures et Appliquées, 98, 1 (2012) · Zbl 1259.35151
[18] V. Calvez, Prion dynamics with size dependency-strain phenomena,, J. Biol. Dyn., 4, 28 (2010) · Zbl 1315.92039
[19] A. M. De Roos, Demographic analysis of continuous-time life-history models,, Ecol. Lett., 11, 1 (2008)
[20] A. M. De Roos, PSPManalysis, 2014,, <a href=
[21] A. M. de Roos, Numerical equilibrium analysis for structured consumer resource models,, Bulletin of Mathematical Biology, 72, 259 (2010) · Zbl 1185.92088
[22] O. Diekmann, Steady-state analysis of structured population models,, Theoretical Population Biology, 63, 309 (2003) · Zbl 1098.92062
[23] C. A. Dorao, Application of the least-squares method for solving population balance problems in Rd+1,, Chemical Engineering Science, 61, 5070 · Zbl 1110.65124
[24] C. A. Dorao, A least squares method for the solution of population balance problems,, Computers & Chemical Engineering, 30, 535 · Zbl 1110.65124
[25] C. A. Dorao, Least-squares spectral method for solving advective population balance problems,, Journal of Computational and Applied Mathematics, 201, 247 (2007) · Zbl 1110.65124
[26] M. Doumic-Jauffret, Eigenelements of a general aggregation-fragmentation model,, <a href= (2009) · Zbl 1201.35086
[27] K. Engelborghs, Numerical Bifurcation Analysis of Delay Differential Equations Using DDE-BIFTOOL,, ACM Trans Math Softw, 28, 1 (2002) · Zbl 1070.65556
[28] J. Z. Farkas, Stability and regularity results for a size-structured population model,, J. Math. Anal. Appl., 328, 119 (2007) · Zbl 1114.35046
[29] J. Z. Farkas, Steady states in hierarchical structured populations with distributed states at birth,, Discrete Contin. Dyn. Syst. - Ser. B, 17, 2671 (2012) · Zbl 1395.92123
[30] S. A. Gourley, Spatiotemporal Patterns of Disease Spread: Interaction of Physiological Structure, Spatial Movements, Disease Progression and Human Intervention,, In Magal, 165 (2008)
[31] M. J. Hounslow, A discretized population balance for continuous systems at steady state,, AIChE J., 36, 106 (1990)
[32] K. Ito, The Trotter-Kato theorem and approximation of PDEs,, Math Comp, 67, 21 (1998) · Zbl 0893.47025
[33] F. Kappel, Spline Approximations for Neutral Functional Differential Equations,, SIAM J. Numer. Anal., 18, 1058 (1981) · Zbl 0511.65053
[34] N. Kato, A Principle of Linearized Stability for Nonlinear Evolution Equations,, Transactions of the American Mathematical Society, 347, 2851 (1995) · Zbl 0843.47037 · doi:10.1090/S0002-9947-1995-1290722-8
[35] T. Kato, Remarks on pseudo-resolvents and infinitesimal generators of semi-groups,, Proc. Japan Acad., 35, 467 (1959) · Zbl 0095.10502 · doi:10.3792/pja/1195524254
[36] T. Kato, <em>Perturbation Theory for Linear Operators,</em>, Classics in Mathematics. Springer Berlin Heidelberg (1976) · Zbl 0342.47009
[37] M. A. Kirkilionis, Numerical continuation of equilibria of physiologically structured population models I: Theory,, Math. Models Methods Appl. Sci., 11, 1101 (2001) · Zbl 1013.92036
[38] P. Laurencot, Steady states for a coagulation-fragmentation equation with volume scattering,, SIAM Journal on Mathematical Analysis, 37, 531 (2005) · Zbl 1105.45007
[39] J. Makino, On the mass distribution of planetesimals in the early runaway stage,, New Astronomy, 3, 411 (1998)
[40] S. A. Matveev, A fast numerical method for the Cauchy problem for the Smoluchowski equation,, Journal of Computational Physics, 282, 23 (2015) · Zbl 1352.65250
[41] G. Menon, Dynamical scaling in Smoluchowski’s coagulation equations: Uniform convergence,, SIAM Rev., 48, 745 (2006) · Zbl 1117.70018
[42] I. Mirzaev, Steady state approximation, 2015., <a href=
[43] I. Mirzaev, Criteria for linearized stability for a size-structured population model,, <a href= (2015)
[44] I. Mirzaev, Stability of steady states for a class of flocculation equations with growth and removal,, <a href= (2015)
[45] M. Nicmanis, Finite-element methods for steady-state population balance equations,, AIChE J., 44, 2258 (1998)
[46] M. Nicmanis, Error estimation and control for the steady state population balance equation: 1. An a posteriori error estimate,, Chemical Engineering Science, 57, 2253 (2002)
[47] H.-S. Niwa, School size statistics of fish,, J. Theor. Biol., 195, 351 (1998)
[48] M. Powell, A hybrid method for nonlinear equations,, In Rabinowitz, 87 (1970) · Zbl 0277.65028
[49] H. R. Pruppacher, <em>Microphysics of Clouds and Precipitation: Reprinted 1980,</em>, Springer Science & Business Media (2012)
[50] D. Ramkrishna, <em>Population Balances: Theory and Applications to Particulate Systems in Engineering,</em>, Academic Press (2000)
[51] S. J. Schreiber, Invasion speeds for structured populations in fluctuating environments,, Theor. Ecol., 4, 423 (2011)
[52] J. W. Sinko, A New Model For Age-Size Structure of a Population,, Ecology, 48, 910 (1967)
[53] H. F. Trotter, Approximation of semi-groups of operators,, Pacific J. Math., 8, 887 (1958) · Zbl 0099.10302 · doi:10.2140/pjm.1958.8.887
[54] J. A. Wattis, An introduction to mathematical models of coagulation fragmentation processes: A discrete deterministic mean-field approach,, Phys. Nonlinear Phenom., 222, 1 (2006) · Zbl 1113.35145 · doi:10.1016/j.physd.2006.07.024
[55] G. F. Webb, <em>Theory of Nonlinear Age-Dependent Population Dynamics</em>,, CRC Press (1985) · Zbl 0555.92014
[56] R. M. Ziff, Kinetics of polymer gelation,, J. Chem. Phys., 73, 3492 (1980)
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.