Search

In this work, we introduce a new semi-Lagrangian numerical method proposed to solve a cell population balance model which describes cell dwarfism, by allowing cell division at any size. We analyse its convergence and derive an optimal rate. Numerical experiments are reported to demonstrate the predicted accuracy of the scheme. Finally, the good behaviour of the numerical method is exhibited in stressed conditions that can provoke lack of smoothness and simulations are included to show how the increase in the division rate of small size-cells promotes dynamics in which nonfunctional dwarf cells are saturating the total cell population. [ABSTRACT FROM AUTHOR]

We deeply researched into the asymptotic behaviour of a numerical method adapted for the solution of mathematical model of hematopoiesis which describes the dynamics of a stem cell population. We investigated the stationary solutions of the original model by their numerical approximation: we proved the existence of a numerical stationary solution that provides a good approximation to the nontrivial equilibrium solution of the problem. Also, we presented a numerical simulation which confirms this behaviour. [ABSTRACT FROM PUBLISHER]

In this paper, we analyze the convergence of a second-order numerical method for the approximation of a size-structured population model whose dependency on the environment is managed by the evolution of a vital resource. Optimal convergence rate is derived. Numerical experiments are also reported to demonstrate the predicted accuracy of the scheme. Also, it is applied to solve a problem that describes the dynamics of a Daphnia magna population, paying attention to the unstable case. [ABSTRACT FROM AUTHOR]

We present an efficiency study for autonomous and non-autonomous non-linear size-structured population models. The study considers three numerical methods: a characteristics scheme, the Lax-Wendroff method and the box method, which are completely described in the paper. Five test problems are considered with diverse degree of complexity: non-trivial equilibrium, periodic solutions and diverse growth functions. The study of the efficiency takes into account the properties of the numerical schemes (such as the stability) and uses a multiple regression analysis to determine the constants of the leading terms of the corresponding global errors. We show how the tables of errors and cpu-times can be used to explain the meaning of the efficiency results. In addition, we present the convergence analysis of the box method.

We formulate explicit second-order finite difference schemes for the numerical integration of non-linear age-dependent population models. These methods have been designed by means of a representation formula for the theoretical solution of the integro-differential equation joint with open quadrature formulae for the numerical approximation of non-local terms. The schemes are analyzed and some numerical experiments are also reported in order to show numerically their accuracy.

We formulate schemes for the numerical solution of size-dependent population models. Such schemes discretize size by means of a natural grid, which introduces a discrete dynamics. The schemes are analysed and optimal rates of convergence are derived. Some numerical experiments are also reported to demonstrate the predicted accuracy of the schemes.