Through boundary integral simulations and asymptotic analysis, we investigate the effect of a finite Navier slip length on the rheological proprieties of a dilute two-dimensional suspension of plate-like particles in the creeping flow limit. Specifically, we study the effects of Navier slip, particle thickness and Péclet number on the effective shear viscosity and average normal stress difference of an isolated two-dimensional plate-like particle in an unbounded shear flow field. We find that Navier slip induces a significant reduction in the effective viscosity and increases the average normal stress difference. The effect of slip becomes more enhanced as the thickness of the particle decreases and as the Péclet number increases. Remarkably, the analysis suggests that it is theoretically possible for a dilute suspension of slip plate-like particles at high Péclet numbers to have a shear viscosity smaller than that of the suspending fluid.

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.

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.

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.

Correct prediction of particle transport by surface waves is crucial in many practical applications such as search and rescue or salvage operations and pollution tracking and clean-up efforts. Recent results by Deike et al. (J. Fluid Mech., vol. 829, 2017, pp. 364–391) and Pizzo et al. (J. Phys. Oceanogr., vol. 49, no. 4, 2019, pp. 983–992) have indicated transport by deep-water breaking waves is enhanced compared with non-breaking waves. To model particle transport in irregular waves, some of which break, we develop a stochastic differential equation describing both mean particle transport and its uncertainty. The equation combines a Brownian motion, which captures non-breaking drift-diffusion effects, and a compound Poisson process, which captures jumps in particle positions due to breaking. From the corresponding Fokker–Planck equation for the evolution of the probability density function for particle position, we obtain closed-form expressions for its first three moments. We corroborate these predictions with new experiments, in which we track large numbers of particles in irregular breaking waves. For breaking and non-breaking wave fields, our experiments confirm that the variance of the particle position grows linearly with time, in accordance with Taylor's single-particle dispersion theory. For wave fields that include breaking, the compound Poisson process increases the linear growth rate of the mean and variance and introduces a finite skewness of the particle position distribution.

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

A thin liquid droplet spreads on a soft viscoelastic substrate with arbitrary rheology. Lubrication theory is applied to the governing field equations in the liquid and solid domains, which are coupled through the free boundary at the solid–liquid interface, to derive a set of reduced equations that describe the spreading dynamics. Fourier transform techniques and the finite difference method are used to construct a solution for the dynamic liquid–gas and solid–liquid interface shapes, as well as the macroscopic contact angle. Substrate properties affect the spreading dynamics through the contact angle and internal droplet flow fields, and these mechanisms are revealed. Increased substrate softness increases the spreading rate, whereas increased viscoelasticity decreases the spreading rate. For the case of a purely elastic substrate, the spreading power-law exponent recovers Tanner's law in the rigid limit and increases with substrate softness.

Eigenvectors of the observability and controllability Gramians represent responsive and receptive flow structures that enjoy a well-established connection to resolvent forcing and response modes. However, whereas resolvent modes have demonstrated great potential to guide sensor and actuator placement, observability and controllability modes have been leveraged exclusively in the context of model reduction via input and output projections. In this work, we introduce interpolatory, rather than orthogonal, input and output projections, that can be leveraged for sensor and actuator placement and open-loop control design. An interpolatory projector is an oblique projector with the property of preserving certain entries in the vector being projected. We review the connection between the resolvent operator and the Gramians, and present several numerical examples where we perform both orthogonal and interpolatory input and output projections onto the dominant forcing and response subspaces. Input projections are used to identify dynamically relevant disturbances, place sensors to measure disturbances, and place actuators for feedforward control in the linearized Ginzburg–Landau equation. Output projections are used to identify coherent structures and place sensors aiming at state reconstruction in the turbulent flow in a minimal channel at $Re_{\tau }=185$. The framework does not require data snapshots and relies only on knowledge of the steady or mean flow.

The whistling induced by a low-Mach turbulent flow through a deep axisymmetric cavity in a duct is investigated theoretically and experimentally. The experiments include acoustic measurements and stereoscopic particle image velocimetry (PIV). The paper focuses on the effect of a mean swirl on the dynamics of the azimuthal aeroacoustic modes. The mean swirl in the cavity has two origins: one component is imposed by a controlled tangential air injection upstream of the cavity, and the other component spontaneously arises under the action of the self-sustained azimuthal aeroacoustic mode, as explained in the companion paper, Part 1 (Faure-Beaulieu, Xiong, Pedergnana & Noiray, J. Fluid Mech., vol. 971, 2023, A21). Experiments show that the dynamics of the aeroacoustic wave is influenced by the imposed swirl. In particular, the spinning wave propagating against the swirl is promoted. To explain this, a linear perturbation analysis is performed around an incompressible mean swirling flow obtained from RANS simulations. It reveals that the dominant shear layer modes of azimuthal order 1 and −1 involved in the whistling phenomenon are helical modes winding respectively with and against the swirl, and spinning respectively in counterswirl and co-swirl directions. The counterswirl hydrodynamic mode is the least damped of the two, which is in agreement with the experimental observations. Finally, a low-order model based on the wave equation is derived. With only a few parameters, it fully reproduces the experimental observations for a wide range of imposed swirl intensity in the duct flow, and it allows us to disentangle the mechanisms responsible for this complex aeroacoustic instability.

In this paper, we study the aeroacoustic instability which occurs in a deep axisymmetric cavity in a turbulent pipe flow. This phenomenon is the axisymmetric counterpart of the classical whistling of a rectangular deep cavity subject to a grazing flow. The whistling of such axisymmetric cavity originates from the interaction of the coherent fluctuations of the vorticity at the cavity's opening with one of its trapped azimuthal or radial acoustic modes. We focus here on the situation involving the first pure azimuthal mode, which is trapped in the cavity. As a consequence of the rotational symmetry of the configuration, azimuthal modes are actually pairs of degenerate eigenmodes, or almost degenerate in the presence of small asymmetries. Therefore, the aeroacoustic instabilities exhibit more complex mechanisms than in the case of a rectangular deep cavity. In particular, we show that self-sustained spinning modes induce a symmetry breaking of the mean flow and we will elucidate the details of this phenomenon. To that end, simultaneous acoustic and time-resolved stereoscopic particle image velocimetry (PIV) measurements are performed. They reveal that when large-amplitude aeroacoustic waves spin around the cavity, a quasi-steady mean flow starts whirling slowly in the opposite direction to the wave propagation. A linear perturbation analysis around an axisymmetric mean flow confirms the experimental observations: although the incoming pipe flow is not swirling, the hydrodynamic component of the aeroacoustic wave induces such whirling motion of the mean flow because of the forcing from the steady part of the coherent Reynolds stress tensor.

To assess how the presence of surfactant in lung airways alters the flow of mucus that leads to plug formation and airway closure, we investigate the effect of insoluble surfactant on the instability of a viscoplastic liquid coating the interior of a cylindrical tube. Evolution equations for the layer thickness using thin-film and long-wave approximations are derived that incorporate yield-stress effects and capillary and Marangoni forces. Using numerical simulations and asymptotic analysis of the thin-film system, we quantify how the presence of surfactant slows growth of the Rayleigh–Plateau instability, increases the size of initial perturbation required to trigger instability and decreases the final peak height of the layer. When the surfactant strength is large, the thin-film dynamics coincide with the dynamics of a surfactant-free layer but with time slowed by a factor of four and the capillary Bingham number, a parameter proportional to the yield stress, exactly doubled. By solving the long-wave equations numerically, we quantify how increasing surfactant strength can increase the critical layer thickness for plug formation to occur and delay plugging. The previously established effect of the yield stress in suppressing plug formation (Shemilt et al., J. Fluid Mech., vol. 944, 2022, A22) is shown to be amplified by introducing surfactant. We discuss the implications of these results for understanding the impact of surfactant deficiency and increased mucus yield stress in obstructive lung diseases.

We report an experimental study of the motion of a clapping body consisting of two flat plates pivoted at the leading edge by a torsion spring. Clapping motion and forward propulsion of the body are initiated by the sudden release of the plates, initially held apart at an angle $2\theta _o$. Results are presented for the clapping and forward motions, and for the wake flow field for 24 cases, where depth-to-length ratio ($d^* = 1.5,1\text { and }0.5$), spring stiffness per unit depth ($Kt$), body mass ($m_b$) and initial separation angle ($2\theta _o = 45^{\circ }\text { and }60^{\circ }$) are varied. The body initially accelerates rapidly forward, then slowly retards to nearly zero velocity. Whereas the acceleration phase involves a complex interaction between plate and fluid motions, the retardation phase is simply fluid dynamic drag slowing the body. The wake consists of either a single axis-switching elliptical vortex loop (for $d^* = 1\text { and }1.5$) or multiple vortex loops (for $d^* = 0.5$). The body motion is nearly independent of $d^*$ and most affected by variation in $\theta _o$ and $Kt$. Using conservation of linear momentum and conversion of spring strain energy into kinetic energy in the fluid and body, we obtain a relation for the translation velocity of the body in terms of the various parameters. Approximately 80 % of the initial stored energy is transferred to the fluid, only 20 % to the body. The experimentally obtained cost of transport lies between 2 and $8\ \mathrm {J}\ \mathrm {kg}^{-1}\ \mathrm {m}^{-1}$.

We present a statistical characterization of the interaction between a planar shock and a finite-diameter, cylindrical column of dense gas based on three-dimensional, large-eddy simulation results. In the simulation, the column of gas is initially inclined at an angle $\alpha _0$ with respect to the shock plane. Effects of the initial column angle on the mixing characteristics are examined at Mach number 2.0 for column incline angles $1^\circ$, $5^\circ$, $10^\circ$ and $30^\circ$. Mean velocity profiles show that the column angle affects the gas velocity components in the vertical plane, but not in the spanwise direction. The gas undergoes higher initial upward acceleration at larger initial column incline angles. With time, the gas motion tends to become one-dimensional in the streamwise direction. Initially, velocity fluctuations are most intense within the interior of the column, but concentrate near the column leading edge over time. At high wavenumbers $\kappa$, the turbulent kinetic energy spectra follow a power-law scaling of $\kappa ^{-1}$. The structure functions of the mass fraction do not clearly demonstrate power-law scaling except at early times for $\alpha _0=30^\circ$, manifesting overall trends very similar to those observed in earlier experiments. Probability distributions of the mass fraction show independence of the mean and the standard deviation of the mixed gas on $\alpha _0$. The column angle was also found to have little effect on the mixing efficiency characterized by the molecular mixedness. Velocity components in the streamwise and transverse directions tend towards a bimodal distribution for larger $\alpha _{0}$.

The shock wave–turbulent boundary layer interaction over a compression corner is studied using global stability analysis (GSA) and resolvent analysis based on a separation of scales between the low-frequency, large-scale motions and the turbulent fluctuations. The GSA identifies a leading stationary mode, which becomes globally unstable as the ramp angle is beyond a critical value. For globally stable flows, the resolvent analysis captures two-dimensional and three-dimensional local maxima in optimal gain, both of which are due to modal resonance between the forcing and the leading global mode. Notably, the frequency-premultiplied optimal gain associated with two-dimensional disturbances peaks at a low frequency. For different interaction strengths, the peak frequencies collapse onto a universal value of 0.015 when non-dimensionalized using the length of the separation region and the free-stream velocity. A numerical simulation perturbed with the corresponding optimal forcing reveals that the response is in the form of a back-and-forth shock motion.

The classical model of evaporation of liquids hinges on Maxwell's assumption that the air near the liquid's surface is saturated. It allows one to find the evaporative flux without considering the interface separating liquid and air. Maxwell's hypothesis is based on an implicit assumption that the vapour-emission capacity of the interface exceeds the throughput of air (i.e. its ability to pass the vapour on to infinity). If this is indeed so, then the air adjacent to the liquid would get quickly saturated, justifying Maxwell's hypothesis. In the present paper, the so-called diffuse-interface model is used to account for the interfacial physics and thus derive a generalised version of Maxwell's boundary condition for the near-interface vapour density. It is then applied to a spherical drop floating in air. It turns out that the vapour-emission capacity of the interface exceeds the throughput of air only if the drop's radius is $r_{d}\gtrsim 10\ \mathrm {\mu } \mathrm {m}$, but for $r_{d}\approx 2\ \mathrm {\mu } {\rm m}$, the two are comparable. For $r_{d} \lesssim 1\ \mathrm {\mu } {\rm m}$, evaporation is interface-driven, and the resulting evaporation rate is noticeably smaller than that predicted by the classical model.

We carry out direct numerical simulations of flow in a plane open channel at friction Reynolds number up to ${{Re}}_{\tau } \approx 6000$. We find solid evidence for the presence of universal large-scale organization in the outer layer, with eddies that are larger and stronger than in the closed channel flow. As a result, velocity fluctuations are found to be stronger than in closed channels, throughout the depth. The inner-layer peak of the streamwise velocity variance is observed to grow logarithmically, as in Townsend's attached-eddy model (Townsend, The Structure of Turbulent Shear Flow, 2nd edn, Cambridge University Press, 1976), but saturation of the growth cannot be discarded based on the present data. Although we do not observe a clear outer peak of the streamwise velocity variance, we present substantial evidence that such a peak should emerge at a Reynolds number barely higher than achieved herein. The most striking feature of the flow is the presence of an extended logarithmic layer, with associated Kármán constant asymptoting to $k \approx 0.375$, in line with observations made in shear-free Couette–Poiseuille flow (Coleman et al., Flow Turbul. Combust., vol. 99, issue 3, 2017, pp. 553–564). The virtual absence of a wake region and of corrective terms to the log law in the present flow leads us to conclude that deviations from the log law observed in internal flows are likely due to the effects of the opposing walls, rather than the presence of a driving pressure gradient.

Submerged vegetation plays a subtle role in exchanging the fluid mass and energy in the vegetated flow zone, where the swaying motions of flexible plants are the important source of turbulent kinetic energy production. Flume experiments were conducted to study the modes, characteristics and factors of swaying of individual submerged flexible plants. A modified plant model in a new form, representing the highly flexible vegetation with clustered leaves, was employed. A ‘rigid-like’ synchronous swaying mode and a ‘whip-like’ asynchronous flapping mode are found to appear alternately for the individual plants. The interaction between these modes depends on the resulting local flow structure affected by the plants. Compared with a plant in isolation with the same flow Reynolds number, the swaying motions of a plant within the vegetation patch are less frequent but more prone to the synchronous mode. The eigen frequency of the motions increases linearly with an increase in flow Reynolds number in the range of 2 × 104–5 × 104, but the normalised amplitude reaches a saturation at a high flow Reynolds number. Moreover, the in-line and spanwise motions have a 2 : 1 frequency ratio for an ‘8’ shaped trajectory on the horizontal plane and a 1 : 1 ratio for a ‘0’ shaped circular trajectory, or a combination of both.

An experimental investigation of two-dimensional dispersively focused laboratory breaking waves is presented. We describe the bandwidth effect on breaking wave energetics, including spectral energy evolution, characteristic group velocity, energy dissipation and its rate, and breaking strength parameter, $b$. To evaluate the role of bandwidth, three definitions of wave group steepness are adopted where $S_s$ and $S_n$ are bandwidth-dependent and $S_p$ remains constant when bandwidth is changed. Our data show two regimes of spectral energy evolution in breaking wave groups, with both regimes bandwidth-dependent: energy dissipation and gain occur at $f > 0.95f_p$ ($f_p$ is the peak frequency) and $f < 0.95f_p$, respectively. The characteristic group velocity, which is used in energy dissipation calculations, increases by up to 7 % after wave breaking, being larger for higher bandwidth breaking waves. An unambiguous bandwidth dependence is found between $S_p$ and both the fractional and absolute wave energy dissipation. Wave groups of larger bandwidth break at a lower value of $S_p$ and consequently lose relatively more energy. The energy dissipation rate depends on the breaking duration which itself is bandwidth dependent. Consequently, no clear bandwidth effect is observed in energy dissipation rate when compared with either $S_p$ or $S_s$. However, there is a systematic bandwidth dependence in the variation of $b$ when parameterised in terms of $S_p$, with their relationship becoming increasingly nonlinear as bandwidth increases. When parameterised with $S_s$, $b$ shows a markedly reduced bandwidth dependence. Finally, the numerical breaking onset and relationship between $b$ and $S_s$ in the numerical study of Derakhti & Kirby (J. Fluid Mech., vol. 790, 2016, pp. 553–581) is validated experimentally.

In this work we investigate the effect of vertical confinement and inertia on the flow past thin ellipses in a Hele-Shaw cell (with centre line velocity $U_c$ and height 2$h$) with different aspect ratios for symmetrical flows and at an angle of attack, using asymptotic methods and numerical simulations. A Stokes region is identified at the ellipse vertices which results in flow different to flow past bluff bodies. Comparison with asymptotic analysis indicates close agreement over the ‘flat’ portion of the ellipse, for $\delta =(b/a)=0.05$, where $a$ and $b$ are the semi-major and -minor ellipse axes, respectively. Two flow conditions are investigated for ellipses at an angle of attack of 10$^\circ$ for a fixed $\delta =0.05$. Firstly, for $\varLambda =(U_ca/\nu )(h/a)^2 \ll 1$, the effect of increasing the vertical confinement of the Hele-Shaw cell results in the rear stagnation point (RSP) moving from close to the potential-flow prediction when $\epsilon =h/a$ is very small to the two-dimensional Stokes-flow prediction when $\epsilon$ is large. Secondly, for a fixed $\epsilon \ll 1$, when inertia is increased past $\varLambda ={O}(\epsilon )$ the RSP moves towards the trailing edge and is located there for $\varLambda ={O}(1)$. Under these conditions an attached exponentially decaying shear layer or ‘viscous tail’ is formed. A modified Bernoulli equation of the depth-averaged flow, together with the Kutta–Joukowski theorem is used to predict the drag and lift coefficients on the ellipse, which include a linear and a nonlinear contribution, corresponding to a Hele-Shaw and circulation component, respectively. Close agreement is found up to $\varLambda ={O}(1)$.