Defect turbulence in a dense suspension of polar, active swimmers

We study the effects of inertia in dense suspensions of polar swimmers. The hydrodynamic velocity field and the polar order parameter field describe the dynamics of the suspension. We show that a dimensionless parameter $R$ (ratio of the swimmer self-advection speed to the active stress invasion speed) controls the stability of an ordered swimmer suspension. For $R$ smaller than a threshold $R_1$, perturbations grow at a rate proportional to their wave number $q$. Beyond $R_1$, we show that the growth rate is $\mathcal{O}(q^2)$ until a second threshold $R=R_2$ is reached. The suspension is stable for $R>R_2$. We perform direct numerical simulations to investigate the steady state properties and observe defect turbulence for $R<R_2$. An investigation of the spatial organisation of defects unravels a hidden transition: for small $R\approx 0$ defects are uniformly distributed and cluster as $R\to R_1$. Beyond $R_1$, clustering saturates and defects are arranged in nearly string-like structures.