90°弯管内气固两相湍流流动的数值研究

利用两相湍流KET模型对90°弯管内气固两相湍流流动进行了数值模拟,得到了弯管内两相流动的一些规律,并提出用颗粒动理学压力来定性表征弯管内磨损严重部位,为管道抗磨损设计提供了一定的理论依据。     

Direct numerical simulation of particle-laden plane turbulent wall jet and the influence of Stokes number

The distribution and motion of inertial particles in plane turbulent wall jet are investigated using direct numerical simulation, under the assumption of one-way coupling. To our knowledge, this appears to be the first direct numerical simulation of a particle-laden plane turbulent wall jet. It is shown that, in outer part of the wall jet, the behaviour of particles closely resembles that of a free plane jet. Due to the streamwise decay of particle Stokes number, the particle streaks formed in the near wall region of the wall jet are characterized by their intensity variation, which differs significantly from those in the channel flow. The streamwise growth of the particle velocity half-width is approximately equal to that of the fluid velocity half-width and the maximum velocity of particles decays slower than that of fluid due to inertia. The outer scaling can collapse the mean particle velocity in both the inner and outer region for heavier particles. In the buffer region, the particle–fluid velocity difference can be negative or positive depending on the Stokes number since there are two competing effects, namely the memory effect and turbophoresis. In the viscous region, the larger particles are on average faster than fluid and the velocity difference is found to be self-similar depending on outer Stokes number. The near-wall distribution of velocity difference is significantly correlated with the presence of high-momentum particles which are entrained by vortical structures generated in the outer region of the wall jet. These results are useful for environmental and engineering applications.     

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.     

Effects of the particle Stokes number on wind turbine airfoil erosion

Under natural conditions, wind turbines are inevitably eroded by the action of sand-wind flow. To further investigate the effects of dust drift on the erosion of the wind turbine blades in sand-wind environments, the effects of the wind velocity, particle diameter, and particle density on the erosion of wind turbine airfoils are studied, and the effects of the particle Stokes number on the airfoil erosion are discussed. The results show that, when the angle of attack(AOA) is 6.1°, there will be no erosion on the airfoil surface if the particle Stokes number is lower than 0.013 5, whereas erosion will occur if the particle Stokes number is higher than 0.015 1. Therefore, there exists a critical range for the particle Stokes number. When the particle Stokes number is higher than the maximum value in the critical range, airfoil erosion will occur. The result is further confirmed by changing the particle diameter, particle density, and inflow speed. It is shown that the erosion area on the airfoil and the maximum erosion rate are almost equal under the same particle Stokes number and AOA. The extent of airfoil erosion increases when the particle Stokes number increases, and the critical particle Stokes number increases when the AOA increases. Moreover, the geometric shape of the airfoil pressure surface greatly affects the airfoil erosion, especially at the curvature near the leading edge.     

Simulation of particles released near the wall in a turbulent boundary layer

A direct numerical simulation was used along with a Lagrangian particle tracking technique to study particle motion in a horizontal, spatially developing turbulent boundary layer along an upper-wall (with terminal velocity directed away from the wall). The objective of the research was to study particle diffusion, dispersion, reflection, and mean velocity in the context of two parametric studies: one investigated the effect of the drift parameter (the ratio of particle terminal velocity to fluid friction velocity) for a fixed and finite particle inertia, and the second varied the drift parameter and particle inertia by the same amount (i.e. for a constant Froude number). A range of drift parameters from 10⁻⁴ to 10⁰ were considered for both cases. The particles were injected into the simulation at a height of four wall units for several evenly distributed points across the span and a perfectly elastic wall collision was specified at one wall unit.Statistics collected along the particle trajectories demonstrated a transition in particle movement from one that is dominated by diffusion to one that is dominated by gravity. For small and intermediate sized particles (i.e. ones with outer Stokes numbers and drift parameters much less than unity) transverse diffusion away from the wall dominated particle motion. However, preferential concentration is seen near the wall for intermediate-sized particles due to inhomogeneous turbulence effects (turbophoresis), consistent with previous channel flow studies. Particle–wall collision statistics indicated that impact velocities tended to increase with increasing terminal velocity for small and moderate inertias, after which initial conditions become important. Finally, high relative velocity fluctuations (compared to terminal velocity) were found as particle inertia increased, and were well described with a quasi-one-dimensional fluctuation model.     

Statistics of Particle Accumulation in Spatially Developing Turbulent Boundary Layers

We present the results of a Direct Numerical Simulation of a particle-laden spatially developing turbulent boundary layer up to Re _θ ?=?2500. Two main features differentiate the behavior of inertial particles in a zero-pressure-gradient turbulent boundary layer from the more commonly studied case of a parallel channel flow. The first is the variation along the streamwise direction of the local dimensionless parameters defining the fluid-particle interactions. The second is the coexistence of an irrotational free-stream and a near-wall rotational turbulent flow. As concerns the first issue, an inner and an outer Stokes number can be defined using inner and outer flow units. The inner Stokes number governs the near-wall behavior similarly to the case of channel flow. To understand the effect of a laminar-turbulent interface, we examine the behavior of particles initially released in the free stream and show that they present a distinct behavior with respect to those directly injected inside the boundary layer. A region of minimum concentration occurs inside the turbulent boundary layer at about one displacement thickness from the wall. Its formation is due to the competition between two transport mechanisms: a relatively slow turbulent diffusion towards the buffer layer and a fast turbophoretic drift towards the wall.     

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.     

Numerical investigation of particle lateral migration in straight channel flows using a direct-forcing immersed boundary method

Inertia-induced cross-stream migration has been recently exploited for precise position of particles in confined channel flows. In this work, a three-dimensional finite volume based immersed boundary method has been developed to study the lateral migration and hydrodynamic self-assembly of neutrally-buoyant particles in pressure-driven flows. Simulation results show that, in 2D channel flows, the equilibrium position for a circular particle is closer to the centreline for larger particle Reynolds number due to the increasing flow rate, while in 3D square duct flow, the equilibrium position for a spherical particle is near a face centre and is closer to the wall for larger particle Reynolds number. Self-assembly of a pair of particles is observed in 3D square duct flows but not in 2D channel flows. Mechanisms for the self-assembly are discussed.     

环形截面螺旋管道内二次流动特性的研究

从曲线柱坐标系下的N-S方程出发,以曲率和挠率为小参数,采用摄动法求解了环形截面螺旋管道内的黏性流动,给出了完全二阶摄动解,结果表明：当挠率为零时,二次流表现为上下对称的四个涡;当挠率不为零,涡的对称性遭到破坏,二次涡的强度和个数受De数和环形截面内外径之比δ的影响,轴向速度最大值在De数较小时靠近管道的内侧,随着De数的增加,其最大值向外侧移动。     

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.     

Direct numerical simulation of particle behaviour in the wall region of turbulent flows in horizontal channels

The motion of small particles in the wall region of turbulent channel flows has been investigated using direct numerical simulation. It is assumed that the particle concentration is low enough to allow the use of one-way coupling in the calculations, i.e. the fluid moves the particles but there is no feedback from the particles on the fluid motion. The velocity of the fluid is calculated by using a pseudospectral, direct solution of the Navier-Stokes equations. The calculations indicate that particles tend to segregate into the low-speed regions of the fluid motion near the wall. The segregation tendency depends on the time constant of the particle made non-dimensional with the wall shear velocity and kinematic viscosity. For very small and very large time constants, the particles are distributed more uniformly. For intermediate time constants (of the order 3), the segregation into the low-speed fluid regions is the highest. The finding that segregation occurs for a range of particle time constants is supported by experimental results. The findings regarding the more uniform distributions, however, still remain to be verified against experimental data which is not yet available. For horizontal channel flows, it is also found that particles are resuspended by ejections (of portions of the low-speed streaks) from the wall and are, therefore, primarily associated with low-speed fluid. The smaller particles are flung further upwards and, as they fall back towards the wall, they tend to be accelerated close to the fluid velocity. The larger particles have greater inertia and, consequently, accelerate to lower velocities giving higher relative velocities. This velocity difference, as a function of wall-normal distance, follows the same trend as in experiments but is always somewhat smaller in the calculations. This appears to be due to the Reynolds number for the numerical simulation being smaller than that in the experiment. It is concluded that the average particle velocity depends not only on the wall variables for scaling, but also on outer variables associated with the mean fluid velocity and fluid depth in the channel. This is because fluid depth in combination with the wall shear velocity determines how much time a particle, of a given size and density, spends in the outer flow and, hence, how close it gets to the local fluid velocity.     

Towards high-fidelity erosion prediction: On time-accurate particle tracking in turbomachinery

Erosion and fouling caused by ingested particles causes performance degradation and safety issues in turbomachinery components. Simulating these processes is a complex multiphysics and multiscale problem which has not reached a satisfactory level of maturity yet. The current state of the art approach is based on RANS solutions, which provide an averaged carrier phase on which the particles are advanced in an a posteriori manner. Upon wall impact, the particle quantities are then fed to the erosion and rebound models. In this work, we present as an alternative to this approach an Eulerian/Lagrangian simulation framework of high-order accuracy in space and time for the time-resolved prediction of particle motion in complex flows. We apply it to the LES of a particle-laden flow over a T106C low pressure turbine linear cascade. We then contrast time-averaged and time-accurate flow fields as carrier phases for the particles and highlight how the associated modeling assumptions influence the solution. Based on the particle Stokes number, we identify characteristic regimes and their interaction with the flow phase. By a detailed comparison of the particle statistics, we highlight the effects of turbulent small scale behavior and define the modeling challenges associated with finding accurate particle closure models for time-averaged simulations. This framework constitutes a first step towards high-fidelity erosion prediction for turbomachinery applications.     

Numerical study of particle motion near a charged collector

The behavior of particles impacting the surface of a charged droplet involves adhesion, rebound, and submersion. In the present study, a numerical model for simulating particle impacts on charged droplets is presented that takes into account the various impact modes. With the droplet considered as a solid boundary, the criterion for rebounding is that the particle’s impact angle is <85°. The simulated trajectories of the particles are verified by comparing with experimental data for low-velocity particles to assess the reliability of the model. For impact angles >85°, particles undergo three distinct modes depending on normal impact velocities. The critical velocity of adhesion/rebound and rebound/submersion is used to identify the mode that the particles are undergoing. The criteria are also verified by comparing with analytical data. The results show that the impact angle of particles increases with increasing Coulomb number and decreases dramatically with increasing Stokes number, both of which lead to a high probability for particle rebound.     

Pressure drop characteristics in a rotating smooth square channel with a sharp 180° bend

The results of an experimental study to investigate the local pressure drop characteristics in a square cross-sectioned smooth channel with a sharp 180° bend rotating about an axis normal to the free-stream direction are reported here. The sharp 180° turn was obtained by dividing a rectangular passage into two channels using a divider wall with a rounded tip at the location where the flow negotiates the turn. The study was carried out for three ratios of divider wall thickness to hydraulic diameter (W/D), namely, 0.24, 0.37 and 0.73 all with a rounded tip divider wall and only for a bend with a W/D ratio of 0.37, the influence of a sharp tip divider wall was studied. Experiments were conducted for a Reynolds number varying from 10 000 to 17 000 with the rotation number (ωD/V) varying from 0 to 0.38. The pressure drop distribution, normalized with the mainstream fluid dynamic pressure head, is presented for the leading, trailing and the outer surfaces. The results indicate that the local pressure drop characteristics in the bend region are significantly affected by a change in the rotation number but the influence of the Reynolds number is weak. The friction factor is less sensitive to rotation for the bend with a W/D ratio of 0.24 when compared to bends with W/D ratios of 0.37 and 0.73. Friction factor correlations are presented which fit the experimental data within 10% for the range of parameters studied.     

Evaluation of erosion in a two‐way coupled fluid–particle system

In this study, the effects of flow turbulence intensity, temperature, particle sizes and impinging velocity on erosion by particle impact are demonstrated numerically. Underlying turbulent flow on an Eulerian frame is described by the compressible Reynolds averaged Navier–Stokes equations with a RNG k–ε turbulence model. The particle trajectories and particle–wall interactions are evaluated by a Eulerian–Lagrangian approach in a two‐way coupling system. An erosion model considering material weight removal from surfaces is used to predict erosive wear. Computational validation against measured data is demonstrated satisfactorily. The analysis of erosion shows that the prevention of erosion is enhanced by increasing the effects of flow temperature and turbulence intensity and reducing particle inertial momentum. Copyright © 2001 John Wiley & Sons, Ltd.     

Cross-stream forces and velocities of fixed and freely suspended particles in viscoelastic Poiseuille flow: Perturbation and numerical analyses

The cross-stream migration of a circular particles (or infinitely long cylinder) in two dimensional, inertia-less viscoelastic pressure-driven flows is examined through complementary finite element simulations and second-order fluid perturbation analyses for small Deborah number (De), where De is defined as the fluid relaxation time divided by the characteristic flow time. A neutrally buoyant, freely suspended particle migrates toward the center of the channel for all particle sizes and cross-stream positions due to the coupled effects of the linear and quadratic variations of the imposed velocity. A particle that is held at a fixed position, in contrast, experiences a cross-stream force directed toward the wall as a result of the coupled effects of the local shear flow and the flow relative to the particle.     

Turbulent heat-transfer characteristics along the heated convex wall of a rectangular cross-sectional return bend

The turbulent heat-transfer characteristics along the heated convex wall of a return bend which has rectangular cross section with large ratio have been examined for various clearances of the duct in detail. The experiments are performed under condition that the convex wall is heated at uniform heat flux while the concave wall is insulated. Water as a working fluid is utilized. Using four kinds of clearances of 15, 40, 60 and 80 mm, the Reynolds number in the turbulent range is varied from 8×10³ to 8×10⁴ with Prandtl number ranging from 6.5 to 8.5. In consequence, it is found that both the local and the mean heat-transfer rates are always smaller than those for straight parallel plates or for the straight duct. It is also found that the local heat-transfer characteristics in the outlet region of the return bend are more sensitively influenced by the variation of duct clearance than those in the inlet region.     

Hydrodynamic force and torque models for a particle moving near a wall at finite particle Reynolds numbers

This research work is aimed at proposing models for the hydrodynamic force and torque experienced by a spherical particle moving near a solid wall in a viscous fluid at finite particle Reynolds numbers. Conventional lubrication theory was developed based on the theory of Stokes flow around the particle at vanishing particle Reynolds number. In order to account for the effects of finite particle Reynolds number on the models for hydrodynamic force and torque near a wall, we use four types of simple motions at different particle Reynolds numbers. Using the lattice Boltzmann method and considering the moving boundary conditions, we fully resolve the flow field near the particle and obtain the models for hydrodynamic force and torque as functions of particle Reynolds number and the dimensionless gap between the particle and the wall. The resolution is up to 50 grids per particle diameter. After comparing numerical results of the coefficients with conventional results based on Stokes flow, we propose new models for hydrodynamic force and torque at different particle Reynolds numbers. It is shown that the particle Reynolds number has a significant impact on the models for hydrodynamic force and torque. Furthermore, the models are validated against general motions of a particle and available modeling results from literature. The proposed models could be used as sub-grid scale models where the flows between particle and wall can not be fully resolved, or be used in Lagrangian simulations of particle-laden flows when particles are close to a wall instead of the currently used models for an isolated particle.     

Particle transport in Poiseuille flow in narrow channels

Particle-tracking experiments were performed to validate a model [Staben, M.E., Zinchenko, A.Z., Davis, R.H., 2003. Motion of a particle between two parallel plane walls in low-Reynolds-number Poiseuille flow. Phys. Fluids 15, 1711–1733] for neutrally buoyant spherical particles convected by a Poiseuille flow in a thin microchannel for particles as large as d_p/H = 0.95, where d_p is the particle diameter and H is the channel width (narrow dimension). The measured and predicted velocities agree within experimental error and show that a particle’s velocity is more retarded when it is larger and/or closer to a channel wall. The particle distribution across the channel for a blunt entrance shows a focusing of small particles away from the walls and towards the center of the channel, whereas the particle distribution for an offset-angled entrance is slightly skewed towards the wall encountered first in the entrance region. As a result, the average particle velocities for the blunt entrance exceed those for the angled entrance. Moreover, due to the depletion of particles from the slow-moving region within one radius of the wall, the average particle velocity exceeds the average fluid velocity unless the particle diameter exceeds about 80% of the channel width.     

Solids behavior in dilute zone of a CFB riser under turbulent conditions

The behavior of the solid phase in the upper zone of a circulating fluidized bed riser was studied using a phase Doppler anemometer. Glass particles of mean diameter 107 μm and superficial gas velocities U_g covering the turbulent and the beginning of the fast fluidization regime were investigated. Three static bed heights were tested. Ascending and descending particles were found co-existing under all operating conditions tested, and at all measurement locations. Superficial gas velocity proved/happened to have a larger effect on descending particles at the wall and on ascending particles in the central region. Transversal particle velocities in both directions (toward the center and toward the wall) behaved relatively equivalently, with only slight difference observed at the wall. However, observation of the number of particles moving in either transversal direction showed a change in bed structure when increasing U_g. Furthermore, a balance was constantly observed between the core zone and the annulus zone where the mutual mass transfer between these two zones occurred continuously. Transition from a slow to a fast particle motion was accompanied by a transition to high levels of velocity fluctuations, and was found corresponding to the appearance of significant solid particle flow rate.