Solids mixing in a fluidized bed riser

The solids mixing in a riser with a height of 10 m and 0.186 m inner diameter was investigated by using pneumatic phosphor tracer technique. Considering the shielding effect of the bed material on the light emitted from the phosphor tracer particle, a modified method for the phosphor tracer measurement is proposed. And then the curves of particle residence time distribution were obtained. The experimental results show that the particle diffusion mechanism can be explained by the dispersions of dispersed particles and particle clusters in the axial direction, and as well the core-annulus nonuniform distribution of the solids fraction in the radial direction of the riser. Moreover, based on the experimental results, a two-dimensional dispersion model was established to predict the solids axial and radial diffusion. Furthermore, the effects of superficial gas velocity and solids circulating flux on the axial and radial Peclet number of the particles were discussed; two empirical correlation formulas about the axial and the radial Peclet numbers were given; the calculated values agree well with the experimental results.     

内构件对于提升管中颗粒混合行为的影响

对加设了钝体式内构件的提升管内颗粒的轴径向混合为进行了研究。研究结果表明,内构件的存在并不能明显减少提升管内颗粒的轴向返混,这与提升管内稳定的微观两相结构密切相关,但颗粒的径向混合能力却可大大加强,而且在越高的气速和固含条件下,这种加强作用越明显,其原因是内构件的存在加强了边壁区颗粒的脉动,有利于破坏提升管边壁区的颗粒浓环,促进颗粒的径向交换与混合。     

Mixing of a lateral gas stream in a two-dimensional riser of a circulating fluidized bed

A cold model of a circulating fluidized bed having a two-dimensional riser, with a 12 × 120 mm section and a 6.4 m height, was equipped with a device to inject a lateral gas stream along the riser. The apparatus was operated under conditions ranging from those characteristic of combustors to those of gas-conversion processes. Flow structures in the interaction region between the rising gas—solids suspension and the lateral gas stream were studied by means of a motion analysis system. Three main configurations were identified. The ratio between the momentum of the lateral gas stream and that of the rising suspension was found to be the parameter able to discriminate among the three configurations. A satisfactory agreement was found with mixing data obtained using a cylindrical riser.     

Segregation and mixing effects in the riser of a circulating fluidized bed

Segregation and mixing effects of binary mixtures of particles having difference in sizes and densities were studied in 0.1016 m-diameter riser of a circulating fluidized bed at gas velocities between 2.01 and 4.681 m/s and solids circulation rate between 12.5 and 50 kg/m² s. Two groups of bed materials (three quartz sand-spent fcc catalyst mixtures with different initial mass % of sand and two coal-iron mixtures, one with almost same sizes but with different densities and the other having both different sizes and densities) were used. Using local axial mass % of heavier/coarser particles and their mean sizes the extent of segregation was evaluated. The influence of operating conditions like superficial gas velocity and solids circulation rate on segregation was examined and found that with their increase segregation effects generally tend to decrease and a uniform mixture conforming to initial composition of the mixture results. Using the data available in the literature and those of the present authors an empirical correlation to obtain the extent of segregation in CFBs has been proposed.     

Study of cracking FCC naphtha in a secondary riser of the FCC unit for maximum propylene production

To satisfy the increasing propylene demand, reprocessing FCC naphtha in a secondary riser of the FCC unit was investigated. To this aim, a full range FCC naphtha was cracked over a mixture of two kinds of commercial equilibrium FCC catalysts, which contained 95 t.% Y zeolite-based catalyst and 5 wt.% ZSM-5 zeolite-based additive. The effects of operating parameters such as reaction temperature (temperature of the riser outlet), catalyst-to-oil ratio and residence time on FCC naphtha cracking were studied in a continuous pilot plant. This work demonstrates that FCC naphtha requires high operating severities to crack, and approximately 12–19 wt.% FCC naphtha can be transformed into propylene. The conversion and yield of propylene showed a rapid increase with increasing reaction temperature, and the increase of catalyst-to-oil ratio also enhanced FCC naphtha cracking, even at high reaction temperature. However, at high catalyst-to-oil reactions, hydrogen-transfer reactions constrain further increases in light olefin yields. At these high operating severities, shortening residence time is an appropriate way to obtain high yields of propylene combined with (i) lower yields of dry gas and (ii) a lower apparent hydrogen-transfer coefficient.     

Numerical simulation of flow behavior of particles and clusters in riser using two granular temperatures

Cluster in CFB riser significantly affects performance of circulating fluidized beds. To model hydrodynamic behavior in CFB risers, three phase flows were assumed in the riser, the gas phase, the dispersed particle phase, and the clusters phase. The gas-solid multi-fluid model is extended to give the macroscopic averaged equations with constitutive equations for both particle phases from kinetic theory of granular flow. The clusters and the dispersed particles have their own fluctuating energy or two individual granular temperatures. Interactions between the cluster and its surrounding dispersed particles were obtained from kinetic theory of granular flow. Drag force for gas to dispersed particles and the clusters are empirically determined. The momentum exchange between dispersed particles and clusters is modeled using the concept of molecular dynamics. Cluster properties are predicted with the cluster-based approach. The distributions of volume fractions and velocities of gas, dispersed particles and clusters are predicted. Computed solid mass fluxes and volume fractions agree with Manyele et al. [S.V. Manyele, J.H. Parssinen, J.X. Zhu, Characterizing particle aggregates in a high-density and high-flux CFB riser, Chemical Engineering Journal, 88 (2002) 151-161.] and Knowlton [T.M. Knowlton, Modelling benchmark exercise. Workshop at the Eighth Engineering Foundation Conference on Fluidization, Tours, France, 1995.] experimental data.     

DSMC prediction of granular temperatures of clusters and dispersed particles in a riser

Flow behavior of gas and particles in a CFB riser is simulated by a Euler-Lagrangian approach. Collisional interaction of particles is modeled by a direct simulation Monte Carlo (DSMC) method. The coefficient of restitution depends on the relative velocity between two particles by taking into account both the viscoelastic and plastic deformations of particles. The Newtonian equations of motion are solved for each simulated particle in the system. The interaction between gas phase and simulated particles is determined by means of Newtonian third law. A criterion proposed by Soong et al. [C.H. Soong, K. Tuzla, J.C. Chen, Experimental determination of clusters size and velocity in circulating fluidized beds, in: J.F. Large, C. Laguerie Eds., Fluidization VIII Engineering Foundation, New York, 1995, pp. 219-227.] is used to obtain information of clusters. A model for the determination of the number of particles inside a cluster is proposed from kinetic theory of granular flow. The flow behavior of clusters and dispersed particles in the riser is numerically predicted. Distributions of granular temperature of clusters and dispersed particles are computed from velocity distributions. The upward moving clusters give a high granular temperature in comparison to the downward moving clusters. The granular temperature of the dispersed particles is an order of magnitude larger than that of the clusters. The computed collision frequencies of dispersed particles and clusters from DSMC are lower than those from kinetic theory of granular flow.     

Solids mixing in a down-flow circulating fluidized bed of 0.418-m in diameter

The axial and lateral solids mixing in a down-flow circulating fluidized bed of 0.418-m diameter was investigated by a pneumatic injection phosphor tracer technique (PIPTT). The axial and lateral solids dispersion were determined by measuring the solids RTD at same axial but different lateral positions using point sources for tracer injection. A two-dimensional dispersion model described the measured RTD curves satisfactorily. The results were compared to those obtained in the small scale downers and the scale-up effect was investigated. The axial solids Peclet number Pe_a is around 110 and invariable with changing U_g, G_s and ?_s, while the lateral solids Peclet number Pe_r is linearly increasing with ?_s. And Pe_r is found to decrease with the square root of inner diameter (ID) in comparison with the results obtained in small ID downers. Correlation of Pe_r, Pe_r = (15 + 70.7 ?_s)D^− 0.5, is proposed.     

The residence time distribution and mixing of the gas phase in the riser of a circulating fluidized bed

CFBs are increasingly used for both gas-catalytic and gas-solid reactions. The conversion is a function of the gas hydrodynamics, subject of the present research.Available literature on the gas mixing in the riser of a CFB shows contradictory results: some investigators neglected back-mixing of gas, whereas others report a considerable amount of back-mixing in CFB risers. The present paper reports experimental findings obtained in a 0.1 m I.D. riser, for a wide range of combined superficial gas velocity (U) and solid circulation flux (G). The gas flow mode (plug vs. mixed) is strongly affected by the operating conditions, however with a dominant mode within a specific (U, G)-range. Sand was used as bed material. The superficial gas velocity was varied from 5.5 to 8.3 m/s, the solids circulation flux was between 40 and 170 kg/m² s. A tracer pulse response technique was used with a pulse of propane injected at the bottom and detected at the riser exit. The cumulative response curves, F(t), define (i) an average residence time (t₅₀) obtained for F(t) = 0.5; and (ii) the slope of the curves (a steeper one corresponding with more pronounced plug flow) and expressed in terms of a span, σ. These parameters (t₅₀ and σ) define the gas flow mode. A quantitative comparison of experimental results with literature RTD-models is inconclusive although the occurrence of both mixed flow and plug flow is evident, and (U, G)-dependent. The experimental results are expressed in empirical design equations, and the comparison of predicted and experimental results is fair: low values of σ determine the plug flow regimes, whereas back-mixing is more pronounced at higher value of σ. Experiments with similar systems might favor plug flow or mixing as function of the combined (U, G)-values. The introduction of the RTD-function in reaction rate equations can improve the prediction of the gas-conversion in a riser-reactor.     

Kinetic theory based CFD simulation of turbulent fluidization of FCC particles in a riser

The turbulent fluidization regime is characterized by the co-existence of a dense, bottom region and a dilute, top bed. A kinetic theory based CFD code with a drag corrected for clusters captured the basic features of this flow regime: the dilute and dense regions, high dispersion coefficients and a strong anisotropy. The computed energy spectrum captures the observed gravity wave and the Kolmogorov -5/3 law at high frequencies. The computed turbulent kinetic energy is close to the measurements for FCC particles. The CFD simulations compared reasonably well with the measured core-annular flow experiments at very high solid fluxes. The computed granular temperatures, solids pressures, FCC viscosities and frequencies of oscillations were close to measurements reported in the literature. The computations suggest that unlike for the flow of group B particles, the oscillations for the FCC particles in the center of the riser are primarily due to the oscillations of clusters and not due to oscillations of individual particles. Hence mixing is not on the level of individual particles.     

循环流化床提升管中固体颗粒停留时间的分布

在内径１４０ｍｍ，高１０ｍ的循环流化床提升管中，采用磷光颗粒示踪法对床内固体颗粒的停留时间分布进行了测定。在气速１．５～９．０ｍ／ｓ，固体循环量１０～１４０ｋｇ／ｍ２ｓ的范围内，实验测得的停留时间分布曲线均有明显的双峰分布。这种双峰分布是由于提升管中弥散颗粒和颗粒团共同作用的结果。本文提出的一维两组分扩散叠加模型可较好地描述提升管中固体的混合行为。考察了在实验条件下，操作条件对固体混合行为的影响。发现：气速及固含对颗粒的轴向Ｐｅｃｌｅｔ数影响不大，提升管中颗粒的返混主要是由于颗粒团引起的。将研究结果与近期文献报道的其他研究进行了对比     

Modeling of vaporization and cracking of liquid oil injected in a gas-solid riser

Vaporization and cracking of liquid oil injected in a gas-solid riser (fluid catalytic cracking riser reactor) was computationally studied in this work. Evaporation of a single drop injected in a stream of gas-solid mixture was analyzed first. A model for simulating evaporation of a drop considering heat transfer from the gas phase as well as from the solid particles was developed. The model relates the evaporation rate of droplet with rate of collisions of solid particles, specific heat capacities of solid and liquid, latent heat of vaporization, relative velocity of gas and liquid and temperatures of three phases. The understanding gained from such a model was then extended to simulate evaporation of liquid drops injected in FCC risers. The Eulerian-Lagrangian approach was used to simulate simultaneous evaporation and cracking reactions occurring in FCC riser reactors. A commercial CFD code, FLUENT (of Fluent Inc., USA) was used. Four and ten lump models were used for simulating cracking reactions. Appropriate user defined functions were developed to implement heterogeneous kinetics and heat transfer models in FLUENT. A special algorithm was developed to calculate accumulated coke on catalyst particles. A boiling point range was considered for simulating realistic oil feedstock. The model was first evaluated by comparing predicted results with published industrial data. The simulations were then carried out to understand influence of key design and operating parameters on performance of FCC riser reactors. The parameters studied included; initial oil droplet distribution, catalyst inlet temperature, catalyst to oil ratio and thermal cracking. The approach, model and results presented here would be useful for optimization of FCC operation, cost to benefit analysis of new FCC nozzles and related decision-making.     

Mixing and segregation in a liquid fluidized bed of particles with different size and density

Experimental studies have been carried out on fluidization of irregular particle mixtures of different size and density. The mixing and segregation phenomena could be interpreted on the basis of the diffusion model of Kennedy and Bretton. The dependence of computed particle dispersion coefficient on liquid velocity, particle density and size has been discussed.     

Species segregation of binary mixtures and a continuous size distribution of Group B particles in riser flow

Experiments involving a gas–solid, pilot-scale circulating fluidized bed (CFB) have been carried out, with a focus on species segregation measurements in a riser. Three mixtures were considered: (i) a binary mixture with particles of different sizes (d_ave) but same material density (ρ_s), (ii) a binary mixture with particles of different material densities (ρ_s) but same size (d_ave), and (iii) a continuous particle size distribution (PSD). Local measurements of the composition (i.e., species segregation) of each mixture were obtained over a range of operating conditions. Similar to previous works, the results show that the more massive species (i.e., greater d_ave or ρ_s) preferentially segregates toward the wall in all cases. Several new trends were also observed. First, for the binary mixtures, composition of the more massive species increases with riser height at the wall under some operating conditions. The operating conditions that cause this phenomenon are mutually exclusive for the size-difference and density-difference systems. Second, for the continuous PSD, radial segregation is observed even when there is a net positive flux in the annular region, contrary to previous findings which indicated segregation only for conditions leading to a net downward flux in the annular region. Finally, two qualitative differences between the binary and continuous mixtures were noted: (i) a monotonic decrease in species segregation is observed for the binary mixtures with an increase in the solid loading (m), while a non-monotonic trend is observed for the continuous PSD, and (ii) while the shape of the radial segregation profile is flattest at the riser bottom for the binary mixtures, the flattest radial profile is at the riser top for the continuous PSD.     

Experimental profiles of lateral mixing of feed particles in a three‐dimensional fluidized bed

Solids mixing data of high quality is one of the most crucial steps for quantitative studies, but it is a difficult task to obtain in a fluidized bed especially with a 3D configuration. Therefore a novel sampling technique is developed with bed collapse method, for measuring lateral mixing of feed particles in a 3D fluidized bed. The sampling tool is designed using a “bottom‐to‐top sampling” idea. Its development, configuration and measurement repetition are discussed in detail. The effects of mixing time, fluidizing gas velocity, and particle size of bed material on the tracer distribution are investigated. A quantitative comparison of lateral dispersion coefficient shows that our results agree fairly well with measurements and predictions of correlations for lab‐scale fluidized systems in previous studies. The presented 2D profiles of the lateral mixing can be used to validate fundamental solids mixing models or verifying convenient measurement techniques. © 2010 American Institute of Chemical Engineers AIChE J, 2011     

逆向气体射流对下行床颗粒混合的影响

下行床入口结构的研究一直被人们所重视。今在内径为0．192m的下行床中颗粒达到均匀分布的部位,沿床四周均布了三个45度方向逆流场气体射流入口,在此处设立气体入口可以使颗粒分散与气固快速接触不再同时进行。采用磷光颗粒示踪技术对下行床有逆流场射流气体存在时颗粒的轴径向混合行为进行了研究。这种逆流场射流气体对下行床颗粒的轴向混合行为无明显影响,在各操作条件下下行床内颗粒均能以接近平推流的方式运动;但该射流气体可以大大加强颗粒的径向混合,有利于气固接触,在下行床颗粒径向混合越差的操作条件下,射流气体对颗粒径向混合的影响效果越明显,下行床的这种入口结构具有良好的应用前景。     

Simulation of particles and gas flow behavior in a riser using a filtered two-fluid model

Flow behavior of gas and particles is predicted by a filtered two-fluid model by taking into the effect of particle clustering on the interphase momentum-transfer account. The filtered gas–solid two-fluid model is proposed on the basis of the kinetic theory of granular flow. The subgrid closures for the solid pressure and drag coefficient (Andrews et al., 2005) and the solid viscosity (Riber et al., 2009) are used in the filtered two-fluid model. The model predicts the heterogeneous particle flow structure, and the distributions of gas and particle velocities and turbulent intensities. Simulated solids concentration and mass fluxes are in agreement with experimental data. Predicted effective solid phase viscosity and pressure increase with the increase of model constant c_g and c_s. At the low concentration of particles, simulations indicate that the anisotropy is obvious in the riser. Simulations show the subgrid closures for viscosity of gas phase and solid phase led to a qualitative change in the simulation results.     

A two-scale second-order moment particle turbulence model and simulation of dense gas-particle flows in a riser

A two-scale second-order moment particle turbulence model is developed, based on the concept of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision. It is then used to simulate dense gas-particle flows in a riser. The predicted results are in reasonable agreement with the experimental results and show the core-annulus flow structure observed in experiments. It is seen that the two-scale model is somewhat better than k_f-ε_f-k_p-ε_p-θ two-phase turbulence model and the single-scale second-order moment two-phase turbulence model (USM-θ model).     

Mixing dynamics in bubbling fluidized beds

Solids mixing affects thermal and concentration gradients in fluidized bed reactors and is, therefore, critical to their performance. Despite substantial effort over the past decades, understanding of solids mixing continues to be lacking because of technical limitations of diagnostics in large pilot and commercial‐scale reactors. This study is focused on investigating mixing dynamics and their dependence on operating conditions using computational fluid dynamics simulations. Toward this end, fine‐grid 3D simulations are conducted for the bubbling fluidization of three distinct Geldart B particles (1.15 mm LLDPE, 0.50 mm glass, and 0.29 mm alumina) at superficial gas velocities U/Umf = 2–4 in a pilot‐scale 50 cm diameter bed. The Two‐Fluid Model (TFM) is employed to describe the solids motion efficiently while bubbles are detected and tracked using MS3DATA. Detailed statistics of the flow‐field in and around bubbles are computed and used to describe bubble‐induced solids micromixing: solids upflow driven in the nose and wake regions while downflow along the bubble walls. Further, within these regions, the hydrodynamics are dependent only on particle and bubble characteristics, and relatively independent of the global operating conditions. Based on this finding, a predictive mechanistic, analytical model is developed which integrates bubble‐induced micromixing contributions over their size and spatial distributions to describe the gross solids circulation within the fluidized bed. Finally, it is shown that solids mixing is affected adversely in the presence of gas bypass, or throughflow, particularly in the fluidization of heavier particles. This is because of inefficient gas solids contacting as 30–50% of the superficial gas flow escapes with 2–3× shorter residence time through the bed. This is one of the first large‐scale studies where both the gas (bubble) and solids motion, and their interaction, are investigated in detail and the developed framework is useful for predicting solids mixing in large‐scale reactors as well as for analyzing mixing dynamics in complex reactive particulate systems. © 2017 American Institute of Chemical Engineers AIChE J, 63: 4316–4328, 2017     

Investigating the dynamics of cylindrical particles in a rotating drum using multiple radioactive particle tracking

The behavior of granular flows inside rotating drums is an ongoing area of research. Only a few studies have investigated non‐spherical particles despite the fact that particle shape is known to have a significant impact on flow behavior. In addition, the experimental techniques limit the interpretation of the results of these studies. In this work, we compared the flow behavior of cylindrical and spherical particles using the multiple radioactive particle tracking technique to capture the positions and orientations of cylindrical particles simultaneously. We analyzed two important components of the transverse flow dynamics, that is, the boundary between the active and passive layers, and the velocity profile on the free surface. For the cylindrical particles, two general models are proposed to calculate the velocity profiles on the free surface and the effective particle sizes in the active and passive layers. © 2016 American Institute of Chemical Engineers AIChE J, 62: 2622–2634, 2016