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.     

液固两相槽道流中湍流调制的数值研究

结合雷诺应力模型 (Reynolds Stress Model, RSM) 和混合模型 (Mixture Model) 对槽道湍流向下流动中的颗粒调制湍流问题进行了研究．该模型考虑了颗粒流的动能理论和颗粒对湍流的反馈作用．着重分析了颗粒对湍流的调制作用,以及颗粒参数变量（如颗粒密度和质量载荷）对湍流调制的影响．结果表明：(1)在颗粒抑制湍流的范围内,当颗粒密度小于载流体密度时,湍流强度的改变量与颗粒密度成反比;当颗粒密度大于载流体密度时,湍流强度的改变量与颗粒密度成正比;(2)在一定范围内,颗粒抑制湍流的能力随颗粒质量载荷增加而变强．     

Unsteady Laminar to Turbulent Flow in a Spacer-Filled Channel

A combined numerical and experimental investigation has been carried out to study the flow behaviour in a spacer-filled channel, representative of those used in spiral-wound membrane modules. Direct numerical simulation and particle image velocimetry were used to investigate the fluid flow characteristics inside a 2 × 2 cell at Reynolds numbers that range between 100 and 1000. It was found that the flow in this geometry moves parallel to and also rotates between the spacer filaments and that the rate of rotation increases with Reynolds number. The flow mechanisms, transition process and onset of turbulence in a spacer-filled channel are investigated including the use of the velocity spectra at different Reynolds numbers. It is found that the flow is steady for Re < 200 and oscillatory at Re ～ 250 and increasingly unsteady with further increases in Re before the onset of turbulent flow at Re ～ 1000.     

Effects of Bubbles, Droplets or Particles on Heat Transfer in Turbulent Channel Flows

Adding small immiscible substances of low volume to turbulent channel water flow modifies turbulence and turbulent heat transfer. We examine the effects on the modifications of the ratios of the physical properties of substances to the physical properties of the carrier fluid, the particle Reynolds number, the Stokes number and the Weber number. We discuss the applicability and importance of the local, instantaneous values of these non-dimensional numbers. In particular, the maxima and minima for the time change in the Stokes number are found to correspond to the minima and maxima respectively for the time change of the wall-normal distance of a heavy, solid particle in turbulent flow in a vertical channel. We also investigate the effectiveness of the recent developments of particle-tracking velocimetry and particle-resolved direct numerical simulation, in increasing the understanding of changes in turbulence.     

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.     

Comparison of the RANS and PDF methods for air-particle flows

Two simulation methods, namely Reynolds-Averaged Navier–Stokes (RANS) equations, and Probability Distribution Function (PDF) are currently widely used for the modeling of multiphase flows. These two approaches are supplemented with appropriate closure equations that take into account all the pertinent forces and interaction effects on the solid particles, such as: particle–turbulence interactions; turbulence modulation; particle–particle interactions; particle–wall interactions; gravitation, drag and lift forces. The two methods have been used in order to simulate the turbulent particulate flow in upward pipes. The flow domain in all cases was a cylindrical pipe and the computations were carried for upward pipe flow. Monodisperse as well as polydisperse mixtures of particles have been considered. In general, the average velocity results obtained from the two methods are in close agreement, because the methods predict well the average velocity distribution of the carrier fluid as well as the solids. Thus, the differences in the average axial velocities predicted by the methods are not substantial. Differences in the turbulence intensity are more significant. A comparison of the numerical results obtained shows the relative importance of retaining the diffusion terms in both the axial and radial directions in the RANS method. Also the comparisons of the results show the relative effect of the lift forces in the distribution of solid particles.     

Lattice–Boltzmann simulations for analysing the detachment of micron-sized spherical particles from surfaces with large-scale roughness structures

Fully resolved numerical simulations of a micron-sized spherical particle residing on a surface with large-scale roughness are performed by using the Lattice–Boltzmann method. The aim is to investigate the influence of surface roughness on the detachment of fine drug particles from larger carrier particles for transporting fine drug particles in a DPI (dry powder inhaler). Often the carrier surface is modified by mechanical treatments for modifying the surface roughness in order to reduce the adhesion force of drug particles. Therefore, drug particle removal from the carrier surface is equivalent to the detachment of a sphere from a rough plane surface. Here a sphere with a diameter of 5 μm at a particle Reynolds number of 1.0, 3.5 and 10 are considered. The surface roughness is described as regularly spaced semi-cylindrical asperities (with the axes oriented normal to the flow direction) on a smooth surface. The influence of asperity distance and size ratio (i.e. the radius of the semi-cylinder to the particle radius, R_c/R_d) on particle adhesion and detachment are studied. The asperity distance is varied in the range 1.2 < L/R_d < 2 and the semi-cylinder radius between 0.5 < R_c/R_d < 0.75. The required particle resolution and domain size are appropriately selected based on numerical studies, and a parametric analysis is performed to investigate the relationship between the contact distance (i.e. half the distance between the particle contact points on two neighbouring semi-cylinders), the asperity distance, the size ratio, and the height of the particle centroid from the plane wall. The drag, lift and torque acting on the spherical particle are measured for different particle Reynolds numbers, asperity distances and sizes or diameters. The detachment of particles from rough surfaces can occur through lift-off, sliding and rolling, and the corresponding detachment models are constructed for the case of rough surfaces. These studies will be the basis for developing Lagrangian detachment models that eventually should allow the optimisation of dry powder inhaler performance through computational fluid dynamics.     

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.     

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.     

Experimental study on flow control of the turbulent boundary layer with micro-bubbles

The effect of micro-bubbles on the turbulent boundary layer in the channel flow with Reynolds numbers (Re) ranging from \(0.87\times 10 ^{5}\) to \(1.23\times 10^{5}\) is experimentally studied by using particle image velocimetry (PIV) measurements. The micro-bubbles are produced by water electrolysis. The velocity profiles, Reynolds stress and instantaneous structures of the boundary layer, with and without micro-bubbles, are measured and analyzed. The presence of micro-bubbles changes the streamwise mean velocity of the fluid and increases the wall shear stress. The results show that micro-bubbles have two effects, buoyancy and extrusion, which dominate the flow behavior of the mixed fluid in the turbulent boundary layer. The buoyancy effect leads to upward motion that drives the fluid motion in the same direction and, therefore, enhances the turbulence intense of the boundary layer. While for the extrusion effect, the presence of accumulated micro-bubbles pushes the flow structures in the turbulent boundary layer away from the near-wall region. The interaction between these two effects causes the vorticity structures and turbulence activity to be in the region far away from the wall. The buoyancy effect is dominant when the Re is relatively small, while the extrusion effect plays a more important role when Re rises.     

Mean and fluctuating components of drag and lift forces on an isolated finite-sized particle in turbulence

We perform fully resolved direct numerical simulations of an isolated particle subjected to free-stream turbulence in order to investigate the effect of turbulence on the drag and lift forces at the level of a single particle, following Bagchi and Balachandar’s work (Bagchi and Balachandar in Phys Fluids 15:3496–3513, 2003). The particle Reynolds numbers based on the mean relative particle velocity and the particle diameter are Re?=?100, 250 and 350, which covers three different regimes of wake evolution in a uniform flow: steady axisymmetric wake, steady planar symmetric wake, and unsteady planar symmetric vortex shedding. At each particle Reynolds number, the turbulent intensity is 5–10% of the mean relative particle velocity, and the corresponding diameter of the particle is comparable to or larger than the Kolmogorov scale. The simulation results show that standard drag values determined from uniform flow simulations can accurately predict the drag force if the turbulence intensity is sufficiently weak (5% or less compared to the mean relative velocity). However, it is shown that for finite-sized particles, flow non-uniformity, which is usually neglected in the case of the small particles, can play an important role in determining the forces as the relative turbulence intensity becomes large. The influence of flow non-uniformity on drag force could be qualitatively similar to the Faxen correction. In addition, finite-sized particles at sufficient Reynolds number are inherently subjected to stochastic forces arising from their self-induced vortex shedding in addition to lift force arising from the local ambient flow properties (vorticity and strain rate). The effect of rotational and strain rate of the ambient turbulence seen by the particle on the lift force is explored based on the conditional averaging using the generalized representation of the quasi-steady force proposed by Bagchi and Balachandar (J Fluid Mech 481:105–148, 2003). From the present study, it is shown that at Re?=?100, the lift force is mainly influenced by the surrounding turbulence, but at Re = 250 and 350, the lift force is affected by the wake structure as well as the surrounding turbulence. Thus, for a finite-sized particle of sufficient Reynolds number supporting self-induced vortex shedding, the lift force will not be completely correlated with the ambient flow. Therefore, it appears that in order to reliably predict the motion of a finite-sized particle in turbulence, it is important to incorporate both a deterministic component and a stochastic component in the force model. The best deterministic contribution is given by the conditional average. The influence of ambient turbulence at the scale of the particle, which are not accounted for in the deterministic contribution, can be considered in stochastic manner. In the modeling of lift force, additional stochastic contribution arising from self-induced vortex shedding must also be included.     

Effect of an unsteady swirled turbulent flow on the motion of a single solid particle

An unsteady swirled turbulent flow between two rotating flat disks is modeled. The flow is directed along the radius toward the rotation axis. A quasi-steady character of the turbulent flow, caused by oscillations of the radial velocity at the entrance to the gap between the disks with a period close to the time of dynamic relaxation of the particle, is studied with the use of the known two-equation Wilcox’s k-ω model of turbulence. The influence of the Stokes number and the frequency and amplitude of oscillations of the carrier medium on the motion of single particles in the field of centrifugal and aerodynamic forces is considered.     

The simulation of turbulent particle‐laden channel flow by the Lattice Boltzmann method

We perform direct numerical simulation of three‐dimensional turbulent flows in a rectangular channel, with a lattice Boltzmann method, efficiently implemented on heavily parallel general purpose graphical processor units. After validating the method for a single fluid, for standard boundary layer problems, we study changes in mean and turbulent properties of particle‐laden flows, as a function of particle size and concentration. The problem of physical interest for this application is the effect of water droplets on the turbulent properties of a high‐speed air flow, near a solid surface. To do so, we use a Lagrangian tracking approach for a large number of rigid spherical point particles, whose motion is forced by drag forces caused by the fluid flow; particle effects on the latter are in turn represented by distributed volume forces in the lattice Boltzmann method. Results suggest that, while mean flow properties are only slightly affected, unless a very large concentration of particles is used, the turbulent vortices present near the boundary are significantly damped and broken down by the turbulent motion of the heavy particles, and both turbulent Reynolds stresses and the production of turbulent kinetic energy are decreased because of the particle effects. We also find that the streamwise component of turbulent velocity fluctuations is increased, while the spanwise and wall‐normal components are decreased, as compared with the single fluid channel case. Additionally, the streamwise velocity of the carrier (air) phase is slightly reduced in the logarithmic boundary layer near the solid walls. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

Transport and deposition of neutral particles in magnetohydrodynamic turbulent channel flows at low magnetic Reynolds numbers

The effect of Lorentz force on particle transport and deposition is studied by using direct numerical simulation of turbulent channel flow of electrically conducting fluids combined with discrete particle simulation of the trajectories of uncharged, spherical particles. The magnetohydrodynamic equations for fluid flows at low magnetic Reynolds numbers are adopted. The particle motion is determined by the drag, added mass, and pressure gradient forces. Results are obtained for flows with particle ensembles of various densities and diameters in the presence of streamwise, wall-normal or spanwise magnetic fields. It is found that the particle dispersion in the wall-normal and spanwise directions is decreased due to the changes of the underlying fluid turbulence by the Lorentz force, while it is increased in the streamwise direction. The particle accumulation in the near-wall region is diminished in the magnetohydrodynamic flows. In addition, the tendency of small inertia particles to concentrate preferentially in the low-speed streaks near the walls is strengthened with increasing Hartmann number. The particle transport by turbophoretic drift and turbulent diffusion is damped by the magnetic field and, consequently, particle deposition is reduced.     

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

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

The effect of particles on fluid turbulence in a turbulent boundary layer over a cylinder

The turbulent fluid and particle interaction in the turbulent boundary layer for cross flow over a cylinder has been experimentally studied. A phase-Doppler anemometer was used to measure the mean and fluctuating velocities of both phases. Two size ranges of particles (30μm–60μm and 80μm–150μm) at certain concentrations were used for considering the effects of particle sizes on the mean velocity profiles and on the turbulent intensity levels. The measurements clearly demonstrated that the larger particles damped fluid turbulence. For the smaller particles, this damping effect was less noticeable. The measurements further showed a delay in the separation point for two phase turbulent cross flow over a cylinder. The project supported by the National Natural Science Foundation of China     

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.     

A Unifying Explanation for the Damping of Turbulence by Additives or External Forces

Recent calculations for the turbulent flow of a suspension of solid spheres in a gas were carried out by solving a Navier Stokes equation for the fluid that recognizes the presence of particles as external point forces. These show that a strong damping of the fluid turbulence can be realized at remarkably small volume fractions. The suggestion is made that the presence of point forces, pseudo-point forces, or added point shear stresses could provide a general explanation for turbulence suppression caused by additives or by the acceleration of a turbulent boundary-layer.