Agglomeration of inertial particles in a random rotating symmetric straining flow

This paper is concerned with the development and validation of a simple Lagrangian model for particle agglomeration in a turbulent flow involving the collision of particles in a sequence of correlated straining and vortical structures which simulate the Kolmogorov small scales of motion of the turbulence responsible for particle pair dispersion and collision. In this particular study we consider the collision rate of monodisperse spherical particles in a symmetric (pure) straining flow which is randomly rotated to create an isotropic flow. The model is similar to the classical model of Saffman and Turner (S&T) (1956) for the collision (agglomeration) of tracer particles suspended in a turbulent flow. However unlike S&T, the straining flow is not frozen in time persisting only for timescales ∼Kolmogorov timescale. Furthermore, we consider the collision of inertial particles as well as tracer particles, and study their behavior not only at the collision boundary but also in its vicinity. In the simulation, particles are injected continuously at the boundaries of the straining flow, the size of the straining region being typical of the Kolmogorov length scale η_K of the turbulence. For steady state conditions, we calculate the flux of particles colliding with a test particle at the centre of the straining flow and consider its dependence on the inertia of the colliding particles (characterized by the particle Stokes number, St). The model replicates the segregation and accumulation observed in DNS and in particular the maximum segregation for St ∼ 1 (where St is the ratio of the particle response time to the Kolmogorov timescale). We also calculate the contributions of the various turbulent forces in the momentum balance equation for satellite particles and show for instance that for small Stokes number, there is a balance between turbulent diffusion and turbophoresis (gradient of kinetic stresses) which in turn is responsible for the build-up of concentration at the collision boundary. As found in previous studies, for the case of inertialess tracer particles, the collision rate turns out to be significantly smaller than the S&T prediction due to a lowering of the concentration at the collision boundary compared to the fully mixed value. The increase in collision rate for St ∼ 0.5 is shown to be a combination of particle segregation (build-up of concentration near the collision boundary) and the decorrelation of the relative velocity between the local fluid and a colliding particle. The difference from the S&T value for the agglomeration kernel is shown to be a consequence of the choice of perfectly absorbing boundary conditions at collision and the influence of the time scale of the turbulence (eddy lifetime). We draw the analogy between turbulent agglomeration and particle deposition in a fully developed turbulent boundary layer.     

Dynamics of disk-like particles in turbulent vertical channel flow

The dynamical behavior of inertial disk-like particles in turbulent vertical channel flow is investigated by an Eulerian–Lagrangian point-particle approach. Gravity effects on distribution, translation, rotation and orientation statistics of non-spherical particles modeled as oblate spheroids have been studied both in an upward and a downward flow and compared with results obtained in the absence of gravity. Altogether 12 different particle classes have been studied, with inertia and shape parameterized by means of Stokes number St and aspect ratio λ ≤ 1. The St = 5 disk-like particles distribute more evenly across the channel in upward than in downward flow. The gravity effect on the particle concentration diminishes with large inertia and the spheroid shape has only a modest influence. Although the gravity significantly affects the streamwise and wall-normal mean slip velocities with increasing inertia, the particle shape rarely has any impact on the translational motion, except for the mean wall-normal velocity. The fluctuations of the velocity of disk-like particles are mainly ascribed to inertia, whereas the gravity and shape only have marginal effects. The presence of gravity is moreover found to have a negligible effect on the particles’ orientation and rotation, in spite of the striking effect of λ on the orientation and rotation seen in the near-wall region. The tendency of the disks to align their symmetry axis orthogonal to the fluid vorticity in the channel center is stronger for particles with modest inertia. In the near-wall region, however, oblate spheroids preferentially align with the fluid vorticity for St >> 1. The observed behavior is believed to be caused by the influence of the gravity force on the turbophoresis; i.e. that inertial particles move towards low-turbulence regions.     

A model for tracking inertial particles in a lattice Boltzmann turbulent flow simulation

A computationally inexpensive model for tracking inertial particles through a turbulent flow is presented and applied to the turbulent flow through a square duct having a friction Reynolds number of Re_τ = 300. Prior to introducing particles into the model, the flow is simulated using a lattice Boltzmann computation, which is allowed to evolve until a steady state turbulent flow is achieved. A snapshot of the flow is then stored, and the trajectories of particles are computed through the flow domain under the influence of this static probability field. Although the flow is not computationally evolving during the particle tracking simulation, the local velocity is obtained stochastically from the local probability function, thus allowing the dynamics of the turbulent flow to be resolved from the point of view of the suspended particles. Particle inertia is modeled by using a relaxation parameter based on the particle Stokes number that allows for a particle velocity history to be incorporated during each time step. Wall deposition rates and deposition patterns are obtained and exhibit a high level of agreement with previously obtained DNS computational results and experimental results for a wide range of particle inertia. These results suggest that accurate particle tracking through complex turbulent flows may be feasible given a suitable probability field, such as one obtained from a lattice Boltzmann simulation. This in turn presents a new paradigm for the rapid acquisition of particle transport statistics without the need for concurrent computations of fluid flow evolution.     

Direct numerical simulation of particle interaction with ejections in turbulent channel flows

Direct numerical simulations (DNS) of incompressible turbulent channel flows coupled with Lagrangian particle tracking are performed to study the characteristics of ejections that surround solid particles. The behavior of particles in dilute turbulent channel flows, without particle collisions and without feedback of particles on the carrier fluid, is studied using high Reynolds number DNS (Re = 12,500). The results show that particles moving away from the wall are surrounded by ejections, confirming previous studies on this issue. A threshold value separating ejections with only upward moving particles is established. When normalized by the square root of the Stokes number and the square of the friction velocity, the threshold profiles follow the same qualitative trends, for all the parameters tested in this study, in the range of the experiments. When compared to suspension thresholds proposed by other studies in the Shields diagram, our simulations predict a much larger value because of the measure used to characterize the fluid and the criterion chosen to decide whether particles are influenced by the surrounding fluid. However, for intermediate particle Reynolds numbers, the threshold proposed here is in fair agreement with the theoretical criterion proposed by Bagnold (1966) [Bagnold, R., 1966. Geological Survey Professional Paper, vol. 422-1]. Nevertheless, further studies will be conducted to understand the normalization of the threshold.     

Turbulent particle dispersion in arbitrary wall-bounded geometries: A coupled CFD-Langevin-equation based approach

A Lagrangian continuous random walk (CRW) model is developed to predict turbulent particle dispersion in arbitrary wall-bounded flows with prevailing anisotropic, inhomogeneous turbulence. The particle tracking model uses 3D mean flow data obtained from the Fluent CFD code, as well as Eulerian statistics of instantaneous quantities computed from DNS databases. The turbulent fluid velocities at the current time step are related to those of the previous time step through a Markov chain based on the normalized Langevin equation which takes into account turbulence inhomogeneities. The model includes a drift velocity correction that considerably reduces unphysical features common in random walk models. It is shown that the model satisfies the well-mixed criterion such that tracer particles retain approximately uniform concentrations when introduced uniformly in the domain, while their deposition velocity is vanishingly small, as it should be. To handle arbitrary geometries, it is assumed that the velocity rms values in the boundary layer can locally be approximated by the DNS data of fully developed channel flows. Benchmarks of the model are performed against particle deposition data in turbulent pipe flows, 90° bends, as well as more complex 3D flows inside a mouth-throat geometry. Good agreement with the data is obtained across the range of particle inertia.     

A stochastic Langevin model of turbulent particle dispersion in the presence of thermophoresis

Direct numerical simulation (DNS) and experimental data have shown that inertial particles exhibit concentration peaks in isothermal turbulent boundary layers, whereas tracer-like particles remain well mixed in the domain. It is therefore expected that the interactions between turbulence and thermophoresis will be strong in particle-laden flows where walls and carrier fluid are at significantly different temperatures. To capture turbulent particle dispersion with active thermophoresis, a coupled CFD-Lagrangian continuous random walk (CRW) model is developed. The model uses 3D mean flow velocities obtained from the Fluent 6.3 CFD code, to which are added turbulent fluid velocities derived from the normalized Langevin equation which accounts for turbulence inhomogeneities. The mean thermophoretic force is included as a body force on the particle following the Talbot formulation. Validation of the model is performed against recent integral thermophoretic deposition data in long pipes as well as the TUBA TT28 test with its detailed local deposition measurements. In all cases, the agreement with the data is very good. In separate parametric studies in a hypothetical cooled channel flow, it is found that turbulence strongly enhances thermophoretic deposition of particles with dimensionless relaxation times τ⁺ of order 1 or more. On the other hand, the thermophoretic deposition of very small inertia particles (τ⁺ < 0.2) in the asymptotic region far from the injection point tends to that which characterizes stagnant flow conditions, in agreement with the DNS results of Thakurta et al.     

Filtered particle tracking in isotropic turbulence and stochastic modeling of subgrid-scale dispersion

A numerical study based on the Eulerian–Lagrangian formulation is performed for dispersed phase motion in a turbulent flow. The effect of spatial filtering, commonly employed in large-eddy simulations, and the role of the subgrid scale turbulence on the statistics of heavy particles, including preferential concentration, are studied through a priori analysis of DNS of particle-laden forced isotropic turbulence. In simulations where the subgrid scale kinetic energy attains 30–35% of the total we observe the impact of residual fluid motions on particles of a smaller inertia. It is shown that neglecting the influence of subgrid scale fluctuations has a significant effect on the preferential concentration of those particles. A stochastic Langevin model is proposed to reconstruct the residual (or subgrid scale) fluid velocity along particle trajectories. The computation results for a selection of particle inertia parameters are performed to appraise the model through comparisons of particle turbulent kinetic energy and the statistics of preferential concentrations.     

Modeling droplet dispersion and interphase turbulent kinetic energy transfer using a new dual-timescale Langevin model

Dispersion of spray droplets and the modulation of turbulence in the ambient gas by the dispersing droplets are two coupled phenomena that are closely linked to the evolution of global spray characteristics, such as the spreading rate of the spray and the spray cone angle. Direct numerical simulations (DNS) of turbulent gas flows laden with sub-Kolmogorov size particles, in the absence of gravity, report that dispersion statistics and turbulent kinetic energy (TKE) evolve on different timescales. Furthermore, each timescale behaves differently with Stokes number, a non-dimensional flow parameter (defined in this context as the ratio of the particle response time to the Kolmogorov timescale of turbulence) that characterizes how quickly a particle responds to turbulent fluctuations in the carrier or gas phase. A new dual-timescale Langevin model (DLM) composed of two coupled Langevin equations for the fluctuating velocities, one for each phase, is proposed. This model possesses a unique feature that the implied TKE and velocity autocorrelation in each phase evolve on different timescales. Consequently, this model has the capability of simultaneously predicting the disparate Stokes number trends in the evolution of dispersion statistics, such as velocity autocorrelations, and TKE in each phase. Predictions of dispersion statistics and TKE from the new model show good agreement with published DNS of non-evaporating and evaporating droplet-laden turbulent flow.     

Direct numerical simulation of a particle-laden low Reynolds number turbulent round jet

Direct numerical simulation method is used for the investigating of particle-laden turbulent flows in a spatially evolution of low Reynolds number axisymmetric jet, and the Eulerian–Lagrangian point-particle approach is employed in the simulation. The simulation uses an explicit coupling scheme between particles and the fluid, which considers two-way coupling between the particle and the fluid. The DNS results are compared well with experimental data with equal Reynolds number (R_e = 1700). Our objects are: (i) to investigate the correlation between the particle number density and the fluctuating of fluid streamwise velocity; (ii) to examine whether the three-dimensional vortex structures in the particle-laden jet are the same as that in the free-air jet and how the particles modulate the thee-dimensional vortex structures and turbulence properties with different Stokes number particles; (iii) to discover the particle circumferential dispersion with different Stokes number particles. Our findings: (i) all the particles, regardless of their particle size, tend to preferentially accumulate in the region with large-than-mean fluid streamwise velocity; (ii) the small Stokes number particles take an important part in the modulation of three-dimensional vortex structures, but for the intermediate and larger sized particles, this modulation effect seems not so apparent; (iii) the particle circumferential dispersion is more effective for the smaller and intermediate sized particles, especially for the intermediate sized particles.     

Large Scale Accumulation Patterns of Inertial Particles in Wall-Bounded Turbulent Flow

Turbulent internal flow in channel and pipe geometry with a diluted second phase of inertial particles is studied numerically. Direct numerical simulations (DNS) are performed at moderate Reynolds number (Re _?? ????200) in pipe and two channels??a smaller one similar in size to previous studies and a 3?×?3-times larger one??and Eulerian statistics pertaining to the particle concentration are evaluated. This simulation box constitutes the largest domain used for particle-laden flows so far. The resulting two-point correlations of the particle concentration show that in the smaller channel the particles organize in thin, streamwise elongated patterns which are very regular and long. The spanwise spacing of these structures is 120 and 160 plus units for the channel and pipe, respectively. Only in the larger box, the streamwise extent is long enough for the particle streaks to decorrelate, thus allowing the particles to move more freely. The influence of the box size on the characteristics of the turbophoresis is clearly shown; a 10% increase of the near-wall correlation is observed for particles with Stokes number St ^?+??=?50. It is thus shown that the box dimensions are an important factor in correctly assessing the motion of inertial particles, and their relation to the underlying velocity field. In addition the binning size effects on the correlation statistics of particle concentration are exploited. In particular the spanwise correlation peak values appear very sensitive to the adopted binning size, although the position of these peaks is found almost independent. Hence to allow a significant comparison between data of different configurations it is necessary to adopt the same binning spacing in inner variable.     

A Reynolds stress model for turbulent flow of homogeneous polymer solutions

Using a priori analyses of direct numerical simulation (DNS) data, a Reynolds stress model (RSM) is developed to account for the influence of polymer additives on turbulent flow over a wide range of flow conditions. The Finitely Extensible Nonlinear Elastic-Peterlin (FENE-P) rheological constitutive model is utilized to evaluate the polymer contribution to the stress tensor. Thirteen DNS data sets are used to analyze the budgets of elastic stress–velocity gradient correlations as well as Reynolds stress and dissipation transport. Closures are developed in the framework of the RSM model for all the required unknown and non-linear terms. The polymer stresses, velocity profiles, turbulent flow statistics and the percentage of friction drag reduction predicted by the RSM model are in good agreement with present and those obtained from independent DNS data over a wide range of rheological and flow parameters.     

Numerical modeling of a two-dimensional vertical turbulent two-phase jet

The results of numerically modeling two-dimensional two-phase flow of the “gas-solid particles” type in a vertical turbulent jet are presented for three cases of its configuration, namely, descending, ascending, and without account of gravity. Both flow phases are modeled on the basis of the Navier-Stokes equations averaged within the framework of the Reynolds approximation and closed by an extended k-? turbulence model. The averaged two-phase flow parameters (particle and gas velocities, particle concentration, turbulent kinetic energy, and its dissipation) are described using the model of mutually-penetrating continua. The model developed allows for both the direct effect of turbulence on the motion of disperse-phase particles and the inverse effect of the particles on turbulence leading to either an increase or a decrease in the turbulent kinetic energy of the gas. The model takes account for gravity, viscous drag, and the Saffman lift. The system of equations is solved using a difference method. The calculated results are in good agreement with the corresponding experimental data which confirms the effect of solid particles on the mean and turbulent characteristics of gas jets.     

Dispersion of particles in wall-bounded particle-laden turbulent flows with high wall permeability

A series of numerical simulations were performed to investigate the distribution and deposition properties of particles in turbulent flows bounded by permeable walls using the Large Eddy Simulation (LES) with a Lagrangian trajectory approach. The wall permeation speeds were taken from 10⁻⁴ to 10⁻² of the bulk velocity. The directions of the permeation speed were the same at both walls, and they were inward on one wall but outward on the other wall to reserve the fluid mass. Particles with Stokes number (respecting viscous time scale) around 0.1, 1 and 10 were released in the fully developed turbulent channel flow. The particle–particle interaction and the retroaction from particles to the fluid were neglected. The fluid-phase turbulence statistical properties and particle's transport characteristics by vortexes were then analyzed in details. If the wall permeation exists, the turbulence intensities will be depressed close to the outward permeable wall but increased near the inward permeable wall. Not influenced by the wall permeation, the suspended particles with St⁺ ∼O(1) tend to accumulate in the less vortical zones away from the wall, while those particles in the flow regions near the outward permeable wall will distribute disregarding of the vorticity. The turbulence structures near the outward permeable wall are found to exert promotional effects on the particle deposition rate, but such effects are different for particles with various Stokes number. A distribution tendency of streamwise streaks for the deposited particles is also found on the wall imposed by the high outward permeation speed and the clustering deposition pattern is more obvious with increasing particle Stokes number.     

Stochastic modeling of particle diffusion in a turbulent boundary layer

Several Continuous Random Walk (CRW) models were constructed to predict turbulent particle diffusion based on Eulerian statistics that can be obtained with Reynolds-Averaged Navier Stokes (RANS) solutions. The test conditions included a wide range of particle inertias (Stokes numbers) with a near-wall injection (y⁺ = 4) in a turbulent boundary layer that is strongly anisotropic and inhomogeneous. To assess the performance of the models, the CRW results were compared to particle diffusion statistics gathered from a Direct Numerical Simulation (DNS). In particular, comparisons were made with transverse concentration profiles, root-mean-square of particle trajectory coordinates, and mean transverse particle velocity away from the wall.The results showed that accurate simulation required a modified (non-dimensionalized) Markov chain to handle the large gradients in turbulence near the wall as shown by simulations with fluid-tracer particles. For finite-inertia particles, an incremental drift correction for the Markov chain developed herein to account for Stokes number effects was critical to avoiding non-physical particle collection in low-turbulence regions. In both cases, inclusion of anisotropy in the turbulence model was found to be important, but the influence of off-diagonal terms was found to be weak. The results were generally good, especially for long-time and large inertia particles.     

Effect of particles on turbulent thermal field of channel flow with different Prandtl numbers

The direct numerical simulation(DNS) of heat transfer in a fully developed non-isothermal particle-laden turbulent channel flow is performed.The focus of this paper is on the modulation of the particles on turbulent thermal statistics in the particle-laden flow with three Prandtl numbers(P r = 0.71,1.5,and 3.0) and a shear Reynolds number(Reτ = 180).Some typical thermal statistics,including normalized mean temperature and their fluctuations,turbulent heat fluxes,Nusselt number and so on,are analyzed.The results show that the particles have less effects on turbulent thermal fields with the increase of Prandtl number.Two reasons can explain this.First,the correlation between fluid thermal field and velocity field decreases as the Prandtl number increases,and the modulation of turbulent velocity field induced by the particles has less influence on the turbulent thermal field.Second,the heat exchange between turbulence and particles decreases for the particle-laden flow with the larger Prandtl number,and the thermal feedback of the particles to turbulence becomes weak.     

Statistics of particle dispersion in direct numerical simulations of wall-bounded turbulence: Results of an international collaborative benchmark test

In this paper the results of an international collaborative test case relative to the production of a direct numerical simulation and Lagrangian particle tracking database for turbulent particle dispersion in channel flow at low Reynolds number are presented. The objective of this test case is to establish a homogeneous source of data relevant to the general problem of particle dispersion in wall-bounded turbulence. Different numerical approaches and computational codes have been used to simulate the particle-laden flow and calculations have been carried on long enough to achieve a statistically steady condition for particle distribution. In such stationary regime, a comprehensive database including both post-processed statistics and raw data for the fluid and for the particles has been obtained. The complete datasets can be downloaded from the web at http://cfd.cineca.it/cfd/repository/. In this paper the most relevant velocity statistics (for both phases) and particle distribution statistics are discussed and benchmarked by direct comparison between the different numerical predictions.     

各向同性湍流内颗粒碰撞率的直接模拟研究

对 Re_{\lambda } 约为51均匀各向同性湍流内 St_{k}(=\tau_{p}/\tau_{k}) 为 0 ~10.0 的 有限惯性颗粒的碰撞行为进行了直接数值模拟，以研究湍流对有限惯性 颗粒碰撞的影响. 结果表明，具有一定惯性颗粒的湍流碰撞率完全不同于零惯性的轻颗粒 (St_{k}=0) 和可忽略湍流作用的重颗粒 (St{k} \to \infty) , 其变化趋势极其复杂： 在Stk为 0~1.0 之间，颗粒的碰撞率随 St 的增加而近乎线性地剧烈增长，在 Stk≈1.0 3.0(对应的StE=τp/Te≈0.5)附近，颗粒碰撞率出现两个峰值，在Stk>3.0以后，颗粒的碰撞率随惯性增 大而逐渐趋向于重颗粒极限；在峰值处，有限惯性颗粒的平均碰撞率的峰值较轻颗粒增强了 30倍左右. 为进一步分析湍流作用下颗粒碰撞率的影响因素，分别使用可能发生碰撞 的颗粒对的径向分布函数和径向相对速度来量化颗粒的局部富集效应和湍流掺混效应，表明 St_{k} \approx 1.0 时局部富集效应最为强烈，使得颗粒的碰撞率出现第1个峰值； 湍流掺混效应则随着颗粒Stk的增大而渐近增大；局部富集和湍流掺混联合作用的结果， 使得颗粒碰撞率在 St_{k} \approx 3.0 附近出现另一个峰值.     

Direct numerical simulation of a particle-laden mixing layer with a chemical reaction

A direct numerical simulation (DNS) is applied to a particle-laden turbulent mixing layer with a chemical reaction, and the effects of particles on turbulence and chemical species’ diffusion and reaction in both zero and finite gravity cases are investigated. Unreactive particles, whose response time, τ_P, is smaller than the Kolmogorov time scale, τ_K [τ_P/τ_K = O(10⁻¹)], are uniformly injected into the high-speed side of the mixing layer. Two reactive chemical species are separately introduced through different sides. The results show that although laden particles generally depress turbulent intensities, they begin to enhance turbulent intensities downstream as the particle size decreases provided that the inlet particle volume fraction is fixed. This is because that the small particles with small particle response time easily accumulate at the circumference of coherent vortices and act to suppress the growth of primitive small-scale coherent vortices upstream but enhance that of relatively developing large-scale ones downstream. Also, since the small-scale turbulence, which promotes the chemical reaction, is suppressed by the laden particles in the entire region, chemical product decreases overall. Furthermore, the presence of finite gravity on the particles acts to depress the turbulent intensities, but its effects on chemical species’ diffusion and reaction are quite small.     

Measurement of inertial particle clustering and relative velocity statistics in isotropic turbulence using holographic imaging

We present the first measurements of relative velocity statistics of inertial particles in a homogeneous isotropic turbulent flow with three-dimensional holographic particle image velocimetry (holographic PIV). From the measurements we are able to obtain the radial relative velocity probability density function (PDF) conditioned on the interparticle separation distance, for distances on the order of the Kolmogorov length scale. Together with measurements of the three-dimensional radial distribution function (RDF) in our turbulence chamber, these statistics, in principle, can be used to determine interparticle collision rates via the formula derived by Sundaram and Collins (1997). In addition, we show temporal development of the RDF, which reveals the existence of an extended quasi-steady-state regime in our facility. Over this regime the measured two-particle statistics are compared to direct numerical simulations (DNS) with encouraging qualitative agreement. Statistics at the same Reynolds number but different Stokes numbers demonstrate the ability of the experiment to correctly capture the trends associated with particles of different inertia. Our results further indicate that even at moderate Stokes numbers turbulence may enhance collision rates significantly. Such experimental investigations may prove valuable in validating, guiding and refining numerical models of particle dynamics in turbulent flows.