scholarly journals Modelling of turbulence modulation in particle- or droplet-laden flows

2012 ◽  
Vol 706 ◽  
pp. 251-273 ◽  
Author(s):  
Daniel W. Meyer
Keyword(s):  
Isotropic Turbulence ◽  
Density Ratio ◽  
Joint Probability ◽  
Stokes Number ◽  
Navier Stokes ◽  
Grid Turbulence ◽  
Heavy Particles ◽  
Turbulence Modulation ◽  
Modelling Framework ◽  
Liquid Flows

AbstractAddition of particles or droplets to turbulent liquid flows or addition of droplets to turbulent gas flows can lead to modulation of turbulence characteristics. Corresponding observations have been reported for very small particle or droplet volume loadings ${\Phi }_{v} $ and therefore may be important when simulating such flows. In this work, a modelling framework that accounts for preferential concentration and reproduces isotropic and anisotropic turbulence attenuation effects is presented. The framework is outlined for both Reynolds-averaged Navier–Stokes (RANS) and joint probability density function (p.d.f.) methods. Validations are performed involving a range of particle and flow-field parameters and are based on the direct numerical simulation (DNS) study of Boivin, Simonin & Squires (J. Fluid Mech., vol. 375, 1998, pp. 235–263) dealing with heavy particles suspended in homogeneous isotropic turbulence (Stokes number $\mathit{St}= O(1{\unicode{x2013}} 10)$, particle/fluid density ratio ${\rho }_{p} / \rho = 2000$, ${\Phi }_{v} = O(1{0}^{- 4} )$) and the experimental investigation of Poelma, Westerweel & Ooms (J. Fluid Mech., vol. 589, 2007, pp. 315–351) involving light particles ($\mathit{St}= O(0. 1)$, ${\rho }_{p} / \rho \gtrsim 1$, ${\Phi }_{v} = O(1{0}^{- 3} )$) settling in grid turbulence. The development in this work is restricted to volume loadings where particle or droplet collisions are negligible.

Validation of Particle-Laden Turbulent Flow Simulations Including Turbulence Modulation

2015 ◽  
Vol 137 (7) ◽  
Author(s):  
Yvonne Reinhardt ◽  
Leonhard Kleiser
Keyword(s):  
Numerical Study ◽  
Density Ratio ◽  
Particle Diameter ◽  
Stokes Number ◽  
Navier Stokes ◽  
Mean Flow ◽  
Turbulent Length Scale ◽  
Turbulence Modulation ◽  
Flow Simulations ◽  
Wide Range

The objective of the present numerical study is the validation of wall-bounded, turbulent particle-laden air flow simulations for a wide range of flow and particle parameters (i.e., flow and particle Reynolds numbers, Stokes number, particle-to-fluid density ratio, ratio of particle diameter to turbulent length scale) covering the one-, two- and four-way coupling regimes. The applied computational fluid dynamics (CFD) model follows the Eulerian two-fluid approach in a Reynolds-Averaged Navier–Stokes (RANS) context and is based on the kinetic theory of granular flow (KTGF) for closures concerning the particulate phase. The fluid turbulence is modeled applying a low-Reynolds-number k–epsilon turbulence model. The main focus is put on the modeling of turbulence coupling between the fluid and the particle phase. Different from common practice, the choice of a model accounting for turbulence modulation is made dependent on the prevailing coupling regime. For the case of four-way coupling, a new modulation model is suggested that well predicts turbulence augmentation and attenuation. The predictive capabilities of the present approach are evaluated by comparing simulation results to experimental benchmark data of various pipe and channel flows. Very good agreement with reference data is obtained for the mean flow and turbulence profiles of both phases.

An Investigation of Droplet Diffusion in Isotropic Turbulence With Digital Holographic PIV

Volume 3 ◽  
2004 ◽  
Author(s):  
Edwin Malkiel ◽  
Jian Sheng ◽  
David Garber ◽  
Joseph Katz
Keyword(s):  
Time Series ◽  
Isotropic Turbulence ◽  
Diffusion Rate ◽  
Density Ratio ◽  
Digital Camera ◽  
Stokes Number ◽  
Sample Volume ◽  
Tracer Particles ◽  
Fluid Velocities ◽  
Digital Holograms

In-line digital holography is utilized to measure the Lagrangian trajectory of droplets in locally isotropic turbulence. The objective of these measurements is to determine the diffusion rate of these droplets as a function of density ratio between the continuous and dispersed phases, Stokes number and turbulence level relative to the quiescent settling/rise velocity of the droplets. The present experiments are conducted using diesel fuel with diameters of 0.5–2 mm, specific gravity of 0.85 and Stokes number in the 0.2 to 5 range. The droplets are injected into a 50 × 50 × 50 mm sample volume located in the center of a 160 1 tank. The turbulence is generated by four spinning grids, located symmetrically in the corners of the tank. Planar PIV has been used to characterize the turbulence prior to the experiments. A time series of in-line digital holograms is recorded at 2000 frames per second using a 1000×1000 pixel digital camera by back illuminating the sample volume with a collimated laser beam. Numerical reconstruction generates a time series of high-resolution images of the droplets and tracer particles throughout the sample volume. Subsequent analysis is used to obtain the velocity along the droplet trajectory. Lagrangian correlations can then be used for calculating the diffusion rate of these droplets. In a smaller sample volume, with a 15×15 mm cross section, and by using localized seeding, we can also simultaneously measure the droplet velocity along with the velocity of the fluid in the vicinity of this droplet. The results provide statistics on the correlations between the droplet and fluid velocities.

scholarly journals Tempered fractional LES modeling

2021 ◽  
Vol 932 ◽  
Author(s):  
Mehdi Samiee ◽  
Ali Akhavan-Safaei ◽  
Mohsen Zayernouri
Keyword(s):  
Isotropic Turbulence ◽  
Stokes Equations ◽  
A Priori ◽  
Correlation Coefficients ◽  
Scale Model ◽  
Navier Stokes ◽  
Model Parameters ◽  
Subgrid Scale ◽  
Modelling Framework ◽  
Non Local

The presence of non-local interactions and intermittent signals in the homogeneous isotropic turbulence grant multi-point statistical functions a key role in formulating a new generation of large-eddy simulation (LES) models of higher fidelity. We establish a tempered fractional-order modelling framework for developing non-local LES subgrid-scale models, starting from the kinetic transport. We employ a tempered Lévy-stable distribution to represent the source of turbulent effects at the kinetic level, and we rigorously show that the corresponding turbulence closure term emerges as the tempered fractional Laplacian, $(\varDelta +\lambda )^{\alpha } (\cdot )$ , for $\alpha \in (0,1)$ , $\alpha \neq \frac {1}{2}$ and $\lambda >0$ in the filtered Navier–Stokes equations. Moreover, we prove the frame invariant properties of the proposed model, complying with the subgrid-scale stresses. To characterize the optimum values of model parameters and infer the enhanced efficiency of the tempered fractional subgrid-scale model, we develop a robust algorithm, involving two-point structure functions and conventional correlation coefficients. In an a priori statistical study, we evaluate the capabilities of the developed model in fulfilling the closed essential requirements, obtained for a weaker sense of the ideal LES model (Meneveau, Phys. Fluids, vol. 6, issue 2, 1994, pp. 815–833). Finally, the model undergoes the a posteriori analysis to ensure the numerical stability and pragmatic efficiency of the model.

Direct numerical simulation of turbulence modulation by particles in compressible isotropic turbulence

2017 ◽  
Vol 832 ◽  
pp. 438-482 ◽  
Author(s):  
Qi Dai ◽  
Kun Luo ◽  
Tai Jin ◽  
Jianren Fan
Keyword(s):  
Kinetic Energy ◽  
Turbulent Kinetic Energy ◽  
Rotational Motion ◽  
Isotropic Turbulence ◽  
Volume Fraction ◽  
Stokes Number ◽  
Systematic Investigation ◽  
Turbulence Modulation ◽  
Physical Mechanisms ◽  
Preferential Concentration

In this paper, a systematic investigation of turbulence modulation by particles and its underlying physical mechanisms in decaying compressible isotropic turbulence is performed by using direct numerical simulations with the Eulerian–Lagrangian point-source approach. Particles interact with turbulence through two-way coupling and the initial turbulent Mach number is 1.2. Five simulations with different particle diameters (or initial Stokes numbers, $St_{0}$) are conducted while fixing both their volume fraction and particle densities. The underlying physical mechanisms responsible for turbulence modulation are analysed through investigating the particle motion in the different cases and the transport equations of turbulent kinetic energy, vorticity and dilatation, especially the two-way coupling terms. Our results show that microparticles ($St_{0}\leqslant 0.5$) augment turbulent kinetic energy and the rotational motion of fluid, critical particles ($St_{0}\approx 1.0$) enhance the rotational motion of fluid, and large particles ($St_{0}\geqslant 5.0$) attenuate turbulent kinetic energy and the rotational motion of fluid. The compressibility of the turbulence field is suppressed for all the cases, and the suppression is more significant if the Stokes number of particles is close to 1. The modifications of turbulent kinetic energy, the rotational motion and the compressibility are all related with the particle inertia and distributions, and the suppression of the compressibility is attributed to the preferential concentration and the inertia of particles.

scholarly journals Inertial particle trapping in an open vortical flow

2014 ◽  
Vol 744 ◽  
pp. 183-216 ◽  
Author(s):  
Jean-Régis Angilella ◽  
Rafael D. Vilela ◽  
Adilson E. Motter
Keyword(s):  
Vortex Pair ◽  
Stokes Number ◽  
Basins Of Attraction ◽  
Critical Distance ◽  
Vortical Flow ◽  
Navier Stokes ◽  
Particle Trapping ◽  
Order Unity ◽  
Heavy Particles ◽  
Open Flow

AbstractRecent numerical results on advection dynamics have shown that particles denser than the fluid can remain trapped indefinitely in a bounded region of an open fluid flow. Here, we investigate this counterintuitive phenomenon both numerically and analytically to establish the conditions under which the underlying particle-trapping attractors can form. We focus on a two-dimensional open flow composed of a pair of vortices and its specular image, which is a system we represent as a vortex pair plus a wall along the symmetry line. Considering particles that are much denser than the fluid, referred to as heavy particles, we show that two attractors form in the neighbourhood of the vortex pair provided that the particle Stokes number is smaller than a critical value of order unity. In the absence of the wall, the attractors are fixed points in the frame rotating with the vortex pair, and the boundaries of their basins of attraction are smooth. When the wall is present, the point attractors describe counter-rotating ellipses in this frame, with a period equal to half the period of one isolated vortex pair. The basin boundaries remain smooth if the distance from the vortex pair to the wall is large. However, these boundaries are shown to become fractal if the distance to the wall is smaller than a critical distance that scales with the inverse square root of the Stokes number. This transformation is related to the breakdown of a separatrix that gives rise to a heteroclinic tangle close to the vortices, which we describe using a Melnikov function. For an even smaller distance to the wall, we demonstrate that a second separatrix breaks down and a new heteroclinic tangle forms farther away from the vortices, at the boundary between the open and closed streamlines. Particles released in the open part of the flow can approach the attractors and be trapped permanently provided that they cross the two separatrices, which can occur under the effect of flow unsteadiness. Furthermore, the trapping of heavy particles from the open flow is shown to be robust to the presence of viscosity, noise and gravity. Navier–Stokes simulations for large flow Reynolds numbers show that viscosity does not destroy the attracting points until vortex merging takes place, while simulation of thermal noise shows that particle trapping persists for extended periods provided that the Péclet number is large. The presence of a gravitational field does not alter the permanent trapping by the attracting points if the settling velocities are not too large. For larger settling velocities, however, gravity can also give rise to a limit-cycle attractor next to the external separatrix and to a new form of trapping from the open flow that is not mediated by a heteroclinic tangle.

scholarly journals Intermittency in the relative separations of tracers and of heavy particles in turbulent flows

2014 ◽  
Vol 757 ◽  
pp. 550-572 ◽  
Author(s):  
L. Biferale ◽  
A. S. Lanotte ◽  
R. Scatamacchia ◽  
F. Toschi
Keyword(s):  
Reynolds Number ◽  
Turbulent Flows ◽  
Isotropic Turbulence ◽  
Exit Time ◽  
Three Dimensional ◽  
Stokes Number ◽  
Inertial Range ◽  
Relative Dispersion ◽  
Heavy Particles ◽  
Self Similar

AbstractResults from direct numerical simulations (DNS) of particle relative dispersion in three-dimensional homogeneous and isotropic turbulence at Reynolds number $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}{\mathit{Re}}_{\lambda } \sim 300$ are presented. We study point-like passive tracers and heavy particles, at Stokes number $\mathit{St}=0.6,1$ and 5. Particles are emitted from localised sources, in bunches of thousands, periodically in time, allowing an unprecedented statistical accuracy to be reached, with a total number of events for two-point observables of the order of ${10^{11}}$. The right tail of the probability density function (PDF) for tracers develops a clear deviation from Richardson’s self-similar prediction, pointing to the intermittent nature of the dispersion process. In our numerical experiment, such deviations are manifest once the probability to measure an event becomes of the order of – or rarer than – one part over one million, hence the crucial importance of a large dataset. The role of finite-Reynolds-number effects and the related fluctuations when pair separations cross the boundary between viscous and inertial range scales are discussed. An asymptotic prediction based on the multifractal theory for inertial range intermittency and valid for large Reynolds numbers is found to agree with the data better than the Richardson theory. The agreement is improved when considering heavy particles, whose inertia filters out viscous scale fluctuations. By using the exit-time statistics we also show that events associated with pairs experiencing unusually slow inertial range separations have a non-self-similar PDF.

Coherent clusters of inertial particles in homogeneous turbulence

2017 ◽  
Vol 833 ◽  
pp. 364-398 ◽  
Author(s):  
Lucia Baker ◽  
Ari Frankel ◽  
Ali Mani ◽  
Filippo Coletti
Keyword(s):  
Isotropic Turbulence ◽  
Stokes Number ◽  
Inertial Particles ◽  
Gravitational Acceleration ◽  
Self Similarity ◽  
Integral Scale ◽  
Heavy Particles ◽  
Diagram Method ◽  
Vorticity Vector ◽  
Definition Of

Despite the widely acknowledged significance of turbulence-driven clustering, a clear topological definition of particle cluster in turbulent dispersed multiphase flows has been lacking. Here we introduce a definition of coherent cluster based on self-similarity, and apply it to distributions of heavy particles in direct numerical simulations of homogeneous isotropic turbulence, with and without gravitational acceleration. Clusters show self-similarity already at length scales larger than twice the Kolmogorov length, as indicated by the fractal nature of their surface and by the power-law decay of their size distribution. The size of the identified clusters extends to the integral scale, with average concentrations that depend on the Stokes number but not on the cluster dimension. Compared to non-clustered particles, coherent clusters show a stronger tendency to sample regions of high strain and low vorticity. Moreover, we find that the clusters align themselves with the local vorticity vector. In the presence of gravity, they tend to align themselves vertically and their fall speed is significantly different from the average settling velocity: for moderate fall speeds they experience stronger settling enhancement than non-clustered particles, while for large fall speeds they exhibit weakly reduced settling. The proposed approach for cluster identification leverages the Voronoï diagram method, but is also compatible with other tessellation techniques such as the classic box-counting method.

scholarly journals Correlation between Buoyancy Flux, Dissipation and Potential Vorticity in Rotating Stratified Turbulence

Atmosphere ◽  
2021 ◽  
Vol 12 (2) ◽  
pp. 157
Author(s):  
Duane Rosenberg ◽  
Annick Pouquet ◽  
Raffaele Marino
Keyword(s):  
Energy Dissipation ◽  
Kinetic Energy ◽  
Potential Vorticity ◽  
Isotropic Turbulence ◽  
Joint Probability ◽  
Buoyancy Flux ◽  
Turbulent Regime ◽  
Stratified Turbulence ◽  
Low Dissipation ◽  
Linear And Nonlinear

We study in this paper the correlation between the buoyancy flux, the efficiency of energy dissipation and the linear and nonlinear components of potential vorticity, PV, a point-wise invariant of the Boussinesq equations, contrasting the three identified regimes of rotating stratified turbulence, namely wave-dominated, wave–eddy interactions and eddy-dominated. After recalling some of the main novel features of these flows compared to homogeneous isotropic turbulence, we specifically analyze three direct numerical simulations in the absence of forcing and performed on grids of 10243 points, one in each of these physical regimes. We focus in particular on the link between the point-wise buoyancy flux and the amount of kinetic energy dissipation and of linear and nonlinear PV. For flows dominated by waves, we find that the highest joint probability is for minimal kinetic energy dissipation (compared to the buoyancy flux), low dissipation efficiency and low nonlinear PV, whereas for flows dominated by nonlinear eddies, the highest correlation between dissipation and buoyancy flux occurs for weak flux and high localized nonlinear PV. We also show that the nonlinear potential vorticity is strongly correlated with high dissipation efficiency in the turbulent regime, corresponding to intermittent events, as observed in the atmosphere and oceans.

scholarly journals Refinement of a Previous Hypothesis of the Lyapunov Analysis of Isotropic Turbulence

2013 ◽  
Vol 2013 ◽  
pp. 1-4 ◽  
Author(s):  
Nicola de Divitiis
Keyword(s):  
Structure Function ◽  
Isotropic Turbulence ◽  
Stokes Equations ◽  
Longitudinal Velocity ◽  
Navier Stokes ◽  
Lyapunov Analysis ◽  
Velocity Difference ◽  
Navier Stokes Equations ◽  
Previous Hypothesis ◽  
Kolmogorov Theory

The purpose of this paper is to improve a hypothesis of the previous work of N. de Divitiis (2011) dealing with the finite-scale Lyapunov analysis of isotropic turbulence. There, the analytical expression of the structure function of the longitudinal velocity differenceΔuris derived through a statistical analysis of the Fourier transformed Navier-Stokes equations and by means of considerations regarding the scales of the velocity fluctuations, which arise from the Kolmogorov theory. Due to these latter considerations, this Lyapunov analysis seems to need some of the results of the Kolmogorov theory. This work proposes a more rigorous demonstration which leads to the same structure function, without using the Kolmogorov scale. This proof assumes that pair and triple longitudinal correlations are sufficient to determine the statistics ofΔurand adopts a reasonable canonical decomposition of the velocity difference in terms of proper stochastic variables which are adequate to describe the mechanism of kinetic energy cascade.

Snowflakes in the atmospheric surface layer: observation of particle–turbulence dynamics

2017 ◽  
Vol 814 ◽  
pp. 592-613 ◽  
Author(s):  
Andras Nemes ◽  
Teja Dasari ◽  
Jiarong Hong ◽  
Michele Guala ◽  
Filippo Coletti
Keyword(s):  
Large Scale ◽  
Isotropic Turbulence ◽  
Volume Fraction ◽  
Response Times ◽  
Field Measurements ◽  
Stokes Number ◽  
Inertial Particles ◽  
Natural Setting ◽  
Scale Particle ◽  
Particle Imaging

We report on optical field measurements of snow settling in atmospheric turbulence at $Re_{\unicode[STIX]{x1D706}}=940$. It is found that the snowflakes exhibit hallmark features of inertial particles in turbulence. The snow motion is analysed in both Eulerian and Lagrangian frameworks by large-scale particle imaging, while sonic anemometry is used to characterize the flow field. Additionally, the snowflake size and morphology are assessed by digital in-line holography. The low volume fraction and mass loading imply a one-way interaction with the turbulent air. Acceleration probability density functions show wide exponential tails consistent with laboratory and numerical studies of homogeneous isotropic turbulence. Invoking the assumption that the particle acceleration has a stronger dependence on the Stokes number than on the specific features of the turbulence (e.g. precise Reynolds number and large-scale anisotropy), we make inferences on the snowflakes’ aerodynamic response time. In particular, we observe that their acceleration distribution is consistent with that of particles of Stokes number in the range $St=0.1{-}0.4$ based on the Kolmogorov time scale. The still-air terminal velocities estimated for the resulting range of aerodynamic response times are significantly smaller than the measured snow particle fall speed. This is interpreted as a manifestation of settling enhancement by turbulence, which is observed here for the first time in a natural setting.

Sign in / Sign up

Export Citation Format

Share Document