We combine flow visualisation techniques and particle image velocimetry to experimentally investigate the higher-order transition to elastoinertial turbulence of Boger fluids ($El = 0.11\unicode{x2013}0.34$) in Taylor–Couette flows. The observed route to turbulence is associated with the appearance of chaotic inflow jets, termed flame patterns, for increasing inertia, and stable structures of solitons, known as diwhirls, for decreasing inertia. We also report for the first time spatially and temporally resolved flow fields in the meridional plane for the three characteristic viscoelastic flow regimes (diwhirls, flame patterns and elastoinertial turbulence). We observe in all cases coherent structures of dynamically independent solitary vortex pairs. The stability of these coherent structures is jet-dominated and can be mainly ascribed to the high extension of the polymer chains in the inflow boundaries in the $r$–$z$ plane. Solitary pairs are self-sustained when created through random events and do not split; instead, they merge when moving sufficiently close and annihilate when hoop stresses are not sufficient to sustain them. The highly localised and random events result in highly fluctuating, chaotic flow states. We estimate the decay exponent of spatial power spectral density, illustrating a universal scaling of $-2.5$ for elastoinertial turbulence. Based on our observations and in an effort to unify and combine precedent theories with our results, we suggest a mechanism for the origins of elastoinertial instabilities, accounting for both the effect of elasticity on the vortex formation and the effect of increasing/decreasing inertia on the flow dynamics.

Studying liquid jet impacts on a liquid pool is crucial for various engineering and environmental applications. During jet impact, the free surface of the pool deforms and a cavity is generated. Simultaneously, the free surface of the cavity extends radially outward and forms a rim. Eventually the cavity collapses by means of gas inertia and surface tension. Our numerical investigation using an axisymmetric model in Basilisk C explores cavity collapse dynamics under different impact velocities and gas densities. We validate our model against theory and experiments across a previously unexplored parameter range. Our results show two distinct regimes in the cavity collapse mechanism. By considering forces pulling along the interface, we derive scaling arguments for the time of closure and maximum radius of the cavity, based on the Weber number. For jets with uniform constant velocity from tip to tail and $We \leqslant 150$, the cavity closure is capillary-dominated and happens below the surface (deep seal). In contrast, for $We \geqslant 180$ the cavity closure happens above the surface (surface seal) and is dominated by the gas entrainment and the pressure gradient that it causes. Additionally, we monitor gas velocity and pressure throughout the impact process. This analysis reveals three critical moments of maximum gas velocity: before impact, at the instant of cavity collapse and during droplet ejection following cavity collapse. Our results provide information for understanding pollutant transport during droplet impacts on large bodies of water, and other engineering applications, like additive manufacturing, lithography and needle-free injections.

A lattice Boltzmann method is used to explore the effect of surfactants on the unequal volume breakup of a droplet in a T-junction microchannel, and the asymmetry due to fabrication defects in real-life microchannels is modelled as the pressure difference between the two branch outlets ($\Delta {P^\ast }$). We first study the effect of the surfactants on the droplet dynamics at different dimensionless initial droplet lengths ($l_0^\ast $) and capillary numbers (Ca) under symmetric boundary conditions ($\Delta {P^\ast } = 0$). The results indicate that the presence of surfactants promotes droplet deformation and breakup at small and moderate $l_0^\ast $ values, while the surfactant effect is weakened at large $l_0^\ast $ values. When the branch channels are completely blocked by the droplet, a linear relationship is observed between the dimensionless droplet length ($l_d^\ast $) and dimensionless time (${t^\ast }$), and two formulas are proposed for predicting the evolution of $l_d^\ast $ with ${t^\ast }$ for the two systems. We then investigate the effect of the surfactants on the droplet breakup at different values of $\Delta {P^\ast }$ and bulk surfactant concentrations (${\psi _b}$) under asymmetric boundary conditions ($\Delta {P^\ast } \ne 0$). It is observed that, as $\Delta {P^\ast }$ increases, the volume ratio of the generated droplets (${V_1}/{V_2}$) decreases to 0 in both systems, while the rate of decrease is higher in the clean system, i.e. the presence of surfactants could cause a decreased pressure difference between the droplet tips. As ${\psi _b}$ increases, ${V_1}/{V_2}$ first increases rapidly, then remains almost constant and finally decreases slightly. We thus establish a phase diagram that describes the ${V_1}/{V_2}$ variation with $\Delta {P^\ast }$ and ${\psi _b}$.

The transient processes of laminar separation bubble (LSB) formation and bursting on a rectangular NACA 0018 wing are studied experimentally. A two-dimensional airfoil model is used as a baseline for the assessment of finite wing effects. The models are subjected to ramp changes in free-stream velocity causing the flow to switch between a state where an LSB forms and a state without reattachment. Lift force and particle image velocimetry measurements are used to relate the flow development to the aerodynamic loading. The lift coefficient of the airfoil exhibits substantial hysteresis, and the duration of the lift transients range from $10$ to $22$ convective time scales for bubble formation and $22$ to $30$ convective time scales for bursting. In contrast, the transient lift coefficients of the wing change gradually, with less hysteresis. The wing tip causes greater three-dimensionality in the separation bubble, whose thickness increases near the midspan where bursting is initiated. During bubble formation, the region of separated flow contracts towards the midspan. The gradual change in lift of the wing is linked to slower spanwise expansion and contraction of the separated flow region relative to the airfoil. On both models, the wavenumbers and amplitudes of disturbances in the separated shear layer rapidly change when reattachment initiates or ceases. Applying the bursting criterion of Gaster (Tech. Rep. Reports and Memoranda 3595. Aeronautical Research Council, London, 1967) to the bubble on the wing shows that bursting of the bubble at a single spanwise location is insufficient to cause complete spanwise failure of reattachment, and that the relationship between bursting parameters depends on spanwise position.

Biological microorganisms and artificial micro-swimmers often locomote in heterogeneous viscous environments consisting of networks of obstacles embedded into viscous fluid media. In this work, we use the squirmer model and present a numerical investigation of the effects of shape on swimming in a heterogeneous medium. Specifically, we analyse the microorganism's propulsion speed as well as its energetic cost and swimming efficiency. The analysis allows us to probe the general characteristics of swimming in a heterogeneous viscous environment in comparison with the case of a purely viscous fluid. We found that a spheroidal microorganism always propels faster, expends less energy and is more efficient than a spherical microorganism in either a homogeneous fluid or a heterogeneous medium. Moreover, we determined that above a critical eccentricity, a spheroidal microorganism in a heterogeneous medium can swim faster than a spherical microorganism in a homogeneous fluid. Based on an analysis of the forces acting on the squirmer, we offer an explanation for the decrease in the squirmer's speed observed in heterogeneous media compared with homogeneous fluids.

Birds employ rapid pitch-up motions close to the ground for different purposes: perching birds use this motion to decelerate and come to a complete stop while hunting birds, such as bald eagles, employ it to catch prey and swiftly fly away. Motivated by these observations, our study investigates how natural flyers accomplish diverse flying objectives by rapidly pitching their wings while decelerating near ground. We conducted experimental and analytical investigations focusing on rapidly pitching plates during deceleration in close proximity to the ground to explore the impact of ground proximity on the unsteady dynamics. Initially, we executed synchronous pitch-up motion, where both pitching and deceleration have the same motion duration, at different ground heights. Experimental results demonstrate that as the pitching wing approaches the ground, the instantaneous lift increases by approximately $38\,\%$ compared with a far-from-ground case, while the initial peak drag force remains relatively unchanged. Our analytical model conforms to this trend, predicting an increase in lift force as the wing approaches the ground, indicating enhanced added mass and circulatory lift force due to the ground effect. Next, we examined asynchronous pitch-up motion cases, where rapid pitching motions were initiated at different stages of deceleration. The results reveal that initiating the wing pitch early in the deceleration leads to the formation of larger counter-rotating vortices at the early stage of the manoeuvre. These vortices generate stronger dipole jets that orient backward in the later stages of the manoeuvre after impinging with the ground surface, which hunting birds utilize to accelerate after catching prey. Conversely, when the wing pitch is delayed, smaller vortices form, but their growth is postponed until late in the manoeuvre. This delayed vortex growth produces lift and drag force at the end phase of the manoeuvre that facilitates a smooth landing or perching. Thus, through strategic tuning of a rapid pitch-up motion with deceleration, natural flyers, such as birds, achieve diverse flying objectives.

The wall-attached structure characteristics of flow and dust concentration fields are investigated in this study based on high-Reynolds-number ($Re_{\tau } \sim O(10^{6})$) synchronous multiphase observations from the Qingtu Lake observation array site (recorded by Liu et al., J. Fluid Mech., vol. 957, 2023, A14). The results show that not only the particle-free flow field but also the particle-laden flow and dust concentration field contain wall-attached structures. The linear coherence spectrum, as a data-driven filter, is adopted to separate the wall-attached portions from the premultiplied spectra. The decomposed spectra for wall-attached structures in streamwise velocity fluctuations under particle-free and particle-laden conditions show obvious outer-scaling and wall-scaling at large-scale and medium-scale ranges, respectively, but the spectra of dust concentration exhibit only wall-scaling rather than outer-scaling. The streamwise length of the most significant wall-attached dust clustering structures in the logarithmic region is approximately five times the boundary layer thickness, and does not change significantly with height. Furthermore, the streamwise turbulence intensity for wall-attached portions follows the universal Townsend–Perry constant, without particle mass loading effect. Correspondingly, the wall-attached portions of the dust concentration also exhibit universal logarithmic decay slope. The remaining non-attached portions for the turbulent velocity and dust concentration have significant dependence on the particle mass loading.

We study the instantaneous inference of an unbounded planar flow from sparse noisy pressure measurements. The true flow field comprises one or more regularized point vortices of various strength and size. We interpret the true flow's measurements with a vortex estimator, also consisting of regularized vortices, and attempt to infer the positions and strengths of this estimator assuming little prior knowledge. The problem often has several possible solutions, many due to a variety of symmetries. To deal with this ill posedness and to quantify the uncertainty, we develop the vortex estimator in a Bayesian setting. We use Markov-chain Monte Carlo and a Gaussian mixture model to sample and categorize the probable vortex states in the posterior distribution, tailoring the prior to avoid spurious solutions. Through experiments with one or more true vortices, we reveal many aspects of the vortex inference problem. With fewer sensors than states, the estimator infers a manifold of equally possible states. Using one more sensor than states ensures that no cases of rank deficiency arise. Uncertainty grows rapidly with distance when a vortex lies outside of the vicinity of the sensors. Vortex size cannot be reliably inferred, but the position and strength of a larger vortex can be estimated with a much smaller one. In estimates of multiple vortices their individual signs are discernible because of the nonlinear coupling in the pressure. When the true vortex state is inferred from an estimator of fewer vortices, the estimate approximately aggregates the true vortices where possible.

Direct numerical simulations are performed to investigate the characteristics of a turbulent core–annular flow with water-lubricated high viscosity oil in a horizontal pipe. Six different superficial velocity ratios ($\,j_w/j_o = 0.057\unicode{x2013}0.41$) are examined by changing the water superficial velocity $j_w$ for a fixed oil superficial velocity $j_o$. The pressure drops in the pipe and the shapes of the phase interface agree well with those from previous experiments. The oil core flow is almost a plug flow, and the gaps between the phase interface and pipe wall are narrow and wide near the upper and lower surfaces of the pipe, respectively, due to the buoyancy. Within a narrow gap, water is confined mostly in a valley region of the wavy phase interface and hardly goes through its crest. On the other hand, water near the phase interface at a wide gap convects downstream almost at the core speed and the flow near the wall is similar to that of single-phase wall-bounded turbulent flow. The annular flow is characterized by three different regimes depending on the clearance Reynolds number ($Re_c$) based on the core velocity and local gap size: laminar Couette flow driven by the core for $Re_c \le 600$, transitional flow for $600 < Re_c < 2500$ and turbulent flow for $Re_c \ge 2500$. The minimum pressure drop occurs at $j_w/j_o = 0.11$ in the early transition regime. For all $j_w/j_o$ considered, the negative lift force acting on the core comes from the pressure force which balances the buoyancy.

Using as a starting point conservation of momentum, a multiphase mechanical energy balance equation is derived that accounts for multiple material phases and interfaces present within a moving control volume. This balance is applied to a control volume that is anchored to a three-phase contact line as it advances continuously over the surface of a rough and chemically homogeneous and inert solid. Using semi-quantitative models for the material behaviour occurring within the control volume, an order of magnitude analysis is performed to neglect insignificant terms, producing an equation for predicting contact-angle hysteresis from a knowledge of the interface dynamics occurring around the three-phase contact line. It is shown that the viscous energy dissipation that occurs during the ‘stick–slip’ motion of the three-phase contact line, being the cause of contact-angle hysteresis on rough surfaces, can be calculated from changes in intermediate equilibrium interface states. The balance is applied to the Wenzel, Cassie–Baxter and Fakir (super-hydrophobic) wetting states, showing for the Fakir case that significant dissipation occurs during both interface advance and recede, and relating these dissipations to interfacial area changes that occur around the ‘stick–slip’ events.

We investigate theoretically and numerically the impact of an elastic sphere on a rigid wall in a viscous fluid. Our focus is on the dynamics of the contact, employing the soft lubrication model in which the sphere is separated from the wall by a thin liquid film. In the limit of large sphere inertia, the sphere bounces and the dynamics is close to the Hertz theory. Remarkably, the film thickness separating the sphere from the wall exhibits non-trivial self-similar properties that vary during the spreading and retraction phases. Leveraging these self-similar properties, we establish the energy budget and predict the coefficient of restitution for the sphere. The general framework derived here opens many perspectives to study the lubrication film in impact problems.

A wide range of environmental, energy, medical and biological processes rely on dispersive transport through complex media. Yet, because of the stagnant and opaque nature of the microscopic system, the role of disordered flow and structure in the dispersive transport of solutes remains poorly understood. Here, we use a circular porous microfluidic system to investigate the radial dispersion in porous media driven by non-Newtonian fluids with strong advection rate (or at high Péclet number) and low-to-moderate Reynolds numbers. We observe for the first time the presence of diffusion ‘blind zones’ in the microstructure for high solution injection velocities. More specifically, an in-depth analysis uncovers that the circumferential flow frame, coformed by obstacles and vortices especially the ‘twin-vortex’ with same rotation direction, is responsible for the diffusion ‘blind zones’ and transport heterogeneity. The vortices are induced by the coupling of microfluidics and porous structures, and correlated to inertial flow-induced instabilities. The trade-off between diffusion efficiency and quality/completeness with respect to the high Péclet number (or high inlet velocity) serves to enhance our comprehension of intricate fluid dynamics and affords a set of principles to aid a diverse range of practical implementations.

We investigate the influence of vortices remote from the boundary on the near-wall flow dynamics in wall-bounded flows. A vortex ring with precisely controlled local twist is introduced into the outer layer of a channel flow at a moderate Reynolds number. We find that the minimum vorticity flux for triggering the transition to turbulence is significantly reduced from the initial disturbance of an untwisted vortex ring to that of a twisted ring. In particular, the latter disturbance can cause vortex bursting in the early transitional stage. The impact of vortex bursting on the transition process is characterised by the near-wall, wall-normal velocity with the rapid distortion theory. The wall-normal velocity grows during vortex bursting, and leads to streak formation and then the transition to turbulence. The notable wall-normal velocity is induced by the large di-vorticity generated in vortex bursting. We model the growing radial component of the di-vorticity in terms of the local twist, and demonstrate that its surge is due to the generation of highly twisted vortex lines in vortex bursting. Then, we derive that the generation of the di-vorticity in the outer layer enhances the wall-normal velocity in the inner layer via the Poisson equation with the image method and the multipole expansion. Thus, we elucidate that the vortex bursting can have an effect on the transition process.

The classical paper by Lighthill (Commun. Pure Appl. Maths, vol. 109, 1952, p. 118) on the propulsion of ciliated microorganisms has become the reference against which many modern studies on swimming in low Reynolds number are compared. However, Lighthill's study was limited to propulsion in a uniform flow, whereas several biologically relevant microorganisms experience non-uniform flows. Here we propose a benchmark for ciliary propulsion in paraboloidal flows. We first consider the axisymmetric problem, with the microorganisms on the centreline of the background flow, and derive exact analytical solutions for the flow field. Our results reveal flow features, swimming characteristics and performance metrics markedly different from those generated in a uniform flow. In particular, the background paraboloidal flow introduces a Stokes quadrupole singularity at the leading-order flow field, generating vortices. Moreover, we determine the necessary conditions on the strength of the background flow for optimal power dissipation and swimming efficiency. We then consider the more general case of a microorganism off the centreline of the background flow. In this case, the squirmer experiences a paraboloidal, linear shear and uniform flows due to its position relative to the flow's centreline. Our findings show that while the linear shear flow does not affect the translational and rotational velocities of the squirmer, it does influence the velocity field and, therefore, the power dissipation.

Numerical simulations are carried out on the vortex-induced rotations of a freely rotatable rigid square cylinder in a two-dimensional uniform cross-flow. A range of Reynolds numbers between 40 and 150 and density ratios between 0.1 and 10 are considered. Results show eight different characteristic regimes, expanding the classification of Ryu & Iaccarino (J. Fluid Mech., vol. 813, 2017, pp. 482–507). New regimes include the transition and wavy rotation regimes; in the ${\rm \pi}$-limited oscillation regime we observe multipeak subregimes. Moment-generating mechanisms of these regimes and subregimes are further elucidated. A phenomenon related to the influence of density ratio is the tooth-like shape of the ${\rm \pi} /2$-limit oscillation regime observed in the regime map, which is explained as a result of the imbalance relation between the main frequencies of rotation response and the vortex shedding frequency. In addition, existence of multiple regimes and multistable states are discussed, indicating multiple stable attractive structures in phase space.

Snap-through is a buckling instability that allows slender objects, including those in plant and biological systems, to generate rapid motion that would be impossible if they were to use their internal forces exclusively. In microfluidic devices, such as micromechanical switches and pumps, this phenomenon has practical applications for manipulating fluids at small scales. The onset of this elastic instability often drives the surrounding fluid into motion – a process known as snap-induced flow. To analyse the complex dynamics resulting from the interaction between a sheet and a fluid, we develop a prototypical model of a thin sheet that is compressed between the two sides of a closed channel filled with an inviscid fluid. At first, the sheet bends towards the upstream direction and the system is at rest. However, once the pressure difference in the channel exceeds a critical value, the sheet snaps to the opposite side and drives the fluid dynamics. We formulate an analytical model that combines the elasticity of thin sheets with the hydrodynamics of inviscid fluids to explore how external pressure differences, material properties and geometric factors influence the system's behaviour. To analyse the early stages of the evolution, we perform a linear stability analysis and obtain the growth rate and the critical pressure difference for the onset of the instability. A weakly nonlinear analysis suggests that the system can exhibit a pressure spike in the vicinity of the inverted configuration.

The flow-induced vibrations (FIVs) of two identical tandem square cylinders with mass ratio m* = 3.5 at Reynolds number Re = 150 are investigated through two-dimensional direct numerical simulation (DNS) and linear stability analysis over a parameter range of spacing ratio 1.5 ≤ L* ≤ 5 and reduced velocity 3 ≤ Ur ≤ 34. Three kinds of FIV responses, namely vortex-induced vibration (VIV), biased oscillation (BO) and galloping (GA), are identified. The FIVs are then further classified into the branches of initial VIV (IV), resonant VIV (RV and RV′), flutter-induced VIV (FV), desynchronized VIV (DV), VIV developing from GA (GV), transitional state between VIV and GA (TR), BO and GA based on the characteristics of the vibration responses. The transitions among different FIV branches are examined by combining the DNS with linear stability analysis, where the transition boundaries among the VIV, BO and GA branches over the concerned parameters are identified on the branch maps. The transition from IV to RV or RV′ is found to be related to the unstable wake mode, while the FV, transiting from RV or RV′, is induced by the unstable structural factor in the wake-structure mode. The structural instability is considered as the physical origin of GA, whereas the mode competition between unstable wake and structure leads to DV, GV and TR, and thus delays the appearance of GA. The transition from DV to BO with biased equilibrium position, accompanied by the even-order harmonic frequencies, is essentially induced by the symmetry breaking bifurcation.

The parameterization of fluxes associated with representing unresolved dynamics in turbulent flows, especially in the atmosphere and ocean (which have a vast range of scales), remains a challenging task. This is especially true for Earth system models including complex biogeochemistry and requiring very long simulations. The problem of representing the dependence of the mean flux of a passive tracer in terms of the mean has a very long history; in this study, we take a somewhat different approach. We use a formalism showing that the mean flux will be a functional of the mean gradients, a formalism that can be used to calculate the structure of the functional which is non-local in both space and time. Two-dimensional turbulent simulations are used to explore the weight of nearby (in space or time) gradients. We also use stochastic velocities and iterated maps to show that the results are similar. The functional formalism provides an understanding of when non-locality needs to be considered and when a local eddy diffusivity can be a reasonably good approximation. Furthermore, the formalism provides guidance for the development of data-driven parameterizations.

Magnetic soap films under the forcing of an inhomogeneous magnetic field are governed by a wide range of interconnected physics. The study of magnetic soap films requires the development of comprehensive models to support experimental observations. In this study, the thin film approximation is applied to the Navier–Stokes equations to derive a model for the film thickness of magnetic soap films that incorporates the effects of interfacial mobility, surfactant transport and magnetite nanoparticle (NP) transport. This derived model consists of a coupled system of equations for the film thickness, interfacial velocity, interfacial surfactant concentration and magnetite NP concentration. Simulations are performed for both soap films and magnetic soap films by solving the system of equations using the Galerkin finite element method, and results are compared with experiments. Simulation results highlight that interfacial flows can dominate the rate of film thinning and that accounting for the dependence of the magnetisation on the local magnetite NP concentration can influence the predicted speed of magnetically driven flows. Furthermore, simulation results demonstrate that the model is able to predict marginal regeneration in qualitative agreement with the experiments for soap films; the model also predicts the same flow pattern as seen in the experiments for magnetic soap films. Overall, this study advances the state of soap film and magnetic soap film modelling and will contribute to acquiring control over the drainage and stability of magnetic soap films in the long term.

Direct numerical simulations of spiral Poiseuille flows in a narrow gap geometry are performed with the aim of identifying the mechanisms governing the dynamics of the axial friction coefficient. The investigation has explored a small portion of the Reynolds number–Taylor number phase space ($600 \leq Re \leq 5766$ and $1500 \leq Ta \leq 5000$), for which reference experimental results are available. The study is focused on the mechanism leading to the enhancement of the axial friction coefficient with the Taylor number when the Reynolds number is kept constant. The analysis of the spatial distribution of the Reynolds stress tensor and of the turbulent energy budget has evidenced the key role of the pressure–strain correlation in the energy transfer from the azimuthal to the axial component. The latter eventually determines the increase of the axial friction coefficient through the enhanced radial mixing of axial momentum. Data have also shown that the flow dynamics is heavily dependent on the $Ta/Re$ ratio, and different regimes develop (ranging from laminar to turbulent), each with peculiar behaviours.