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.

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.

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 review some of the processes leading to dispersion and mixing in porous media, exploring the differences between the travel time distribution of fluid particles within a pore throat and between pore throats of different size within the porous layer. A recent paper of Liu et al. has combined a model of these travel time distributions with a continuous time random walk to quantify the dispersion as a function of the Peclet number. We describe some further problems relating to dispersive mixing of tracer which may be amenable to this approach, including dispersion caused by macroscopic lenses of different permeability, dispersion of tracer which partitions between the fluid and matrix and the effects of buoyancy on mixing.

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.

This study investigates the spatial distribution of inertial particles in turbulent Taylor–Couette flow. Direct numerical simulations are performed using a one-way coupled Eulerian–Lagrangian approach, with a fixed inner-wall Reynolds number of 2500 for the carrier flow, while the particle Stokes number ($St$) varies from 0.034 to 1 for the dispersed phase. We first examine the issue of preferential concentration of particles near the outer-wall region. Employing two-dimensional Voronoï analysis, we observe a pronounced particle clustering with increasing $St$, particularly evident in regions of low fluid velocity. Additionally, we investigate the concentration balance equation, inspired by the work of Johnson et al. (J. Fluid Mech., vol. 883, 2020, A27), to examine the particle radial distribution. We discern the predominant sources of influence, namely biased sampling, turbophoresis and centrifugal effects. Across all cases, centrifugal force emerges as the primary driver, causing particle migration toward the outer wall. Biased sampling predominantly affects smaller inertial particles, driving them toward the inner wall due to sampling within Taylor rolls with inward radial velocity. Conversely, turbophoresis primarily impacts larger inertial particles, inducing migration towards both walls where turbulent intensity is weaker compared with the bulk. With the revealed physics, our work provides a basis for predicting and controlling particle movement and distribution in industrial applications.

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.

The flow-induced oscillation of a transversely clamped buckled flexible filament in a uniform flow was explored using the penalty immersed boundary method. Both inverted and conventional configurations were analysed. The effects of bending rigidity, filament length and Reynolds number were examined. As these parameters were varied, four distinct modes were identified: conventional transverse oscillation mode, deflected oscillation mode, inverted transverse oscillation mode and structurally steady mode. The filament exhibited a 2S wake pattern under the conventional transverse oscillation mode and the small-amplitude inverted transverse oscillation mode, a P wake pattern under the deflected oscillation mode and a 2S + 2P wake pattern for the large-amplitude inverted transverse oscillation mode. Irrespective of their initial conditions, all of the filaments converged to the conventional transverse oscillation mode under low bending rigidity. Multistability was observed in the transversely clamped buckled flexible filament under moderate bending rigidity. The deflection in the oscillation mode increased with increasing filament length. The inverted buckled filament was sensitive to the Reynolds number, unlike the conventional buckled filament. The transverse oscillation mode demonstrated superior energy-harvesting performance.

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.

Acoustic resonances in cascade structures may cause structural damage and instability problems in aero-engines and other industrial plants; thus, developing corresponding prediction methods is important. However, works published in the open literature mostly focus on the special case of the stationary Parker modes and provide little knowledge into the rotating resonances in annular cascades, especially in the presence of non-zero background mean flows. This paper develops a three-dimensional semi-analytic model to study the acoustic resonances in an annular cascade in the presence of axial mean flow. The model applies an unsteady cascade response based on the three-dimensional lifting surface method to construct a matrix equation. Characteristic frequencies are solved in the complex domain by numerically searching for singular points. Both the oscillation frequency and the growth rate of the three-dimensional resonance modes are theoretically calculated for the first time under non-zero mean flow conditions. The results reveal an organised distribution with varying inter-blade phase angle and show obvious change with the background flow speed. It is found that the unsteady vortex shedding from the trailing edges of the cascade is a key factor influencing the dissipation rate of the resonance modes. In addition, the important effects of acoustic scattering by the cascade during resonances are examined, which qualitatively corroborate some previous experimental observations.

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.

A numerical study supplemented with theoretical analysis is made, to analyse the electrophoresis of highly charged soft particles in electrolytes with trivalent counterions. The electrokinetic model is devised under the continuum hypothesis, which incorporates the ion–ion electrostatic correlations, hydrodynamic steric interactions of finite sized ions and ion–solvent interactions. The governing equations for ion transport and electric field are derived from the volumetric free energy of the system, which includes the first-order correction for the non-local electrostatic correlations of interacting ions, excess electrochemical potential due to finite ion size as well as the Born energy difference of ions due to dielectric permittivity variation. The electrolyte viscosity is considered to be a function of the local volume fraction of finite-sized ions, which causes the diffusivity of ions to vary locally. The occurrence of mobility reversal of a soft particle having the same polarity of its core and soft shell charge and formation of a coion-dominated zone in the soft layer is elaborated through this study. This can explain the mechanisms for the attraction between like-charged soft particles, as seen in the condensation of DNAs. The impact of ion–ion correlations and ion–solvent interactions of finite-sized ions are analysed by comparing them with the results based on the standard model. At a higher range of the core charge density, the ion–ion correlations induce a condensed layer of counterions on the outer surface of the core, which draws coions in the electric double layer, leading to an inversion in polarity of the charge density and mobility reversal. The dielectric decrement and ion steric interactions create a saturation in ion distribution and hence, modify the condensed layer of counterions. The enhanced fixed charge density of the polyelectrolyte layer diminishes the ion correlations due to the stronger screening effects and prevents the formation of a coion dominated zone in the Debye layer. The impact of the counterion size and the mixture of monovalent and trivalent counterions on mobility is analysed.

We analysed the performance of convolutional autoencoders in generating reduced-order representations of the temperature field of two-dimensional Rayleigh–Bénard flows at $\textit{Pr} =1$ and Rayleigh numbers extending from $10^6$ to $10^8$, capturing the range where the flow transitions to turbulence. We present a way of estimating the minimum number of dimensions needed by the autoencoders to capture all the relevant physical scales of the data that is more apt for highly multiscale flows than previous criteria applied to lower-dimensional systems. We compare our architecture with two regularized variants as well as with linear methods, and find that manually fixing the dimension of the latent space produces the best results. We show how the estimated minimum dimension presents a sharp increase around $Ra\sim 10^7$, when the flow starts to transition to turbulence. Furthermore, we show how this dimension does not follow the same scaling as the physically relevant scales, such as the dissipation length scale and the thermal boundary layer.

We studied flow organization and heat transfer properties in mixed turbulent convection within Poiseuille–Rayleigh–Bénard channels subjected to temporally modulated sinusoidal wall temperatures. Three-dimensional direct numerical simulations were performed for Rayleigh numbers in the range $10^6 \leqslant Ra \leqslant 10^8$, a Prandtl number $Pr = 0.71$ and a bulk Reynolds number $Re_b \approx 5623$. We found that high-frequency wall temperature oscillations had minimal impact on flow structures, while low-frequency oscillations induced adaptive changes, forming stable stratified layers during cooling. Proper orthogonal decomposition (POD) analysis revealed a dominant streamwise unidirectional shear flow mode. Large-scale rolls oriented in the streamwise direction appeared as higher POD modes and were significantly influenced by lower-frequency wall temperature variations. Long-time-averaged statistics showed that the Nusselt number increased with decreasing frequency by up to 96 %, while the friction coefficient varied by less than 15 %. High-frequency modulation predominantly influenced near-wall regions, enhancing convective effects, whereas low frequencies reduced these effects via stable stratified layer formation. Phase-averaged statistics showed that high-frequency modulation resulted in phase-stable streamwise velocity and temperature profiles, while low-frequency modulation caused significant variations due to weakened turbulence. Turbulent kinetic energy (TKE) profiles remained high near the wall during both heating and cooling at high frequency, but decreased during cooling at low frequencies. A TKE budget analysis revealed that during heating, TKE production was dominated by shear near the wall and by buoyancy in the bulk region; while during cooling, the production, distribution and dissipation of TKE were all nearly zero.

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.

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 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.

This study investigates how the spatial configuration of submerged three-dimensional patches of vegetation impacts turbulence. Laboratory experiments were conducted in a channel with submerged patches of model vegetation configured with different patch area densities ($\phi _{p}$), representing the bed area fraction occupied by patches, ranging from 0.13 to 0.78, and different spatial patterns transitioning from two dimensional (channel-spanning patches) to three dimensional (laterally unconfined patches). These configurations produced a range of flow regimes within the canopy, from wake interference flow to skimming flow. At low area density ($\phi _{p}\lt0.5$), turbulence within the canopy increased with increasing $\phi _{p}$ regardless of spatial configuration, while at high area density ($\phi _{p}\gt0.5$), the relationship between turbulence and $\phi _{p}$ depended on the spatial configuration of the patches. For the same patch area density, the configuration with smaller lateral gaps generated stronger turbulence within the canopy. The relative contributions of wake and shear production also varied with the spatial configuration of the patches. At low area densities, wake production dominated over shear production, while at high area densities, shear production was more dominant due to an enhanced shear layer at the top of the canopy and reduced mean velocity within the canopy. A new predictive model for the channel-averaged turbulent kinetic energy (TKE) was developed as a function of channel-averaged velocity, canopy geometry, and patch area density, which showed good agreement with the measured TKE.

We find the optimally time-dependent (OTD) orthogonal modes about a time-varying flow generated by a strong gust vortex impacting a NACA 0012 airfoil. This OTD analysis reveals the amplification characteristics of perturbations about the unsteady base flow and their amplified spatiotemporal structures that evolve over time. We consider four time-varying laminar base flows in which a vortex with a strength corresponding to the gust ratio $G$ of $\{-1,-0.5,0.5,1\}$ impinges on the leading edge of the airfoil at an angle of attack of $12^\circ$. In these cases, the impingement of the strong gust vortex causes massive separation and the generation of large-scale vortices around the airfoil within two convective time units. As these flow structures develop around the airfoil on a short time scale, the airfoil experiences large transient vortical lift variations in the positive and negative directions that are approximately five to ten times larger than the baseline lift. The highly unsteady nature of these vortex–airfoil interactions necessitates an advanced analytical technique capable of capturing the transient perturbation dynamics. For each of the considered gust ratios, the OTD analysis identifies the most amplified region to perturbations, the location of which changes as the wake evolves differently. For interactions between a moderate positive vortex gust ($G=0.5$) and the airfoil, the area where perturbations are amplified transitions from the leading-edge vortex (LEV) sheet to the forming LEV. Later, this most amplified structure becomes supported in the airfoil wake directly behind the trailing edge. In contrast, a strong vortex gust ($G=\pm 1$) encountered by the airfoil shows the most amplified OTD mode to appear around the core of the shed vortices. This study provides an analysis technique and fundamental insights into the broader family of unsteady aerodynamic problems.