Eigenelements of a General Aggregation-Fragmentation Model

We consider a linear integro-differential equation which arises to describe both aggregation-fragmentation processes and cell division. We prove the existence of a solution (\lb,\U,\phi) to the related eigenproblem. Such eigenelements are useful to study the long time asymptotic behaviour of solutions as well as the steady states when the equation is coupled with an ODE. Our study concerns a non-constant transport term that can vanish at $x=0,$ since it seems to be relevant to describe some biological processes like proteins aggregation. Non lower-bounded transport terms bring difficulties to find $a\ priori$ estimates. All the work of this paper is to solve this problem using weighted-norms

Biomixing by chemotaxis and efficiency of biological reactions: the critical reaction case

Many phenomena in biology involve both reactions and chemotaxis. These processes can clearly influence each other, and chemotaxis can play an important role in sustaining and speeding up the reaction. In continuation of our earlier work, we consider a model with a single density function involving diffusion, advection, chemotaxis, and absorbing reaction. The model is motivated, in particular, by the studies of coral broadcast spawning, where experimental observations of the efficiency of fertilization rates significantly exceed the data obtained from numerical models that do not take chemotaxis (attraction of sperm gametes by a chemical secreted by egg gametes) into account. We consider the case of the weakly coupled quadratic reaction term, which is the most natural from the biological point of view and was left open. The result is that similarly to higher power coupling, the chemotaxis plays a crucial role in ensuring efficiency of reaction. However, mathematically, the picture is quite different in the quadratic reaction case and is more subtle. The reaction is now complete even in the absence of chemotaxis, but the timescales are very different. Without chemotaxis, the reaction is very slow, especially for the weak reaction coupling coefficient. With chemotaxis, the timescale and efficiency of reaction are independent of the coupling parameter.Comment: 10 pages. arXiv admin note: text overlap with arXiv:1101.244

On the Calibration of a Size-Structured Population Model from Experimental Data

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

The flashing ratchet and unidirectional transport of matter

We study the flashing ratchet model of a Brownian motor, which consists in cyclical switching between the Fokker-Planck equation with an asymmetric ratchet-like potential and the pure diffusion equation. We show that the motor really performs unidirectional transport of mass, for proper parameters of the model, by analyzing the attractor of the problem and the stationary vector of a related Markov chain.Comment: 11 page

The one-dimensional Keller-Segel model with fractional diffusion of cells

We investigate the one-dimensional Keller-Segel model where the diffusion is replaced by a non-local operator, namely the fractional diffusion with exponent $0<\alpha\leq 2$. We prove some features related to the classical two-dimensional Keller-Segel system: blow-up may or may not occur depending on the initial data. More precisely a singularity appears in finite time when $\alpha<1$ and the initial configuration of cells is sufficiently concentrated. On the opposite, global existence holds true for $\alpha\leq1$ if the initial density is small enough in the sense of the $L^{1/\alpha}$ norm.Comment: 12 page

Three-points interfacial quadrature for geometrical source terms on nonuniform grids

International audienceThis paper deals with numerical (finite volume) approximations, on nonuniform meshes, for ordinary differential equations with parameter-dependent fields. Appropriate discretizations are constructed over the space of parameters, in order to guarantee the consistency in presence of variable cells' size, for which $L^p$-error estimates, $1\le p < +\infty$, are proven. Besides, a suitable notion of (weak) regularity for nonuniform meshes is introduced in the most general case, to compensate possibly reduced consistency conditions, and the optimality of the convergence rates with respect to the regularity assumptions on the problem's data is precisely discussed. This analysis attempts to provide a basic theoretical framework for the numerical simulation on unstructured grids (also generated by adaptive algorithms) of a wide class of mathematical models for real systems (geophysical flows, biological and chemical processes, population dynamics)

Phase appearance or disappearance in two-phase flows

This paper is devoted to the treatment of specific numerical problems which appear when phase appearance or disappearance occurs in models of two-phase flows. Such models have crucial importance in many industrial areas such as nuclear power plant safety studies. In this paper, two outstanding problems are identified: first, the loss of hyperbolicity of the system when a phase appears or disappears and second, the lack of positivity of standard shock capturing schemes such as the Roe scheme. After an asymptotic study of the model, this paper proposes accurate and robust numerical methods adapted to the simulation of phase appearance or disappearance. Polynomial solvers are developed to avoid the use of eigenvectors which are needed in usual shock capturing schemes, and a method based on an adaptive numerical diffusion is designed to treat the positivity problems. An alternate method, based on the use of the hyperbolic tangent function instead of a polynomial, is also considered. Numerical results are presented which demonstrate the efficiency of the proposed solutions

On interfaces between cell populations with different mobilities

Partial differential equations describing the dynamics of cell population densities from a fluid mechanical perspective can model the growth of avascular tumours. In this framework, we consider a system of equations that describes the interaction between a population of dividing cells and a population of non-dividing cells. The two cell populations are characterised by different mobilities. We present the results of numerical simulations displaying two-dimensional spherical waves with sharp interfaces between dividing and non-dividing cells. Furthermore, we numerically observe how different ratios between the mobilities change the morphology of the interfaces, and lead to the emergence of finger-like patterns of invasion above a threshold. Motivated by these simulations, we study the existence of one-dimensional travelling wave solutions

Two classes of nonlocal Evolution Equations related by a shared Traveling Wave Problem

We consider reaction-diffusion equations and Korteweg-de Vries-Burgers (KdVB) equations, i.e. scalar conservation laws with diffusive-dispersive regularization. We review the existence of traveling wave solutions for these two classes of evolution equations. For classical equations the traveling wave problem (TWP) for a local KdVB equation can be identified with the TWP for a reaction-diffusion equation. In this article we study this relationship for these two classes of evolution equations with nonlocal diffusion/dispersion. This connection is especially useful, if the TW equation is not studied directly, but the existence of a TWS is proven using one of the evolution equations instead. Finally, we present three models from fluid dynamics and discuss the TWP via its link to associated reaction-diffusion equations

A General Inverse Problem for the Growth-Fragmentation Equation

The growth-fragmentation equation arises in many different contexts, ranging from cell division, protein polymerization, biopolymers, neurosciences etc. Direct observation of temporal dynamics being often difficult, it is of main interest to develop theoretical and numerical methods to recover reaction rates and parameters of the equation from indirect observation of the solution. Following the work done in (Perthame, Zubelli, 2006) and (Doumic, Perthame, Zubelli, 2009) for the specific case of the cell division equation, we address here the general question of recovering the fragmentation rate of the equation from the observation of the time-asymptotic solution, when the fragmentation kernel and the growth rates are fully general. We give both theoretical results and numerical methods, and discuss the remaining issues