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

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

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.     

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.     

Particle statistics in a gas–solid coaxial strongly swirling flow: A direct numerical simulation

A direct numerical simulation of a strongly coaxial swirling particle-laden flow is conducted with reference to a previous experiment. The carrier phase is simulated as a coaxial swirling flow through a short nozzle injecting into a large container. The particle phase is carried by the primary jet, and simulated in the Lagrangian approach. The drag force, slip-shear force and slip-rotation force experienced by particles are calculated. A partial validation of the results is followed. The results are analyzed in Eulerian approach focusing on the statistical behavior of particle motion. The relative importance of the drag, slip-shear and slip-rotation forces under different Stokes numbers is indicated quantitatively. The particle velocity profiles, fluctuations, Reynolds stress, and turbulence intensity are demonstrated and analyzed respectively. An important “choke” behavior for large particles within the mainstream is found and interpreted. Additionally, the patterns of particle distribution and the helical structures of particle motion under different Stokes numbers are demonstrated qualitatively and analyzed quantitatively.     

土体颗粒物流动物质点法模拟的弹塑性和非牛顿流体本构模型比较研究

土体颗粒物流动是一种典型的大变形破坏,具有非牛顿流体的流动特征。准确模拟土体颗粒物的流动及冲击过程,对滑坡和泥石流等地质灾害的防治具有重要意义。物质点法是一种无网格粒子类方法,已在各类大变形问题中得到了广泛应用。以往土体颗粒物流动的模拟,通常采用弹塑性本构模型,但缺乏对非牛顿本构模型的模拟分析。本文引入非牛顿本构模型的模拟分析,旨在为土体颗粒物流动模拟提供一种新的方法与思路。非牛顿本构模型的模拟分析是将非牛顿广义Cross模型引入三维物质点法,通过人工阻尼力模拟颗粒间的摩擦力,对土体颗粒物的坍塌、沿斜面滑动以及冲击障碍物等问题进行了动态模拟,研究了其运动全过程,并与弹塑性本构模型的模拟结果进行了对比验证。结果表明,基于非牛顿流体本构模型的物质点法可以较好地模拟土体颗粒物加速、减速到再次稳定的流动全过程及其对障碍物的冲击效应。     

Two- and Four-Way Coupled Euler–Lagrangian Large-Eddy Simulation of Turbulent Particle-Laden Channel Flow

Large-eddy simulations (LES) of a vertical turbulent channel flow laden with a very large number of solid particles are performed. The motivation for this research is to get insight into fundamental aspects of co-current turbulent gas-particle flows, as encountered in riser reactors. The particle volume fraction equals about 1.3%, which is relatively high in the context of modern LES of two-phase flows. The channel flow simulations are based on large-eddy approximations of the compressible Navier–Stokes equations in a porous medium. The Euler–Lagrangian method is adopted, which means that for each individual particle an equation of motion is solved. The method incorporates four-way coupling, i.e., both the particle-fluid and particle–particle interactions are taken into account. The results are compared to single-phase channel flow in order to investigate the effect of the particles on turbulent statistics. The present results show that due to particle–fluid interactions the mean fluid profile is flattened and the boundary layer is thinner. Compared to single-phase turbulent flow, the streamwise turbulence intensity of the gas phase is increased, while the normal and spanwise turbulence intensities are reduced. This finding is generally consistent with existing experimental data. The four-way coupled simulations are also compared with two-way coupled simulations, in which the inelastic collisions between particles are neglected. The latter comparison clearly demonstrates that the collisions have a large influence on the main statistics of both phases. In addition, the four-way coupled simulations contain stronger coherent particle structures. It is thus essential to include the particle–particle interactions in numerical simulations of two-phase flow with volume fractions around one percent.     

Effect of temperature gradient within a solid particle on the rotation and oscillation modes in solid-dispersed two-phase flows

Liquid–solid two-phase flow with heat transfer is simulated, and the effect of temperature gradient within a solid particle on the particle behaviour and heat transfer is studied. The interaction between fluid and particles is considered with our original immersed solid approach on a rectangular grid system. The local heat flux at the fluid–solid interface is described with an anisotropic heat conductivity matrix, and the governing equation of temperature is time-updated with an implicit treatment for the diffusion term. The method is applied to a 2-D natural convection flow of a relatively low Rayleigh number including multiple particles. Heat transfer and particle behaviours are studied for different solid heat conductivities (ratio to the fluid conductivity ranging between 10⁻³ and 10³) and solid volume fractions. Under a condition of relatively low heat conductivity ratio, the particles show a simple circulating flow. By increasing the heat conductivity ratio, a transition of the particulate flow is observed to oscillation mode around the domain centre due to the buoyancy force as a restitution force. The oscillation period is found to vary with the heat conductivity ratio, and it is related to the time scales for the heat transfer via fluid and solid.     

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.     

Modeling fluid-particle interaction in dilute-phase turbulent liquid-particle flow simulation

The present work examines the predictive capability of a two-fluid CFD model that is based on the kinetic theory of granular flow in simulating dilute-phase turbulent liquid-particle pipe flows in which the inter-stitial fluid effect on the particle fluctuating motion is significant.The impacts of employing different drag correlations and turbulence closure models to describe the fluid-particle interactions(i.e.drag force and long-range interaction)are examined at both the mean and fluctuating velocity levels.The model pre-dictions are validated using experimental data of turbulent liquid-particle flows in a vertical pipe at different particle Reynolds numbers(ReP > 400 and ReP < 400),which characterize the importance of the vortex shedding phenomenon in the fluid-phase turbulence modulation.The results indicate that(1)the fluctuating velocity level predictions at different ReP are highly sensitive to the drag correlation selec-tion and(2)different turbulence closure models must be employed to accurately describe the long-range fluid-particle interaction in each phase.In general,good agreement is found between the model predic-tions and the experimental data at both the mean and fluctuating velocity levels provided that appropriate combinations of the drag correlation and the turbulence closure model are selected depending on Rep.     

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.     

气固两相流介尺度LBM-DEM模型

LBM-DEM耦合方法通常是指一种颗粒流体系统直接数值模拟算法,即是一种不引入经验曳力模型的计算方法,颗粒尺寸通常比计算网格的长度大一个量级,颗粒的受力通过表面的粘性力与压力积分获得,其优点是能描述每个颗粒周围的详细流场,产生详细的颗粒-流体相互作用的动力学信息,可以探索颗粒流体界面的流动、传递和反应的详细信息及两相相互作用的本构关系,但其缺点是计算量巨大,无法应用于真实流化床过程模拟。本文针对气固流化床中的流体以及固体颗粒间的多相流体力学行为,建立了一种稠密气固两相流的介尺度LBMDEM模型,即LBM-DEM耦合的离散颗粒模型,实现在颗粒尺度上流化床的快速离散模拟。该耦合模型采用格子玻尔兹曼方法(LBM)描述气相的流动和传递行为,离散单元法(DEM)用于描述颗粒相的运动,并利用能量最小多尺度(EMMS)曳力解决气固耦合不成熟问题,以提高其模拟精度。通过经典快速流态化的模拟,验证了介尺度LBM-DEM耦合模型的有效性。模拟结果表明介尺度LBM-DEM模型是一种探索实验室规模气固系统的有力手段。     

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.     

Experimental investigation on turbulence modification in a horizontal channel flow at relatively low mass loading

Particle-laden flows in a horizontal channel were investigated by means of a two-phase particle image velocimetry (PIV) technique. Experiments were performed at a Reynolds number of 6 826 and the flow is seeded with polythene beads of two sizes, 60 μm and 110 μm. One was slightly smaller than and the other was larger than the Kolmogorov length scale. The particle loadings were relatively low, with mass loading ratio ranging from 5×10⁻⁴ to 4×10⁻² and volume fractions from 6×10⁻⁷ to 4.8×10⁻⁵, respectively. The results show that the presence of particles can dramatically modify the turbulence even under the lowest mass loading ratio of 5×10⁻⁴. The mean flow is attenuated and decreased with increasing particle size and mass loading. The turbulence intensities are enhanced in all the cases concerned. With the increase of the mass loading, the intensities vary in a complicated manner in the case of small particles, indicating complicated particle-turbulence interactions; whereas they increase monotonously in the case of large particles. The particle velocities and concentrations are also given. The particles lag behind the fluid in the center region but lead in the wall region, and this trend is more prominent for the large particles. The streamwise particle fluctuations are larger than the gas fluctuations for both sizes of particles, however their varying trend with the mass loadings is not so clear. The wall-normal fluctuations increase with increasing mass loadings. They are smaller in the 60 μm particle case but larger in the 110 μm particle case than those of the gas phase. It seems that the small particles follow the fluid motion to certain extent while the larger particles are more likely dominated by their own inertia. Finally, remarkable non-uniform distributions of particle concentration are observed, especially for the large particles. The inertia of particles is proved to be very important for the turbulence modification and particles behaviors and thus should be considered in horizontal channels. The project supported by the National Natural Science Foundation of China (50276021), and Program for New Century Excellent Talents in University, Ministry of Education (NCET-04-0708) The English text was polished by Yunming Chen.     

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.     

Large eddy simulation of gas-particle turbulent channel flow with momentum exchange between the phases

This paper presents results of a large eddy simulation (LES) combined with Lagrangian particle tracking and a point-force approximation for the feedback effect of particles on the downward turbulent gaseous flow in a vertical channel. The LES predictions are compared with the results obtained by direct numerical simulation (DNS) of a finer computational mesh. A parametric study is conducted for particles with two response times in simulations with and without streamwise gravitational settling and elastic, binary interparticle collisions. It is shown that the classical and the dynamic Smagorinsky turbulence models adequately predict the particle-induced changes in the mean streamwise velocity and the Reynolds stresses of the carrier phase for the range of parameters studied. However, the largest discrepancies between the LES and DNS results are found in the cases of particle-laden flows. Conditional sampling of the instantaneous resolved flow fields indicates that the mechanisms by which particles directly oppose the production of momentum and vorticity of the organized fluid motions are also observed in the LES results. However, the geometric features of the near-wall quasistreamwise vortices are overestimated by the use of both turbulence models compared to the DNS predictions.     

The effect of two-way coupling and inter-particle collisions on turbulence modulation in a vertical channel flow

Turbulence modulation due to its interaction with dispersed solid particles in a downward fully developed channel flow was studied. The Eulerian framework was used for the gas-phase, whereas the Lagrangian approach was used for the particle-phase. The steady-state equations of conservation of mass and momentum were used for the gas-phase, and the effect of turbulence on the flow-field was included via the standard k–ε model. The particle equation of motion included the drag, the Saffman lift and the gravity forces. Turbulence dispersion effect on the particles was simulated as a continuous Gaussian random field. The effects of particles on the flow were modeled by appropriate source terms in the momentum, k and ε equations. Particle–particle collisions and particle–wall collisions were accounted for in these simulations. Gas-phase velocities and turbulence kinetic energy in the presence of 2–100% mass loadings of two particle classes (50 μm glass and 70 μm copper) were evaluated, and the results were compared with the available experimental data and earlier numerical results. The simulation results showed that when the inter-particle collisions were important and was included in the computational model, the fluid turbulence was attenuated. The level of turbulence attenuation increased with particle mass loading, particle Stokes number, and the distance from the wall. When the inter-particle collisions were negligible and/or was neglected in the model, the fluid turbulence was augmented for the range of particle sizes considered.     

Modulation of homogeneous turbulence seeded with finite size bubbles or particles

The dynamics of homogeneous, isotropic turbulence seeded with finite sized particles or bubbles is investigated in a series of numerical simulations, using the force-coupling method for the particle phase and low wavenumber forcing of the flow to sustain the turbulence. Results are given on the modulation of the turbulence due to massless bubbles, neutrally buoyant particles and inertial particles of specific density 1.4 at volumetric concentrations of 6%. Buoyancy forces due to gravity are excluded to emphasize finite size and inertial effects for the bubbles or particles and their interactions with the turbulence. Besides observing the classical entrapment of bubbles and the expulsion of inertial particles by vortex structures, we analyze the Lagrangian statistics for the velocity and acceleration of the dispersed phase. The turbulent fluctuations are damped at mid-range wavenumbers by the bubbles or particles while the small-scale kinetic energy is significantly enhanced. Unexpectedly, the modulation of turbulence depends only slightly on the dispersion characteristics (bubble entrapment in vortices or inertial sweeping of the solid particles) but is closely related to the stresslet component (finite size effect) of the flow disturbances. The pivoting wavenumber characterizing the transition from damped to enhanced energy content is shown to vary with the size of the bubbles or particles. The spectrum for the energy transfer by the particle phase is examined and the possibility of representing this, at large scales, through an additional effective viscosity is discussed.     

Assessment of force models on finite-sized particles at finite Reynolds numbers

Finite-sized inertial spherical particles are fully-resolved with the immersed boundary projection method(IBPM) in the turbulent open-channel flow by direct numerical simulation(DNS). The accuracy of the particle surface force models is investigated in comparison with the total force obtained via the fully-resolved method. The results show that the steady-state resistance only performs well in the streamwise direction, while the fluid acceleration force, the added-mass force, and the shear-induc...