Three-dimensional turbulent flows enhance velocity gradients via strong nonlinear interactions of the rate-of-strain tensor with the vorticity vector, and with itself. For statistically homogeneous flows, their total contributions to gradient production are related to each other by conservation of mass, and so are the total enstrophy and total dissipation. However, locally, they do not obey this relation and have different (often extreme) values, and for this reason both production mechanisms have been subject to numerous studies, often decomposed into multi-scale interactions. In general lines, their dynamics and contributions to the cascade processes and turbulent kinetic dissipation are different, which poses a difficulty for turbulence modelling. In this paper, we explore the consequence of the ‘Betchov’ relations locally, and show that they implicitly define a length scale. This length scale is found to be approximately three times the size of the turbulent structures and their interactions. It is also found that, while the non-locality of the dissipation and enstrophy at a given scale comes mostly from larger scales that do not cancel, the non-local production of strain and vorticity comes from multi-scale interactions. An important consequence of this work is that isotropic cascade models need not distinguish between vortex stretching and strain self-amplification, but can instead consider both entities as part of a more complex transfer mechanism, provided that their detailed point value is not required and a local average of reasonable size is sufficient.

This study examines the effects of surface topography on the flow and performance of a self-propelled swimming (SPS) body. We consider a thin flat plate with an egg-carton roughness texture undergoing prescribed undulatory swimming kinematics at Strouhal number $0.3$ and tail amplitude to length ratio $0.1$; we use plate Reynolds numbers $Re=6$, 12 and $24\times 10^3$, and focus on $12\,000$. As the roughness wavelength is decreased, we find that the undulation wave speed must be increased to overcome the additional drag from the roughness and maintain SPS. Correspondingly, the extra wave speed raises the power required to maintain SPS, making the swimmer less efficient. To decouple the roughness and the kinematics, we compare the rough plates to equivalent smooth cases by matching the kinematic conditions. We find that all but the longest roughness wavelengths reduce the required swimming power and the unsteady amplitude of the forces when compared to a smooth plate undergoing identical kinematics. Additionally, roughness can enhance flow enstrophy by up to 116 % compared to the smooth cases without a corresponding spike in forces; this suggests that the increased mixing is not due to increased vorticity production at the wall. Instead, the enstrophy is found to peak strongly when the roughness wavelength is approximately twice the boundary layer thickness over the $Re$ range, indicating the roughness induces large-scale secondary flow structures that extend to the edge of the boundary layer. This study reveals the nonlinear interaction between roughness and kinematics beyond a simple increase or decrease in drag, illustrating that roughness studies on static shapes do not transfer directly to unsteady swimmers.

We investigate the nonlinear evolution of pairs of three-dimensional, equal-sized and opposite-signed vortices at finite Froude and Rossby numbers. The two vortices may be offset in the vertical direction. The initial conditions stem from relative equilibria obtained numerically in the quasi-geostrophic regime, for vanishing Froude and Rossby numbers. We first address the linear stability of the quasi-geostrophic opposite-signed pairs of vortices, and show that for all vertical offsets, the vortices are sensitive to an instability when close enough together. In the nonlinear regime, the instability may lead to the partial destruction of the vortices. We then address the nonlinear interaction of the vortices for various values of the Rossby number. We show that as the Rossby number increases, destructive interactions, where the vortices break into pieces, may occur for a larger separation between the vortices, compared to the quasi-geostrophic case. We also show that for well-separated vortices, the interaction is non-destructive, and ageostrophic effects lead to the deviation of the trajectory of the pair of vortices, as the anticyclonic vortex dominates the interaction. Finally, we show that the flow remains remarkably close to a balanced state, emitting only waves containing negligible energy, even when the interaction leads to the destruction of the vortices.

In the context of linear stability analysis, considering unsteady base flows is notoriously difficult. A generalisation of modal linear stability analysis, allowing for arbitrarily unsteady base flows over a finite time, is therefore required. The recently developed optimally time-dependent (OTD) modes form a projection basis for the tangent space. They capture the leading amplification directions in state space under the constraint that they form an orthonormal basis at all times. The present numerical study illustrates the possibility to describe a complex flow case using the leading OTD modes. The flow under investigation is an unsteady case of the Blasius boundary layer, featuring streamwise streaks of finite length and relevant to bypass transition. It corresponds to the state space trajectory initiated by the minimal seed; such a trajectory is unsteady, free from any spatial symmetry and shadows the laminar–turbulent separatrix for a finite time only. The finite-time instability of this unsteady base flow is investigated using the 8 leading OTD modes. The analysis includes the computation of finite-time Lyapunov exponents as well as instantaneous eigenvalues, and of the associated flow structures. The reconstructed instantaneous eigenmodes are all of outer type. They map unambiguously the spatial regions of largest instantaneous growth. Other flow structures, previously reported as secondary, are identified with this method as relevant to streak switching and to streamwise vortical ejections. The dynamics inside the tangent space features both modal and non-modal amplification. Non-normality within the reduced tangent subspace, quantified by the instantaneous numerical abscissa, emerges only as the unsteadiness of the base flow is reduced.

The near wake of a hemisphere immersed in a laminar boundary layer is studied utilizing time-resolved tomographic particle image velocimetry (TPIV). Focus is placed on the three-dimensional vortical structures and the formation details of hairpin vortices before the onset of transition. The three-dimensional instantaneous pressure field of the hemisphere wake is reconstructed for better understanding the flow mechanism. Experiments are carried out with Reynolds number $Re_{r}=1370$, based on the hemisphere radius $R$. Features of periodicity of the near wake are analysed using proper orthogonal decomposition and Fourier transformation. The velocity fluctuation in the wall-normal direction is shown to be crucial to the formation of hairpin vortices in the near wake. By investigating the transport of mass and vorticity, and the correlation between pressure and hairpin vortex strength, the formation mechanism is revealed clearly. Specifically, the main hairpin vortices (MHVs) are formed within the reaction of outer high-speed flow and near-wall flow. The formation of the head portion is followed by the leg portion formation. The shedding of the MHVs is highly correlated with the pressure, as well as the pressure gradient in the wall-normal direction. For the side hairpin vortices (SHVs), the leg portion is formed first, followed by the generation of the head portion thanks to induction of the re-oriented standing vortices. The generation of the SHVs can be regarded as the downstream bridging of the standing vortices, similar to the generation of hairpin vortices due to the connection of streamwise vortices in turbulent boundary layers.

Droplet coalescence is a common phenomenon and plays an important role in multidisciplinary applications. Previous studies mainly consider the coalescence of miscible liquids, even though the coalescence of immiscible droplets on a solid surface is a common process. In this study, we explore the coalescence of two immiscible droplets on a partial wetting surface experimentally and theoretically. We find that the coalescence process can be divided into three stages based on the time scales and force interactions involved, namely (I) the growth of a liquid bridge, (II) the oscillation of the coalescing sessile droplet and (III) the formation of a partially engulfed compound sessile droplet and the subsequent retraction. In stage I, the immiscible interface is found not to affect the scaling of the temporal evolution of the liquid bridge, which follows the same 2/3 power law as that of miscible droplets. In stage II, by developing a new capillary time scale considering both surface and interfacial tensions, we show that the interfacial tension between the two immiscible liquids functions as a non-negligible resistance to the oscillation which decreases the oscillation periods. In stage III, a modified Ohnesorge number is developed to characterize the visco-capillary and inertia-capillary time scales involved during the displacement of water by oil; a new model based on energy balance is proposed to analyse the maximum retraction velocity, highlighting that the viscous resistance is concentrated in a region close to the contact line.

The model eukaryotic microalgae Chlamydomonas reinhardtii is well known for its ability to generate bioconvection flows that are associated to intricate concentration patterns. Recently, it was demonstrated that the propensity of these algae to move toward a light source – a phenomenon termed phototaxis – can be exploited to locally concentrate micro-organisms and induce (photo)-bioconvection in algal suspensions by inducing a localised excess of density. In the present study we show experimentally that a cell population in a thin liquid layer self-organises in the presence of a heterogeneous light field and displays remarkable symmetry-breaking instabilities that are ruled by both the width of the light beam and the photo-bioconvection Rayleigh number. Beside circular stable states, fingers, dendrites and wave instabilities are reported, quantified and classified in a general phase diagram. Next, we use lubrication theory to develop an asymptotic model for bioconvection in a thin liquid layer, that includes the influences of both gyrotaxis and phototaxis. We obtain a single nonlinear anisotropic diffusion–drift equation describing the spatiotemporal evolution of the depth-averaged algal population. Analytical and numerical solutions are presented and show a very good agreement with the experimental results. In particular, we show that the dendrite instability arises as a result of a subtle coupling between the nonlinearity of the phototactic response, the gyrotactic effect and the self-induced bioconvective flows. Such complex flow fields might find applications in photo-bioreactors, through the efficient stirring of the harvested biomass.

We present numerical and analytical predictions of mucociliary clearance based on the continuum description of a viscoelastic mucus film, where momentum transfer from the beating cilia is represented via a Navier-slip boundary condition introduced by Bottier et al. (PLoS Comput. Biol., vol. 13, issue 7, 2017a, e1005552). Mucus viscoelasticity is represented via the Oldroyd-B model, where the relaxation time and the viscosity ratio have been fitted to experimental data for the storage and loss moduli of different types of real mucus, ranging from healthy to diseased conditions. We solve numerically the fully nonlinear governing equations for inertialess flow, and develop analytical solutions via asymptotic expansion in two limits: (i) weak viscoelasticity, i.e. low Deborah number; (ii) low cilia beat amplitude (CBA). All our approaches predict a drop in the mucus flow rate in relation to the Newtonian reference value, as the cilia beat frequency is increased. This relative drop increases with decreasing CBA and slip length. In diseased conditions, e.g. mucus properties characteristic of cystic fibrosis, the drop reaches 30 % in the physiological frequency range. In the case of healthy mucus, no significant drop is observed, even at very high frequency. This contrasts with the deterioration of microorganism propulsion predicted by the low-amplitude theory of Lauga (Phys. Fluids, vol. 19, issue 8, 2007, 083104), and is due to the larger beat amplitude and slip length associated with mucociliary clearance. In the physiological range of the cilia beat frequency, the low-amplitude prediction is accurate for both healthy and diseased conditions. Finally, we find that shear-thinning, modelled via a multi-mode Giesekus model, does not significantly alter our weakly viscoelastic and low-amplitude predictions based on the Oldroyd-B model.

The linear and nonlinear stability of two concentric jets separated by a duct wall is analysed by means of global linear stability and weakly nonlinear analysis. Three governing parameters are considered, the Reynolds number based on the inner jet, the inner-to-outer jet velocity ratio ($\delta _u$) and the length of the duct wall ($L$) separating the jet streams. Global linear stability analysis demonstrates the existence of unsteady modes of inherent convective nature, and symmetry-breaking modes that lead to a new non-axisymmetric steady state with a single or double helix. Additionally, we highlight the existence of multiple steady states, as a result of a series of saddle-node bifurcations and its connection to the changes in the topology of the flow. The neutral lines of stability have been computed for inner-to-outer velocity ratios within the range $0 < \delta _u < 2$ and duct wall distances in the interval $0.5 < L < 4$. They reveal the existence of hysteresis, and mode switching between two symmetry-breaking modes with azimuthal wavenumbers $1:2$. Finally, the mode interaction is analysed, highlighting the presence of travelling waves emerging from the resonant interaction of the two steady states, and the existence of robust heteroclinic cycles that are asymptotically stable.