About this Journal Submit a Manuscript Table of Contents
International Journal of Chemical Engineering
Volume 2012 (2012), Article ID 786982, 13 pages
http://dx.doi.org/10.1155/2012/786982
Research Article

Numerical Studies of the Gas-Solid Hydrodynamics at High Temperature in the Riser of a Bench-Scale Circulating Fluidized Bed

1School of Chemical Engineering, University of Campinas, 500 Albert Einstein Avenue, Campinas 13083-970, SP, Brazil
2Department of Mechanical Engineering, Federal University of Technology of Parana, Monteiro Lobato Avenue, 84016-210 Ponta Grossa, PR, Brazil
3Faculty of Mechanical Engineering, University of Campinas, 200 Mendeleyev Avenue, Campinas, 13083-970, SP, Brazil

Received 8 May 2012; Revised 3 July 2012; Accepted 10 July 2012

Academic Editor: Adrian Schumpe

Copyright © 2012 Maximilian J. Hodapp et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

The hydrodynamics of circulating fluidized beds (CFBs) is a complex phenomenon that can drastically vary depending on operational setup and geometrical configuration. A research of the literature shows that studies for the prediction of key variables in CFB systems operating at high temperature still need to be implemented aiming at applications in energy conversion, such as combustion, gasification, or fast pyrolysis of solid fuels. In this work the computational fluid dynamics (CFD) technique was used for modeling and simulation of the hydrodynamics of a preheating gas-solid flow in a cylindrical bed section. For the CFD simulations, the two-fluid approach was used to represent the gas-solid flow with the k-epsilon turbulence model being applied for the gas phase and the kinetic theory of granular flow (KTGF) for the properties of the dispersed phase. The information obtained from a semiempirical model was used to implement the initial condition of the simulation. The CFD results were in accordance with experimental data obtained from a bench-scale CFB system and from predictions of the semiempirical model. The initial condition applied in this work was shown to be a viable alternative to a more common constant solid mass flux boundary condition.

1. Introduction

Circulating fluidized bed reactors are systems in which gas-solid heterogeneous reactions take place in a fast fluidization regime. One of the most successful applications of this technology is the combustion of low-grade fuels, such as biomass, waste, and coal with high ash content. Fluidized beds provide high specific transfer rates, high solids throughput, and thermal uniformity within the reactor [1].

The total particulate material circulating in the system, the solids inventory, is an important parameter for an efficient design of CFB for combustion applications. Together with the gas superficial velocity it drives the bed particles along the loop. The pressure drop and the particle residence time are strongly related to the solid distribution in the fast bed zone (riser). Moreover, the solid distribution also affects the mass and heat transfer rates. For combustion and gasification purposes, quartz sand is commonly used as inert material due to its relatively low cost and excellent performance at high temperatures. In such systems, the amount of solid inert material can reach up to 97% of the total mass in the solids inventory [2].

Several studies have been carried out to describe the hydrodynamic behavior of the gas-solid flow in CFB systems, based on empirical and theoretical analysis at high or environmental temperatures [37]. The mathematical models developed sought to explain the main characteristics of the CFB system in terms of significant variables, such as geometrical configurations, solids properties and inventory, superficial gas velocity, and solid circulation flux [811]. However, only a few studies reported in the literature describe a straightforward procedure to define the operational variables and geometrical measures of the components in a CFB unit for project purposes [1214].

In a different approach to the mathematical model of CFBs, the literature has shown the successful application of the CFD technique for the simulation of many gas-solid-fluidized bed applications. Great attention has been given to the hydrodynamics in the riser component of the CFB unit. Some works have studied the behavior of FCC catalyst [15], others coal [16], and also sand particles [9].

The Eulerian-Eulerian approach is commonly applied to describe a gas-solid system. Also known as the two-fluid model, it treats the gas as well as the particulate phases as fluids, which introduces some variables that are difficult to determine for the solid phase. Thus, the kinetic theory of granular flow has been widely applied [1719] as an analogy to the kinetic theory of gases, thereby providing closure relations for the transport equation. The two-fluid approach has been applied to different system configurations, also producing reliable results for industrial scale CFBs [20].

Cases where the available experimental data provide only the total solids inventory of the CFB loop are inappropriate for the simulation of individual components, such as risers. For the determination of initial and boundary conditions, many authors [2123] have alternatively applied constant values for the solid mass flux entering the riser. Although this approach is valid, these are mean fluxes, not representing the variation in mass flux history over time that is observed in actual experimental units [24]. Also with a constant solid mass flux boundary condition sufficient simulation time must be computed in order to reach a pseudo-steady state, in which the actual solids inventory and thus the correct multiphase hydrodynamics are achieved in the riser. Depending on the operational and geometrical configuration, this may take too much computational effort to be viable for a tridimensional and transient CFD approach. Knowledge of the actual amount of solid material in each section of the CFB gives a further option for setting the boundary and initial conditions in the simulation of risers. Therefore the aim of this work was to use a semiempirical model to determine the required solids inventory information and verify the validity of an alternative solid reflux boundary condition.

2. Experimental Unit

A schematic view of the bench-scale CFB unit installed at University of Campinas is provided in Figure 1. The loop is composed primarily of a riser and a cyclone for solid separation, followed by a standpipe and L-valve to close the circuit.

786982.fig.001
Figure 1: Experimental assembly setup, adapted from Behainne [25].

Preheated air is supplied to the riser by a screw compressor, while the solids are fed in by a screw feeder located at the bottom of the solids hopper. Solids that are fluidized in the riser and dragged by the gas phase are collected by a tangential cyclone, which is responsible for releasing the solids into the standpipe. The L-valve controls the solid reinjection into the riser, creating a circulating system. Fine particles escaping the unit are collected by a bag filter located before the stack. Temperature measurement points are placed along the riser. Table 1 shows the main characteristics of the system as well as some of the physical and fluidization properties of the bed particles (quartz sand).

tab1
Table 1: Geometrical and physical properties of the system.

Yang [26] revised the particle classification of the Geldart groups for other temperatures and pressures, using the dimensionless density and the Archimedes number to determine the group divisions. Under the CFB operational conditions of the present work the solids, which were initially classified in group B under atmospheric conditions, are still considered to be on the same group.

A set of several experimental tests were conducted on this unit, as reported by Hory et al. [27]. For this study only the experimental data presented in Table 2, which represent the mean values of each variable tested in three runs, were chosen.

tab2
Table 2: Operational conditions of the CFB unit.

3. Numerical Studies

Two approaches were employed in the numerical investigation of the CFB unit, each with a specific aim. Initially a semiempirical model was assembled to determine key operational variables of a CFB system. This model was based on mathematical correlations obtained from the literature. The results produced by the semiempirical model were applied to determine boundary and initial conditions of the second approach, a CFD study of the hydrodynamics of the gas-solid flow in the unit’s riser section. Both model strategies are presented in more detail as follows.

3.1. Semiempirical Model

The practical hydrodynamic model presented by Basu [2] was adapted to determine the bed solids inventory as well as other basic operational and geometrical parameters for a bench-scale CFB system. An overview of the correlations used in the semiempirical model is given in Table 3, which is divided into three stages, one for each component group of the CFB unit (riser, cyclone, and standpipe/solid recycle valve). The model requires input information such as riser diameter and properties of the fluidizing air and the bed particles to be initially defined. Then the solids inventory and other key conditions required for stable operation of the loop are calculated.

tab3
Table 3: Semiempirical CFB mathematical model.

The total solids inventory is calculated as the sum of individual values for the riser and the standpipe/L-valve, since the particles of the CFB system are mainly concentrated in those components.

3.2. CFD Model

For a better understanding of the behavior of the gas-solid flow in the riser, a CFD study was carried out. CFD can produce detailed information on multiphase flows, which is not always readily obtained on experimental units. The gas-solid flow information produced can also be used to aid in the improvement of the system design.

The governing equations for the momentum, turbulence, and energy transport were used to describe the phenomena involved in the multiphase flow inside the riser. The gas-solid system was represented by the Eulerian-Eulerian approach, which considers each phase as a fluid. For describing the turbulence in the gas phase, the realizable k-epsilon model was applied and the solid phase was modeled by the KTGF, following the general approach adopted by many authors for multiphase flow in risers. The main equations in the model are presented in Table 4.

tab4
Table 4: Hydrodynamic and thermal models for the two-fluid gas-solid flow approach.

The solid stress tensor was described through the KTGF, which accounts for the solid viscosity and pressure terms. The granular temperature was related to the kinetic turbulent energy of the particle and the transport equation derived from it in the kinetic theory is simplified to an algebraic formulation. As for the thermal balance in the riser, a transport equation was solved for each phase. The pressure work, kinetic energy, and viscous heating were neglected as the flow occurs at low Mach numbers [16]. The heat exchange between the two phases was expressed by the Ranz and Marshall [40] correlation.

Only the riser was considered for the CFD simulation, since this is the section of interest in most of industrial applications and where experimental data was collected for validation. This simplification also avoids extracomputational efforts for simulating the very dense regions in the L-valve, through which the solids return to the riser. A hexahedral mesh containing approximately 400 thousand control volumes was created for the riser. A set of finer and coarser mesh was also created, whereas the refinement of the one selected was sufficient for producing grid independent results for the nonisothermal multiphase flow. The mean equivalent cell length to particle diameter for the selected grid was around 37. Details of the geometry and numerical grid are shown in Figure 2.

786982.fig.002
Figure 2: Details of (a) the generated geometry, (b) numerical mesh of the top section, and (c) numerical mesh of the bottom section of the riser for the CFD simulation.

Atmospheric air represented the continuum phase and quartz sand the dispersed phase. The boundary conditions are described in Table 5. At the main air entrance (bottom zone), a constant gas mass flow rate was specified, at the temperature measured at the exit of the GLP burner. At the secondary inlet, where the solid returns to the riser, air from the L-valve enters at a constant mass flow rate. The temperature boundary condition of the secondary entrance was determined by experimental measurements at the L-valve exit and was considered the same for both phases as they reenter the riser. A mathematical function was implemented to inject sand at the same mass flow rate as that calculated at the riser exit at each time step. In order to apply this condition it was necessary to determine the mass of sand that would be present only in the riser, since considering the total solids inventory of the entire loop would greatly overestimate the solid recirculation rate, which in turn would modify the system’s fluid dynamic behavior. This information was obtained from the semiempirical model. At the walls the no-slip boundary condition was applied for the gas phase and free slip for the solid phase, a simplified but close to the correct low-friction condition for the dispersed phase, according to [19]. The heat loss to the environment was also considered as a constant negative heat flux.

tab5
Table 5: Boundary conditions for the CFD riser simulation.

In order to reproduce the behavior of the gas-solid flow in a vertical riser, transient simulations were considered. The numerical formulation followed the finite volume method by means of the commercial software Ansys Fluent 12.0. For the interpolation of the discretized partial differential equations, the first-order upwind scheme was used. The convergence criterion was established as 10−3 in absolute values. Reduced underrelaxation factors were chosen to improve the stability of the solution. The time step of 10−5 s was initially adopted and gradually increased to a maximum value of 10−4 s.

4. Results and Discussion

As the main purpose of the semiempirical model was to ease the preliminary system setup for the CFB unit, a validation was carried out comparing its results to those obtained from the experimental data shown in Table 2. Furthermore, the internal solids reflux ratio was set at the appropriate value to guarantee the same riser height of the experimental CFB combustor. Table 6 contains the input data of the model based on the experimental unit conditions.

tab6
Table 6: Input data used in the semiempirical model based on Hory et al. [27].

The results obtained from the semiempirical model for each section of the CFB unit are presented in Table 7, and the comparison with experimental data for superficial gas velocity, solids inventory, and aeration mass flow rate in L-valve is shown in Figure 3.

tab7
Table 7: Results obtained from the semiempirical model.
fig3
Figure 3: Comparison of the semiempirical model results for (a) superficial gas velocity (), (b) solids inventory (), (c) aeration mass flow rate (), and (d) deviation for each variable.

In Figure 3 it is shown that the predictions for , , and were consistent with the experimental values. For the total solids inventory, the result obtained in the simulation was overestimated by approximately 26%. The gas superficial velocity and the aeration rate were very close to the plant operational conditions. In general, these results are considered satisfactory for preliminary design of CFB loops.

According to these results, the amount of solids in the riser corresponds to approximately 43% of the total solids inventory. This information produced by the semiempirical model was used for setting the initial condition for the CFD simulation of the riser. Thus an initial packed bed was set at the base of the riser corresponding to the stipulated solid mass.

The simulation was carried out for 10 s, during which time a quasi-steady state was achieved. At the initial seconds of simulation, the inert bed was observed to begin a rapid fluidization, reaching the riser’s exit and starting to return to the riser. Nevertheless, only after 5 seconds of simulation, transient statistics averages were collected for the main variables in order to analyze the riser characteristics.

The multiphase flow in fluidized beds (risers) is known for its chaotic behavior, as can be observed by the solid phase mass flux leaving the riser (Figure 4(a)). The CFD simulation captured this behavior, although the solid residence time in the cyclone and standpipe was neglected. The calculated at the exit at every time step was applied to the recirculation boundary condition. The mass flux fluctuations at the recirculation entry were computed without any delay. The average solid recirculation mass flux () predicted by the CFD simulation was 59 kg/m2s, which agreed well with the semiempirical estimation, an over-predicted simulated result of 7.5%.

fig4
Figure 4: CFD simulation results of solid phase mass flux leaving the riser: (a) time history fluctuations and (b) power spectral density.

The oscillatory characteristic of the mass flux at the exit can be explained by the continuous disappearance and formation of clusters along the riser, which determines the solid radial distribution and the vertical displacement of particles into the column, also to be observed in Figure 7. The higher mass flux peaks occur as a consequence of this behavior. The slug frequency at the riser exit is around 0.36 Hz (2.8 s time period), as verified by the power spectrum density (PSD) analysis of solids mass flux fluctuation, shown in Figure 4(b). The frequency of the largest peak was within the same range as reported by other research works [21, 41]. This analysis was carried out with a sampling time of 0.01 s.

In Figure 5 comparison of the pressure drop with the values predicted by the CFD and semiempirical models is shown. The total pressure drop for the full height of the riser had a higher value with the semiempirical model; however, considering the region above of the solids return level, the results were similar. This behavior can be better understood by analyzing Figure 6, which shows the mean axial pressure and solid volume fraction profiles in the riser. The time-averaged volume fraction profile of Figure 6(b) was calculated on several consecutive cross-section planes along the height of the riser, according to the following expression:

786982.fig.005
Figure 5: Pressure drop in the riser as given by the semiempirical model and the CFD simulations.
fig6
Figure 6: Time-averaged CFD axial profiles for (a) pressure at the center line and (b) solid phase volume fraction calculated as an area average at several cross section planes along the height of the riser.
786982.fig.007
Figure 7: Instantaneous solid-phase volume fraction at 8, 9, and 10 seconds and the mean time average on an axial plane of the riser.

The CFD simulation indicates that the voidage fraction is not strongly segregated into a bottom dense region and an upper dilute one, but these regions are still clearly identified. Higher solid volume fraction zones were observed to occur close to the recirculation inlet height and at the T-shaped exit. This differs from the semiempirical model assumption of a much denser region at the bottom, a mean volume fraction of 0.1 against 0.04 from the simulation, thus explaining the lower pressure at the base of the riser.

Although the CFD model showed only a small difference between a denser region at the bottom and a more dilute above, the core-annulus behavior of the gas-solid flow found in the fast fluidization regime was verified. This configuration is commonly found in fast fluidization gas-solid flow in risers [19]. In this pattern a dilute core flows upwards and a denser ring falls near the walls. The core-annulus is more accentuated in the upper-half region of the riser, where the influence of the recirculation entry is less intense. This tendency is also observed at the column base where the riser’s diameter decreases to prevent the solids from falling further down and out off the loop, Figure 7.

Karri and Knowlton [42] remarked that the annulus structure can move cocurrently or countercurrently depending on the gas velocity and particulate flow rate. In Figure 8 the velocity and volume fractions profiles of the solid phase at three different positions in the riser are shown. A mostly counter-current pattern was identified, indicating internal recirculation. This characteristic helps to build a more uniform temperature profile along the height of the riser, as vertical mixing of the particulate phase occurs.

fig8
Figure 8: CFD simulation solid phase mean radial profiles of (a) volume fraction; (b) axial velocity component at 1.5 m, 2.5 m, and 3.5 m.

Furthermore, the thermal behavior of the gas-solid flow along the riser was analyzed. Figure 9 provides the gas phase temperature profile in the axial direction, showing reasonable agreement with the experimental data.

786982.fig.009
Figure 9: Time-averaged gas phase temperature profile at the axial central position in the riser compared with the experimental data.

As described previously, fluidizing gas heated in a GLP burner is introduced by the main inlet flow. The system loses thermal energy through the walls and by heating the solids returning to the system through the L-valve. To account for the heat loss, the heat flux at the wall was estimated considering only the external natural convection for a vertical cylinder with a constant mean external wall temperature of 338 K and an atmospheric temperature of 303 K. The total amount of heat lost to the environment under this condition was around 175 W/m2, which was set as the thermal boundary condition at the walls. This approximation was considered more realistic than an adiabatic wall. From the base of the riser to the lateral inlet, a temperature reduction was observed due to the return of the recirculated particles and the injection of auxiliary gas. In that region, the simulated results showed some deviation from the experimental measurement. However, the difference was lower than 22 K for the second and third points. Above that region, where the flow is more developed, a better agreement was found.

5. Conclusions

A practical semiempirical model based on hydrodynamic correlations from experimental data was proposed for determination of the main characteristics of a bench-scale CFB unit. This model predicts the solids inventory in components of the system, thereby providing important information for more detailed simulations through a CFD approach. The semiempirical model was validated with experimental information obtained from a bench-scale CFB unit operating under precombustion conditions. Deviation for key operational parameters, such as superficial gas velocity, solids inventory, and aeration mass flow rate in the recirculating valve were below 30%.

CFD simulations of the riser of the CFB unit were also carried out. In these simulations the realizable k-epsilon model described the turbulence of the continuous phase and the KTGF was applied to the description of the dispersed phase flow in the two-fluid model of the riser section of the CFB. The hydrodynamics verified the core-annulus model, which is also confirmed by most of the studies found in literature on fast fluidization beds. The solid circulation rate showed good agreement with the semiempirical model predictions. The simulated thermal profile in the riser showed small underpredicted values with respect to the experimental data, not larger then 3%.

The study also showed the viability of applying the semiempirical correlations to determine the solids inventory as an initial condition to a CFD simulation. By this approach less computational time is required to attain the pseudo-steady state than it would be to initialize the same simulation with an empty riser and constant . Because the semiempirical model was based on operational data of several CFB units, some deviation will always be found in the values predicted for a specific unit. Nevertheless, in the absence of experimental information, the approach presented in this paper can be a useful alternative approximation.

Nomenclature

Decay constant
:Area,
Ar:Archimedes number
:Solid-to-gas mass ratio at the entrance of the cyclone
Cd:Drag coefficient
:Diameter, m
:Friction coefficient
:Gravity, m/
:Radial distribution
:Circulation mass flux, kg/.s
:Height of the secondary air injection, m
:Total height of the riser, m
:Inventory, kg
:Thermal conductivity, W/m.K
:Model factor
:Length, m
:Mass rate, kg/h
Nu:Nusselt number
:Pressure, Pa
:Internal recirculation rate
Re:Reynolds number
:Time, s
:Temperature, K
Velocity, m/s
:Velocity, m/s
:Volume, .
Greek Letters
Model parameter
Model parameter
:Interphase momentum exchange coefficient, kg/·s
:Volume fraction
:Sphericity
:Solid bulk viscosity, P·s
:Shear viscosity, Pa·s
:Granular temperature, /
:Density, kg/
:Shear stress, Pa
:Velocity, m/s
:Stress tensor, Pa
:Solid stress tensor, Pa.
Subscripts
:Cyclone
ch:Choking
:Entry
:Falling
:Gas phase
:Phase index
lv:L-valve
mf:Minimum fluidization
ns:Nonspherical
:Outlet
:Riser
:Solid phase
sp:Standpipe
sr:Solids return level
:Terminal
tr:Transport
:Wall.

Acknowledgment

This work was supported by the national agency CNPq.

References

  1. E. Hartge, Y. Li, and J. Werther, “Flow structure in fast fluidized beds,” in Fluidization V, K. Ostergaard and A. Sorensen, Eds., pp. 345–352, Engineering Foundation, New York, NY, USA, 1986.
  2. P. Basu, Combustion and Gasification in Fluidized Beds, Taylor & Francis, Boca Raton, Fla, USA, 2006.
  3. T. Knowlton, “Nonmechanical solids feed and recycle devices for circulating fluidized beds,” in Circulating Fluidized Bed Technology II, P. Basu and J. Large, Eds., pp. 31–41, Pergamon Press, Oxford, UK, 1988.
  4. X. Wang, M. Rhodes, and B. Gibbs, “Solids flux distribution in a CFB riser operating at elevated temperatures,” in Fluidization VIII, F. Large and C. Laguérie, Eds., pp. 236–244, Engineering Foundation, New York, NY, USA, 1996.
  5. J. F. Davidson, “Circulating fluidised bed hydrodynamics,” Powder Technology, vol. 113, no. 3, pp. 249–260, 2000. View at Publisher · View at Google Scholar · View at Scopus
  6. H. Tong, H. Li, X. Lu, and Q. Zheng, “Hydrodynamic modeling of the L-valve,” Powder Technology, vol. 129, no. 1–3, pp. 8–14, 2003. View at Publisher · View at Google Scholar · View at Scopus
  7. A. Gungor and N. Eskin, “Effects of operational parameters on emission performance and combustion efficiency in small-scale CFBCs,” Journal of the Chinese Institute of Chemical Engineers, vol. 39, no. 6, pp. 541–556, 2008. View at Publisher · View at Google Scholar · View at Scopus
  8. M. J. Rhodes and P. Laussmann, “Study of the pressure balance around the loop of a circulating fluidized bed,” Canadian Journal of Chemical Engineering, vol. 70, no. 4, pp. 625–630, 1992. View at Scopus
  9. M. L. Mastellone and U. Aren, “The effect of particle size and density on solids distribution along the riser of a circulating fluidized bed,” Chemical Engineering Science, vol. 54, no. 22, pp. 5383–5391, 1999. View at Scopus
  10. S. W. Kim and S. D. Kim, “Effects of particle properties on solids recycle in loop-seal of a circulating fluidized bed,” Powder Technology, vol. 124, no. 1-2, pp. 76–84, 2002. View at Publisher · View at Google Scholar · View at Scopus
  11. X. Qi, J. Zhu, and W. Huang, “Hydrodynamic similarity in circulating fluidized bed risers,” Chemical Engineering Science, vol. 63, no. 23, pp. 5613–5625, 2008. View at Publisher · View at Google Scholar · View at Scopus
  12. Y. Chong, D. O'Dea, L. Leung, and D. Nicklin, “Design of standpipe and non-mechanical V-valve for a circulating fluidized bed,” in Circulating Fluidized Bed Technology II, P. Basu and J. Large, Eds., pp. 493–500, Pergamon Press, Oxford, UK, 1988.
  13. X. L. Yin, C. Z. Wu, S. P. Zheng, and Y. Chen, “Design and operation of a CFB gasification and power generation system for rice husk,” Biomass and Bioenergy, vol. 23, no. 3, pp. 181–187, 2002. View at Publisher · View at Google Scholar · View at Scopus
  14. R. Dewil, J. Baeyens, and B. Caerts, “CFB cyclones at high temperature: operational results and design assessment,” Particuology, vol. 6, no. 3, pp. 149–156, 2008. View at Publisher · View at Google Scholar · View at Scopus
  15. G. C. Lopes, L. M. Rosa, M. Mori, J. R. Nunhez, and W. P. Martignoni, “Three-dimensional modeling of fluid catalytic cracking industrial riser flow and reactions,” Computers and Chemical Engineering, vol. 35, pp. 2159–2168, 2011. View at Publisher · View at Google Scholar · View at Scopus
  16. L. Yu, J. Lu, X. Zhang, and S. Zhang, “Numerical simulation of the bubbling fluidized bed coal gasification by the kinetic theory of granular flow (KTGF),” Fuel, vol. 86, no. 5-6, pp. 722–734, 2007. View at Publisher · View at Google Scholar · View at Scopus
  17. D. Gidaspow, Multiphase Flow and Fluidization, Academic Press, London, UK, 1994.
  18. B. G. M. van Wachen, Derivation, implementation and validation of computer simulation models for gas-solid fluidized beds [Ph.D. thesis], Delft University, 1998.
  19. S. Benyahia, M. Syamlal, and T. J. O'Brien, “Particle Technology and Fluidization. Study of the ability of multiphase continuum models to predict core-annulus flow,” AIChE Journal, vol. 53, no. 10, pp. 2549–2568, 2007. View at Publisher · View at Google Scholar · View at Scopus
  20. J. C. S. C. Bastos, L. M. Rosa, M. Mori, F. Marini, and W. P. Martignoni, “Modelling and simulation of a gas-solids dispersion flow in a high-flux circulating fluidized bed (HFCFB) riser,” Catalysis Today, vol. 130, no. 2–4, pp. 462–470, 2008. View at Publisher · View at Google Scholar · View at Scopus
  21. V. Jiradilok, D. Gidaspow, S. Damronglerd, W. J. Koves, and R. Mostofi, “Kinetic theory based CFD simulation of turbulent fluidization of FCC particles in a riser,” Chemical Engineering Science, vol. 61, no. 17, pp. 5544–5559, 2006. View at Publisher · View at Google Scholar · View at Scopus
  22. X. Lan, C. Xu, G. Wang, L. Wu, and J. Gao, “CFD modeling of gas-solid flow and cracking reaction in two-stage riser FCC reactors,” Chemical Engineering Science, vol. 64, no. 17, pp. 3847–3858, 2009. View at Publisher · View at Google Scholar · View at Scopus
  23. X. Wang, B. Jin, W. Zhong, and R. Xiao, “Modeling on the hydrodynamics of a high-flux circulating fluidized bed with geldart group a particles by kinetic theory of granular flow,” Energy and Fuels, vol. 24, no. 2, pp. 1242–1259, 2010. View at Publisher · View at Google Scholar · View at Scopus
  24. A. Samuelsberg and B. H. Hjertager, “An experimental and numerical study of flow patterns in a circulating fluidized bed reactor,” International Journal of Multiphase Flow, vol. 22, no. 3, pp. 575–591, 1996. View at Publisher · View at Google Scholar · View at Scopus
  25. J. J. Ramirez-Behainne, Assessment of mercury emissions from Brazilian coal combustion in fast fluidized bed [Ph.D. thesis], UNICAMP, Campinas, Brazil, 2007.
  26. W. C. Yang, “Modification and re-interpretation of Geldart’s classification of powders,” Powder Technology, vol. 171, no. 2, pp. 69–74, 2007.
  27. R. I. Hory, J. J. Ramirez-Behainne, A. A. B. Pécora, and L. Goldstein Jr, “An empirical model to predict the mass flow rate of solids in a high temperature circulating fluidized bed system,” in Proceedings of the 11th Brazilian Congress of Thermal Sciences and Engineering, ENCIT, 2006.
  28. C. Y. Wen and Y. U. Yu, “A generalized method for predicting the minimum fluidization velocity,” American Institute of Chemical Engineers Journal, vol. 12, pp. 610–612, 1966.
  29. D. Bai, Y. Jin, and Z. Yu, “Flow regimes in circulating fluidized beds,” Chemical Engineering and Technology, vol. 16, no. 5, pp. 307–313, 1993. View at Scopus
  30. D. Kunii and O. Levenspiel, “Entrainment of solids from fluidized beds I. Hold-up of solids in the freeboard II. Operation of fast fluidized beds,” Powder Technology, vol. 61, no. 2, pp. 193–206, 1990. View at Scopus
  31. D. Kunii and O. Levenspiel, “Flow modeling of fast fluidized beds,” in Circulating Fluidized Bed Technology III, P. Basu, M. Hasatani, and M. Horio, Eds., pp. 91–98, Pergamon Press, Oxford, UK, 1991.
  32. W. C. Yang, “Criteria for choking in vertical pneumatic conveying lines,” Powder Technology, vol. 35, no. 2, pp. 143–150, 1983. View at Scopus
  33. J. Perales, T. Coll, M. Llop, L. Puigjaner, J. Arnaldos, and J. Casal, “On the transition from bubbling to fast fluidization regimes,” in Circulating Fluidized Bed Technology III, P. Basu, M. Hasatani, and M. Horio, Eds., pp. 73–78, Pergamon Press, Oxford, UK, 1991.
  34. J. J. Ramirez-Behainne and G. Martins, “A semi-empirical model to calculate the inert solid inventory and the main dimensions of a CFB reactor aiming biomass thermochemical conversion,” in Proceedings of the 20th International Congress of Mechanical Engineering, 2009.
  35. E. Muschelknautz and V. Greif, “Cyclone and other gas-solid separators,” in Circulating Fluidized Beds, J. Grace, A. Avidan, and T. Knowlton, Eds., pp. 181–213, Blackie Academic & Professional, London, UK, 1997.
  36. T. Knowlton, “Standpipes and return systems,” in Circulating Fluidized Beds, J. Grace, A. Avidan, and T. Knowlton, Eds., pp. 214–260, Blackie Academic & Professional, London, UK, 1997.
  37. D. Geldart and P. Jones, “The behaviour of L-valves with granular powders,” Powder Technology, vol. 67, no. 2, pp. 163–174, 1991. View at Scopus
  38. S. Ogawa, A. Umemura, and N. Oshima, “On the equations of fully fluidized granular materials,” Journal of Applied Mathematics and Physics, vol. 31, no. 4, pp. 483–493, 1980. View at Publisher · View at Google Scholar · View at Scopus
  39. C. K. K. Lun, S. B. Savage, D. J. Jeffrey, and N. Chepurniy, “Kinetic theories for granular flow: inelastic particles in Couette flow and slightly inelastic particles in a general flowfield,” Journal of Fluid Mechanics, vol. 140, pp. 223–256, 1984. View at Scopus
  40. W. E. Ranz and W. R. Marshall Jr., “Evaporation and drops—part I and II,” Chemical Engineering Prog, vol. 48, pp. 173–180, 1952.
  41. S. Benyahia, H. Arastoopour, T. M. Knowlton, and H. Massah, “Simulation of particles and gas flow behavior in the riser section of a circulating fluidized bed using the kinetic theory approach for the particulate phase,” Powder Technology, vol. 112, no. 1-2, pp. 24–33, 2000. View at Publisher · View at Google Scholar · View at Scopus
  42. S. B. R. Karri and T. M. Knowlton, “A comparison of annulus solids flow direction and radial solids mass flux profiles at low and high mass fluxes in a riser,” in Proceedings of the 6th International Conference on Circulating Fluidized Beds, pp. 71–76, 1999.