Particle and droplet deposition in turbulent swirled pipe flow

In this work we study deposition of particles and droplets in non-rotating swirled turbulent pipe flow. We aim at verifying whether the capability of swirl to enhance particle separation from the core flow and the capability of turbulence to efficiently trap particles at the wall can co-exist to optimize collection efficiency in axial separators. We perform an Eulerian–Lagrangian study based on Direct Numerical Simulation (DNS) of turbulence, considering the effect of different swirl intensities on turbulence structures and on particle transfer at varying particle inertia. We show that, for suitably-chosen flow parameters, swirl may be superimposed to the base flow without disrupting near-wall turbulent structures and their regeneration mechanisms. We also quantify collection efficiency demonstrating for the first time that an optimal synergy between swirl and wall turbulence can be identified to promote separation of particles and droplets.     

Statistical models for predicting the effect of bidisperse particle collisions on particle velocities and stresses in homogeneous anisotropic turbulent flows

The purpose of this paper is to present and compare two statistical models for predicting the effect of collisions on particle velocities and stresses in bidisperse turbulent flows. These models start from a kinetic equation for the probability density function (PDF) of the particle velocity distribution in a homogeneous anisotropic turbulent flow. The kinetic equation describes simultaneously particle–turbulence and particle–particle interactions. The paper is focused on deriving the collision terms in the governing equations of the PDF moments. One of the collision models is based on a Grad-like expansion for the PDF of the velocity distributions of two particles. The other model stems from a Grad-like expansion for the joint fluid–particle PDF. The validity of these models is explored by comparing with Lagrangian simulations of particle tracking in uniformly sheared and isotropic turbulent flows generated by LES. Notwithstanding the fact that the fluid turbulence may be isotropic, the particle velocity fluctuations are anisotropic due to the impact of gravitational settling. Comparisons of the model predictions and the numerical simulations show encouraging agreement.     

Interactions of particles and microbubbles with turbulence

Dilute, dispersed two-phase flows arise in many contexts ranging from solid particles or droplets in gas flows to bubbles in liquids. Many of the flows of interest are turbulent, which presents a complex problem to analyze or to determine the dominant physical processes contributing to the observed phenomena. Advances in experimental techniques have made it possible to measure directly turbulent and particle velocity fluctuations in dilute systems. This has provided a counterpart to advances in computational and analytical models and a basis on which to test these models. Three specific areas are considered: the fluctuating forces on an individual particle in an unsteady flow, the response of a solid particle to a turbulent air flow, and the corresponding response of a small bubble in turbulent liquid flows. Results from direct numerical simulations are presented for each of these, including the nonuniform distribution of particles generated by local instantaneous features of the flow. The issue of turbulence modulation at low to moderate void fractions is discussed.     

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.     

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.     

Two-phase micro- and macro-time scales in particle-laden turbulent channel flows

The micro-and macro-time scales in two-phaseturbulent channel flows are investigated using the direct numerical simulation and the Lagrangian particle trajectorymethods for the fluid-and the particle-phases,respectively.Lagrangian and Eulerian time scales of both phases are calculated using velocity correlation functions.Due to flowanisotropy,micro-time scales are not the same with the theoretical estimations in large Reynolds number(isotropic) turbulence.Lagrangian macro-time scales of particle-phaseand of fluid-phase seen by particles are both dependent onparticle Stokes number.The fluid-phase Lagrangian integral time scales increase with distance from the wall,longerthan those time scales seen by particles.The Eulerian integral macro-time scales increase in near-wall regions but decrease in out-layer regions.The moving Eulerian time scalesare also investigated and compared with Lagrangian integraltime scales,and in good agreement with previous measurements and numerical predictions.For the fluid particles themicro Eulerian time scales are longer than the Lagrangianones in the near wall regions,while away from the walls themicro Lagrangian time scales are longer.The Lagrangianintegral time scales are longer than the Eulerian ones.Theresults are useful for further understanding two-phase flowphysics and especially for constructing accurate predictionmodels of inertial particle dispersion.     

Modeling of dynamics, heat transfer, and combustion in two-phase turbulent flows: 1. Isothermal flows

This paper presents a review of authors' collective works in the field of two-phase flow modeling done in the past few decades. The paper is aimed at the construction of mathematical models for simulation of particle-laden turbulent flows. A kinetic equation was obtained for the probability density function (PDF) of the particle velocity distribution in turbulent flows. The proposed kinetic equation describes both the interaction of particles with turbulent eddies of the carrier phase and particle-particle collisions. This PDF equation is used for the derivation of different schemes describing turbulent momentum transfer in the dispersed particle phase. The turbulent characteristics of the gaseous phase are calculated on the basis of the k - turbulence model with a modulation effect of particles on the turbulence.

The constructed models have been applied to the calculation of various two-phase gas-particle turbulent flows in jets and channels as well as particle deposition in tubes and separators. For validating the theoretical and numerical results, a wide range of comparisons with experimental data from Russian and foreign sources has been done.     

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.     

Transport and deposition of colliding particles in turbulent channel flows

The purpose of this paper is twofold: (i) to present statistical models that describe particle–turbulence interactions as well as particle–particle collisions and (ii) to gain a better understanding of the effect of inter-particle collisions on transport, deposition, and preferential concentration of heavy particles in turbulent channel flows. The models presented are based on a kinetic equation for the probability density function of the particle velocity distribution in anisotropic turbulent flow. The model predictions compare reasonable well with numerical simulations and properly reproduce the crucial trends of computations.     

Particle drag in a dilute turbulent two-phase suspension flow

The drag between phases plays an important role in the study of a turbulent two-phase suspension flow and its physical understanding will greatly promote progress in theoretical treatments of a whole range of important industrial and technical problems involving such a flow. The conventional practice of using the results of measurements based on a single particle in a laminar stream for the case of a turbulent flow of a dilute suspension is questioned. An analysis of the results of local measurements of upward turbulent flows of a solid particle-air two-phase suspension leads to the determination of the realistic particle drag coefficient over a wide range of flow conditions. It is established that the particle drag can be described by the simple Stokes law, based on an apparent turbulent viscosity of the fluid for the particles in the suspension flow. A correlation is provided for this apparent turbulent viscosity in terms of the particle size and concentration in the suspension, the local flow turbulence Reynolds number and the particle-to-fluid density ratio.     

Turbulent collision of inertial particles: Point-particle based,hybrid simulations and beyond

Point-particle based direct numerical simulation (PPDNS) has been a productive research tool for studying both single-particle and particle-pair statistics of inertial particles suspended in a turbulent carrier flow. Here we focus on its use in addressing particle-pair statistics relevant to the quantification of turbulent collision rate of inertial particles. PPDNS is particularly useful as the interaction of particles with small-scale (dissipative) turbulent motion of the carrier flow is mostly relevant. Furthermore, since the particle size may be much smaller than the Kolmogorov length of the background fluid turbulence, a large number of particles are needed to accumulate meaningful pair statistics. Starting from the relative simple Lagrangian tracking of so-called ghost particles, PPDNS has significantly advanced our theoretical understanding of the kinematic formulation of the turbulent geometric collision kernel by providing essential data on dynamic collision kernel, radial relative velocity, and radial distribution function. A recent extension of PPDNS is a hybrid direct numerical simulation (HDNS) approach in which the effect of local hydrodynamic interactions of particles is considered, allowing quantitative assessment of the enhancement of collision efficiency by fluid turbulence. Limitations and open issues in PPDNS and HDNS are discussed. Finally, on-going studies of turbulent collision of inertial particles using large-eddy simulations and particle-resolved simulations are briefly discussed.     

Particle-turbulence interaction in a boundary layer

Particle-turbulence interaction in wall turbulent flows has been studied. A series of experiments varying particle size, particle density, particle loading and flow Re has been conducted. The results show that the larger polystyrene particles (1100 μm) cause an increase in the number of wall ejections, giving rise to an increase in the measured values of the turbulence intensities and Reynolds stresses. On the other hand, the smaller polystyrene particles (120 μm) bring about a decrease in the number of wall ejections, causing a decrease in the measured intensities and Reynolds stresses. These effects are enhanced as the particle loading is increased. It was also found that the heavier glass particles (88 μm) do not bring about any significant modulation of turbulence. In addition, measurements of the burst frequency and the mean streak-spacing show no significant change with increase in particle loading. Based on these observations, a mechanism of particle transport in wall turbulent flows has been proposed, in which the particles are transported (depending on their size, density and flow Re) by the bursting events of the wall regions.     

Experimental and numerical investigation of turbulent cross-flow in a staggered tube bundle

This paper presents the results of measurements and numerical predictions of turbulent cross-flow in a staggered tube bundle. The bundle consists of transverse and longitudinal pitch-to-diameter ratios of 3.8 and 2.1, respectively. The experiments were conducted using a particle image velocimetry technique, in a flow of water in a channel at a Reynolds number of 9300 based on the inlet velocity and the tube diameter. A commercial CFD code, ANSYS CFX V10.0, is used to predict the turbulent flow in the bundle. The steady and isothermal Reynolds–Averaged Navier–Stokes (RANS) equations were used to predict the turbulent flow using each of the following four turbulence models: a k-epsilon, a standard k-omega, a k-omega-based shear stress transport, and an epsilon-based second moment closure. The epsilon-based models used a scalable wall function and the omega-based models used a wall treatment that switches automatically between low-Reynolds and standard wall function formulations.

The experimental results revealed extremely high levels of turbulence production by the normal stresses, as well as regions of negative turbulence production. The convective transport by mean flow and turbulent diffusion were observed to be significantly higher than in classical turbulent boundary layers. As a result, turbulence production is generally not in equilibrium with its dissipation rate. In spite of these characteristics, it was observed that the Reynolds normal stresses approximated from the k-based two-equation models were in a closer agreement with experiments than values obtained from the second moment closure. The results show that none of the turbulence models was able to consistently reproduce the mean and turbulent quantities reasonably well. The omega-based models predicted the mean velocities better in the developing region while the epsilon-based models gave better results in the region where the flow is becoming spatially periodic.     

Application of Lattice-Boltzmann Method for Analysing Detachment of Micron-Sized Particles from Carrier Particles in Turbulent Flows

Numerical calculations based on the Lattice-Boltzmann method were performed for a particle cluster consisting of a large spherical carrier particle covered with hundreds of small spherical drug particles. This cluster, fixed in space within a cubic computational domain, was exposed to turbulent plug airflow with predefined intensity. Such a situation is found in dry powder inhalers where carrier particles blended with fine drug powder are dispersed in a highly turbulent flow with the objective of detaching the drug powder for pulmonary delivery. Turbulence was generated by a digital filtering technique applied to the inflow velocity boundary condition. This technique was first validated by analysing the turbulence intensity at 15 fluid nodes along the stream-wise direction of the computational domain. The size ratio between the drug and carrier particle was 5 μ m/100 μ m, and the coverage degree of the carrier by the small particles was 50%, which is a typical value for carrier particle blending. The range of carrier particle Reynolds numbers considered was between 80 and 200, typical values found in inhaler devices. Exemplarily, at Re = 200 turbulence intensity was varied from 0.3% to 9.0%. The systematic increase of the mean flow (i.e. 80 < Re <200) resulted in varying turbulence intensities from 20 to 9%. These simulations provided the temporal evolution of the fluid dynamic forces on the drug particles in dependence of their angular position on the carrier in order to estimate the possibility of drug particle detachment. For turbulent conditions (i.e. Re = 200 and I = 9.0%) the maximum fluid forces on the drug particles were found to be about 10-times larger than found in laminar flow. The fluctuations in the forces were found to be higher than the flow velocity fluctuations due to the modification of the boundary layer around the cluster and instabilities triggered by the turbulent flow. There are three possibilities for detaching the drug powder, namely, through lift-off and sliding or rolling. Lift-off was found to be of minor importance due to the observed small normal fluid forces even at Re = 200 and I = 9.0%. The probability of sliding and rolling detachment in dependence of the angular position was estimated based on measured adhesion properties, i.e. van der Waals force, adhesion surface energy and friction coefficient. The remarkable rise of detachment probability for both effects due to the action of turbulence is an important finding of this study. In accordance with laminar flow, rolling detachment occurs before sliding, however in turbulent conditions over the entire carrier particle. The present studies improve the understanding of drug particle detachment from carrier particles in an inhaler device. The results will be the basis for developing Lagrangian detachment models that eventually should allow the optimisation of dry powder inhalators through computational fluid dynamics.     

Validation of a numerical model for the simulation of an electrostatic powder coating process

Numerical modeling of a complete powder coating process is carried out to understand the gas-particle two-phase flow field inside a powder coating booth and results of the numerical simulations are compared with experimental data to validate the numerical results. The flow inside the coating booth is modeled as a three-dimensional turbulent continuous gas flow with solid powder particles as a discrete phase. The continuous gas flow is predicted by solving Navier–Stokes equations using a standard k–ε turbulence model with non-equilibrium wall functions. The discrete phase is modeled based on a Lagrangian approach. In the calculation of particle propagation, a particle size distribution obtained through experiments is applied. The electrostatic field, including the effect of space charge due to free ions, is calculated with the use of the user defined scalar transport equations and user defined scalar functions in the software package, FLUENT, for the electrostatic potential and charge density.     

Lagrangian particle dispersion in turbulent flow over a wall mounted obstacle

Large-eddy simulations (LES) of particle-laden turbulent flows are presented in order to investigate the effects of particle response time on the dispersion patterns of a space developing flow with an obstruction, where solid particles are injected inside the wake of an obstacle [Vincont, J.Y., Simoens, S., Ayrault M., Wallace, J.M., 2000. Passive scalar dispersion in a turbulent boundary layer from a line source at the wall and downstream of an obstacle. J. Fluid Mech. 424, 127–167]. The numerical method is based on a fully explicit fractional step approach and finite-differences on Cartesian grids, using the immersed boundary method (IBM) to represent the existence of solid obstacles. Two different turbulence models have been tested, the classical Smagorinsky turbulence model and the filtered structure function model. The dispersed phase was modelled either by an Eulerian approach or a Lagrangian particle tracking scheme of solid particles with Stokes numbers in the range St = 0–25, assuming one-way coupling between the two phases. A very good agreement was observed between the Lagrangian and Eulerian approaches. The effect of particle size was found to significantly differentiate the dispersion pattern for the inhomogeneous flow over the obstacle. Although in homogeneous flows like particle-laden turbulent channels near-wall particle clustering increases monotonically with particle size, for the examined flow over an obstacle, preferential concentration effects were stronger only for an intermediate range of Stokes numbers.     

Numerical modelling of complex turbulent free‐surface flows with the SPH method: an overview

The gridless smoothed particle hydrodynamics (SPH) method is now commonly used in computational fluid dynamics (CFD) and appears to be promising in predicting complex free‐surface flows. However, increasing flow complexity requires appropriate approaches for taking account of turbulent effects, whereas some authors are still working without any turbulence closure in SPH. A review of recently developed turbulence models adapted to the SPH method is presented herein, from the simplistic point of view of a one‐equation model involving mixing length to more sophisticated (and thus realistic) models like explicit algebraic Reynolds stress models (EARSM) or large eddy simulation (LES). Each proposed model is tested and validated on the basis of schematic cases for which laboratory data, theoretical or numerical solutions are available in the general field of turbulent free‐surface incompressible flows (e.g. open‐channel flow and schematic dam break). They give satisfactory results, even though some progress should be made in the future in terms of free‐surface influence and wall conditions. Recommendations are given to SPH users to apply this method to the modelling of complex free‐surface turbulent flows. Copyright © 2006 John Wiley & Sons, Ltd.     

Thermo-mechanical modeling of turbulent heat transfer in gas–solid flows including particle collisions

A thermo-mechanical turbulence model is developed and used for predicting heat transfer in a gas–solid flow through a vertical pipe with constant wall heat flux. The new four-way interaction model makes use of the thermal k_θ–τ_θ equations, in addition to the hydrodynamic k–τ transport, and accounts for the particle–particle and particle–wall collisions through a Eulerian/Lagrangian formulation. The simulation results indicate that the level of thermal turbulence intensity and the heat transfer are strongly affected by the particle collisions. Inter-particle collisions attenuate the thermal turbulence intensity near the wall but somewhat amplify the temperature fluctuations in the pipe core region. The hydrodynamic-to-thermal times-scale ratio and the turbulent Prandtl number in the region near the wall increase due to the inter-particle collisions. The results also show that the use of a constant or the single-phase gas turbulent Prandtl number produces error in the thermal eddy diffusivity and thermal turbulent intensity fields. Simulation results also indicate that the inter-particle contact heat conduction during collision has no significant effect in the range of Reynolds number and particle diameter studied.     

明渠湍流边界层中颗粒的运动与分布

利用直接数值模拟、点球浸入边界法和颗粒离散元法相结合的方法, 模拟了颗粒在明渠湍流边界层中的运动, 并对颗粒的瞬时位置进行了Voronoi 分析, 定量研究了颗粒在湍流边界层中的运动和分布规律. 研究发现:颗粒的输运对湍流的统计特征有影响, 其运动与近壁区湍流拟序结构密切相关, 在"喷发"结构作用下被带离壁面, 在"扫掠" 结构和自身重力作用下回到壁面; 在湍流边界层中, 颗粒倾向于聚集在低流速带, 呈条带状分布;颗粒在大部分时间处于"簇"状态, 偶尔跳跃到"空" 状态, 但能够很快回到邻近低速区域.      

Experimental studies on particle behaviour and turbulence modification in horizontal channel flow with different wall roughness

Detailed measurements in a developed particle-laden horizontal channel flow (length 6 m, height 35 mm, the length is about 170 channel heights) are presented using phase-Doppler anemometry for simultaneous determination of air and particle velocity. The particles were spherical glass beads with mean diameters in the range of 60 µm-1 mm. The conveying velocity could be varied between about 10 m/s and 25 m/s, and the particle mass loading could reach values of about 2 (the mass loading is defined as the ratio of particle to gas phase mass flow rates), depending on particle size. For the first time, the degree of wall roughness could be modified by exchanging the wall plates. The influence of these parameters and the effect of inter-particle collisions on the profiles of particle mean and fluctuating velocities and the normalised concentration in the developed flow were examined. It was shown that wall roughness decreases the particle mean velocity and enhances fluctuating velocities due to irregular wall bouncing and an increase in wall collision frequency, i.e. reduction in mean free path. Thereby, the larger particles are mainly more uniformly distributed across the channel, and gravitational settling is reduced. Both components of the particle velocity fluctuation were reduced with increasing mass loading due to inter-particle collisions and the momentum loss involved. Moreover, the effect of the particles on the air flow and the turbulent fluctuations was studied on the basis of profiles in the developed flow and turbulence spectra determined for the streamwise velocity component. In addition to the effect of particle size and mass loading on turbulence modulation, the influence of wall roughness was analysed. It was clearly shown that increasing wall roughness also results in a stronger turbulence dissipation due to two-way coupling.