Abstract
First-order hyperbolic partial differential equations with two internal variables have been used to model biological and epidemiological problems with two physiological structures, such as chronological age and infection age in epidemic models, age and another physiological character (maturation, size, stage) in population models, and cell-age and molecular content (cyclin content, maturity level, plasmid copies, telomere length) in cell population models. In this paper, we study nonlinear double physiologically structured population models with two internal variables by applying integrated semigroup theory and non-densely defined operators. We consider first a semilinear model and then a nonlinear model, use the method of characteristic lines to find the resolvent of the infinitesimal generator and the variation of constant formula, apply Krasnoselskii’s fixed point theorem to obtain the existence of a steady state, and study the stability of the steady state by estimating the essential growth bound of the semigroup. Finally, we generalize the techniques to investigate a nonlinear age-size structured model with size-dependent growth rate.
Similar content being viewed by others
References
Bekkal Brikci, F., Clairambault, J., Ribba, B., Perthame, B.: An age-and-cyclin-structured cell population model for healthy and tumoral tissues. J. Math. Biol. 57, 91–110 (2008)
Bénilan, P., Crandall, M.G., Pazy, A.: “Bonnes solutions” d’un problème d’évolution semi-linéaire. C. R. Acad. Sci. Paris Sér. I Math. 306, 527–530 (1988)
Bernard, S., Pujo-Menjouet, L., Mackey, M.C.: Analysis of cell kinetics using a cell division marker: mathematical modeling of experimental data. Biophys. J. 84(5), 3414–3424 (2003)
Burie, J.-B., Ducrot, A., Mbengue, A.A.: Asymptotic behavior of an age and infection age structured model for the propagation of fungal diseases in plants. Discrete Contin. Dyn. Syst. Ser. B 22(7), 2879–2905 (2017)
Busenberg, S.N., Iannelli, M., Thieme, H.R.: Global behavior of an age-structured epidemic model. SIAM J. Math. Anal. 22(4), 1065–1080 (1991)
Calsina, A., Saldana, J.: A model of physiologically structured population dynamics with a nonlinear individual growth rate. J. Math. Biol. 33, 335–364 (1995)
Chu, J., Ducrot, A., Magal, P., Ruan, S.: Hopf bifurcation in a size-structured population dynamic model with random growth. J. Differ. Equ. 247, 956–1000 (2009)
Chu, J., Magal, P.: Hopf bifurcation for a size-structured model with resting phase. Discrete Contin. Dyn. Syst. 33, 4891–4921 (2013)
Cushing, J.M.: Equilibria in structured populations. J. Math. Biol. 23, 15–39 (1985)
Cushing, J.M.: Equilibria in systems of interacting structured populations. J. Math. Biol. 24, 627–649 (1987)
Cushing, J.M.: A competition model for size-structured species. SIAM. J. Appl. Math. 49, 838–858 (1989)
Da Prato, G., Sinestrari, E.: Differential operators with non dense domain. Ann. Sc. Norm. Sup. Pisa Cl. Sci. 14(2), 285–344 (1987)
Doumic, M.: Analysis of a population model structured by the cells molecular content. Math. Model. Nat. Phenom. 2(3), 121–152 (2007)
Dyson, J., Villella-Bressan, R., Webb, G.: A nonlinear age and maturity structured model of population dynamics: I. Basic theory. J. Math. Anal. Appl. 242(1), 93–104 (2000)
Dyson, J., Villella-Bressan, R., Webb, G.: A nonlinear age and maturity structured model of population dynamics: II. Chaos. J. Math. Anal. Appl. 242(2), 255–270 (2000)
Farkas, J.Z., Hagen, T.: Stability and regularity results for a size-structured population model. J. Math. Anal. Appl. 328, 119–136 (2007)
Farkas, J.Z., Green, D.M., Hinow, P.: Semigroup analysis of structured parasite populations. Math. Model. Nat. Phenom. 5(3), 94–114 (2010)
Gripenberg, G.: On a nonlinear integral equation modelling an epidemic in an age-structured population. J. Reine Angew. Math. 341, 54–67 (1893)
Gwiazda, P., Lorenz, T., Marciniak-Czochra, A.: A nonlinear structured population model: Lipschitz continuity of measure-valued solutions with respect to model ingredients. J. Differ. Equ. 248, 2703–2735 (2010)
Gyllenberg, M., Webb, G.F.: Age-size structure in populations with quiescence. Math. Biol. 86, 67–95 (1987)
Heijmans, H.J.A.M.: The dynamical behaviour of the age-size-distribution of a cell population. In: Metz, J.A.J., Diekmann, O. (eds.) The Dynamics of Physiologically Structured Populations. Lecture Notes in Biomathematics, vol. 68, pp. 185–202. Springer, Berlin (1986)
Hethcote, H.W.: Optimal ages of vaccination for measles. Math. Biosci. 89(1), 29–52 (1988)
Hethcote, H.W.: An age-structured model for pertussis transmission. Math. Biosci. 145(2), 89–136 (1997)
Hethcote, H.W.: Simulations of pertussis epidemiology in the United States: effects of adult booster vaccinations. Math. Biosci. 158(1), 47–73 (1999)
Hoppensteadt, F.: An age dependent epidemic model. J. Frankl. Inst. 297, 325–333 (1974)
Inaba, H.: Threshold and stability results for age-structured epidemic model. J. Math. Biol. 28, 411–434 (1990)
Inaba, H.: Endemic threshold results in an age-duration-structured population model for HIV infection. Math. Biosci. 201(1–2), 15–47 (2006)
Inaba, H.: On pandemic threshold theorem of the early Kermack–McKendrick model with individual heterogeneity. Math. Popul. Stud. 21, 95–111 (2014)
Inaba, H.: Endemic threshold analysis for the Kermack–McKendrick reinfection model. Josai Math. Monogr. 9, 105–133 (2016)
Inaba, H.: Age-Structured Population Dynamics in Demography and Epidemiology. Springer, New York (2017)
Kang, H., Huo, X., Ruan, S.: On first-order hyperbolic partial differential equations with two internal variables modeling population dynamics of two physiological structures. Ann. Mat. Pura Appl. (2020). https://doi.org/10.1007/s10231-020-01001-5
Kapitanov, G.: A mathematical model of cancer stem cell lineage population dynamics with mutation accumulation and telomere length hierarchies. Math. Model. Nat. Phenom. 7(1), 136–165 (2012)
Koijman, S.A.L.M., Metz, J.A.J.: On the dynamics of chemically stressed populations: the deduction of population consequences from effects on individuals. Ecotoxicol. Environ. Saf. 8, 254–274 (1984)
Krasnoselskii, M.A.: Positive Solutions of Operator Equations. Noordhoff, Groningen (1964)
Laroche, B., Perasso, A.: Threshold behaviour of a SI epidemiological model with two structuring variables. J. Evol. Equ. 16, 293–315 (2016)
Magal, P., Ruan, S.: Theory and Applications of Abstract Semilinear Cauchy Problems. Springer, New York (2018)
Marek, I.: Frobenius theory of positive operators: comparison theorems and applications. SIAM J. Appl. Math. 19(3), 607–628 (1970)
Matucci, S.: Existence, uniqueness and asymptotic behavior for a multi stage evolution problem. Math. Models Methods Appl. Sci. 8, 1013–1041 (1995)
McNair, J.N., Goulden, C.E.: The dynamics of age-structured population with a gestation period: density-independent growth and egg ratio methods for estimating the birth-rate. Theor. Popul. Biol. 39, 1–29 (1991)
Sawashima, I.: On spectral properties of some positive operators. Nat. Sci. Rep. Ochanomizu Univ. 15(2), 53–64 (1964)
Sinko, J.W., Streifer, W.: A new model for age-size structure of a population. Ecology 48, 910–918 (1967)
Stadler, E.: Eigensolutions and spectral analysis of a model for vertical gene transfer of plasmids. J. Math. Biol. 78(5), 1299–1330 (2019)
Thieme, H.R.: Well-posedness of physiologically structured population models for Daphnia magna—how biological concepts can benefit by abstract mathematical analysis. J. Math. Biol. 26(3), 299–317 (1988)
Thieme, H.R.: Semiflows generated by Lipschitz perturbations of non-densely defined operators. Differ. Integral Equ. 3(6), 1035–1066 (1990)
Tucker, S.L., Zimmerman, S.O.: A nonlinear model of population dynamics containing an arbitrary number of continuous structure variables. SIAM. J. Appl. Math. 48, 549–591 (1988)
Webb, G.F.: Theory of Nonlinear Age-Dependent Population Dynamics. Marcel Dekker, New York (1984)
Webb, G.F.: Dynamics of populations structured by internal variables. Math. Z. 189, 319–335 (1985)
Webb, G.F.: Population models structured by age, size, and spatial position. In: Magal, P., Ruan, S. (eds.) Structured Population Models in Biology and Epidemiology. Lecturer Notes in Mathematics, vol. 1936, pp. 1–49. Springer, Berlin (2008)
Acknowledgements
We would like to thank the three anonymous reviewers for their helpful comments and suggestions which helped us to improve the presentation of the paper.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Anthony Bloch.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Research was partially supported by National Science Foundation (DMS-1853622).
A Appendix: Positive Operators
A Appendix: Positive Operators
In this “Appendix”, we recall some definitions and results of positive operator theory on ordered Banach spaces from Inaba (2006). For more complete exposition, we refer to Heijmans (1986); Marek (1970), and Sawashima (1964).
Let E be a real or complex Banach space and \(E^*\) be its dual (the space of all linear functionals on E). Write the value of \(f\in E^*\) at \(\psi \in E\) as \(\langle f, \psi \rangle \). A nonempty closed subset \(E_+\) is called a cone if the following hold: (i) \(E_++E_+\subset E_+\), (ii) \(\lambda E_+\subset E_+\) for \(\lambda \ge 0\), (iii) \(E_+\cap (-E_+)=\{0\}\). Define the order in E such that \(x\le y\) if and only if \(y-x\in E_+\) and \(x<y\) if and only if \(y-x\in E_+{\setminus }\{0\}\). The cone \(E_+\) is called total if the set \(\{\psi -\phi : \psi , \phi \in E_+\}\) is dense in E. The dual cone \(E^*_+\) is the subset of \(E^*\) consisting of all positive linear functionals on E; that is, \(f\in E^*_+\) if and only if \(\langle f, \psi \rangle \ge 0\) for all \(\psi \in E_+\). \(\psi \in E_+\) is called a quasi-interior point if \(\langle f, \psi \rangle >0\) for all \(f\in E^*_+{\setminus }\{0\}\). \(f\in E^*_+\) is said to be strictly positive if \(\langle f, \psi \rangle >0\) for all \(\psi \in E_+{\setminus }\{0\}\). The cone \(E_+\) is called generating if \(E=E_+-E_+\) and is called normal if \(E^*=E^*_+-E^*_+\).
An ordered Banach space \((E, \le )\) is called a Banach lattice if (i) any two elements \(x, y \in E\) have a supremum \(x \vee y=\sup \{x, y\}\) and an infimum \(x \wedge y=\inf \{x, y\}\) in E; and (ii) \(|x| \le |y|\) implies \(\Vert x\Vert \le \Vert y\Vert \) for \(x, y \in E,\) where the modulus of x is defined by \(|x|=x \vee (-x).\)
Let B(E) be the set of bounded linear operators from E to E. \(T\in B(E)\) is said to be positive if \(T(E_+)\subset E_+\). For \(T, S\in B(E)\), we say \(T\ge S\) if \((T-S)(E_+)\subset E_+\). A positive operator \(T\in B(E)\) is called semi-nonsupporting if for every pair \(\psi \in E_+{\setminus }\{0\}, f\in E^*_+{\setminus }\{0\}\), there exists a positive integer \(p=p(\psi , f)\) such that \(\langle f, T^p\psi \rangle >0\). A positive operator \(T\in B(E)\) is called nonsupporting if for every pair \(\psi \in E_+{\setminus }\{0\}, f\in E^*_+{\setminus }\{0\}\), there exists a positive integer \(p=p(\psi , f)\) such that \(\langle f, T^n\psi \rangle >0\) for all \(n\ge p\). The spectral radius of \(T\in B(E)\) is denoted by r(T). \(\sigma (T)\) denotes the spectrum of T and \(\sigma _P(T)\) denotes the point spectrum of T. If there exists a nonzero \(x\in E\) which satisfies \(Tx=\lambda x\), \(\lambda \) is called a proper value and x a proper vector corresponding to \(\lambda \).
From results in Sawashima (1964) and Inaba (2006), we state the following propositions.
Proposition A.1
Let E be a Banach space, and let \(T\in B(E)\) be compact and semi-nonsupporting. Then, the following statements hold:
-
(i)
\(r(T)\in \sigma _P(T){\setminus }\{0\}\) and r(T) is a simple pole of the resolvent \(\lambda I-T;\) that is, r(T) is an algebraically simple eigenvalue of T;
-
(ii)
The eigenspace of T corresponding to r(T) is one-dimensional and the corresponding eigenvector \(\psi \in E_+\) is a quasi-interior point. The relation \(T\phi =\mu \phi \) with \(\phi \in E_+\) implies that \(\phi =c\psi \) for some constant c;
-
(iii)
The eigenspace of \(T^*\) corresponding to r(T) is also a one-dimensional subspace of \(E^*\) spanned by a strictly positive functional \(f\in E^*_+\).
Proposition A.2
Let E be a Banach space with positive cone \(E_+\) which is total. Let \(T\in B(E)\) be positive and have the resolvent \(\lambda I-T\) with the point r(T) as its pole. Then, T is a semi-nonsupporting operator if and only if \(r(T)>0\) and T satisfies (A), where
-
(A)
Every proper eigenvector corresponding to the proper eigenvalue r(T) lying in \(E_+\) is a quasi-interior point of \(E_+\) and every proper eigenvector corresponding to r(T) lying in \(E_+^*\) is strictly positive.
Rights and permissions
About this article
Cite this article
Kang, H., Huo, X. & Ruan, S. Nonlinear Physiologically Structured Population Models with Two Internal Variables. J Nonlinear Sci 30, 2847–2884 (2020). https://doi.org/10.1007/s00332-020-09638-5
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00332-020-09638-5
Keywords
- Physiological structure
- Cauchy problem with non-dense domain
- Integrated semigroups
- Infinitesimal generator
- Spectrum theory
- Stability