Source author record

Jose A. Carrillo

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

ResearcherUnclaimed source record

math.AP math.NA Numerical Analysis math-ph math.DS math.MP

Catalog footprint

What is connected

10works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Convergence of a Fully Discrete and Energy-Dissipating Finite-Volume Scheme for Aggregation-Diffusion Equations

We study an implicit finite-volume scheme for non-linear, non-local aggregation-diffusion equations which exhibit a gradient-flow structure, recently introduced by Bailo, Carrillo, and Hu (2020). Crucially, this scheme keeps the dissipation property of an associated fully discrete energy, and does so unconditionally with respect to the time step. Our main contribution in this work is to show the convergence of the method under suitable assumptions on the diffusion functions and potentials involved.

preprint2022arXiv

Noise-driven bifurcations in a neural field system modelling networks of grid cells

The activity generated by an ensemble of neurons is affected by various noise sources. It is a well-recognised challenge to understand the effects of noise on the stability of such networks. We demonstrate that the patterns of activity generated by networks of grid cells emerge from the instability of homogeneous activity for small levels of noise. This is carried out by analysing the robustness of network activity patterns with respect to noise in an upscaled noisy grid cell model in the form of a system of partial differential equations. Inhomogeneous network patterns are numerically understood as branches bifurcating from unstable homogeneous states for small noise levels. We show that there is a phase transition occurring as the level of noise decreases. Our numerical study also indicates the presence of hysteresis phenomena close to the precise critical noise value.

preprint2021arXiv

Primal dual methods for Wasserstein gradient flows

Combining the classical theory of optimal transport with modern operator splitting techniques, we develop a new numerical method for nonlinear, nonlocal partial differential equations, arising in models of porous media, materials science, and biological swarming. Our method proceeds as follows: First, we discretize in time, either via the classical JKO scheme or via a novel Crank-Nicolson type method we introduce. Next, we use the Benamou-Brenier dynamical characterization of the Wasserstein distance to reduce computing the solution of the discrete time equations to solving fully discrete minimization problems, with strictly convex objective functions and linear constraints. Third, we compute the minimizers by applying a recently introduced, provably convergent primal dual splitting scheme for three operators [Yan 2018]. By leveraging the PDEs' underlying variational structure, our method overcomes stability issues present in previous numerical work built on explicit time discretizations, which suffer due to the equations' strong nonlinearities and degeneracies. Our method is also naturally positivity and mass preserving and, in the case of the JKO scheme, energy decreasing. We prove that minimizers of the fully discrete problem converge to minimizers of the spatially continuous, discrete time problem as the spatial discretization is refined. We conclude with simulations of nonlinear PDEs and Wasserstein geodesics in one and two dimensions that illustrate the key properties of our approach, including higher order convergence our novel Crank-Nicolson type method, when compared to the classical JKO method.

preprint2021arXiv

Vortex formation for a non-local interaction model with Newtonian repulsion and superlinear mobility

We consider density solutions for gradient flow equations of the form $u_t = \nabla \cdot ( γ(u) \nabla \mathrm N(u))$, where $\mathrm N$ is the Newtonian repulsive potential in the whole space $\mathbb R^d$ with the nonlinear convex mobility $γ(u)=u^α$, and $α>1$. We show that solutions corresponding to compactly supported initial data remain compactly supported for all times leading to moving free boundaries as in the linear mobility case $γ(u)=u$. For linear mobility it was shown that there is a special solution in the form of a disk vortex of constant intensity in space $u=c_1t^{-1}$ supported in a ball that spreads in time like $c_2t^{1/d}$, thus showing a discontinuous leading front or shock. Our present results are in sharp contrast with the case of concave mobilities of the form $γ(u)=u^α$, with $0<α<1$ studied in [9]. There, we developed a well-posedness theory of viscosity solutions that are positive everywhere and moreover display a fat tail at infinity. Here, we also develop a well-posedness theory of viscosity solutions that in the radial case leads to a very detail analysis allowing us to show a waiting time phenomena. This is a typical behavior for nonlinear degenerate diffusion equations such as the porous medium equation. We will also construct explicit self-similar solutions exhibiting similar vortex-like behaviour characterizing the long time asymptotics of general radial solutions under certain assumptions. Convergent numerical schemes based on the viscosity solution theory are proposed analysing their rate of convergence. We complement our analytical results with numerical simulations ilustrating the proven results and showcasing some open problems.

preprint2020arXiv

A particle method for the homogeneous Landau equation

We propose a novel deterministic particle method to numerically approximate the Landau equation for plasmas. Based on a new variational formulation in terms of gradient flows of the Landau equation, we regularize the collision operator to make sense of the particle solutions. These particle solutions solve a large coupled ODE system that retains all the important properties of the Landau operator, namely the conservation of mass, momentum and energy, and the decay of entropy. We illustrate our new method by showing its performance in several test cases including the physically relevant case of the Coulomb interaction. The comparison to the exact solution and the spectral method is strikingly good maintaining 2nd order accuracy. Moreover, an efficient implementation of the method via the treecode is explored. This gives a proof of concept for the practical use of our method when coupled with the classical PIC method for the Vlasov equation.

preprint2020arXiv

Lagrangian schemes for Wasserstein gradient flows

This paper reviews different numerical methods for specific examples of Wasserstein gradient flows: we focus on nonlinear Fokker-Planck equations,but also discuss discretizations of the parabolic-elliptic Keller-Segel model and of the fourth order thin film equation. The methods under review are of Lagrangian nature, that is, the numerical approximations trace the characteristics of the underlying transport equation rather than solving the evolution equation for the mass density directly. The two main approaches are based on integrating the equation for the Lagrangian maps on the one hand, and on solution of coupled ODEs for individual mass particles on the other hand.

preprint2016arXiv

Sharp conditions to avoid collisions in singular Cucker-Smale interactions

We consider the Cucker-Smale flocking model with a singular communication weight $ψ(s) = s^{-α}$ with $α> 0$. We provide a critical value of the exponent $α$ in the communication weight leading to global regularity of solutions or finite-time collision between particles. For $α\geq 1$, we show that there is no collision between particles in finite time if they are placed in different positions initially. For $α\geq 2$ we investigate a version of the Cucker-Smale model with expanded singularity i.e. with weight $ψ_δ(s) = (s-δ)^{-α}$, $δ\geq 0$. For such model we provide a uniform with respect to the number of particles estimate that controls the $δ$-distance between particles. In case of $δ= 0$ it reduces to the estimate of non-collisioness.

preprint2011arXiv

An Asymptotic Preserving Scheme for the Diffusive Limit of Kinetic systems for Chemotaxis

In this work we numerically study the diffusive limit of run & tumble kinetic models for cell motion due to chemotaxis by means of asymptotic preserving schemes. It is well-known that the diffusive limit of these models leads to the classical Patlak-Keller-Segel macroscopic model for chemotaxis. We will show that the proposed scheme is able to accurately approximate the solutions before blow-up time for small parameter. Moreover, the numerical results indicate that the global solutions of the kinetic models stabilize for long times to steady states for all the analyzed parameter range. We also generalize these asymptotic preserving schemes to two dimensional kinetic models in the radial case. The blow-up of solutions is numerically investigated in all these cases.

preprint2010arXiv

On the Dynamics of a Fluid-Particle Interaction Model: The Bubbling Regime

This article deals with the issues of global-in-time existence and asymptotic analysis of a fluid-particle interaction model in the so-called bubbling regime. The mixture occupies the physical space $Ω\subset \mathbb{R}^3$ which may be unbounded. The system under investigation describes the evolution of particles dispersed in a viscous compressible fluid and is expressed by the conservation of fluid mass, the balance of momentum and the balance of particle density often referred as the Smoluchowski equation. The coupling between the dispersed and dense phases is obtained through the drag forces that the fluid and the particles exert mutually by action-reaction principle. We show that solutions exist globally in time under reasonable physical assumptions on the initial data, the physical domain, and the external potential. Furthermore, we prove the large-time stabilization of the system towards a unique stationary state fully determined by the masses of the initial density of particles and fluid and the external potential.

preprint2008arXiv

Strong Convergence towards homogeneous cooling states for dissipative Maxwell models

We show the propagation of regularity, uniformly in time, for the scaled solutions of the inelastic Maxwell model for small inelasticity. This result together with the weak convergence towards the homogenous cooling state present in the literature implies the strong convergence in Sobolev norms and in the $L^1$ norm towards it depending on the regularity of the initial data. The strategy of the proof is based on a precise control of the growth of the Fisher information for the inelastic Boltzmann equation. Moreover, as an application we obtain a bound in the $L^1$ distance between the homogeneous cooling state and the corresponding Maxwellian distribution vanishing as the inelasticity goes to zero.