Yongzheng Li ,Xiaolai Zhang ,Guangwei Zhai,Haitao Zhang,Tao Li, *,Qiwen Sun,Weiyong Ying
1 Engineering Research Center of Large Scale Reactor Engineering and Technology,Ministry of Education,State Key Laboratory of Chemical Engineering,East China University of Science and Technology,Shanghai 200237,China
2 School of Chemistry and Chemical Engineering,Shandong University,Jinan 250100,China
3 State Key Laboratory of Coal Liquefaction and Coal Chemical Technology,Shanghai 201203,China
Keywords:Annular stripper Flow Hydrodynamics Particle velocity Laser Doppler velocimeter
ABSTRACT Particle descent velocities in an annular stripper were measured by a laser Doppler velocimetry(LDV)system.In the radial direction,particle descent velocity was relatively constant in the mid-region of the stripper and increased towards the walls on both sides,exhibiting an anti-U-shaped distribution.Particle descent velocity in the radial mid-region increased with the increase of superficial gas velocity,and the maximum in the outer wall region increased significantly with the increase of solid mass flux.Superficial stripping gas velocity had stronger effect on particle velocity distributions near the stripper gas distributor,and such effect weakened with the increase of the distance from the distributor.Local particle velocity and its radial profiles could be adjusted by changing the superficial stripping gas velocity.Empirical formulas were established to describe the relationships between the local particle velocity and cross-sectional averaged velocity based on the effects of operating conditions and measuring positions.The result showed that the predicted data was in good agreement with the experimental value.
In a free vertical movement,particles and gas usually move either in an upward or downward direction when there is no distribution plate in the reactor.Therefore,the gas-solid two-phase flow in the pipelines has three different conditions[1]:particles move upward with gas(such as a riser),gas-solid phases flow co-currently downward(such as a downer),and gas travels upward while solid travels downward(such as a stripper).In the past decades,many hydrodynamics researches[2-6]about riser and downer have been carried out due to their wide use in commercial applications,but the investigation on complex gas-solid counter-current movement in the stripper was relatively deficient.
When particles transfer from a fluidized bed into another one through standpipe,the interstitial hydrocarbon gases move into another fluidized bed with the particles.If particles are porous solids,the entrained gas in holes should be replaced by the stripping gas[1].For example,in FCC units,catalysts require continuous regeneration due to rapid deactivation.Hydrocarbons,which present in the interstitial gas,and adsorb on the catalyst surface[7],should be stripped before catalysts enter regenerator.During the downward movement of solid catalysts along the tripper,the countercurrent steam strips the interstitial gases from the surface or holes of the particles[8].High efficiency stripper not only enhances the yield of light oil,but also reduces the stream consumption and the load of the regenerator,which is important for saving energy and maintaining heat balance of the units.However,for a commercial stripper,flooding and bridging can happen as undesired phenomena.These operating problems impede solid circulation,reduce gas-solid contact efficiency,and may even stop the particle circulation and then interrupt the process[8].Therefore,proper understanding of the hydrodynamics in strippers is essential in order to choose the appropriate operating parameters and optimize the design process.
The hydrodynamics and stripping characteristics have been found to depend strongly on superficial gas velocity and solid mass flux[9-14].Concentration and velocity profiles of gas-particle two-phases were unevenly distributed in the stripper[15].Moreover,the shape and position of internal components greatly influenced the solids holdup profiles of the stripper[16]and appropriate baffles intensified the gas-solid contact and efficiently improved the stripping efficiency[11,17,18].Various configurations for strippers have been studied[19-21].Zhang et al.[22]found that a two segment annular stripper had higher stripping efficiency than an empty cylinder stripper and disk-donut stripper under the same operation condition.Cui et al.[23]studied the effects of stripping steam injection on stripping efficiency and explored the steam injection configurations.Li et al.[24]proposed a new stripping process,which significantly reduced the energy consumption and simplified the overhead vapor route.
However,among these various investigations about different inlet structures and internals of stripper,few studies were conducted on particle velocity distribution in strippers.Particle velocity and its distribution of the gas-solid flow are important hydrodynamic characteristics,which affect the heat and mass transfer,contacting efficiency,and catalyst residence time distribution.It is very important for appropriate industrial design,accurate operation,successful scale-up and process modification to fully understand the particle flow structure in stripper.
In this study,laser Doppler velocimeter(LDV)measurement technique was used to research the particle velocity characteristic in a large scale annular stripper.Without probe inserted into the fluid,this technique does not interfere with the flow pattern[25].The effects of superficial gas velocity,solid mass flux,and superficial stripping gas velocity on particle velocity were discussed.The purpose of this paper is to deepen the understanding of particle velocity distribution in the stripper and provide useful data for the design and operation of the stripper.

Fig.1.Schematic diagram of large-scale cold model experimental apparatus.1.Riser;2.Stripper;3.Precipitator;4.Feed tank;5.Catalyst regenerator;6.Branched pipe distributor;7.Perforated plate;8.Stripper gas distributor;9-14.Ball valve.
Experiments were conducted in a Plexiglas circulating fluidized bed(CFB)(see Fig.1),which consists of riser,stripper,catalyst regenerator,precipitator and feed tank.The stripper was an annular device coaxial to the riser.The riser was 4.8 m long and 0.15 m i.d.,and the stripper was 1.4 m long and 0.43 m i.d.A branched pipe distributor(see Fig.2(a))with 92 nozzles was equipped at an altitude of 0.2 m in the riser.To make the gas distribution more uniform,a perforated plate(see Fig.2(b))with 208 holes of 0.002 m in diameter was installed above the gas distributor.
Glass beads with a mean diameter of 89 μm,2410 kg·m-3particle density and 1521 kg·m-3bulk density were used as fluidized solids.After being measured by a glass rotameter,the treated compressed air(water,oil,and impurities were removed)was fed into the branched pipe gas distributor.Initially,the solid particles were piled up at the bottom of the stripper.When the valves 9 and 10 were opened,solids traveled into the riser reactor above the perforated plate.Due to gassolid drag,the particles moved upward along the riser with the fluidized gas.At the top of the riser,most of the particles bounced back into the stripper after hitting the hat shaped structure and passed through the stripper from top to bottom.A ring type distributor(see Fig.2(c))with 34 nozzles was installed at the bottom of the stripper(H=2.65 m).In the annular stripper,stripping gas traveled upward,which was contrary to the motion of the particles.Solid particles passed through the stripping section and returned to the CFB riser via a circulating pipe of 0.08 m in diameter.Solid mass flux was adjusted with ball valves installed on the circulating pipe.During the experiment,only a very small amount of particles following with the gas left the reactor and entered a bag filter whose aim was for gas-solid separation.The separated particles were returned periodically to the fluidize bed system.Since this paper only aimed to study the particle velocity characteristic in the stripper,valves 12 and 14 were closed because the particles did not need to enter the catalyst regeneration system.
LDV offers an accurate and applicable method to research particle velocities in a sufficiently dilute fluidization system[26-30].In this study,particle velocity in an annular stripper was investigated by using a TSI LDV system.An illustration of the LDV system is provided in Fig.3,which mainly included laser,multicolor beam separator,transmitting/receiving optics,photo detector module,signal processor,and data analysis system.A 5 W argon-ion laser was employed in LDV to produce a coherent beam that illuminated a particle moving in the flow field.And this particle produced a Doppler shift which was proportional to its velocity.Then,particle velocity was calculated using Doppler frequency,laser wavelength,and the angle between beams.The detailed measurement principle has been mentioned in previous published articles[27,31,32].Unfortunately,this technique can only be used in the low solids fraction.In order to determine the available range of LDV,local solids holdups were obtained with the fiber-optic concentration probe at each axial height.As shown in Fig.4,the solids holdups were higher in both wall(inner wall and outer wall)regions of the annular stripper,and the maximum of local solids holdup was about 0.016,which was higher than the solids holdup obtained under experimental conditions.Thus,the solids holdup in the stripper was low enough to guarantee the accuracy of the LDV system.In addition,the accuracy of particle velocity could be determined by data rate,and higher data rate collected can obtain more stable and accurate measurement data.In our study,the data rate was large enough and the sampling time was typically over 60 s or the number of sampled particle reached 20000 at each measurement position,which was enough to ensure the accuracy of the measurement data.

Fig.2.(a)Branched pipe distributor;(b)perforated plate;(c)stripping gas distributor.
In the experiments,various operating conditions were considered.The superficial gas velocity(based on cross-sectional area of riser),Ug,mean solid mass flux(based on cross-sectional area of riser),Gs,and stripping gas velocity(based on cross-sectional area of stripper),Ug′,were 1.57-2.20 m·s-1,4.42-12.61 kg·m-2·s-1,and 0.000-0.088 m·s-1,respectively.H represents the distance from the measured levels to the perforated plate.The local particle velocity in the stripper was measured at 12 radial positions(r/R=0.42,0.47,0.51,0.56,0.60,0.65,0.70,0.74,0.79,0.84,0.88,0.93)on 6 axial levels(H=2.72,2.82,2.92,3.02,3.12,3.22 m).The top view of the annular stripper and the radial measuring positions(filled circles)are showed in Fig.5(a).In order to eliminate the random error and ensure the precision and reliability of the measured data,the measurement of each point was repeated several times.In addition,the local particle velocities were measured on six different radial measuring paths around the circumference of the annular stripper(as shown in Fig.5(a))to investigate the influence of different measuring directions on the measurement results.The angle between the two adjacent paths was 60°.Fig.5(b)shows the local particle velocities on six different measuring paths when Ug=2.20 m·s-1,Gs=4.42 kg·m-2·s-1,Ug′=0.000 m·s-1,and H=2.92 m.After the gas-solid two-phase in the annular stripper reached steady fluidization state,there was no significant difference in the measured values of local particle velocity along each path.Wang et al.[33]have also uncovered a similar phenomenon:the standard deviations among measurements(at the same radial position)in different directions were less than 0.15%.Therefore,radial symmetry was assumed in this paper and all particle velocities were obtained along one of the measuring paths(green filled circles in Fig.5(a)).In addition,some relevant solids holdups were measured at different operating conditions with the fiber-optic concentration probe to better explain the particle velocity distribution.

Fig.3.Schematic diagram of laser Doppler velocimeter.

Fig.4.Radial profile of solids holdup in the annular stripper.
Fig.6 provides the radial distributions of particle velocity at 6 axial levels.The superficial gas velocity,solid mass flux,and stripping gas velocity were 1.57 m·s-1,4.42 kg·m-2·s-1,and 0.044 m·s-1,respectively.Overall,the particle velocity profiles were non-uniform along radial direction.All the particle velocities were less than zero,which meant the direction of movement of the particles was the same as the direction of the gravity.It should be noted that smaller values denoted larger particle descent velocities.
As shown in Fig.6,radial profiles of particle velocity on each level(except for H=2.72 m)exhibited a relatively stable mid-region(r/R=0.56-r/R=0.84),and then gradually decreased towards the inner wall and outer wall,which was contrary to the radial distribution trend of solids holdup.According to this study(e.g.Fig.4)and previous research[6,34-37],particles tended to accumulate and form a dense annulus near the wall.By comparison,the radial position of the maximum particle descent velocity(r/R=0.42 or r/R=0.93)was almost consistent with that of the maximum solids holdup.

Fig.5.The top view of the annular stripper(including radial measuring positions and measuring paths)and the effect of measuring path on local particle velocity.

Fig.6.Radial profiles of particle velocity in the large scale stripper(Ug=1.57 m·s-1,Gs=4.42 kg·m-2·s-1,Ug′=0.044 m·s-1).
In fact,particle descent velocity was greater than gas descent velocity within the measurement range(H=2.72 m-H=3.22 m)of the current experiments.In the process of particle falling,gravity was the dominant force,while the gas-solid drag force became the resistance.Therefore,the formation reason of the radial particle velocity distribution could be explained as follows:in the wall region on both sides of the stripper,the solids holdup was higher,forming two relatively denser annular layers,which may be related to the wall effect and/or the special hat shaped structure on the top of the riser.A large number of particles rebounded back into the stripper with certain angles and moved downward along the inner and outer walls.The outer particles of the dense annulus had shielding effect on the gas-solid drag force,causing a decrease of upward force applied to the inner particles[28].Therefore,the total upward drag force on particles in this region was reduced,causing a higher slip velocity and a larger descent velocity[6].Meanwhile,the number of particles in the radial mid-region was relatively small due to the particle aggregation in the wall region.According to the energy-minimization multi-scale(EMMS)model[38,39],stripping gas was more inclined to move upward in the mid-region of the stripper with minimal resistance rather than through the dense phase in the wall region.In addition,the shielding effect of outer particles in the mid-region became weaker,so the particles encountered greater upward resistance in the process of falling.It was evident that clustered particles moved downward faster than the dispersed particles.Therefore,the particle descent velocities exhibited an anti-U-shaped distribution along the radial direction of the stripper.
The axial development of particle velocity radial profiles was also displayed in Fig.6.There was no obvious development trend in the radial profile shape along the gas-solid flow direction from H=3.22 to H=2.82 m,which was probably because the axial measurement range was too narrow to observe the whole axial development process of particle velocities.Another reason was that particle velocity might have been fully developed.The variation in particle velocity along the direction of motion could be neglected,which meant acceleration force(gravity)and deceleration force(air drag force and wall friction)on particles were in equilibrium.However,at the bottom of the stripper(H=2.72 m),the radial profile became more uniform,which was attributable to the fact that particles in this region were closer to the stripping distributor and were strongly influenced by the stripping gas.The effects of stripping gas velocity on the radial profiles of particle velocity will be discussed in Section 3.4.
Fig.7(a)showed that the superficial gas velocity clearly affected the local particle velocity.Under the same solid mass flux and stripping gas velocity,increasing superficial gas velocity led to larger particle descent velocity at almost every radial position.It was probably because the particles,which were bounced into the stripper,got more energy both in the wall and mid-region to impel particles downward at higher superficial gas velocity.
In Fig.7(a),with the increment of superficial gas velocity,the particle descent velocity increased obviously in the radial mid-region of the stripper,while small changes occurred in the wall region.To better describe this change,the Probability Density Functions(PDFs)distribution of particle velocity was examined,as shown in Fig.7(b).In the wall regions on both sides of the stripper,with the increase of superficial gas velocity,the shape and peak value(the occurrence probability of the particle velocity was the largest)of the probability distribution curve of particle instantaneous velocity did not change much.In the radial mid-region,the PDFs curve and its peak value were shifted to the left,indicating that the particle descent velocity gradually increased as the superficial gas velocity increased.
In addition,particles in the dispersed state had an apparently random velocity in a direction perpendicular to the mean velocity gradient[40].Therefore,the particles had a horizontal velocity during the falling process in the stripper.When the horizontal motion of particles collided with the wall,the particles rotated and moved away from the wall towards the mid-region.But in the meanwhile,particles in the radial mid-region were subjected to a greater drag force because of higher gas velocity,which caused the particles to move from the mid-region to the walls[36,41].With the increase of superficial gas velocity,particles obtained more energy to move away from the wall while the air drag force was increased and made more particles move towards the walls[36].If these two opposing trends were in equilibrium,the solids holdup would remain relatively stable in the dense wall region.This phenomenon was confirmed by the variation of the radial profiles of solids holdup in Fig.7(c)as a function of superficial gas velocity.As the superficial gas velocity increased,the solids holdup in wall regions on both sides of the stripper decreased slightly,so the particle descent velocity change in these regions was not obvious.In the radial mid-region of the stripper,as the superficial gas velocity increased,the solids holdup decreased significantly,but the effect of the superficial gas velocity(downward in stripper)on the downward particle was greater than that of the upward drag force.Finally,the particle descent velocity increased with superficial gas velocity.
In addition to superficial gas velocity,solid mass flux is another important factor influencing the particle velocity profiles of the gas-solid system.Radial profiles of particle velocity under different solid mass flux were plotted in Fig.8(a).Particle descent velocity in the outer wall region obviously increased with solid mass flux,but altered slightly in other regions of the stripper.This implied that the difference between the maximum and minimum velocities was greater under higher solid mass flux than that under the lower solid mass flux.As shown in Fig.8(b),as the solid mass flux increased,the peak of the particle velocity probability density curve shifted significantly towards the left side of coordinate axis(representing the particle descent velocity increased)at the r/R=0.93 position where the maximum particle descent velocity appeared.However,in the inner wall region and the radial mid-region,the peak value of the PDFs changed slightly to the right(representing the particle descent velocity decrease).This phenomenon could be explained by the radial profiles of solids holdup shown in Fig.8(c):with the increase of solid mass flux,the crosssectional averaged solids holdup increased,which made more particles move horizontally from the radial mid-region towards the wall under the action of gas drag force.At the same time,the wall effect was aggravated,and more particles in the outer wall region were moved towards the mid-region by the wall shear force[36].However,particles from the wall region had low kinetic energy and had difficulty in entering into the mid-region where the energy was relatively high[42].Therefore,under the combined effect of the two opposite trends,a large number of particles gathered in the region near the outer wall(r/R=0.93)and eventually formed a dense ring.Besides,compared with dispersion in dilute phase,particles tended to aggregate in dense clusters[39].As the solid mass flux increased,the solids holdup in the wall region became higher and the shielding effect of the outer particles enhanced.Therefore,particle descent velocity increased obviously in the dense annular region.In addition,particle descent velocity reduced weakly in the mid-region as the stripping gas tended to move upward from the radial mid-region.
As shown in Fig.6,the effect of stripping gas velocity on particle velocity radial profiles gradually decreased along the stripping gas flow direction.It was mainly because the effect on the particles in upper space was hindered by the gas-solid suspension segment.In the process of passing through the stripper from the bottom to the upper section,stripping gas was bound to lose energy.In the near wall region,suspended particle concentration was higher.Therefore,the energy loss was more serious than that in the mid-region because of the shielding effect of the outer particles during the stripping gas ascent process.Taking the radial profiles at Ug′=0.044 m·s-1as an example,the particles in the upper levels were less affected by the stripping gas velocity and the overall radial profiles remained obviously non-uniform.
Fig.9,which compared the radial profiles of particle velocity among different stripping gas velocities in the upper section of stripper,indicted that radial profiles were similar while the non-uniformity increased with stripping gas velocity.Particle descent velocity,especially in the radial mid-region,decreased obviously with increasing stripping gas velocity,which was mainly because the more stripping gas tended to move upward in the mid-region with minimal resistance.In addition,higher solids holdup near the wall made the particle velocity insensitive to the change of stripping gas velocity.
In Fig.10,the effects of the stripping gas velocity on the radial profiles of particle velocity at H=2.72 m were provided.With the increase of the stripping gas velocity,radial profiles became uniform and the maximum particle descent velocity disappeared under a fixed solid mass flux and superficial gas velocity.As the stripping gas velocity further increased,the U-shaped distribution arose possibly because of the special geometry of the stripping gas distributor(see Fig.2(c))at the bottom section.Moreover,under higher solid mass flux,larger striping velocity was needed to make the maximum of particle descent velocity disappear.
Fig.11 further showed the effects of superficial stripping gas velocity on particle velocity in the lower section of the stripper(H=2.82 m).Particle descent velocity always decreased with the increase of the stripping gas velocity,but this trend was particularly apparent in the wall regions on the both sides of the stripper,which was mainly because particle aggregation in these regions became difficult so that the dense annulus could not be formed due to the increase of superficial stripping gas velocity(as shown in Fig.12).Eventually,the shielding effect of the outer particles on the drag force mentioned above was reduced or even disappeared.
Fig.13 showed the effect of superficial stripping gas velocity on the instantaneous particle velocity distribution and its occurrence frequency at different radial positions.As the superficial stripping gas velocity increased,the probability density curves of the particle velocity moved to the right,which meant that the descent velocity of the particles at each local position was reduced.The specific performance was that the variation of particle velocity distribution curve in the radial middle area of the stripping section was weak,while that in the inner and outer wall regions changed significantly.Taking the PDFs of particle velocity at r/R=0.88 and r/R=0.42 as examples,not only the peak of PDFs curve shifted to the right obviously,but also the velocity distribution was more concentrated.
Fig.10 also showed that when solid mass flux and stripping gas velocity were constant,particle descent velocities always increased with the increase of Ug.This is similar to the situation described previously in Section 3.2.Besides,the radial profiles of particle velocities became more uneven as the solid mass flux increased.Taking the radial profiles at Ug′=0.044 m·s-1as an example,when solid mass flux went up to 12.61 kg·m-2·s-1,the maximum descent velocity occurred at r/R=0.93.When the stripping gas velocity was less than or equal to 0.044 m·s-1,the variation of particle velocity radial profiles with solid mass flux was similar to what has been discussed in Section 3.3.However,the particle descent velocities increased in all radial regions with increasing solid mass flux when stripping gas velocity was higher than 0.044 m·s-1.This could be ascribed to the cross-sectional averaged solids holdup increased with the enhancement of solid mass flux.The particles,which should have been clustered in the wall region,were more or less dispersed by upflow with the increase of stripping gas velocity because of the special structure of the stripper gas distributor and the increase of superficial stripping gas velocity.When stripping gas velocity went up to some extent,particles were difficult to aggregate in the wall region and tended to move to the radial mid-region and ultimately resulted in higher solids holdup(as shown in Fig.12),so that particle descent velocity in this region increased with the increase of solid mass flux.

Fig.7.Effect of superficial gas velocity on(a)radial profiles of particle velocity,and(b)probability density function of particle velocity,and(c)radial profiles of solids holdup.

Fig.8.Effect of solid mass flux on(a)radial profiles of particle velocity,and(b)probability density function of particle velocity,and(c)radial profiles of solids holdup.

Fig.9.Radial particle velocity profiles under different stripping gas velocities.(Ug=1.57 m·s-1,Gs=4.42 kg·m-2·s-1).
Therefore,the stripping gas velocity has a noticeable impact on the macro gas-solid flow.It is reasonable to expect that superficial stripping gas velocity can be used to adjust the residence time distribution of catalysts in the stripper of the industrial reactor according to the actual need.

Fig.10.Radial particle velocity profiles at the axial level of H=2.72 m.

Fig.11.Effects of stripping gas velocity on particle velocity at different radial positions.
The other purpose of this paper was to propose a suitable empirical formula for predicting the particle velocity in the annular stripper.From the experimental results in Figs.6-11,we may conclude that the superficial gas velocity,solid mass flux,and stripping gas velocity jointly affect the radial profiles of particle velocity in the annular stripper.Therefore,particle velocity can be considered as a function of the particle/gas properties,stripper structure,measuring positions,and operation conditions,i.e.

To make the Vpterm dimensionless,it was modified into the form asSo Eq.(1)could be written as


Fig.12.Effects of stripping gas velocity on radial profiles of solids holdup.(Ug=1.89 m·s-1,Gs=16.8 kg·m-2·s-1).

Fig.13.Effect of superficial stripping gas velocity on probability density function of particle velocity at different radial positions(Ug=1.57 m·s-1,Gs=4.42 kg·m-2·s-1,H=2.72 m).

To test the predictive abilities of these correlations,Eqs.(3)-(5)were verified against independent data obtained under different opera-tion conditions given in Table 2.By comparing the calculated values with the experimental data used to establish the empirical formula and the independent data,Fig.14 showed that Eqs.(3)-(5)had an average relative error of 12.47%which implied great performance in predicting the dimensionless quantity

Table 1 Operating conditions for fitting correlations

Table 2 Operating conditions used in verifying experiment
The radial profiles of particle velocity in a large scale annular stripper were studied with the help of a laser Doppler velocimetry(LDV)system.The particle movement directions at almost all of the radial positions were the same as the direction of gravity within the experimental operating range.Particle velocity radial profiles exhibited a relatively stable mid-region(r/R=0.56-0.84),while in the wall region on both sides of the stripper(r/R<0.56 or r/R>0.84)the particle descent velocity gradually increased towards the stripper walls.
The effect of different operating conditions on particle velocity distribution was obvious.When superficial gas velocity enlarged from 1.57 to 2.20 m·s-1(under a constant solid mass flux of 4.42 kg·m-2·s-1),particle descent velocity increased at every radial position,especially in the mid-region of the stripper.Moreover,with solid mass flux ranging from 4.42 to 12.61 kg·m-2·s-1,the maximum of particle descent velocity in the outer wall region increased obviously and the radial profiles became more non-uniform.Furthermore,stripping gas velocity had a tremendous influence on particle descent velocity and particle velocity radial profiles.The effect of stripping gas velocity on the radial profiles near the stripper gas distributor was quite obvious and gradually diminished away from the distributor.Because of the special geometry of the stripping distributor and the increase of the stripping gas velocity,uniform or U-shaped radial distribution was observed in the lower levels of the stripper.
Based on effects of the superficial gas velocity,solid mass flux,and stripping gas velocity,empirical formulas were developed for predicting thein the annular stripper.Upon examination,the calculated results of these formulas agreed well with the experimental data.
Chinese Journal of Chemical Engineering2019年10期