Model reduction of a kinetic swarming model by operator projection
This provides a new numerical method for simulating collective behavior in swarming models, but the contribution is incremental as it extends existing moment methods to a specific kinetic model.
The authors derive a globally hyperbolic moment system for the Vicsek swarming model using operator projection, proving its hyperbolicity, rotational invariance, and mass conservation. Numerical tests show convergence to the kinetic equation and capture of rarefaction, shock waves, and vortex formation.
This paper derives the arbitrary order globally hyperbolic moment system for a non-linear kinetic description of the Vicsek swarming model by using the operator projection. It is built on our careful study of a family of the complicate Grad type orthogonal functions depending on a parameter (angle of macroscopic velocity). We calculate their derivatives with respect to the independent vari- able, and projection of those derivatives, the product of velocity and basis, and collision term. The moment system is also proved to be hyperbolic, rotational invariant, and mass-conservative. The relationship between Grad type expansions in different parameter is also established. A semi-implicit numerical scheme is presented to solve a Cauchy problem of our hyperbolic moment system in order to verify the convergence behavior of the moment method. It is also compared to the spectral method for the kinetic equation. The results show that the solutions of our hyperbolic moment system converge to the solutions of the kinetic equation for the Vicsek model as the order of the moment system increases, and the moment method can capture key features such as rarefaction and shock waves, contact discontinuity, and vortex formation.