The concept of statistical stability is central to Malkus's 1956 attempt to predict the mean profile in shear flow turbulence. Here we discuss how his original attempt to assess this – an Orr–Sommerfeld (OS) analysis on the mean profile – can be improved by considering a cumulant expansion of the Navier–Stokes equations. Focusing on the simplest non-trivial closure (commonly referred to as CE2) that corresponds to the quasilinearized Navier–Stokes equations, we develop an extended OS analysis that also incorporates information about the fluctuation field. A more practical version of this – minimally extended OS analysis – is identified and tested on a number of statistically steady and, therefore, statistically stable turbulent channel flows. Beyond the concept of statistical stability, this extended stability analysis should also improve the popular approach of mean flow linear analysis in time-dependent shear flows by including more information about the underlying flow in its predictions as well as for other flows with additional physics such as convection.

Chiral fluids – such as fluids under rotation or a magnetic field as well as synthetic and biological active fluids – flow in a different way than ordinary ones. Due to symmetries broken at the microscopic level, chiral fluids may have asymmetric stress and viscosity tensors, for example giving rise to a hydrostatic torque or non-dissipative (odd) and parity-violating viscosities. In this article, we investigate the motion of rigid bodies in such an anisotropic fluid in the incompressible Stokes regime through the mobility matrix, which encodes the response of a solid body to forces and torques. We demonstrate how the form of the mobility matrix, which is usually determined by particle geometry, can be analogously controlled by the symmetries of the fluid. By computing the mobility matrix for simple shapes in a three-dimensional (3-D) anisotropic chiral fluid, we predict counterintuitive phenomena such as motion at an angle to the direction of applied forces and spinning under the force of gravity.

The present experimental study shows that a nozzle with optimal flexibility can enhance the impulse and entrainment of a pulsed jet. Near the nozzle exit, vortex rings emanating from the flexible nozzle move faster because of the timely release of the elastic energy (stored during the expansion) to the jet, which is maximized at the structural stiffness that needs to be optimally tuned to the jet acceleration. The total circulation, hydrodynamic impulse and entrained fluid volume are enhanced substantially. Interestingly, we find that the same condition for optimal flexibility to maximize the hydrodynamic impulse and circulation of the primary vortex ring of the continuous jet (Choi & Park, J. Fluid Mech., vol. 949, 2022, A39) holds universally for the pulsed jet, indicating that abrupt jet termination is irrelevant to the impulse augmentation mechanism. Compared to the rigid counterpart, increments of the impulse (${\sim }400\,\%$) and entrainment (${\sim }220\,\%$) of a pulsed jet in the present study are considerably larger than those ($200\,\%$ and $50\,\%$, respectively) in a continuous jet from previous studies, which is attributed to the significant suppression of negative pressure at the nozzle exit by the collapsing motion of the flexible nozzle in the phase with the jet-driven upstream propagation of the surface wave on the nozzle. This universal mechanism provides a guideline for a novel jet propulsor using a flexible nozzle, for example, for small-scale underwater robots.

Based on a generalized local Kolmogorov–Hill equation expressing the evolution of kinetic energy integrated over spheres of size $\ell$ in the inertial range of fluid turbulence, we examine a possible definition of entropy and entropy generation for turbulence. Its measurement from direct numerical simulations in isotropic turbulence leads to confirmation of the validity of the fluctuation relation (FR) from non-equilibrium thermodynamics in the inertial range of turbulent flows. Specifically, the ratio of probability densities of forward and inverse cascade at scale $\ell$ is shown to follow exponential behaviour with the entropy generation rate if the latter is defined by including an appropriately defined notion of ‘temperature of turbulence’ proportional to the kinetic energy at scale $\ell$.

In this paper, we study the spreading of droplets of density-matched granular suspensions on the surface of a solid. Bidispersity of the particle size distribution enriches the conclusions drawn from monodisperse experiments by highlighting key elements of the wetting dynamics. In all cases, the relation between the dynamic contact angle and the velocity of the contact line is similar to that for a simple fluid, despite the complexity introduced by the presence of particles. We extract from this relation an apparent wetting viscosity of the suspensions that differs from that measured in the bulk. Dimensional analysis supported by experimental measurements yields an estimate of the size of the region inside the droplet where the value of the dynamic contact angle depends on a balance of viscous dissipation and capillary stresses. How particle size compares with this viscous cut-off length seems crucial in determining the value of the apparent wetting viscosity. With bimodal blends, the particle size ratio can be used to show the effects of the local structure and volume fraction at the contact line, both impacting the value of the corresponding wetting viscosity.

We study the effect of insoluble surfactants on the impact of surfactant-free droplets onto surfactant-laden thin liquid films via a fully three-dimensional direct numerical simulation approach that employs a hybrid interface-tracking/level-set method, and by taking into account surfactant-induced Marangoni stresses due to gradients in interfacial surfactant concentration. Our numerical predictions for the temporal evolution of the surfactant-free crown are validated against the experimental work by Che & Matar (Langmuir, vol. 33, 2017, pp. 12140–12148). We focus on the ‘crown-splash regime’, and we observe that the crown dynamics evolves through various stages: from the growth of linear modes (through a Rayleigh–Plateau instability) to the development of nonlinearities leading to primary and secondary breakup events (through droplet shedding modulated by an end-pinching mechanism). We show that the addition of surfactants does not affect the wave selection via the Rayleigh–Plateau instability. However, the presence of surfactants plays a key role in the late stages of the dynamics as soon as the ligaments are driven out from the rim. Surfactant-induced Marangoni stresses delay the end-pinching mechanisms to result in longer ligaments prior to their capillary singularity. Our results indicate that Marangoni stresses bridge the gap between adjacent protrusions promoting the adjacent protrusions' collision and the merging of ligaments. Finally, we demonstrate that the addition of surfactants leads to surface rigidification and consequently to the retardation of the flow dynamics.

Density segregation of multi-component granular mixtures in a dense, gravity-driven flow over a rough and bumpy periodic chute surface is studied using theory and simulations. An existing theoretical model for predicting the steady-state concentration field of each species in a binary mixture using the forces acting on the particles is generalised for multi-component mixtures in this work. In addition, the rheological model for binary mixtures is also extended to multi-component mixtures. In contrast to the percolation velocity-based empirical segregation models that do not account for the rheology and need prior knowledge of the velocity field, the present approach accounts for the inter-coupling of rheology with segregation. The momentum balance equations are solved, along with the mixture rheological model as well as the multi-component density segregation model, to obtain concentration fields using an iterative numerical method. The theoretical predictions are compared with discrete element method (DEM) simulations for ternary, quaternary and quinary granular mixtures differing in density. The steady-state profiles for the concentration of different species as well as other flow properties predicted from the theory are in excellent agreement with the DEM simulation results for a variety of compositions over a range of inclination angles for different density ratios. Through this work, we take the first essential step towards proposing a generalised particle force-based segregation theory for multi-component mixtures differing in size and/or density.

We investigate how inhomogeneity influences the $k^{-5/3}$ inertial range scaling of turbulent kinetic energy spectra (with $k$ the wavenumber). For weak statistical inhomogeneity, the energy spectrum can be described as an equilibrium spectrum plus a perturbation. Theoretical arguments suggest that this latter contribution scales as $k^{-7/3}$. This prediction is assessed using direct numerical simulations of three-dimensional Kolmogorov flow.

Flexible canopy flows are often encountered in natural scenarios, e.g. when crops sway in the wind or when submerged kelp forests are agitated by marine currents. Here, we provide a detailed characterisation of the turbulent flow developed above and between the flexible filaments of a fully submerged dense canopy and we describe their dynamical response to the turbulent forcing. We investigate a wide range of flexibilities, encompassing the case in which the filaments are completely rigid and standing upright as well as that where they are fully compliant to the flow and deflected in the streamwise direction. We are thus able to isolate the effect of the canopy flexibility on the drag and on the inner–outer flow interactions, as well as the two flapping regimes of the filaments already identified for a single fibre. Furthermore, we offer a detailed description of the Reynolds stresses throughout the wall-normal direction resorting to the Lumley triangle formalism, and we show the multi-layer nature of turbulence inside and above the canopy. The relevance of our investigation is thus twofold: the fundamental physical understanding developed here paves the way towards the investigation of more complex and realistic scenarios, while we also provide a thorough characterisation of the turbulent state that can prove useful in the development of accurate turbulence models for RANS and LES.

Heat transfer by large deformable drops in a turbulent flow is a complex and rich-in-physics system, in which drop deformation, breakage and coalescence influence the transport of heat. We study this problem by coupling direct numerical simulation (DNS) of turbulence with a phase-field method for the interface description. Simulations are run at fixed-shear Reynolds and Weber numbers. To evaluate the influence of microscopic flow properties, like momentum/thermal diffusivity, on macroscopic flow properties, like mean temperature or heat transfer rates, we consider four different values of the Prandtl number, which is the momentum to thermal diffusivity ratio: $Pr=1$, $Pr=2$, $Pr=4$ and $Pr=8$. The drop volume fraction is $\varPhi \simeq 5.4\,\%$ for all cases. Drops are initially warmer than the turbulent carrier fluid and release heat at different rates depending on the value of $Pr$, but also on their size and on their own dynamics (topology, breakage, drop–drop interaction). Computing the time behaviour of the drops and carrier fluid average temperatures, we clearly show that an increase of $Pr$ slows down the heat transfer process. We explain our results by a simplified phenomenological model: we show that the time behaviour of the drop average temperature is self-similar, and a universal behaviour can be found upon rescaling by $t/Pr^{2/3}$. Accordingly, the heat transfer coefficient $\mathcal {H}$ (respectively its dimensionless counterpart, the Nusselt number $Nu$) scales as $\mathcal {H}\sim Pr^{-2/3}$ (respectively $Nu\sim Pr^{1/3}$) at the beginning of the simulation, and tends to $\mathcal {H}\sim Pr^{-1/2}$ (respectively $Nu\sim Pr^{1/2}$) at later times. These different scalings can be explained via the boundary layer theory and are consistent with previous theoretical/numerical predictions.

The statistics of breaking wave fields are characterised within a novel multi-layer framework, which generalises the single-layer Saint-Venant system into a multi-layer and non-hydrostatic formulation of the Navier–Stokes equations. We simulate an ensemble of phase-resolved surface wave fields in physical space, where strong nonlinearities, including directional wave breaking and the subsequent highly rotational flow motion, are modelled, without surface overturning. We extract the kinematics of wave breaking by identifying breaking fronts and their speed, for freely evolving wave fields initialised with typical wind wave spectra. The $\varLambda (c)$ distribution, defined as the length of breaking fronts (per unit area) moving with speed $c$ to $c+{\rm d}c$ following Phillips (J. Fluid Mech., vol. 156, 1985, pp. 505–531), is reported for a broad range of conditions. We recover the $\varLambda (c) \propto c^{-6}$ scaling without wind forcing for sufficiently steep wave fields. A scaling of $\varLambda (c)$ based solely on the root-mean-square slope and peak wave phase speed is shown to describe the modelled breaking distributions well. The modelled breaking distributions are in good agreement with field measurements and the proposed scaling can be applied successfully to the observational data sets. The present work paves the way for simulations of the turbulent upper ocean directly coupled to a realistic breaking wave dynamics, including Langmuir turbulence, and other sub-mesoscale processes.

Over the last decade, substantial progress has been made in understanding the topology of quasi-two-dimensional (2-D) non-equilibrium fluid flows driven by ATP-powered microtubules and microorganisms. By contrast, the topology of three-dimensional (3-D) active fluid flows still poses interesting open questions. Here, we study the topology of a spherically confined active flow using 3-D direct numerical simulations of generalized Navier–Stokes (GNS) equations at the scale of typical microfluidic experiments. Consistent with earlier results for unbounded periodic domains, our simulations confirm the formation of Beltrami-like bulk flows with spontaneously broken chiral symmetry in this model. Furthermore, by leveraging fast methods to compute linking numbers, we explicitly connect this chiral symmetry breaking to the entanglement statistics of vortex lines. We observe that the mean of linking number distribution converges to the global helicity, consistent with the asymptotic result by Arnold [In Vladimir I. Arnold – Collected Works (ed. A.B. Givental, B.A. Khesin, A.N. Varchenko, V.A. Vassiliev & O.Y. Viro), pp. 357–375. Springer]. Additionally, we characterize the rate of convergence of this measure with respect to the number and length of observed vortex lines, and examine higher moments of the distribution. We find that the full distribution is well described by a k-Gamma distribution, in agreement with an entropic argument. Beyond active suspensions, the tools for the topological characterization of 3-D vector fields developed here are applicable to any solenoidal field whose curl is tangent to or cancels at the boundaries of a simply connected domain.

Low Reynolds number turbulence in wall-bounded shear flows en route to laminar flow takes the form of oblique, spatially intermittent turbulent structures. In plane Couette flow, these emerge from uniform turbulence via a spatio-temporal intermittent process in which localised quasi-laminar gaps randomly nucleate and disappear. For slightly lower Reynolds numbers, spatially periodic and approximately stationary turbulent–laminar patterns predominate. The statistics of quasi-laminar regions, including the distributions of space and time scales and their Reynolds-number dependence, are analysed. A smooth, but marked transition is observed between uniform turbulence and flow with intermittent quasi-laminar gaps, whereas the transition from gaps to regular patterns is more gradual. Wavelength selection in these patterns is analysed via numerical simulations in oblique domains of various sizes. Via lifetime measurements in minimal domains, and a wavelet-based analysis of wavelength predominance in a large domain, we quantify the existence and nonlinear stability of a pattern as a function of wavelength and Reynolds number. We report that the preferred wavelength maximises the energy and dissipation of the large-scale flow along laminar–turbulent interfaces. This optimal behaviour is due primarily to the advective nature of the large-scale flow, with turbulent fluctuations playing only a secondary role.

We derive the evolution equation of the average uncertainty energy for periodic/homogeneous incompressible Navier–Stokes turbulence and show that uncertainty is increased by strain rate compression and decreased by strain rate stretching. We use three different direct numerical simulations (DNS) of non-decaying periodic turbulence and identify a similarity regime where (a) the production and dissipation rates of uncertainty grow together in time, (b) the parts of the uncertainty production rate accountable to average strain rate properties on the one hand and fluctuating strain rate properties on the other also grow together in time, (c) the average uncertainty energies along the three different strain rate principal axes remain constant as a ratio of the total average uncertainty energy, (d) the uncertainty energy spectrum's evolution is self-similar if normalised by the uncertainty's average uncertainty energy and characteristic length and (e) the uncertainty production rate is extremely intermittent and skewed towards extreme compression events even though the most likely uncertainty production rate is zero. Properties (a), (b) and (c) imply that the average uncertainty energy grows exponentially in this similarity time range. The Lyapunov exponent depends on both the Kolmogorov time scale and the smallest Eulerian time scale, indicating a dependence on random large-scale sweeping of dissipative eddies. In the two DNS cases of statistically stationary turbulence, this exponential growth is followed by an exponential of exponential growth, which is, in turn, followed by a linear growth in the one DNS case where the Navier–Stokes forcing also produces uncertainty.

The root systems of mangroves, a tree species found in intertidal tropical and subtropical coastal zones, provide a natural barrier that dissipates wave energy effectively and reduces sediment erosion. Here, we use a combination of experiments and numerical simulations to examine the wake and drag characteristics of porous arrays of cylinders, which serve as simplified models of mangrove root networks. Optimal arrangements of the arrays are obtained by coupling Navier–Stokes simulations with a multi-objective optimization algorithm, which seeks configurations that minimize wake enstrophy and maximize drag on the porous structure. These optimal configurations are investigated using particle image velocimetry, and the internal and external flows around the porous arrays are analysed using a combination of proper orthogonal decomposition and Lagrangian particle tracking. Large variations in drag and enstrophy are observed by varying the relative positions of the cylinders, which indicates that the geometrical arrangement of porous arrays plays a prominent role in determining wake and drag characteristics. A sensitivity analysis suggests that enstrophy is more sensitive than drag to specific cylinder placement, and depends on distinctive flow patterns that develop in the interior due to interactions among neighbouring cylinders. Arrays with higher drag involve a combination of larger projected frontal area and minimal flux through the interior, leading to increased wake enstrophy, which is unfavourable for particle deposition and erosion. Based on the analysis of characteristics associated with the optimal arrays, several manually designed arrays are tested, and they display the expected behaviour with regard to drag and enstrophy.

It is known that inertial lift forces can lead to particle focusing in channel flows; yet oscillatory straining effects have also been suggested as a mechanism for particle focusing in wavy channels. To explore the synergy between these two mechanisms, we analytically and experimentally investigate the focusing behaviour of rigid neutrally buoyant particles in a wavy channel. We decompose the particle-free channel flow into a primary Poiseuille flow and secondary eddies induced by the waviness. We calculate the perturbation of the particle on the particle-free flow and the resulting lateral lift force exerted on the particle using the method of matched asymptotic expansions. This yields a zeroth-order lift force arising from the Poiseuille flow and a first-order lift force due to the waviness of the channel. We further incorporate the inertial lift force into the Maxey–Riley equation and simulate particle trajectories in wavy channels. The interactions between the zeroth-order lift force and the particle-free flow largely determine the focusing locations. Experiments in wavy channels with varying amplitudes at channel Reynolds numbers ranging from 5 to 250 are consistent with the predictions of the focusing locations, which are mainly governed by the channel Reynolds number as well as the competition between the inertial lift and the oscillatory straining effects.

Advective dispersion of solutes in long thin axisymmetric channels is important to the analysis and design of a wide range of devices, including chemical separation systems and microfluidic chips. Despite extensive analysis of Taylor dispersion in various scenarios, most studies focus on long-term dispersion behaviour and cannot capture the transient evolution of the solute zone across the spatial variations in the channel. In the current study, we analyse the Taylor–Aris dispersion for arbitrarily shaped axisymmetric channels. We derive an expression for solute dynamics in terms of two coupled ordinary differential equations, which allow prediction of the time evolution of the mean location and axial (standard deviation) width of the solute zone as a function of the channel geometry. We compare and benchmark our predictions with Brownian dynamics simulations for a variety of cases, including linearly expanding/converging channels and periodic channels. We also present an analytical description of the physical regimes of transient positive versus negative axial growth of solute width. Finally, to further demonstrate the utility of the analysis, we demonstrate a method to engineer channel geometries to achieve desired solute width distributions over space and time. We apply the latter analysis to generate a geometry that results in a constant axial width and a second geometry that results in a sinusoidal axial variance in space.

In this paper a fully coupled particle-size segregation model for granular flows (Barker et al., J. Fluid Mech., vol. 909, 2021, p. A22) is used to simulate the development of the patterns in a triangular rotating drum. The results are compared with the experimental patterns formed with bidisperse and tridisperse granular mixtures, and with varying compositions and fill heights. In all cases the agreement between the simulations and experiments is remarkably good. The experimental patterns are generated in a narrow gap between transparent front and back sidewalls. These prevent three-dimensional motion, but also impose friction on the flow, making it thinner and faster than it would otherwise be. This promotes segregation, as it simultaneously increases the shear rate and reduces the local pressure. To obtain the correct flow dynamics and segregation, width-averaged sidewall friction is incorporated into the two-dimensional simulations, which are performed in OpenFOAM$^{\circledR}$. The free-surface avalanche forms a boundary layer within which all the segregation occurs. Material in the lower reach of the avalanche is continuously deposited into an underlying solid body of grains, which rotates with the drum, and is eventually re-entrained into the avalanche along its upper reach. The changing geometry of the granular region (as the drum rotates) implies that the avalanche is constantly adjusting its length, position and depth. This generates a complex quasi-periodic flow, which when combined with particle-size segregation generates amazing patterns in the solid rotating granular body after only two drum rotations.

In 2017, Brosseau & Vlahovska (Phys. Rev. Lett, vol. 119, no. 3, 2017, p. 034501) found that, in a strong electric field, a weakly conductive, low-viscosity droplet immersed in a highly conductive, high-viscosity medium formed a lens shape, and liquid rings continuously detached from its equatorial plane and subsequently broke up into satellite droplets. This fascinating multiphase electrohydrodynamic (EHD) phenomenon is known as droplet equatorial streaming. In this paper, based on the unified lattice Boltzmann method framework proposed by Luo et al. (Phil. Trans. R. Soc. A Math. Phys. Engng Sci, vol. 379, no. 2208, 2021, p. 20200397), a novel lattice Boltzmann (LB) model is constructed for multiphase EHD by coupling the Allen–Cahn type of multiphase LB model and two new LB equations to solve the Poisson equation of the electric field and the conservation equation of the surface charge. Using the proposed LB model, we successfully reproduced, for the first time, the complete process of droplet equatorial streaming, including the continuous ejection and breakup of liquid rings on the equatorial plane. In addition, it is found that, under conditions of high electric field strength or significant electrical conductivity contrast, droplets exhibit fingering equatorial streaming that was unknown before. A power-law relationship is discovered for droplet total charge evolution and a theoretical model is then proposed to describe the droplet radius and height over time. The breakup of liquid rings is found to be dominated by capillary instability, while the breakup of liquid fingers is governed by the end-pinching mechanism. Finally, a phase diagram is constructed for fingering equatorial streaming and ring equatorial streaming, and a criterion equation is established for the phase boundary.

We present a suite of large-eddy simulations (LES) of a wind farm operating in conventionally neutral boundary layers. A fixed 1.6 GW wind farm is considered for 40 different atmospheric stratification conditions to investigate effects on wind-farm efficiency and blockage, as well as related gravity-wave excitation. A tuned Rayleigh damping layer and a wave-free fringe-region method are used to avoid spurious excitation of gravity waves, and a domain-size study is included to evaluate and minimize effects of artificial domain blockage. A fully neutral reference case is also considered, to distinguish between a case with hydrodynamic blockage only, and cases that include hydrostatic blockage induced by the air column above the boundary layer and the excitation of gravity waves therein. We discuss in detail the dependence of gravity-wave excitation, flow fields and wind-farm blockage on capping-inversion height, strength and free-atmosphere lapse rate. In all cases, an unfavourable pressure gradient is present in front of the farm, and a favourable pressure gradient in the farm, with hydrostatic contributions arising from gravity waves at least an order of magnitude larger than hydrodynamic effects. Using respectively non-local and wake efficiencies $\eta _{nl}$ and $\eta _{w}$, we observe a strong negative correlation between the unfavourable upstream pressure rise and $\eta _{nl}$, and a strong positive correlation between the favourable pressure drop in the farm and $\eta _{w}$. Using a simplified linear gravity-wave model, we formulate a simple scaling for the ratio $(1-\eta _{nl})/\eta _{w}$, which matches reasonably well with the LES results.