We examine the evaporation-induced coalescence of two droplets undergoing freezing by conducting numerical simulations employing the lubrication approximation. When two sessile drops undergo freezing in close vicinity over a substrate, they interact with each other through the gaseous phase and the simultaneous presence of evaporation/condensation. In an unsaturated environment, the evaporation flux over the two volatile sessile drops is asymmetric, with lower evaporation in the region between the two drops. This asymmetry in the evaporation flux generates an asymmetric curvature in each drop, which results in a capillary flow that drives the drops closer to each other, eventually leading to their coalescence. This capillary flow, driven by evaporation, competes with the upward movement of the freezing front, depending on the relative humidity in the surrounding environment. We found that higher relative humidity reduces the evaporative flux, delaying capillary flow and impeding coalescence by restricting contact line motion. For a constant relative humidity, the substrate temperature governs the coalescence phenomenon and the resulting condensation can accelerate this process. Interestingly, lower substrate temperatures are observed to facilitate faster propagation of the freezing front, which, in turn, restricts coalescence.

Manipulation of small-scale particles across streamlines is the elementary task of microfluidic devices. Many such devices operate at very low Reynolds numbers and deflect particles using arrays of obstacles, but a systematic quantification of relevant hydrodynamic effects has been lacking. Here, we explore an alternative approach, rigorously modelling the displacement of force-free spherical particles in vortical Stokes flows under hydrodynamic particle–wall interaction. Certain Moffatt-like eddy geometries with broken symmetry allow for systematic deflection of particles across streamlines, leading to particle accumulation at either Faxen field fixed points or limit cycles. Moreover, particles can be forced onto trajectories approaching channel walls exponentially closely, making possible quantitative predictions of particle capture (sticking) by short-range forces. This rich, particle-size-dependent behaviour suggests the versatile use of inertia-less flow in devices with a long particle residence time for concentration, sorting or filtering.

We perform direct numerical simulations of turbulent channel flows. Secondary motions are produced by applying a streamwise-homogeneous, spanwise-heterogeneous roughness pattern of spanwise period $\Lambda _s$ to the walls of the channel; their time evolution is observed. Notice that, owing to the geometry, the secondary motions are streamwise-invariant at any instant of time, so that no spatial development is seen. Once the secondary motions reach a statistically steady state, the roughness pattern is suddenly removed, so that the secondary motions decay. The time needed for the secondary motions to vanish is then measured; in doing so, we distinguish between the streamwise-momentum pathways and the cross-sectional circulatory motions that compose the secondary motions. Larger values of $\Lambda _s$ are generally associated with a longer time scale for the decay of the momentum pathways, although this might not hold true for $\Lambda _s/h\gt 4$ (where $h$ is the channel half-height). The value of such a time scale for the circulatory motions, instead, saturates for $\Lambda _s/h \geqslant 2$; this may be related to the observed spatial confinement of said circulatory motions. For specific values of $\Lambda _s$ ($2 \leqslant \Lambda _s/h \leqslant 4$), the volume-averaged energy associated with the momentum pathways undergoes an unexpected transient growth with respect to its value at the beginning of the decay. This might indicate that structures of such a specific size are able to self-sustain as postulated by Townsend (The Structure of Turbulent Shear Flow, 2nd edition, 1976, ch. 7.19); the evidence we gather in this respect is however inconclusive. Finally, the present data suggest that most of the energy of the momentum pathways is produced by the circulatory motions transporting the mean (spanwise-averaged) velocity.

Particle image velocimetry is used to study the control of swirl momentum, delivered through an orifice formed by a physically rotating tube of finite length, relevant to the evolution of vortex rings produced at a Reynolds number ${Re}\approx 1000$ based on the average discharge velocity, for swirl numbers ${S} \in [0, 1]$. Experiments without discharge, reinforced with complimentary numerical predictions, reveal the presence of an intriguing secondary flow pattern in the rotating tube, preventing attainment of a solid-body-like swirl distribution. Nevertheless, it is found that fully established rings produced in this way, following discharge once conditions in the tube have reached a steady state, exhibit similar characteristics to rings formed by an otherwise solid-body rotating initial condition as explored computationally by Ortega-Chavez et al. (2023, J. Fluid Mech. 967, A16). Namely, opposite-signed vorticity forms due to vortex tilting, which subsequently interacts with the ring, promoting vorticity cancellation and vortex ring breakdown. A key feature of the experimental work is that partially established vortex rings, produced before a steady-state rotating tube condition is reached, show unique characteristics. Their creation, a short time after the onset of tube rotation: (i) facilitates more efficient delivery of swirl momentum to the vortex core area; (ii) maintains a low level of swirl in the ring bubble’s central region which would otherwise promote the formation of opposite-signed vorticity and vortex breakdown.

We investigate the natural oscillations of sessile drops with a central trapped bubble on a plane using linear potential flow theory, considering both free and pinned contact lines. The system is governed by the contact angle $\alpha$ and the ratio $\tau$ of inner to outer contact line radii. For bubble-containing (BC) hemispherical drops with free contact lines (referred to as free BC semi-drops), the modes mirror half of those in concentric spherical BC drops due to plane symmetry. These modes are labelled ‘plus’ (with greater inner surface deformation) and ‘minus’ (with greater outer surface deformation). As $\tau \to 0$, minus modes converge to those of bubble-free drops. Results show that varying $\alpha$ from $90^\circ$ or pinning the contact line in free BC semi-drops alters the topology of spectral lines, turning original crossings of spectral lines between minus and plus modes into avoided crossings. This shift causes minus and plus modes to form spectral trends with avoided crossings, maintaining their original spectral shapes. In an avoided crossing, two coupled modes cannot be classified as plus or minus due to their comparable inner and outer surface deformations, resulting in mode beating when both are excited, as confirmed by our direct numerical simulations. This study on the impact of inner bubbles on the spectrum may help in predicting bubble size in opaque sessile drops.

Inertia–gravity waves are scattered by background flows as a result of Doppler shift by a non-uniform velocity. In the Wentzel–Kramers–Brillouin regime, the scattering process reduces to a diffusion in spectral space. Other inhomogeneities that the waves encounter, such as density variations, also cause scattering and spectral diffusion. We generalise the spectral diffusion equation to account for these inhomogeneities. We apply the result to a rotating shallow-water system, for which height inhomogeneities arise from velocity inhomogeneities through geostrophy, and to the Boussinesq system for which buoyancy inhomogeneities arise similarly. We compare the contributions that height and buoyancy variations make to the spectral diffusion with the contribution of the Doppler shift. In both systems, we find regimes where all contributions are significant. We support our findings with exact solutions of the diffusion equation and with ray tracing simulations in the shallow-water case.

We investigate fully developed turbulent flow in curved channels to explore the interaction between turbulence and curvature-driven coherent structures. By focusing on two cases of mild and strong curvature, we examine systematically the effects of the Reynolds number through a campaign of direct numerical simulations, spanning flow regimes from laminar up to the moderately high Reynolds number – based on bulk velocity and channel height – of $87\,000$. Our analysis highlights the influence of curvature on the friction coefficient, showing that flow transition is anticipated by concave curvature and delayed by convex curvature. In the case of mild curvature, a frictional drag reduction compared with plane channel flow is found in the transitional regime. Spectral analysis reveals that the near-wall turbulence regeneration cycle is maintained in mildly curved channels, while it is absent or severely inhibited on the convex wall of strongly curved channels. Streamwise large-scale structures resembling Dean vortices are found to be weakly dependent on the Reynolds number and strongly affected by curvature: increasing curvature shifts these vortices towards the outer wall and reduces their size and coherence, limiting their contribution to streamwise velocity fluctuations and momentum transport. In the case of strong curvature, spanwise large-scale structures are also detected. These structures are associated with large pressure fluctuations and the suppression of turbulent stresses near the convex wall, where a region with negative turbulence production is observed and characterised via quadrant analysis.

We introduce a wall model for large-eddy simulation (WMLES) applicable to rough surfaces with Gaussian and non-Gaussian distributions for both the transitionally and fully rough regimes. The model is applicable to arbitrary complex geometries where roughness elements are assumed to be underresolved, i.e. subgrid-scale roughness. The wall model is implemented using a multi-hidden-layer feedforward neural network, with the mean geometric properties of the roughness topology and near-wall flow quantities serving as input. The optimal set of non-dimensional input features is identified using information theory, selecting variables that maximize information about the output while minimizing redundancy among inputs. The model also incorporates a confidence score based on Gaussian process modelling, enabling the detection of potentially low model performance for untrained rough surfaces. The model is trained using a direct numerical simulation (DNS) roughness database comprising approximately 200 cases. The roughness geometries for the database are selected from a large repository through active learning. This approach ensures that the rough surfaces incorporated into the database are the most informative, achieving higher model performance with fewer DNS cases compared with passive learning techniques. The performance of the model is evaluated both apriori and aposteriori in WMLES of turbulent channel flows with rough walls. Over 550 channel flow cases are considered, including untrained roughness geometries, roughness Reynolds numbers and grid resolutions for both transitionally and fully rough regimes. Our rough-wall model offers higher accuracy than existing models, generally predicting wall shear stress within an accuracy range of 1%–15 %. The performance of the model is also assessed on a high-pressure turbine blade with two different rough surfaces. We show that the new wall model predicts the skin friction and the mean velocity deficit induced by the rough surface on the blade within 1%–10 % accuracy except the region with transition or shock waves. This work extends the building-block flow wall model (BFWM) introduced by Lozano-Durán & Bae (2023. J. Fluid Mech. 963, A35) for smooth walls, expanding the BFWM framework to account for rough-wall scenarios.

Large-eddy simulations are analysed to determine the influence of suspended canopies, such as those formed in macroalgal farms, on ocean mixed layer (OML) deepening and internal wave generation. In the absence of a canopy, we show that Langmuir turbulence, when compared with wind-driven shear turbulence, results in a deeper OML and more pronounced internal waves beneath the OML. Subsequently, we examine simulations with suspended canopies of varying densities located in the OML, in the presence of a background geostrophic current. Intensified turbulence occurs in the shear layer at the canopy’s bottom edge, arising from the interaction between the geostrophic current and canopy drag. Structures resembling Kelvin–Helmholtz (KH) instability emerge as the canopy shear layer interacts with the underlying stratification, radiating internal waves beneath the OML. Both intensified turbulence and lower-frequency motions associated with KH-type structures are critical factors in enhancing mixing. Consequently, the OML depth increases by up to a factor of two compared with cases without a canopy. Denser canopies and stronger geostrophic currents lead to more pronounced KH-type structures and internal waves, stronger turbulence and greater OML deepening. Additionally, vertical nutrient transport is enhanced as the OML deepens due to the presence of the canopy. Considering that the canopy density investigated in this study closely represents offshore macroalgal farms, these findings suggest a mechanism for passive nutrient entrainment conducive to sustainable farming. Overall, this study reveals the intricate interactions between the suspended canopy, turbulent mixing and stratification, underscoring their importance in reshaping OML characteristics.

A collection of secondary instability calculations in streaky boundary layers is presented. The data are retrieved from well-resolved numerical simulations of boundary layers forced by free-stream turbulence (FST), considering different geometries and FST conditions. The stability calculations are performed before streak breakdown, taking place at various $Rey_x$ the Reynolds number based on the streamwise coordinate. Despite the rich streak population of various sizes, it is found that breaking streaks have similar aspect ratios, independently of the streamwise position where they appear. This suggests that wider streaks will break down further downstream than thinner ones, making the appearance of secondary instabilities somewhat independent of the streak’s wavelength. Moreover, the large difference in the integral length scale among the simulations suggests that this aspect ratio is also independent of the FST scales. An explanation for this behaviour is provided by showing that these breaking streaks are in the range of perturbations that can experience maximum transient growth according to optimal disturbance theory. This could explain why, at a given streamwise position, there is a narrow spanwise wavelength range where streak breakdown is more likely to occur.

Wave propagation in channels with area changes is a topic of significant practical interest that involves a rich set of coupled physics. While the acoustic wave problem has been studied extensively, the shock propagation problem has received less attention. In addition to its practical significance, this problem also introduces deep fundamental issues associated with how energy in propagating large-amplitude disturbances is redistributed upon interaction with inhomogeneities. This paper presents a study of shock scattering and entropy and vorticity coupling for shock wave propagation through discrete area changes. It compares results from computational fluid dynamics to those of one-dimensional quasi-steady calculations. The solution space is naturally divided into five ‘regimes’ based upon the incident shock strength and area ratio. This paper also presents perturbation methods to quantify the dimensionless scaling of physical effects associated with wave reflection/transmission and energy transfer to other disturbances. Finally, it presents an analysis of the ‘energetics’ of the interaction, quantifying how energy that initially resides in dilatational disturbances and propagates at the shock speed is redistributed into finite-amplitude reflected and transmitted waves as well as convecting vortical and entropy disturbances.

The formation process of a vortex pair generated by a two-dimensional starting jet has been investigated numerically over a range of Reynolds numbers from 500 to 2000. The effects of stroke ratio and nozzle configuration are examined. Only a single vortex pair can be observed in the vorticity field generated by small stroke ratios less than 10 while the leading vortex pair formed by larger stroke ratios eventually disconnects from the trailing jet. The formation numbers (13.6 and 9.3) for a straight nozzle and an orifice nozzle have been identified by the circulation criterion and they are further analysed by four other criteria. Using the contraction coefficient, formation numbers can be transformed into a universal value at about 16.5 for both nozzles. The effect of Reynolds number on the formation number is found to be within 12 % for parallel flow cases but it will increase up to 27 % for non-parallel flow cases due to shear-layer instability. A modified contraction-based slug model is proposed, and it can accurately predict the total invariants (e.g. circulation, hydrodynamic impulse and kinetic energy) shedding from the nozzle edge. Analytical estimation of the formation number is further conducted by matching the predicted total invariants to the Pierrehumbert model of steady vortex pairs. By assuming that pinch-off starts when the vortex pair achieves the steady state, two analytical models are proposed in terms of vortex impulse and translational velocity. The latter appears to be more appropriate to predict the formation number for two-dimensional flows.

Turbulent transonic buffet is an aerodynamic instability causing periodic (albeit, often irregular) oscillations of lift/drag in aerospace applications. Involving complex coupling between inviscid and viscous effects, buffet is characterised by shock wave oscillations and flow separation/reattachment. Previous studies have identified both two-dimensional (2-D) chordwise shock-oscillation and three-dimensional (3-D) buffet-/stall-cell modes. While the 2-D instability has been studied extensively, investigations of 3-D buffet have been limited to only low-fidelity simulations or experiments. Due to computational cost, almost all high-fidelity studies to date have been limited to narrow span-widths around 5 % of aerofoil chord length (aspect ratio, ), which is insufficiently wide to observe large-scale three-dimensionality. In this work, high-fidelity simulations are performed up to , on an infinite unswept NASA Common Research Model (CRM) wing profile at $Re=5\times 10^{5}$. At , intermittent 3-D separation bubbles are observed at buffet conditions. While previous Reynolds-averaged Navier–Stokes (RANS)/stability-based studies predict quasi-simultaneous onset of 2-D- and 3-D-buffet, a case that remains essentially 2-D is identified here. Strongest three-dimensionality was observed near low-lift phases of the buffet cycle at maximum flow separation, reverting to essentially 2-D behaviour during high-lift phases. Buffet was found to become 3-D when extensive mean flow separation was present. At , multiple 3-D separation bubbles form in a spanwise wavelength range $\lambda =1c$ to $1.5c$. Spectral proper orthogonal decomposition (SPOD) was applied to analyse the spatio/temporal structure of 3-D buffet-cells. In addition to the 2-D chordwise shock-oscillation mode (Strouhal number $St \approx 0.07-0.1$), 3-D modal structures were observed at the shock wave/boundary layer interaction at $St \approx 0.002-0.004$.

Experiments are presented to explore the non-axisymmetric instabilities of spreading films of aqueous suspensions of Carbopol and Xanthan gum floating on a bath of perfluoropolyether oil. The experimental observations are compared against theoretical predictions exploiting a shallow-film model in which the viscoplastic rheology is captured by the Herschel–Bulkley constitutive law. With this model, we construct axisymmetric base states that evolve from the moment that the film floats onto the bath, out towards long times at which spreading becomes self-similar, and then test their linear stability towards non-axisymmetric perturbations. In the geometry of a thinning expanding film, we find that shear thinning does not drive a loss of axisymmetry at early times (when the degree of expansion is small), but when the film has expanded in radius by a factor of two or so, shear-thinning hoop stresses drive non-axisymmetric instabilities. Unstable modes possess relatively low angular wavenumber, and the loss of symmetry is not particularly dramatic. When the oil in the bath is replaced by salty water, the experiments are completely different, with dramatic non-axisymmetric patterns emerging from interfacial effects.

We report the first measurement of turbulent mixing developing from the convergent Richtmyer–Meshkov (RM) instability using time-resolved planar laser-induced fluorescence in a semi-annular convergent shock tube. A membraneless yet sharp interface with random short-wavelength perturbations, but controllable long-wavelength perturbations, is created by an automatically retractable plate, enhancing the reproducibility and reliability of RM turbulence experiments. The cylindrical air/SF$_6$ interface formed is first subjected to a convergent shock, then to its reflected shock and subsequently transits to turbulent mixing. It is found that the mixing width after reshock has a linear growth rate more than twice the rate in planar geometry. Also, the mixing width does not present power-law growth at late stages as in a planar geometry. However, the scalar spectrum and transition criterion obtained are similar to their planar counterparts. These findings indicate that the geometric constraint greatly affects the large scales of the flow, while having a weaker effect on the small scales. It is also found that the reflected shock significantly increases both scale separation and Reynolds number, explaining the rapid transition to turbulence following reshock.

The transport of particles in elastoviscoplastic (EVP) fluids is of significant interest across various industrial and scientific domains. However, the physical mechanisms underlying the various particle distribution patterns observed in experimental studies remain inadequately understood in the current literature. To bridge this gap, we perform interface-resolved direct numerical simulations to study the collective dynamics of spherical particles suspended in a pressure-driven EVP duct flow. In particular, we investigate the effects of solid volume fraction, yield stress, inertia, elasticity, shear-thinning viscosity, and secondary flows on particle migration and formation of plug regions in the suspending fluid. Various cross-streamline migration patterns are observed depending on the rheological parameters of the carrier fluid. In EVP fluids with constant plastic viscosity, particles aggregate into a large cluster at the duct centre. Conversely, EVP fluids with shear-thinning plastic viscosity induce particle migration towards the duct walls, leading to formation of particle trains at the corners. Notably, we observe significant secondary flows ($O(10^{-2})$ compared to the mean velocity) in shear-thinning EVP suspensions, arising from the interplay of elasticity, shear-thinning viscosity and particle presence, which further enhances corner-ward particle migration. We elucidate the physical mechanism by which yield stress augments the first normal stress difference, thereby significantly amplifying elastic effects. Furthermore, through a comprehensive analysis of various EVP suspensions, we identify critical thresholds for elasticity and yield stress necessary to achieve particle focusing at the duct corners.

In binary mixtures, the lifetimes of surface bubbles can be five orders of magnitude longer than those in pure liquids because of slightly different compositions of the bulk and the surfaces, leading to a thickness-dependent surface tension of thin films. Taking advantage of the resulting simple surface rheology, we derive the equations describing the thickness, flow velocity and surface tension of a single liquid film, using thermodynamics of ideal solutions and thin film mechanics. Numerical resolution shows that, after a first step of tension equilibration, the Laplace-pressure-driven flow is associated with a flow at the interface driven by an induced Marangoni stress. The resulting parabolic flow with mobile interfaces in the film further leads to its pinching, eventually causing its rupture. Our model paves the way for a better understanding of the rupture dynamics of liquid films of binary mixtures.

Linear and nonlinear mechanisms governing the growth of second-mode waves are analysed using a newly derived disturbance energy conservation equation that highlights the physical processes responsible for fluctuation energy production, flux-based transport and destruction. Axisymmetric direct numerical simulations (DNS) data from a Mach 6 hypersonic boundary layer, simulated over a $3^\circ$ half-angle sharp cone at zero angle of attack, is used as a reference. A Legendre polynomial-based forcing methodology is used to trigger transition in the DNS over a range of various amplitude levels and different forcing frequency content. Closure of the disturbance energy budgets is demonstrated numerically using the DNS data. The terms responsible for the amplification of the disturbance are identified, and nonlinear attenuation effects are discussed. We show that the interaction between the entropy/velocity fluctuations and the base temperature gradient governs the second-mode growth in the linear regime. Energy production occurs in the critical layer due to non-isentropic processes and accumulates in acoustic form below the relative sonic line through downward transport. At higher forcing amplitudes, nonlinear spectral broadening is observed, with simultaneous thermoviscous diffusion attenuating the disturbance energy growth. This effect is responsible for the non-monotonic streamwise variation of the wall-pressure spectrum. Phase speed and growth rate analyses, informed by linear stability theory (LST), reveal wave steepening effects preceding this nonlinear attenuation effect. The disturbance energy is observed to match the LST predictions at lower forcing amplitudes, deviating, as expected, at higher amplitudes.

Skin-friction drag reduction (DR) in a turbulent boundary layer (TBL) using plasma-generated streamwise vortices (PGSVs) is governed by plasma-induced spanwise wall-jet velocity $W$, the distance $L$ between the positive electrodes of two adjacent plasma actuators (PAs) and the friction Reynolds number $Re_\tau$. It is found experimentally that DR increases logarithmically with the growing maximum spanwise mean velocity $\overline {W}_{max}^+$ but decreases with rising $L^+$ and $Re_\tau$, where superscript ‘+’ denotes normalization by the inner scales. It is further found from theoretical and empirical scaling analyses that the dimensionless drag variation $\Delta F = g_1 (\overline {W}_{max}^+, L^+, {Re_\tau })$ may be reduced to $\Delta F = g_2 (\xi )$, where $g_1$ and $g_2$ are different functions and the scaling factor $\xi = [k_{2} \log _{10} (k_{1} \overline {W}_{max }^{+} ) ] / (L^{+} Re_{\tau } )$ ($k_{2}$ and $k_{1}$ are constants) is physically the circulation of the PGSVs. Discussion is conducted based on $\Delta F = g_2 (\xi )$, which provides important insight into the physics of TBL control based on PAs.

Three-dimensional eddy-resolving simulations are used to study the structure of turbulent boundary layers developing in an open channel where an array of sparse roughness elements in the form of partially burrowed mussels is placed on the smooth, flat channel bed starting at a certain streamwise location. The identical mussels are oriented with their major axis parallel to the mean flow. Their positions are randomised while making sure that the mussels are close to uniformly distributed inside the array. The turbulent flow approaching the leading edge of the array is fully developed. The protruding parts of the mussels play the role of sparse roughness elements and generate a rough-bed, internal boundary that is characterised by zero net flow exchange but non-zero local flow exchange due to active filtering. A double-averaging technique is used to obtain an equivalent, width- and time-averaged, boundary layer over a ‘flat’ rough surface containing no mussels. The paper discusses the effects of varying the mussel array density, protruding height of the mussels and filtering discharge on the spatial growth of the two-dimensional boundary layer. With proper scaling, the profiles of the (double-averaged) streamwise velocity are close to self-similar inside the inertial layer (e.g. for h $\lt$ z $\lt$ $\delta$, where h is the height of the protruding part of the mussels and $\delta$ is the boundary layer thickness) starting some distance from the leading edge of the array. The scaled turbulent kinetic energy and concentration profiles associated with the scalar advected through the excurrent syphons are also found to be self-similar above the vertical location where the maximum value is reached. An analytical model containing three subzones is proposed for the streamwise velocity in between the bed (z $=$ 0) and z $=$ $\delta$. The velocity profile inside the inertial region contains a law-of-the-wall component supplemented by a law-of-the-wake component. The scaling coefficient of the law-of-the-wake component is found to be larger than typical values used to describe velocity variation in turbulent boundary layers developing in a surrounding flow with close-to-uniform free stream velocity. The equivalent roughness height for this particular type of boundary layers developing over sparse roughness elements increases monotonically with h and the mussel array density, $\rho$N. The paper also discusses the effect of the mussel bed density on the average refiltration fraction and the phytoplankton removal efficiency of the mussel bed.