Source author record

J. A. Carrillo

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

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

24 published item(s)

preprint2022arXiv

Nonlinear Aggregation-Diffusion Equations: Radial Symmetry and Long Time Asymptotics

We analyze under which conditions equilibration between two competing effects, repulsion modeled by nonlinear diffusion and attraction modeled by nonlocal interaction, occurs. This balance leads to continuous compactly supported radially decreasing equilibrium configurations for all masses. All stationary states with suitable regularity are shown to be radially symmetric by means of continuous Steiner symmetrization techniques. Calculus of variations tools allow us to show the existence of global minimizers among these equilibria. Finally, in the particular case of Newtonian interaction in two dimensions they lead to uniqueness of equilibria for any given mass up to translation and to the convergence of solutions of the associated nonlinear aggregation-diffusion equations towards this unique equilibrium profile up to translations as $t\to\infty$.

preprint2019arXiv

A Lipschitz metric for the Camassa--Holm equation

We analyze stability of conservative solutions of the Cauchy problem on the line for the Camassa--Holm (CH) equation. Generically, the solutions of the CH equation develop singularities with steep gradients while preserving continuity of the solution itself. In order to obtain uniqueness, one is required to augment the equation itself by a measure that represents the associated energy, and the breakdown of the solution is associated with a complicated interplay where the measure becomes singular. The main result in this paper is the construction of a Lipschitz metric that compares two solutions of the CH equation with the respective initial data. The Lipschitz metric is based on the use of the Wasserstein metric.

preprint2019arXiv

Long-time behaviour and phase transitions for the McKean--Vlasov equation on the torus

We study the McKean-Vlasov equation \[ \partial_t \varrho= β^{-1} Δ\varrho + κ\nabla \cdot (\varrho \nabla (W \star \varrho)) \, , \] with periodic boundary conditions on the torus. We first study the global asymptotic stability of the homogeneous steady state. We then focus our attention on the stationary system, and prove the existence of nontrivial solutions branching from the homogeneous steady state, through possibly infinitely many bifurcations, under appropriate assumptions on the interaction potential. We also provide sufficient conditions for the existence of continuous and discontinuous phase transitions. Finally, we showcase these results by applying them to several examples of interaction potentials such as the noisy Kuramoto model for synchronisation, the Keller--Segel model for bacterial chemotaxis, and the noisy Hegselmann--Krausse model for opinion dynamics.

preprint2016arXiv

Numerical Study of a Particle Method for Gradient Flows

We study the numerical behaviour of a particle method for gradient flows involving linear and nonlinear diffusion. This method relies on the discretisation of the energy via non-overlapping balls centred at the particles. The resulting scheme preserves the gradient flow structure at the particle level, and enables us to obtain a gradient descent formulation after time discretisation. We give several simulations to illustrate the validity of this method, as well as a detailed study of one-dimensional aggregation-diffusion equations.

preprint2015arXiv

Efficient numerical calculation of drift and diffusion coefficients in the diffusion approximation of kinetic equations

In this paper we study the diffusion approximation of a swarming model given by a system of interacting Langevin equations with nonlinear friction. The diffusion approximation requires the calculation of the drift and diffusion coefficients that are given as averages of solutions to appropriate Poisson equations. We present a new numerical method for computing these coefficients that is based on the calculation of the eigenvalues and eigenfunctions of a Schrödinger operator. These theoretical results are supported by numerical simulations showcasing the efficiency of the method.

preprint2014arXiv

Exponential Convergence Towards Stationary States for the 1D Porous Medium Equation with Fractional Pressure

We analyse the asymptotic behaviour of solutions to the one dimensional fractional version of the porous medium equation introduced by Caffarelli and Vázquez, where the pressure is obtained as a Riesz potential associated to the density. We take advantage of the displacement convexity of the Riesz potential in one dimension to show a functional inequality involving the entropy, entropy dissipation, and the Euclidean transport distance. An argument by approximation shows that this functional inequality is enough to deduce the exponential convergence of solutions in self-similar variables to the unique steady states.

preprint2014arXiv

Local well-posedness of the generalized Cucker-Smale model

In this paper, we study the local well-posedness of two types of generalized Cucker-Smale (in short C-S) flocking models. We consider two different communication weights, singular and regular ones, with nonlinear coupling velocities $v|v|^{β-2}$ for $β> \frac{3-d}{2}$. For the singular communication weight, we choose $ψ^1(x) = 1/|x|^α$ with $α\in (0,d-1)$ and $β\geq 2$ in dimension $d > 1$. For the regular case, we select $ψ^2(x) \geq 0$ belonging to $(L_{loc}^\infty \cap \mbox{Lip}_{loc})(\mathbb{R}^d)$ and $β\in (\frac{3-d}{2},2)$. We also remark the various dynamics of C-S particle system for these communication weights when $β\in (0,3)$.

preprint2014arXiv

Regularity of local minimizers of the interaction energy via obstacle problems

The repulsion strength at the origin for repulsive/attractive potentials determines the regularity of local minimizers of the interaction energy. In this paper, we show that if this repulsion is like Newtonian or more singular than Newtonian (but still locally integrable), then the local minimizers must be locally bounded densities (and even continuous for more singular than Newtonian repulsion). We prove this (and some other regularity results) by first showing that the potential function associated to a local minimizer solves an obstacle problem and then by using classical regularity results for such problems.

preprint2014arXiv

Renyi entropy and improved equilibration rates to self-similarity for nonlinear diffusion equations

We investigate the large-time asymptotics of nonlinear diffusion equations $u_t = Δu^p$ in dimension $n \ge 1$, in the exponent interval $p > n/(n+2)$, when the initial datum $u_0$ is of bounded second moment. Precise rates of convergence to the Barenblatt profile in terms of the relative Rényi entropy are demonstrated for finite-mass solutions defined in the whole space when they are re-normalized at each time $t> 0$ with respect to their own second moment. The analysis shows that the relative Rényi entropy exhibits a better decay, for intermediate times, with respect to the standard Ralston-Newton entropy. The result follows by a suitable use of the so-called concavity of Rényi entropy power.

preprint2014arXiv

Single to Double Mill Small Noise Transition via Semi-Lagrangian Finite Volume Methods

We show that double mills are more stable than single mills under stochastic perturbations in swarming dynamic models with basic attraction-repulsion mechanisms. In order to analyse accurately this fact, we will present a numerical technique for solving kinetic mean field equations for swarming dynamics. Numerical solutions of these equations for different sets of parameters will be presented and compared to microscopic and macroscopic results. As a consequence, we numerically observe a phase transition diagram in term of the stochastic noise going from single to double mill for small stochasticity fading gradually to disordered states when the noise strength gets larger. This bifurcation diagram at the inhomogeneous kinetic level is shown by carefully computing the distribution function in velocity space.

preprint2013arXiv

A new interaction potential for swarming models

We consider a self-propelled particle system which has been used to describe certain types of collective motion of animals, such as fish schools and bird flocks. Interactions between particles are specified by means of a pairwise potential, repulsive at short ranges and attractive at longer ranges. The exponentially decaying Morse potential is a typical choice, and is known to reproduce certain types of collective motion observed in nature, particularly aligned flocks and rotating mills. We introduce a class of interaction potentials, that we call Quasi-Morse, for which flock and rotating mills states are also observed numerically, however in that case the corresponding macroscopic equations allow for explicit solutions in terms of special functions, with coefficients that can be obtained numerically without solving the particle evolution. We compare thus obtained solutions with long-time dynamics of the particle systems and find a close agreement for several types of flock and mill solutions.

preprint2013arXiv

Contractivity of the Wasserstein metric for the kinetic Kuramoto equation

We present synchronization and contractivity estimates for the kinetic Kuramoto model obtained from the Kuramoto phase model in the mean-field limit. For identical Kuramoto oscillators, we present an admissible class of initial data leading to time-asymptotic complete synchronization, that is, all measure valued solutions converge to the traveling Dirac measure concentrated on the initial averaged phase. If two initial Radon measures have the same natural frequency density function and strength of coupling, we show that the Wasserstein p-distance between corresponding measure valued solutions is exponentially decreasing in time. This contraction principle is more general than previous L1-contraction properties of the Kuramoto phase model.

preprint2013arXiv

Nonlinear stability of flock solutions in second-order swarming models

In this paper we consider interacting particle systems which are frequently used to model collective behavior in animal swarms and other applications. We study the stability of orientationally aligned formations called flock solutions, one of the typical patterns emerging from such dynamics. We provide an analysis showing that the nonlinear stability of flocks in second-order models entirely depends on the linear stability of the first-order aggregation equation. Flocks are shown to be nonlinearly stable as a family of states under reasonable assumptions on the interaction potential. Furthermore, we numerically verify that commonly used potentials satisfy these hypotheses and investigate the nonlinear stability of flocks by an extensive case-study of uniform perturbations.

preprint2013arXiv

On the analysis of a coupled kinetic-fluid model with local alignment forces

This paper studies global existence, hydrodynamic limit, and large-time behavior of weak solutions to a kinetic flocking model coupled to the incompressible Navier-Stokes equations. The model describes the motion of particles immersed in a Navier-Stokes fluid interacting through local alignment. We first prove the existence of weak solutions using energy and $L^p$ estimates together with the velocity averaging lemma. We also rigorously establish a hydrodynamic limit corresponding to strong noise and local alignment. In this limit, the dynamics can be totally described by a coupled compressible Euler - incompressible Navier-Stokes system. The proof is via relative entropy techniques. Finally, we show a conditional result on the large-time behavior of classical solutions. Specifically, if the mass-density satisfies a uniform in time integrability estimate, then particles align with the fluid velocity exponentially fast without any further assumption on the viscosity of the fluid.

preprint2013arXiv

Splitting-Particle Methods for Structured Population Models: Convergence and Applications

We propose a new numerical scheme designed for a wide class of structured population models based on the idea of operator splitting and particle approximations. This scheme is related to the Escalator Boxcar Train (EBT) method commonly used in biology, which is in essence an analogue of particle methods used in physics. Our method exploits the split-up technique, thanks to which the transport step and the nonlocal integral terms in the equation can be separately considered. The order of convergence of the proposed method is obtained in the natural space of finite nonnegative Radon measures equipped with the flat metric. This convergence is studied even adding reconstruction and approximation steps in the particle simulation to keep the number of approximation particles under control. We validate our scheme in several test cases showing the theoretical convergence error. Finally, we use the scheme in situations in which the EBT method does not apply showing the flexibility of this new method to cope with the different terms in general structured population models.

preprint2013arXiv

Stability Analysis of Flock and Mill rings for 2nd Order Models in Swarming

We study the linear stability of flock and mill ring solutions of two individual based models for biological swarming. The individuals interact via a nonlocal interaction potential that is repulsive in the short range and attractive in the long range. We relate the instability of the flock rings with the instability of the ring solution of the first order model. We observe that repulsive-attractive interactions lead to new configurations for the flock rings such as clustering and fattening formation. Finally, we numerically explore mill patterns arising from this kind of interactions together with the asymptotic speed of the system.

preprint2013arXiv

The derivation of Swarming models: Mean-Field Limit and Wasserstein distances

These notes are devoted to a summary on the mean-field limit of large ensembles of interacting particles with applications in swarming models. We first make a summary of the kinetic models derived as continuum versions of second order models for swarming. We focus on the question of passing from the discrete to the continuum model in the Dobrushin framework. We show how to use related techniques from fluid mechanics equations applied to first order models for swarming, also called the aggregation equation. We give qualitative bounds on the approximation of initial data by particles to obtain the mean-field limit for radial singular (at the origin) potentials up to the Newtonian singularity. We also show the propagation of chaos for more restricted set of singular potentials.

preprint2012arXiv

Asymptotic Fixed-Speed Reduced Dynamics for Kinetic Equations in Swarming

We perform an asymptotic analysis of general particle systems arising in collective behavior in the limit of large self-propulsion and friction forces. These asymptotics impose a fixed speed in the limit, and thus a reduction of the dynamics to a sphere in the velocity variables. The limit models are obtained by averaging with respect to the fast dynamics. We can include all typical effects in the applications: short-range repulsion, long-range attraction, and alignment. For instance, we can rigorously show that the Cucker-Smale model is reduced to the Vicsek model without noise in this asymptotic limit. Finally, a formal expansion based on the reduced dynamics allows us to treat the case of diffusion. This technique follows closely the gyroaverage method used when studying the magnetic confinement of charged particles. The main new mathematical difficulty is to deal with measure solutions in this expansion procedure.

preprint2012arXiv

Confinement for Repulsive-Attractive Kernels

We investigate the confinement properties of solutions of the aggregation equation with repulsive-attractive potentials. We show that solutions remain compactly supported in a large fixed ball depending on the initial data and the potential. The arguments apply to the functional setting of probability measures with mildly singular repulsive-attractive potentials and to the functional setting of smooth solutions with a potential being the sum of the Newtonian repulsion at the origin and a smooth suitably growing at infinity attractive potential.

preprint2012arXiv

Uniqueness for Keller-Segel-type chemotaxis models

We prove uniqueness in the class of integrable and bounded nonnegative solutions in the energy sense to the Keller-Segel (KS) chemotaxis system. Our proof works for the fully parabolic KS model, it includes the classical parabolic-elliptic KS equation as a particular case, and it can be generalized to nonlinear diffusions in the particle density equation as long as the diffusion satisfies the classical McCann displacement convexity condition. The strategy uses Quasi-Lipschitz estimates for the chemoattractant equation and the above-the-tangent characterizations of displacement convexity. As a consequence, the displacement convexity of the free energy functional associated to the KS system is obtained from its evolution for bounded integrable initial data.

preprint2011arXiv

A mass-transportation approach to a one dimensional fluid mechanics model with nonlocal velocity

We consider a one dimensional transport model with nonlocal velocity given by the Hilbert transform and develop a global well-posedness theory of probability measure solutions. Both the viscous and non-viscous cases are analyzed. Both in original and in self-similar variables, we express the corresponding equations as gradient flows with respect to a free energy functional including a singular logarithmic interaction potential. Existence, uniqueness, self-similar asymptotic behavior and inviscid limit of solutions are obtained in the space $\mathcal{P}_{2}(\mathbb{R})$ of probability measures with finite second moments, without any smallness condition. Our results are based on the abstract gradient flow theory developed in \cite{Ambrosio}. An important byproduct of our results is that there is a unique, up to invariance and translations, global in time self-similar solution with initial data in $\mathcal{P}_{2}(\mathbb{R})$, which was already obtained in \textrm{\cite{Deslippe,Biler-Karch}} by different methods. Moreover, this self-similar solution attracts all the dynamics in self-similar variables. The crucial monotonicity property of the transport between measures in one dimension allows to show that the singular logarithmic potential energy is displacement convex. We also extend the results to gradient flow equations with negative power-law locally integrable interaction potentials.

preprint2011arXiv

Nonlocal interactions by repulsive-attractive potentials: radial ins/stability

In this paper, we investigate nonlocal interaction equations with repulsive-attractive radial potentials. Such equations describe the evolution of a continuum density of particles in which they repulse each other in the short range and attract each other in the long range. We prove that under some conditions on the potential, radially symmetric solutions converge exponentially fast in some transport distance toward a spherical shell stationary state. Otherwise we prove that it is not possible for a radially symmetric solution to converge weakly toward the spherical shell stationary state. We also investigate under which condition it is possible for a non-radially symmetric solution to converge toward a singular stationary state supported on a general hypersurface. Finally we provide a detailed analysis of the specific case of the repulsive-attractive power law potential as well as numerical results. We point out the the conditions of radial ins/stability are sharp.