Modeling straight and circle swimmers: from single swimmer to collective motion.

We propose a simple numerical model for the motion of microswimmers based on the immersed boundary method. The swimmer, either pusher or puller, is represented by a distribution of point forces corresponding to the body and the flagellum. We study in particular the minimal model consisting of only three beads (two for the body and one for the flagellum) connected by rigid, inextensible links.

Immiscible Rayleigh-Taylor turbulence: Implications for bacterial degradation in oil spills.

An unstable density stratification between two fluids mixes spontaneously under the effect of gravity, a phenomenon known as Rayleigh-Taylor (RT) turbulence. If the two fluids are immiscible, for example, oil and water, surface tension prevents intermixing at the molecular level. However, turbulence fragments one fluid into the other, generating an emulsion in which the typical droplet size decreases over time as a result of the competition between the rising kinetic energy and the surface energy density.

Flocking turbulence of microswimmers in confined domains.

We extensively study the Toner-Tu-Swift-Hohenberg model of motile active matter by means of direct numerical simulations in a two-dimensional confined domain. By exploring the space of parameters of the model we investigate the emergence of a new state of active turbulence which occurs when the aligning interactions and the self-propulsion of the swimmers are strong. This regime of flocking turbulence is characterized by a population of few strong vortices, each surrounded by an island of coherent flocking motion.

Steering self-organisation through confinement.

Self-organisation is the spontaneous emergence of spatio-temporal structures and patterns from the interaction of smaller individual units. Examples are found across many scales in very different systems and scientific disciplines, from physics, materials science and robotics to biology, geophysics and astronomy. Recent research has highlighted how self-organisation can be both mediated and controlled by confinement.

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.

Incompressible Rayleigh-Taylor mixing in circular and spherical geometries.

We present a numerical investigation of the turbulent evolution of the mixing layer developing from the Rayleigh-Taylor instability for miscible incompressible fluids in circular (in two dimensions) and in spherical (in three dimensions) geometries in the Boussinesq approximation. We show that the main difference caused by the convergent geometry with respect to the planar case is that the center of the mixing layer drifts toward the center of the domain during the evolution of the mixing layer. A similar effect is observed for the radial profile of the density flux.

Dimensional effects in Rayleigh-Taylor mixing.

We study the effects of dimensional confinement on the evolution of incompressible Rayleigh-Taylor mixing both in a bulk flow and in porous media by means of numerical simulations of the transport equations. In both cases, the confinement to two-dimensional flow accelerates the mixing process and increases the speed of the mixing layer. Dimensional confinement also produces stronger correlations between the density and the velocity fields affecting the efficiency of the mass transfer, quantified by the dependence of the Nusselt number on the Rayleigh number.

Conformal invariance of weakly compressible two-dimensional turbulence.

We study conformal invariance of vorticity clusters in weakly compressible two-dimensional turbulence at low Mach numbers. On the basis of very high resolution direct numerical simulation we demonstrate the scaling invariance of the inverse cascade with scaling close to Kolmogorov prediction. In this range of scales, the statistics of zero-vorticity isolines are found to be compatible with those of critical percolation, thus generalizing the results obtained in incompressible Navier-Stokes turbulence.

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).

Alignment of Nonspherical Active Particles in Chaotic Flows.

We study the orientation statistics of spheroidal, axisymmetric microswimmers, with shapes ranging from disks to rods, swimming in chaotic, moderately turbulent flows. Numerical simulations show that rodlike active particles preferentially align with the flow velocity. To explain the underlying mechanism, we solve a statistical model via the perturbation theory.

Turbulence induces clustering and segregation of non-motile, buoyancy-regulating phytoplankton.

Turbulence plays a major role in shaping marine community structure as it affects organism dispersal and guides fundamental ecological interactions. Below oceanographic mesoscale dynamics, turbulence also impinges on subtle physical-biological coupling at the single cell level, setting a sea of chemical gradients and determining microbial interactions with profound effects on scales much larger than the organisms themselves. It has been only recently that we have started to disentangle details of this coupling for swimming microorganisms.

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

The collective motion of microswimmers in suspensions induce patterns of vortices on scales that are much larger than the characteristic size of a microswimmer, attaining a state called bacterial turbulence. Hydrodynamic turbulence acts on even larger scales and is dominated by inertial transport of energy. Using an established modification of the Navier-Stokes equation that accounts for the small-scale forcing of hydrodynamic flow by microswimmers, we study the properties of a dense suspension of microswimmers in two dimensions, where the conservation of enstrophy can drive an inverse cascade through which energy is accumulated on the largest scales.

Suppression of Rayleigh-Taylor turbulence by time-periodic acceleration.

The dynamics of Rayleigh-Taylor turbulence convection in the presence of an alternating, time-periodic acceleration is studied by means of extensive direct numerical simulations of the Boussinesq equations. Within this framework, we discover a mechanism of relaminarization of turbulence: the alternating acceleration, which initially produces a growing turbulent mixing layer, at longer times suppresses turbulent fluctuation and drives the system toward an asymptotic stationary configuration. Dimensional arguments and linear stability theory are used to predict the width of the mixing layer in the asymptotic state as a function of the period of the acceleration.

Gyrotactic phytoplankton in laminar and turbulent flows: A dynamical systems approach.

Gyrotactic algae are bottom heavy, motile cells whose swimming direction is determined by a balance between a buoyancy torque directing them upwards and fluid velocity gradients. Gyrotaxis has, in recent years, become a paradigmatic model for phytoplankton motility in flows. The essential attractiveness of this peculiar form of motility is the availability of a mechanistic description which, despite its simplicity, revealed predictive, rich in phenomenology, easily complemented to include the effects of shape, feedback on the fluid and stochasticity (e.

Buoyancy-Driven Flow through a Bed of Solid Particles Produces a New Form of Rayleigh-Taylor Turbulence.

Rayleigh-Taylor (RT) fluid turbulence through a bed of rigid, finite-size spheres is investigated by means of high-resolution direct numerical simulations, fully coupling the fluid and the solid phase via a state-of-the-art immersed boundary method. The porous character of the medium reveals a totally different physics for the mixing process when compared to the well-known phenomenology of classical RT mixing. For sufficiently small porosity, the growth rate of the mixing layer is linear in time (instead of quadratical) and the velocity fluctuations tend to saturate to a constant value (instead of linearly growing).

Time irreversibility in reversible shell models of turbulence.

Turbulent flows governed by the Navier-Stokes equations (NSE) generate an out-of-equilibrium time irreversible energy cascade from large to small scales. In the NSE, the energy transfer is due to the nonlinear terms that are formally symmetric under time reversal. As for the dissipative term: first, it explicitly breaks time reversibility; second, it produces a small-scale sink for the energy transfer that remains effective even in the limit of vanishing viscosity.

Point-particle method to compute diffusion-limited cellular uptake.

We present an efficient point-particle approach to simulate reaction-diffusion processes of spherical absorbing particles in the diffusion-limited regime, as simple models of cellular uptake. The exact solution for a single absorber is used to calibrate the method, linking the numerical parameters to the physical particle radius and uptake rate. We study the configurations of multiple absorbers of increasing complexity to examine the performance of the method by comparing our simulations with available exact analytical or numerical results.

Chaos and Predictability of Homogeneous-Isotropic Turbulence.

We study the chaoticity and the predictability of a turbulent flow on the basis of high-resolution direct numerical simulations at different Reynolds numbers. We find that the Lyapunov exponent of turbulence, which measures the exponential separation of two initially close solutions of the Navier-Stokes equations, grows with the Reynolds number of the flow, with an anomalous scaling exponent, larger than the one obtained on dimensional grounds. For large perturbations, the error is transferred to larger, slower scales, where it grows algebraically generating an "inverse cascade" of perturbations in the inertial range.

Scale-dependent colocalization in a population of gyrotactic swimmers.

We study the small scale clustering of gyrotactic swimmers transported by a turbulent flow, when the intrinsic variability of the swimming parameters within the population is considered. By means of extensive numerical simulations, we find that the variety of the population introduces a characteristic scale R^{*} in its spatial distribution. At scales smaller than R^{*} the swimmers are homogeneously distributed, while at larger scales an inhomogeneous distribution is observed with a fractal dimension close to what observed for a monodisperse population characterized by mean parameters.

Irreversibility of the two-dimensional enstrophy cascade.

We study the time irreversibility of the direct cascade in two-dimensional turbulence by looking at the time derivative of the square vorticity along Lagrangian trajectories, a quantity called metenstrophy. By means of extensive direct numerical simulations we measure the time irreversibility from the asymmetry of the probability density function of the metenstrophy and we find that it increases with the Reynolds number of the cascade, similarly to what is found in three-dimensional turbulence. A detailed analysis of the different contributions to the enstrophy budget reveals a remarkable difference with respect to what is observed for the energy cascade, in particular the role of the statistics of the forcing to determine the degree of irreversibility.

Centripetal focusing of gyrotactic phytoplankton.

A suspension of gyrotactic microalgae Chlamydomonas augustae swimming in a cylindrical water vessel in solid-body rotation is studied. Our experiments show that swimming algae form an aggregate around the axis of rotation, whose intensity increases with the rotation speed. We explain this phenomenon by the centripetal orientation of the swimming direction towards the axis of rotation.

Statistical conservation law in two- and three-dimensional turbulent flows.

Particles in turbulence live complicated lives. It is nonetheless sometimes possible to find order in this complexity. It was proposed in Falkovich et al.

Universal profile of the vortex condensate in two-dimensional turbulence.

An inverse turbulent cascade in a restricted two-dimensional periodic domain creates a condensate-a pair of coherent system-size vortices. We perform extensive numerical simulations of this system and carry out theoretical analysis based on momentum and energy exchanges between the turbulence and the vortices. We show that the vortices have a universal internal structure independent of the type of small-scale dissipation, small-scale forcing, and boundary conditions.

Turbulent channel without boundaries: the periodic Kolmogorov flow.

The Kolmogorov flow provides an ideal instance of a virtual channel flow: It has no boundaries, but it possesses well defined mean flow in each half wavelength. We exploit this remarkable feature for the purpose of investigating the interplay between the mean flow and the turbulent drag of the bulk flow. By means of a set of direct numerical simulations at increasing Reynolds number, we show the dependence of the bulk turbulent drag on the amplitude of the mean flow.

Dimensional transition in rotating turbulence.

In this work we investigate, by means of direct numerical hyperviscous simulations, how rotation affects the bidimensionalization of a turbulent flow. We study a thin layer of fluid, forced by a two-dimensional forcing, within the framework of the "split cascade" in which the injected energy flows both to small scales (generating the direct cascade) and to large scale (to form the inverse cascade). It is shown that rotation reinforces the inverse cascade at the expense of the direct one, thus promoting bidimensionalization of the flow.