This empirical study aims to characterize the dynamical behaviour of sprays using time-series analysis of the size–velocity data acquired using a phase Doppler particle analyser. The prime motivation of this analysis is to capture the spatio-temporal correlations using time-series modelling paradigms that provide valuable new insights into spray dynamics. As a first step, we study long-held assumptions, especially on stationarity and time unsteadiness. We show that air-blast sprays have increased drop size as well as velocity ordering near the edge of the spray. Analysis of the inter-particle time of the droplets shows non-Poisson behaviour where droplets that are closely spaced in time are also closely spaced in the size and velocity coordinates. Temporal auto-correlation and partial auto-correlation calculations reveal the presence of inherent correlated features in the spray. This correlation is stronger and short lived in an air-blast spray and weaker but more persistent in a pressure swirl spray. These correlations render the probability density function (p.d.f.) estimate obtained from standard methods inaccurate; therefore, we propose a technically correct way of estimating the p.d.f. using a suitable downsampling and averaging method. Statistical analysis of residuals (from appropriate autoregressive integrated moving average time-series models) uncovers an interesting feature of spray data pertaining to heteroskedasticity (stochastically changing variance) of the diameter series. In order to account for heteroskedasticity, appropriate generalized autoregressive conditional heteroskedasticity models are developed. Finally, we present a utilitarian view of these results as an empirically consistent boundary condition implementation tool for computational fluid dynamics (CFD).

Hydrodynamic approaches that treat granular materials as a continuum via the homogenization of discrete flow properties have become viable options for efficient predictions of bulk flow behaviours. However, simplified boundary conditions in computational fluid dynamics are often adopted, which have difficulty in describing the complex stick–slip phenomenon at the boundaries. This paper extends the lattice Boltzmann method for granular flow simulations by incorporating a novel frictional boundary condition. The wall slip velocity is first calculated based on the shear rate limited by the Coulomb friction, followed by the reconstruction of unknown density distribution functions through a modified bounce-back scheme. Validation is performed against a unique plane Couette flow configuration, and the analytical solutions for the flow velocity profile and the wall slip velocity, as functions of the friction coefficient, are reproduced by the numerical model. The transition between no-slip and partial-slip regimes is captured well, but the convergence rate drops from second order to first order when slip occurs. The rheological parameters and the basal friction coefficient are calibrated further against the discrete element simulation of a square granular column collapsing over a horizontal bottom plane. It is found that the calibrated continuum model can predict other granular column collapses with different initial aspect ratios and slope inclination angles, including the basal slip and the complex internal flow structures, without any further adjustments to the model parameters. This highlights the generalization ability of the numerical model, which has a wide range of application in granular flow predictions and controls.

A cylinder immersed in a current and free to translate along a circular arc is considered to investigate the impact of path curvature on the flow-induced vibrations (FIV) occurring without structural restoring force. Path curvature magnitude ($\kappa$, inverse of path radius non-dimensionalized by the body diameter $D$) is varied from $0$ (transverse rectilinear path) to $20$, over a wide range of values of the structure to displaced fluid mass ratio, $m^\star \in [0.05,10]$. The exploration is carried out numerically at subcritical and postcritical values of the Reynolds number ($Re$, based on $D$ and the inflow velocity), i.e. below and above the critical value $47$ for the onset of flow unsteadiness when the body is fixed, up to $100$. Path curvature triggers a desynchronized regime of the flow–body system in addition to the synchronized regime typical of vortex-induced vibrations, and alters the composition of fluid forcing. The most prominent effect uncovered here is, however, a global enhancement of FIV, with three principal results: (i) vibrations and flow unsteadiness are found to arise at lower subcritical $Re$ along a curved path, down to $19.5$ versus $31$ for $\kappa =0$; (ii) the $m^\star$ range where substantial responses develop is considerably extended and encompasses the entire interval under study, which contrasts with the narrow band of low $m^\star$ identified for $\kappa =0$; (iii) the vibrations are amplified, $+45\,\%$ relative to the peak amplitude measured along a rectilinear path at $Re =100$.