Solid atmospheric particles, such as ice crystals, pollen, dust, ash and microplastics, strongly influence Earth’s climate, ecosystems and air quality. Previous studies have typically relied on analytical models valid only for very small particles or experiments in liquids, where the particle-to-fluid density ratio $R$ is much lower than values encountered in the atmosphere. We combine a novel experimental set-up with particle-resolved direct numerical simulations to study the settling of sub-millimetre ellipsoids in still air. Particle shapes span elongation and flatness values $ 0.2 \leqslant {\textit{EL}}, {\textit{FL}} \leqslant 1.0$ at a density ratio $ R = 1000$ and particle Reynolds numbers $ 2.1 \lt {\textit{Re}}_{\!p} \lt 4.5$, a regime well below the onset of wake-induced instabilities. Nonetheless, we observe unexpectedly rich dynamics: all non-spherical particles exhibit damped oscillatory motion, and some triaxial ellipsoids follow fully three-dimensional, non-planar trajectories due to rotation about all three axes. Simulations at lower density ratios ($ R = 10, 100$) confirm that these behaviours are driven by strong lateral forces happening only at $R=1000$. Key settling characteristics exhibit nonlinear and non-trivial dependencies on shape. In the two-dimensional phase space of elongation and flatness, settling velocity is symmetric about the principal diagonal ($ {\textit{EL}} = {\textit{FL}}$), while oscillation frequency and damping rate show symmetry about the anti-diagonal. Flatness strongly influences pressure drag, while elongation governs lateral drift and swept volume, which can reach up to ten times the particle diameter and four times the volume-equivalent sphere, respectively.

We investigate the energetics of mixing induced by a continuously supplied dense current (density $\rho _0$) propagating beneath a lighter ambient fluid (density $\rho _a$) along a horizontal rigid boundary within a rectangular domain. The flow fields are computed using direct numerical simulations (DNS) performed with the Nek5000 spectral element solver. Mixing is quantified through the temporal evolution of the background potential energy, which exhibits a linear increase over time. This linear trend enables the definition of a dimensionless mixing parameter $\gamma$, representing the rate of background potential energy growth. The value of $\gamma$ depends on the initial density contrast for a fixed volumetric discharge at the source, characterised by the dimensionless source Froude number. The results reveal a non-monotonic dependence of $\gamma$ on the source Froude number, highlighting a complex interaction between flow forcing and mixing efficiency. We find that, under the assumption of uniform mixing along the current’s length, a fraction $\gamma /2$ of the total supplied energy is invested in mixing along a horizontal distance equal to the height of the inlet.

The present article investigates the stability of Rayleigh–Bénard convection in a composite system consisting of a horizontal fluid layer overlying a fluid-saturated Darcy porous layer subjected to a time-periodic temperature distribution. The bottom surface is heated periodically with time, whereas a Biot number-dependent thermal boundary condition represents the heat transfer at the upper surface. The Beavers–Joseph–Saffman–Jones condition describes the ‘slip’ at the interface of the domains, and the Lions interface condition governs the normal force balance, incorporating a dynamic pressure term. The Chebyshev tau method and Fourier analysis are utilised to obtain linear instability bounds, which are compared with strong global and asymptotic limits derived from the nonlinear analysis using the energy method. Four deliberately chosen configurations of superposed fluid- and porous-layer systems are investigated. Two configurations validate the analysis through the limiting cases of the classical Darcy–Bénard and Rayleigh–Bénard systems obtained by setting the fluid-to-porous depth ratio $(\hat {d})$ to zero and infinity, respectively. The other two configurations involve layers with equal depths $(\hat {d} =1)$ and a shallow fluid layer overlying a porous layer $(\hat {d} \sim 0.1)$. For these cases, modulation substantially influences the onset of convection. In the last case, the linear theory points out that modulation parameters can control the dominant convective mode (fluid/porous). Furthermore, unlike the previously reported studies, the nonlinear stability bounds are found to be significantly lower than the linear instability bounds, indicating the possibility of subcritical instabilities in the presence of modulation. The region of subcritical instabilities increases with modulation amplitude.

Supersonic jets impinging on a ground plane produce a highly unsteady jet shear layer, often resulting in extremely high noise level. The widely accepted mechanism for this jet resonance involves a feedback loop consisting of downstream-travelling coherent structures and upstream-propagating acoustic waves. Despite the importance of coherent structures, often referred to as disturbances, that travel downstream, a comprehensive discussion on the disturbance convection velocity has been limited due to the challenges posed by non-intrusive measurement requirements. To determine the convection velocity of disturbances in the jet shear layer, a high-speed schlieren flow visualisation is carried out, and phase-averaged wave diagrams are constructed from the image sets. The experiments are conducted using a Mach 1.5 jet under various nozzle pressure ratios and across a range of impingement distances. A parametric analysis is performed to examine the influence of nozzle pressure ratio on the convection velocity and phase lead/lag at specific impingement distances. The results reveal that impingement tonal frequency is nearly independent of the disturbance convection velocity, except in cases of staging behaviour. They also demonstrated that slower downstream convection velocity of the disturbance corresponds to larger coherent structures, resulting in increased noise levels. Based on the observation of acoustic standing waves, an acoustic speed-based frequency model has been proposed. With the help of the allowable frequency range calculated from the vortex-sheet model, this model can provide a good approximation for the majority of axisymmetric impingement tonal frequencies.

One-degree-of-freedom flow-induced vibration (FIV) and energy harvesting through FIV of an elastically mounted circular cylinder with mechanically coupled rotation were investigated numerically for low Reynolds number 100, mass ratio 8 and a wide range of reduced velocities. The aims of this study are to investigate the effect of the flow direction angle $\beta$ on the vibration and energy harvesting through FIV. Two types of lock-in are found: vortex-induced vibration (VIV) and galloping. The response amplitude increases with the increase of $\beta$ in both regimes. Both VIV response and galloping regimes are found for $\beta$ = 45° to $\beta$ = 90°. For $\beta$ = −90° to $\beta$ = 0°, only VIV response regimes are found. The fluid force and fluid torque play different roles in exciting/damping the vibration. In the high-amplitude gallop regime, the fluid force excites the vibration, and the torque damps the vibration. Energy harvesting at flow direction angle 90° is investigated as this flow direction has the maximum galloping amplitude. The energy harvesting is achieved by a linear electric damping coefficient in the numerical model. The maximum harvestable power in the galloping regime is significantly greater than that in the VIV regime, and it increases with the increase of the reduced velocity. When the reduced velocity is 20, the harvested power is over 20 times that in the VIV regime, and can further increase if reduced velocity further increases. The maximum efficiency over all simulated parameters is 0.424, occurring when the reduced velocity is 20, and electric damping factor is 0.04.