Drop deformation in fluid flows is investigated as an exchange between the kinetic energy of the fluid and the surface energy of the drop. We show analytically that this energetic exchange is controlled only by the stretching (or compression) of the drop surface by the rate-of-strain tensor. This mechanism is analogous to the stretching of the vorticity field in turbulence. By leveraging the non-local nature of turbulence dynamics, we introduce a new decomposition that isolates the energetic exchange arising from local drop-induced surface effects from the non-local action of turbulent fluctuations. We perform direct numerical simulations of single inertial drops in isotropic turbulence and show that an important contribution to the increments of the surface energy arises from the non-local stretching of the fluid–fluid interface by eddies far from the drop surface (outer eddies). We report that this mechanism is dominant and independent of surface dynamics in a range of Weber numbers in which drop breakup occurs. These findings shed new light on drop deformation and breakup in turbulent flows, and opens the possibility for the improvement and simplification of breakup models.

The three-dimensional flow topology of a laminar separation bubble forming on the suction surface of a semispan wing with an aspect ratio of $2.5$ and NACA 0018 airfoil section is characterised experimentally using surface pressure measurements and particle image velocimetry at a chord Reynolds number of $125\ 000$. In the inboard region of the wing, the separation bubble is essentially two-dimensional, and the transition process in the separated shear layer leads to periodic vortex shedding, which dominates the bubble dynamics, similar to two-dimensional separation bubbles. However, progressive spanwise changes in the mean structure and vortex dynamics occur near the wingtip, leading to an open separation and eventual suppression of the bubble. In the immediate proximity of the wingtip, the boundary layer remains attached, no vortex shedding occurs and the flow remains laminar, terminating separation bubble formation. Despite variations in the mean separation bubble topology and vortex dynamics along the span, the fundamental shedding characteristics remain nearly invariant across the portion of the wing where vortex shedding occurs, and the flow appears to lock onto a common instability mode across the span, leading to minimal changes in the mean bubble characteristics despite notable changes in the effective angle of attack along the span. A comparison with available surface flow visualisations from previous studies indicates that the observed changes to the mean bubble footprint along the span of the wing are similar across different geometries and flow characteristics, suggesting similarities in the three-dimensional bubble topology and dynamics on finite wings.

In this paper, we study the fluid–structure interaction of a three-dimensional (3-D) flexible membrane immersed in an unsteady separated flow at moderate Reynolds numbers. We employ a body-conforming variational fluid–structure interaction solver based on the recently developed partitioned iterative scheme for the coupling of turbulent fluid flow with nonlinear structural dynamics. Of particular interest is to understand the flow-excited instability of a 3-D flexible membrane as a function of the non-dimensional mass ratio ($m^{*}$), Reynolds number ($Re$) and aeroelastic number ($Ae$). For a wide range of parameters, we examine two distinct stability regimes of the fluid–membrane interaction: deformed steady state (DSS) and dynamic balance state (DBS). We propose stability phase diagrams to demarcate the DSS and DBS regimes for the parameter space of mass ratio versus Reynolds number ($m^{*}$-$Re$) and mass ratio versus aeroelastic number ($m^{*}$-$Ae$). With the aid of the global Fourier mode decomposition technique, the distinct dominant vibrational modes are identified from the intertwined membrane responses in the parameter space of $m^{*}$-$Re$ and $m^{*}$-$Ae$. Compared to the deformed steady membrane, the flow-excited vibration produces relatively longer attached leading-edge vortices which improve the aerodynamic performance when the coupled system is near the flow-excited instability boundary. The optimal aerodynamic performance is achieved for lighter membranes with higher $Re$ and larger flexibility. Based on the global aeroelastic mode analysis, we observe a frequency lock-in phenomenon between the vortex-shedding frequency and the membrane vibration frequency causing self-sustained vibrations in the dynamic balance state. To characterize the origin of the frequency lock-in, we propose an approximate analytical formula for the nonlinear natural frequency by considering the added mass effect and employing a large deflection theory for a simply supported rectangular membrane. Through our systematic high-fidelity numerical investigation, we find that the onset of the membrane vibration and the mode transition has a direct dependence on the frequency lock-in between the natural frequency of the tensioned membrane and the vortex-shedding frequency or its harmonics. These findings on the fluid-elastic instability of membranes have implications for the design and development of control strategies for membrane wing-based unmanned systems and drones.

Vortex-induced vibration (VIV) of an elastically mounted sphere placed close to or piercing a free surface (FS) was investigated numerically. The submergence depth ($h$) was systematically varied between $1$ and $-$0.75 sphere diameters ($D$) and the response simulated over the reduced velocity range $U^*\in [3.5,14]$. The incompressible flow was coupled with the sphere motion modelled by a spring–mass–damper system, treating the free-surface boundary as a slip wall. In line with the previous experimental findings, as the submergence depth was decreased from $h^* = h/D =1$, the maximum response amplitude of the fully submerged sphere decreased; however, as the sphere pierced the FS, the amplitude increased until $h^* = -0.375$, and then decreased beyond that point. The fluctuating components of the lift and drag coefficients also followed the same pattern. The variation of the near-wake vortex dynamics over this submergence range was examined in detail to understand the effects of $h^*$ on the VIV response. It was found that $h^* = 1$ is a critical submergence depth, beyond which, as $h^*$ is decreased, the vortical structures in the wake vary significantly. For a fully submerged sphere, the influence of the stress-free condition on the VIV response was dominant over the kinematic constraint preventing flow through the surface. For piercing sphere cases, two previously unseen vortical recirculations were formed behind the sphere near times of maximal displacement, enhancing the VIV response. These were strongest at $h^* = -0.375$, and much weaker for small submergence depths, explaining the observed response-amplitude variation.

We study the stability of unsteady particle-laden flows in long, tilted water columns in batch settling mode, where the quasi-steady assumption of base flow no longer holds for the fast settling of particles. For this purpose, we introduce a settling time scale in the momentum and transport equations to solve the unsteady base flow, and utilise non-modal analysis to examine the stability of the disturbance flow field. The base flow increases in magnitude as the settling speed decreases and attains its maximum value when the settling speed becomes infinitesimal. The time evolution of the disturbance flow energy experiences an algebraic growth caused by the lift-up mechanism of the wall-normal disturbance, followed by an exponential growth owing to the shear instability of the base flow. The streamwise and spanwise wavenumbers corresponding to the peak energy gain are identified for both stages. In particular, the flow instability is enhanced as the Prandtl number increases, which is attributed to the sharpening of the particle-laden interface. On the other hand, the flow instability is suppressed by the increase in settling speed, because less disturbance energy can be extracted from the base flow. There exists an optimal tilted angle for efficient sedimentation, where the particle-laden flow is relatively stable and is accompanied by a smaller energy gain of the disturbance.

Although the behaviour of fluid-filled vesicles in steady flows has been extensively studied, far less is understood regarding the shape dynamics of vesicles in time-dependent oscillatory flows. Here, we investigate the nonlinear dynamics of vesicles in large amplitude oscillatory extensional (LAOE) flows using both experiments and boundary integral (BI) simulations. Our results characterize the transient membrane deformations, dynamical regimes and stress response of vesicles in LAOE in terms of reduced volume (vesicle asphericity), capillary number (${Ca}$, dimensionless flow strength) and Deborah number (${De}$, dimensionless flow frequency). Results from single vesicle experiments are found to be in good agreement with BI simulations across a wide range of parameters. Our results reveal three distinct dynamical regimes based on vesicle deformation: pulsating, reorienting and symmetrical regimes. We construct phase diagrams characterizing the transition of vesicle shapes between pulsating, reorienting and symmetrical regimes within the two-dimensional Pipkin space defined by ${De}$ and ${Ca}$. Contrary to observations on clean Newtonian droplets, vesicles do not reach a maximum length twice per strain rate cycle in the reorienting and pulsating regimes. The distinct dynamics observed in each regime result from a competition between the flow frequency, flow time scale and membrane deformation time scale. By calculating the particle stresslet, we quantify the nonlinear relationship between average vesicle stress and strain rate. Additionally, we present results on tubular vesicles that undergo shape transformation over several strain cycles. Broadly, our work provides new information regarding the transient dynamics of vesicles in time-dependent flows that directly informs bulk suspension rheology.

Air filaments and cavities in plunging breaking waves, generically cylinders, produce bubbles through an interface instability. The effects of gravity, surface tension and surface curvature on cylinder breakup are explored. A generalized dispersion relation is obtained that spans the Rayleigh–Taylor and Plateau–Rayleigh instabilities as cylinder radius varies. The analysis provides insight into the role of surface tension in the formation of bubbles from filaments and cavities. Small filaments break up into bubbles through a Plateau–Rayleigh instability driven through the action of surface tension. Large air cavities produce bubbles through a Rayleigh–Taylor instability driven by gravity and moderated by surface tension, which has a stabilizing effect. Surface tension, interface curvature and gravity are all important for cases between these two extremes. Predicted unstable mode wavenumber and bubble size show good agreement with direct numerical simulations of plunging breaking waves and air cylinders.

Electrophoresis of a tightly fitting sphere of radius $a$ along the centreline of a liquid-filled circular cylinder of radius $R$ is studied for a gap width $h_0=R-a\ll a$. We assume a Debye length $\kappa ^{-1}\ll h_0$, so that surface conductivity is negligible for zeta potentials typically seen in experiments, and the Smoluchowski slip velocity is imposed as a boundary condition at the solid surfaces. The pressure difference between the front and rear of the sphere is determined. If the cylinder has finite length $L$, this pressure difference causes an additional volumetric flow of liquid along the cylinder, increasing the electrophoretic velocity of the sphere, and an analytic prediction for this increase is found when $L\gg R$. If $N$ identical, well-spaced spheres are present, the electrophoretic velocity of the spheres increases with $N$, in agreement with the experiments of Misiunas & Keyser (Phys. Rev. Lett., vol. 122, 2019, 214501).