The impact of intrinsic compressibility effects – changes in fluid volume due to pressure variations – on high-speed wall-bounded turbulence has often been overlooked or incorrectly attributed to mean property variations. To quantify these intrinsic compressibility effects unambiguously, we perform direct numerical simulations of compressible turbulent channel flows with nearly uniform mean properties. Our simulations reveal that intrinsic compressibility effects yield a significant upward shift in the logarithmic mean velocity profile that can be attributed to the reduction in the turbulent shear stress. This reduction stems from the weakening of the near-wall quasi-streamwise vortices. In turn, we attribute this weakening to the spontaneous opposition of sweeps and ejections from the near-wall expansions and contractions of the fluid, and provide a theoretical explanation for this mechanism. Our results also demonstrate that intrinsic compressibility effects play a crucial role in the increase in inner-scaled streamwise turbulence intensity in compressible flows, as compared with incompressible flows, which was previously regarded to be an effect of mean property variations alone.

Coherent structures over two distinct, organized wall perturbations – a transverse sinusoidal bump with and without small-scale longitudinal grooves – are studied using direct numerical simulations. Large-scale spanwise rollers (SRs) form via shear layer rollup past the bump peak, enveloping a large separation bubble (SB) for both a smooth wall (SW) and a grooved wall (GW). In a GW, small-scale alternatingly spinning jets emanating from the crests’ corners merge with the shear layer, altering the SRs compared with SRs in a SW. The underlying coherence of the highly turbulent SRs is educed via phase-locked ensemble averaging. Coherent vorticity contours of SRs are ellipses tilted downward, hence causing co-gradient Reynolds stress. The limited streamwise length of SB precludes SR tumbling, unlike in a free shear layer. The coherent field reveals minibubbles attached to the bump’s downstream wall with circulation opposite to that of the SB – they are larger, stronger and more numerous in GW than in SW – reducing skin friction. Compared with SW, the swirling jets in GW increase coherent production while decreasing incoherent production. Additionally, the jets push the SRs to travel faster and farther before reattachment. The SB experiences two different modes of oscillation due to high-frequency advection of the shear layer SR and low-frequency breathing of the SB, where the former dominates in GW and the latter in SW. Negative production is caused by counter-rotating vortex dipoles inducing flow ejections (for both SW and GW) and single vortices penetrating the grooves – both occurring in the region of flow acceleration.

We study the effect of turbulence on collisions between a finite-size bubble and small inertial particles based on interface-resolved simulations. Our results show that the interaction with the flow field around the bubble remains the dominant effect. Nonlinear dependencies in this process can enhance the turbulent collision rate by up to 100 % compared to quiescent flow. Fluctuations in the bubble slip velocity during the interaction with the particle additionally increase the collision rate. We present a frozen-turbulence model that captures the relevant effects providing a physically consistent framework to model collisions of small inertial particles with finite-sized objects in turbulence.

Deep reinforcement learning (DRL) is employed to develop control strategies for drag reduction in direct numerical simulations of turbulent channel flows at high Reynolds numbers. The DRL agent uses near-wall streamwise velocity fluctuations as input to modulate wall blowing and suction velocities. These DRL-based strategies achieve significant drag reduction, with maximum rates $35.6\,\%$ at $Re_{\tau }\thickapprox 180$, $30.4\,\%$ at $Re_{\tau }\thickapprox 550$, and $27.7\,\%$ at $Re_{\tau }\thickapprox 1000$, outperforming traditional opposition control methods. An expanded range of wall actions further enhances drag reduction, although effectiveness decreases at higher Reynolds numbers. The DRL models elevate the virtual wall through blowing and suction, aiding in drag reduction. However, at higher Reynolds numbers, the amplitude modulation of large-scale structures significantly increases the residual Reynolds stress on the virtual wall, diminishing the drag reduction. Analysis of budget equations provides a systematic understanding of the underlying drag reduction dynamics. The DRL models reduce skin friction by inhibiting the redistribution of wall-normal turbulent kinetic energy. This further suppresses the wall-normal velocity fluctuations, reducing the production of Reynolds stress, thereby decreasing skin friction. This study showcases the successful application of DRL in turbulence control at high Reynolds numbers, and elucidates the nonlinear control mechanisms underlying the observed drag reduction.

The resonance mechanism in the initial of wind-wave generation proposed by Phillips (1957. J. Fluid Mech. 2, 417–445) is a foundation of wind-wave generation theory, but a precise theoretical quantification of wave energy growth in this initial stage has not been obtained yet after more than six decades of research. In this study, we aim to address this knowledge gap by developing an analytical approach based on a novel complex analysis method to theoretically investigate the temporal evolution of the wave energy in the Phillips initial stage. We quantitatively derive and analyse the growth behaviour of the surface wave energy and obtain an analytical solution for its upper bound. Our result highlights the crucial effects of surface tension. Because the phase velocity of gravity–capillary waves has a minimal value at a critical wavenumber, gravity–capillary waves and gravity waves (which neglect surface tension) exhibit distinct resonance curve properties and wave energy growth behaviours. For gravity waves, the resonance curve extends indefinitely; for gravity–capillary waves, it either forms a finite-length curve or does not exist, depending on the wind speed. The leading-order term of the upper-bound solution of the energy of gravity waves increases linearly over time, while for gravity–capillary waves, the term increases linearly over time under strong wind conditions but remains finite under weak wind conditions. This theoretical study provides an analytical framework for the generation of wind-waves in the Phillips initial stage, which may inspire further theoretical, numerical and experimental research.

This study introduces a novel approach to investigate the Reynolds analogy in complex flow scenarios. It is shown that the total mechanical energy $\mathit {B}$, viz. the sum of kinetic energy and pressure work, and the field $\Gamma =\theta ^2/2$ (where $\theta$ is the transported passive scalar) are governed by two equations that are similar in form, when time-averaged for statistically stationary flows. For fully developed channel flows the integral energy balance links the mean bulk velocity and scalar with the volume averages of the respective dissipation rates, allowing the assessment of the Reynolds analogy in terms of the dissipation fields. This approach is tested on direct numerical simulation data of rough-wall turbulent channel flow at two different roughness Reynolds numbers, namely $k^+=15$ and $k^+=90$. For a unit Prandtl number, the same qualitative behaviour is observed for the mean wall-normal distributions of the budget-equation terms of $B$ and $\Gamma$, the latter being larger than the corresponding terms in the mechanical-energy budget. The Reynolds decomposition of the flow into temporal mean and stochastic parts reveals that roughness primarily affects the mean-flow dissipation. For the $k^+=90$ case, the analysis shows that attached-flow and high-shear regions dominate the integral mean scalar and momentum transfer and exhibit the greatest differences between the mean mechanical and scalar dissipation rates. In contrast, well-mixed regions, sheltered by large roughness elements, contribute similarly and minimally to the integral scalar and momentum transfer.

We investigate the linear instability of flows that are stable according to Rayleigh’s criterion for rotating fluids. Using Taylor–Couette flow as a primary test case, we develop large-Reynolds-number-matched asymptotic expansion theories. Our theoretical results not only aid in detecting instabilities previously reported by Deguchi (Phys. Rev. E, vol 95, 2017, p. 021102(R)) across a wide parameter range, but also clarify the physical mechanisms behind this counterintuitive phenomenon. Instability arises from the interaction between large-scale inviscid vortices and the viscous flow structure near the wall, which is analogous to Tollmien–Schlichting waves. Furthermore, our asymptotic theories and numerical computations reveal that similar instability mechanisms occur in boundary layer flows over convex walls.

The turbulent boundary layer is a region where both preferential dissipation of energy and the production of significant vorticity arises as a consequence of the strong velocity gradients. Previous work has shown that, following a Reynolds decomposition, the purely fluctuating component of the enstrophy production is the dominant term. Near the wall this varies in a complex manner with height. In this study, we additionally decompose the strain rate and vorticity terms into normal and non-normal components using a Schur decomposition and are able to explain all these features in terms of contributions at different heights from constituents involving different combinations of normal and non-normal quantities. What is surprising about our results is that, while the mean shear and the action of larger-scale structures should mean that non-normal effects are of over-riding importance at the wall, the most important individual term involves the fluctuating normal strain rate in the transverse direction. In part, this is because of a strong correlation between this term and the non-normal vorticity with a transverse axis, but it is also the case that individual components of the purely non-normal enstrophy production are negative in the mean. Hence, a local strain rate that is orthogonal to the direction of the dominant mean and fluctuating shear plays a crucial role in amplifying vorticity that is yet to have developed a local component. These conclusions support the emphasis in the control literature on the transverse velocity components at the wall.

Understanding wave kinematics is crucial for analysing the thermodynamic effects of sloshing, which can lead to pressure drops in non-isothermal cryogenic fuel tanks. In the research reported here, Faraday waves in a horizontal circular tank (partially filled with water) under vertical excitation are investigated. The tank geometry is referred to as a horizontal circular tank throughout, with its circular face oriented perpendicular to the horizontal plane. Firstly, this paper addresses the eigenvalue problem through linear potential flow theory, in order to provide theoretical evidence of Faraday waves in horizontal circular tanks, the impact the density ratio has on the eigenvalues is then considered. Secondly, an experimental investigation testing multiple liquid fill levels is conducted. A soft-spring nonlinear response is demonstrated throughout the parameter space. The results showed larger sloshing amplitudes for low fill levels and smaller sloshing amplitudes for high fill levels. Asymmetry between anti-nodes at the container sidewalls and through the tank centreline are evident for low fill levels. Moreover, the sloshing wave amplitude at which breaking waves occur is smaller for high fill level conditions. Finally, period tripling was observed for all fill levels tested, confirming nonlinear mode interactions before the onset to wave breaking.

Inspired by laboratory experiments showing internal waves generated by a plume impinging upon a stratified fluid layer (Ansong & Sutherland. 2010 J. Fluid Mech. 648, 405–434), we perform large eddy simulations in three dimensions to examine the structure and source of internal waves emanating from the top of a plume that rises vertically into stratification whose strength ranges over two orders of magnitude between different simulations. Provided the plume is sufficiently energetic to penetrate into the stratified layer, internal waves are generated with frequencies in a relatively narrow band moderately smaller than the buoyancy frequency. Through adaptations of ray theory including viscosity and use of dynamic mode decomposition, we show that the waves originate from within the turbulent flow rather than at the turbulent/non-turbulent interface between the fountain top and the surrounding stratified fluid.

We report the first shock-tube experiments on Richtmyer–Meshkov instability at a single-mode light–heavy interface accelerated by a strong shock wave with Mach number higher than 3.0. Under the proximity effect of the transmitted shock and its induced secondary compression effect, the interface profile is markedly different from that in weakly compressible flows. For the first time, the validity of the compressible linear theory and the failure of the impulsive model in predicting the linear amplitude evolution in highly compressible flows are verified through experiments. Existing nonlinear and modal models fail to accurately describe the perturbation evolution, as they do not account for the shock proximity and secondary compression effects on interface evolution. The shock proximity effect manifests mainly in the early stages when the transmitted shock remains close to the interface, while the effect of secondary compression manifests primarily at the period when interactions of transverse shocks occur at the bubble tips. Based on these findings, we propose an empirical model capable of predicting the bubble evolution in highly compressible flows.

Evolution of solitary waves and an undular bore intruding through an abrupt transition from a wide basin into a narrow channel with opposing current is investigated. The laboratory experiments are performed in a wave tank that is crafted to achieve a steady and symmetrical shallow-water jet in the basin. The channel has a breadth comparable to the wave lengths, and the flow has Froude number approximately 0.1. The opposing current amplifies and slows the incoming waves on the jet in the basin, but the propagation speed is faster than the local Doppler effect of the current due to the influence of the wave propagating in the flank of the jet. At the channel mouth, the wave amplitude is enhanced due to the waveform altered by the current in the basin, although the amplification in the upstream channel is similar with and without the current. The longer incident waves have greater amplification into the channel. The leading wave of the undular bore is impacted by the opposing flow and transition similarly to the solitary waves. In contrast, the subsequent waves of the undular bore have a complex phase interference on the jet that causes disconnection in the lateral wave formation across the breadth of the jet. At the transition, the subsequent waves exhibit greater amplification than the leading one due to accumulated wave energy at the channel mouth. The intrusion of the undular bore against the current further enhances a rise in mean water level in the channel.

Ocean turbulence at meso- and submesocales affects the propagation of surface waves through refraction and scattering, inducing spatial modulations in significant wave height (SWH). We develop a theoretical framework that relates these modulations to the current that induces them. We exploit the asymptotic smallness of the ratio of typical current speed to wave group speed to derive a linear map – the U2H map – between surface current velocity and SWH anomaly. The U2H map is a convolution, non-local in space, expressible as a product in Fourier space by a factor independent of the magnitude of the wavenumber vector. Analytic expressions of the U2H map show how the SWH responds differently to the vortical and divergent parts of the current, and how the anisotropy of the wave spectrum is key to large current-induced SWH anomalies. We implement the U2H map numerically and test its predictions against WAVEWATCH III numerical simulations for both idealised and realistic current configurations.

In this experimental and numerical study, we revisit the question of the onset of the elastic regime in viscoelastic pinch-off. This is relevant to all modern filament thinning techniques, which aim to measure the extensional properties of low-viscosity polymer solutions. Examples are the slow retraction method (SRM) for capillary breakup extensional rheometry (CaBER), or the dripping method, in which a drop detaches from a nozzle. As part of these techniques, a stable liquid bridge is brought slowly to its stability threshold, where capillary-driven thinning starts. This thinning slows down dramatically at a critical radius $h_1$, marking the onset of the elasto-capillary regime, characterised by a filament of nearly uniform radius. While a theoretical scaling exists for this transition in the case of the classical step-strain CaBER protocol, where polymer chains stretch without relaxing during the fast plate separation, we show that this theory is not necessarily valid for a slow protocol such as the SRM. In that case, polymer chains start stretching (beyond their equilibrium coiled configuration) only when the bridge thinning rate becomes comparable to the inverse of their relaxation time. We derive a universal scaling for $h_1$, valid for both low- and high-viscosity polymer solutions. This scaling is validated by CaBER experiments with a slow plate separation protocol using different polymer solutions, plate diameters and sample volumes, as well as by numerical simulations using the FENE-P model.

Sea ice is a mushy layer, a porous material whose properties depend on the relative proportions of solid and liquid. The growth of sea ice is governed by heat transfer through the ice together with appropriate boundary conditions at the interfaces with the atmosphere and ocean. The salinity of sea ice has a major effect on its thermal properties so might naïvely be expected to have a major effect on its growth rate. However, previous studies observed a low sensitivity throughout the winter growth season. The goal of this study is to identify the controlling physical mechanisms that explain this observation. We develop a simplified quasi-static framework by applying a similarity transformation to the underlying heat equation and neglecting the explicit time dependence. We find three key processes controlling the sensitivity of growth rate to salinity. First, the trade-off between thermal conductivity and (latent) heat capacity leads to low sensitivity to salinity even at moderately high salinity and brine volume fraction. Second, the feedback on the temperature profile reduces the sensitivity relative to models that assume a linear profile, such as zero-layer Semtner models. Third, thicker ice has the opposite sensitivity of growth rate to salinity compared with thinner ice, sensitivities that counteract each other as the ice grows. Beyond its use in diagnosing these sensitivities, we show that the quasi-static approach offers a valuable sea-ice model of intermediate complexity between zero-layer Semtner models and full partial-differential-equation-based models such as Maykut–Untersteiner/Bitz–Lipscomb and mushy-layer models.

In this work, the stability and transition to turbulence over blunt flat plates with different leading-edge radii are investigated computationally. The benchmark experimental work for comparative studies is conducted by Borovoy et al. (AIAA J., vol. 60, 2022, pp. 497–507). The freestream Mach number is 5, the unit Reynolds number is $6\times 10^7$ m$^{-1}$, and the maximum nose-tip radius 3 mm exceeds the experimental reversal value. High-resolution numerical simulation and stability analysis are performed. Three-dimensional broadband perturbation is added on the far field boundary to initiate the transition. The highlight of this work is that the complete physical process is considered, including the three-dimensional receptivity, linear and nonlinear instabilities, and transition. The experimental reversal phenomenon is reproduced favourably in the numerical simulation for the first time. Linear stability analysis shows that unstable first and second modes are absent in the blunt-plate flows owing to the presence of the entropy layer, although these modes are evident in the sharp-leading-edge case. Therefore, the transition on the blunt plate is due to non-modal instabilities. Numerical results for all the blunt-plate cases reveal the formation of streamwise streaky structures downstream of the nose (stage I) and then the presence of intermittent turbulent spots in the transitional region (stage II). In stage I, a preferential spanwise wavelength approximately 0.9 mm is selected for all the nose-tip radii, and low-frequency components are dominant. In stage II, high-frequency secondary instabilities appear to grow, which participate in the eventual breakdown. By contrast, leading-edge streaks are not remarkable in the sharp-leading-edge case, where transition is induced by oblique first and Mack second modes. The transition reversal beyond the critical nose-tip radius arises from an increasing magnitude of the streaky response in the early stage, while the transition mechanism stays similar qualitatively.

Understanding the mechanism of hydrodynamic cloud cavitation is crucial to reducing noise, vibration and wear. Recent studies have clarified the physics of two distinct formation mechanisms of cloud cavitation. Ganesh et al. (J. Fluid Mech., vol. 802, 2016, pp. 37–78) identified the propagation of bubbly shockwaves as a cloud detachment mechanism. Pelz et al. (J. Fluid Mech., vol. 817, 2017, pp. 439–454) explained the influence of Reynolds number and cavitation number on asymptotic growth of the cavity sheet and its periodic shedding caused by re-entrant flow. In this paper the two mechanisms are set in relation to each other. For this, we show firstly that the transition from re-entrant flow to shockwave-driven cloud cavitation is given by a kinematic condition, namely the asymptotic sheet length equates to the chord length, $\hat {a}=L$. For $\hat {a}>L$ shockwave-driven cloud cavitation dominates. For $\hat {a}< L$ re-entrant flow-driven cloud cavitation dominates. As the cavitation number decreases, the closure region of the cavity sheet reaches the trailing edge of the hydrofoil, identifying the trailing edge as a trigger for condensation shockwaves, particularly as re-entrant flow-driven cavitation diminishes. Since the sheet length is an implicit function of the cavitation number, the kinematic condition $\hat {a}/L=1$ results in a critical cavitation number ${\sigma _\mathrm {II,III}}$ that is calculated analytically and validated by experiments. Secondly, we derive the relationship between the Strouhal number and the asymptotic sheet length for re-entrant flow-driven cloud cavitation. The model presented here is thoroughly validated by experiments.

We demonstrate the first successful non-invasive stabilisation of nonlinear travelling waves in a straight cylindrical pipe using time-delayed feedback control working in various symmetric subspaces. By using an approximate linear stability analysis and by analysing the frequency-domain effect of the control using transfer functions, we find that solutions with well-separated unstable eigenfrequencies can have narrow windows of stabilising time delays. To mitigate this issue we employ a ‘multiple time-delayed feedback’ approach, where several control terms are included to attenuate a broad range of unstable eigenfrequencies. We implement a gradient descent method to dynamically adjust the gain functions in order to reduce the need for tuning a high-dimensional parameter space. This results in a novel control method where the properties of the target state are not needed in advance, and speculative guesses can result in robust stabilisation. This enables travelling waves to be stabilised from generic turbulent states and unknown travelling waves to be obtained in highly symmetric subspaces.

The instantaneous structure of a turbulent boundary layer (TBL) subjected to freestream turbulence (FST) is investigated at several streamwise locations downstream of an active turbulence-generating grid. Using planar particle image velocimetry, three grid sequences are tested at four streamwise locations with FST intensities up to 10.9 %. A low-turbulence reference case is included for comparison. A novel method is proposed to separate the instantaneous TBL and FST flows by identifying a distinct interface for each realisation using probability density functions of the vorticity field. Two alternative approaches are used to define the interfaces, based on either constant velocity contour lines or constant vorticity magnitude contour lines. The former is found to highlight the momentum events in the velocity fields, whereas the latter outlines the vortical features of the flow. Regardless of the interface choice, when faced with FST, the interface moves closer to the wall on average, and its location fluctuates more. When FST is present, the shear and mean spanwise vorticity magnitudes increase on the TBL side of the interface. Uniform momentum zones (UMZs) beneath the velocity interfaces are identified. In the presence of FST, UMZs located closer to the wall appear to be compressed, resulting in fewer identified UMZs. Moving downstream, the FST intensity decays while the TBL develops. As a result, many characteristics of the TBL recover to an undisturbed state, with the interface moving away from the wall, vorticity and turbulent fluctuations returning to their natural state undisturbed by FST and the number of detected UMZs increasing.