The paper presents a simulation of the turbulent flow over and through a submerged aquatic canopy composed of 672 long, slender ribbons modelled as Cosserat rods. It is characterized by a bulk Reynolds number of 20 000, and a friction Reynolds number of 2638. Compared with a smooth turbulent channel at the same bulk Reynolds number, the canopy increases drag by a factor of 12. The ribbons are highly flexible, with a Cauchy number of 25 000, slightly buoyant, and densely packed. Their length exceeds the channel height by a factor of 1.6, while their average reconfigured height is only a quarter of the channel height. Different from lower-Cauchy-number cases, the movement of the ribbons, characterized by the motion of their tips, is very pronounced in the vertical direction, and even more in the spanwise direction, with root-mean-square fluctuations of the spanwise tip position 1.5 times the vertical ones. A canopy hull is defined to analyse the collective motion of the canopy and its interaction with the outer flow. Dominant spanwise wavelengths at this interface measure approximately one channel height, corresponding to twice the spacing of adjacent high- and low-speed streaks identified in two-point correlations of fluid velocity fluctuations. Conditional averages associated with troughs and ridges in the topography of the hull reveal streamwise-oriented counter-rotating vortices. They are reminiscent of the head-down structures related to the monami phenomenon in lower-Cauchy-number cases.

Predicting particle segregation has remained challenging due to the lack of a general model for the segregation velocity that is applicable across a range of granular flow geometries. Here, a segregation-velocity model for dense granular flows is developed by exploiting force balance and recent advances in particle-scale modelling of the segregation driving and drag forces over the entire particle concentration range, size ratios up to 3 and inertial numbers as large as 0.4. This model is shown to correctly predict particle segregation velocity in a diverse set of idealised and natural granular flow geometries simulated using the discrete element method. When incorporated in the well-established advection–diffusion–segregation formulation, the model has the potential to accurately capture segregation phenomena in many relevant industrial applications and geophysical settings.

Resolvent-based modelling and estimation is critically dependent on the nonlinear forcing input and hence understanding its role in the flow response is of great significance. This study quantifies the nonlinear forcing input in the resolvent formulation and investigates its characteristics for compressible turbulent boundary layers at Mach number 5.86 and friction Reynolds number 420 subject to adiabatic- and cold-wall conditions. Results show that, with the addition of the eddy viscosity to the resolvent operator, the cross-spectral density (CSD) of the forcing tends to exhibit a spatially uncorrelated distribution, which suggests that the spatial cross-coherence may be neglected and makes the modelling of the forcing input potentially easier. Aiming to quantify the different importance of each forcing component in generating turbulent fluctuations, contributions of the eddy-viscosity-corrected forcing to the flow responses are investigated through reduced-order analysis and matrix decomposition. The streamwise motions are almost insensitive to the temperature-related forcing, and can be oppositely influenced by the wall-normal and spanwise forcing components. By retaining only the diagonal components in the CSD of the forcing input, the assumption of forcing decorrelation in space and among components is also examined in the input–output framework. It is found that this simplified input is able to capture the dominant turbulence features and the local forcing is observed to cause inner-layer responses. That is, present results suggest adequate modelling of the CSD of the forcing can be achieved retaining only its diagonal components. On the basis of the current findings, the forcing input in the resolvent-based framework is thus modelled, with the wall-normal dependence and amplitude ratio between forcing components designed for compressible turbulent boundary layers. Through an algebraic Lyapunov equation, improved estimations of the statistical spectral densities of velocity and temperature fluctuations are finally obtained, in contrast to the results by simply assuming the forcing CSD to be an identity matrix.

Direct numerical simulations are carried out to investigate the underlying mechanism of the low-frequency unsteadiness of a transitional shock reflection with separation at $M=1.5$. To clarify the nonlinear mechanisms, the incoming laminar boundary layer is forced with two different arrangements of oblique unstable modes. Each wave arrangement is given by a combination of two unstable waves such that their difference in frequency falls in a low-frequency range corresponding to a Strouhal number (based on the length of interaction) of 0.04. This deterministic forcing allows the introduction of nonlinearities, and high-order statistical tools are used to identify the properties of quadratic couplings. It is found that the low-frequency unsteadiness and the transition to turbulence are decoupled problems. On the one hand, the unstable modes of the boundary layer interact nonlinearly such that energy cascades to higher frequencies, initiating the turbulent cascade process, and to lower frequencies. On the other hand, the low-frequency quadratic coupling of the oblique modes is found to be responsible for low-frequency unsteadiness affecting the separation point. The direction of the quadratic interactions is extracted and it is shown that, in the presence of low-frequency unsteadiness, these interactions enter the separated zone just before reattachment and travel both downstream and upstream, extending beyond the separation point, hence feeding the low-frequency bubble response. In addition to the two main arrangements of oblique modes, two other combinations are analysed, including multiple oblique waves and streaks. Interestingly, their inclusion did not alter the low-frequency unsteadiness phenomenon. Furthermore, the effect of the forcing difference frequency is examined and it is shown that the breathing phenomenon is sensitive to the range of frequencies present in the system due to a low-pass filter effect.

Mass transport induced by group-forced subharmonic waves (infragravity waves) is investigated in the present study. A theoretical solution for subharmonic waves’ kinematic contributions to fourth-order mass transport and drift velocity has been proposed for any depth and bandwidth for the first time. This model is validated using particle-tracking simulations driven by the flow field generated by the SWASH. The subharmonic-induced mass transport solution is a weighted sum of the subharmonic velocity variance spectrum and velocity skewness bispectrum due to the triad-difference interaction among two primary and one subharmonic components. For narrow-banded waves with long wave group relative to depth, the weightings become independent of spectral components, and the solution is recovered in the time domain. Two mechanisms contributing to mass transport were identified: a forward drift resulting from self-interaction similar to Stokes drift, and a depth-decaying backward drift induced by negative subharmonic velocity skewness due to the anti-phase coupling between subharmonics and wave groups. For narrow-banded waves the forward transport surpasses the backward transport for kh< 0.72, where k is the short wave wavenumber and h is the water depth. For other waves, the critical kh for this phenomenon decreases with increasing wave period and bed slope and decreasing bandwidth. At greater depths or steeper bed slopes, near-surface backward transport predominates over forward transport; at shallower depths or gentler slopes, forward transport is dominant throughout the water column. Although smaller than Stokes transport by short waves, the subharmonic wave-induced mass transport can affect the long-term trajectory of a floating and suspended particle. This study provides the first evidence and insight for the influences of group-forced subharmonics on vertically varying mass transport from the ocean surface to seabed in coastal environments.

Convection in planetary environments is often modelled using stress-free boundary conditions, with diffusion-free geostrophic turbulence scalings frequently assumed. However, key questions remain about whether rotating convection with stress-free boundary conditions truly achieves the diffusion-free geostrophic turbulence regime. Here, we investigated the scaling behaviours of the Nusselt number ($Nu$), Reynolds number (${Re}$) and dimensionless convective length scale ($\ell /H$, where $H$ is the height of the domain) in rotating Rayleigh–Bénard convection under stress-free boundary conditions within a Boussinesq framework. Using direct numerical simulation data for Ekman number $Ek$ down to $5\times 10^{-8}$, Rayleigh number $Ra$ up to $5\times 10^{12}$, and Prandtl number $Pr = 1$, we show that the diffusion-free scaling of the heat transfer $Nu - 1 \sim Ra^{3/2}\, Pr^{-1/2}\, Ek^2$ alone does not necessarily imply that the flow is in a geostrophic turbulence regime. Under the stress-free conditions, ${Re}$ and $\ell /H$ deviate from the diffusion-free scalings, indicating a dependence on molecular diffusivity. We propose new non-diffusion-free scaling relations for this diffusion-free heat transfer regime with stress-free boundary conditions: $\ell /H \sim Ra^{1/8}\, Pr^{-1/8}\, Ek^{1/2}$ and ${Re} \sim Ra^{11/8}\, Pr^{-11/8}\, Ek^{3/2}$. Our findings highlight the need to assess both thermal and dynamic characteristics to confirm geostrophic turbulence.

In turbulent pipe flows, drag-reducing polymers are commonly used to reduce skin-friction drag; however, predicting this reduction in industry applications, such as crude oil pipelines, remains challenging. The skin-friction coefficient ($C_f$) of polymer drag-reduced turbulent pipe flows can be related to three dimensionless parameters: the solvent Reynolds number ($Re_s$), the Weissenberg number ($Wi$) and the ratio of solvent viscosity ($\eta _s$) to zero-shear-rate viscosity ($\eta _0$), denoted as $\beta$. The function that relates these four dimensionless numbers was determined using experiments of various pipe diameters ($D$), flow velocities ($U$) and drag-reducing polyacrylamide solutions. The experiments included measurements of streamwise pressure drop ($\Delta P$) for determining $C_f$, and measurements of shear viscosity ($\eta$) and elastic relaxation time ($\lambda$). This experimental campaign involved 156 flow conditions, each characterised by distinct values for $C_f$, $Re_s$, $Wi$ and $\beta$. Experimental results demonstrated good agreement with the relationship: $C_f^{-1/2} = \widehat {A}\log _{10}(Re_sC_f^{1/2})+\widehat {B}$, where $\widehat {A} = 27.6(Wi \beta )^{0.346}$ and $\widehat {B} = 122/15-58.9(Wi \beta )^{0.346}$. Based on this relationship, onset and maximum drag reduction are predicted to occur when $Wi \beta$ equals $3.76 \times 10^{-3}$ and $3.40 \times 10^{-1}$, respectively. This function can predict $C_f$ of dilute polyacrylamide solutions based on predefined parameters (bulk velocity, pipe diameter, density, solvent viscosity) and two measurable rheological properties of the solution (shear viscosity and elastic relaxation time) with an accuracy of $\pm 9.36$ %.

Turbulence closures are essential for predictive fluid flow simulations in both natural and engineering systems. While machine learning offers promising avenues, existing data-driven turbulence models often fail to generalise beyond their training datasets. This study identifies the root cause of this limitation as the conflation of generalisable flow physics and dataset-specific behaviours. We address this challenge using symbolic regression, which yields interpretable, white-box expressions. By decomposing the learned corrections into inner-layer, outer-layer and pressure-gradient components, we isolate universal physics from flow-specific features. The model is trained progressively using high-fidelity datasets for plane channel flows, zero-pressure-gradient turbulent boundary layers (ZPGTBLs), and adverse pressure-gradient turbulent boundary layers (PGTBLs). For example, direct application of a model trained on channel flow data to ZPGTBLs results in incorrect skin friction predictions. However, when only the generalisable inner-layer component is retained and combined with an outer-layer correction specific to ZPGTBLs, predictions improve significantly. Similarly, a pressure-gradient correction derived from PGTBL data enables accurate modelling of aerofoil flows with both favourable and adverse pressure gradients. The resulting symbolic corrections are compact, interpretable, and generalise across configurations – including unseen geometries such as aerofoils and Reynolds numbers outside the training set. The models outperform baseline Reynolds-averaged Navier–Stokes closures (e.g. the Spalart–Allmaras and shear stress transport models) in both a priori and a posteriori tests. These results demonstrate that explicit identification and retention of generalisable components is key to overcoming the generalisation challenge in machine-learned turbulence closures.

Compressible jets impinging on a perpendicular surface can produce high-intensity, discrete-frequency tones. The character of these tones is a function of nozzle shape, jet Mach number, impingement-plate geometry, and the distance between nozzle and plate. Though it has long been recognised that these tones are associated with a resonance cycle, the exact mechanism by which they are generated has remained a topic of some debate. In this work, we present evidence for a number of distinct tone-generation mechanisms, reconciling some of the different findings of prior authors. We demonstrate that the upstream-propagating waves that close resonance can be confined within the jet, or external to it. These waves can be either weak and relatively linear, or strong and nonlinear from their inception. The waves can undergo coalescence or merging, and in some configurations, pairs of waves rather than singletons appear. We discuss both historical and new evidence for multiple distinct processes by which upstream-propagating waves are produced: direct vortex sound, shock leakage, wall-jet-boundary fluctuations, and wall-jet shocklets. We link these various mechanisms to the disparate collection of upstream-propagating waves observed in the data. We also demonstrate that multiple mechanisms can be provoked by a single vortex, providing an explanation as to why sometimes pairs of waves or merging waves are observed. Through this body of work, we demonstrate that rather than being in opposition, the various pieces of past research on this topic were simply identifying different mechanisms that can support resonance.

The paper by Pružina et al. (2025) J. Fluid Mech. 1009, sheds new light on the physical processes responsible for the formation of distinct layers in double-diffusive convection. Towards this end, it discusses direct numerical simulation results within the framework of sorted buoyancy coordinates. In particular, it demonstrates that the eddy diffusivity is negative everywhere, including in the interior of the well-mixed layers. This approach holds promise for analysing other, closely related, flow configurations that give rise to the emergence of pronounced layering features.

We present a flexible, multilayer fabric strain sensor composed of a carbon fabric layer sandwiched between elastic bands. The sensor achieved a gauge factor of 3.4 and maintained its durability up to 635% strain. Its uniform graphite layer enabled reliable fabrication and easy integration into wearable formats. Performing well on commercial gloves and bands, the sensor effectively captured strain variations during body movement and enabled wireless transmission for real-time monitoring. Distinct resistance patterns were recorded for various body motions such as walking, jogging, jumping, and knee bending with a clear separation between high- and low-intensity activities. The overall design supports scalable fabrication and practical integration into wearable systems.

Dot array deposition through electrohydrodynamic (EHD) printing is widely used for high resolution and material utilization advantages. However, the conventional printing method is subject to a printing frequency limit known as the capillary frequency of the meniscus oscillation, where the jet directly contacts the substrate. This makes the printing frequency of EHD printing maintain at a low level and that is difficult to improve. In this work, a method for high-frequency EHD printing through continuous pinch-off is proposed. The characteristic frequency is broken through. A model is established to reveal the printing mechanism by combining the Poisson–Nernst–Planck equation and the phase field method. The unreal charge leakage is prevented by constructing a transition function for the fluid’s properties. The stability of the Taylor cone’s deformation and the droplets’ generation is studied. The measurement criterion for printing frequency is determined. The suitable printing height that can prevent the jet from directly contacting the substrate is obtained by investigating its influence on the printing states and frequency. The phase diagram considering the liquid’s conductivity and viscosity is presented to distinguish whether the printing is based on the end-pinching or Rayleigh–Plateau instability. The influence of the conductivity, viscosity, flow rate and printing voltage on the printing frequencies is studied quantitatively. Finally, scaling laws for printing frequency are proposed by theoretical analyses and summarizing the numerical data. This work could be beneficial for further enhancing the printing frequency of EHD printing.

Breaking wave impacts on rigid structures have been extensively studied, yet the role of structural elasticity in shaping the impact and response remains insufficiently understood. In this study, we experimentally investigate the hydroelastic behaviour of a vertical cantilever plate subjected to multimodal solitary breaking wave impacts. The plate is mounted near the still water level on a 1 : 10 sloping beach, and the wave height-to-depth ratio ($H/h$) is varied from 0.15 to 0.40 to systematically control the impact type from non-breaking to highly aerated wave impacts. We show that aeration significantly affects hydroelastic impacts. The spatio-temporal extent of the impact pressure on the elastic plate increases with air entrapment, while the peak pressure becomes highly sensitive as the wave approaches the flip-through regime. Pressure oscillations associated with bubble formation induce high-frequency structural vibrations, particularly under low-aeration conditions. Furthermore, we find that the elasticity has a limited effect on the peak pressure, impact duration and impulse, but increases the maximum quasi-hydrostatic force on the plate for the scenarios investigated. Following the impact, two distinct free-top deflections are identified, i.e. a deflection $\Delta x_{\textit{imp}}$ with high acceleration induced by the impact pressure and a deflection $\Delta x_{{hp}}$ with high magnitude caused by the maximum quasi-hydrostatic pressure. These deflections scale with the Cauchy number as $\Delta x_{\textit{imp}}/l \sim Ca_{\textit{imp}}/6$ and $\Delta x_{{hp}}/l \sim Ca_{{hp}}/12$ (where l is the plate length), exhibiting parabolic and linear trends with $H/h$, respectively. This work presents a benchmark dataset and introduces a predictive law for structural deflection, providing practical insights into hydroelastic effects across various impact regimes.

The effect of Stokes number on turbulence modulation in particle-laden channel flow is investigated through four-way coupled point-particle direct numerical simulations, with the mass loading fixed at 0.6 and the friction Stokes number $St^+$ varying from 3 to 300. A full transition pathway is observed, from a drag-enhanced to a drag-reduced regime, eventually approaching the single-phase state as $St^+$ increases towards 300. A set of transport equations for the particle phase is derived analytically to characterise the interphase coupling, within the framework of the point-based statistical description of particle-laden turbulence. By virtue of this, two dominant mechanisms are identified and quantitatively characterised: a positive, particle-induced extra transport that decreases monotonically with increasing $St^+$, and a negative, particle-induced extra dissipation that varies non-monotonically with $St^+$. The coupling of these two mechanisms leads to a direct contribution of the particle phase to the shear stress balance, the turbulent kinetic energy budgets and the Reynolds stress budgets. Consequently, as $St^+$ increases, the self-sustaining cycle of near-wall turbulence transitions from being augmented to being suppressed and, eventually, returns to the single-phase state. This gives rise to an indirect effect, manifested as a non-monotonic modulation of Reynolds shear stress and turbulence production rate. Taken together, complex interplays between particle-modified turbulent transport, particle-induced extra transport and extra dissipation are analysed and summarised, providing a holistic physical picture composed of consistent interpretations of turbulence modulation induced by small heavy particles.