Inertial particle relative velocity statistics in homogeneous isotropic turbulence

2012 ◽  
Vol 696 ◽  
pp. 45-66 ◽  
Author(s):  
Juan P. L. C. Salazar ◽  
Lance R. Collins
Keyword(s):  
Structure Function ◽  
Relative Velocity ◽  
Velocity Structure ◽  
Structure Functions ◽  
Inertial Subrange ◽  
Inertial Particles ◽  
Inertial Particle ◽  
Scaling Exponents ◽  
Velocity Statistics ◽  
Dissipation Range

AbstractIn the present study, we investigate the scaling of relative velocity structure functions, of order two and higher, for inertial particles, both in the dissipation range and the inertial subrange using direct numerical simulations (DNS). Within the inertial subrange our findings show that contrary to the well-known attenuation in the tails of the one-point acceleration probability density function (p.d.f.) with increasing inertia (Bec et al., J. Fluid Mech., vol. 550, 2006, pp. 349–358), the opposite occurs with the velocity structure function at sufficiently large Stokes numbers. We observe reduced scaling exponents for the structure function when compared to those of the fluid, and correspondingly broader p.d.f.s, similar to what occurs with a passive scalar. DNS allows us to isolate the two effects of inertia, namely biased sampling of the velocity field, a result of preferential concentration, and filtering, i.e. the tendency for the inertial particle velocity to attenuate the velocity fluctuations in the fluid. By isolating these effects, we show that sampling is playing the dominant role for low-order moments of the structure function, whereas filtering accounts for most of the scaling behaviour observed with the higher-order structure functions in the inertial subrange. In the dissipation range, we see evidence of so-called ‘crossing trajectories’, the ‘sling effect’ or ‘caustics’, and find good agreement with the theory put forth by Wilkinson et al. (Phys. Rev. Lett., vol. 97, 2006, 048501) and Falkovich & Pumir (J. Atmos. Sci., vol. 64, 2007, 4497) for Stokes numbers greater than 0.5. We also look at the scaling exponents within the context of the model proposed by Bec et al. (J. Fluid Mech., vol. 646, 2010, pp. 527–536). Another interesting finding is that inertial particles at low Stokes numbers sample regions of higher kinetic energy than the fluid particle field, the converse occurring at high Stokes numbers. The trend at low Stokes numbers is predicted by the theory of Chun et al. (J. Fluid Mech., vol. 536, 2005, 219–251). This work is relevant to modelling the particle collision rate (Sundaram & Collins, J. Fluid Mech., vol. 335, 1997, pp. 75–109), and highlights the interesting array of phenomena induced by inertia.

Investigation of sub-Kolmogorov inertial particle pair dynamics in turbulence using novel satellite particle simulations

2013 ◽  
Vol 720 ◽  
pp. 192-211 ◽  
Author(s):  
Baidurja Ray ◽  
Lance R. Collins
Keyword(s):  
Structure Function ◽  
Power Law ◽  
Particle Motion ◽  
Relative Velocity ◽  
Velocity Structure ◽  
Inertial Particles ◽  
Particle Clustering ◽  
Particle Pair ◽  
Velocity Statistics ◽  
Particle Simulations

AbstractClustering (or preferential concentration) of weakly inertial particles suspended in a homogeneous isotropic turbulent flow is driven primarily by the smallest eddies at the so-called Kolmogorov scale. In particle-laden large-eddy simulations (LES), these small scales are not resolved by the grid and hence their effect on both the resolved flow scales and the particle motion have to be modelled. In order to predict clustering in a particle-laden LES, it is crucial that the subgrid model for the particles captures the mechanism by which the subgrid scales affect the particle motion (Ray & Collins, J. Fluid Mech., vol. 680, 2011, pp. 488–510). In this paper, we describe novel satellite particle simulations (SPS), in which we study the clustering and relative velocity statistics of inertial particles at separation distances well below the Kolmogorov length scale. SPS is designed to isolate pairwise interactions of particles, and is therefore well suited for developing two-particle models. We show that the power-law dependence of the radial distribution function (RDF), a statistical measure of clustering, is predicted by the SPS in excellent agreement with direct numerical simulations (DNS) for Stokes numbers up to 3, implying that no explicit information from the inertial range is required to accurately describe particle clustering. This result further explains our successful prediction of the RDF power using the drift-diffusion model of Chun et al. (J. Fluid Mech., vol. 536, 2005, pp. 219–251) for $\mathit{St}\leq 0. 4$. We also consider the second-order longitudinal relative velocity structure function for the particles; we show that the SPS is able to capture its power-law exponent for $\mathit{St}\leq 0. 5$ and attribute the disagreement at larger $\mathit{St}$ to the effect of the larger scales of motion not captured by the SPS. Further, the SPS is able to capture the ‘caustic activation’ of the structure function at zero separation and predict the critical $\mathit{St}$ and rate of activation in agreement with the DNS (Salazar & Collins, J. Fluid. Mech., vol. 696, 2012, pp. 45–66). We show comparisons between filtered DNS and equivalently filtered SPS, and the findings are similar to the unfiltered case. Overall, SPS is an efficient and accurate computational tool for investigating particle pair dynamics at small separations, as well as an interesting platform for developing LES subgrid models designed to accurately reproduce particle clustering.

scholarly journals The effect of Reynolds number on inertial particle dynamics in isotropic turbulence. Part 1. Simulations without gravitational effects

2016 ◽  
Vol 796 ◽  
pp. 617-658 ◽  
Author(s):  
Peter J. Ireland ◽  
Andrew D. Bragg ◽  
Lance R. Collins
Keyword(s):  
Reynolds Number ◽  
Relative Velocity ◽  
Isotropic Turbulence ◽  
Reynolds Numbers ◽  
Inertial Particles ◽  
Inertial Particle ◽  
Physical Mechanisms ◽  
Preferential Sampling ◽  
Velocity Statistics ◽  
Particle Statistics

In this study, we analyse the statistics of both individual inertial particles and inertial particle pairs in direct numerical simulations of homogeneous isotropic turbulence in the absence of gravity. The effect of the Taylor microscale Reynolds number, $R_{{\it\lambda}}$, on the particle statistics is examined over the largest range to date (from $R_{{\it\lambda}}=88$ to 597), at small, intermediate and large Kolmogorov-scale Stokes numbers $St$. We first explore the effect of preferential sampling on the single-particle statistics and find that low-$St$ inertial particles are ejected from both vortex tubes and vortex sheets (the latter becoming increasingly prevalent at higher Reynolds numbers) and preferentially accumulate in regions of irrotational dissipation. We use this understanding of preferential sampling to provide a physical explanation for many of the trends in the particle velocity gradients, kinetic energies and accelerations at low $St$, which are well represented by the model of Chun et al. (J. Fluid Mech., vol. 536, 2005, pp. 219–251). As $St$ increases, inertial filtering effects become more important, causing the particle kinetic energies and accelerations to decrease. The effect of inertial filtering on the particle kinetic energies and accelerations diminishes with increasing Reynolds number and is well captured by the models of Abrahamson (Chem. Engng Sci., vol. 30, 1975, pp. 1371–1379) and Zaichik & Alipchenkov (Intl J. Multiphase Flow, vol. 34 (9), 2008, pp. 865–868), respectively. We then consider particle-pair statistics, and focus our attention on the relative velocities and radial distribution functions (RDFs) of the particles, with the aim of understanding the underlying physical mechanisms contributing to particle collisions. The relative velocity statistics indicate that preferential sampling effects are important for $St\lesssim 0.1$ and that path-history/non-local effects become increasingly important for $St\gtrsim 0.2$. While higher-order relative velocity statistics are influenced by the increased intermittency of the turbulence at high Reynolds numbers, the lower-order relative velocity statistics are only weakly sensitive to changes in Reynolds number at low $St$. The Reynolds-number trends in these quantities at intermediate and large $St$ are explained based on the influence of the available flow scales on the path-history and inertial filtering effects. We find that the RDFs peak near $St$ of order unity, that they exhibit power-law scaling for low and intermediate $St$ and that they are largely independent of Reynolds number for low and intermediate $St$. We use the model of Zaichik & Alipchenkov (New J. Phys., vol. 11, 2009, 103018) to explain the physical mechanisms responsible for these trends, and find that this model is able to capture the quantitative behaviour of the RDFs extremely well when direct numerical simulation data for the structure functions are specified, in agreement with Bragg & Collins (New J. Phys., vol. 16, 2014a, 055013). We also observe that at large $St$, changes in the RDF are related to changes in the scaling exponents of the relative velocity variances. The particle collision kernel closely matches that computed by Rosa et al. (New J. Phys., vol. 15, 2013, 045032) and is found to be largely insensitive to the flow Reynolds number. This suggests that relatively low-Reynolds-number simulations may be able to capture much of the relevant physics of droplet collisions and growth in the adiabatic cores of atmospheric clouds.

scholarly journals On the universality of anomalous scaling exponents of structure functions in turbulent flows

2018 ◽  
Vol 837 ◽  
pp. 657-669 ◽  
Author(s):  
E.-W. Saw ◽  
P. Debue ◽  
D. Kuzzay ◽  
F. Daviaud ◽  
B. Dubrulle
Keyword(s):  
Turbulent Flows ◽  
Large Scale ◽  
Velocity Structure ◽  
Swirling Flow ◽  
Spatial Scales ◽  
Structure Functions ◽  
Inertial Subrange ◽  
Scaling Exponents ◽  
Von Karman ◽  
Von Kármán

All previous experiments in open turbulent flows (e.g. downstream of grids, jets and the atmospheric boundary layer) have produced quantitatively consistent values for the scaling exponents of velocity structure functions (Anselmet et al., J. Fluid Mech., vol. 140, 1984, pp. 63–89; Stolovitzky et al., Phys. Rev. E, vol. 48 (5), 1993, R3217; Arneodo et al., Europhys. Lett., vol. 34 (6), 1996, p. 411). The only measurement of scaling exponents at high order (${>}6$) in closed turbulent flow (von Kármán swirling flow) using Taylor’s frozen flow hypothesis, however, produced scaling exponents that are significantly smaller, suggesting that the universality of these exponents is broken with respect to change of large scale geometry of the flow. Here, we report measurements of longitudinal structure functions of velocity in a von Kármán set-up without the use of the Taylor hypothesis. The measurements are made using stereo particle image velocimetry at four different ranges of spatial scales, in order to observe a combined inertial subrange spanning approximately one and a half orders of magnitude. We found scaling exponents (up to ninth order) that are consistent with values from open turbulent flows, suggesting that they might be in fact universal.

Scale-dependent alignment, tumbling and stretching of slender rods in isotropic turbulence

2018 ◽  
Vol 860 ◽  
pp. 465-486 ◽  
Author(s):  
Nimish Pujara ◽  
Greg A. Voth ◽  
Evan A. Variano
Keyword(s):  
Isotropic Turbulence ◽  
Fluid Velocity ◽  
Inertial Range ◽  
Small Scale ◽  
Strain Rate Tensor ◽  
Scaling Exponents ◽  
Slender Body Theory ◽  
Velocity Statistics ◽  
Rigid Rods ◽  
Dissipation Range

We examine the dynamics of slender, rigid rods in direct numerical simulation of isotropic turbulence. The focus is on the statistics of three quantities and how they vary as rod length increases from the dissipation range to the inertial range. These quantities are (i) the steady-state rod alignment with respect to the perceived velocity gradients in the surrounding flow, (ii) the rate of rod reorientation (tumbling) and (iii) the rate at which the rod end points move apart (stretching). Under the approximations of slender-body theory, the rod inertia is neglected and rods are modelled as passive particles in the flow that do not affect the fluid velocity field. We find that the average rod alignment changes qualitatively as rod length increases from the dissipation range to the inertial range. While rods in the dissipation range align most strongly with fluid vorticity, rods in the inertial range align most strongly with the most extensional eigenvector of the perceived strain-rate tensor. For rods in the inertial range, we find that the variance of rod stretching and the variance of rod tumbling both scale as $l^{-4/3}$, where $l$ is the rod length. However, when rod dynamics are compared to two-point fluid velocity statistics (structure functions), we see non-monotonic behaviour in the variance of rod tumbling due to the influence of small-scale fluid motions. Additionally, we find that the skewness of rod stretching does not show scale invariance in the inertial range, in contrast to the skewness of longitudinal fluid velocity increments as predicted by Kolmogorov’s $4/5$ law. Finally, we examine the power-law scaling exponents of higher-order moments of rod tumbling and rod stretching for rods with lengths in the inertial range and find that they show anomalous scaling. We compare these scaling exponents to predictions from Kolmogorov’s refined similarity hypotheses.

scholarly journals Solar wind low-frequency magnetohydrodynamic turbulence: extended self-similarity and scaling laws

1996 ◽  
Vol 3 (4) ◽  
pp. 247-261 ◽  
Author(s):  
V. Carbone ◽  
P. Veltri ◽  
R. Bruno
Keyword(s):  
Solar Wind ◽  
Scaling Laws ◽  
Velocity Structure ◽  
Fluid Flows ◽  
Structure Functions ◽  
Point Of View ◽  
Self Similarity ◽  
Extended Self ◽  
Scaling Exponents ◽  
Magnetohydrodynamic Turbulence

Abstract. In this paper we review some of the work done in investigating the scaling properties of Magnetohydrodynamic turbulence, by using velocity fluctuations measurements performed in the interplanetary space plasma by the Helios spacecraft. The set of scaling exponents ξq for the q-th order velocity structure functions, have been determined by using the Extended Self-Similarity hypothesis. We have found that the q-th order velocity structure function, when plotted vs. the 4-th order structure function, displays a range of self-similarity which extends over all the lengths covered by measurements, thus allowing for a very good determination of ξq. Moreover the results seem to show that the scaling exponents are the same regardless the various observation periods considered. The obtained scaling exponents have been compared with the results of some intermittency models for Kraichnan's turbulence, derived in the framework of infinitely divisible fragmentation processes, showing the good agreement between these models and our observations. Finally, on the basis of the actually available data sets, we show that scaling laws in Solar Wind turbulence seem to be different from turbulent scaling laws in the ordinary fluid flows. This is true for high-order velocity structure functions, while low-order velocity structure functions show the same scaling laws. Since our measurements involve length scales which extend over many order of magnitude where dissipation is practically absent, our results show that Solar Wind turbulence can be regarded as a testing bench for the investigation of general scaling behaviour in turbulent flows. In particular our results strongly support the point of view which attributes a key role to the inertial range dynamics in determining the intermittency characteristics in fluid flows, in contrast with the point of view which attributes intermittency to a finite Reynolds number effect.

scholarly journals Analysis of the forward and backward in time pair-separation probability density functions for inertial particles in isotropic turbulence

2017 ◽  
Vol 830 ◽  
pp. 63-92 ◽  
Author(s):  
Andrew D. Bragg
Keyword(s):  
Probability Density ◽  
Isotropic Turbulence ◽  
Probability Density Functions ◽  
Density Functions ◽  
Inertial Particles ◽  
Mean Square ◽  
Inertial Particle ◽  
Particle Pair ◽  
Pair Separation ◽  
Dissipation Range

In this paper we investigate, using theory and direct numerical simulations (DNS), the forward in time (FIT) and backward in time (BIT) probability density functions (PDFs) of the separation of inertial particle pairs in isotropic turbulence. In agreement with our earlier study (Bragg et al., Phys. Fluids, vol. 28, 2016, 013305), where we compared the FIT and BIT mean-square separations, we find that inertial particles separate much faster BIT than FIT, with the strength of the irreversibility depending upon the final/initial separation of the particle pair and their Stokes number $St$. However, we also find that the irreversibility shows up in subtle ways in the behaviour of the full PDF that it does not in the mean-square separation. In the theory, we derive new predictions, including a prediction for the BIT/FIT PDF for $St\geqslant O(1)$, and for final/initial separations in the dissipation regime. The prediction shows how caustics in the particle relative velocities in the dissipation range affect the scaling of the pair-separation PDF, leading to a PDF with an algebraically decaying tail. The predicted functional behaviour of the PDFs is universal, in that it does not depend upon the level of intermittency in the underlying turbulence. We also analyse the pair-separation PDFs for fluid particles at short times, and construct theoretical predictions using the multifractal formalism to describe the fluid relative velocity distributions. The theoretical and numerical results both suggest that the extreme events in the inertial particle-pair dispersion at the small scales are dominated by their non-local interaction with the turbulent velocity field, rather than due to the strong dissipation range intermittency of the turbulence itself. In fact, our theoretical results predict that for final/initial separations in the dissipation range, when $St\gtrsim 1$, the tails of the pair-separation PDFs decay faster as the Taylor Reynolds number $Re_{\unicode[STIX]{x1D706}}$ is increased, the opposite of what would be expected for fluid particles.

scholarly journals Equations relating structure functions of all orders

2001 ◽  
Vol 434 ◽  
pp. 379-388 ◽  
Author(s):  
REGINALD J. HILL
Keyword(s):  
Velocity Structure ◽  
Arbitrary Order ◽  
Structure Functions ◽  
Navier Stokes ◽  
Isotropic Case ◽  
Incompressibility Condition ◽  
Scaling Exponents ◽  
Navier Stokes Equation ◽  
Local Isotropy ◽  
Locally Homogeneous

Exact equations are given that relate velocity structure functions of arbitrary order with other statistics. ‘Exact’ means that no approximations are used except that the Navier–Stokes equation and incompressibility condition are assumed to be accurate. The exact equations are used to determine the structure function equations of all orders for locally homogeneous but anisotropic turbulence as well as for the locally isotropic case. The uses of these equations for investigating the approach to local homogeneity as well as to local isotropy and the balance of the equations and identification of scaling ranges are discussed. The implications for scaling exponents and investigation of intermittency are briefly discussed.

scholarly journals Experimental study of the influence of anisotropy on the inertial scales of turbulence

2012 ◽  
Vol 692 ◽  
pp. 464-481 ◽  
Author(s):  
Kelken Chang ◽  
Gregory P. Bewley ◽  
Eberhard Bodenschatz
Keyword(s):  
Structure Function ◽  
Velocity Structure ◽  
Transverse Velocity ◽  
Scaling Exponent ◽  
Mean Flow ◽  
Scaling Exponents ◽  
Velocity Fluctuations ◽  
Velocity Structure Function ◽  
Central Volume ◽  
Kolmogorov Constant

AbstractWe ask whether the scaling exponents or the Kolmogorov constants depend on the anisotropy of the velocity fluctuations in a turbulent flow with no shear. According to our experiment, the answer is no for the Eulerian second-order transverse velocity structure function. The experiment consisted of 32 loudspeaker-driven jets pointed toward the centre of a spherical chamber. We generated anisotropy by controlling the strengths of the jets. We found that the form of the anisotropy of the velocity fluctuations was the same as that in the strength of the jets. We then varied the anisotropy, as measured by the ratio of axial to radial root-mean-square (r.m.s.) velocity fluctuations, between 0.6 and 2.3. The Reynolds number was approximately constant at around ${R}_{\lambda } = 481$. In a central volume with a radius of 50 mm, the turbulence was approximately homogeneous, axisymmetric, and had no shear and no mean flow. We observed that the scaling exponent of the structure function was $0. 70\pm 0. 03$, independent of the anisotropy and regardless of the direction in which we measured it. The Kolmogorov constant, ${C}_{2} $, was also independent of direction and anisotropy to within the experimental error of 4 %.

Finite Reynolds number effect on the scaling range behaviour of turbulent longitudinal velocity structure functions

2017 ◽  
Vol 820 ◽  
pp. 341-369 ◽  
Author(s):  
S. L. Tang ◽  
R. A. Antonia ◽  
L. Djenidi ◽  
L. Danaila ◽  
Y. Zhou
Keyword(s):  
Reynolds Number ◽  
Structure Function ◽  
Large Scale ◽  
Velocity Structure ◽  
Longitudinal Velocity ◽  
Structure Functions ◽  
Third Order ◽  
Reynolds Number Effect ◽  
Rate Of Increase ◽  
Velocity Structure Function

The effect of large-scale forcing on the second- and third-order longitudinal velocity structure functions, evaluated at the Taylor microscale $r=\unicode[STIX]{x1D706}$, is assessed in various turbulent flows at small to moderate values of the Taylor microscale Reynolds number $R_{\unicode[STIX]{x1D706}}$. It is found that the contribution of the large-scale terms to the scale by scale energy budget differs from flow to flow. For a fixed $R_{\unicode[STIX]{x1D706}}$, this contribution is largest on the centreline of a fully developed channel flow but smallest for stationary forced periodic box turbulence. For decaying-type flows, the contribution lies between the previous two cases. Because of the difference in the large-scale term between flows, the third-order longitudinal velocity structure function at $r=\unicode[STIX]{x1D706}$ differs from flow to flow at small to moderate $R_{\unicode[STIX]{x1D706}}$. The effect on the second-order velocity structure functions appears to be negligible. More importantly, the effect of $R_{\unicode[STIX]{x1D706}}$ on the scaling range exponent of the longitudinal velocity structure function is assessed using measurements of the streamwise velocity fluctuation $u$, with $R_{\unicode[STIX]{x1D706}}$ in the range 500–1100, on the axis of a plane jet. It is found that the magnitude of the exponent increases as $R_{\unicode[STIX]{x1D706}}$ increases and the rate of increase depends on the order $n$. The trend of published structure function data on the axes of an axisymmetric jet and a two-dimensional wake confirms this dependence. For a fixed $R_{\unicode[STIX]{x1D706}}$, the exponent can vary from flow to flow and for a given flow, the larger $R_{\unicode[STIX]{x1D706}}$ is, the closer the exponent is to the value predicted by Kolmogorov (Dokl. Akad. Nauk SSSR, vol. 30, 1941a, pp. 299–303) (hereafter K41). The major conclusion is that the finite Reynolds number effect, which depends on the flow, needs to be properly accounted for before determining whether corrections to K41, arising from the intermittency of the energy dissipation rate, are needed. We further point out that it is imprudent, if not incorrect, to associate the finite Reynolds number effect with a consequence of the modified similarity hypothesis introduced by Kolmogorov (J. Fluid Mech., vol. 13, 1962, pp. 82–85) (K62); we contend that this association has misled the vast majority of post K62 investigations of the consequences of K62.

Sign in / Sign up

Export Citation Format

Share Document