Finite difference approximations for a size-structured population model with distributed states in the recruitment

In this paper we consider a size-structured population model where individuals may be recruited into the population at different sizes. First and second order finite difference schemes are developed to approximate the solution of the mathematical model. The convergence of the approximations to a unique weak solution with bounded total variation is proved. We then show that as the distribution of the new recruits become concentrated at the smallest size, the weak solution of the distributed states-at-birth model converges to the weak solution of the classical Gurtin-McCamy-type size-structured model in the weak$^*$ topology. Numerical simulations are provided to demonstrate the achievement of the desired accuracy of the two methods for smooth solutions as well as the superior performance of the second-order method in resolving solution-discontinuities. Finally we provide an example where supercritical Hopf-bifurcation occurs in the limiting single state-at-birth model and we apply the second-order numerical scheme to show that such bifurcation occurs in the distributed model as well.

Steady states in hierarchical structured populations with distributed states at birth

We investigate steady states of a quasilinear first order hyperbolic partial integro-differential equation. The model describes the evolution of a hierarchical structured population with distributed states at birth. Hierarchical size-structured models describe the dynamics of populations when individuals experience size-specific environment. This is the case for example in a population where individuals exhibit cannibalistic behavior and the chance to become prey (or to attack) depends on the individuals size. The other distinctive feature of the model is that individuals are recruited into the population at arbitrary size. This amounts to an infinite rank integral operator describing the recruitment process. First we establish conditions for the existence of a positive steady state of the model. Our method uses a fixed point result of nonlinear maps in conical shells of Banach spaces. Then we study stability properties of steady states for the special case of a separable growth rate using results from the theory of positive operators on Banach lattices.