Numerical study of elbow erosion due to sand particles under annular flow considering liquid entrainment

Sand particle erosion is always a challenge in natural gas production. In particular, the erosion in gas–liquid–solid annular flow is more complicated. In this study, a three-phase flow numerical model that couples the volume of fluid multiphase flow model and the discrete phase model was developed for prediction of erosion in annular flow. The ability of the numerical model to simulate the gas–liquid annular flow is validated through comparison with the experimental data. On the basis of the above numerical model, the phase distribution in the pipe was analyzed. The liquid entrainment behavior was reasonably simulated through the numerical model, which guaranteed the accuracy of predicting the particle erosion. Additionally, four erosion prediction models were used for the erosion calculation, among them, the Zhang et al. erosion model predicted the realistic results. Through the analysis of the particle trajectory and the particle impact behavior on the elbow, the cushion effect of the liquid film on the particles and the erosion morphology generation at the elbow were revealed.     

Determination of the Parameters of a Radial Turbine Operating on Freon in the Quasi-One-Dimensional Approximation

The study contains a mathematical model for calculating the main parameters of a gas flow in the working channel of a turbine rotor in the quasi-one-dimensional approximation, the formulation of the problemof determining the geometric parameters of the turbine, and the results of calculations intended for an analysis of peculiarities occurring in modeling freon flows in the working channel of this radial turbine.     

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 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.     

Inertial-microfluidic radial migration in solid/liquid two-phase flow through a microcapillary: Particle equilibrium position

Although equilibrium of spherical particles under radial migration has been extensively investigated, mostly in macroscale flows with characteristic lengths on the order of centimeters, it is not fully characterized at relatively small Reynolds numbers, 1 ≤ Re ≤ 100. This paper experimentally studies “inertial microfluidic” radial migration of spherical particles in circular Poiseuille flow through a microcapillary. Microparticle tracking experiments are performed to obtain the spatial distribution of the particles by adopting a depth-resolved measurement technique. Through the analysis of the radial distribution of particles, inertial microfluidic circular Poiseuille flow is shown to induce a strong radial migration of particles at substantially small Re, which is quite in contrast to the pipe flows at large Re previously reported. This particle migration phenomenon is so prominent that particle equilibrium positions are formed even at small Re. However, it turns out that there exists a certain critical Re below which particle equilibrium position is almost fixed, but above which it seems to drift toward the channel wall.     

Gas turbulence modulation in the pneumatic conveying of massive particles in vertical tubes

The present work is concerned with the interaction between large particles and gas phase turbulence. Gas turbulence modulation in these systems is considered to be dominated by a generation mechanism which arises due to the presence of wakes behind particles. Following a recent proposal, a closure for gas turbulence modulation accounting for the effect of wakes is employed within the context of a mathematical model for particle-laden, turbulent flows. The model accounts for particle particle and particle-wall interactions associated with larger particles based on concepts from gas kinetic theory. It is shown that due to the significant flattening of the mean gas velocity profile with the addition of particles, and the corresponding decrease in turbulent energy production, a generation mechanism must be present in order to produce gas velocity fluctuation predictions which are consistent with the experimental measurements, even in the case where the experimental results indicate a net suppression of gas phase turbulence in the presence of particles.     

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.     

Dual characters of interaction between particles and fluid

The unique characteristics of gas-solids two-phase flow and fluidization in terms of the flow structures and the apparent behavior of particles and fluid-particle interactions are closely linked to physical properties of the particles, operating conditions and bed configurations. Fluidized beds behave quite differently when solid properties, gas velocities or vessel geometries are varied. An understanding of hydrodynamic changes and how they, in turn, influence the transfer and reaction characteristics of chemical and thermal operations by variations in gas-solid contact, residence time, solid circulation and mixing and gas distribution is very important for the proper design and scale-up of fluidized bed reactors. In this paper, rather than attempting a comprehensive survey, we concentrate on examining some important positive and negative impacts of particle sizes, bubbles, clusters and column walls on the physical and chemical aspects of chemical reactor performance from the engineering application point of view with the aim of forming an adequate concept for guiding the design of multiphase fluidized bed chemical reactors.One unique phenomenon associated with particle size is that fluidized bed behavior does not always vary monotonically with changing the average particle size. Different behaviors of particles with difference sizes can be well understood by analyzing the relationship between particle size and various forces. For both fine and coarse particles, too narrow a distribution is generally not favorable for smooth fluidization. A too wide size distribution, on the other hand, may lead to particle segregation and high particle elutriation. Good fluidization performance can be established with a proper size distribution in which inter-particle cohesive forces are reduced by the lubricating effect of fine particles on coarse particles for Type A, B and D particles or by the spacing effect of coarse particles or aggregates for Type C powders.Much emphasis has been paid to the negative impacts of bubbles, such as gas bypassing through bubbles, poor bubble-to-dense phase heat & mass transfer, bubble-induced large pressure fluctuations, process instabilities, catalyst attrition and equipment erosion, and high entrainment of particles induced by erupting bubbles at the bed surface. However, it should be noted that bubble motion and gas circulation through bubbles, together with the motion of particles in bubble wakes and clouds, contribute to good gas and solids mixing. The formation of clusters can be attributed to the movement of trailing particles into the low-pressure wake region of leading particles or clusters. On one hand, the existence of down-flowing clusters induces strong solid back-mixing and non-uniform radial distributions of particle velocities and holdups, which is undesirable for chemical reactions. On the other hand, the formation of clusters creates high solids holdups in the riser by inducing internal solids circulations, which are usually beneficial for increasing concentrations of solid catalysts or solid reactants.Wall effects have widely been blamed for complicating the scale-up and design of fluidized-bed reactors. The decrease in wall friction with increasing the column diameter can significantly change the flow patterns and other important characteristics even under identical operating conditions with the same gas and particles. However, internals, which can be considered as a special wall, have been used to improve the fluidized bed reactor performance.Generally, desirable and undesirable dual characteristics of interaction between particles and fluid are one of the important natures of multiphase flow. It is shown that there exists a critical balance between those positive and negative impacts. Good fluidization quality can always be achieved with a proper choice of right combinations of particle size and size distribution, bubble size and wall design to alleviate the negative impacts.     

A nonlinear model for aerodynamic configuration of wake behind horizontal-axis wind turbine

Determination of the aerodynamic configuration of wake is the key to analysis and evaluation of the rotor aerodynamic characteristics of a horizontal-axis wind turbine.According to the aerodynamic configuration, the real magnitude and direction of the onflow velocity at the rotor blade can be determined, and subsequently, the aerodynamic force on the rotor can be determined. The commonly employed wake aerodynamic models are of the cylindrical form instead of the actual expanding one. This is because the influence of the radial component of the induced velocity on the wake configuration is neglected. Therefore, this model should be called a "linear model". Using this model means that the induced velocities at the rotor blades and aerodynamic loads on them would be inexact. An approximately accurate approach is proposed in this paper to determine the so-called "nonlinear" wake aerodynamic configuration by means of the potential theory,where the influence of all three coordinate components of the induced velocity on wake aerodynamic configuration is taken into account to obtain a kind of expanding wake that approximately looks like an actual one. First, the rotor aerodynamic model composed of axial(central), bound, and trailing vortexes is established with the help of the finite aspect wing theory. Then, the Biot-Savart formula for the potential flow theory is used to derive a set of integral equations to evaluate the three components of the induced velocity at any point within the wake. The numerical solution to the integral equations is found,and the loci of all elementary trailing vortex filaments behind the rotor are determined thereafter. Finally, to formulate an actual wind turbine rotor, using the nonlinear wake model, the induced velocity everywhere in the wake, especially that at the rotor blade,is obtained in the case of various tip speed ratios and compared with the wake boundary in a neutral atmospheric boundary layer. Hereby, some useful and referential conclusions are offered for the aerodynamic computation and design of the rotor of the horizontal-axis wind turbine.     

Numerical investigation of full scale coal combustion model of tangentially fired boiler with the effect of mill ducting

In this paper a full scale combustion model incorporating upstream mill ducting of a large tangentially fired boiler with flue gas recirculation was examined numerically. Lagrangian particle tracking was used to determine the coal particle paths and the Eddy Dissipation Model for the analysis of the gas phase combustion. Moreover volatiles and gaseous char products, given off by the coal particles were modelled by Arrhenius single phase reactions and a transport equation was solved for each material given off by the particles. Thermal, prompt, fuel and reburn NO _x models with presumed probability density functions were used to model NO _x production and the discrete transfer radiation model was used to model radiation heat transfer. Generally, the findings indicated reasonable agreement with observed qualitative and quantitative data of incident heat flux on the walls. The model developed here could be used for a range of applications in furnace design and optimisation of gas emissions of coal fired boiler plants.     

Reliable plate–plate MRF magnetorheometry based on validated radial magnetic flux density profile simulations

The apparent shear stress from plate–plate magnetorheometry, using the commercial magnetocell MRD180/1T (Physica/Anton Paar) in standard configuration, is distinctly overestimated. The effect is due to a flux density maximum near the sample rim and radial migration of iron particles towards the rim. Radial magnetic flux density profiles were investigated by finite element simulations using the Maxwell^®2D code and by direct Hall probe measurements. The reliability of the finite element method results, both for the empty magnetocell and with magnetorheological fluid (MRF) in the measuring gap, allows conclusions on the true flux density within the MRF, which cannot be accessed by Hall probe measurements. If the MRF sits on top of the bottom yoke (standard configuration), the flux density maximum reaches twice the plateau value (0.74 T for 3 A coil current and 0.3 mm gap height of the investigated MRF). This yields a higher effective flux density and causes radial iron particle migration, resulting in an additional magnetic flux increase near the rim due to augmented MRF magnetisation. As a result, the rotor torque at constant rotary speed increases with time. Reliable results are achieved by a modification of the magnetocell, such that the MRF sits on a non-magnetic Hall disc of 1.5 mm thickness, allowing an online flux density measurement by a FW Bell Hall probe. In this configuration, the radial flux density profile near the rim remains sufficiently flat and no iron particle migration is detected.     

杆状弯曲行波型超声电机运动机理研究

为了深入研究杆状弯曲行波型超声电机的运动机理，本文分析了定子表面质点椭圆运动轨迹的形成，并将椭圆运动分解为产生驱动力的有效椭圆运动与使定、转子产生径向滑动的径向直线运动；然后基于有效椭圆运动变化情况，分析了在两相振幅不平衡和任意相位差情况下对转子运动平稳性的影响，为这类电机的设计和性能控制提供理论依据。     

基于五阶WENO格式的燃气在药床中流动过程二维两相流研究

为研究内弹道初始阶段中心点火管燃气在膛内药床中的流动特性和传播规律,设计了可视化点传火实验平台,并进行了膛内假药床的点传火实验。基于加权本质无震荡(weighted essentially non-oscillatory, WENO)格式,构造了膛内轴对称二维内弹道两相流模型,对膛内燃气在假药床中的流动过程进行数值模拟。计算结果与可视化实验结果符合较好,全局压力平均误差为5.35%。表明数值计算准确地描述了燃气流动特性,完整地呈现了点火管燃气在假药床中的发展过程。在点火初始阶段,膛内压力径向效应明显,气相沿径向传播较快,药床药粒基本不会发生运动;随着燃气逐渐在膛内传播,膛内压力呈现径向一致、轴向梯度分布的特征,在压力梯度作用下,气相轴向速度开始占据主导,径向速度在膛底和中部区域减小为零,而固相速度随气相速度变化而变化;气相在到达弹底前,由于固相颗粒的壅塞,会提前出现速度反向波动现象。     

Study of a new macro-particle model for the low-temperature pyrolysis of dried wood chips

 This paper presents a transient one-dimensional mathematical model which simulates the pyrolysis of a single dried wood particle. The porous wood particle is considered as a two-phase system: the solid phase consisting of wood and char and the gas phase consisting of volatiles and tar. Conservation equations for mass, momentum and energy are developed for each phase. Chemical processes are described through an existing one-stage three-reactions scheme, leading separately to char, tar and volatiles evolution dynamics. The variation of transport and physical properties with temperature and with composition is presented by algebraic equations. The model presented in this paper is more advanced than the current models found in literature, since it contains physical effects not included in past models, such as cross diffusion, differing solid and gas phase temperatures and a transient gas phase momentum equation incorporating the wall friction experienced by a fluid flowing through a porous medium. Furthermore, an adequate reference system for enthalpy, based on temperature dependent reaction heats, is used. The mathematical equations with initial and boundary conditions are solved numerically by means of a commercial CFD code (PHOENICS). The validity of the pyrolysis kinetics scheme is examined through comparison with experimental data. Furthermore, the macro-particle model is validated by (1) examining the limitations and importance of the newly-modelled effects (different solid phase and gas phase temperature, cross diffusion and wall friction) and (2) making a comparison between predicted and experimental results for large particles. Received on 18 December 2000     

Solid‐particle erosion in the tube end of the tube sheet of a shell‐and‐tube heat exchanger

Erosion is one of the major problems in many industrial processes, and in particular, in heat exchangers. The effects of flow velocity and sand particle size on the rate of erosion in a typical shell‐and‐tube heat exchanger were investigated numerically using the Lagrangian particle‐tracking method. Erosion and penetration rates were obtained for sand particles of diameters ranging from 10 to 500 µm and for inlet flow velocities ranging from 0.197 to 2.95 m/s. A flow visualization experiment was conducted with the objective of verifying the accuracy of the continuous phase calculation procedure. Comparison with available experimental data of penetration rates was also conducted. These comparisons resulted in a good agreement. The results show that the location and number of eroded tubes depend mainly on the particle size and velocity magnitude at the header inlet. The rate of erosion depends exponentially on the velocity. The particle size shows negligible effect on the erosion rate at high velocity values and the large‐size particles show less erosion rates compared to the small‐size particles at low values of inlet flow velocities. The results indicated that the erosion and penetration rates are insignificant at the lower end of the velocity range. However, these rates were found to increase continuously with the increase of the inlet flow velocity for all particle sizes. The particle size creating the highest erosion rate was found to depend on the flow velocity range. Copyright © 2005 John Wiley & Sons, Ltd.     

复杂气固两相系统的微观结构

流化床中的气固两相流动是一个高度复杂的非线性混沌系统。本文利用激光粒子动态分析仪（ＰＤＡ）得到的循环流化床中颗粒脉动速度信号,采用ＦＦＴ分析了脉动信号的宽频谱特征,在此基础上应用小波法分析了脉动信号的动态特征,得到了颗粒脉动速度的微观结构,指出颗粒脉动速度的非线性特性是流化床具有混沌特性的根源,且在不同的尺度上颗粒脉动速度表现出各向异性的特征。     

A numerical investigation into the solids phase chromatography using a combined continuous and discrete approach

Solids phase chromatography for particle classification is based on different retention times of particles with different properties when they are elutriated through a confined geometry.This work aims at a fundamental understanding of such a technology by using the combined continuous and discrete method.A packed bed is employed as the model confined geometry.The numerical method is compared first with experimental observations,followed by a parametric analysis of the effects on the flow hydrodynamics and solids behaviour of various parameters including the number of injected particles,the superficial gas velocity,the contact stiffness and the diameter ratio of the packed column to the packed particles.The results show that the modelling captures some important features of the flow of an injected pulse of fine particles in a packed bed. An increase in the number of injected particles or the superficial gas velocity reduces the retention time,whereas the contact stiffness does not show much effect over the range of 5×102 to5×104 N/m.It is also found that the effect on the retention time of the diameter ratio of the packed column to the packed particles seems complex showing a non-monotonous dependence.     

A numerical investigation into the solids phase chromatography using a combined continuous and discrete approach

Solids phase chromatography for particle classification is based on different retention times of particles with different properties when they are elutriated through a confined geometry. This work aims at a fundamental understanding of such a technology by using the combined continuous and discrete method. A packed bed is employed as the model confined geometry. The numerical method is compared first with experimental observations, followed by a parametric analysis of the effects on the flow hydrodynamics and solids behaviour of various parameters including the number of injected particles, the superficial gas velocity, the contact stiffness and the diameter ratio of the packed column to the packed particles. The results show that the modelling captures some important features of the flow of an injected pulse of fine particles in a packed bed. An increase in the number of injected particles or the superficial gas velocity reduces the retention time, whereas the contact stiffness does not show much effect over the range of 5 × 10^2 to 5× 10^4 N/m. It is also found that the effect on the retention time of the diameter ratio of the packed column to the packed particles seems complex showing a non-monotonous dependence.     

蒙特卡洛方法数值研究大气颗粒物动力学效应和辐射传输性质

爆发性增强的雾天,空气污染严重能见度低,这与大气边界层湍流性质、悬浮颗粒的动力学及散射性质密切相关.文中基于颗粒群平衡方程和Mie理论,采取加权蒙特卡洛方法,自行开发了Fortran程序.文中计算所得的颗粒尺度分布函数、颗粒散射性质与实验值、理论解一致,验证了数值模型和方法的正确性.此外,数值研究了雾爆发性增强阶段雾滴谱拓宽、能见度降低的机理,讨论湍流输运和颗粒局部聚集效应下颗粒间的碰并过程,并耦合颗粒散射性质,数值分析雾发展中湍流耗散率对颗粒对径向相对速度、系统透过率的影响;以及颗粒对径向相对速度与系统透过率、颗粒尺度的关系.研究结果表明:随着湍流耗散率的增大,颗粒的径向相对速度呈现先缓慢而后快速增大的变化趋势.1000s时刻,湍流的耗散率为1.0×10-2m2/s3,颗粒径向相对速度(无量纲)为0.0969;对于0.6μm的可见光,雾环境颗粒系统的透过率为0.47.此外,雾发展中雾滴易与气溶胶碰并,系统的散射性质与水组成的雾滴系统不同,天气的能见度明显降低.