This study examines the cross-flow vortex-induced vibration (VIV) of a circular cylinder in combined current–oscillatory inflows, revealing a distinct multi-frequency response characterised by beat-like modulation. Systematic water-channel experiments were conducted across a range of reduced velocities, inflow oscillation intensities and frequency ratios to investigate the synchronisation mechanisms among inflow velocity variations, cylinder motion and hydrodynamic loading. Results show that the presence of oscillatory inflow can lead to significant deviations of vibration amplitudes from quasi-steady predictions within the upper-branch regime. At a given reduced velocity, the cylinder motion is dominated by a primary frequency component similar to that observed in steady flow, but accompanied by two secondary components. The contributions of these supplementary frequencies increase with inflow oscillation intensity but diminish as the oscillation frequency rises. Analysis of time-varying hydrodynamic forces reveals that, in the upper-branch regime, the vortex-force phase angle deviates substantially from quasi-steady estimation based on instantaneous reduced velocity, which is associated with non-quasi-steady vortex-shedding patterns. Particle image velocimetry measurements reveal that when the minimum vortex-force phase angle lies between 0$^\circ$ and 180$^\circ$ over the inflow oscillation cycle, a mixed vortex-shedding mode emerges. This mode is characterised by a vortex sequence resembling the ‘2P’ (two-pair) shedding pattern but with negligible secondary vortices, occurring predominantly during intervals of low inflow velocity. A theoretical framework incorporating nonlinear damping and excitation coefficients assuming quasi-steady response well predicts VIV amplitudes and elucidates the influence of inflow oscillation intensity and frequency on the emergence of supplementary vibration frequencies.

This study experimentally investigates bubble size evolution and void fraction redistribution in an unexplored, coalescence-dominated regime of a decaying turbulent bubbly flow. The flow is generated downstream of a regenerative pump in a duct, with bulk Reynolds number (Re) $\sim \mathcal{O}(10^5)$, Taylor-scale Reynolds number (Re$_\lambda$) $\sim \mathcal{O}(10^3)$ and void fraction ($\phi$) $\sim \mathcal{O}(1\,\%)$, where the inlet turbulence is extremely intense (turbulence intensity $\gt 30\,\%$) but decays rapidly along the duct. Shadowgraph imaging and particle shadow velocimetry are used for measurements. The experimentally obtained turbulent dissipation in the duct flow decays as $\varepsilon \sim \mathcal{L}^{-2}$, where $\mathcal{L}$ is the axial position, in close agreement with the homogeneous isotropic turbulence prediction of $\varepsilon \sim \mathcal{L}^{-2.2}$. High-speed imaging and statistical analysis reveal that bubble coalescence dominates over breakup across most of the domain, leading to monotonic growth in the Sauter mean diameter ($d_{32}$) and progressive broadening of the bubble size distribution. The normalised extreme-to-mean diameter ratio ($\mathcal{D}$) increases axially and asymptotically from ${\sim} 1.9$ (breakup regime) and saturates at ${\sim} 2.2$ (coalescence regime), indicating the emergence of a quasi-self-similar bubble size distribution. The probability density function of the bubble diameter exhibits a dual power-law tail with exponents $-10/3$ and $-3/2$ near the duct inlet. However, after a few hydraulic diameters, a single $-3/2$ power-law scaling emerges, indicating a regime of pure coalescence in which all bubbles are smaller than the Hinze scale. The cumulative distribution plotted against $d/d_{32}$ shows that the slope decreases and the distribution width increases with both axial position and void fraction $(\phi )$. Although classical Hinze scaling gives $d_{\textit{H}} \propto \mathcal{L}^{0.9}$, our theory for $d_{32}$ and $d_{99.8}$ (99.8th percentile bubble diameter) in a pure-coalescence regime predicts the slower law $\propto \mathcal{L}^{0.5}$, which our experimental results confirm – indicating negligible breakup and sub-Hinze growth. Concurrently, in contrast to current models, transient $\phi$ profiles evolve from nearly uniform to sharply core-peaked Gaussian distributions in the developing regime, with increasing centreline values and decreasing near-wall values, due to lift-force reversal. These results provide the first spatially resolved characterisation of coalescence-dominated bubbly flows at high Re, advancing the design of industrial systems as in nuclear cooling and multiphase forming processes (e.g. paper manufacturing, chemical reactors).

In this work, we revisit the Generalised Navier Boundary Condition (GNBC) introduced by Qian et al. in the sharp interface volume-of-fluid context. We replace the singular uncompensated Young stress by a smooth function with a characteristic width $\varepsilon \gt 0$ that is understood as a physical parameter of the model. Therefore, we call the model the ‘contact region GNBC’ (CR-GNBC). We show that the model is consistent with the fundamental kinematics of the contact angle transport described by Fricke, Köhne and Bothe. We implement the model in the geometrical volume-of-fluid solver Basilisk using a ‘free angle’ approach. This means that the dynamic contact angle is not prescribed, but reconstructed from the interface geometry and subsequently applied as an input parameter to compute the uncompensated Young stress. We couple this approach to the two-phase Navier–Stokes solver and study the withdrawing tape problem with a receding contact line. It is shown that the model allows for grid-independent solutions and leads to a full regularisation of the singularity at the moving contact line, which is in accordance with the thin film equation subject to this boundary condition. In particular, it is shown that the curvature at the moving contact line is finite and mesh converging. As predicted by the fundamental kinematics, the parallel shear stress component vanishes at the moving contact line for quasi-stationary states (i.e. for $\dot \theta _d=0$), and the dynamic contact angle is determined by a balance between the uncompensated Young stress and an effective contact line friction. Furthermore, a nonlinear generalisation of the model is proposed, which aims at reproducing the molecular kinetic theory of Blake and Haynes for quasi-stationary states.

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.