In this paper, we systematically derive a fourth-order continuum theory capable of reproducing mesoscale turbulence in a three-dimensional suspension of microswimmers. We start from overdamped Langevin equations for a generic microscopic model (pushers or pullers), which include hydrodynamic interactions on both small length scales (polar alignment of neighboring swimmers) and large length scales, where the solvent flow interacts with the order parameter field. The flow field is determined via the Stokes equation supplemented by an ansatz for the stress tensor. In addition to hydrodynamic interactions, we allow for nematic pair interactions stemming from excluded-volume effects. The results here substantially extend and generalize earlier findings [S. Heidenreich et al., Phys. Rev. E 94, 020601 (2016)2470-004510.1103/PhysRevE.94.020601], in which we derived a two-dimensional hydrodynamic theory. From the corresponding mean-field Fokker-Planck equation combined with a self-consistent closure scheme, we derive nonlinear field equations for the polar and the nematic order parameter, involving gradient terms of up to fourth order. We find that the effective microswimmer dynamics depends on the coupling between solvent flow and orientational order. For very weak coupling corresponding to a high viscosity of the suspension, the dynamics of mesoscale turbulence can be described by a simplified model containing only an effective microswimmer velocity.

Recent experiments have shown that the complex spatio-temporal vortex structures emerging in active fluids are susceptible to weak geometrical constraints. This observation poses the fundamental question of how boundary effects stabilize a highly ordered pattern from seemingly turbulent motion. Here we show, by a combination of continuum theory and experiments on a bacterial suspension, how artificial obstacles guide the flow profile and reorganize topological defects, which enables the design of bacterial vortex lattices with tunable properties. To this end, the continuum model is extended by appropriate boundary conditions. Beyond the stabilization of square and hexagonal lattices, we also provide a striking example of a chiral, antiferromagnetic lattice exhibiting a net rotational flow, which is induced by arranging the obstacles in a Kagome-like array.

This paper studies the influence of orienting external fields on pattern formation, particularly mesoscale turbulence, in microswimmer suspensions. To this end, we apply a hydrodynamic theory that can be derived from a microscopic microswimmer model (Reinken et al 2018 Phys. Rev. E 97, 022613). The theory combines a dynamic equation for the polar order parameter with a modified Stokes equation for the solvent flow. Here, we extend the model by including an external field that exerts an aligning torque on the swimmers (mimicking the situation in chemo-, photo-, magneto-or gravitaxis). Compared to the field-free case, the external field breaks the rotational symmetry of the vortex dynamics and leads instead to strongly asymmetric, traveling stripe patterns, as demonstrated by numerical solution and linear stability analysis. We further analyze the emerging structures using a reduced model which involves only an (effective) microswimmer velocity field. This model is significantly easier to handle analytically, but still preserves the main features of the anisotropic pattern formation. We observe an underlying transition between a square vortex lattice and a traveling stripe pattern. These structures can be well described in the framework of weakly nonlinear analysis, provided the strength of nonlinear advection is sufficiently weak.

No abstract

We investigate the occurrence of shear banding in nematogenic fluids under planar Couette flow, based on mesoscopic dynamical equations for the orientational order parameter and the shear stress. We focus on parameter values where the sheared homogeneous system exhibits regular oscillatory orientational dynamics, whereas the equilibrium system is either isotropic (albeit close to the isotropic-nematic transition) or deep in its nematic phase. The numerical calculations are restricted to spatial variations in shear gradient direction. We find several new types of shear-banded states characterized by regions with regular oscillatory orientational dynamics. In all cases shear banding is accompanied by a non-monotonicity of the flow curve of the homogeneous system; however, only in the case of the initially isotropic system this curve has the typical S-like shape. We also analyze the influence of different orientational boundary conditions and of the spatial correlation length.