To recover the behavior of active particles as demonstrated by [1] in the hydrodynamic limit one might implement the so-called Toner-Tu equations in a Lattice Boltzmann Method [2]. This repository contains a D2Q9 Python implementation of the modified LBM which embeds the Contact Inhibition of Locomotion (CIL) model as proposed by [2].
[1] Tamás Vicsek et al. 1995, Novel Type of Phase Transition in a System of Self-Driven Particles
[2] David Nesbitt et al. 2021, Unconvering novel phase transitions in dense dry polar fluids using a lattice Boltzmann method