Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userOn the Calibration of a Size-Structured Population Model from Experimental Data
117
0
0.0
(
0
)
Authors
Marie Doumic Jauffret
Ask ChatGPT about the research
No Arabic abstract
The aim of this work is twofold. First, we survey the techniques developed in (Perthame, Zubelli, 2007) and (Doumic, Perthame, Zubelli, 2008) to reconstruct the division (birth) rate from the cell volume distribution data in certain structured population models. Secondly, we implement such techniques on experimental cell volume distributions available in the literature so as to validate the theoretical and numerical results. As a proof of concept, we use the data reported in the classical work of Kubitschek [3] concerning Escherichia coli in vitro experiments measured by means of a Coulter transducer-multichannel analyzer system (Coulter Electronics, Inc., Hialeah, Fla, USA.) Despite the rather old measurement technology, the reconstructed division rates still display potentially useful biological features.
rate research
Read More
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.
We consider a size-structured model for cell division and address the question of determining the division (birth) rate from the measured stable size distribution of the population. We propose a new regularization technique based on a filtering approach. We prove convergence of the algorithm and validate the theoretical results by implementing numerical simulations, based on classical techniques. We compare the results for direct and inverse problems, for the filtering method and for the quasi-reversibility method proposed in [Perthame-Zubelli].
We consider a size-structured aggregation and growth model of phytoplankton community proposed by Ackleh and Fitzpatrick [2]. The model accounts for basic biological phenomena in phytoplankton community such as growth, gravitational sedimentation, predation by zooplankton, fecundity, and aggregation. Our primary goal in this paper is to investigate the long-term behavior of the proposed aggregation and growth model. Particularly, using the well-known principle of linearized stability and semigroup compactness arguments, we provide sufficient conditions for local exponential asymptotic stability of zero solution as well as sufficient conditions for instability. We express these conditions in the form of an easy to compute characteristic function, which depends on the functional relationship between growth, sedimentation and fecundity. Our results can be used to predict long-term phytoplankton dynamic
We introduce and analyze several aspects of a new model for cell differentiation. It assumes that differentiation of progenitor cells is a continuous process. From the mathematical point of view, it is based on partial differential equations of transport type. Specifically, it consists of a structured population equation with a nonlinear feedback loop. This models the signaling process due to cytokines, which regulate the differentiation and proliferation process. We compare the continuous model to its discrete counterpart, a multi-compartmental model of a discrete collection of cell subpopulations recently proposed by Marciniak-Czochra et al. in 2009 to investigate the dynamics of the hematopoietic system. We obtain uniform bounds for the solutions, characterize steady state solutions, and analyze their linearized stability. We show how persistence or extinction might occur according to values of parameters that characterize the stem cells self-renewal. We also perform numerical simulations and discuss the qualitative behavior of the continuous model vis a vis the discrete one.
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 also be used to produce 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 Smoluchowksi 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).
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
