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

Pierre Degond

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

ResearcherUnclaimed source record
math-phmath.MPmath.APmath.NAnlin.AOphysics.flu-dynphysics.comp-phmath.PRphysics.plasm-phphysics.soc-phPopulations and EvolutionCell Behaviorcond-mat.softmath.OCq-fin.GNcond-mat.stat-mechnlin.CGNumerical AnalysisQuantitative Methods

Catalog footprint

What is connected

74works
19topics
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

74 published item(s)

preprint2022arXiv

Bulk topological states in a new collective dynamics model

In this paper, we demonstrate the existence of topological states in a new collective dynamics model. This individual-based model (IBM) describes self-propelled rigid bodies moving with constant speed and adjusting their rigid-body attitude to that of their neighbors. In previous works, a macroscopic model has been derived from this IBM in a suitable scaling limit. In the present work, we exhibit explicit solutions of the macroscopic model characterized by a non-trivial topology. We show that these solutions are well approximated by the IBM during a certain time but then the IBM transitions towards topologically trivial states. Using a set of appropriately defined topological indicators, we reveal that the breakage of the non-trivial topology requires the system to go through a phase of maximal disorder. We also show that similar but topologically trivial initial conditions result in markedly different dynamics, suggesting that topology plays a key role in the dynamics of this system.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multi-species viscous models for tissue growth: incompressible limit and qualitative behaviour

We introduce two 2D mechanical models reproducing the evolution of two viscous tissues in contact. Their main property is to model the swirling cell motions while keeping the tissues segregated, as observed during vertebrate embryo elongation. Segregation is encoded differently in the two models: by passive or active segregation (based on a mechanical repulsion pressure). We compute the incompressible limits of the two models, and obtain strictly segregated solutions. The two models thus obtained are compared. A striking feature in the active segregation model is the persistence of the repulsion pressure at the limit: a ghost effect is discussed and confronted to the biological data. Thanks to a transmission problem formulation at the incompressible limit, we show a pressure jump at the tissues' boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Topological states and continuum model for swarmalators without force reciprocity

Swarmalators are systems of agents which are both self-propelled particles and oscillators. Each particle is endowed with a phase which modulates its interaction force with the other particles. In return, relative positions modulate phase synchronization between interacting particles. In the present model, there is no force reciprocity: when a particle attracts another one, the latter repels the former. This results in a pursuit behavior. In this paper, we derive a hydrodynamic model of this swarmalator system and show that it has explicit doubly-periodic travelling-wave solutions in two space dimensions. These special solutions enjoy non-trivial topology quantified by the index of the phase vector along a period in either dimension. Stability of these solutions is studied by investigating the conditions for hyperbolicity of the model. Numerical solutions of both the particle and hydrodynamic models are shown. They confirm the consistency of the hydrodynamic model with the particle one for small times or large phase-noise but also reveal the emergence of intriguing patterns in the case of small phase-noise.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A mean field game model of firm--level innovation

Knowledge spillovers occur when a firm researches a new technology and that technology is adapted or adopted by another firm, resulting in a social value of the technology that is larger than the initially predicted private value. As a result, firms systematically under--invest in research compared with the socially optimal investment strategy. Understanding the level of under--investment, as well as policies to correct it, is an area of active economic research. In this paper, we develop a new model of spillovers, taking inspiration from the available microeconomic data. We prove existence and uniqueness of solutions to the model, and we conduct some initial simulations to understand how indirect spillovers contribute to the productivity of a sector.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Fokker-Planck approach to the study of robustness in gene expression

We study several Fokker-Planck equations arising from a stochastic chemical kinetic system modeling a gene regulatory network in biology. The densities solving the Fokker-Planck equations describe the joint distribution of the messenger RNA and micro RNA content in a cell. We provide theoretical and numerical evidences that the robustness of the gene expression is increased in the presence of micro RNA. At the mathematical level, increased robustness shows in a smaller coefficient of variation of the marginal density of the messenger RNA in the presence of micro RNA. These results follow from explicit formulas for solutions. Moreover, thanks to dimensional analyses and numerical simulations we provide qualitative insight into the role of each parameter in the model. As the increase of gene expression level comes from the underlying stochasticity in the models, we eventually discuss the choice of noise in our models and its influence on our results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Large-scale dynamics of self-propelled particles moving through obstacles: model derivation and pattern formation

We model and study the patterns created through the interaction of collectively moving self-propelled particles (SPPs) and elastically tethered obstacles. Simulations of an individual-based model reveal at least three distinct large-scale patterns: travelling bands, trails and moving clusters. This motivates the derivation of a macroscopic partial differential equations model for the interactions between the self-propelled particles and the obstacles, for which we assume large tether stiffness. The result is a coupled system of non-linear, non-local partial differential equations. Linear stability analysis shows that patterning is expected if the interactions are strong enough and allows for the predictions of pattern size from model parameters. The macroscopic equations reveal that the obstacle interactions induce short-ranged SPP aggregation, irrespective of whether obstacles and SPPs are attractive or repulsive.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A Markov jump process modelling animal group size statistics

We translate a coagulation-framentation model, describing the dynamics of animal group size distributions, into a model for the population distribution and associate the \blue{nonlinear} evolution equation with a Markov jump process of a type introduced in classic work of H.~McKean. In particular this formalizes a model suggested by H.-S. Niwa [J.~Theo.~Biol.~224 (2003)] with simple coagulation and fragmentation rates. Based on the jump process, we develop a numerical scheme that allows us to approximate the equilibrium for the Niwa model, validated by comparison to analytical results by Degond et al. [J.~Nonlinear Sci.~27 (2017)], and study the population and size distributions for more complicated rates. Furthermore, the simulations are used to describe statistical properties of the underlying jump process. We additionally discuss the relation of the jump process to models expressed in stochastic differential equations and demonstrate that such a connection is justified in the case of nearest-neighbour interactions, as opposed to global interactions as in the Niwa model.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Phase transitions and macroscopic limits in a BGK model of body-attitude coordination

In this article we investigate the phase transition phenomena that occur in a model of self-organisation through body-attitude coordination. Here, the body-attitude of an agent is modelled by a rotation matrix in $\mathbb{R}^3$ as in [Degond, Frouvelle, Merino-Aceituno, 2017]. The starting point of this study is a BGK equation modelling the evolution of the distribution function of the system at a kinetic level. The main novelty of this work is to show that in the spatially homogeneous case, self-organisation may appear or not depending on the local density of agents involved. We first exhibit a connection between body-orientation models and models of nematic alignment of polymers in higher dimensional space from which we deduce the complete description of the possible equilibria Then, thanks to a gradient-flow structure specific to this BGK model, we are able to prove the stability and the convergence towards the equilibria in the different regimes. We then derive the macroscopic models associated to the stable equilibria in the spirit of [Degond, Frouvelle, Merino-Aceituno, 2017] and [Degond, Frouvelle, Liu, 2015].

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Finite Volume approximations of the Euler system with variable congestion

We are interested in the numerical simulations of the Euler system with variable congestion encoded by a singular pressure. This model describes for instance the macroscopic motion of a crowd with individual congestion preferences. We propose an asymptotic preserving (AP) scheme based on a conservative formulation of the system in terms of density, momentum and density fraction. A second order accuracy version of the scheme is also presented. We validate the scheme on one-dimensional test-cases and extended here to higher order accuracy. We finally carry out two dimensional numerical simulations and show that the model exhibit typical crowd dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A new flocking model through body attitude coordination

We present a new model for multi-agent dynamics where each agent is described by its position and body attitude: agents travel at a constant speed in a given direction and their body can rotate around it adopting different configurations. In this manner, the body attitude is described by three orthonormal axes giving an element in $SO(3)$ (rotation matrix). Agents try to coordinate their body attitudes with the ones of their neighbours. In the present paper, we give the Individual Based Model (particle model) for this dynamics and derive its corresponding kinetic and macroscopic equations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Continuum dynamics of the intention field under weakly cohesive social interactions

We investigate the long-time dynamics of an opinion formation model inspired by a work by Borghesi, Bouchaud and Jensen. Firstly, we derive a Fokker-Planck type equation under the assumption that interactions between individuals produce little consensus of opinion (grazing collision approximation). Secondly, we study conditions under which the Fokker-Planck equation has non-trivial equilibria and derive the macroscopic limit (corresponding to the long-time dynamics and spatially localized interactions) for the evolution of the mean opinion. Finally, we compare two different types of interaction rates: the original one given in the work of Borghesi, Bouchaud and Jensen (symmetric binary interactions) and one inspired from works by Motsch and Tadmor (non-symmetric binary interactions). We show that the first case leads to a conservative model for the density of the mean opinion whereas the second case leads to a non-conservative equation. We also show that the speed at which consensus is reached asymptotically for these two rates has fairly different density dependence.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Damped Arrow-Hurwicz algorithm for sphere packing

We consider algorithms that, from an arbitrarily sampling of $N$ spheres (possibly overlapping), find a close packed configuration without overlapping. These problems can be formulated as minimization problems with non-convex constraints. For such packing problems, we observe that the classical iterative Arrow-Hurwicz algorithm does not converge. We derive a novel algorithm from a multi-step variant of the Arrow-Hurwicz scheme with damping. We compare this algorithm with classical algorithms belonging to the class of linearly constrained Lagrangian methods and show that it performs better. We provide an analysis of the convergence of these algorithms in the simple case of two spheres in one spatial dimension. Finally, we investigate the behaviour of our algorithm when the number of spheres is large.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Kinetic theory of particle interactions mediated by dynamical networks

We provide a detailed multiscale analysis of a system of particles interacting through a dynamical network of links. Starting from a microscopic model, via the mean field limit, we formally derive coupled kinetic equations for the particle and link densities, following the approach of [Degond et al., M3AS, 2016]. Assuming that the process of remodelling the network is very fast, we simplify the description to a macroscopic model taking the form of single aggregation-diffusion equation for the density of particles. We analyze qualitatively this equation, addressing the stability of a homogeneous distribution of particles for a general potential. For the Hookean potential we obtain a precise condition for the phase transition, and, using the central manifold reduction, we characterize the type of bifurcation at the instability onset.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the asymptotic limit of the three dimensional Vlasov-Poisson system for large magnetic field : formal derivation

This paper establishes the long time asymptotic limit of the three dimensional Vlasov-Poisson equation with strong external magnetic field. The guiding center approximation is investigated in the three dimensional case with a non-constant magnetic field. In the long time asymptotic limit, the motion can be split in two parts : one stationary flow along the lines of the magnetic field and the guiding center motion in the orthogonal plane of the magnetic field where classical drift velocities are recovered. We discuss in particular the effect of nonconstant external magnetic fields.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Self-Organized Hydrodynamics with nonconstant velocity

Motivated by recent experimental and computational results that show a motility-induced clustering transition in self-propelled particle systems, we study an individual model and its corresponding Self-Organized Hydrodynamic model for collective behaviour that incorporates a density-dependent velocity, as well as inter-particle alignment. The modal analysis of the hydrodynamic model elucidates the relationship between the stability of the equilibria and the changing velocity, and the formation of clusters. We find, in agreement with earlier results for non-aligning particles, that the key criterion for stability is $(ρv(ρ))'> 0$, i.e. a non-rapid decrease of velocity with density. Numerical simulation for both the individual and hydrodynamic models with a velocity function inspired by experiment demonstrates the validity of the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A continuum model for nematic alignment of self-propelled particles

A continuum model for a population of self-propelled particles interacting through nematic alignment is derived from an individual-based model. The methodology consists of introducing a hydrodynamic scaling of the corresponding mean-field kinetic equation. The resulting perturbation problem is solved thanks to the concept of generalized collision invariants. It yields a hyperbolic but non-conservative system of equations for the nematic mean direction of the flow and the densities of particles flowing parallel or anti-parallel to this mean direction. Diffusive terms are introduced under a weakly non-local interaction assumption and the diffusion coefficient is proven to be positive. An application to the modeling of myxobacteria is outlined.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A multi-layer model for self-propelled disks interacting through alignment and volume exclusion

We present an individual-based model describing disk-like self-propelled particles moving inside parallel planes. The disk directions of motion follow alignment rules inside each layer. Additionally, the disks are subject to interactions with those of the neighboring layers arising from volume exclusion constraints. These interactions affect the disk inclinations with respect to the plane of motion. We formally de-rive a macroscopic model composed of planar Self-Organized Hydrodynamic (SOH) models describing the transport of mass and evolution of mean direction of motion of the disks in each plane, supplemented with transport equations for the mean disk inclination. These planar models are coupled due to the interactions with the neighboring planes. Numerical comparisons between the individual-based and macroscopic models are carried out. These models could be applicable, for instance, to describe sperm-cell collective dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An Asymptotic-Preserving Method for a Relaxation of the Navier-Stokes-Korteweg Equations

The Navier-Stokes-Korteweg (NSK) equations are a classical diffuse-interface model for compressible two-phase flow. As direct numerical simulations based on the NSK system are quite expensive and in some cases even impossible, we consider a relaxation of the NSK system, for which robust numerical methods can be designed. However, time steps for explicit numerical schemes depend on the relaxation parameter and therefore numerical simulations in the relaxation limit are very inefficient. To overcome this restriction, we propose an implicit-explicit asymptotic-preserving finite volume method. We prove that the new scheme provides a consistent discretization of the NSK system in the relaxation limit and demonstrate that it is capable of accurately and efficiently computing numerical solutions of problems with realistic density ratios and small interfacial widths.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Asymptotic-preserving Particle-In-Cell methods for the Vlasov-Maxwell system near quasi-neutrality

In this article, we design Asymptotic-Preserving Particle-In-Cell methods for the Vlasov-Maxwell system in the quasi-neutral limit, this limit being characterized by a Debye length negligible compared to the space scale of the problem. These methods are consistent discretizations of the Vlasov-Maxwell system which, in the quasi-neutral limit, remain stable and are consistent with a quasi-neutral model (in this quasi-neutral model, the electric field is computed by means of a generalized Ohm law). The derivation of Asymptotic-Preserving methods is not straightforward since the quasi-neutral model is a singular limit of the Vlasov-Maxwell model. The key step is a reformulation of the Vlasov-Maxwell system which unifies the two models in a single set of equations with a smooth transition from one to another. As demonstrated in various and demanding numerical simulations, the Asymptotic-Preserving methods are able to treat efficiently both quasi-neutral plasmas and non-neutral plasmas, making them particularly well suited for complex problems involving dense plasmas with localized non-neutral regions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coagulation-fragmentation model for animal group-size statistics

We study coagulation-fragmentation equations inspired by a simple model proposed in fisheries science to explain data for the size distribution of schools of pelagic fish. Although the equations lack detailed balance and admit no $H$-theorem, we are able to develop a rather complete description of equilibrium profiles and large-time behavior, based on recent developments in complex function theory for Bernstein and Pick functions. In the large-population continuum limit, a scaling-invariant regime is reached in which all equilibria are determined by a single scaling profile. This universal profile exhibits power-law behavior crossing over from exponent $-\frac23$ for small size to $-\frac32$ for large size, with an exponential cut-off.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Continuum model for linked fibers with alignment interactions

We introduce an individual-based model for fiber elements having the ability to cross-link or unlink each other and to align with each other at the cross links. We first formally derive a kinetic model for the fiber and cross-links distribution functions. We then consider the fast linking/unlinking regime in which the model can be reduced to the fiber distribution function only and investigate its diffusion limit. The resulting macroscopic model consists of a system of nonlinear diffusion equations for the fiber density and mean orientation. In the case of a homogeneous fiber density, we show that the model is elliptic.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Topological interactions in a Boltzmann-type framework

We consider a finite number of particles characterised by their positions and velocities. At random times a randomly chosen particle, the follower, adopts the velocity of another particle, the leader. The follower chooses its leader according to the proximity rank of the latter with respect to the former. We study the limit of a system size going to infinity and, under the assumption of propagation of chaos, show that the limit equation is akin to the Boltzmann equation. However, it exhibits a spatial non-locality instead of the classical non-locality in velocity space. This result relies on the approximation properties of Bernstein polynomials.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Boltzmann model for rod alignment and schooling fish

We consider a Boltzmann model introduced by Bertin, Droz and Greegoire as a binary interaction model of the Vicsek alignment interaction. This model considers particles lying on the circle. Pairs of particles interact by trying to reach their mid-point (on the circle) up to some noise. We study the equilibria of this Boltzmann model and we rigorously show the existence of a pitchfork bifurcation when a parameter measuring the inverse of the noise intensity crosses a critical threshold. The analysis is carried over rigorously when there are only finitely many non-zero Fourier modes of the noise distribution. In this case, we can show that the critical exponent of the bifurcation is exactly 1/2. In the case of an infinite number of non-zero Fourier modes, a similar behavior can be formally obtained thanks to a method relying on integer partitions first proposed by Ben-Naim and Krapivsky.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A hierarchy of models related to nanoflows and surface diffusion

In last years a great interest was brought to molecular transport problems at nanoscales, such as surface diffusion or molecular flows in nano or sub-nano-channels. In a series of papers V. D. Borman, S. Y. Krylov, A. V. Prosyanov and J. J. M. Beenakker proposed to use kinetic theory in order to analyze the mechanisms that determine mobility of molecules in nanoscale channels. This approach proved to be remarkably useful to give new insight on these issues, such as density dependence of the diffusion coefficient. In this paper we revisit these works to derive the kinetic and diffusion models introduced by V. D. Borman, S. Y. Krylov, A. V. Prosyanov and J. J. M. Beenakker by using classical tools of kinetic theory such as scaling and systematic asymptotic analysis. Some results are extended to less restrictive hypothesis.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Evolution of wealth in a nonconservative economy driven by local Nash equilibria

We develop a model for the evolution of wealth in a non-conservative economic environment, extending a theory developed earlier by the authors. The model considers a system of rational agents interacting in a game theoretical framework. This evolution drives the dynamic of the agents in both wealth and economic configuration variables. The cost function is chosen to represent a risk averse strategy of each agent. That is, the agent is more likely to interact with the market, the more predictable the market, and therefore the smaller its individual risk. This yields a kinetic equation for an effective single particle agent density with a Nash equilibrium serving as the local thermodynamic equilibrium. We consider a regime of scale separation where the large scale dynamics is given by a hydrodynamic closure with this local equilibrium. A class of generalized collision invariants (GCIs) is developed to overcome the difficulty of the non-conservative property in the hydrodynamic closure derivation of the large scale dynamics for the evolution of wealth distribution. The result is a system of gas dynamics-type equations for the density and average wealth of the agents on large scales. We recover the inverse Gamma distribution, which has been previously considered in the literature, as a local equilibrium for particular choices of the cost function.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Local stability of perfect alignment for a spatially homogeneous kinetic model

We prove the nonlinear local stability of Dirac masses for a kinetic model of alignment of particles on the unit sphere, each point of the unit sphere representing a direction. A population concentrated in a Dirac mass then corresponds to the global alignment of all individuals. The main difficulty of this model is the lack of conserved quantities and the absence of an energy that would decrease for any initial condition. We overcome this difficulty thanks to a functional which is decreasing in time in a neighborhood of any Dirac mass (in the sense of the Wasserstein distance). The results are then extended to the case where the unit sphere is replaced by a general Riemannian manifold.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Macroscopic models of collective motion with repulsion

We study a system of self-propelled particles which interact with their neighbors via alignment and repulsion. The particle velocities result from self-propulsion and repulsion by close neighbors. The direction of self-propulsion is continuously aligned to that of the neighbors, up to some noise. A continuum model is derived starting from a mean-field kinetic description of the particle system. It leads to a set of non conservative hydrodynamic equations. We provide a numerical validation of the continuum model by comparison with the particle model. We also provide comparisons with other self-propelled particle models with alignment and repulsion.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Meanfield games and model predictive control

Mean-Field Games are games with a continuum of players that incorporate the time-dimension through a control-theoretic approach. Recently, simpler approaches relying on the Best Reply Strategy have been proposed. They assume that the agents navigate their strategies towards their goal by taking the direction of steepest descent of their cost function (i.e. the opposite of the utility function). In this paper, we explore the link between Mean-Field Games and the Best Reply Strategy approach. This is done by introducing a Model Predictive Control framework, which consists of setting the Mean-Field Game over a short time interval which recedes as time moves on. We show that the Model Predictive Control offers a compromise between a possibly unrealistic Mean-Field Game approach and the sub-optimal Best Reply Strategy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-Conforming Multiscale Finite Element Method for Stokes Flows in Heterogeneous Media. Part I: Methodologies and Numerical Experiments

The Multiscale Finite Element Method (MsFEM) is developed in the vein of Crouzeix-Raviart element for solving viscous incompressible flows in genuine heterogeneous media. Such flows are relevant in many branches of engineering, often at multiple scales and at regions where analytical representations of the microscopic features of the flows are often unavailable. Full accounts to these problems heavily depend on the geometry of the system under consideration and are computationally expensive. Therefore, a method capable of solving multiscale features of the flow without confining itself to fine scale calculations is sought after. The approximation of boundary condition on coarse element edges when computing the multiscale basis functions critically influences the eventual accuracy of any MsFEM approaches. The weakly enforced continuity of Crouzeix - Raviart function space across element edges leads to a natural boundary condition for the multiscale basis functions which relaxes the sensitivity of our method to complex patterns of obstacles exempt from the needs of implementing any oversampling techniques. Additionally, the application of penalization method makes it possible to avoid complex unstructured domain and allows extensive use of simpler Cartesian meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase transitions, hysteresis, and hyperbolicity for self-organized alignment dynamics

We provide a complete and rigorous description of phase transitions for kinetic models of self-propelled particles interacting through alignment. These models exhibit a competition between alignment and noise. Both the alignment frequency and noise intensity depend on a measure of the local alignment. We show that, in the spatially homogeneous case, the phase transition features (number and nature of equilibria, stability, convergence rate, phase diagram, hysteresis) are totally encoded in how the ratio between the alignment and noise intensities depend on the local alignment. In the spatially inhomogeneous case, we derive the macroscopic models associated to the stable equilibria and classify their hyperbolicity according to the same function.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-organized Hydrodynamics in an Annular Domain: Modal Analysis and Nonlinear Effects

The Self-Organized Hydrodynamics model of collective behavior is studied on an annular domain. A modal analysis of the linearized model around a perfectly polarized steady-state is conducted. It shows that the model has only pure imaginary modes in countable number and is hence stable. Numerical computations of the low-order modes are provided. The fully non-linear model is numerically solved and nonlinear mode-coupling is then analyzed. Finally, the efficiency of the modal decomposition to analyze the complex features of the nonlinear model is demonstrated.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Crouzeix-Raviart MsFEM with Bubble Functions for Diffusion and Advection-Diffusion in Perforated Media

The adaptation of Crouzeix - Raviart finite element in the context of multiscale finite element method (MsFEM) is studied and implemented on diffusion and advection-diffusion problems in perforated media. It is known that the approximation of boundary condition on coarse element edges when computing the multiscale basis functions critically influences the eventual accuracy of any MsFEM approaches. The weakly enforced continuity of Crouzeix - Raviart function space across element edges leads to a natural boundary condition for the multiscale basis functions which relaxes the sensitivity of our method to complex patterns of perforations. Another ingredient to our method is the application of bubble functions which is shown to be instrumental in maintaining high accuracy amid dense perforations. Additionally, the application of penalization method makes it possible to avoid complex unstructured domain and allows extensive use of simpler Cartesian meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evolution of the distribution of wealth in an economic environment driven by local Nash equilibria

We present and analyze a model for the evolution of the wealth distribution within a heterogeneous economic environment. The model considers a system of rational agents interacting in a game theoretical framework, through fairly general assumptions on the cost function. This evolution drives the dynamic of the agents in both wealth and economic configuration variables. We consider a regime of scale separation where the large scale dynamics is given by a hydrodynamic closure with a Nash equilibrium serving as the local thermodynamic equilibrium. The result is a system of gas dynamics-type equations for the density and average wealth of the agents on large scales. We recover the inverse gamma distribution as an equilibrium in the particular case of quadratic cost functions which has been previously considered in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Flow on sweeping networks

We introduce a cellular automaton model coupled with a transport equation for flows on graphs. The direction of the flow is described by a switching process where the switching probability dynamically changes according to the value of the transported quantity in the neighboring cells. A motivation is pedestrian dynamics in a small corridor where the propagation of people in a part of the corridor can be either left or rightgoing. Under the assumptions of propagation of chaos and mean-field limit, we derive a master equation and the corresponding meanfield kinetic and macroscopic models. Steady--states are computed and analyzed analytically and exhibit the possibility of multiple meta-stable states and hysteresis.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hydrodynamics of the Kuramoto-Vicsek model of rotating self-propelled particles

We consider an Individual-Based Model for self-rotating particles interacting through local alignment and investigate its macroscopic limit. This model describes self-propelled particles moving in the plane and trying to synchronize their rotation motion with their neighbors. It combines the Kuramoto model of synchronization and the Vicsek model of swarm formation. We study the mean-field kinetic and hydrodynamic limits of this system within two different scalings. In the small angular velocity regime, the resulting model is a slight modification of the 'Self-Organized Hydrodynamic' model which has been previously introduced by the first author. In the large angular velocity case, a new type of hydrodynamic model is obtained. A preliminary study of the linearized stability is proposed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Macroscopic models of collective motion and self-organization

In this paper, we review recent developments on the derivation and properties of macroscopic models of collective motion and self-organization. The starting point is a model of self-propelled particles interacting with its neighbors through alignment. We successively derive a mean-field model and its hydrodynamic limit. The resulting macroscopic model is the Self-Organized Hydrodynamics (SOH). We review the available existence results and known properties of the SOH model and discuss it in view of its possible extensions to other kinds of collective motion.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Nash equilibrium macroscopic closure for kinetic models coupled with Mean-Field Games

We introduce a new mean field kinetic model for systems of rational agents interacting in a game theoretical framework. This model is inspired from non-cooperative anonymous games with a continuum of players and Mean-Field Games. The large time behavior of the system is given by a macroscopic closure with a Nash equilibrium serving as the local thermodynamic equilibrium. An application of the presented theory to a social model (herding behavior) is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A note on phase transitions for the Smoluchowski equation with dipolar potential

In this note, we study the phase transitions arising in a modified Smoluchowski equation on the sphere with dipolar potential. This equation models the competition between alignment and diffusion, and the modification consists in taking the strength of alignment and the intensity of the diffusion as functions of the order parameter. We characterize the stable and unstable equilibrium states. For stable equilibria, we provide the exponential rate of convergence. We detail special cases, giving rise to second order and first order phase transitions, respectively. We study the hysteresis diagram, and provide numerical illustrations of this phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phase transition and diffusion among socially interacting self-propelled agents

We consider a hydrodynamic model of swarming behavior derived from the kinetic description of a particle system combining a noisy Cucker-Smale consensus force and self-propulsion. In the large self-propulsion force limit, we provide evidence of a phase transition from disordered to ordered motion which manifests itself as a change of type of the limit model (from hyperbolic to diffusive) at the crossing of a critical noise intensity. In the hyperbolic regime, the resulting model, referred to as the `Self-Organized Hydrodynamics (SOH)', consists of a system of compressible Euler equations with a speed constraint. We show that the range of SOH models obtained by this limit is restricted. To waive this restriction, we compute the Navier-Stokes diffusive corrections to the hydrodynamic model. Adding these diffusive corrections, the limit of a large propulsion force yields unrestricted SOH models and offers an alternative to the derivation of the SOH using kinetic models with speed constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Self-Organized Hydrodynamics with congestion and path formation in crowds

A continuum model for self-organized dynamics is numerically investigated. The model describes systems of particles subject to alignment interaction and short-range repulsion. It consists of a non-conservative hyperbolic system for the density and velocity orientation. Short-range repulsion is included through a singular pressure which becomes infinite at the jamming density. The singular limit of infinite pressure stiffness leads to phase transitions from compressible to incompressible dynamics. The paper proposes an Asymptotic-Preserving scheme which takes care of the singular pressure while preventing the breakdown of the CFL stability condition near congestion. It relies on a relaxation approximation of the system and an elliptic formulation of the pressure equation. Numerical simulations of impinging clusters show the efficiency of the scheme to treat congestions. A two-fluid variant of the model provides a model of path formation in crowds.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Traffic Instabilities in Self-Organized Pedestrian Crowds

In human crowds as well as in many animal societies, local interactions among individuals often give rise to self-organized collective organizations that offer functional benefits to the group. For instance, flows of pedestrians moving in opposite directions spontaneously segregate into lanes of uniform walking directions. This phenomenon is often referred to as a smart collective pattern, as it increases the traffic efficiency with no need of external control. However, the functional benefits of this emergent organization have never been experimentally measured, and the underlying behavioral mechanisms are poorly understood. In this work, we have studied this phenomenon under controlled laboratory conditions. We found that the traffic segregation exhibits structural instabilities characterized by the alternation of organized and disorganized states, where the lifetime of well-organized clusters of pedestrians follow a stretched exponential relaxation process. Further analysis show that the inter-pedestrian variability of comfortable walking speeds is a key variable at the origin of the observed traffic perturbations. We show that the collective benefit of the emerging pattern is maximized when all pedestrians walk at the average speed of the group. In practice, however, local interactions between slow- and fast-walking pedestrians trigger global breakdowns of organization, which reduce the collective and the individual payoff provided by the traffic segregation. This work is a step ahead toward the understanding of traffic self-organization in crowds, which turns out to be modulated by complex behavioral mechanisms that do not always maximize the group's benefits. The quantitative understanding of crowd behaviors opens the way for designing bottom-up management strategies bound to promote the emergence of efficient collective behaviors in crowds.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An Asymptotic-Preserving all-speed scheme for the Euler and Navier-Stokes equations

We present an Asymptotic-Preserving 'all-speed' scheme for the simulation of compressible flows valid at all Mach-numbers ranging from very small to order unity. The scheme is based on a semi-implicit discretization which treats the acoustic part implicitly and the convective and diffusive parts explicitly. This discretization, which is the key to the Asymptotic-Preserving property, provides a consistent approximation of both the hyperbolic compressible regime and the elliptic incompressible regime. The divergence-free condition on the velocity in the incompressible regime is respected, and an the pressure is computed via an elliptic equation resulting from a suitable combination of the momentum and energy equations. The implicit treatment of the acoustic part allows the time-step to be independent of the Mach number. The scheme is conservative and applies to steady or unsteady flows and to general equations of state. One and Two-dimensional numerical results provide a validation of the Asymptotic-Preserving 'all-speed' properties.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An Asymptotic-Preserving method for highly anisotropic elliptic equations based on a micro-macro decomposition

The concern of the present work is the introduction of a very efficient Asymptotic Preserving scheme for the resolution of highly anisotropic diffusion equations. The characteristic features of this scheme are the uniform convergence with respect to the anisotropy parameter $0<\eps <<1$, the applicability (on cartesian grids) to cases of non-uniform and non-aligned anisotropy fields $b$ and the simple extension to the case of a non-constant anisotropy intensity $1/\eps$. The mathematical approach and the numerical scheme are different from those presented in the previous work [Degond et al. (2010), arXiv:1008.3405v1] and its considerable advantages are pointed out.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Asymptotic-Preserving scheme for a bi-fluid Euler-Lorentz model

The present work is devoted to the simulation of a strongly magnetized plasma considered as a mixture of an ion fluid and an electron fluid. For the sake of simplicity, we assume that the model is isothermal and described by Euler equations coupled with a term representing the Lorentz force. Moreover we assume that both Euler systems are coupled through a quasi-neutrality constraint. The numerical method which is described in the present document is based on an Asymptotic-Preserving semi-discretization in time of a variant of this two-fluid Euler-Lorentz model with a small perturbation of the quasi-neutrality constraint. Firstly, we present the two-fluid model and the motivations for introducing a small perturbation into the quasi-neutrality equation, then we describe the time semi-discretization of the perturbed model and a fully-discrete finite volume scheme based on it. Finally, we present some numerical results which have been obtained with this method.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Asymptotic-Preserving Schemes for Fluid Models of Plasmas

These notes summarize a series of works related to the numerical approximation of plasma fluid problems. We construct so-called 'Asymptotic-Preserving' schemes which are valid for a large range of values (from very small to order unity) of the dimensionless parameters that appear in plasma fluid models. Specifically, we are interested in two parameters, the scaled Debye length which quantifies how close to quasi-neutrality the plasma is, and the scaled cyclotron period, which is inversely proportional to the magnetic field strength. We will largely focus on the ideas, in order to enable the reader to apply these concepts to other situations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hydrodynamic models of self-organized dynamics: derivation and existence theory

This paper is concerned with the derivation and analysis of hydrodynamic models for systems of self-propelled particles subject to alignment interaction and attraction-repulsion. The starting point is the kinetic model considered in earlier work of Degond & Motsch with the addition of an attraction-repulsion interaction potential. Introducing different scalings than in Degond & Motsch, the non-local effects of the alignment and attraction-repulsion interactions can be kept in the hydrodynamic limit and result in extra pressure, viscosity terms and capillary force. The systems are shown to be symmetrizable hyperbolic systems with viscosity terms. A local-in-time existence result is proved in the 2D case for the viscous model and in the 3D case for the inviscid model. The proof relies on the energy method.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hydrodynamics of self-alignment interactions with precession and derivation of the Landau-Lifschitz-Gilbert equation

We consider a kinetic model of self-propelled particles with alignment interaction and with precession about the alignment direction. We derive a hydrodynamic system for the local density and velocity orientation of the particles. The system consists of the conservative equation for the local density and a non-conservative equation for the orientation. First, we assume that the alignment interaction is purely local and derive a first order system. However, we show that this system may lose its hyperbolicity. Under the assumption of weakly non-local interaction, we derive diffusive corrections to the first order system which lead to the combination of a heat flow of the harmonic map and Landau-Lifschitz-Gilbert dynamics. In the particular case of zero self-propelling speed, the resulting model reduces to the phenomenological Landau-Lifschitz-Gilbert equations. Therefore the present theory provides a kinetic formulation of classical micromagnetization models and spin dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kinetic formulation and global existence for the Hall-Magneto-hydrodynamics system

This paper deals with the derivation and analysis of the the Hall Magneto-Hydrodynamic equations. We first provide a derivation of this system from a two-fluids Euler-Maxwell system for electrons and ions, through a set of scaling limits. We also propose a kinetic formulation for the Hall-MHD equations which contains as fluid closure different variants of the Hall-MHD model. Then, we prove the existence of global weak solutions for the incompressible viscous resistive Hall-MHD model. We use the particular structure of the Hall term which has zero contribution to the energy identity. Finally, we discuss particular solutions in the form of axisymmetric purely swirling magnetic fields and propose some regularization of the Hall equation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kinetic hierarchy and propagation of chaos in biological swarm models

We consider two models of biological swarm behavior. In these models, pairs of particles interact to adjust their velocities one to each other. In the first process, called 'BDG', they join their average velocity up to some noise. In the second process, called 'CL', one of the two particles tries to join the other one's velocity. This paper establishes the master equations and BBGKY hierarchies of these two processes. It investigates the infinite particle limit of the hierarchies at large time-scale. It shows that the resulting kinetic hierarchy for the CL process does not satisfy propagation of chaos. Numerical simulations indicate that the BDG process has similar behavior to the CL process.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kinetic limits for pair-interaction driven master equations and biological swarm models

We consider a class of stochastic processes modeling binary interactions in an N-particle system. Examples of such systems can be found in the modeling of biological swarms. They lead to the definition of a class of master equations that we call pair interaction driven master equations. We prove a propagation of chaos result for this class of master equations which generalizes Mark Kac's well know result for the Kac model in kinetic theory. We use this result to study kinetic limits for two biological swarm models. We show that propagation of chaos may be lost at large times and we exhibit an example where the invariant density is not chaotic.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Macroscopic limits and phase transition in a system of self-propelled particles

We investigate systems of self-propelled particles with alignment interaction. Compared to previous work, the force acting on the particles is not normalized and this modification gives rise to phase transitions from disordered states at low density to aligned states at high densities. This model is the space inhomogeneous extension of a previous work by Frouvelle and Liu in which the existence and stability of the equilibrium states were investigated. When the density is lower than a threshold value, the dynamics is described by a non-linear diffusion equation. By contrast, when the density is larger than this threshold value, the dynamics is described by a hydrodynamic model for self-alignment interactions previously derived in Degond and Motsch. However, the modified normalization of the force gives rise to different convection speeds and the resulting model may lose its hyperbolicity in some regions of the state space.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical approximation of the Euler-Maxwell model in the quasineutral limit

We derive and analyze an Asymptotic-Preserving scheme for the Euler-Maxwell system in the quasi-neutral limit. We prove that the linear stability condition on the time-step is independent of the scaled Debye length $λ$ when $λ\to 0$. Numerical validation performed on Riemann initial data and for a model Plasma Opening Switch device show that the AP-scheme is convergent to the Euler-Maxwell solution when $Δx/ λ\to 0$ where $Δx$ is the spatial discretization. But, when $λ/Δx \to 0$, the AP-scheme is consistent with the quasi-neutral Euler-Maxwell system. The scheme is also perfectly consistent with the Gauss equation. The possibility of using large time and space steps leads to several orders of magnitude reductions in computer time and storage.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Trail formation based on directed pheromone deposition

We propose an Individual-Based Model of ant-trail formation. The ants are modeled as self-propelled particles which deposit directed pheromones and interact with them through alignment interaction. The directed pheromones intend to model pieces of trails, while the alignment interaction translates the tendency for an ant to follow a trail when it meets it. Thanks to adequate quantitative descriptors of the trail patterns, the existence of a phase transition as the ant-pheromone interaction frequency is increased can be evidenced. Finally, we propose both kinetic and fluid descriptions of this model and analyze the capabilities of the fluid model to develop trail patterns. We observe that the development of patterns by fluid models require extra trail amplification mechanisms that are not needed at the Individual-Based Model level.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Two-way multi-lane traffic model for pedestrians in corridors

We extend the Aw-Rascle macroscopic model of car traffic into a two-way multi-lane model of pedestrian traffic. Within this model, we propose a technique for the handling of the congestion constraint, i.e. the fact that the pedestrian density cannot exceed a maximal density corresponding to contact between pedestrians. In a first step, we propose a singularly perturbed pressure relation which models the fact that the pedestrian velocity is considerably reduced, if not blocked, at congestion. In a second step, we carry over the singular limit into the model and show that abrupt transitions between compressible flow (in the uncongested regions) to incompressible flow (in congested regions) occur. We also investigate the hyperbolicity of the two-way models and show that they can lose their hyperbolicity in some cases. We study a diffusive correction of these models and discuss the characteristic time and length scales of the instability.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A hierarchy of models related to nanoflows and surface diffusion

In last years a great interest was brought to molecular transport problems at nanoscales, such as surface diffusion or molecular flows in nano or sub-nano-channels. In a series of papers V. D. Borman, S. Y. Krylov, A. V. Prosyanov and J. J. M. Beenakker proposed to use kinetic theory in order to analyze the mechanisms that determine mobility of molecules in nanoscale channels. This approach proved to be remarkably useful to give new insight on these issues, such as density dependence of the diffusion coefficient. In this paper we revisit these works to derive the kinetic and diffusion models introduced by V. D. Borman, S. Y. Krylov, A. V. Prosyanov and J. J. M. Beenakker by using classical tools of kinetic theory such as scaling and systematic asymptotic analysis. Some results are extended to less restrictive hypothesis.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A macroscopic model for a system of swarming agents using curvature control

In this paper, we study the macroscopic limit of a new model of collective displacement. The model, called PTWA, is a combination of the Vicsek alignment model and the Persistent Turning Walker (PTW) model of motion by curvature control. The PTW model was designed to fit measured trajectories of individual fish. The PTWA model (Persistent Turning Walker with Alignment) describes the displacements of agents which modify their curvature in order to align with their neighbors. The derivation of its macroscopic limit uses the non-classical notion of generalized collisional invariant. The macroscopic limit of the PTWA model involves two physical quantities, the density and the mean velocity of individuals. It is a system of hyperbolic type but is non-conservative due to a geometric constraint on the velocity. This system has the same form as the macroscopic limit of the Vicsek model (the 'Vicsek hydrodynamics') but for the expression of the model coefficients. The numerical computations show that the numerical values of the coefficients are very close. The 'Vicsek Hydrodynamic model' appears in this way as a more generic macroscopic model of swarming behavior as originally anticipated.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Degenerate anisotropic elliptic problems and magnetized plasma simulations

This paper is devoted to the numerical approximation of a degenerate anisotropic elliptic problem. The numerical method is designed for arbitrary space-dependent anisotropy directions and does not require any specially adapted coordinate system. It is also designed to be equally accurate in the strongly and the mildly anisotropic cases. The method is applied to the Euler-Lorentz system, in the drift-fluid limit. This system provides a model for magnetized plasmas.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Diffusion in a continuum model of self-propelled particles with alignment interaction

In this paper, we provide the $O(ε)$ corrections to the hydrodynamic model derived by Degond and Motsch from a kinetic version of the model by Vicsek & coauthors describing flocking biological agents. The parameter $ε$ stands for the ratio of the microscopic to the macroscopic scales. The $O(ε)$ corrected model involves diffusion terms in both the mass and velocity equations as well as terms which are quadratic functions of the first order derivatives of the density and velocity. The derivation method is based on the standard Chapman-Enskog theory, but is significantly more complex than usual due to both the non-isotropy of the fluid and the lack of momentum conservation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Duality-based Asymptotic-Preserving method for highly anisotropic diffusion equations

The present paper introduces an efficient and accurate numerical scheme for the solution of a highly anisotropic elliptic equation, the anisotropy direction being given by a variable vector field. This scheme is based on an asymptotic preserving reformulation of the original system, permitting an accurate resolution independently of the anisotropy strength and without the need of a mesh adapted to this anisotropy. The counterpart of this original procedure is the larger system size, enlarged by adding auxiliary variables and Lagrange multipliers. This Asymptotic-Preserving method generalizes the method investigated in a previous paper [arXiv:0903.4984v2] to the case of an arbitrary anisotropy direction field.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fluid Simulations with Localized Boltzmann Upscaling by Direct Simulation Monte-Carlo

In the present work, we present a novel numerical algorithm to couple the Direct Simulation Monte Carlo method (DSMC) for the solution of the Boltzmann equation with a finite volume like method for the solution of the Euler equations. Recently we presented in [14],[16],[17] different methodologies which permit to solve fluid dynamics problems with localized regions of departure from thermodynamical equilibrium. The methods rely on the introduction of buffer zones which realize a smooth transition between the kinetic and the fluid regions. In this paper we extend the idea of buffer zones and dynamic coupling to the case of the Monte Carlo methods. To facilitate the coupling and avoid the onset of spurious oscillations in the fluid regions which are consequences of the coupling with a stochastic numerical scheme, we use a new technique which permits to reduce the variance of the particle methods [11]. In addition, the use of this method permits to obtain estimations of the breakdowns of the fluid models less affected by fluctuations and consequently to reduce the kinetic regions and optimize the coupling. In the last part of the paper several numerical examples are presented to validate the method and measure its computational performances.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Kinetic models for dilute solutions of dumbbells in non-homogeneous flows revisited

We propose a two fluid theory to model a dilute polymer solution assuming that it consists of two phases, polymer and solvent, with two distinct macroscopic velocities. The solvent phase velocity is governed by the macroscopic Navier-Stokes equations with the addition of a force term describing the interaction between the two phases. The polymer phase is described on the mesoscopic level using a dumbbell model and its macroscopic velocity is obtained through averaging. We start by writing down the full phase-space distribution function for the dumbbells and then obtain the inertialess limits for the Fokker-Planck equation and for the averaged friction force acting between the phases from a rigorous asymptotic analysis. The resulting equations are relevant to the modelling of strongly non-homogeneous flows, while the standard kinetic model is recovered in the locally homogeneous case.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical approximation of the Euler-Poisson-Boltzmann model in the quasineutral limit

This paper analyzes various schemes for the Euler-Poisson-Boltzmann (EPB) model of plasma physics. This model consists of the pressureless gas dynamics equations coupled with the Poisson equation and where the Boltzmann relation relates the potential to the electron density. If the quasi-neutral assumption is made, the Poisson equation is replaced by the constraint of zero local charge and the model reduces to the Isothermal Compressible Euler (ICE) model. We compare a numerical strategy based on the EPB model to a strategy using a reformulation (called REPB formulation). The REPB scheme captures the quasi-neutral limit more accurately.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical simulations of the Euler system with congestion constraint

In this paper, we study the numerical simulations for Euler system with maximal density constraint. This model is developed in [1, 3] with the constraint introduced into the system by a singular pressure law, which causes the transition of different asymptotic dynamics between different regions. To overcome these difficulties, we adapt and implement two asymptotic preserving (AP) schemes originally designed for low Mach number limit [2,4] to our model. These schemes work for the different dynamics and capture the transitions well. Several numerical tests both in one dimensional and two dimensional cases are carried out for our schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Multiscale Kinetic-Fluid Solver with Dynamic Localization of Kinetic Effects

This paper collects the efforts done in our previous works [P. Degond, S. Jin, L. Mieussens, A Smooth Transition Between Kinetic and Hydrodynamic Equations, J. Comp. Phys., 209 (2005) 665--694.],[P.Degond, G. Dimarco, L. Mieussens, A Moving Interface Method for Dynamic Kinetic-fluid Coupling, J. Comp. Phys., Vol. 227, pp. 1176-1208, (2007).],[P. Degond, J.G. Liu, L. Mieussens, Macroscopic Fluid Model with Localized Kinetic Upscaling Effects, SIAM Multi. Model. Sim. 5(3), 940--979 (2006)] to build a robust multiscale kinetic-fluid solver. Our scope is to efficiently solve fluid dynamic problems which present non equilibrium localized regions that can move, merge, appear or disappear in time. The main ingredients of the present work are the followings ones: a fluid model is solved in the whole domain together with a localized kinetic upscaling term that corrects the fluid model wherever it is necessary; this multiscale description of the flow is obtained by using a micro-macro decomposition of the distribution function [P. Degond, J.G. Liu, L. Mieussens, Macroscopic Fluid Model with Localized Kinetic Upscaling Effects, SIAM Multi. Model. Sim. 5(3), 940--979 (2006)]; the dynamic transition between fluid and kinetic descriptions is obtained by using a time and space dependent transition function; to efficiently define the breakdown conditions of fluid models we propose a new criterion based on the distribution function itself. Several numerical examples are presented to validate the method and measure its computational efficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

All speed scheme for the low mach number limit of the Isentropic Euler equation

An all speed scheme for the Isentropic Euler equation is presented in this paper. When the Mach number tends to zero, the compressible Euler equation converges to its incompressible counterpart, in which the density becomes a constant. Increasing approximation errors and severe stability constraints are the main difficulty in the low Mach regime. The key idea of our all speed scheme is the special semi-implicit time discretization, in which the low Mach number stiff term is divided into two parts, one being treated explicitly and the other one implicitly. Moreover, the flux of the density equation is also treated implicitly and an elliptic type equation is derived to obtain the density. In this way, the correct limit can be captured without requesting the mesh size and time step to be smaller than the Mach number. Compared with previous semi-implicit methods, nonphysical oscillations can be suppressed. We develop this semi-implicit time discretization in the framework of a first order local Lax-Friedrich (LLF) scheme and numerical tests are displayed to demonstrate its performances.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

An asymptotic preserving scheme for strongly anisotropic elliptic problems

In this article we introduce an asymptotic preserving scheme designed to compute the solution of a two dimensional elliptic equation presenting large anisotropies. We focus on an anisotropy aligned with one direction, the dominant part of the elliptic operator being supplemented with Neumann boundary conditions. A new scheme is introduced which allows an accurate resolution of this elliptic equation for an arbitrary anisotropy ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Congestion in a macroscopic model of self-driven particles modeling gregariousness

We analyze a macroscopic model with a maximal density constraint which describes short range repulsion in biological systems. This system aims at modeling finite-size particles which cannot overlap and repel each other when they are too close. The parts of the fluid where the maximal density is reached behave like incompressible fluids while lower density regions are compressible. This paper investigates the transition between the compressible and incompressible regions. To capture this transition, we study a one-dimensional Riemann problem and introduce a perturbation problem which regularizes the compressible-incompressible transition. Specific difficulties related to the non-conservativity of the problem are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Kinetic models for polymers with inertial effects

Novel kinetic models for both Dumbbell-like and rigid-rod like polymers are derived, based on the probability distribution function $f(t, x, n, \dot n)$ for a polymer molecule positioned at $x$ to be oriented along direction $n$ while embedded in a $\dot n$ environment created by inertial effects. It is shown that the probability distribution function of the extended model, when converging, will lead to well accepted kinetic models when inertial effects are ignored such as the Doi models for rod like polymers, and the Finitely Extensible Non-linear Elastic (FENE) models for Dumbbell like polymers.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Moment Guided Monte Carlo Method

In this work we propose a new approach for the numerical simulation of kinetic equations through Monte Carlo schemes. We introduce a new technique which permits to reduce the variance of particle methods through a matching with a set of suitable macroscopic moment equations. In order to guarantee that the moment equations provide the correct solutions, they are coupled to the kinetic equation through a non equilibrium term. The basic idea, on which the method relies, consists in guiding the particle positions and velocities through moment equations so that the concurrent solution of the moment and kinetic models furnishes the same macroscopic quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

An Asymptotic Preserving Scheme for the Euler equations in a strong magnetic field

This paper is concerned with the numerical approximation of the isothermal Euler equations for charged particles subject to the Lorentz force. When the magnetic field is large, the so-called drift-fluid approximation is obtained. In this limit, the parallel motion relative to the magnetic field direction splits from perpendicular motion and is given implicitly by the constraint of zero total force along the magnetic field lines. In this paper, we provide a well-posed elliptic equation for the parallel velocity which in turn allows us to construct an Asymptotic-Preserving (AP) scheme for the Euler-Lorentz system. This scheme gives rise to both a consistent approximation of the Euler-Lorentz model when epsilon is finite and a consistent approximation of the drift limit when epsilon tends to 0. Above all, it does not require any constraint on the space and time steps related to the small value of epsilon. Numerical results are presented, which confirm the AP character of the scheme and its Asymptotic Stability.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Continuum limit of self-driven particles with orientation interaction

We consider the discrete Couzin-Vicsek algorithm (CVA), which describes the interactions of individuals among animal societies such as fish schools. In this article, we propose a kinetic (mean-field) version of the CVA model and provide its formal macroscopic limit. The final macroscopic model involves a conservation equation for the density of the individuals and a non conservative equation for the director of the mean velocity and is proved to be hyperbolic. The derivation is based on the introduction of a non-conventional concept of a collisional invariant of a collision operator.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Large scale dynamics of the Persistent Turning Walker model of fish behavior

This paper considers a new model of individual displacement, based on fish motion, the so-called Persistent Turning Walker (PTW) model, which involves an Ornstein-Uhlenbeck process on the curvature of the particle trajectory. The goal is to show that its large time and space scale dynamics is of diffusive type, and to provide an analytic expression of the diffusion coefficient. Two methods are investigated. In the first one, we compute the large time asymptotics of the variance of the individual stochastic trajectories. The second method is based on a diffusion approximation of the kinetic formulation of these stochastic trajectories. The kinetic model is a Fokker-Planck type equation posed in an extended phase-space involving the curvature among the kinetic variables. We show that both methods lead to the same value of the diffusion constant. We present some numerical simulations to illustrate the theoretical results.

0 citations0 reviewsNo code linkedOpen access