Phase Transition to Large Scale Coherent Structures in Two-Dimensional Active Matter Turbulence

Turbulent cascade arrests and the formation of intermediate-scale condensates

Energy cascades lie at the heart of the dynamics of turbulent flows. In a recent study of turbulence in fluids with odd viscosity X. M. de Wit et al. [Nature (London) 627, 515 (2024)0028-083610.1038/s41586-024-07074-z], the two dimensionalization of the flow at small scales leads to the arrest of the energy cascade and selection of an intermediate scale, between the forcing and the viscous scales. To demonstrate the generality of the phenomenon and its existence for a wide class of turbulent systems, we study a shell model that is carefully constructed to have three-dimensional turbulent dynamics at small wave numbers and two-dimensional turbulent dynamics at large wave numbers. The large scale separation that we can achieve in our shell model allows us to examine clearly the interplay between these dynamics, which leads to an arrest of the energy cascade at a transitional wave number and an associated accumulation of energy at the same scale. Such pile-up of energy around the transitional wave number is reminiscent of the formation of condensates in two-dimensional turbulence, but, in contrast, it occurs at intermediate wave numbers instead of the smallest wave number.

Nonequilibrium hyperuniform states in active turbulence

We demonstrate that the complex spatiotemporal structure in active fluids can feature characteristics of hyperuniformity. Using a hydrodynamic model, we show that the transition from hyperuniformity to nonhyperuniformity and antihyperuniformity depends on the strength of active forcing and can be related to features of active turbulence without and with scaling characteristics of inertial turbulence. Combined with identified signatures of Levy walks and nonuniversal diffusion in these systems, this allows for a biological interpretation and the speculation of nonequilibrium hyperuniform states in active fluids as optimal states with respect to robustness and strategies of evasion and foraging.

Novel turbulence and coarsening arrest in active-scalar fluids

We uncover a new type of turbulence - activity-induced homogeneous and isotropic turbulence - in a model that has been employed to investigate motility-induced phase separation (MIPS) in a system of microswimmers. The active Cahn-Hilliard-Navier-Stokes (CHNS) equations, also called active model H, provide a natural theoretical framework for our study. In this CHNS model, a single scalar order parameter ϕ, positive (negative) in regions of high (low) microswimmer density, is coupled with the velocity field u. The activity of the microswimmers is governed by an activity parameter ζ that is positive for extensile swimmers and negative for contractile swimmers. With extensile swimmers, this system undergoes complete phase separation, which is similar to that in binary-fluid mixtures. By carrying out pseudospectral direct numerical simulations (DNSs), we show, for the first time, that (a) this model develops an emergent nonequilibrium, but statistically steady, state (NESS) of active turbulence, for the case of contractile swimmers, if ζ is sufficiently large and negative, and (b) this turbulence arrests the phase separation. We quantify this suppression by showing how the coarsening-arrest length scale does not grow indefinitely, with time t, but saturates at a finite value at large times. We characterise the statistical properties of this active-scalar turbulence by employing energy spectra and fluxes and the spectrum of ϕ. For sufficiently high Reynolds numbers, the energy spectrum (k) displays an inertial range, with a power-law dependence on the wavenumber k. We demonstrate that, in this range, the flux Π(k) assumes a nearly constant, negative value, which indicates that the system shows an inverse cascade of energy, even though energy injection occurs over a wide range of wavenumbers in our active-CHNS model.

Pattern formation by turbulent cascades

Fully developed turbulence is a universal and scale-invariant chaotic state characterized by an energy cascade from large to small scales at which the cascade is eventually arrested by dissipation1-6. Here we show how to harness these seemingly structureless turbulent cascades to generate patterns. Pattern formation entails a process of wavelength selection, which can usually be traced to the linear instability of a homogeneous state7. By contrast, the mechanism we propose here is fully nonlinear. It is triggered by the non-dissipative arrest of turbulent cascades: energy piles up at an intermediate scale, which is neither the system size nor the smallest scales at which energy is usually dissipated. Using a combination of theory and large-scale simulations, we show that the tunable wavelength of these cascade-induced patterns can be set by a non-dissipative transport coefficient called odd viscosity, ubiquitous in chiral fluids ranging from bioactive to quantum systems8-12. Odd viscosity, which acts as a scale-dependent Coriolis-like force, leads to a two-dimensionalization of the flow at small scales, in contrast with rotating fluids in which a two-dimensionalization occurs at large scales4. Apart from odd viscosity fluids, we discuss how cascade-induced patterns can arise in natural systems, including atmospheric flows13-19, stellar plasma such as the solar wind20-22, or the pulverization and coagulation of objects or droplets in which mass rather than energy cascades23-25.

Active turbulence and spontaneous phase separation in inhomogeneous extensile active gels

We report numerical results for the hydrodynamics of inhomogeneous lyotropic and extensile active nematic gels. By simulating the coupled Cahn-Hilliard, Navier-Stokes, and Beris-Edwards equation for the evolution of the composition, flow and orientational order of an active nematic, we ask whether composition variations are important to determine its emergent physics. As in active gels of uniform composition, we find that increasing either activity or nematic tendency (e.g., overall active matter concentration) triggers a transition between an isotropic passive phase and an active nematic one. We show that composition inhomogeneities are important in the latter phase, where we find three types of possible dynamical regimes. First, we observe regular patterns with defects and vortices: these exist close to the passive-active transition. Second, for larger activity, or deeper in the nematic phase, we find active turbulence, as in active gels of uniform composition, but with exceedingly large composition variation. In the third regime, which is uniquely associated with inhomogeneity and occurs for large nematic tendency and low activity, we observe spontaneous microphase separation into active and passive domains. The microphase separated regime is notable in view of the absence of an explicit demixing term in the underlying free energy which we use, and we provide a theoretical analysis based on the common tangent construction which explains its existence. We hope this regime can be probed experimentally in the future.

Coexistence of Active and Hydrodynamic Turbulence in Two-Dimensional Active Nematics

In active nematic liquid crystals, activity is able to drive chaotic spatiotemporal flows referred to as active turbulence. Active turbulence has been characterized through theoretical and experimental work as a low Reynolds number phenomenon. We show that, in two dimensions, the active forcing alone is able to trigger hydrodynamic turbulence leading to the coexistence of active and inertial turbulence. This type of flow develops for sufficiently active and extensile flow-aligning nematics. We observe that the combined effect of an extensile nematic and large values of the flow-aligning parameter leads to a broadening of the elastic energy spectrum that promotes a growth of kinetic energy able to trigger an inverse energy cascade.

Giant vortex dynamics in confined bacterial turbulence

We report the numerical evidence of a new state of bacterial turbulence in confined domains. By means of extensive numerical simulations of the Toner-Tu-Swift-Hohenberg model for dense bacterial suspensions in circular geometry, we discover the formation a stable, ordered state in which the angular momentum symmetry is broken. This is achieved by self-organization of a turbulent-like flow into a single, giant vortex of the size of the domain. The giant vortex is surrounded by an annular region close to the boundary, characterized by small-scale, radial vorticity streaks. The average radial velocity profile of the vortex is found to be in agreement with a simple analytical prediction. We also provide an estimate of the temporal and spatial scales of a suitable experimental setup comparable with our numerical findings.

Exact relations for energy transfer in simple and active binary fluid turbulence

Inertial range energy transfer in three-dimensional fully developed binary fluid turbulence is studied under the assumption of statistical homogeneity. Using two-point statistics, exact relations corresponding to the energy cascade are derived in terms of (i) two-point increments and (ii) two-point correlators. Despite having some apparent resemblances, the exact relation in binary fluid turbulence is found to be different from that of the incompressible magnetohydrodynamic turbulence [H. Politano and A. Pouquet, Geophys. Res. Lett. 25, 273 (1998)]0094-827610.1029/97GL03642. Besides the usual direct cascade of energy, under certain situations, an inverse cascade of energy is also speculated depending upon the strength of the activity parameter and the interplay between the two-point increments of the fluid velocity and the composition gradient fields. An alternative form of the exact relation is also derived in terms of the "upsilon" variables and a subsequent phenomenology is proposed predicting a k^{-3/2} law for the turbulent energy spectrum.

Nematic order condensation and topological defects in inertial active nematics

Living materials at different length scales manifest active nematic features such as orientational order, nematic topological defects, and active nematic turbulence. Using numerical simulations we investigate the impact of fluid inertia on the collective pattern formation in active nematics. We show that an incremental increase in inertial effects due to reduced viscosity results in gradual melting of nematic order with an increase in topological defect density before a discontinuous transition to a vortex-condensate state. The emergent vortex-condensate state at low enough viscosities coincides with nematic order condensation within the giant vortices and the drop in the density of topological defects. We further show flow field around topological defects is substantially affected by inertial effects. Moreover, we demonstrate the strong dependence of the kinetic energy spectrum on the inertial effects, recover the Kolmogorov scaling within the vortex-condensate phase, but find no evidence of universal scaling at higher viscosities. The findings reveal complexities in active nematic turbulence and emphasize the important cross-talk between active and inertial effects in setting flow and orientational organization of active particles.

Langevin based turbulence model and its relationship with Kappa distributions

Kappa distributions (or \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa $$\end{document}κ-like distributions) represent a robust framework to characterize and understand complex phenomena with high degrees of freedom, as turbulent systems, using non-extensive statistical mechanics. Here we consider a coupled map lattice Langevin based model to analyze the relation of a turbulent flow, with its spatial scale dynamic, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa $$\end{document}κ-like distributions. We generate the steady-state velocity distribution of the fluid at each scale, and show that the generated distributions are well fitted by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa $$\end{document}κ-like distributions. We observe a robust relation between the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa $$\end{document}κ parameter, the scale, and the Reynolds number of the system, Re. In particular, our results show that there is a closed scaling relation between the level of turbulence and the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa $$\end{document}κ parameter; namely \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa \sim \text {Re}\,k^{-5/3}$$\end{document}κ∼Rek-5/3. We expect these results to be useful to characterize turbulence in different contexts, and our numerical predictions to be tested by observations and experimental setups.

Fluctuation-Induced Interaction in Turbulent Flows

Fluctuation-induced forces are observed in numerous physical systems spanning from quantum to macroscopic scale. However, there is as yet no experimental report of their existence in hydrodynamic turbulence. Here, we present evidence of an attraction force mediated via turbulent fluctuations by using two walls locally confining 2D turbulence. This long-range interaction is a function of the wall separation and the energy injection rate in the turbulent flow. As the wall spacing decreases, the confined flow becomes less energetic and more anisotropic in the bounded domain, producing stronger attraction. The mechanism of force generation is rooted in a nontrivial fluid-wall coupling where coherent flow structures are guided by the cavity walls. For the narrowest cavities studied, a resonance phenomenon at the flow forcing scale leads to a complex short-range interaction. The results could be relevant to problems encountered in a range of fields from industrial multiphase flows to modeling of planetary formation.

Exact Coherent Structures and Phase Space Geometry of Preturbulent 2D Active Nematic Channel Flow

Confined active nematics exhibit rich dynamical behavior, including spontaneous flows, periodic defect dynamics, and chaotic "active turbulence." Here, we study these phenomena using the framework of exact coherent structures, which has been successful in characterizing the routes to high Reynolds number turbulence of passive fluids. Exact coherent structures are stationary, periodic, quasiperiodic, or traveling wave solutions of the hydrodynamic equations that, together with their invariant manifolds, serve as an organizing template of the dynamics. We compute the dominant exact coherent structures and connecting orbits in a preturbulent active nematic channel flow, which enables a fully nonlinear but highly reduced-order description in terms of a directed graph. Using this reduced representation, we compute instantaneous perturbations that switch the system between disparate spatiotemporal states occupying distant regions of the infinite-dimensional phase space. Our results lay the groundwork for a systematic means of understanding and controlling active nematic flows in the moderate- to high-activity regime.

Exploring order in active turbulence: Geometric rule and pairing order transition in confined bacterial vortices

Ordered collective motion emerges in a group of autonomously motile elements (known as active matter) as their density increases. Microswimmers, such as swimming bacteria, have been extensively studied in physics and biology. A dense suspension of bacteria forms seemingly chaotic turbulence in viscous fluids. Interestingly, this active turbulence driven by bacteria can form a hidden ensemble of many vortices. Understanding the active turbulence in a bacterial suspension can provide physical principles for pattern formation and insight into the instability underlying biological phenomena. This review presents recent findings regarding ordered structures causing active turbulence and discusses a physical approach for controlling active turbulence via geometric confinement. When the active matter is confined in a compartment with a size comparable to the correlation length of the collective motion, vortex-like rotation appears, and the vortex pairing order is indicated by the patterns of interacting vortices. Additionally, we outline the design principle for controlling collective motions via the geometric rule of the vortex pairing, which may advance engineering microdevices driven by a group of active matter. This article is an extended version of the Japanese article, Ordered Structure and Geometric Control of Active Matter in Dense Bacterial Suspensions, published in SEIBUTSU BUTSURI Vol. 60, p. 13–18 (2020).

Role of Advective Inertia in Active Nematic Turbulence

Suspensions of active agents with nematic interactions exhibit complex spatiotemporal dynamics such as mesoscale turbulence. Since the Reynolds number of microscopic flows is very small on the scale of individual agents, inertial effects are typically excluded in continuum theories of active nematic turbulence. Whether active stresses can collectively excite inertial flows is currently unclear. To address this question, we investigate a two-dimensional continuum theory for active nematic turbulence. In particular, we compare mesoscale turbulence with and without the effects of advective inertia. We find that inertial effects can influence the flow already close to the onset of the turbulent state and, moreover, give rise to large-scale fluid motion for strong active driving. A detailed analysis of the kinetic energy budget reveals an energy transfer to large scales mediated by inertial advection. While this transfer is small in comparison to energy injection and dissipation, its effects accumulate over time. The inclusion of friction, which is typically present in experiments, can compensate for this effect. The findings suggest that the inclusion of inertia and friction may be necessary for dynamically consistent theories of active nematic turbulence.

Imaging the emergence of bacterial turbulence: Phase diagram and transition kinetics

We experimentally study the emergence of collective bacterial swimming, a phenomenon often referred to as bacterial turbulence. A phase diagram of the flow of 3D Escherichia coli suspensions spanned by bacterial concentration, the swimming speed of bacteria, and the number fraction of active swimmers is systematically mapped, which shows quantitative agreement with kinetic theories and demonstrates the dominant role of hydrodynamic interactions in bacterial collective swimming. We trigger bacterial turbulence by suddenly increasing the swimming speed of light-powered bacteria and image the transition to the turbulence in real time. Our experiments identify two unusual kinetic pathways, i.e., the one-step transition with long incubation periods near the phase boundary and the two-step transition driven by long-wavelength instabilities deep inside the turbulent phase. Our study provides not only a quantitative verification of existing theories but also insights into interparticle interactions and transition kinetics of bacterial turbulence.

Understanding contagion dynamics through microscopic processes in active Brownian particles

Together with the universally recognized SIR model, several approaches have been employed to understand the contagion dynamics of interacting particles. Here, Active Brownian particles (ABP) are introduced to model the contagion dynamics of living agents that perform a horizontal transmission of an infectious disease in space and time. By performing an ensemble average description of the ABP simulations, we statistically describe susceptible, infected, and recovered groups in terms of particle densities, activity, contagious rates, and random recovery times. Our results show that ABP reproduces the time dependence observed in traditional compartmental models such as the Susceptible-Infected-Recovery (SIR) models and allows us to explore the critical densities and the contagious radius that facilitates the virus spread. Furthermore, we derive a first-principles analytical expression for the contagion rate in terms of microscopic parameters, without considering free parameters as the classical SIR-based models. This approach offers a novel alternative to incorporate microscopic processes into analyzing SIR-based models with applications in a wide range of biological systems.

Coarsening in the two-dimensional incompressible Toner-Tu equation: Signatures of turbulence

We investigate coarsening dynamics in the two-dimensional, incompressible Toner-Tu equation. We show that coarsening proceeds via vortex merger events, and the dynamics crucially depend on the Reynolds number Re. For low Re, the coarsening process has similarities to Ginzburg-Landau dynamics. On the other hand, for high Re, coarsening shows signatures of turbulence. In particular, we show the presence of an enstrophy cascade from the intervortex separation scale to the dissipation scale.

Soft channel formation and symmetry breaking in exotic active emulsions

We use computer simulations to study the morphology and rheological properties of a bidimensional emulsion resulting from a mixture of a passive isotropic fluid and an active contractile polar gel, in the presence of a surfactant that favours the emulsification of the two phases. By varying the intensity of the contractile activity and of an externally imposed shear flow, we find three possible morphologies. For low shear rates, a simple lamellar state is obtained. For intermediate activity and shear rate, an asymmetric state emerges, which is characterized by shear and concentration banding at the polar/isotropic interface. A further increment in the active forcing leads to the self-assembly of a soft channel where an isotropic fluid flows between two layers of active material. We characterize the stability of this state by performing a dynamical test varying the intensity of the active forcing and shear rate. Finally, we address the rheological properties of the system by measuring the effective shear viscosity, finding that this increases as active forcing is increased-so that the fluid thickens with activity.

Condensate formation and multiscale dynamics in two-dimensional active suspensions

The collective effects of microswimmers in active suspensions result in active turbulence, a spatiotemporally chaotic dynamics at mesoscale, which is characterized by the presence of vortices and jets at scales much larger than the characteristic size of the individual active constituents. To describe this dynamics, Navier-Stokes-based one-fluid models driven by small-scale forces have been proposed. Here, we provide a justification of such models for the case of dense suspensions in two dimensions (2D). We subsequently carry out an in-depth numerical study of the properties of one-fluid models as a function of the active driving in view of possible transition scenarios from active turbulence to large-scale pattern, referred to as condensate, formation induced by the classical inverse energy cascade in Newtonian 2D turbulence. Using a one-fluid model it was recently shown [M. Linkmann et al., Phys. Rev. Lett 122, 214503 (2019)10.1103/PhysRevLett.122.214503] that two-dimensional active suspensions support two nonequilibrium steady states, one with a condensate and one without, which are separated by a subcritical transition. Here, we report further details on this transition such as hysteresis and discuss a low-dimensional model that describes the main features of the transition through nonlocal-in-scale coupling between the small-scale driving and the condensate.

Lattice Boltzmann methods and active fluids

We review the state of the art of active fluids with particular attention to hydrodynamic continuous models and to the use of Lattice Boltzmann Methods (LBM) in this field. We present the thermodynamics of active fluids, in terms of liquid crystals modelling adapted to describe large-scale organization of active systems, as well as other effective phenomenological models. We discuss how LBM can be implemented to solve the hydrodynamics of active matter, starting from the case of a simple fluid, for which we explicitly recover the continuous equations by means of Chapman-Enskog expansion. Going beyond this simple case, we summarize how LBM can be used to treat complex and active fluids. We then review recent developments concerning some relevant topics in active matter that have been studied by means of LBM: spontaneous flow, self-propelled droplets, active emulsions, rheology, active turbulence, and active colloids.