Influence of gravity and lift on particle velocity statistics and transfer rates in turbulent vertical channel flow

The present study reports detailed statistics for velocity and transfer rates of heavy particles dispersed in turbulent boundary layers. Statistics have been extracted from a homogeneous source of data covering a large target parameter space and are used here to analyze the effects of gravity and lift on particle dispersion and deposition in a systematic way. Datasets were obtained performing Direct Numerical Simulation (DNS) of particle-laden turbulent upward/downward flow in a vertical channel. Six values for the particle timescale (the particle Stokes number, St) ranging three orders of magnitude were considered to analyze the deposition process as the controlling mechanism was shifting from turbulent diffusion to inertia-moderated crossing trajectories. For the particle timescales examined, gravity and lift do not influence the qualitative behavior of particles even though velocity profiles and deposition coefficients are modified in a non-monotonic fashion, reaching an optimum for St ? 15. Physical mechanisms for the different behavior are discussed. Raw data and statistics obtained from the present DNS are made available at http://cfd.cineca.it (mirror site: http://158.110.32.35/download/database) and can be used to test simple models and closure equations for multiphase RANS and Large Eddy simulations.     

Direct numerical simulation of particle dispersion in a turbulent jet considering inter-particle collisions

To investigate the behaviour of inter-particle collision and its effects on particle dispersion, direct numerical simulation of a three-dimensional two-phase turbulent jet was conducted. The finite volume method and the fractional-step projection algorithm were used to solve the governing equations of the gas phase fluid and the Lagrangian method was applied to trace the particles. The deterministic hard-sphere model was used to describe the inter-particle collision. In order to allow an analysis of inter-particle collisions independent of the effect of particles on the flow, two-way coupling was neglected. The inter-particle collision occurs frequently in the local regions with higher particle concentration of the flow field. Under the influence of the local accumulation and the turbulent transport effects, the variation of the average inter-particle collision number with the Stokes number takes on a complex non-linear relationship. The particle distribution is more uniform as a result of inter-particle collisions, and the lateral and the spanwise dispersion of the particles considering inter-particle collision also increase. Furthermore, for the case of particles with the Rosin–Rammler distribution (the medial particle size is set d₅₀ = 36.7 μm), the collision number is significantly larger than that of the particles at the Stokes number of 10, and their effects on calculated results are also more significant.     

Effect of Particle Clusters on Carrier Flow Turbulence: A Direct Numerical Simulation Study

Experiments indicate that particle clusters that form in fluidized–bed risers can enhance gas-phase velocity fluctuations. Direct numerical simulations (DNS) of turbulent flow past uniform and clustered configurations of fixed particle assemblies at the same solid volume fraction are performed to gain insight into particle clustering effects on gas-phase turbulence, and to guide model development. The DNS approach is based on a discrete-time, direct-forcing immersed boundary method (IBM) that imposes no-slip and no-penetration boundary conditions on each particle’s surface. Results are reported for mean flow Reynolds number Re _p ?=?50 and the ratio of the particle diameter d _p to Kolmogorov scale is 5.5. The DNS confirm experimental observations that the clustered configurations enhance the level of fluid-phase turbulent kinetic energy (TKE) more than the uniform configurations, and this increase is found to arise from a lower dissipation rate in the clustered particle configuration. The simulations also reveal that the particle-fluid interaction results in significantly anisotropic fluid-phase turbulence, the source of which is traced to the anisotropic nature of the interphase TKE transfer and dissipation tensors. This study indicates that when particles are larger than the Kolmogorov scale (d _p ?>?η), modeling the fluid-phase TKE alone may not be adequate to capture the underlying physics in multiphase turbulence because the Reynolds stress is anisotropic. It also shows that multiphase turbulence models should consider the effect of particle clustering in the dissipation model.     

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.     

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

对 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 附近出现另一个峰值.     

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.     

Temporal stability of a particle-laden jet

The temporal instability of a particle-laden jet was investigated numerically which took into consideration the parametric effects of jet parameter, B, jet Reynolds number, Re_j, particle mass loading, Z and Stokes number, St. The linear stability theory was used to derive the instability equations of a viscous particle-laden jet flow. The single-phase instability of a top-hat jet was then calculated and compared with the available analytical theories. The numerical results agree well with the analytical results for both the axisymmetric (n = 0) and first azimuthal (n = 1) modes. The results show that the first azimuthal mode disturbance is usually more unstable than that of the axisymmetric mode. But the axisymmetric mode disturbance can be more unstable when Z is high enough (i.e., Z ? 0.1). The higher B and Re_j are, the more unstable the particle-laden jet will be. The existence of particles enhances the flow stability. With the increasing of Z, the jet flow will grow more stable. The inviscid single-phase jet is the most unstable. The wave amplification, c_i first decreases with the increasing of St and then increases afterwards. There exist certain values of St, at which the jet is the most stable.     

Turbulent transport of particles in a straight square duct

Turbulent flow through a duct of square cross-section gives rise to off-axis secondary flows, which are known to transfer momentum between fluid layers thereby flattening the velocity profile. The aim of this study is to investigate the role of the secondary flows in the transport and dispersion of particles suspended in a turbulent square duct flow. We have numerically simulated a flow through a square duct having a Reynolds number of Re_τ = 300 through discretization of the Navier–Stokes equations, and followed the trajectories of a large number of passive tracers and finite-inertia particles under a one-way coupling assumption. Snapshots of particle locations and statistics of single-particle and particle pair dispersion were analyzed. It was found that lateral mixing is enhanced for passive tracers and low-inertia particles due to the lateral advective transport that is absent in straight pipe and channels flows. Higher inertia particles accumulate close to the wall, and thus tend to mix more efficiently in the streamwise direction since a large number of the particles spend more time in a region where the mean fluid velocity is small compared to the bulk. Passive tracers tend to remain within the secondary swirling flows, circulating between the core and boundary of the duct.     

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.     

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.     

Numerical study of gas–solid two-phase flow and erosion in a cavity with a slope

Gate valve is mainly used to turn on or turn off the pipeline in pneumatic conveying. When the gate valve is fully open, the particles are easy to collide with the cavity rear wall and enter into the cavity, resulting in particles’ accumulation in the cavity. The particles in cavity will accumulate between the cavity bottom and the flashboard bottom wall and prevent the gate from turning off normally. Meanwhile, the particles’ collision with cavity rear wall will cause serious erosion. Both the particles’ accumulation and erosion will cause the poor sealing of the gate valve, further resulting in the leakage of the pipeline system. To reduce the particles’ accumulation in cavity and erosion on cavity when the gate valve is fully open, we simplify the gate valve into a cavity structure and study it. We find that adding a slope upstream the cavity can effectively reduce the particles’ accumulation in the cavity and the erosion on the cavity rear wall. In this work, Eulerian–Lagrangian method in commercial code (FLUENT) was used to study the gas–solid two-phase flow and erosion characteristics of a cavity with a slope. The particle distribution shows that the particles with Stokes number St = 1.3 and St = 13 cannot enter the cavity due to the slope, but the particles with St = 0.13 enter the cavity following the gas. For St = 13, the particles collide with the wall many times in the ideal cavity. Erosion results show that the slope can transfer the erosion on cavity rear wall to the slope and reduce the maximum erosion rate of the wall near the cavity to some degrees.     

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.     

Bubble migration in a turbulent boundary layer

The wall void peaking distribution observed in an upward turbulent bubbly boundary layer along a flat plate is generated by bubbles that move towards the plate, come into contact with the wall and then slide along it. This transverse ‘migration’ has been studied using flow visualization, high speed video and particle tracking techniques to measure the trajectories of mono-disperse air bubbles at very low void fractions. Investigations have been performed at four Reynolds numbers in the range 280 < Re_θ < 3000, covering both the laminar and turbulent regimes, with mono-disperse bubbles of mean equivalent diameter between 2 mm and 6 mm. Lagrangian statistics calculated from hundreds of trajectories show that the migration only occurs in the turbulent regime and for bubble diameters below some critical value: 3.5 mm < d_eqcrit < 4 mm. Above this size (We > 3), the interface deformation is such that bubbles do not remain at the wall, even when they are released at the surface. Also, bubble migration is shown to be non-systematic, to have a non-deterministic character in the sense that trajectories differ significantly, to increase with Reynolds number and to take place on a short time scale. A series of experiments with isolated bubbles demonstrates that these results are not influenced by bubble–bubble interactions and confirm that two-way coupling in the flow is limited. Flow visualizations show that the migration originates with the capture of bubbles inside the large turbulent structures of the boundary layer (‘bulges’). The bubbles begin to move towards the wall as they cross these structures, and the point at which they reach the wall is strongly correlated with the position of the deep ‘valleys’ which separate the turbulent ‘bulges’. The analysis of the mean Lagrangian trajectories of migrating bubbles confirms these observations. Firstly, the average time of migration calculated from these trajectories coincides with the mean transit time of the bubbles across the structures. Secondly, once the trajectories have been scaled by this transit time and the boundary layer thickness δ, they all have the same form in the region y/δ < 0.4, independent of the Reynolds number.     

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.     

Experimental investigation of heat transfer in vertical upward and downward supercritical CO2 flow in a circular tube

An experimental investigation of turbulent heat transfer in vertical upward and downward supercritical CO₂ flow was conducted in a circular tube with an inner diameter of 4.5 mm. The experiments were performed for bulk fluid temperatures from 29 to 115 °C, pressures from 74.6 to 102.6 bar, local wall heat fluxes from 38 to 234 kW/m², and mass fluxes from 208 to 874 kg/m² s. At a moderate wall heat flux and low mass flux, the wall temperature had a noticeable peak value for vertical upward flow, but increased monotonically along the flow direction without a peak value for downward flow. The ratios of the experimental Nusselt number to the value obtained from a reference correlation were compared with Bo^* and q⁺ distributions to observe the buoyancy and flow-acceleration effects on heat transfer. In the experimental range of this study, the flow acceleration predominantly affected the heat-transfer phenomena. Based on an analysis of the shear-stress distribution in the turbulent boundary layer and the significant variation of the specific heat across the turbulent boundary layer, a new heat-transfer correlation for vertical upward and downward flow of supercritical pressurized fluid was developed; this correlation agreed with various experimental datasets within ±30%.     

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.     

SECOND-ORDER MOMENT MODEL FOR DENSE TWO-PHASE TURBULENT FLOW OF BINGHAM FLUID WITH PARTICLES

The USM-θmodel of Bingham fluid for dense two-phase turbulent flow was developed, which combines the second-order moment model for two-phase turbulence with the particle kinetic theory for the inter-particle collision. In this model, phases interaction and the extra term of Bingham fluid yield stress are taken into account. An algorithm for USM-θmodel in dense two-phase flow was proposed, in which the influence of particle volume fraction is accounted for. This model was used to simulate turbulent flow of Bingham fluid single-phase and dense liquid-particle two-phase in pipe. It is shown USM-θmodel has better prediction result than the five-equation model, in which the particle-particle collision is modeled by the particle kinetic theory, while the turbulence of both phase is simulated by the two-equation turbulence model. The USM-θmodel was then used to simulate the dense two-phase turbulent up flow of Bingham fluid with particles. With the increasing of the yield stress, the velocities of Bingham and particle decrease near the pipe centre. Comparing the two-phase flow of Bingham-particle with that of liquid-particle, it is found the source term of yield stress has significant effect on flow.     

An analytical approach for the closure equations of gas–solid flows with inter-particle collisions

This work examines the effect of inter-particle collisions on the motion of solid particles in two-phase turbulent pipe and channel flows. Two mechanisms for the particle–particle collisions are considered, with and without friction sliding. Based on these collision mechanisms, the correlations of the various velocity components of colliding particles are obtained analytically by using an averaging procedure. This takes into account three collision coordinates, two angles and the distance between the centers of colliding particles. The various stress tensor components are obtained and then introduced in the mass, linear momentum and angular momentum equations of the dispersed phase. The current approach applies to particle–particle collisions that result from both the average velocity difference and the turbulent velocity fluctuations. In order to close the governing equations of the dispersed phase, the pseudo-viscosity coefficients are defined and determined by the time of duration of the inter-particle collision process. The model is general enough to apply to both polydisperse and monodisperse particulate systems and has been validated by comparisons with experimental data.