Numerical simulations are conducted to investigate particle suspension and deposition within turbidity currents. Utilizing Lagrangian particle tracking and a discrete element model, our numerical approach enables a detailed examination of autosuspension, deposition and bulk behaviours of turbidity current. We specifically focus on flow regimes where particle settling and buoyancy-induced hydrodynamics play equally important roles. Our discussion is divided into three parts. Firstly, we examine the main body of the current formed by suspended particles, revealing a temporal evolution consisting of initial slumping, propagation and dissipation stages. Our particle calculation allows for the tracking of autosuspended particles, enabling a deeper understanding of the connection between autosuspension and current propagation through energy budget analysis. In the second part, we delve into particle deposition, highlighting transverse and longitudinal variations. Transverse variations arise from lobe-and-cleft (LC) flow features, while longitudinal variations result from vortex detachment, particularly notable with large-sized particles. We observe that as particle size increases, leading to a particle Stokes number greater than 0.1, rapid particle settling suppresses the LC flow structure, resulting in wider lobes at the deposition height. Lastly, we propose a new scaling law for the propagation speed and current length. Our simulation results demonstrate close agreement with this new scaling law, providing valuable insights into turbidity current dynamics.

We investigate the drag reduction effect of the streamwise travelling wave-like wall deformation in a high-Reynolds-number turbulent channel flow by large-eddy simulation (LES). First, we assess the validity of subgrid-scale models in uncontrolled and controlled flows. For friction Reynolds numbers $Re_\tau = 360$ and $720$, the Smagorinsky and wall-adapting local eddy-viscosity (WALE) models with a damping function can reproduce well the mean velocity profile obtained by direct numerical simulation (DNS) in both the uncontrolled and controlled flows, leading to a small difference in drag reduction rate between LES and DNS. The LES with finer grid resolution can reproduce well the key structures observed in the DNS of the controlled flow. These results show that the high-fidelity LES is valid for appropriately predicting the drag reduction effect. In addition, a small computational domain is sufficient for reproducing the turbulence statistics, key structures and drag reduction rate obtained by DNS. Subsequently, to investigate the trend of drag reduction rate at higher Reynolds numbers, we utilize the WALE model with the damping function to investigate the control effect at higher Reynolds numbers up to $Re_\tau = 3240$. According to the analyses of turbulence statistics and instantaneous flow fields, the drag reduction at higher Reynolds numbers occurs basically through the same mechanism as that at lower Reynolds numbers. In addition, the drag reduction rate obtained by the present LES approaches that predicted using the semi-empirical formula (Nabae et al., Intl J. Heat Fluid Flow, vol. 82, 2020, 108550) as the friction Reynolds number increases, which supports the high predictability of the semi-empirical formula at significantly high Reynolds numbers.

The collapse of an initially spherical cavitation bubble near a free surface leads to the formation of two jets: a downward jet into the liquid, and an upward jet penetrating the free surface. In this study, we examine the surprising interaction of a bubble trapped in a stable cavitating vortex ring approaching a free surface. As a result, a single fast and tall liquid jet forms. We find that this jet is observed only above critical Froude numbers ($Fr$) and Weber numbers ($We$) when ${Fr}^2 (1.6-2.73/{We}) > 1$, illustrating the importance of inertia, gravity and surface tension in accelerating this novel jet and thereby reaching heights several hundred times the radius of the vortex ring. Our experimental results are supported by numerical simulations, revealing that the underlying mechanism driving the vortex ring acceleration is the disruption of the equilibrium of high-pressure regions at the front and rear of the vortex ring caused by the free surface. Quantitative analysis based on the energy relationships elucidates that the velocity ratio between the maximum velocity of the free-surface jet and the translational velocity of the vortex ring is relatively stable yet is attenuated by surface tension when the jet is mild.

Ekman pumping is a phenomenon induced by no-slip boundary conditions in rotating fluids. In the context of Rayleigh–Bénard convection, Ekman pumping causes a significant change in the linear stability of the system compared with when it is not present (that is, stress-free). Motivated by numerical solutions to the marginal stability problem of the incompressible Navier–Stokes equation (iNSE) system, we seek analytical asymptotic solutions which describe the departure of the no-slip solution from the stress-free one. The substitution of normal modes into a reduced asymptotic model yields a linear system for which we explore analytical solutions for various scalings of wavenumber. We find very good agreement between the analytical asymptotic solutions and the numerical solutions to the iNSE linear stability problem with no-slip boundary conditions.