# American Institute of Mathematical Sciences

August  2017, 14(4): 933-952. doi: 10.3934/mbe.2017049

## A numerical framework for computing steady states of structured population models and their stability

 Department of Applied Mathematics, University of Colorado, Boulder, CO, 80309-0526, United States

Received  January 2016 Accepted  January 2016 Published  February 2017

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 \emph{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).

Citation: Inom Mirzaev, David M. Bortz. A numerical framework for computing steady states of structured population models and their stability. Mathematical Biosciences & Engineering, 2017, 14 (4) : 933-952. doi: 10.3934/mbe.2017049
##### References:

show all references

##### References:
Results-of-the Results of the numerical simulations. a) $a,\,b\,,c$ values satisfying the necessary condition (12), form a 3D surface (blue surface). Steady states of the Sinko-Streifer model only exist on the red line b) Comparison of exact stationary solution (for the point marked with red star in Figure 1a) with approximate stationary solution for $n=100$. c) Absolute error between exact stationary solution and approximate stationary solution decays linearly as the dimension of approximate subspaces $\mathcal{X}_{n}$ increase.
Existence and stability regions for the steady states of the PBE a) Existence region for the steady states of the PBE forms a wedge like shape. b) Stability region for $b=0.1$, $a\in[0,\,15]$ and $c\in[0,\,5]$. c) Stability region for $b=0.5$, $a\in[0,\,15]$ and $c\in[0,\,5]$. d) Stability region for $b=1.0$, $a\in[0,\,15]$ and $c\in[0,\,5]$. Color bar represents the real part of rightmost eigenvalue of the Jacobian matrix evaluated at each steady state. Yellow regions represents the region for which a positive steady state does not exists.
An example steady-state solution of the PBE for $b=0.5,\,a=c=1$. b) Steady states for increasing renewal rate and $b=c=1$
Time evolution of the flocculation model with arbitrary initial conditions. a) Four different initial conditions are chosen close to the steady state. b) Solution of the PBE for those initial conditions at $t=10$. c) Evolution of the total number $M_{0}(t)$ of the flocs for $t\in[0,\,10]$. d) Evolution of the total mass $M_{1}(t)$ of the flocs for $t\in[0,\,10]$.
Change in zeroth and first moments with increasing dimension of the approximate space $\mathcal{X}_{n}$. a) Change in the total number and the total mass of the flocs with respect to increasing dimension $n$. Dashed red lines and dotted green lines corresponds to the total number and the total mass of the flocs of the steady state for $n=1000$, respectively. b) Steady state solution for $n=100$ and $n=500$.
Eigenvalues of the Jacobian $J_{\mathcal{F}}(\alpha)$ multiplied by $\Delta x$ for the steady state illustrated in Figure 3a. a) Eigenvalues of the Jacobian plotted in the complex plane for $n=20$. b) Eigenvalues of the Jacobian plotted in the complex plane for $n=50$. c) Eigenvalues of the Jacobian plotted in the complex plane for $n=200$. d) Change in the rightmost eigenvalue for increasing $n$. Dashed black line corresponds to the rightmost eigenvalue of the Jacobian for $n=1000$.
 [1] Xianlong Fu, Dongmei Zhu. Stability analysis for a size-structured juvenile-adult population model. Discrete & Continuous Dynamical Systems - B, 2014, 19 (2) : 391-417. doi: 10.3934/dcdsb.2014.19.391 [2] L. M. Abia, O. Angulo, J.C. López-Marcos. Size-structured population dynamics models and their numerical solutions. Discrete & Continuous Dynamical Systems - B, 2004, 4 (4) : 1203-1222. doi: 10.3934/dcdsb.2004.4.1203 [3] Xianlong Fu, Dongmei Zhu. Stability results for a size-structured population model with delayed birth process. Discrete & Continuous Dynamical Systems - B, 2013, 18 (1) : 109-131. doi: 10.3934/dcdsb.2013.18.109 [4] Dongxue Yan, Xianlong Fu. Asymptotic analysis of a spatially and size-structured population model with delayed birth process. Communications on Pure & Applied Analysis, 2016, 15 (2) : 637-655. doi: 10.3934/cpaa.2016.15.637 [5] Dongxue Yan, Yu Cao, Xianlong Fu. Asymptotic analysis of a size-structured cannibalism population model with delayed birth process. Discrete & Continuous Dynamical Systems - B, 2016, 21 (6) : 1975-1998. doi: 10.3934/dcdsb.2016032 [6] Dongxue Yan, Xianlong Fu. Asymptotic behavior of a hierarchical size-structured population model. Evolution Equations & Control Theory, 2018, 7 (2) : 293-316. doi: 10.3934/eect.2018015 [7] Keng Deng, Yixiang Wu. Extinction and uniform strong persistence of a size-structured population model. Discrete & Continuous Dynamical Systems - B, 2017, 22 (3) : 831-840. doi: 10.3934/dcdsb.2017041 [8] Azmy S. Ackleh, Vinodh K. Chellamuthu, Kazufumi Ito. Finite difference approximations for measure-valued solutions of a hierarchically size-structured population model. Mathematical Biosciences & Engineering, 2015, 12 (2) : 233-258. doi: 10.3934/mbe.2015.12.233 [9] Qihua Huang, Hao Wang. A toxin-mediated size-structured population model: Finite difference approximation and well-posedness. Mathematical Biosciences & Engineering, 2016, 13 (4) : 697-722. doi: 10.3934/mbe.2016015 [10] Blaise Faugeras, Olivier Maury. An advection-diffusion-reaction size-structured fish population dynamics model combined with a statistical parameter estimation procedure: Application to the Indian Ocean skipjack tuna fishery. Mathematical Biosciences & Engineering, 2005, 2 (4) : 719-741. doi: 10.3934/mbe.2005.2.719 [11] Rajesh Kumar, Jitendra Kumar, Gerald Warnecke. Convergence analysis of a finite volume scheme for solving non-linear aggregation-breakage population balance equations. Kinetic & Related Models, 2014, 7 (4) : 713-737. doi: 10.3934/krm.2014.7.713 [12] Jianquan Li, Zhien Ma. Stability analysis for SIS epidemic models with vaccination and constant population size. Discrete & Continuous Dynamical Systems - B, 2004, 4 (3) : 635-642. doi: 10.3934/dcdsb.2004.4.635 [13] Jacek Banasiak, Wilson Lamb. Coagulation, fragmentation and growth processes in a size structured population. Discrete & Continuous Dynamical Systems - B, 2009, 11 (3) : 563-585. doi: 10.3934/dcdsb.2009.11.563 [14] József Z. Farkas, Thomas Hagen. Asymptotic analysis of a size-structured cannibalism model with infinite dimensional environmental feedback. Communications on Pure & Applied Analysis, 2009, 8 (6) : 1825-1839. doi: 10.3934/cpaa.2009.8.1825 [15] Guangrui Li, Ming Mei, Yau Shu Wong. Nonlinear stability of traveling wavefronts in an age-structured reaction-diffusion population model. Mathematical Biosciences & Engineering, 2008, 5 (1) : 85-100. doi: 10.3934/mbe.2008.5.85 [16] Rong Liu, Feng-Qin Zhang, Yuming Chen. Optimal contraception control for a nonlinear population model with size structure and a separable mortality. Discrete & Continuous Dynamical Systems - B, 2016, 21 (10) : 3603-3618. doi: 10.3934/dcdsb.2016112 [17] Azmy S. Ackleh, Kazufumi Ito. An approximation scheme for a nonlinear size-dependent population model. Conference Publications, 1998, 1998 (Special) : 1-6. doi: 10.3934/proc.1998.1998.1 [18] Azmy S. Ackleh, H.T. Banks, Keng Deng, Shuhua Hu. Parameter Estimation in a Coupled System of Nonlinear Size-Structured Populations. Mathematical Biosciences & Engineering, 2005, 2 (2) : 289-315. doi: 10.3934/mbe.2005.2.289 [19] Diène Ngom, A. Iggidir, Aboudramane Guiro, Abderrahim Ouahbi. An observer for a nonlinear age-structured model of a harvested fish population. Mathematical Biosciences & Engineering, 2008, 5 (2) : 337-354. doi: 10.3934/mbe.2008.5.337 [20] H. L. Smith, X. Q. Zhao. Competitive exclusion in a discrete-time, size-structured chemostat model. Discrete & Continuous Dynamical Systems - B, 2001, 1 (2) : 183-191. doi: 10.3934/dcdsb.2001.1.183

2018 Impact Factor: 1.313