This work reports an experimental study of the turbulent entrainment into the planar wake of a circular cylinder, exposed to various turbulent backgrounds, from the near- to the far-field. The background turbulence features independently varying turbulence intensity and integral length scale, thereby rendering different turbulent/turbulent interfaces (TTIs) between the background and the primary flow (wake). Combined, simultaneous particle image velocimetry and planar laser induced fluorescence measurements were conducted to quantify the entrainment characteristics across these various TTIs at an inlet Reynolds number of 3800. The primary focus was on understanding how turbulent entrainment evolves spatially in conjunction with the rapid development of the large-scale coherent vortices in the planar wake, and how such evolution is affected by the background turbulence. It is found that TTIs can establish two layers when the background turbulence is sufficiently intense, which distinguishes TTIs from the turbulent/non-turbulent interface (TNTI). The two layers are underpinned by different physical mechanisms but have the same thickness and appear to scale with the local Kolmogorov length scale after the wake spreading transition position (Chen & Buxton, J. Fluid Mech., vol. 969, 2023, A4). It is also found that the probability density functions of the entrainment velocity for both TTIs and a TNTI display power law tails, which are associated with extremely large entrainment velocities occurring more frequently than for a Gaussian process. These intermittent, extreme entrainment velocities make a remarkable contribution to the mean entrainment velocity, particularly in the near wake, which leads to a much higher mean entrainment velocity than farther downstream, for both a TNTI and the TTIs. Conditionally averaged analysis reveals that these extreme events of the entrainment velocity are directly associated with intense enstrophy structures close to the interface.

Understanding the solutal convection is a crucial step towards accurate prediction of CO$_2$ sequestration in deep saline aquifers. In this study, pore-scale resolved direct numerical simulations (DNS) are performed to analyse the scaling laws of the solutal convection in porous media. The porous media studied are composed of uniformly distributed square or circular elements. The Rayleigh numbers in the range $426 \le Ra \le 80\,000$, the Darcy numbers in the range $1.7\times 10^{-8} \le Da \le 8.8\times 10^{-6}$ and the Schmidt numbers in the range $250 \le Sc \le 10^4$ are considered in the DNS. The main time, length and velocity scales affecting the solutal convection are classified as the diffusive scales ($t_I$, $l_I$ and $u_I$), the convective scales ($t_{II}$, $l_{II}$ and $u_{II}$) and the shut-down scales ($t_{III}$, $l_{III}$ and $u_{III}$). These scales determine the pore-scale Rayleigh number $Ra_K$ and the Rayleigh number $Ra$. Based on the DNS results, the scaling laws for the transient dissolution flux are proposed in the different regimes of convection. It is found that the onset time of convection ($t_{oc}$) and the period of the flux-growth regime ($\Delta t_{fg}$) vary linearly with the convective time scale $t_{II}$. The merging of the original plumes and the re-initiation of the new plumes start in the same period, resulting in the merging re-initiation regime. Horizontal dispersion plays an important role in plume merging. The dissolution flux $F$ and the time since the onset of convection $t^{\ast }$ have an $F / u_{II} \sim (t^{\ast }/ t_{II})^{-0.2}$ scaling in the later stage of the merging re-initiation regime. This scaling is caused by the merging of the low-wavenumber plumes. It becomes more pronounced with decreasing porosity and leads to the nonlinear relationship between the Sherwood number $Sh$ and $Ra$ when the domain is not high enough for the plumes to fully develop. According to the DNS results, a regime is expected that is independent of both $Ra$ and $Ra_K$, while the dimensionless constant flux $F_{cf}/u_{II}$ in this regime decreases with decreasing porosity. When the mean solute concentration reaches approximately 30 %, convection enters the shut-down regime. For large $Ra$, the dimensionless flux in the shut-down regime follows the scaling law $F/u_{III}\sim (t/t_{III})^{-1.2}$ at large porosity ($\phi =0.56$) and $F/u_{III}\sim (t/t_{III})^{-1.5}$ at small porosity ($\phi =0.36$ or $0.32$). The study shows the significant pore-scale effect on the convection. The DNS cases in this study are mainly for simplified geometries and large $Ra_K$. This can lead to uncertainties of the obtained scaling coefficients. It is important to determine the scaling coefficients for real geological formations to predict a real CO$_2$ sequestration process.

Developing a model to describe the shock-accelerated cylindrical fluid layer with arbitrary Atwood numbers is essential for uncovering the effect of Atwood numbers on the perturbation growth. The recent model (J. Fluid Mech., vol. 969, 2023, p. A6) reveals several contributions to the instability evolution of a shock-accelerated cylindrical fluid layer but its applicability is limited to cases with an absolute value of Atwood numbers close to $1$, due to the employment of the thin-shell correction and interface coupling effect of the fluid layer in vacuum. By employing the linear stability analysis on a cylindrical fluid layer in which two interfaces separate three arbitrary-density fluids, the present work generalizes the thin-shell correction and interface coupling effect, and thus, extends the recent model to cases with arbitrary Atwood numbers. The accuracy of this extended model in describing the instability evolution of the shock-accelerated fluid layer before reshock is confirmed via direct numerical simulations. In the verification simulations, three fluid-layer configurations are considered, where the outer and intermediate fluids remain fixed and the density of the inner fluid is reduced. Moreover, the mechanisms underlying the effect of the Atwood number at the inner interface on the perturbation growth are mainly elucidated by employing the model to analyse each contribution. As the Atwood number decreases, the dominant contribution of the Richtmyer–Meshkov instability is enhanced due to the stronger waves reverberated inside the layer, leading to weakened perturbation growth at initial in-phase interfaces and enhanced perturbation growth at initial anti-phase interfaces.

The resolvent analysis reveals the worst-case disturbances and the most amplified response in a fluid flow that can develop around a stationary base state. The recent work by Padovan et al. (J. Fluid Mech., vol. 900, 2020, A14) extended the classical resolvent analysis to the harmonic resolvent analysis framework by incorporating the time-varying nature of the base flow. The harmonic resolvent analysis can capture the triadic interactions between perturbations at two different frequencies through a base flow at a particular frequency. The singular values of the harmonic resolvent operator act as a gain between the spatiotemporal forcing and the response provided by the singular vectors. In the current study, we formulate the harmonic resolvent analysis framework for compressible flows based on the linearized Navier–Stokes equation (i.e. operator-based formulation). We validate our approach by applying the technique to the low-Mach-number flow past an airfoil. We further illustrate the application of this method to compressible cavity flows at Mach numbers of 0.6 and 0.8 with a length-to-depth ratio of $2$. For the cavity flow at a Mach number of 0.6, the harmonic resolvent analysis reveals that the nonlinear cross-frequency interactions dominate the amplification of perturbations at frequencies that are harmonics of the leading Rossiter mode in the nonlinear flow. The findings demonstrate a physically consistent representation of an energy transfer from slow-evolving modes toward fast-evolving modes in the flow through cross-frequency interactions. For the cavity flow at a Mach number of 0.8, the analysis also sheds light on the nature of cross-frequency interaction in a cavity flow with two coexisting resonances.

Modular floating solar farms exhibit periodic open surface coverages due to the strip configuration of floating modules that support the photovoltaic (PV) panels on top. The associated modulations in the surface boundary layer and its turbulence characteristics are investigated in the present study under fully developed open channel flows. Different coverage percentages of 100 % (i.e. full cover), 60 %, 30 % and 0 % (i.e. open surface) were tested and measurements were obtained using particle image velocimetry. The results showed that the turbulence statistics are similar when the coverage decreases from 100 % to 60 %. However, with 30 %, both the turbulence intensities and Reynolds stresses increase substantially, reaching up to 50 % higher compared with the 100 % coverage, and the boundary layer thickness increases by more than 25 %. The local skin friction beneath the openings increases by 50 %. Analysis of spanwise vortices and premultiplied spectra indicates that the periodic coverage elongates the hairpin vortex packets and reduces their inclination angle, imposing limitations on sustainable coherent structures. At 30 %, flow detachment and smaller-scale vortices become dominant, reducing the mean velocities and increasing the turbulence intensities. Decreasing coverage percentage with flow detachment also shifts the energy transfer to higher wavenumbers, increasing energy dissipation and decreasing the bulk flow velocity. The kinetic energy and Reynolds stress carried by very large-scale motions decreases from 40 %–50 % with the 100 % and 60 % coverage to around 30 %–40 % with the 30 % coverage. Further research studies involving spanwise heterogeneity, higher Reynolds number and varying submergence of PV modules are needed for environmental considerations.

We propose a theoretical method to decompose the solution of a Stokes flow past a body immersed in a confined fluid into two simpler problems, related separately to the two geometrical elements of these systems: (i) the body immersed in the unbounded fluid (represented by its Faxén operators); and (ii) the domain of the confinement (represented by its Stokesian multipoles). Specifically, by using a reflection method, and assuming linear and reciprocal boundary conditions (Procopio & Giona, Phys. Fluids, vol. 36, issue 3, 2024, 032016), we provide the expression for the velocity field, the forces, torques and higher-order moments acting on the body in terms of: (i) the volume moments of the body in the unbounded ambient flow; (ii) the multipoles in the domain of the confinement; (iii) the collection of all the volumetric moments on the body immersed in all the regular parts of the multipoles considered as ambient flows. A detailed convergence analysis of the reflection method is developed. In light of the practical applications, we estimate the truncation error committed by considering only the lower-order moments (thus, truncating the matrices) and the errors associated with the approximated expressions available in the literature for force and torques. We apply the theoretical results to the archetypal hydrodynamic system of a sphere with Navier-slip boundary conditions near a plane wall with no-slip boundary conditions, to determine forces and torques on a translating and rotating sphere as a function of the slip length and of the distance of the sphere from the plane. The hydromechanics of a spheroid is also addressed.

Non-helical turbulence within a linear shear flow has demonstrated efficient amplification of large-scale magnetic fields in numerical simulations, but its precise mechanism remains elusive. The incoherent $\alpha$ mechanism proposes that a zero-mean fluctuating transport coefficient $\alpha$ (linked to kinetic helicity) in the shear flow is a candidate driver. Previous renovating-flow models have proposed that the correlation time of helicity fluctuations must be sufficiently extended to overcome turbulent magnetic diffusivity, yet only empirical validation of this concept has been obtained. In this study, we conduct direct numerical simulations of weakly compressible non-helical hydrodynamic turbulence. We scrutinize the correlation times of velocity and kinetic helicity fluctuations in distinct flow configurations, including rotation, shearing and Keplerian flows, as well as the shearing burgulence counterpart. Our findings indicate that rotation contributes to a prolonged correlation time of helicity compared with velocity, particularly notable in auto-correlations of both volume-averaged quantities and individual Fourier modes due to the formation of large-scale vortices. In contrast, moderate shear strength does not exhibit significant scale separation, with shear flows elongating vortices in the shear direction. Shearing burgulence, characterized by shorter helicity correlation times, appears less conducive to hosting the incoherent $\alpha$ effect. Notably, at modest shear rates, only Keplerian flows exhibit sufficiently coherent helicity fluctuations, in contrast to shearing flows. However, the relative strength of helicity fluctuations compared with turbulent diffusivity is significantly lower, raising doubts about the viability of the incoherent $\alpha$ effect as a potential dynamo driver in the subsonic flows examined in this study.

In this work we study features of inertia-gravity wave turbulence in the rotating shallow water equations. On examining the dynamics of waves with varying rotation rates, we find that the turbulent cascade of waves is strongest at low rotation rates, forming a $k^{-2}$ energy spectrum, and a rich distribution of shocks in physical space. At high rotation rates, the forward cascade of waves weakens along with a steeper energy spectra and vanishing of shocks in physical space. The wave cascade is seen to be scale-local, resulting in a noticeable time interval for energy to get transferred from domain scale to dissipative scale. Furthermore, we find that the vortical flow has a non-negligible effect on the wave cascade, especially at high rotation rates. The vortical flow assists in the forward cascade of waves and shock formation at high rotation rates, while the waves by themselves in the absence of the vortical flow lack a forward cascade and shock formation at such high rotation rates. On investigating the physical space structures in the vortical flow and their connections to the wave cascade, we find that strain-dominant regions, that are located around the boundaries of coherent vortices, are the physical space regions that contribute majorly to the forward cascade of waves. Our results in general highlight intriguing features of dispersive inertia-gravity wave turbulence that are qualitatively similar to those seen in three-dimensional homogeneous isotropic turbulence and are beyond the predictions of asymptotic resonant wave interaction theory.

This paper presents advances towards the data-based control of periodic oscillator flows, from their fully developed regime to their equilibrium stabilized in closed loop, with linear time-invariant (LTI) controllers. The proposed approach directly builds upon the iterative method of Leclercq et al. (J. Fluid Mech., vol. 868, 2019, pp. 26–65) and provides several improvements for an efficient online implementation, aimed at being applicable in experiments. First, we use input–output data to construct an LTI mean transfer functions of the flow. The model is subsequently used for the design of an LTI controller with linear quadratic Gaussian synthesis, which is practical to automate online. Then, using the controller in a feedback loop, the flow shifts in phase space and oscillations are damped. The procedure is repeated until equilibrium is reached, by stacking controllers and performing balanced truncation to deal with the increasing order of the compound controller. In this article, we illustrate the method for the classic flow past a cylinder at Reynolds number $Re=100$. Care has been taken such that the method may be fully automated and hopefully used as a valuable tool in a forthcoming experiment.

The results of an experimental investigation of smooth-body adverse pressure gradient (APG) turbulent boundary layer flow separation and reattachment over a two-dimensional ramp are presented. These results are part of a larger archival smooth-body flow separation data set acquired in partnership with NASA Langley Research Center and archived on the NASA Turbulence Modeling Resource website. The experimental geometry provides initial canonical turbulent boundary layer growth under nominally zero pressure gradient conditions prior to encountering a smooth, two-dimensional, backward facing ramp geometry onto which a streamwise APG that is fully adjustable is imposed. Detailed surface and off-surface flow field measurements are used to fully characterize the smooth-body APG turbulent boundary layer separation and reattachment at multiple spanwise locations over the ramp geometry. Unsteady aspects of the flow separation are characterized. It is shown that the first and second spatial derivatives of the streamwise static surface pressure profile are sufficient to determine key detachment and reattachment locations. The imposed streamwise APG gives rise to inflectional mean velocity profiles and the associated formation of an embedded shear layer, which is shown to play a dominant role in the subsequent flow development. Similarity scaling is developed for both the mean velocity and turbulent stresses that is found to provide self-similar collapse of profiles for different regions of the ramp flow. Despite the highly non-equilibrium flow environment, a new similarity scaling proved capable of providing self-similar turbulent stress profiles over the full streamwise extent of flow separation and downstream reattachment.

The causal relevance of local flow conditions in open-channel turbulence is analysed using ensembles of interventional experiments in which the effect of perturbing the flow within a small cell is monitored at some future time. When this is done using the relative amplification of the perturbation energy, causality depends on the flow conditions within the cell before it is perturbed, and can be used as a probe of the flow dynamics. The key scaling parameter is the ambient shear, which is also the dominant diagnostic variable for wall-attached perturbations. Away from the wall, the relevant variables are the streamwise and wall-normal velocities. Causally significant cells are associated with sweeps that carry the perturbation towards the stronger shear near the wall, whereas irrelevant ones are associated with ejections that carry it towards the weaker shear in the outer layers. Causally significant and irrelevant cells are themselves organised into structures that share many characteristics with classical sweeps and ejections, such as forming spanwise pairs whose dimensions and geometry are similar to those of classical quadrants. At the wall, this is consistent with causally significant configurations in which a high-speed streak overtakes a low-speed one, and causally irrelevant ones in which the two streaks pull apart from each other. It is argued that this is probably associated with streak meandering.

A model is formulated of a two-dimensional migrating, or swimming, inviscid bubble in a viscous fluid whose unsteady displacement is caused by the spreading over its surface of an initial distribution of insoluble surfactant. Assuming small capillary and Reynolds numbers, and a linear equation of state giving the surface tension as a function of surfactant concentration, the quasi-steady Stokes flow around the bubble is found analytically and explicit formulas are determined for the time-dependent bubble speed and its final overall displacement. At infinite surface Péclet number this is done using a complex version of the method of characteristics to solve a complex partial differential equation of Burgers type. For a finite non-zero surface Péclet number, the problem is shown to be linearizable by a complex variant of the classical Cole–Hopf transformation. The formulation allows general statements to be made on the bubble speed and its total net displacement in terms of the initial surfactant distribution. A weak finite-time singularity in the surface activity associated with an isolated clean point on the bubble surface is also identified and studied in detail.

This study investigates the flow structures and combustion regimes in an axisymmetric cavity-based scramjet combustor with a total temperature of 1800 K and a high Reynolds number of approximately 1 × 107. The hydroxyl planar laser-induced fluorescence technique, along with the broadband flame emission and CH* chemiluminescence, is employed to visualize the instantaneous flame structure in the optically accessible cavity. The jet-wake flame stabilization mode is observed, with intense heat release occurring in the jet wake upstream of the cavity. A hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation approach is performed for the 0.18-equivalent-ratio case with a pressure-corrected flamelet/progress variable model. The combustion regime is identified mainly in the corrugated or wrinkled flamelet regime (approximately 102 < Da < 104, 103 < Ret < 105 where $Da$ is the Damköhler number and $Re_t$ is the turbulent Reynolds number). The combustion process is jointly dominated by supersonic combustion (which accounts for approximately 58 %) and subsonic combustion, although subsonic combustion has a higher heat release rate (peak value exceeding 1 × 109 J (m3s)−1). A partially premixed flame is observed, where the diffusion flame packages a considerable quantity of twisted premixed flame. The shockwave plays a critical role in generating vorticity by strengthening the volumetric expansion and baroclinic torque term, and it can facilitate the chemical reaction rates through the pressure and temperature surges, thereby enhancing the combustion. Combustion also shows a remarkable effect on the overall flow structures, and it drives alterations in the vorticity of the flow field. In turn, the turbulent flow facilitates the combustion and improves the flame stabilization by enhancing the reactant mixing and increasing the flame surface area.

During oscillatory wetting, a phase retardation emerges between contact angle variation and contact line velocity, presenting as a hysteresis loop in their correlation – an effect we term dynamic hysteresis. This phenomenon is found to be tunable by modifying the surface with different molecular layers. A comparative analysis of dynamic hysteresis, static hysteresis and contact line friction coefficients across diverse substrates reveals that dynamic hysteresis is not a result of dissipative effects but is instead proportionally linked to the static hysteresis of the surface. In the quest for appropriate conditions to model oscillatory contact line motion, we identify the generalized Hocking's linear law and modified generalized Navier boundary condition as alternative options for predicting realistic dynamic hysteresis.

Wind tunnel measurements of the incident turbulent velocity fields and axial forces on a horizontal axis turbine and porous disc analogues are reported. The models were tested in both a simulated atmospheric boundary layer (ABL) and in grid turbulence, allowing for a range of turbulence length scale to rotor diameter ratios to be considered. A theoretical framework to account for the combined effect of distortion and potential flow blocking in the induction zone is presented. In the case of very large length-scale turbulence to diameter ratios, where distortion effects are minimal, a quasi-steady approach is adopted for the effect of blocking. For the small length-scale ratio limit, the method is developed from the classical analyses for rapid distortion of turbulence and blockage from flow through a porous sheet of resistance. For general length-scale ratios, an efficient prediction method based on interpolation between the two length-scale ratio extremes is established. For very large length-scale ratios, a quasi-steady theory without distortion is appropriate for a rotor or disc in a simulated ABL. The small length-scale theory is applicable for tests conducted in grid turbulence. The results of the study can inform the prediction and interpretation of typical measurements of turbulence within the induction zone and the fluctuating loads on a rotor, at both prototype and full scale. This is of particular importance to fatigue load assessments.

Experiments of transitional shock wave–boundary layer interactions (SBLIs) over 6$^\circ$ and 10$^\circ$ compression ramps were performed at Mach number 1.65. The unit Reynolds number was varied by a factor of two between 5.6 million per metre and 11 million per metre. Schlieren flow visualization was performed, and mean flow measurements were made using Pitot probes. Free interaction theory was verified from pressure measurements for all operating conditions. A new non-dimensional parameter was developed for scaling the strength of the imposed shock, which was based on the pressure required to separate a boundary layer. The validity of this new scaling was supported by the reconciliation of large discrepancies in a diverse collection of experimental results on the length scales of transitional interactions. This non-dimensional scaling was also applied to turbulent interactions, where different models were used to determine the pressure required to separate a turbulent boundary layer. Finally, a direct comparison between transitional and turbulent SBLIs was made, which revealed new insights into the evolution of length scales based on the state of the boundary layer.

Linear instability analysis of a viscous swirling liquid jet surrounded by ambient gas is carried out by considering the significant influence of axial shear effect. The jet azimuthal flow is assumed as a Rankine vortex, and the non-uniform velocity distribution in the jet axial direction is approximated by parabolic and error functions. The enhancement of jet rotation is found to promote the predominant mode with larger azimuthal wavenumbers, and the mode transition is decided by the competition between centrifugal force and axial shear stress. Subsequently, the influence of the axial shear effect is examined through changing the degree of shear stress and the thickness of the gas velocity boundary layer. It is found that an increase of jet average velocity or surface velocity in the axial direction leads to the predominant mode transition to smaller azimuthal wavenumbers, due to the combined effects of shear stress and gas pressure perturbation. A larger velocity difference between ambient gas and liquid jet also promotes the predominant modes with smaller azimuthal wavenumbers, and the physical mechanism is attributed to gas pressure perturbation. Phase diagrams of different azimuthal modes are given and compared with the study of Kubitschek & Weidman (J. Fluid Mech., vol. 572, 2007, pp. 261–286), where a static swirling column without axial shear stress was considered. The strengthened axial shear effect is found to delay the transition of predominant modes with the increase of angular velocity. Experimental studies considering the swirling jets with different axial velocities are further carried out, which validate the theoretical findings. Different instability mechanisms and their transition rules are also identified through energy budget analysis. This study is expected to give scientific guidance on understanding the instability mechanisms of the swirling jets that widely exist in natural phenomena and engineering applications.

This paper explores active wake-flow control on a notchback Ahmed body using genetically inspired optimization. Hotwire and particle image velocimetry measurements record velocity data and flow structures in the wake. Pulsed jets at four actuation slots (two at the roof trailing edge, two at the side trailing edges) dynamically control the wake to minimize aerodynamic drag. The study achieves up to 9.2 % (without consideration of energy consumption) drag reduction, primarily by manipulating vortices from the roof rear end. The paper elucidates the underlying flow mechanism and evaluates various actuation strategies, highlighting how optimal control leads to reattachment of wake separation at the rear slant, diminishing the slant bubble and promoting downstream reattachment for enhanced drag reduction.