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.

Wall modelling in large-eddy simulation (LES) is of high importance to allow scale resolving simulations of industrial applications. Numerous models were developed and validated for incompressible flows, including a simple quasi-analytical model based on Reichardt's formula that approximates the law of the wall. In this paper, a scaling is proposed to generalize this wall model to highly compressible flows. First, the results of wall-resolved LES (wrLES) of adiabatic compressible turbulent channel flows at $Re_\tau = 1000$ and at centreline Mach numbers of $M_c= 0.76$ and $1.5$ are presented. Then, three potential scalings of the incompressible wall model are proposed, and their a priori performance is evaluated : (i) the Howarth–Stewartson scaling, (ii) an improved Van Driest scaling and (iii) a new scaling obtained from a blending of those two. The results of wall-modelled LES (wmLES) of compressible channel flows using these three models are compared with the reference wrLES data, showing the superior accuracy of the hybrid scaling. The consistency of the new wall model at low Mach numbers is also verified by comparing the results of a wmLES at $M_c= 0.25$ with those of reference incompressible DNS data at $Re_\tau = 1000$ and $5200$. Finally, the proposed wall model is also applied to a turbulent channel flow at $M_c=1.5$ and $Re_\tau =5200$.

First predicted by Richtmyer in 1960 and experimentally confirmed by Meshkov in 1969, the Richtmyer–Meshkov instability (RMI) is crucial in fields such as physics, astrophysics, inertial confinement fusion and high-energy-density physics. These disciplines often deal with strong shocks moving through condensed materials or high-pressure plasmas that exhibit non-ideal equations of state (EoS), thus requiring theoretical models with realistic fluid EoS for accurate RMI simulations. Approximate formulae for asymptotic growth rates, like those proposed by Richtmyer, are helpful but rely on heuristic prescriptions for compressible materials. These prescriptions can sometimes approximate the RMI growth rate well, but their accuracy remains uncertain without exact solutions, as the fully compressible RMI growth rate is influenced by both vorticity deposited during shock refraction and multiple sonic wave refractions. This study advances previous work by presenting an analytic, fully compressible theory of RMI for reflected shocks with arbitrary EoS. It compares theoretical predictions with heuristic prescriptions using ideal gas, van der Waals gas and three-term constitutive equations for simple metals, the latter being analysed with detailed and simplified ideal-gas-like EoS. We additionally offer an alternative explicit approximate formula for the asymptotic growth rate. The comprehensive model also incorporates the effects of constant-amplitude acoustic waves at the interface, associated with the D'yakov–Kontorovich instability in shocks.

The composition of the Earth's mantle, built by impacts of metallic core planetesimals, depends on the time offered to metal–silicate equilibration (concentration equalization by mixing) as the impactor fragments settle down. Smaller fragments equilibrate while bigger ones remain segregated and accumulate in the – therefore iron rich – Earth's core. To understand the primary instability setting the interpenetration depth and wavelength between the impactor and the mantle, Lherm et al. (J. Fluid Mech., vol. 937, 2022, A20) analyse the intermingled phenomena occurring when a drop impacts a substrate of a pool of another liquid, digging a crater. The density difference between the phases triggers a Rayleigh–Taylor instability when the drop is denser than the pool. The overall planet chemical composition relies on a tradeoff between fragmentation and mixing; this study provides elements to decipher their relative importance.

We present a deep probabilistic convolutional neural network (PCNN) model for predicting local values of small-scale mixing properties in stratified turbulent flows, namely the dissipation rates of turbulent kinetic energy and density variance, $\varepsilon$ and $\chi$. Inputs to the PCNN are vertical columns of velocity and density gradients, motivated by data typically available from microstructure profilers in the ocean. The architecture is designed to enable the model to capture several characteristic features of stratified turbulence, in particular the dependence of small-scale isotropy on the buoyancy Reynolds number $Re_b:=\varepsilon /(\nu N^2)$, where $\nu$ is the kinematic viscosity and $N$ is the background buoyancy frequency, the correlation between suitably locally averaged density gradients and turbulence intensity and the importance of capturing the tails of the probability distribution functions of values of dissipation. Empirically modified versions of commonly used isotropic models for $\varepsilon$ and $\chi$ that depend only on vertical derivatives of density and velocity are proposed based on the asymptotic regimes $Re_b\ll 1$ and $Re_b\gg 1$, and serve as an instructive benchmark for comparison with the data-driven approach. When trained and tested on a simulation of stratified decaying turbulence which accesses a range of turbulent regimes (associated with differing values of $Re_b$), the PCNN outperforms assumptions of isotropy significantly as $Re_b$ decreases, and additionally demonstrates improvements over the fitted empirical models. A differential sensitivity analysis of the PCNN facilitates a comparison with the theoretical models and provides a physical interpretation of the features enabling it to make improved predictions.

Direct numerical simulations (DNS) were carried out to investigate flow control for transition delay using steady blowing/suction strips at the wall of a flared cone at Mach 6 and zero angle of attack. For the numerical investigations of the transition control strategy, the flared cone geometry and the flow conditions of the experiments in the Boeing/Air Force Office of Scientific Research (AFOSR) Mach 6 Quiet Tunnel (BAM6QT) at Purdue University were chosen. For the DNS, transition was initiated by introducing random disturbances at the inflow of the computational domain, emulating ‘natural’ transition in wind-tunnel experiments caused by free-stream noise. In both wind-tunnel experiments and numerical simulations, streamwise ‘hot’ streaks were found on the surface of the flared cone, which are caused by a nonlinear interaction of an axisymmetric second-mode wave and a pair of oblique waves of the same frequency (‘fundamental resonance’). The objective of the flow control strategy proposed here is to delay the transition onset, and thus mitigate the negative consequences associated with the nonlinear transition stages, i.e. the development of hot streaks and large wall-pressure amplitudes that were observed in experiments and DNS. Our previous so-called ‘controlled’ transition simulations have shown that flow control using steady blowing and suction strips can lead to a significant delay of the hot streak development on the surface of the flared cone. The simulation results presented in this paper show that this flow control strategy remains effective, even in a natural transition scenario characterized by broadband disturbances.

The influence of outer large-scale motions (LSMs) on near-wall structures in compressible turbulent channel flows is investigated. To separate the compressibility effects, velocity fluctuations are decomposed into solenoidal and dilatational components using the Helmholtz decomposition method. Solenoidal velocity fluctuations manifest as near-wall streaks and outer large-scale structures. The spanwise drifting of near-wall solenoidal streaks is found to be driven by the outer LSMs, while LSMs have a trivial influence on the spanwise density of solenoidal streaks, consistent with the outer LSM impacts found in incompressible flows (Zhou et al., J. Fluid Mech., vol. 940, 2022, p. A23). Dilatational motions are characterized by the near-wall small-scale travelling-wave packets and the large-scale parts in the outer region. The streamwise advection velocity of the near-wall structures remains at $16 \sim 18u_{\tau }$, hardly influenced by Mach numbers, Reynolds numbers and wall temperatures. The spanwise drifting of near-wall dilatational structures, quantified by the particle image velocimetry method, follows a mechanism distinct from solenoidal streaks. This drifting velocity is notably larger than those of the solenoidal streaks, and the influence of outer LSMs is not the primary trigger for this drifting.

Despite elegant theoretical descriptions and numerous potential applications in nature, the various features of wave turbulence and its associated transfers across scales are still debated. One reason is the difficulty, both experimentally and numerically, of studying a chaotic nonlinear wavy system with sufficient precision, over a sufficient duration, and with sufficient separation between the injection scale, the system scale and the dissipation scale, to dwell in detail on the description of a statistically converged state. The numerical study of Thomas & Ding (J. Fluid Mech., 2023, this issue) sheds new light on one key aspect of this global problem: the upscale transfer of waves in rotating shallow water equations, which is of relevance to atmosphere and ocean dynamics. Their results first confirm the theoretical predictions, with a robust inverse wave cascade, a predominant role of quartic resonances and a slope value of the wave spectrum in agreement with the expected range. But their deeper analysis of the results also highlights some weakness in the theoretical foundations, exhibiting strong intermittency and departure from ideal Gaussian statistics. This work thus calls for improved modelling, acknowledging that important large-scale climatic features could actually result from the piling up of small-scale waves, unresolved by typical Earth system models. Beyond rotating shallow water systems, this study more generically resonates with open questions in the field of wave turbulence and its geophysical applications, including recent research in deep rotating and stably stratified systems.

We investigate systematically the free settling of a single Platonic polyhedron in an unbounded domain filled with an otherwise quiescent Newtonian fluid. We consider a particle–fluid density mimicking a rock in water. Five Platonic polyhedrons of increasing sphericity are studied for a range of Galileo numbers $10 \leqslant \mathcal {G}a \leqslant 300$. We construct a regime map in the parameter space of Galileo number and particle volume fraction ($\mathcal {G}a, \phi$), highlighting how the angularity of the Platonic polyhedron impacts its settling path and the onset of instabilities. We find that the initial angular position solely affects the transient settling process. All the Platonic polyhedrons maintain a stable settling angular position at low $\mathcal {G}a$. Higher angularity leads to path unsteadiness at lower $\mathcal {G}a$. Path instability progresses from steady vertical to unsteady vertical, followed by oblique settling observed for highly spherical particles, but helical settling (HS) for more angular particles. The particle autorotation is found to be the pivotal factor influencing path instability and the regime transition of angular particles. Beginning in the unsteady oblique and helical regimes, particle autorotation becomes more prevalent, escalating further in the chaotic regime as $\mathcal {G}a$ increases. The particle angular velocity vector is shown to be predominantly situated in the horizontal plane. A thorough force balance in the horizontal plane reveals that the Magnus force is the primary driving force of the HS regime. Additionally, we establish two new empirical correlations to predict the particle settling velocity and the disturbed wake length that solely require the physical properties of the system ($\mathcal {G}a$ and $\phi$). Our numerical results suggest that an increase of the density ratio from $2$ to $3$ exerts only a marginal impact on the path instability of the most angular particle, the settling tetrahedron.

This study investigates the coherence of turbulent fluctuations in a turbulent vertical natural convection boundary layer immersed in a stably stratified medium (turbulent buoyancy layer). A turbulent buoyancy layer of a fluid having a Prandtl number of $0.71$ at a Reynolds number of $800$ is numerically simulated using direct numerical simulation. The two-point correlations reveal that the streamwise velocity fluctuations are coherent over large streamwise distances, with the length scale of the streamwise coherence being greater than the boundary layer thickness. This is due to large-scale motions (LSMs), similar to the LSMs observed in canonical wall-bounded turbulence despite the stark differences in flow dynamics. Both high-speed (positive) and low-speed (negative) streamwise velocity fluctuations form LSMs, with their streamwise length scales increasing with increasing wall-normal distance. High-speed LSMs are composed of upwash flow with high temperatures, while low-speed LSMs are composed of downwash flow with low temperatures. Both high-speed and low-speed LSMs meander appreciably in the streamwise direction, with the degree of meandering being correlated with the sign of the spanwise velocity fluctuations. The LSMs exhibit coherence across significant wall-normal distances and contribute significantly to the turbulence production in the outer layer. Examining the one-dimensional energy spectra of the turbulent buoyancy layer shows that the LSMs are the dominant energy-containing motions, implying that the length scale of the energy-containing range is of the order of boundary layer thickness. Notably, wall-normal velocity, spanwise velocity and buoyancy fluctuations do not form LSMs with streamwise length scales comparable to streamwise velocity fluctuations.

Particulate flows at moderate particle Reynolds numbers are important in critical engineering and geological applications. This experimental study explores neutrally buoyant suspensions in an outer-rotating coaxial rheometer for solid fractions, $\phi$, from 0.1 to 0.5, and particle Reynolds number, $Re$, from 0.5 to 800, covering laminar, transitional and turbulent regimes; $Re$ is defined in terms of the square of the particle diameter and the shear rate. For $0.1 < \phi < 0.4$ and $0.5 < Re <10$, the direct torque measurements normalised by the laminar flow torque, $M/M_{lam}$, are independent of $Re$, but depend on $\phi$. For the same range of $\phi$ and for $10< Re<100$, the normalised torques depend on both $\phi$ and $Re$, and show an increasing dependence on $Re$. As $Re$ increases, the flow transitions to turbulence. Small particles delay the turbulent transition for $\phi \leqslant 0.3$, while large particles augment the transition. A modified Reynolds number, $Re^\prime$, that depends linearly on the particle diameter and the maximum velocity, $U_{o}$, is introduced for both laminar and turbulent flows and shows a better correlation of the results as compared with $Re$. For $\phi = 50\,\%$, the normalised torque minus the torque at zero rotational speed is nearly independent of $Re^\prime$. Rheological models based on $Re^\prime$ and the Krieger–Dougherty relative viscosity are proposed in the laminar regime for $10< Re^\prime <500$; in the turbulent regime, a correlation is proposed in terms of $Re^\prime$ and $\phi$ for $1000< Re^\prime < 6000$.

The stability and postcritical behaviour of a horizontal flag undergoing gravity-induced deformation and periodic contact with a nearby horizontal rigid wall are experimentally investigated. The results elucidate the combined effects of gravity and contact on flutter, and reveal design principles for application to triboelectric energy harvesting. By varying the free-stream velocity, flag thickness and distance between the flagpole and the wall, the dynamics of the flag are classified into quasistatic equilibrium, flutter, partial contact and saturated contact modes. Considering the significance of gravitational effects, a new dimensionless flow velocity is proposed to identify the distribution of the dynamic modes, and its definition varies according to whether the wall is placed above or below the flag. The critical conditions for transitions between the dynamic modes are determined from the balance of fluid dynamic and gravitational effects. The distance from the flagpole to the wall is found to be more critical for transitions in the lower-wall configuration than in the upper-wall configuration. The peak contact force as well as the oscillation amplitude and frequency at postequilibrium exhibits remarkably different trends depending on the location of the wall. The peak contact force imposed on the wall by the fluttering flag weakens as the distance to the wall increases in the case of an upper wall, whereas it becomes stronger in the case of a lower wall.

A computational study was performed of the transport of both non-adhesive and adhesive particles in a porous bed with a body-centred cubic (BCC) structure. Pore-scale simulation of the flow within the porous bed was achieved through combining the immersed boundary method and the lattice Boltzmann method. Particle transport is computed using an adhesive discrete-element method based on a multi-time-scale soft-sphere model. The fluid flow results are validated by comparison with experimental data for dimensionless permeability of flow in a porous bed of spheres. For computations with non-adhesive particles, the particles are observed to drift to the centre of ‘channels’ in the BCC array, within which most of the fluid flow occurs. The mechanism of this inward drift was found to be related to the phenomenon of oscillatory clustering, which is an inertial drift mechanism observed for particles in a corrugated channel. A measure for particle drift into these channels was developed, and the time rate of change of this measure was found to compare closely with an approximate theoretical prediction based on oscillatory clustering theory. The drift measure was observed to be limited at long time by hold-up of outlier particles caught in long-duration collisions with the fixed bed particles in regions of low fluid velocity magnitude. Simulations with adhesive particles exhibited marked increase in collision duration, as well as inhibition of the tendency to drift toward the flow channels due to adhesive hold-up.

Fragmentation of a fluid body into droplets underlies many contamination and disease transmission processes where pathogens are transported in a liquid phase. An important class of such processes involves formation of a fluid ligament and its destabilization into droplets. Inertial detachment (Gilet & Bourouiba, J. R. Soc. Interface, vol. 12, 2015, 20141092) is one of these modes: upon impact on a sufficiently compliant substrate, the substrate's motion can transfer its impulse to a contaminated sessile drop residing on it. The fragmentation of the sessile drop is efficient at producing contaminated ejected droplets with little dilution. Inertial detachment, particularly from substrates of intermediate wetting, is also interesting as a fundamental fragmentation process on its own merit, involving the asymmetric stretching of the sessile drop under impulsive axial forcing with one-sided pinning due to the substrate's intermediate wetting. Our experiments show that the radius, $R_{tip}$, of the tip drop ejected become insensitive to the Bond number value for $Bo>1$. Here, $Bo$ quantifies the inertial effects via the relative axial impulsive acceleration compared with capillarity. The time, $t_{tip}$, of tip-drop breakup is also insensitive to $Bo$. Combining experiments, theory and validated numerics, we decipher the selection of $R_{tip}$ and its sensitivity to the surface-wetting and substrate foot dynamics. Using asymptotic theory in the large $Bo$ limit for which the thin-film/slender-jet approximations hold, we derive a reduced physical model that predicts $R_{tip}$ consistent with our experiments. Finally, we discuss how pathogen physical properties (e.g. wetting and buoyancy) within the sessile drop determine their distribution in the tip and secondary fragmentation droplets.

We study linear convective instability in a mushy layer formed by solidification of a binary alloy, cooled by either an isothermal perfectly conducting boundary or an imperfectly conducting boundary where the surface temperature depends linearly on the surface heat flux. A companion paper (Hitchen & Wells, J. Fluid Mech., 2025, in press) showed how thermal and salinity conditions impact mush structure. We here quantify the impact on convective instability, described by a Rayleigh number characterising the ratio of buoyancy to dissipative mechanisms. Two limits emerge for a perfectly conducting boundary. When the salinity-dependent freezing-point depression is large versus the temperature difference across the mush, convection penetrates throughout the depth of a high-porosity mush. The other limit, which we will call the Stefan limit, has small freezing-point depression and inhibits convection, which localises at onset to a high-porosity boundary layer near the mush–liquid interface. Scaling arguments characterise variation of the critical Rayleigh number and wavenumber based on the potential energy contained in order-one aspect ratio convective cells over the high-porosity regions. The Stefan number characterises the ratio of latent and sensible heats, and has moderate impact on stability via modification of the background temperature and porosity. For imperfectly conducting boundaries, the changing surface temperature causes stability to decrease over time in the limit of large freezing-point depression, but in the Stefan limit combines with the decreasing porosity to yield non-monotonic variation of the critical Rayleigh number. We discuss the implications for convection in growing sea ice.

We investigate the Faraday instabilities of a three-layer fluid system in a cylindrical container containing low-viscosity liquid metal, sodium hydroxide solution and air by establishing the Mathieu equations with considering the viscous model derived by Labrador et al. (J. Phys.: Conf. Ser., vol. 2090, 2021, 012088). The Floquet analysis, asymptotic analysis, direct numerical simulation and experimental method are adopted in the present study. We obtain the dispersion relations and critical oscillation amplitudes of zigzag and varicose modes from the analysis of the Mathieu equations, which agree well with the experimental result. Furthermore, considering the coupling strength of two interfaces, besides zigzag and varicose modes, we find a beating instability mode that contains two primary frequencies, with its average frequency equalling half of the external excitation frequency in the strongly coupled system. In the weakly coupled system, the $A$-interface instability, $B$-interface instability and $A$&$B$-interface instability are defined. Finally, we obtain a critical wavenumber $k_c$ that can determine the transition from zigzag or varicose modes to the corresponding $A$-interface or $B$-interface instability.

In this paper, we propose a series expansion of the baroclinic torque in low-Mach-number flows, so that the accuracy and universality of any buoyancy term could be examined analytically, and new types of buoyancy terms could be constructed and validated. We first demonstrate that the purpose of introducing a buoyancy term is to approximate the baroclinic torque, and straightforwardly the error of any buoyancy term could be defined with the deviation of its curl from the corresponding baroclinic torque. Then a regular perturbation method is introduced for the elliptic equation of the hydrodynamic pressure in low-Mach-number flows, resulting in a sequence of Poisson equations, whose solutions lead to the series representation of the baroclinic torque and the new types of buoyancy terms. It is found that the frame invariance of the momentum equation is maintained with one of the new types of buoyancy terms. With the error definition of buoyancy terms and the series representation of the baroclinic torque, the validity and accuracy of previous and new buoyancy terms are examined. Finally, numerical simulations confirm that, with a decreasing density variation or an increasing order of our new buoyancy term, the simplified equations can converge to the original low-Mach-number equations.

We investigate the instabilities and transition mechanisms of Boussinesq stratified boundary layers on sloping boundaries when subjected to oscillatory body forcing parallel to the slope. We examine idealized forms of boundary layers on hydraulically smooth abyssal slopes in tranquil mid- to low-latitude regions, where low-wavenumber internal tides gently heave isopycnals up and down adiabatic slopes in the absence of mean flows, high-wavenumber internal tides, shelf breaks, resonant tide–bathymetry interactions (critical slopes) and other phenomena associated with turbulence ‘hot spots’. In non-rotating low-Reynolds-number flow, increased stratification on the downslope phase has a relaminarizing effect, while on the upslope phase we find transition-to-turbulence pathways arise from shear production triggered by gravitational instabilities. When rotation is significant (low slope Burger numbers) we find that boundary layer turbulence is sustained throughout the oscillation period, resembling stratified Stokes–Ekman layer turbulence. Simulation results suggest that oscillating boundary layers on smooth slopes at low Reynolds number ($\textit {Re}\leqslant 840$), unity Prandtl number and slope Burger numbers greater than unity do not cause significant irreversible turbulent buoyancy flux (mixing), and that flat-bottom dissipation rate models derived from the tide amplitude are accurate within an order of magnitude.

Helicity plays a key role in the evolution of vortex structures and turbulent dynamics. The helicity dynamics and vortex structures in streamwise-rotating channel turbulence are discussed in this paper using the helicity budget equation and the differentiated second-order structure function equation of helicity. Generally, rotation and Reynolds numbers exhibit opposing effects on the interscale helicity dynamics and the vortices. Under the buffer layer, the positions of the helicity peaks are proportional to the ratio between the Reynolds and rotation numbers. The mechanism is related to the opposing effects of convection and rotation. Rotation directly affects the helicity balance through the Coriolis term and corresponding pressure term. In the buffer layer, the scale helicity is negative at small scales but positive at large scales, which is mainly induced by the spatial effects (the production and the spatial turbulent convection) but reduced by interscale cascades. Examination of structures reveals the close association between scale helicity and streaks, with streak lift angles exhibiting an increase with rotation and a decrease with Reynolds numbers. In the log-law layer, the Coriolis terms and corresponding pressure terms are proportional to the rotation numbers but remain independent of the Reynolds numbers. The negative scale helicity is forward cascaded towards small scales. Generally, spanwise vortices in the log-law layer are related to sweep events and forward cascades. Our findings indicate that these spanwise vortices are suppressed by rotation but recover with increasing Reynolds numbers, aligning with the effects observed in the scale helicity balance.

Experiments and numerical simulations of inertial particles in underexpanded jets are performed. The structure of the jet is controlled by varying the nozzle pressure ratio, while the influence of particles on emerging shocks and rarefaction patterns is controlled by varying the particle size and mass loading. Ultra-high-speed schlieren and Lagrangian particle tracking are used to experimentally determine the two-phase flow quantities. Three-dimensional simulations are performed using a high-order, low-dissipative discretization of the gas phase while particles are tracked individually in a Lagrangian manner. A simple two-way coupling strategy is proposed to handle interphase exchange in the vicinity of shocks. Velocity statistics of each phase are reported for a wide range of pressure ratios, particle sizes and volume fractions. An upstream shift of the Mach disk in the presence of particles reveals significant two-way coupling even at low mass loading. A semi-analytic model that predicts the extent of the Mach disk shift is presented based on a one-dimensional Fanno flow that takes into account volume displacement by particles and interphase exchange due to drag and heat transfer. The per cent shift in Mach disk is found to scale with the mass loading, nozzle pressure ratio and interphase slip velocity and inversely with the particle diameter.