Patchiness plays a fundamental role in phytoplankton ecology by dictating the rate at which individual cells encounter each other and their predators. The distribution of motile phytoplankton species is often considerably more patchy than that of non-motile species at submetre length scales, yet the mechanism generating this patchiness has remained unknown. Here we show that strong patchiness at small scales occurs when motile phytoplankton are exposed to turbulent flow. We demonstrate experimentally that Heterosigma akashiwo forms striking patches within individual vortices and prove with a mathematical model that this patchiness results from the coupling between motility and shear. When implemented within a direct numerical simulation of turbulence, the model reveals that cell motility can prevail over turbulent dispersion to create strong fractal patchiness, where local phytoplankton concentrations are increased more than 10-fold. This 'unmixing' mechanism likely enhances ecological interactions in the plankton and offers mechanistic insights into how turbulence intensity impacts ecosystem productivity.

This contribution provides a general framework to use Lagrange multipliers for the simulation of low Reynolds number fiber dynamics based on Bead Models (BM). This formalism provides an efficient method to account for kinematic constraints. We illustrate, with several examples, to which extent the proposed formulation offers a flexible and versatile framework for the quantitative modeling of flexible fibers deformation and rotation in shear flow, the dynamics of actuated filaments and the propulsion of active swimmers. Furthermore, a new contact model called Gears Model is proposed and successfully tested. It avoids the use of numerical artifices such as repulsive forces between adjacent beads, a source of numerical difficulties in the temporal integration of previous Bead Models.

The motility of microorganisms is often biased by gradients in physical and chemical properties of their environment, with myriad implications on their ecology. Here we show that fluid acceleration reorients gyrotactic plankton, triggering small-scale clustering. We experimentally demonstrate this phenomenon by studying the distribution of the phytoplankton Chlamydomonas augustae within a rotating tank and find it to be in good agreement with a new, generalized model of gyrotaxis. When this model is implemented in a direct numerical simulation of turbulent flow, we find that fluid acceleration generates multifractal plankton clustering, with faster and more stable cells producing stronger clustering. By producing accumulations in high-vorticity regions, this process is fundamentally different from clustering by gravitational acceleration, expanding the range of mechanisms by which turbulent flows can impact the spatial distribution of active suspensions.

We show that gyrotactic motility within a steady vortical flow leads to tightly clustered aggregations of microorganisms. Two dimensionless numbers, characterizing the relative swimming speed and stability against overturning by vorticity, govern the coupling between motility and flow. Exploration of parameter space reveals a striking array of patchiness regimes. Aggregations are found to form within a few overturning time scales, suggesting that vortical flows might be capable of efficiently separating species with different motility characteristics.

OATAO is an open access repository that collects the work of Toulouse researchers and makes it freely available over the web where possible.

This paper presents a numerical study of the dynamic self-assembly of neutrally buoyant particles rotating in a plane in a viscous fluid. The particles experience simultaneously a magnetic torque that drives their individual spinning motion, a magnetic attraction toward the center of the domain, and flow-induced interactions. A hydrodynamic repulsion balances the centripetal attraction of the magnetized particles and leads to the formation of an aggregate of several particles that rotates with a precession velocity related to the inter-particle distance. This dynamic self-assembly is stable (but not stationary) and the morphology depends on the number of particles. The repulsion force between the particles is shown to be the result of the secondary flow generated by each particle at low but nonzero Reynolds number. Comparisons are made with analogous experiments of spinning disks at a liquid–air interface, where it is found that the variation in the characteristic scales of the aggregate with the rotation rate of individual particles are consistent with the numerical results.

Based on direct numerical simulations of the coupled motion of particles and fluid, this study analyzes the collective hydrodynamic and colloidal effects of flowing microparticles during the formation of different 3D clogging patterns at a pore entrance. Simulations of flowing suspensions through a pore with various simulation conditions show that particle concentration and surface interactions play a major role in the occurrence of the bridging phenomenon (simultaneous adhesion of many particles). In the absence of DLVO repulsive forces, plugging is characterized by the temporal reduction of the bulk permeability when increasing the volume fraction of the flowing suspension up to 20%. Under these conditions, different structures of particle aggregates (from cluster to cake plug) are formed at the pore entrance yielding different evolution rates of hydrodynamic resistance to the flow. Taking into account DLVO repulsive forces in simulations for a particle concentration equal to 10%, we observe the transition from dendritic structures (for low repulsion) to dense aggregates (for high repulsion). At high DLVO repulsive forces, the scenario of pore clogging is controlled by the collective behavior of many interacting particles. It leads to the formation of a jamming phase (Wigner glass phase) with transient clusters of interacting particles at the pore entrance. The network of jammed particles collapses when the force chains among the particles are overcome by the flow stress. The build-up and the collapse of the jammed phase at the pore entrance induce temporal permeability fluctuations. According to the relative magnitude of particle-particle and particle-wall interactions, when the jammed phase is disorganized by the flow, the residual force in the network can accelerate particles and lead to particle adhesion at the wall inducing a pore blockage and a rapid reduction of the bulk permeability.

We investigate bubble dispersion in turbulent Taylor-Couette flow. The aim of this study is to describe the main mechanisms yielding preferential bubble accumulation in near-wall structures of the flow. We first proceed to direct numerical simulation of Taylor-Couette flows for three different geometrical configurations (three radius ratios η = R 1 /R 2 : η = 0.5, η = 0.72, and η = 0.91 with the outer cylinder at rest) and Reynolds numbers corresponding to turbulent regime ranging from 3000 to 8000. The statistics of the flow are discussed using two different averaging procedures that permit to characterize the mean azimuthal velocity, the Taylor vortices contribution and the small-scale turbulent fluctuations. The simulations are compared and validated with experimental and numerical data from literature. The second part of this study is devoted to bubble dispersion. Bubble accumulation is analyzed by comparing the dispersion obtained with the full turbulent flow field to bubble dispersion occurring at lower Reynolds numbers in previous works. Several patterns of preferential accumulation of bubbles have been observed depending on bubble size and the effect of gravity. For the smaller size considered, bubbles disperse homogeneously throughout the gap, while for the larger size they accumulate along the inner wall for the large gap width (η = 0.5). Varying the intensity of buoyancy yields complex evolution of the bubble spatial distribution. For low gravity effect, bubble entrapment is strong leading to accumulation along the inner wall in outflow regions (streaks of low wall shear stress). When buoyancy effect dominates on vortex trapping, bubbles rise through the vortices, while spiral patterns stretched along the inner cylinder are clearly identified. Force balance is analyzed to identify dominating forces leading to this accumulation and accumulation patterns are compared with previous experiments.