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

Nicolas Vauchelet

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

ResearcherUnclaimed source record
math.APmath.NAmath.OCPopulations and EvolutionCell Behaviormath.CAnlin.PSphysics.soc-ph

Catalog footprint

What is connected

26works
8topics
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

26 published item(s)

preprint2026arXiv

Optimal Dirac controls for time-periodic bistable ODEs, application to population replacement

This work addresses an optimal control problem on a dynamics governed by a nonlinear differential equation with a bistable time-periodic nonlinearity. This problem, relevant in population dynamics, models the strategy of replacing a population of A-type individuals by a population of B-type individuals in a time-varying environment, focusing on the evolution of the proportion of B-type individuals among the whole population. The control term accounts for the instant release of B-type individuals. Our main goal, after noting some interesting properties on the differential equation, is to determine the optimal time at which this release should be operated to ensure population replacement while minimizing the release effort. The results establish that the optimal release time appears to be the minimizer of a function involving the carrying capacity of the environment and the threshold periodic solution of the dynamics; they also describe the convergence of the whole optimal release strategy. An application to the biocontrol of mosquito populations using Wolbachia-infected individuals illustrates the relevance of the theoretical results. Wolbachia is a bacterium that helps preventing the transmission of some viruses from mosquitoes to humans, making the optimization of Wolbachia propagation in a mosquito population a crucial issue.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical scheme for kinetic transport equation with internal state *

We investigate the numerical discretization of a two-stream kinetic system with an internal state, such system has been introduced to model the motion of cells by chemotaxis. This internal state models the intracellular methylation level. It adds a variable in the mathematical model, which makes it more challenging to simulate numerically. Moreover, it has been shown that the macroscopic or mesoscopic quantities computed from this system converge to the Keller-Segel system at diffusive scaling or to the velocity-jump kinetic system for chemotaxis at hyperbolic scaling. Then we pay attention to propose numerical schemes uniformly accurate with respect to the scaling parameter. We show that these schemes converge to some limiting schemes which are consistent with the limiting macroscopic or kinetic system. This study is illustrated with some numerical simulations and comparisons with Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the role of the epithelium in a model of sodium exchange in renal tubules

In this study we present a mathematical model describing the transport of sodium in a fluid circulating in a counter-current tubular architecture, which constitutes a simplified model of Henle's loop in a kidney nephron. The model explicitly takes into account the epithelial layer at the interface between the tubular lumen and the surrounding interstitium. In a specific range of parameters, we show that explicitly accounting for transport across the apical and basolateral membranes of epithelial cells, instead of assuming a single barrier, affects the axial concentration gradient , an essential determinant of the urinary concentrating capacity. We present the solution related to the stationary system, and we perform numerical simulations to understand the physiological behaviour of the system. We prove that when time grows large, our dynamic model converges towards the stationary system at an exponential rate. In order to prove rigorously this global asymptotic stability result, we study eigen-problems of an auxiliary linear operator and its dual.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optimal immunity control by social distancing for the SIR epidemic model

Until a vaccine or therapy is found against the SARS-CoV-2 coronavirus, reaching herd immunity appears to be the only mid-term option. However, if the number of infected individuals decreases and eventually fades only beyond this threshold, a significant proportion of susceptible may still be infected until the epidemic is over. A containment strategy is likely the best policy in the worst case where no vaccine or therapy is found. In order to keep the number of newly infected persons to a minimum, a possible strategy is to apply strict containment measures, so that the number of susceptible individuals remains close to herd immunity. Such an action is unrealistic since containment can only last for a finite amount of time and is never total. In this article, using a classical SIR model, we determine the (partial or total) containment strategy on a given finite time interval that maximizes the number of susceptible individuals over an infinite horizon, or equivalently that minimizes the total infection burden during the curse of the epidemic. The existence and uniqueness of the optimal strategy is proved and the latter is fully characterized. If applicable in practice, such a strategy would lead theoretically to an increase by 30% of the proportion of susceptible on an infinite horizon, for a containment level corresponding to the sanitary measures put in place in France from March to May 2020. We also analyze the minimum intervention time to reach a fixed distance from herd immunity, and show the relationship with the previous problem. Simulations are provided that illustrate and validate the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reduction of a model for sodium exchanges in kidney nephron

This work deals with a mathematical analysis of sodium's transport in a tubular architecture of a kidney nephron. The nephron is modelled by two counter-current tubules. Ionic exchange occurs at the interface between the tubules and the epithelium and between the epithelium and the surrounding environment (interstitium). From a mathematical point of view, this model consists of a 5x5 semi-linear hyperbolic system. In the literature similar models neglect the epithelial layers. In this paper, we show rigorously that such models may be obtained by assuming that the permeabilities between lumen and epithelium are large. Indeed we show that when these grow, solutions of the 5x5 system converge in a certain way to solutions of a reduced 3x3 system where no epithelial layer is present. The problem is dened on a bounded spacial domain with initial and boundary data. Establishing BV compactness forces to introduce initial layers and to handle carefully the presence of lateral boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The sterile insect technique used as a barrier control against reinfestation

The sterile insect technique consists in massive release of sterilized males in the aim to reduce the size of mosquitoes population or even eradicate it. In this work, we investigate the feasability of using the sterile insect technique as a barrier against reinvasion. More precisely, we provide some numerical simulations and mathematical results showing that performing the sterile insect technique on a band large enough may stop reinvasion.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convergence order of upwind type schemes for transport equations with discontinuous coefficients

An analysis of the error of the upwind scheme for transport equation with discontinuous coefficients is provided. We consider here a velocity field that is bounded and one-sided Lipschitz continuous. In this framework, solutions are defined in the sense of measures along the lines of Poupaud and Rascle's work. We study the convergence order of the upwind scheme in the Wasserstein distances. More precisely, we prove that in this setting the convergence order is 1/2. We also show the optimality of this result. In the appendix, we show that this result also applies to other "diffusive" "first order" schemes and to a forward semi-Lagrangian scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

From Kinetic Theory of Multicellular Systems to Hyperbolic Tissue Equations: Asymptotic Limits and Computing

This paper deals with the analysis of the asymptotic limit toward the derivation of macroscopic equations for a class of equations modeling complex multicellular systems by methods of the kinetic theory. After having chosen an appropriate scaling of time and space, a Chapman-Enskog expansion is combined with a closed, by minimization, technique to derive hyperbolic models at the macroscopic level. The resulting macroscopic equations show how the macroscopic tissue behavior can be described by hyperbolic systems which seem the most natural in this context. We propose also an asymptotic-preserving well-balanced scheme for the one-dimensional hyperbolic model, in the two dimensional case, we consider a time splitting method between the conservative part and the source term where the conservative equation is approximated by the Lax-Friedrichs scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Reduction to a single closed equation for 2 by 2 reaction-diffusion systems of Lotka-Volterra type

We consider general models of coupled reaction-diffusion systems for interacting variants of the same species. When the total population becomes large with intensive competition, we prove that the frequencies (i.e. proportions) of the variants can be approached by the solution of a simpler reaction-diffusion system, through a singular limit method and a relative compactness argument. As an example of application, we retrieve the classical bistable equation for Wolbachia's spread into an arthropod population from a system modeling interaction between infected and uninfected individuals.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Traveling Pulses for a Two-Species Chemotaxis Model

Mathematical models have been widely used to describe the collective movement of bacteria by chemotaxis. In particular, bacterial concentration waves traveling in a narrow channel have been experimentally observed and can be precisely described thanks to a mathematical model at the macroscopic scale. Such model was derived in [1] using a kinetic model based on an accurate description of the mesoscopic run-and-tumble process. We extend this approach to study the behavior of the interaction between two populations of E. Coli. Separately, each population travels with its own speed in the channel. When put together, a synchronization of the speed of the traveling pulses can be observed. We show that this synchronization depends on the fraction of the fast population. Our approach is based on mathematical analysis of a macroscopic model of partial differential equations. Numerical simulations in comparison with experimental observations show qualitative agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Derivation of the bacterial run-and-tumble kinetic equation from a model with biochemical pathway

Kinetic-transport equations are, by now, standard models to describe the dynamics of populations of bacteria moving by run-and-tumble. Experimental observations show that bacteria increase their run duration when encountering an increasing gradient of chemotactic molecules. This led to a first class of models which heuristically include tumbling frequencies depending on the path-wise gradient of chemotactic signal. More recently, the biochemical pathways regulating the flagellar motors were uncovered. This knowledge gave rise to a second class of kinetic-transport equations, that takes into account an intra-cellular molecular content and which relates the tumbling frequency to this information. It turns out that the tumbling frequency depends on the chemotactic signal, and not on its gradient. For these two classes of models, macroscopic equations of Keller-Segel type, have been derived using diffusion or hyperbolic rescaling. We complete this program by showing how the first class of equations can be derived from the second class with molecular content after appropriate rescaling. The main difficulty is to explain why the path-wise gradient of chemotactic signal can arise in this asymptotic process. Randomness of receptor methylation events can be included, and our approach can be used to compute the tumbling frequency in presence of such a noise.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

One-dimensional aggregation equation after blow up: existence, uniqueness and numerical simulation

The nonlocal nonlinear aggregation equation in one space dimension is investigated. In the so-called attractive case smooth solutions blow up in finite time, so that weak measure solutions are introduced. The velocity involved in the equation becomes discontinuous, and a particular care has to be paid to its definition as well as the formulation of the corresponding flux. When this is done, the notion of duality solutions allows to obtain global in time existence and uniqueness for measure solutions. An upwind finite volume scheme is also analyzed, and the convergence towards the unique solution is proved. Numerical examples show the dynamics of the solutions after the blow up time.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Synchronising and Non-synchronising Dynamics for a Two-species Aggregation Model

This paper deals with analysis and numerical simulations of a one-dimensional two-species hyperbolic aggregation model. This model is formed by a system of transport equations with nonlocal velocities, which describes the aggregate dynamics of a two-species population in interaction appearing for instance in bacterial chemotaxis. Blow-up of classical solutions occurs in finite time. This raises the question to define measure-valued solutions for this system. To this aim, we use the duality method developed for transport equations with discontinuous velocity to prove the existence and uniqueness of measure-valued solutions. The proof relies on a stability result. In addition, this approach allows to study the hyperbolic limit of a kinetic chemotaxis model. Moreover, we propose a finite volume numerical scheme whose convergence towards measure-valued solutions is proved. It allows for numerical simulations capturing the behaviour after blow up. Finally, numerical simulations illustrate the complex dynamics of aggregates until the formation of a single aggregate: after blow-up of classical solutions, aggregates of different species are synchronising or nonsynchronising when collide, that is move together or separately, depending on the parameters of the model and masses of species involved.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Derivation of a Hele-Shaw type system from a cell model with active motion

We formulate a Hele-Shaw type free boundary problem for a tumor growing under the combined effects of pressure forces, cell multiplication and active motion, the latter being the novelty of the present paper. This new ingredient is considered here as a standard diffusion process. The free boundary model is derived from a description at the cell level using the asymptotic of a stiff pressure limit. Compared to the case when active motion is neglected, the pressure satisfies the same complementarity Hele-Shaw type formula. However, the cell density is smoother (Lipschitz continuous), while there is a deep change in the free boundary velocity, which is no longer given by the gradient of the pressure, because some kind of 'mushy region' prepares the tumor invasion.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Equivalence between duality and gradient flow solutions for one-dimensional aggregation equations

Existence and uniqueness of global in time measure solution for a one dimensional nonlinear aggregation equation is considered. Such a system can be written as a conservation law with a velocity field computed through a selfconsistant interaction potential. Blow up of regular solutions is now well established for such system. In Carrillo et al. (Duke Math J (2011)), a theory of existence and uniqueness based on the geometric approach of gradient flows on Wasserstein space has been developped. We propose in this work to establish the link between this approach and duality solutions. This latter concept of solutions allows in particular to define a flow associated to the velocity field. Then an existence and uniqueness theory for duality solutions is developped in the spirit of James and Vauchelet (NoDEA (2013)). However, since duality solutions are only known in one dimension, we restrict our study to the one dimensional case.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Existence and diffusive limit of a two-species kinetic model of chemotaxis

In this paper, we propose a kinetic model describing the collective motion by chemotaxis of two species in interaction emitting the same chemoattractant. Such model can be seen as a generalisation to several species of the Othmer-Dunbar-Alt model which takes into account the run-and-tumble process of bacteria. Existence of weak solutions for this two-species kinetic model is studied and the convergence of its diffusive limit towards a macroscopic model of Keller-Segel type is analysed.

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

Numerical methods for one-dimensional aggregation equations

We focus in this work on the numerical discretization of the one dimensional aggregation equation $\pa_tρ+ \pa_x (vρ)=0$, $v=a(W'*ρ)$, in the attractive case. Finite time blow up of smooth initial data occurs for potential $W$ having a Lipschitz singularity at the origin. A numerical discretization is proposed for which the convergence towards duality solutions of the aggregation equation is proved. It relies on a careful choice of the discretized macroscopic velocity $v$ in order to give a sense to the product $v ρ$. Moreover, using the same idea, we propose an asymptotic preserving scheme for a kinetic system in hyperbolic scaling converging towards the aggregation equation in hydrodynamical limit. Finally numerical simulations are provided to illustrate the results.

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
preprint2014arXiv

Traveling wave solution of the Hele-Shaw model of tumor growth with nutrient

Several mathematical models of tumor growth are now commonly used to explain medical observations and predict cancer evolution based on images. These models incorporate mechanical laws for tissue compression combined with rules for nutrients availability which can differ depending on the situation under consideration, in vivo or in vitro. Numerical solutions exhibit, as expected from medical observations, a proliferative rim and a necrotic core. However, their precise profiles are rather complex, both in one and two dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical simulation of an energy-transport model for partially quantized particles

A coupled quantum-classical model describing the transport of electrons confined in nanoscale semiconductor devices is considered. Using the subband decomposition approach allows to separate the transport directions from the confinement direction. The motion of the gas in the transport direction is assumed to be classical. Then a hierarchy of adiabatic quantum-classical model is obtained, leading to subband SHE and energy-transport models, with explicit expression of the diffusion coefficients. The energy-transport-Schrödinger-Poisson model is then used for the numerical simulation of the transport of the electron gas in an ultra-scaled Double-Gate-MOSFET.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A remark on duality solutions for some weakly nonlinear scalar conservation laws

We investigate existence and uniqueness of duality solutions for a scalar conservation law with a nonlocal interaction kernel. Following the work of Bouchut and James (Comm. Partial Diff. Eq., 24, 1999), a notion of duality solution for such a nonlinear system is proposed, for which we do not have uniqueness. Then we prove that a natural definition of the flux allows to select a solution for which uniqueness holds.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analysis of a diffusive effective mass model for nanowires

We propose in this paper to derive and analyze a self-consistent model describing the diffusive transport in a nanowire. From a physical point of view, it describes the electron transport in an ultra-scaled confined structure, taking in account the interactions of charged particles with phonons. The transport direction is assumed to be large compared to the wire section and is described by a drift-diffusion equation including effective quantities computed from a Bloch problem in the crystal lattice. The electrostatic potential solves a Poisson equation where the particle density couples on each energy band a two dimensional confinement density with the monodimensional transport density given by the Boltzmann statistics. On the one hand, we study the derivation of this Nanowire Drift-Diffusion Poisson model from a kinetic level description. On the other hand, we present an existence result for this model in a bounded domain.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Chemotaxis: from kinetic equations to aggregate dynamics

The hydrodynamic limit for a kinetic model of chemotaxis is investigated. The limit equation is a non local conservation law, for which finite time blow-up occurs, giving rise to measure-valued solutions and discontinuous velocities. An adaptation of the notion of duality solutions, introduced for linear equations with discontinuous coefficients, leads to an existence result. Uniqueness is obtained through a precise definition of the nonlinear flux as well as the complete dynamics of aggregates, i.e. combinations of Dirac masses. Finally a particle method is used to build an adapted numerical scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the hydrodynamical limit for a one dimensional kinetic model of cell aggregation by chemotaxis

The hydrodynamic limit of a one dimensional kinetic model describing chemotaxis is investigated. The limit system is a conservation law coupled to an elliptic problem for which the macroscopic velocity is possibly discontinuous. Therefore, we need to work with measure-valued densities. After recalling a blow-up result in finite time of regular solutions for the hydrodynamic model, we establish a convergence result of the solutions of the kinetic model towards solutions of a problem limit defined thanks to the flux. Numerical simulations illustrate this convergence result.

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