Skip to main content
Log in

Nonlinear Physiologically Structured Population Models with Two Internal Variables

  • Published:
Journal of Nonlinear Science Aims and scope Submit manuscript

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.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

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)

    MathSciNet  MATH  Google Scholar 

  • 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)

    MathSciNet  MATH  Google Scholar 

  • 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)

    Google Scholar 

  • 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)

    MathSciNet  MATH  Google Scholar 

  • 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)

    MathSciNet  MATH  Google Scholar 

  • Calsina, A., Saldana, J.: A model of physiologically structured population dynamics with a nonlinear individual growth rate. J. Math. Biol. 33, 335–364 (1995)

    MathSciNet  MATH  Google Scholar 

  • 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)

    MathSciNet  MATH  Google Scholar 

  • Chu, J., Magal, P.: Hopf bifurcation for a size-structured model with resting phase. Discrete Contin. Dyn. Syst. 33, 4891–4921 (2013)

    MathSciNet  MATH  Google Scholar 

  • Cushing, J.M.: Equilibria in structured populations. J. Math. Biol. 23, 15–39 (1985)

    MathSciNet  MATH  Google Scholar 

  • Cushing, J.M.: Equilibria in systems of interacting structured populations. J. Math. Biol. 24, 627–649 (1987)

    MathSciNet  MATH  Google Scholar 

  • Cushing, J.M.: A competition model for size-structured species. SIAM. J. Appl. Math. 49, 838–858 (1989)

    MathSciNet  MATH  Google Scholar 

  • Da Prato, G., Sinestrari, E.: Differential operators with non dense domain. Ann. Sc. Norm. Sup. Pisa Cl. Sci. 14(2), 285–344 (1987)

    MATH  Google Scholar 

  • Doumic, M.: Analysis of a population model structured by the cells molecular content. Math. Model. Nat. Phenom. 2(3), 121–152 (2007)

    MathSciNet  MATH  Google Scholar 

  • 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)

    MathSciNet  MATH  Google Scholar 

  • 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)

    MathSciNet  MATH  Google Scholar 

  • Farkas, J.Z., Hagen, T.: Stability and regularity results for a size-structured population model. J. Math. Anal. Appl. 328, 119–136 (2007)

    MathSciNet  MATH  Google Scholar 

  • Farkas, J.Z., Green, D.M., Hinow, P.: Semigroup analysis of structured parasite populations. Math. Model. Nat. Phenom. 5(3), 94–114 (2010)

    MathSciNet  MATH  Google Scholar 

  • Gripenberg, G.: On a nonlinear integral equation modelling an epidemic in an age-structured population. J. Reine Angew. Math. 341, 54–67 (1893)

    MathSciNet  MATH  Google Scholar 

  • 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)

    MathSciNet  MATH  Google Scholar 

  • Gyllenberg, M., Webb, G.F.: Age-size structure in populations with quiescence. Math. Biol. 86, 67–95 (1987)

    MathSciNet  MATH  Google Scholar 

  • 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)

    Google Scholar 

  • Hethcote, H.W.: Optimal ages of vaccination for measles. Math. Biosci. 89(1), 29–52 (1988)

    MathSciNet  MATH  Google Scholar 

  • Hethcote, H.W.: An age-structured model for pertussis transmission. Math. Biosci. 145(2), 89–136 (1997)

    MathSciNet  MATH  Google Scholar 

  • Hethcote, H.W.: Simulations of pertussis epidemiology in the United States: effects of adult booster vaccinations. Math. Biosci. 158(1), 47–73 (1999)

    MathSciNet  MATH  Google Scholar 

  • Hoppensteadt, F.: An age dependent epidemic model. J. Frankl. Inst. 297, 325–333 (1974)

    MATH  Google Scholar 

  • Inaba, H.: Threshold and stability results for age-structured epidemic model. J. Math. Biol. 28, 411–434 (1990)

    MathSciNet  MATH  Google Scholar 

  • Inaba, H.: Endemic threshold results in an age-duration-structured population model for HIV infection. Math. Biosci. 201(1–2), 15–47 (2006)

    MathSciNet  MATH  Google Scholar 

  • Inaba, H.: On pandemic threshold theorem of the early Kermack–McKendrick model with individual heterogeneity. Math. Popul. Stud. 21, 95–111 (2014)

    MathSciNet  MATH  Google Scholar 

  • Inaba, H.: Endemic threshold analysis for the Kermack–McKendrick reinfection model. Josai Math. Monogr. 9, 105–133 (2016)

    Google Scholar 

  • Inaba, H.: Age-Structured Population Dynamics in Demography and Epidemiology. Springer, New York (2017)

    MATH  Google Scholar 

  • 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)

    MathSciNet  MATH  Google Scholar 

  • 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)

    Google Scholar 

  • Krasnoselskii, M.A.: Positive Solutions of Operator Equations. Noordhoff, Groningen (1964)

    Google Scholar 

  • Laroche, B., Perasso, A.: Threshold behaviour of a SI epidemiological model with two structuring variables. J. Evol. Equ. 16, 293–315 (2016)

    MathSciNet  MATH  Google Scholar 

  • Magal, P., Ruan, S.: Theory and Applications of Abstract Semilinear Cauchy Problems. Springer, New York (2018)

    MATH  Google Scholar 

  • Marek, I.: Frobenius theory of positive operators: comparison theorems and applications. SIAM J. Appl. Math. 19(3), 607–628 (1970)

    MathSciNet  MATH  Google Scholar 

  • Matucci, S.: Existence, uniqueness and asymptotic behavior for a multi stage evolution problem. Math. Models Methods Appl. Sci. 8, 1013–1041 (1995)

    MathSciNet  MATH  Google Scholar 

  • 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)

    MathSciNet  MATH  Google Scholar 

  • Sawashima, I.: On spectral properties of some positive operators. Nat. Sci. Rep. Ochanomizu Univ. 15(2), 53–64 (1964)

    MathSciNet  MATH  Google Scholar 

  • Sinko, J.W., Streifer, W.: A new model for age-size structure of a population. Ecology 48, 910–918 (1967)

    Google Scholar 

  • Stadler, E.: Eigensolutions and spectral analysis of a model for vertical gene transfer of plasmids. J. Math. Biol. 78(5), 1299–1330 (2019)

    MathSciNet  MATH  Google Scholar 

  • 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)

    MathSciNet  MATH  Google Scholar 

  • Thieme, H.R.: Semiflows generated by Lipschitz perturbations of non-densely defined operators. Differ. Integral Equ. 3(6), 1035–1066 (1990)

    MathSciNet  MATH  Google Scholar 

  • 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)

    MathSciNet  MATH  Google Scholar 

  • Webb, G.F.: Theory of Nonlinear Age-Dependent Population Dynamics. Marcel Dekker, New York (1984)

    Google Scholar 

  • Webb, G.F.: Dynamics of populations structured by internal variables. Math. Z. 189, 319–335 (1985)

    MathSciNet  MATH  Google Scholar 

  • 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)

    Google Scholar 

Download references

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

Authors

Corresponding author

Correspondence to Shigui Ruan.

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:

  1. (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;

  2. (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;

  3. (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

  1. (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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s00332-020-09638-5

Keywords

Mathematics Subject Classification

Navigation