Shock–boundary-layer interactions on hypersonic cone-step flows exhibit a range of intrinsic unsteady behaviours, from shear-layer oscillations to large-scale pulsations. This work investigates the unsteadiness in a cone-step geometry at Mach 6 under quiet-flow conditions at different free-stream Reynolds numbers using time-resolved schlieren imaging and spectral proper orthogonal decomposition. Experimental results are compared with high-fidelity axisymmetric and three-dimensional simulations. Results demonstrate regime transition in the parameter space, across the unsteadiness boundary, all the way from shear-layer breakdown to shock system oscillations and ultimately to large-amplitude pulsations. The dominant mode in the experiments and the simulations corresponds to a Strouhal number St $\approx 0.17$ for small oscillations reducing to St $ \approx 0.13$ for large pulsations. A detailed description of the unsteady shock dynamics and an analysis of the nonlinear limit cycle is presented.

Bed shear stress is a key parameter governing sediment transport and fluxes at the sediment–water interface. In vegetated channels, predicting bed shear stress, especially for rough beds, remains a challenge. This study developed a unified theoretical model for bed shear stress that smoothly spans conditions from bare bed to vegetated bed for both smooth and rough beds. Building on phenomenological turbulence theory, the model relates bed shear stress to the characteristic velocities of the larger energy-containing eddies and the smaller, near-bed eddies, with the new assumption that the bottom boundary layer (BBL) thickness controls the larger, energy-containing eddy length scale. The BBL was defined as the region within which the bed shear stress contributed significantly, compared to vegetation drag, and a force balance predicted that the BBL thickness scales with the ratio of bed shear stress to vegetation drag. In the limit of zero vegetation density, the BBL thickness equals the water depth, and the bed shear stress model reduces to the classical bare bed formulation. With increasing vegetation density (drag), the thickness of the boundary layer decreases, and the bed friction coefficient increases, which is consistent with previous observations. For rough beds, the bed friction coefficient increases with bed roughness, but is not dependent on the mean velocity. In contrast, for smooth beds, the bed friction coefficient decreases with increasing mean velocity. The coupled models for bed shear stress and BBL thickness were compared against 114 physical and numerical experiments from multiple previous studies.

Kinetic theory offers a promising alternative to conventional turbulence modelling by providing a mesoscopic perspective that naturally captures non-equilibrium physics such as non-Newtonian effects. In this work, we present an extension and theoretical analysis of the kinetic model for incompressible turbulent flows developed by Chen et al. (Atmosphere, 2023, vol. 14(7), p. 1109), constructed for unbounded flows. The first extension is to reselect a relaxation time such that the turbulent transport coefficients are obtained consistently and better align with well-established turbulence theory. The Chapman–Enskog (CE) analysis of the kinetic model reproduces the linear eddy-viscosity and gradient diffusion models for Reynolds stress and turbulent kinetic energy flux at the first order, and yields nonlinear eddy-viscosity and closure models at the second order. In particular, a previously unreported CE solution for turbulent kinetic energy flux is obtained. The second extension is to enable the model for wall-bounded turbulent flows with preserved near-wall asymptotic behaviours. This involves developing a low-Reynolds-number model incorporating wall damping effects and viscous diffusion, with boundary conditions enabling both viscous sublayer resolution and wall function application. Comprehensive validation against experimental and direct numerical simulation data for turbulent Couette flow demonstrates excellent agreement in predicting mean velocity profiles, skin friction coefficients and Reynolds shear-stress distributions, although the near-wall-normal stress anisotropy is underestimated. The results show that averaged turbulent flow behaves similarly to rarefied-gas flow at finite Knudsen number, capturing non-Newtonian effects beyond linear eddy-viscosity models. This kinetic model provides a physics-based foundation for turbulence modelling with reduced empirical dependence.

This study elucidates the influence of liquid viscosity on the hydrodynamics of simultaneous and non-simultaneous droplet-pair impacts on solid substrates. Using synchronised high-speed imaging and quantitative analysis, the spreading dynamics of droplet lamellae and their interaction-driven central sheet evolution are examined across a range of viscosities from 1.01 to 91.46 mPa s, representing Ohnesorge numbers of 0.002–0.177, under controlled impact Weber numbers in the range of 81–131 and dimensionless inter-droplet spacings in the range of 1.43–1.85. The findings reveal that increasing viscosity results in thicker lamella fronts, reduced spreading and a lower maximum central sheet height. In addition, the central sheet morphology transitions from ‘semilunar’ sheets to ephemeral liquid bumps, accompanied by suppressed capillary waves and reduced rim instabilities. A novel scaling law is derived for the maximum sheet extension, demonstrating its robust applicability to both simultaneous and non-simultaneous impacts of droplet pairs across varying viscosities and impact conditions. Furthermore, distinct morphological differences emerge between simultaneous and non-simultaneous impacts, primarily governed by lamella–lamella interactions and the momentum transfer dynamics. These findings enhance our understanding of the interplay between viscous and inertial forces in droplet-pair impacts, offering valuable insights for optimising spray-based technologies and multiphase fluid systems.

The active-layer model used to account for mixed-size sediment morphodynamic processes may be ill-posed under certain circumstances. Well-posedness guarantees the existence of a unique solution continuously depending on the problem data. When a model becomes ill-posed, infinitesimal perturbations to a solution grow infinitely fast. Apart from the fact that this behaviour cannot represent a physical process, numerical simulations of an ill-posed model continue to change as the grid is refined. For this reason, ill-posed models cannot be used as predictive tools. There exists a regularisation strategy based on a preconditioning method that guarantees that the one-dimensional active-layer model is well-posed. Here, we show that the extension of this strategy to two dimensions does not regularise the model and we propose a different regularisation strategy based on diffusion that guarantees that both the one-dimensional and two-dimensional active-layer models are well-posed. We implement the strategy in Delft3D Flexible Mesh and show an application.

This study investigates convection in a non-isothermal spherical Taylor–Couette flow (sTC) under the influence of the dielectrophoretic (DEP) force. The convective flow is driven by differential rotation of the inner and outer boundaries rotating with $\varOmega$ and $\Delta {\varOmega }$ in combination of an electric tension applied between both shells to induce thermo-electrohydrodynamic (TEHD) convection. To understand the interaction between DEP force-driven and rotation-driven mechanisms, we first analysed TEHD convection and non-isothermal sTC flow independently. For the TEHD case, we establish scaling relations for heat transport by expressing the Nusselt number, ${\textit{Nu}}$, as a function of the electric Rayleigh number, ${\textit{Ra}}_{{E}}$, and the kinetic energy density, $\tilde {E}_k$. These relations are evaluated against classical models of convection to assess consistency and deviations. A similar approach was applied to the non-isothermal sTC flow in the absence of the DEP force, where we identified axisymmetric and non-axisymmetric flow regimes which were classified by ${\textit{Nu}}$, $\tilde {E}_k$ and $\Delta {\varOmega }$, and developed corresponding scaling relations. When both mechanisms were active, ${\textit{Nu}}$ generally increased, however, the DEP force locally suppressed angular momentum transport, especially near the equator. This interplay revealed three distinct regimes: (A) DEP force-dominated TEHD convection, (C) rotation-dominated non-isothermal sTC flow and (B) a transitional regime with reduced heat transport. A decomposition of a derived inflow Nusselt number, ${\textit{Nu}}^q$, based on conductive and convective contributions, further elucidated the underlying heat transport mechanism.

Dispersion in turbulent flows is of broad interest in engineering and environmental processes, particularly for rivers, lakes and oceanic water bodies. Based on our streamwise dispersion model grounded in a Lagrangian perspective of convection–diffusion dynamics (Guan & Chen, 2024, J. Fluid Mech., vol. 980, A33), this work presents a comprehensive solution that consistently unifies dispersion across the Reynolds number spectrum, bridging laminar and turbulent regimes. The streamwise dispersion mechanism is general across time scales, yet its statistical behaviour cannot be fully described using conventional coarse-grained moments averaged over cross-sections. While classical drift–diffusion models that are effective for long-time asymptotics fail to capture the turbulent dynamics of the pre-asymptotic phase, our analytical model enables a complete spatio-temporal characterisation of concentration, and reveals how local statistics evolve towards their asymptotic, coarse-grained limits. Through asymptotic expansions and eigenfunction analysis, we quantify the time-dependent behaviour of phenomenological dispersion coefficients, and distinguish between local and mean statistics, which diverge significantly during the pre-asymptotic phase. The early regime exhibits robust features, including an overshoot in local dispersivity, asymmetric long tails in mean concentration, and island-shaped solute accumulation near the free surface. Three regimes are identified in the evolution of the local concentration: (i) an initially uniform line source, (ii) a transitional logarithmic profile shaped by vertical shear, and (iii) an emergent Gaussian dispersion regime approaching vertical uniformity. Comparisons of both local and mean concentration demonstrate quantitative agreement with finite difference and Monte Carlo simulations across all regimes. These findings clarify the interplay between shear and turbulent diffusion, laying a foundation for addressing more intricate and physically significant transport problems.

This study extends the two-step stochastic metafrontier approach to multiple-output farms to investigate technology gap in the dryland areas of West Africa using cross-sectional data from randomly selected farmers in Ghana, Mali and Nigeria. Results show that Nigeria has a higher average metatechnical efficiency (MTE) with on average 64% for Nigeria, 46% for Mali and 28% for Ghana with the same trends for the metatechnology ratios (MTR). This implies that, from the dryland farming metafrontier, farmers in Ghana are technically inefficient than those in Mali, and both are technically inefficient than farmers in Nigeria. The comparison of the MTR with the country specific efficiency (TE) shows that, technical inefficiency with respect to the metafrontier in Ghana and Mali comes mainly from the technology gap than operating inefficiency, whereas it comes primarily from production inefficiency in Nigeria. Finally, analysis of the sources of inefficiencies suggests that country-specific policies and programmes that help reduce the gap and targeting the important drivers of efficiency such as household size and access to credit are needed to improve technical efficiency. Overall, this study shows that there is economic and food security gain in enhancing farmer managerial performance strategies in the dryland zone.

In a previous study, we have proposed a mechanism for simultaneous reduction of drag and lift by half-rotation at moderately low Reynolds numbers. The axis of rotation ($z$) is perpendicular to both the drag ($x$) and lift ($y$) directions, i.e. the rotation is transverse to the incoming flow direction. Under laminar flow conditions, force-element analysis indicates that a partially rotating sphere can significantly reduce both drag and lift with suppression of vortex shedding. This study extends investigation of the same mechanism of half-rotating a sphere to the turbulence regime at a Reynolds number $Re = 1 \times 10^4$. Similar to the laminar case, half-rotation of the sphere introduces a significant negative velocity drag term, which effectively counteracts the rapid increase in the volume- and surface-vorticity drag terms. Numerical simulations with delayed detached eddy simulation, aided by direct numerical simulation, show that the drag coefficient decreases monotonically with increasing the non-dimensional rotational speed $\alpha$, even becoming negative at $\alpha =10$, while the lift and side-force coefficients remain small for all $\alpha$. However, in contrast to laminar conditions, the turbulent regime is characterised by an earlier onset of shear-layer instabilities, which accelerates the transition of the wake into a fully turbulent state. The relative importance of volume- and surface-vorticity contributions to the drag and lift is the most outstanding difference between the laminar and turbulent flows. In turbulent flow, simultaneous reduction of drag and lift is more pronounced as the contributions of volume- and surface-vorticity lift terms almost cancel each other exactly. These mechanisms and characteristics are systematically compared with those observed in the flow around a fully rotating sphere at the same Reynolds number in terms of vorticity structures, force elements, pressure distributions as well as surface-vorticity distributions.

The anomalous scaling of passive scalar fluctuations is experimentally investigated in turbulent pipe flow with a Taylor-microscale Péclet number of $\mathcal{O}(10^5)$, where the turbulence is known to deviate from the homogeneous isotropic turbulence. The scalar structure functions and intermittency in the mixing are examined. The experimental results consolidate that the scaling exponents of scalar structure functions saturate at high-order even moments, evidenced previously in homogeneous isotropic turbulence with a Taylor-microscale Péclet number of $\mathcal{O}(10^3)$. The saturation scaling exponent decreases to approach unity as the Taylor-microscale Péclet number increases. This saturation scaling exponent is further corroborated by the fractal codimension of sharp scalar fronts.

An experimental study is conducted investigating the characteristics of tones due to ‘guided jet waves’ in the forward arc of jet noise radiation fields. These spectral peaks disappear in the far acoustic field in the aft and sideline directions, which is why they went unnoticed in decades of jet noise measurements. However, it is clearly shown in the present study that they radiate in the forward arc (shallow upstream polar locations in the approximate range θ < 45°). Jet noise spectra in the forward arc, data on which had been lacking in the literature, are not smooth but are characterised by these peaks. This is found for high subsonic to supersonic jets up to the highest jet Mach number covered in the experiment (MJ ≈ 1.9), for round and rectangular as well as convergent–divergent nozzles. In heated jets, these spectral peaks appear to get weakened especially around transonic conditions, whereas they clearly persist in supersonic conditions.

Elastoviscoplastic (EVP) fluids, characterised by the coexistence of elastic, viscous and yield-stress properties, play a central role in diverse applications, including drug delivery, 3D printing and hydraulic fracturing. These fluids often transport non-spherical particles whose migration dynamics strongly influences flow behaviour. In this work, we employ interface-resolved direct numerical simulations to investigate the migration and orientation dynamics of finite-size spheroidal particles suspended in EVP duct flows across a wide range of governing parameters. Our results show that the equilibrium position and orientation of the particles are influenced significantly by both their aspect ratio and the carrier fluid rheology. In Saramito fluids, spheroidal particles migrate towards the duct centre and align along the duct diagonals in the presence of inertia. At sufficiently high elasticity, they penetrate the central plug and reach the duct core, irrespective of their initial position or shape. At lower elasticities, where larger plug regions persist, interactions with the plug alter the angular dynamics of the particles, leading to unsteady, quasi-periodic tumbling and spinning motions. In contrast, in Saramito–Giesekus fluids, the interplay between inertial forces, shear-thinning plastic viscosity and yield stress drives particles towards the duct corners, aligning them perpendicular to the duct diagonals. In semi-dilute suspensions, flattened particles maintain a greater distance from the walls, whereas their spherical counterparts tend to cluster directly at the corners. These findings reveal complex migration and orientation behaviours unique to EVP media and suggest new opportunities for geometry-based particle separation in microfluidic applications.

The nonlinear evolution of free-stream vortical disturbances entrained in the entrance region of a channel is investigated using asymptotic and numerical methods, building on the linear framework developed by Ricco & Alvarenga (2021 J. Fluid Mech., vol. 927, A18). The focus is on low-frequency disturbances that induce streamwise-elongated structures at Reynolds numbers for which the entrance flow is locally stable according to classical linear stability theory. The perturbation flow along the channel entrance is generated by free-stream vortical disturbances located at the channel inlet. These disturbances are symmetric or antisymmetric with respect to the centreplane and their amplitude is sufficiently intense to provoke nonlinear interactions within the channel. The formation and evolution of the perturbation flow are described by the nonlinear unsteady boundary-region equations. Combined with physically realistic initial conditions, the resulting initial-boundary-value problems are solved numerically using a streamwise integration method. A parametric study is conducted to elucidate how the nonlinear channel flow is influenced by the Reynolds number and the inlet-disturbance properties, i.e. the amplitude and the streamwise, wall-normal and spanwise wavelengths. Nonlinearity is found to stabilise the intense algebraic growth and to drive the formation of elongated channel-entrance structures that span the entire cross-section. These structures, characterised by low- and high-speed regions and streamwise vortices, meander along the streamwise direction and persist even when the base flow is fully developed. They exhibit a half-turn rotational symmetry with respect to the vortex centres. These properties emerge downstream regardless of the symmetry of the initial perturbation flow, provided nonlinear interactions are sufficiently intense. The occurrence of travelling waves is detected sufficiently downstream, and their similarity to those found in the fully developed region by other researchers is discussed. Our results show good agreement with theoretical predictions, numerical results and experimental measurements for both the mean flow and the perturbation flow.

The present work numerically investigates the dynamics of inclined thin flexible plates in oscillatory flows to assess the effects of bending stiffness, wave orbital excursion and inclination angle on plate deflection, reconfiguration, drag force and energy conversion. Four distinct structural response modes are identified, together with their transition conditions. An analytical expression for the lift coefficient of inclined rigid plates in the oscillatory flow is derived. By combining the drag and lift coefficients, we propose a modified Cauchy number, which quantitatively reveals the inclination effect from both geometric and mechanical perspectives. The dynamic behaviours of the plate deflection can be separated into two states. In the fully reconfigured state, a balance is achieved between elastic restoring force and hydrodynamics-driven force. Based on moment and energy balance, we derive a scaling law incorporating the modified Cauchy number, which accurately predicts the variation of tip deflection. In the passive movement state, the flexible plate moves passively along the flow, and its tip deflection saturates to the order of wave orbital excursions. A pronounced drag reduction is induced by the plate reconfiguration, following a $-1$ scaling law with a combined parameter, which is explained by the model of effective plate length. The energy conversion from fluid kinetic energy to structural elastic strain energy first increases and then decreases with increasing flexibility, yielding an optimal energy conversion efficiency. The modified Cauchy-number-based scaling law accurately predicts the averaged elastic energy growth and critical conditions for optimal energy conversion through a time scale competition mechanism.