Analysis of a semi-Lagrangian numerical method in a cell dwarfism model

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.