The relative velocities and positions of monodisperse high-inertia particle pairs in isotropic turbulence are studied using direct numerical simulations (DNS), as well as Langevin simulations (LS) based on a probability density function (PDF) kinetic model for pair relative motion. In a prior study (Rani et al. , J. Fluid Mech. , vol. 756, 2014, pp. 870–902), the authors developed a stochastic theory that involved deriving closures in the limit of high Stokes number for the diffusivity tensor in the PDF equation for monodisperse particle pairs. The diffusivity contained the time integral of the Eulerian two-time correlation of fluid relative velocities seen by pairs that are nearly stationary. The two-time correlation was analytically resolved through the approximation that the temporal change in the fluid relative velocities seen by a pair occurs principally due to the advection of smaller eddies past the pair by large-scale eddies. Accordingly, two diffusivity expressions were obtained based on whether the pair centre of mass remained fixed during flow time scales, or moved in response to integral-scale eddies. In the current study, a quantitative analysis of the (Rani et al. 2014) stochastic theory is performed through a comparison of the pair statistics obtained using LS with those from DNS. LS consist of evolving the Langevin equations for pair separation and relative velocity, which is statistically equivalent to solving the classical Fokker–Planck form of the pair PDF equation. Langevin simulations of particle-pair dispersion were performed using three closure forms of the diffusivity – i.e. the one containing the time integral of the Eulerian two-time correlation of the seen fluid relative velocities and the two analytical diffusivity expressions. In the first closure form, the two-time correlation was computed using DNS of forced isotropic turbulence laden with stationary particles. The two analytical closure forms have the advantage that they can be evaluated using a model for the turbulence energy spectrum that closely matched the DNS spectrum. The three diffusivities are analysed to quantify the effects of the approximations made in deriving them. Pair relative-motion statistics obtained from the three sets of Langevin simulations are compared with the results from the DNS of (moving) particle-laden forced isotropic turbulence for $St_{\unicode[STIX]{x1D702}}=10,20,40,80$ and $Re_{\unicode[STIX]{x1D706}}=76,131$ . Here, $St_{\unicode[STIX]{x1D702}}$ is the particle Stokes number based on the Kolmogorov time scale and $Re_{\unicode[STIX]{x1D706}}$ is the Taylor micro-scale Reynolds number. Statistics such as the radial distribution function (RDF), the variance and kurtosis of particle-pair relative velocities and the particle collision kernel were computed using both Langevin and DNS runs, and compared. The RDFs from the stochastic runs were in good agreement with those from the DNS. Also computed were the PDFs $\unicode[STIX]{x1D6FA}(U|r)$ and $\unicode[STIX]{x1D6FA}(U_{r}|r)$ of relative velocity $U$ and of the radial component of relative velocity $U_{r}$ respectively, both PDFs conditioned on separation $r$ . The first closure form, involving the Eulerian two-time correlation of fluid relative velocities, showed the best agreement with the DNS results for the PDFs.
more »
« less
Particle pair statistics of inertial particles at small separation using stereoscopic particle tracking
Particle pair statistics of inertial particles having average Stokes numbers of 2.1 and 14 are measured in isotropic turbulence at a Reynolds number of Reλ = 240. The radial distribution function (RDF) and mean relative approach velocity are obtained at small separation distances using 2-frame stereoscopic particle tracking velocimetry (stereo-PTV). At small separation distance, the RDF varies by an order of magnitude in the range of Stokes numbers investigated. However, the mean relative approach velocity is found to have a weak dependence on Stokes number. The results are shown to have high accuracy when compared to analogous mono-PTV datasets, and can be used to provide a more reliable estimate of the inter-particle collision rate. The main limitation of the measurement is observed at separation distances less than the laser sheet thickness, where the technique tended to underestimate the mean relative approach velocity.
more »
« less
- Award ID(s):
- 1756068
- NSF-PAR ID:
- 10302024
- Date Published:
- Journal Name:
- 14th International Symposium on Particle Image Velocimetry
- Volume:
- 1
- Issue:
- 1
- ISSN:
- 2769-7576
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Particle ingestion into turbine engines is a widespread problem that can cause significant degradation in engine service life. One primary damage mechanism is deposition of particulate matter in internal cooling passages. Musgrove et al. proposed a compact particle separator that could be installed between the combustor bypass exit and turbine vane cooling passage inlet. The design had small pressure losses but provided limited particle separation, and its performance has proved difficult to replicate in subsequent experiments. Borup et al recently developed a Magnetic Resonance Imaging (MRI) based technique for making full-field, 3D measurements of the mean particle concentration distribution in complex flows. A particle separator based on the Musgrove et all design was fabricated out of plastic using 3D printing. The primary difference from earlier designs was the addition of a drain from the collector, through which 3% of the total flow was extracted. The separator efficiency was measured at two Reynolds numbers, using water as the working fluid and 33 micron titanium microspheres to represent dust particles. Particle Stokes number was shown to play the dominant role in determining efficiency across studies. MRI was used to obtain the 3D volume fraction and 3-component velocity fields. The velocity data showed that flow was poorly distributed between the separator louvers, while the collector flow followed the optimal pattern for particle retention. The particle distribution data revealed that strong swirling flow in the collector centrifuged particles toward the outer wall of the collector and intro a partitioned region of quiescent flow, where they proceeded to exit the collector.more » « less
-
The transport and deposition of firebrand particles is an important fire spread mechanism in wildland fires. These particles can be transported by wind over large distances and can ignite secondary fires upon landing. The transport of firebrands by wind is a complex, multiscale process that is largely controlled by interactions between the firebrand particles and the atmospheric wind. To account for the complex temporal evolution of atmospheric turbulence over large scales, the use of large-eddy simulation (LES) techniques is necessary. However, filtering of subgrid-scale (SGS) turbulence in LES hinders the accuracy of particle transport models. In this work, we employ a Lagrangian SGS model in an LES framework to investigate the effects of small-scale turbulence on the transport of mass- and size-changing firebrand particles. The impact of SGS turbulence was analyzed by comparing landing and trajectory statistics for firebrand and regular (fixed size and mass) particles under different Stokes numbers. It was found that the presence of SGS turbulence modifies the particle transport behavior, which is characterized by smaller spanwise dispersions but larger travel distances along the streamwise direction compared with particles under no SGS turbulence. As expected, the enhanced velocity field produced by the SGS model has larger influence on the statistics of firebrand particles compared with regular particles due to the time-evolving reduction in particle mass and size induced by pyrolysis.more » « less
-
more » « less
Abstract Explosive eruptions expel volcanic gases and particles at high pressures and velocities. Within this multiphase fluid, small ash particles affect the flow dynamics, impacting mixing, entrainment, turbulence, and aggregation. To examine the role of turbulent particle behavior, we conducted an analogue experiment using a particle‐laden jet. We used compressed air as the carrier fluid, considering turbulent conditions at Reynolds numbers from approximately 5,000 to 20,000. Two different particles were examined: 14‐μm diameter solid nickel spheres and 13‐μm diameter hollow glass spheres. These resulted in Stokes numbers between 1 and 35 based on the convective scale. The particle mass percentage in the mixture is varied from 0.3% to more than 20%. Based on a 1‐D volcanic plume model, these Stokes numbers and mass loadings corresponded to millimeter‐scale particle diameters at heights of 4–8 km above the vent during large, sustained eruptions. Through particle image velocimetry, we measured the mean flow behavior and the turbulence statistics in the near‐exit region, primarily focusing on the dispersed phase. We show that the flow behavior is dominated by the particle inertia, with high Stokes numbers reducing the entrainment by more than 40%. When applied to volcanic plumes, these results suggest that high‐density particles can greatly increase the probability of column collapse.
-
This paper presents a theory to obtain the force per unit length acting on a slender filament with a non-circular cross-section moving in a fluid at low Reynolds number. Using a regular perturbation of the inner solution, we show that the force per unit length has $O(1/\ln (2A))+O(\unicode[STIX]{x1D6FC}/\ln ^{2}(2A))$ contributions driven by the relative motion of the particle and the local fluid velocity and an $O(\unicode[STIX]{x1D6FC}/(\ln (2A)A))$ contribution driven by the gradient in the imposed fluid velocity. Here, the aspect ratio ( $A=l/a_{0}$ ) is defined as the ratio of the particle size ( $l$ ) to the cross-sectional dimension ( $a_{0}$ ) and $\unicode[STIX]{x1D6FC}$ is the amplitude of the non-circular perturbation. Using thought experiments, we show that two-lobed and three-lobed cross-sections affect the response to relative motion and velocity gradients, respectively. A two-dimensional Stokes flow calculation is used to extend the perturbation analysis to cross-sections that deviate significantly from a circle (i.e. $\unicode[STIX]{x1D6FC}\sim O(1)$ ). We demonstrate the ability of our method to accurately compute the resistance to translation and rotation of a slender triaxial ellipsoid. Furthermore, we illustrate novel dynamics of straight rods in a simple shear flow that translate and rotate quasi-periodically if they have two-lobed cross-section, and rotate chaotically and translate diffusively if they have a combination of two- and three-lobed cross-sections. Finally, we show the remarkable ability of our theory to accurately predict the motion of rings, retaining great accuracy for moderate aspect ratios ( ${\sim}10$ ) and cross-sections that deviate significantly from a circle, thereby making our theory a computationally inexpensive alternative to other Stokes flow solvers.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.18409/ispiv.v1i1.27
