An experimental investigation of separation bubble shaped control bumps for oblique shock wave–boundary-layer interactions was performed in two supersonic wind tunnel facilities at Mach 2.5 and 2, with incident shock deflection angles of $8^\circ$ and $12^\circ$, respectively, and momentum thickness Reynolds numbers of approximately $1.5 \times 10^4$. Shock control bumps were designed to replicate the time-averaged separation bubble shape, and were placed onto the floor in the separation location. This resulted in almost complete elimination of flow separation. There was also a marked improvement in the downstream boundary-layer state. A low-frequency bubble breathing oscillation was identified in the baseline interaction using high-speed shadowgraphy and particle image velocimetry measurements. This oscillation was strongly suppressed in the controlled interactions. Velocity fluctuations in the downstream boundary layer were also significantly reduced. We propose that the key mechanism by which flow separation is eliminated is by breaking down the overall pressure rise into smaller steps, each of which is below the separation threshold. A key feature is the bump crest expansion fan, located near to where the incident shock terminates, which negates the shock induced pressure jump. Thus, the precise bump geometry is critical for control efficacy and should be designed to manage these pressure rise steps as well as the expansion fan strength and location with respect to the incident shock wave. The length of the bump faces must also be sufficiently long for the boundary layer to recover between successive adverse pressure jumps.

Spectral turbulence models commonly used in the design and certification of wind turbines have only been validated at heights up to 70 m in the atmosphere, but many offshore wind turbines now operate at heights above 150 m. Moreover, there is a lack of measurement data on the spatial structure of turbulence at such heights in the marine atmospheric boundary layer (MBL). Consequently, it is uncertain whether these turbulence models are valid for the design of tall offshore wind turbines. To fill this gap, we present measurements of one-point auto-spectra and two-point spectral coherence at heights of 150–250 m and lateral separations up to 241 m providing lateral coherence of turbulence in the MBL that has never been measured before for these heights and separations. Five light detection and ranging (lidar) instruments were deployed on the west coast of Denmark, and we reconstructed the along-wind and cross-wind components at the lidar beam intersection points. The measurements were compared with the theoretical predictions of auto-spectra and lateral coherence from the Mann model and its extension, the Syed–Mann model. The latter models turbulence down to frequencies of 1 h$^{-1}$ through the $-5/3$ scaling observed in the mesoscale range. The results show that the Mann model did not compare well with the measurements under stable and near-neutral conditions. On the other hand, the Syed–Mann model predicted the lateral coherence for a range of different conditions. However, the lateral coherence was under predicted in about $8\,\%$ of the data, possibly due to gravity waves. We believe that the high coherence from mesoscale turbulence at these heights can influence the loads on floating wind turbines and large offshore wind farms.

The interaction of an object with an unsteady flow is non-trivial and is still far from being fully understood. When an aerofoil or hydrofoil, for example, undergoes time-dependent motion, nonlinear flow phenomena such as dynamic stall can emerge. The present work experimentally investigates the interaction between a hydrofoil and surface gravity waves. The waves impose periodic fluctuations of the velocity magnitude and orientation, causing a steadily translating hydrofoil to be susceptible to dynamic stall at large wave forcing amplitudes. Simultaneous measurement of both the forces acting on the hydrofoil and the flow around it by means of particle image velocimetry (PIV) are performed, to properly characterise the hydrofoil–wave interaction. In an attempt at alleviating the impact of the flow unsteadiness via passive flow control, a bio-inspired tubercle geometry is applied along the hydrofoil leading edge. This geometry is known to delay stall in steady cases but has scarcely been studied in unsteady flow conditions. The vortex structures associated with dynamic stall are identified, and their trajectories, dimension and strength characterised. This analysis is performed for both straight- and tubercled-leading-edge geometries, with tubercles found to qualitatively modify the flow behaviour during dynamic stall. In contrast to previous studies, direct measurements of lift do not evidence any strong modification by tubercles. Drag-driven horizontal force fluctuations, however, which have not previously been measured in this context, are found to be strongly attenuated. This decrease is quantified and a physical model based on the flow observations is finally proposed.

Triply periodic minimal surfaces (TPMS)-based media (a type of metamaterial) are defined by mathematical expressions, which are amenable to additive manufacturing, and are finding increasing practical applications owing to their porous nature. We present experimental pressure drop measurements for a range of velocities spanning laminar to turbulent regimes for three TPMS geometries – gyroid, primitive and body-centred cubic (BCC) – with different porosity, unit cell length and surface finish. Dimensional Darcy and Forchheimer permeabilities are estimated via quadratic fitting for the gyroid geometry, which closely resembles random packed porous media. Subsequently, the non-dimensional drag (${\kern-0.5pt}f$) is plotted against Reynolds number ($Re$) yielding distinct curves for each case. The lack of collapse stems from varying definitions of pore diameter, complicating comparisons across porous media (not just TPMS). Therefore, a method is developed to estimate an equivalent hydraulic diameter $d_{{H\hbox{-}\textit{equ}}}$ from pressure drop data by matching the laminar drag $f$ of packed spheres via the Ergun equation, allowing the collapse of all porous media $f-Re$ curves in the laminar regime. The value of $d_{ {H\hbox{-}\textit{equ}}}$ is related to the ‘true’ Darcy permeability defined strictly in the linear regime (unlike permeability from quadratic fitting). We observe an approximate linear relationship between the $d_{ {H\hbox{-}\textit{equ}}}$ and the hydraulic diameter for self-similar TPMS configurations. The common basis of $d_{ {H\hbox{-}\textit{equ}}}$ allows intercomparison of TPMS geometries, and shows that BCC achieves significant drag reduction compared with packed spheres in the turbulent regime partially because of their open tube-like structure, whereas some configurations show drag increase. Although gyroid can be represented using the traditional quadratic drag law, primitive and BCC show an increase in $f$ with increasing $Re$ immediately before transitioning to fully turbulent regime – akin to rough-wall pipe flows, likely owing to their periodic streamwise elongated open structures.

A linear theory for unsteady aerodynamic effects of the actuator line method (ALM) is developed. This theory is validated using two-dimensional ALM simulations, where we compute the unsteady lift generated by the plunging and pitching motion of a thin aerofoil in uniform flow, comparing the results with Theodorsen’s theory. This comparison elucidates the underlying characteristics and limitations of ALM when applied to unsteady aerodynamics. Numerical simulations were conducted across a range of chord lengths and oscillation frequencies. Comparison of ALM results with theoretical predictions shows consistent accuracy, with all Gaussian parameter choices yielding accurate results at low reduced frequencies. Furthermore, the study indicates that selecting a width parameter ratio of $\varepsilon /c$ (the Gaussian width parameter over the chord length) between 0.33 and 0.4 in ALM yields the closest alignment with analytical results across a broader frequency range. Additionally, a proper definition of angle of attack for a pitching aerofoil is shown to be important for accurate computations. These findings offer valuable guidance for the application of ALM in unsteady aerodynamics and aeroelasticity.

We investigate the influence of shear-thinning and viscoelasticity on turbulent drag reduction in lubricated channel flow – a configuration where a thin lubricating layer of non-Newtonian fluid facilitates the transport of a primary Newtonian fluid. Direct numerical simulations are performed in a channel flow driven by a constant mean pressure gradient at a reference shear Reynolds number $\textit{Re}_\tau = 300$. The interface between the two fluid layers is characterised by Weber number $\textit{We} = 0.5$. The fluids are assumed to have matched densities. In addition to a single-phase reference case, we analyse four configurations: a Newtonian lubrication layer, a shear-thinning Carreau fluid layer, a shear-thinning and viscoelastic FENE-P fluid layer, and a purely viscoelastic FENE-CR fluid layer. Consistent with previous findings (Roccon et al. 2019, J. Fluid Mech., vol. 863, R1), surface tension is found to induce significant drag reduction across all cases. Surprisingly, variations in the lubricating layer viscosity do not yield noticeable drag-reducing effects: the Carreau fluid, despite its lower apparent viscosity, behaves similarly to the Newtonian case. In contrast, viscoelastic effects lead to a further reduction in drag, with both the FENE-P and FENE-CR fluids demonstrating enhanced drag-reducing capabilities.

Previous literature has shown that the introduction of homogeneous perforation on plates and cylinders decreases aerodynamic drag. Here, it is shown that the opposite is true for a sphere; drag can increase with porosity. Hollow porous spheres exposed to a uniform free stream are studied experimentally using force and flow field measurements. The parameter space encompasses moderate to high Reynolds numbers ($5 \times 10^4 \leq \textit{Re} \leq 4 \times 10^5$) and porosities ranging from $0\,\%$ to $80\,\%$. The main conclusion is that drag increases with porosity, at super-critical Reynolds numbers, for all studied porosities. At low porosities (less than $9\,\%$), the effect of porosity on drag can be explained by shifts in the separation point. At higher porosities the drag increase cannot be explained by separation shifts, and instead is explained by two competing forms of kinetic energy dissipation: (i) shear on the macro-scale of the body, and (ii) hole losses from flow through the pores. The former generally decreases with porosity, as bleeding flow passing through the body decreases the characteristic velocity difference in the body-scale wake. In a sphere, hole losses increase with porosity sufficiently fast to overcome decreasing body-scale shear losses, in contrast to plates and cylinders where this is not the case. Relatively weak wake vortex structures, and associated low drag coefficient at zero porosity, for a sphere reduce the impact of wake bleeding. Moreover, fluid entering the fore of a sphere can exit perpendicular to the free stream, further reducing wake bleeding while still contributing to hole losses.

This work aims to complement the description of the atomisation process in a typical commercial pressure-swirl atomiser. Conventional characterisation focuses on the final spray, where established experimental techniques allow for measuring spherical droplets in a dilute regime. However, the early stages of atomisation involve distorted liquid structures with complex interface morphology that challenge both experimental and numerical approaches. While numerical simulations with interface-capturing methods have provided access to this region, they currently remain computationally prohibitive to follow the atomisation process until the formation of the final spherical droplets. To characterise the evolving interface morphology, we propose analysing the curvature distribution obtained from both simulations and two-photon laser-induced fluorescence (2P-LIF) imaging. This curvature-based methodology, recently developed to characterise numerical sprays (Palanti et al. Intl J. Multiphase Flow 147, 2022, 103879; Ferrando et al. Atomiz. Sprays 33, 2023, 1–28), is here extended to experimental data. Both approaches are compared with available phase Doppler anemometry (PDA) measurements performed further downstream on spherical droplets. The morphological evolution of the atomising spray is interpreted through curvature statistics, which provide a unified framework applicable to all atomisation stages. When applied to spherical droplets, the curvature distribution recovers the conventional drop size distribution, linking early interface deformation to the final spray structure. The birth of this final drop size distribution can thus be observed by comparing the three approaches – numerical simulation limited to the early stage of atomisation, curvature derived from 2P-LIF images limited to two-dimensional (2-D) contour analysis, and PDA measurements of the dilute spray. The results show that curvature properties evolve in a way that can be directly representative of the final spray even at early atomisation stages.

A lattice Boltzmann method is adopted to investigate the breakup of surfactant-free and surfactant-laden droplets in both regular and irregular T-junction microchannels. During droplet neck contraction, the neck thinning shifts from inertia dominated to interfacial tension dominated, causing spontaneous rapid neck collapse due to Rayleigh–Plateau instability. For the regular rectangular microchannels, we find that the prerequisite for the spontaneous breakup of a surfactant-free droplet is that the local capillary pressure in the triggering area exceeds the Laplace pressure difference between the inside and outside of the droplet neck. Results show that the critical neck thickness $\delta _\textit{cr}^{*}$ for the droplet spontaneous breakup increases with increasing height-to-width ratio $\chi$ of the microchannel in both surfactant-free and surfactant-laden systems. The presence of surfactants decreases $\delta _\textit{cr}^{*}$ at the identified $\chi$, while the surfactant effects are gradually enhanced as $\chi$ increases. Subsequently, a constriction section is incorporated into the upper microchannel wall to establish an irregular microchannel. As constriction depth (length) increases, $\delta _\textit{cr}^{*}$ linearly decreases (increases) in the surfactant-free system, while $\delta _\textit{cr}^{*}$ exponentially decreases (linearly increases) in the surfactant-laden system. Four empirical formulas are proposed to predict the values of $\delta _\textit{cr}^{*}$ under varying constriction depths and lengths in the two systems.