BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Benoit Perthame

Benoit Perthame appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.APCell Behaviormath.PRnlin.PSPopulations and Evolutionmath.NAmath.OCNeurons and CognitionNumerical AnalysisTissues and Organs

Catalog footprint

What is connected

27works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2026arXiv

Multiscale analysis of a kinetic equation for mechanotaxis

We present a new kinetic equation for cell migration driven by mechanical interactions with the substrate, an effect not previously captured in kinetic models, and essential for explaining observed collective behaviors such as those in bacterial colonies. The model introduces an acceleration term that accounts for the dynamics of motile cells undergoing mechanotaxis, where extracellular signals modulate the forces arising from cell-substrate interactions. From this formulation, we derive a family of macroscopic limit equations and analyze their principal properties. In particular, we examine linear stability and pattern formation ability through theoretical analysis, supported by numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Confined run-and-tumble model with boundary aggregation: long time behavior and convergence to the confined Fokker-Planck model

The motile micro-organisms such as E. coli, sperm, or some seaweed are usually modelled by self-propelled particles that move with the run-and-tumble process. Individual-based stochastic models are usually employed to model the aggregation phenomenon at the boundary, which is an active research field that has attracted a lot of biologists and biophysicists. Self-propelled particles at the microscale have complex behaviors, while characteristics at the population level are more important for practical applications but rely on individual behaviors. Kinetic PDE models that describe the time evolution of the probability density distribution of the motile micro-organisms are widely used. However, how to impose the appropriate boundary conditions that take into account the boundary aggregation phenomena is rarely studied. In this paper, we propose the boundary conditions for a 2D confined run-and-tumble model (CRTM) for self-propelled particle populations moving between two parallel plates with a run-and-tumble process. The proposed model satisfies the relative entropy inequality and thus long-time convergence. We establish the relation between CRTM and the confined Fokker-Planck model (CFPM) studied in [22]. We prove theoretically that when the tumble is highly forward peaked and frequent enough, CRTM converges asymptotically to the CFPM. A numerical comparison of the CRTM with aggregation and CFPM is given. The time evolution of both the deterministic PDE model and individual-based stochastic simulations are displayed, which match each other well.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Structured Model Conserving Biomass for the Size-spectrum Evolution in Aquatic Ecosystems

Mathematical modelling of the evolution of the size-spectrum dynamics in aquatic ecosystems was discovered to be a powerful tool to have a deeper insight into impacts of human- and environmental driven changes on the marine ecosystem. In this article we propose to investigate such dynamics by formulating and investigating a suitable model. The underlying process for these dynamics is given by predation events, causing both growth and death of individuals, while keeping the total biomass within the ecosystem constant. The main governing equation investigated is deterministic and non-local of quadratic type, coming from binary interactions. Predation is assumed to strongly depend on the ratio between a predator and its prey, which is distributed around a preferred feeding preference value. Existence of solutions is shown in dependence of the choice of the feeding preference function as well as the choice of the search exponent, a constant influencing the average volume in water an individual has to search until it finds prey. The equation admits a trivial steady state representing a died out ecosystem, as well as - depending on the parameterregime - steady states with gaps in the size spectrum, giving evidence to the well known cascade effect. The question of stability of these equilibria is considered, showing convergence to the trivial steady state in a certain range of parameters. These analytical observations are underlined by numerical simulations, with additionally exhibiting convergence to the non-trivial equilibrium for specific ranges of parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A Hamilton-Jacobi Approach to Evolution of Dispersal

The evolution of dispersal is a classical question in evolutionary biology, and it has been studied in a wide range of mathematical models. A selection-mutation model, in which the population is structured by space and a phenotypic trait, with the trait acting directly on the dispersal (diffusion) rate, was formulated by Perthame and Souganidis [Math. Model. Nat. Phenom. 11 (2016), 154-166] to study the evolution of random dispersal towards the evolutionarily stable strategy. For the rare mutation limit, it was shown that the equilibrium population concentrates on a single trait associated to the smallest dispersal rate. In this paper, we consider the corresponding evolution equation and characterize the asymptotic behaviors of the time-dependent solutions in the rare mutation limit, under mild convexity assumptions on the underlying Hamiltonian function.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Rare mutations limit of a steady state dispersion trait model

The evolution of dispersal is a classical question in evolutionary ecology, which has been widely studied with several mathematical models. The main question is to define the fittest dispersal rate for a population in a bounded domain, and, more recently, for traveling waves in the full space. In the present study, we reformulate the problem in the context of adaptive evolution. We consider a population structured by space and a genetic trait acting directly on the dispersal (diffusion) rate under the effect of rare mutations on the genetic trait. We show that, as in simpler models, in the limit of vanishing mutations, the population concentrates on a single trait associated to the lowest dispersal rate. We also explain how to compute the evolution speed towards this evolutionary stable distribution. The mathematical interest stems from the asymptotic analysis which requires a completely different treatment of the different variables. For the space variable, the ellipticity leads to the use the maximum principle and Sobolev-type regularity results. For the trait variable, the concentration to a Dirac mass requires a different treatment. This is based on the WKB method and viscosity solutions leading to an effective Hamiltonian (effective fitness of the population) and a constrained Hamilton-Jacobi equation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Free boundary problems for Tumor Growth: a Viscosity solutions approach

The mathematical modeling of tumor growth leads to singular stiff pressure law limits for porous medium equations with a source term. Such asymptotic problems give rise to free boundaries, which, in the absence of active motion, are generalized Hele-Shaw flows. In this note we use viscosity solutions methods to study limits for porous medium-type equations with active motion. We prove the uniform convergence of the density under fairly general assumptions on the initial data, thus improving existing results. We also obtain some additional information/regularity about the propagating interfaces, which, in view of the discontinuities, can nucleate and, thus, change topological type. The main tool is the construction of local, smooth, radial solutions which serve as barriers for the existence and uniqueness results as well as to quantify the speed of propagation of the free boundary propagation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Notes on a PDE System for Biological Network Formation

We present new analytical and numerical results for the elliptic-parabolic system of partial differential equations proposed by Hu and Cai, which models the formation of biological transport networks. The model describes the pressure field using a Darcy's type equation and the dynamics of the conductance network under pressure force effects. Randomness in the material structure is represented by a linear diffusion term and conductance relaxation by an algebraic decay term. The analytical part extends the results of Haskovec, Markowich and Perthame regarding the existence of weak and mild solutions to the whole range of meaningful relaxation exponents. Moreover, we prove finite time extinction or break-down of solutions in the spatially onedimensional setting for certain ranges of the relaxation exponent. We also construct stationary solutions for the case of vanishing diffusion and critical value of the relaxation exponent, using a variational formulation and a penalty method. The analytical part is complemented by extensive numerical simulations. We propose a discretization based on mixed finite elements and study the qualitative properties of network structures for various parameters values. Furthermore, we indicate numerically that some analytical results proved for the spatially one-dimensional setting are likely to be valid also in several space dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A simple derivation of BV bounds for inhomogeneous relaxation systems

We consider relaxation systems of transport equations with heterogeneous source terms and with boundary conditions, which limits are scalar conservation laws. Classical bounds fail in this context and in particular BV estimates. They are the most standard and simplest way to prove compactness and convergence. We provide a novel and simple method to obtain partial BV regularity and strong compactness in this framework. The standard notion of entropy is not convenient either and we also indicate another, but closely related, notion. We give two examples motivated by renal flows which consist of 2 by 2 and 3 by 3 relaxation systems with 2-velocities but the method is more general.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Asymptotic analysis of a selection model with space

Selection of a phenotypical trait can be described in mathematical terms by 'stage structured' equations which are usually written under the form of integral equations so as to express competition for resource between individuals whatever is their trait. The solutions exhibit a concentration effect (selection of the fittest); when a small parameter is introduced they converge to a Dirac mass. An additional space variable can be considered in order to take into account local environmental conditions. Here we assume this environment is a single nutrient which diffuses in the domain. In this framework, we prove that the solution converges to a Dirac mass in the physiological trait which depends on time and on the location in space with Lipschitz continuity. The main difficulties come from the lack of compactness in time and trait variables. Strong convergence can be recovered from uniqueness in the limiting constrained equation after Hopf-Cole change of unknown. Our analysis is motivated by a model of tumor growth introduced by Lorz et al. (preprint) in order to explain emergence of resistance to therapy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Competition and boundary formation in heterogeneous media: Application to neuronal differentiation

We analyze an inhomogeneous system of coupled reaction-diffusion equations representing the dynamics of gene expression during differentiation of nerve cells. The outcome of this developmental phase is the formation of distinct functional areas separated by sharp and smooth boundaries. It proceeds through the competition between the expression of two genes whose expression is driven by monotonic gradients of chemicals, and the products of gene expression undergo local diffusion and drive gene expression in neighboring cells. The problem therefore falls in a more general setting of species in competition within a non-homogeneous medium. We show that in the limit of arbitrarily small diffusion, there exists a unique monotonic stationary solution, which splits the neural tissue into two winner-take-all parts at a precise boundary point: on both sides of the boundary, different neuronal types are present. In order to further characterize the location of this boundary, we use a blow-up of the system and define a traveling wave problem parametrized by the position within the monotonic gradient: the precise boundary location is given by the unique point in space at which the speed of the wave vanishes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Data Assimilation for hyperbolic conservation laws. A Luenberger observer approach based on a kinetic description

Developing robust data assimilation methods for hyperbolic conservation laws is a challenging subject. Those PDEs indeed show no dissipation effects and the input of additional information in the model equations may introduce errors that propagate and create shocks. We propose a new approach based on the kinetic description of the conservation law. A kinetic equation is a first order partial differential equation in which the advection velocity is a free variable. In certain cases, it is possible to prove that the nonlinear conservation law is equivalent to a linear kinetic equation. Hence, data assimilation is carried out at the kinetic level, using a Luenberger observer also known as the nudging strategy in data assimilation. Assimilation then resumes to the handling of a BGK type equation. The advantage of this framework is that we deal with a single "linear" equation instead of a nonlinear system and it is easy to recover the macroscopic variables. The study is divided into several steps and essentially based on functional analysis techniques. First we prove the convergence of the model towards the data in case of complete observations in space and time. Second, we analyze the case of partial and noisy observations. To conclude, we validate our method with numerical results on Burgers equation and emphasize the advantages of this method with the more complex Saint-Venant system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Incompressible limit of mechanical model of tumor growth with viscosity

Various models of tumor growth are available in the litterature. A first class describes the evolution of the cell number density when considered as a continuous visco-elastic material with growth. A second class, describes the tumor as a set and rules for the free boundary are given related to the classical Hele-Shaw model of fluid dynamics. Following the lines of previous papers where the material is described by a purely elastic material, or when active cell motion is included, we make the link between the two levels of description considering the 'stiff pressure law' limit. Even though viscosity is a regularizing effect, new mathematical difficulties arise in the visco-elastic case because estimates on the pressure field are weaker and do not imply immediately compactness. For instance, traveling wave solutions and numerical simulations show that the pressure may be discontinous in space which is not the case for the elastic case.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Long-term analysis of phenotypically structured models

Phenotypically structured equations arise in population biology to describe the interaction of species with their environment that brings the nutrients. This interaction usually leads to selection of the fittest individuals. Models used in this area are highly nonlinear, and the question of long term behaviour is usually not solved. However, there is a particular class of models for which convergence to an Evolutionary Stable Distribution is proved, namely when the quasi-static assumption is made. This means that the environment, and thus the nutrient supply, reacts immediately to the population dynamics. One possible proof is based on a Total Variation bound for the appropriate quantity. We extend this proof to several cases where the nutrient is regenerated with delay. A simple example is the chemostat with a rendering factor, then our result does not use any smallness assumption. For a more general setting, we can treat the case with a fast reaction of nutrient supply to the population dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mathematical Analysis of a System for Biological Network Formation

Motivated by recent physics papers describing rules for natural network formation, we study an elliptic-parabolic system of partial differential equations proposed by Hu and Cai. The model describes the pressure field thanks to Darcy's type equation and the dynamics of the conductance network under pressure force effects with a diffusion rate $D$ representing randomness in the material structure. We prove the existence of global weak solutions and of local mild solutions and study their long term behaviour. It turns out that, by energy dissipation, steady states play a central role to understand the pattern capacity of the system. We show that for a large diffusion coefficient $D$, the zero steady state is stable. Patterns occur for small values of $D$ because the zero steady state is Turing unstable in this range; for $D=0$ we can exhibit a large class of dynamically stable (in the linearized sense) steady states.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scalar conservation laws with rough (stochastic) fluxes; the spatially dependent case

We continue the development of the theory of pathwise stochastic entropy solutions for scalar conservation laws in $\R^N$ with quasilinear multiplicative ''rough path'' dependence by considering inhomogeneous fluxes and a single rough path like, for example, a Brownian motion. Following our previous note where we considered spatially independent fluxes, we introduce the notion of pathwise stochastic entropy solutions and prove that it is well posed, that is we establish existence, uniqueness and continuous dependence in the form of a (pathwise) $L^1$-contraction. Our approach is motivated by the theory of stochastic viscosity solutions, which was introduced and developed by two of the authors, to study fully nonlinear first- and second-order stochastic pde with multiplicative noise. This theory relies on special test functions constructed by inverting locally the flow of the stochastic characteristics. For conservation laws this is best implemented at the level of the kinetic formulation which we follow here.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Transversal instability for the thermodiffusive reaction-diffusion system

The propagation of unstable interfaces is at the origin of remarkable patterns that are observed in various areas of science as chemical reactions, phase transitions, growth of bacterial colonies. Since a scalar equation generates usually stable waves, the simplest mathematical description relies on two by two reaction-diffusion systems. Our interest is the extension of the Fisher/KPP equation to a two species reaction which represents reactant concentration and temperature when used for flame propagation, bacterial population and nutrient concentration when used in biology. We study in which circumstances instabilities can occur and in particular the effect of dimension. It is observed numerically that spherical waves can be unstable depending on the coefficients. A simpler mathematical framework is to study transversal instability, that means a one dimensional wave propagating in two space dimensions. Then, explicit analytical formulas give explicitely the range of paramaters for instability.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On a voltage-conductance kinetic system for integrate and fire neural networks

The voltage-conductance kinetic equation for integrate and fire neurons has been used in neurosciences since a decade and describes the probability density of neurons in a network. It is used when slow conductance receptors are activated and noticeable applications to the visual cortex have been worked-out. In the simplest case, the derivation also uses the assumption of fully excitatory and moderately all-to-all coupled networks; this is the situation we consider here. We study properties of solutions of the kinetic equation for steady states and time evolution and we prove several global a priori bounds both on the probability density and the firing rate of the network. The main difficulties are related to the degeneracy of the diffusion resulting from noise and to the quadratic aspect of the nonlinearity. This result constitutes a paradox; the solutions of the kinetic model, of partially hyperbolic nature, are globally bounded but it has been proved that the fully parabolic integrate and fire equation (some kind of diffusion limit of the former) blows-up in finite time.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Scalar conservation laws with rough (stochastic) fluxes

We develop a pathwise theory for scalar conservation laws with quasilinear multiplicative rough path dependence, a special case being stochastic conservation laws with quasilinear stochastic dependence. We introduce the notion of pathwise stochastic entropy solutions, which is closed with the local uniform limits of paths, and prove that it is well posed, i.e., we establish existence, uniqueness and continuous dependence, in the form of pathwise $L^1$-contraction, as well as some explicit estimates. Our approach is motivated by the theory of stochastic viscosity solutions, which was introduced and developed by two of the authors, to study fully nonlinear first- and second-order stochastic pde with multiplicative noise. This theory relies on special test functions constructed by inverting locally the flow of the stochastic characteristics. For conservation laws this is best implemented at the level of the kinetic formulation which we follow here.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Time fluctuations in a population model of adaptive dynamics

We study the dynamics of phenotypically structured populations in environments with fluctuations. In particular, using novel arguments from the theories of Hamilton-Jacobi equations with constraints and homogenization, we obtain results about the evolution of populations in environments with time oscillations, the development of concentrations in the form of Dirac masses, the location of the dominant traits and their evolution in time. Such questions have already been studied in time homogeneous environments. More precisely we consider the dynamics of a phenotypically structured population in a changing environment under mutations and competition for a single resource. The mathematical model is a non-local parabolic equation with a periodic in time reaction term. We study the asymptotic behavior of the solutions in the limit of small diffusion and fast reaction. Under concavity assumptions on the reaction term, we prove that the solution converges to a Dirac mass whose evolution in time is driven by a Hamilton-Jacobi equation with constraint and an effective growth/death rate which is derived as a homogenization limit. We also prove that, after long-time, the population concentrates on a trait where the maximum of an effective growth rate is attained. Finally we provide an example showing that the time oscillations may lead to a strict increase of the asymptotic population size.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Invasion fronts with variable motility: phenotype selection, spatial sorting and wave acceleration

Invasion fronts in ecology are well studied but very few mathematical results concern the case with variable motility (possibly due to mutations). Based on an apparently simple reaction-diffusion equation, we explain the observed phenomena of front acceleration (when the motility is unbounded) as well as other quantitative results, such as the selection of the most motile individuals (when the motility is bounded). The key argument for the construction and analysis of traveling fronts is the derivation of the dispersion relation linking the speed of the wave and the spatial decay. When the motility is unbounded we show that the position of the front scales as $t^{3/2}$. When the mutation rate is low we show that the canonical equation for the dynamics of the fittest trait should be stated as a PDE in our context. It turns out to be a type of Burgers equation with source term.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Populational adaptive evolution, chemotherapeutic resistance and multiple anti-cancer therapies

Resistance to chemotherapies, particularly to anticancer treatments, is an increasing medical concern. Among the many mechanisms at work in cancers, one of the most important is the selection of tumor cells expressing resistance genes or phenotypes. Motivated by the theory of mutation-selection in adaptive evolution, we propose a model based on a continuous variable that represents the expression level of a resistance gene (or genes, yielding a phenotype) influencing in healthy and tumor cells birth/death rates, effects of chemotherapies (both cytotoxic and cytostatic) and mutations. We extend previous work by demonstrating how qualitatively different actions of chemotherapeutic and cytostatic treatments may induce different levels of resistance. The mathematical interest of our study is in the formalism of constrained Hamilton-Jacobi equations in the framework of viscosity solutions. We derive the long-term temporal dynamics of the fittest traits in the regime of small mutations. In the context of adaptive cancer management, we also analyse whether an optimal drug level is better than the maximal tolerated dose.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stochastic averaging lemmas for kinetic equations

We develop a class of averaging lemmas for stochastic kinetic equations. The velocity is multiplied by a white noise which produces a remarkable change in time scale. Compared to the deterministic case and as far as we work in $L^2$, the nature of regularity on averages is not changed in this stochastic kinetic equation and stays in the range of fractional Sobolev spaces at the price of an additional expectation. However all the exponents are changed; either time decay rates are slower (when the right hand side belongs to $L^2$), or regularity is better when the right hand side contains derivatives. These changes originate from a different space/time scaling in the deterministic and stochastic cases. Our motivation comes from scalar conservation laws with stochastic fluxes where the structure under consideration arises naturally through the kinetic formulation of scalar conservation laws.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Singular Hamilton-Jacobi equation for the tail problem

In this paper we study the long time-long range behavior of reaction diffusion equations with negative square root -type reaction terms. In particular we investigate the exponential behavior of the solutions after a standard hyperbolic scaling. This leads to a Hamilton-Jacobi variational inequality with an obstacle that depends on the solution itself and defines the open set where the limiting solution does not vanish. Counter-examples show a nontrivial lack of uniqueness for the variational inequality depending on the conditions imposed on the boundary of this open set. Both Dirichlet and state constraints boundary conditions play a role. When the competition term does not change sign, we can identify the limit, while, in general, we find lower and upper bounds for the limit. Although models of this type are rather old and extinction phenomena are as important as blow-up, our motivation comes from the so-called "tail problem" in population biology. One way to avoid meaningless exponential tails, is to impose extra-mortality below a given survival threshold. Our study shows that the precise form of this extra-mortality term is asymptotically irrelevant and that, in the survival zone, the population profile is impacted by the survival threshold (except in the very particular case when the competition term is non-positive).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Structured Population Model of Cell Differentiation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Can a traveling wave connect two unstable states? The case of the nonlocal Fisher equation

This note investigates the properties of the traveling waves solutions of the nonlocal Fisher equation. The existence of such solutions has been proved recently in \cite{BNPR} but their asymptotic behavior was still unclear. We use here a new numerical approximation of these traveling waves which shows that some traveling waves connect the two homogeneous steady states $0$ and $1$, which is a striking fact since $0$ is dynamically unstable and $1$ is unstable in the sense of Turing.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modelling the spatial organization of cell proliferation in the developing central nervous system

How far is neuroepithelial cell proliferation in the developing central nervous system a deterministic process? Or, to put it in a more precise way, how accurately can it be described by a deterministic mathematical model? To provide tracks to answer this question, a deterministic system of transport and diffusion partial differential equations, both physiologically and spatially structured, is introduced as a model to describe the spatially organized process of cell proliferation during the development of the central nervous system. As an initial step towards dealing with the three-dimensional case, a unidimensional version of the model is presented. Numerical analysis and numerical tests are performed. In this work we also achieve a first experimental validation of the proposed model, by using cell proliferation data recorded from histological sections obtained during the development of the optic tectum in the chick embryo.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Traveling plateaus for a hyperbolic Keller-Segel system with attraction and repulsion: existence and branching instabilities

How can repulsive and attractive forces, acting on a conservative system, create stable traveling patterns or branching instabilities? We have proposed to study this question in the framework of the hyperbolic Keller-Segel system with logistic sensitivity. This is a model system motivated by experiments on cell communities auto-organization, a field which is also called socio-biology. We continue earlier modeling work, where we have shown numerically that branching patterns arise for this system and we have analyzed this instability by formal asymptotics for small diffusivity of the chemo-repellent. Here we are interested in the more general situation, where the diffusivities of both the chemo-attractant and the chemo-repellent are positive. To do so, we develop an appropriate functional analysis framework. We apply our method to two cases. Firstly we analyze steady states. Secondly we analyze traveling waves when neglecting the degradation coefficient of the chemo-repellent; the unique wave speed appears through a singularity cancelation which is the main theoretical difficulty. This shows that in different situations the cell density takes the shape of a plateau. The existence of steady states and traveling plateaus are a symptom of how rich the system is and why branching instabilities can occur. Numerical tests show that large plateaus may split into smaller ones, which remain stable.

0 citations0 reviewsNo code linkedOpen access