skip to main content


Title: Global hypocoercivity of kinetic Fokker-Planck-Alignment equations

In this note we establish hypocoercivity and exponential relaxation to the Maxwellian for a class of kinetic Fokker-Planck-Alignment equations arising in the studies of collective behavior. Unlike previously known results in this direction that focus on convergence near Maxwellian, our result is global for hydrodynamically dense flocks, which has several consequences. In particular, if communication is long-range, the convergence is unconditional. If communication is local then all nearly aligned flocks quantified by smallness of the Fisher information relax to the Maxwellian. In the latter case the class of initial data is stable under the vanishing noise limit, i.e. it reduces to a non-trivial and natural class of traveling wave solutions to the noiseless Vlasov-Alignment equation.

The main novelty in our approach is the adaptation of a mollified Favre filtration of the macroscopic momentum into the communication protocol. Such filtration has been used previously in large eddy simulations of compressible turbulence and its new variant appeared in the proof of the Onsager conjecture for inhomogeneous Navier-Stokes system. A rigorous treatment of well-posedness for smooth solutions is provided. Lastly, we prove that in the limit of strong noise and local alignment solutions to the Fokker-Planck-Alignment equation Maxwellialize to solutions of the macroscopic hydrodynamic system with the isothermal pressure.

 
more » « less
Award ID(s):
2107956 1813351
NSF-PAR ID:
10329347
Author(s) / Creator(s):
Date Published:
Journal Name:
Kinetic and Related Models
Volume:
15
Issue:
2
ISSN:
1937-5093
Page Range / eLocation ID:
213
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We establish existence of finite energy weak solutions to the kinetic Fokker-Planck equation and the linear Landau equation near Maxwellian, in the presence of specular reflection boundary condition for general domains. Moreover, by using a method of reflection and the \begin{document}$ S_p $\end{document} estimate of [7], we prove regularity in the kinetic Sobolev spaces \begin{document}$ S_p $\end{document} and anisotropic Hölder spaces for such weak solutions. Such \begin{document}$ S_p $\end{document} regularity leads to the uniqueness of weak solutions.

     
    more » « less
  2. We propose an equilibrium-driven deformation algorithm (EDDA) to simulate the inbetweening transformations starting from an initial image to an equilibrium image, which covers images varying from a greyscale type to a colorful type on planes or manifolds. The algorithm is based on the Fokker-Planck dynamics on manifold, which automatically incorporates the manifold structure suggested by dataset and satisfies positivity, unconditional stability, mass conservation law and exponentially convergence. The thresholding scheme is adapted for the sharp interface dynamics and is used to achieve the finite time convergence. Using EDDA, three challenging examples, (I) facial aging process, (II) coronavirus disease 2019 (COVID-19) pneumonia invading/fading process, and (III) continental evolution process are computed efficiently.

     
    more » « less
  3. We consider the well-known Lieb-Liniger (LL) model for \begin{document}$ N $\end{document} bosons interacting pairwise on the line via the \begin{document}$ \delta $\end{document} potential in the mean-field scaling regime. Assuming suitable asymptotic factorization of the initial wave functions and convergence of the microscopic energy per particle, we show that the time-dependent reduced density matrices of the system converge in trace norm to the pure states given by the solution to the one-dimensional cubic nonlinear Schrödinger equation (NLS) with an explict rate of convergence. In contrast to previous work [3] relying on the formalism of second quantization and coherent states and without an explicit rate, our proof is based on the counting method of Pickl [65,66,67] and Knowles and Pickl [44]. To overcome difficulties stemming from the singularity of the \begin{document}$ \delta $\end{document} potential, we introduce a new short-range approximation argument that exploits the Hölder continuity of the \begin{document}$ N $\end{document}-body wave function in a single particle variable. By further exploiting the \begin{document}$ L^2 $\end{document}-subcritical well-posedness theory for the 1D cubic NLS, we can prove mean-field convergence when the limiting solution to the NLS has finite mass, but only for a very special class of \begin{document}$ N $\end{document}-body initial states.

     
    more » « less
  4. Efficient simulation of SDEs is essential in many applications, particularly for ergodic systems that demand efficient simulation of both short-time dynamics and large-time statistics. However, locally Lipschitz SDEs often require special treatments such as implicit schemes with small time-steps to accurately simulate the ergodic measures. We introduce a framework to construct inference-based schemes adaptive to large time-steps (ISALT) from data, achieving a reduction in time by several orders of magnitudes. The key is the statistical learning of an approximation to the infinite-dimensional discrete-time flow map. We explore the use of numerical schemes (such as the Euler-Maruyama, the hybrid RK4, and an implicit scheme) to derive informed basis functions, leading to a parameter inference problem. We introduce a scalable algorithm to estimate the parameters by least squares, and we prove the convergence of the estimators as data size increases.

    We test the ISALT on three non-globally Lipschitz SDEs: the 1D double-well potential, a 2D multiscale gradient system, and the 3D stochastic Lorenz equation with a degenerate noise. Numerical results show that ISALT can tolerate time-step magnitudes larger than plain numerical schemes. It reaches optimal accuracy in reproducing the invariant measure when the time-step is medium-large.

     
    more » « less
  5. In this note we study a new class of alignment models with self-propulsion and Rayleigh-type friction forces, which describes the collective behavior of agents with individual characteristic parameters. We describe the long time dynamics via a new method which allows us to reduce analysis from the multidimensional system to a simpler family of two-dimensional systems parametrized by a proper Grassmannian. With this method we demonstrate exponential alignment for a large (and sharp) class of initial velocity configurations confined to a sector of opening less than \begin{document}$ \pi $\end{document}.

    In the case when characteristic parameters remain frozen, the system governs dynamics of opinions for a set of players with constant convictions. Viewed as a dynamical non-cooperative game, the system is shown to possess a unique stable Nash equilibrium, which represents a settlement of opinions most agreeable to all agents. Such an agreement is furthermore shown to be a global attractor for any set of initial opinions.

     
    more » « less