- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Table of Contents
Advances in Mechanical Engineering
Volume 2013 (2013), Article ID 616435, 14 pages
Particle Dispersion Behaviors of Dense Gas-Particle Flows in Bubble Fluidized Bed
1College of Chemistry and Environmental Engineering, Dongguan University of Technology, Dongguan 523808, China
2School of Electronic and Information Engineering, Dalian Jiaotong University, Dalian 116028, China
3Performance Test Center, Harbin Power System Engineering and Research Institute, Harbin 150046, China
Received 11 January 2013; Revised 26 April 2013; Accepted 20 May 2013
Academic Editor: Guan Heng Yeoh
Copyright © 2013 Sihao Lv 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.
An Euler-Euler two-fluid model incorporating a developed momentum transfer empirical coefficient is developed to study the particle dispersion behaviors of dense gas-particle flows in gas-fluidization reactor. In this model, the four-way couplings among gas-particles, particle-gas, and particle-particle collisions are fully considered based on kinetic theory of granular flows and an improved smooth continuous drag coefficient is utilized. Gas turbulent flow is solved by large eddy simulation. The particle fraction, the time-averaged axial particle velocity, the histogram of particle fluctuation velocity, and the wavelet analysis of pressure signals are obtained. The results are in good agreement with experimental measurements. The mean value and the variance of axial particle velocity are greater than those of radial particle velocities. Particle collision frequencies at bubble vibrant movement regions along axial direction are much higher than those of radial direction and attenuated along height increase. Low-frequency component of pressure signal indicating the bubble movement behaviors in the center of reactor is stronger than wall regions. Furthermore, the negative values represent the passed bubble and positive peak values disclose the continuous motion of single bubble.
Circulating fluidized beds (CFB) have been widely applied in the chemical industry due to their practical applications, such as riser, downer, and gas-fluidized bubble bed . The efficiency of systems relies primarily on their hydrodynamic behaviors, especially for particle dispersion behaviors. Thus, a better understanding of the dynamics of fluidized beds is the key problem to improve the efficiency and scale-up laboratory, which can be achieved through numerical simulation due to the rapid development of computer hardware technique. Except the traditional experimental measurement approach, the computational fluid dynamics (CFD) a powerful tool for the better understanding of physical phenomena involved in the gas-particle fluidization system [2, 3].
Generally, mathematical simulation of gas-particle bubbling fluidized bed can be divided into the Euler-Euler two-fluid model and Euler-Lagrange discrete particle model. In the Eulerian-Lagrangian approach, numerous discrete particles are tracked and interparticle collisions are simulated using a hard-sphere or a soft-sphere model. In the framework of two-fluid model, particles and gas are described as interpenetrating continua and the two sets of conservation equations for mass, momentum, and energy are numerically solved. The Eulerian-Eulerian models require closure laws for particle interactions, interaction between gas and particle phases to close the system of equations, which are generally a set of models, either physical or empirical correlations. Comparing these two approaches, due to the huge computation consumption for discrete particle model in industry application, two-fluid model has successfully applied. So far, the most popular theory for closing particle-particle collision is the kinetic theory of granular flow (KTGF), firstly proposed by Gidaspow [4, 5]. It is analogous with the kinetic theory of gases and the macroscopic behavior of particle phase is described by the equations that account for the energy associated with particles arising out of collisions and fluctuating motions of the particles. The kinetic and collision momentum transfer due to particle collisions can be modeled as a function of the local particle concentration and the fluctuating motion of the particles. Ding and Gidaspow 1990  firstly proposed the concept of granular temperature to indicate particle collisions, which is defined as , where is the random fluctuating velocity of the particulate phase. Thus, the rheologic properties of particles, the solid phase pressure, viscosity of particles, conductive energy flux, and dissipation are considered as a function of granular temperature. Simulations using kinetic theory of granular flow have been reported by many researchers in gas-particle fluidization [7–20]. In KFGT, the particle flow is assumed to be laminar flows instead of large-scale turbulent flows. In order to solve this point, Zhou et al. [21–26] research group of Tsinghua University have successfully proposed a series of kinetic energy equations for gas and particle phases , second-order and unified second-order moment (SOM, USM) models for dilute and dense gas-particle turbulent flows, that is, model, four-equation model, five-equation model, USM- particle temperature model, and subgrid scale USM model to simulate particle turbulent flows based on the Reynolds-averaged Navier-Stokes (RANS) method. In USM model, the anisotropy between gas and particle two-phase stresses and their interaction can be fully considered by establishing the two-phase Reynolds stresses transport equations and their stress correlations.
Besides the effects of particle and particle collision, the other very important factor on gas-particle two-phase turbulent flow behavior is the momentum transfer process between gas and particle phases. This kind of transfer behavior can be typically represented by drag force. Due to the fact that it is difficult to get an accurate value restricted to different Reynolds numbers and packing fractions, as well as in terms of homogeneity, monodispersity, and sphericity of the particles, a large number of the gas-solid drag forces have some limitations due to empirical relations from the measurments. Although the relations by Ergun and Orning , Wen and Yu , and Giblaro et al.  have been the most widely used, there is at present no real consensus as to what the most accurate prediction for the drag force is at given Reynolds numbers and packing fraction. Thus, so far, all correlations are all based on experimental data. In order to assure these correlations more accurately in theory, Huilin et al. 2003  introduced a weighted average two-scale switch function to prevent the discontinuous behaviors when solid volume fraction is less than 0.2. Beetstra et al. 2007  proposed a drag force model based on the kinetic theory of granular flows and lattice Boltzmann data for monodisperse and bidisperse systems. Even if these corrections have successfully got better results in specified cases, they should be further validated for wide applications.
As we have known, the motion of bubbles can produce the pressure fluctuations and the amplitude of the pressure fluctuations is definitely related to the particle dispersion behaviors. Successfully, wavelet analysis approach has identified the multiscale frequency wave components for bubble upward and downward movements with different particle dispersion behaviors. Van Der Schaaf 1998  demonstrated that the pressure fluctuations in a gas-solid fluidized bed could be partitioned into a fractional Brownian motion (FBM) and Gaussian white noise (GWN), caused by the jet and the formation of the small bubbles near the distributor. He et al. 1997  indicated that the upward moving compression waves originate from the formation and coalesce of gas bubbles, and downward moving compression waves are caused by gas bubble eruptions at the fluidized bed surface. Lu and Li 1999  found that the peak frequency of these signals can represent the bubbling frequency, and the average peak value can represent the bubble size. Although many models have been developed to describe these behaviors, a valid method to process the pressure fluctuation signals is to be further developed and validated.
To date, considering the effects of gas-particle moment transfer represented by transfer coefficient coupling KGFT theory and wavelet analysis of pressure fluctuation approaches on dense particle dispersion behaviors in bubble fluidized bed have never been reported. Compared to RANS approach using LES, the instantaneous flow information of two phases and coherent structure of particle dispersion only can be provided in detail rather than the RANS. In this study, coupled with LES, KGFT, and improved transfer coefficient model, as well as wavelet analysis approach, are used to study the particle hydrodynamics and transport characteristics.
2. Conservation Equations of Two-Phase Turbulent Flows
2.1. Governing Equations
The continuity equations, representing the mass conservation for gas and solid phases, respectively, can be written as follows: where is the velocities, is densities, and is the volumetric fractions of gas and solid phases, respectively. Subscript is gas phase and is solid phase. In relation to single-phase continuity equation, the above equations differ by the presence of the phases volumetric fractions. It is defined as follows: The momentum equations for gas and solid phases, respectively, can be written as follows: In (3), is the thermodynamic pressure, is the interface drag coefficient between gas phase and particle phase, is the gas phase viscous stress tensor, is the gravity source term of gas phase, is the dynamic viscosity and assumed to be constant and is set equal to , is gravitational accelerations, is the unit tensor, subscript lam is the laminar flow, and is the turbulence flow.
2.2. Closure for Stress Equations
The Newtonian rheology and the Stokes hypothesis are assumed for both gas phase and particle phase; the viscous stress tensor of gas and particle phase are computed as follows:
The gas is assumed to be incompressible and has a constant density. Gas turbulent flow is solved by large eddy simulation; that is, the larger-scale eddy is directly solved by Navier-Stokes equations and the smaller-scale eddy is modeled by subgrid scale model. The traditional Smagorinsky’s  eddy viscosity has a problem that there is an overlarge eddy viscosity dissipation at wall. In this paper, empirical corrected model presented by Yuu et al. 2000  is utilized; it is as follows:
2.3. Particle Temperature Equations
The conservation equation of particle fluctuating energy or translational granular temperature is given as follows: where is the conductivity coefficient of granular temperature: The translational fluctuation energy dissipation rate is The particle pressure is
where is the particle temperature and is the restitution coefficient of particles.
The particle dynamic viscosity is simulated as follows: where is the particle diameter and is the radial distribution function at contact particle:
where is the particle maximum volume fraction at random packing.
The particle bulk viscosity with the effective coefficient of restitution is computed as follows:
2.4. Momentum Transfer between Gas and Particle Phases
In order to couple the momentum transfer between gas and particle phases, a model for the drag force is required. Ergun and Orning , Wen and Yu , and Giblaro et al.  have presented some correlations based on the experimental results, but there still are some limitations in theory. Although Beetstra et al. 2007  established a drag force model based on the kinetic theory of granular flows and lattice Boltzmann data for monodisperse and bidisperse systems, it still should be further validated for applications. In this paper, based on the Huilin et al. 2003  smooth function when solid volume fraction is less than 0.2 and the McKeen and Pugsley 2003  experimental results, a modified momentum transfer coefficient is presented. where is the interface drag coefficient for a single particle in an infinite medium, calculated by: The Reynolds number Re based on the particle mean diameter is a function of the particle sphericity : The final sooth function is
3. Numerical Strategy, Experiments Setup and Boundary Conditions
3.1. Numerical Strategy
The aforementioned governing equations are discretized using the finite differences approach with second-order upwind scheme and solved by point relaxation technique. Continuity equations are in the form of implicit format and momentum equations are discretized over a staggered mesh. In the momentum equations, the convective terms are treated explicitly and all other terms are treated implicitly. The two-dimensional computational mesh is uniform staggered grid. Scalar quantity parameters such as the volume fraction of particles, density, and turbulent kinetic energy are stored at the main grid points. Staggered grid scheme is used and the velocity components are solved at the volume surfaces. A guessed pressure field is initialized for simulation start and the velocities are calculated from the momentum equations. The particulate phase continuity equation is solved using the updated velocities to compute the particulate phase volume fraction. The gas phase volume fraction is then computed. By using the gas volume fraction, the updated velocities the gas phase mass residue from continuity equation for gas phase is computed, which is used as a convergence criterion. For convergence, the gas pressure is corrected in each cell at a time until convergence is attained or the number of iterations exceeds an inner iteration limit (a number of iterations for one cell at a given time). The computations proceed until the entire computational domain is covered. At the end of such a computational sweep, if a pressure adjustment was necessary in any of the cells, the procedure is repeated until simultaneous convergence in all the cells is obtained.
3.2. Experimental Setup and Boundary Conditions
Experimental setup is chosen from Peltola et al. 2010 . The detailed geometry dimension, gas and particle property parameters are seen in Figure 1 and Table 1, respectively. Height of initial particle packing is set equal to 110 mm. Values of the volume fractions of particles at the inlet are set to 0.4 and of initial jet velocity to 27.5 m/s. At the outlet, the pressure was set to 1 atm. Cartesian coordinate system is used with 40 cells uniform computational mesh in radial direction and 150 cells uniform in axial direction. Cell flag and characters computational cells are shown in Table 2 and Figure 2.
3.3. Boundary Conditions
At the wall, no slip condition is applied, in which the gas tangential and normal velocities were set to zero, and the normal velocity of particles was also set at zero. Pressure boundary condition is set to the zero for the normal gradient of pressure at the inlet and outlet. Velocity boundary condition is to provide the superficial gas velocity at the inlet cell and used to the free slip condition at the outlet profile. For temporal integration was used a fixed time step equal to s based on the CFL stable criteria and the simulation running time is 25 seconds. The turbulent kinetic energy and its dissipation rate are taken by empirical expressions. The fully developed flow conditions of two-phase are taken at the outlet. At the wall, no slip condition is used for gas phase velocity and gas Reynolds stress are determined via production term including the effect of wall function for near wall grid nodes. . At the near-wall grid nodes, the wall-function approximation is used. At the axis, symmetric conditions are adopted for both the two phase. The convergence criteria for gas phase are mass source .
4. Results and Discussion
In order to assure the accuracy of the simulation, the grid independence test for sensitivity of computed axial particle velocity distributions to spatial grid sizes is carried out. Here, coarse (75 × 20 axial by radial), medium , and fine mesh grids are used in simulation process. Results indicated that the velocity distributions have similar trends along both axial directions, especially for medium and fine grids. So, we choose the medium grid size scheme for the simulations in order to achieve good accuracy within manageable computation time. In the meantime, the validation with experimental data also carried out for the statistical time-averaged axial particle velocity, which is defined by the dimensionless time between 100 and 500, is determined as the sample of instantaneous velocity and is expressed by The simulated results of the 0.8 m height are validated by experimental data (see Figure 4). We can see that the predicted results are in good agreement with experimental results from the literature . So, this model and code in this paper have been accurately validated.
Figure 3 shows the bubbles movement behaviors in bubble gas-fluidized bed at 0.5 s, 2.0 s, 4.0 s, 8.0 s, 20.0, and 25.0 s, respectively. The bubble can be defined by the instantaneous local high porosity of particles. Two kinds of bubble behaviors named the faster and the slower are identified as the halation layer around bubble, which is composed of motion circulation gas between bubble and environments. It would be affected by the bubble rising velocity. When bubble rising velocity is greater than gas gap velocity between bubble and environment, it is called the halation bubble. However, if bubble rising velocity is too fast to lead to the decrease of halation layer, the most of gas has to circulate inside bubble. When gas gap velocity is larger than bubble rising, the gas in halation layer will penetrate the bubble; it is called non halation bubble.
When the higher jet gas passed through the interior reactor, the stronger interaction between gas and particle phases occurred. Then, particles located at the jet inlet position are pushed upward. Thus, the local porosity will rapidly increase and the whole system begins to form the unstable stage in the gas-particle system. Once the stable system was lost, the rapidly increasing porosity is produced. At the bottom of bubble, the flow pressure is less than other positions of bubble. It leads to drive particles around bubble entered into the inner space of bubble. Thus, bubble had to leave the jet position and begin to move toward bed top surface and to grow gradually. Result of this kind of motion and development is to produce the individual bubble. When it reaches up to the top surface of bed, it will be broken up. Becaues of motion effects of bubble, they lead to the drastic collision of particles and accelerates the mixed process of particles. So, particle dispersions gradually become vigorous and the particle dispersions of gas-particle system shows on the intensely heterogeneous. From these figures, they clearly show the formation, the movements and the breakup of the bubbles, the movement and circulation of particle flows and dispersion and the interaction of gas and particle phases. Some useful information, is higher particle temperature due to particle collisions, the higher porosity due to bubble movement and two-phase turbulent behaviors, and so forth.
Figures 5, 6, and 7 show the distributions of instantaneous concentration of particles at three different heights of 0.4 m, 0.8 m, and 1.2 m and at of 0.25 and 0.73, respectively. The distribution of the computed concentration of particles shows uniformity in the bubbling fluidized bed, except for low values in the bubble burst region. The concentration of particles is high near the walls and in the center regime of the bed. It is low between the center and the walls. Near the bed surface, the concentrations of particles become low in the center of the bed. So, the uniformity and nonperiodicity are the main characters of particle dispersions.
Figures 8, 9, and 10 show the distributions of particles velocity histogram at axial and radial directions at three different heights of 0.4 m, 0.8 m, and 1.2 m and at of 0.25, 0.47, and 0.73, respectively. We can see the distributions of the mathematical expectation and the variance of particle velocity in radial and axial direction. Mathematical expectation represents the mean value of the random variables and the variance is the fluctuation intensity of random variables. At the height of 0.4 m, the mean values at axial direction are both less than radial directions with different radial positions. But the variance values at axial direction are larger than radial directions. With respect to the those of heights of 0.8 m and 1.2 m, they have the same trend compared to those of 0.4 m. As for the axial direction at the three heights, the mean and the variance at of 0.47 and 0.72 and the middle of the 0.8 m have largest values, while at the bottom height of 0.4 m, they have the largest values at of 0.25. As for the radial direction at three heights, the mean and the variance at all of of 0.47 and 0.72 have the larger values than those of at of 0.25. So, the larger particle collisions are often formed at the more uniform gas-particle fluid system.
To further understand the nonlinear particle dispersion behaviors, nonuniform and multiscale turbulent flows are very important to optimize design and operation for gas-fluidized bed. To great extents, the complexity degrees are all represented by particle concentration and pressure fluctuations. Figures 11, 12, and 13 show the distributions of the wavelet analysis of pressure fluctuations at axial and radial directions at three different heights of 0.4 m, 0.8 m, and 1.2 m and at of 0.25, 0.47, and 0.73, respectively. The primary pressure signal are provided in these figures, where stand for the detailed information of primary signals and are the approximation of signals. Results showed that the higher frequent components gradually disappeared with increasing scale and the lower frequent components are finally lost, which indicated that wavelet analysis could be stopped at the 8th layer. Pressure drops will be affected due to the variation of particle concentration or gas velocity. In the gas particle bubble fluidized bed, pressure drop have different frequencies for time behaviors. The higher frequency is mainly caused by the particle fluctuation and the lower frequency is bubble movements. In other word particle fluctuation can be considered as the higher-faster frequent components and bubble is the lower-slower frequent components.
As seen from these figures, regarding the same radial positions at three heights, the components fluctuations at the scale of the d1, the d2, and the d3 at height of 0.4 m will last a longer time than those of at 0.8 m and 1.2 m, in which small scale with the higher frequent component is a function of time. Particle fluctuation is a kind of the higher frequency and bubble is a kind of the lower frequency. It indicated that the random particle fluctuation will be attenuated with the height increase. As for the different radial positions at the same height, the component signal at d7 is obvious, where the peak values are larger than others. In addition, larger scale shows that the lower frequency of pressure fluctuation will be changed with the time development and the higher frequency of pressure fluctuation bubble reflected the bubble movements. At near-wall region, the bubble movements in the center bed are more vigorous. The variation of peak means the alternative formation of bubble and particle movement, where negative peaks represent the bubble passing and positive peaks are the continuous passing of the single bubble.
The instantaneous and statistical particle dispersions, the nonlinear hydrodynamic characteristic, multiscale flow regions, and heterogeneous characteristics in complex dense gas-particle flow system are obtained using developed Euler-Euler two-fluid model, an improved momentum transfer empirical coefficient, and wavelet analysis methods, as well as LES for solving the gas turbulent flows, which can be considered to be a novelty research approach for this study. The most important conclusions are as follows.(1)Particle collision frequencies at bubble vibrant movement regions along axial direction are much higher than those of radial direction and attenuated along height increase.(2)Representing bubble movement, low frequency components of pressure signal in the center of reactor regions are stronger than wall regions.(3)In wavelet analysis results, the negative values represent the passed bubble and positive peaks are the continuous motion of single bubbles.
|:||Particle restitution coefficient|
|:||Effective particle restitution coefficient|
|:||Radial distribution function|
|:||Kronic-delta unit tensor|
|, , , :||Coordinates directions|
The authors sincerely appreciate the financial support of Projects of National Natural Science Foundation of China under the Grant no. 20906014.
- J. R. Grace, “Reflections on turbulent fluidization and dense suspension upflow,” Powder Technology, vol. 113, no. 3, pp. 242–248, 2000.
- S. Elghobashi, “On predicting particle-laden turbulent flows,” Applied Scientific Research, vol. 52, no. 4, pp. 309–329, 1994.
- C. T. Crowe, M. Sommerfeld, and Y. Tsuji, Multiphase Flows with Droplets and Particles, CRC Press, Boca Raton, Fla, USA, 1998.
- D. Gidaspow, Multiphase Flows and Fluidization, Academic Press, San Diego, Calif, USA, 1994.
- S. Chapman and T. G. Cowling, The Mathematical Theory of Non-Uniform Gases, Cambridge University Press, 1970.
- J. Ding and D. Gidaspow, “Bubbling fluidization model using kinetic theory of granular flow,” AIChE Journal, vol. 36, no. 4, pp. 523–538, 1990.
- J. L. Sinclair and R. Jackson, “Gas-particle flow in a vertical pipe with particle-particle interactions,” AIChE Journal, vol. 35, no. 9, pp. 1473–1486, 1989.
- S. Dasgupta, R. Jackson, and S. Sundaresan, “Turbulent gas-particle flow in vertical risers,” AIChE Journal, vol. 40, no. 2, pp. 215–228, 1994.
- C. M. Hrenya and J. L. Sinclair, “Effects of particle-phase turbulence in gas-solid fows,” AIChE Journal, vol. 43, no. 4, pp. 853–869, 1997.
- V. Mathiesen, T. Solberg, H. Arastoopour, and B. H. Hjertager, “Experimental and computational study of multiphase gas/particle flow in a CFB riser,” AIChE Journal, vol. 45, no. 12, pp. 2503–2518, 1999.
- G. Balzer, “Gas-solid flow modelling based on the kinetic theory of granular media: validation, applications and limitations,” Powder Technology, vol. 113, no. 3, pp. 299–309, 2000.
- 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.
- B. G. M. Van Wachem, J. C. Schouten, C. M. Van den Bleek, R. Krishna, and J. L. Sinclair, “Comparative analysis of CFD models of dense gas-solid systems,” AIChE Journal, vol. 47, no. 5, pp. 1035–1051, 2001.
- K. Agrawal, P. N. Loezos, M. Syamlal, and S. Sundaresan, “The role of meso-scale structures in rapid gas-solid flows,” Journal of Fluid Mechanics, vol. 445, pp. 151–181, 2001.
- M. J. V. Goldschmidt, J. A. M. Kuipers, and W. P. M. Van Swaaij, “Hydrodynamic modelling of dense gas-fluidised beds using the kinetic theory of granular flow: effect of coefficient of restitution on bed dynamics,” Chemical Engineering Science, vol. 56, no. 2, pp. 571–578, 2001.
- J. De Wilde, G. B. Marin, and G. J. Heynderickx, “The effects of abrupt T-outlets in a riser: 3D simulation using the kinetic theory of granular flow,” Chemical Engineering Science, vol. 58, no. 3-6, pp. 877–885, 2003.
- L. Huilin, D. Gidaspow, J. Bouillard, and L. Wentie, “Hydrodynamic simulation of gas-solid flow in a riser using kinetic theory of granular flow,” Chemical Engineering Journal, vol. 95, no. 1, pp. 1–13, 2003.
- J. Leboreiro, G. G. Joseph, and C. M. Hrenya, “Revisiting the standard drag law for bubbling, gas-fluidized beds,” Powder Technology, vol. 183, no. 3, pp. 385–400, 2008.
- N. Reuge, L. Cadoret, C. Coufort-Saudejaud, S. Pannala, M. Syamlal, and B. Caussat, “Multifluid Eulerian modeling of dense gas-solids fluidized bed hydrodynamics: influence of the dissipation parameters,” Chemical Engineering Science, vol. 63, no. 22, pp. 5540–5551, 2008.
- A. Passalacqua and L. Marmo, “A critical comparison of frictional stress models applied to the simulation of bubbling fluidized beds,” Chemical Engineering Science, vol. 64, no. 12, pp. 2795–2806, 2009.
- L. X. Zhou, Theory and Numerical Simulation Modeling of Turbulent Gas-Particle Flows and Combustion, Science Press and Florida, CRC Press, Beijing, China, 1993.
- L. X. Zhou, C. M. Liao, and T. Chen, “Unified second-order-moment two-phase turbulence model for simulating gas-particle flows,” in Proceedings of the ASME Fluids Engineering Division Summer Meeting, vol. 185, pp. 307–313, June 1994.
- L. X. Zhou and T. Chen, “Simulation of swirling gas-particle flows using USM and k-ε-kp two-phase turbulence models,” Powder Technology, vol. 114, no. 1–3, pp. 1–11, 2001.
- Y. Yu, L. X. Zhou, B. T. Wang, and F. Cai, “USM-Θ two-phase turbulence model for simulating dense gas-particle flows,” Acta Mechanica Sinica/Lixue Xuebao, vol. 21, no. 3, pp. 228–234, 2005.
- Y. Liu, L. X. Zhou, and C. X. Xu, “Numerical simulation of instantaneous flow structure of swirling and non-swirling coaxial-jet particle-laden turbulence flows,” Physica A, vol. 389, no. 23, pp. 5380–5389, 2010.
- Y. Liu, X. Liu, S. Kallio, and L. X. Zhou, “Hydrodynamic predictions of dense gas-particle flows using a second-order-moment frictional stress model,” Advanced Powder Technology, vol. 22, no. 4, pp. 504–511, 2011.
- S. Ergun and A. A. Orning, “Fluid flow through packed columns,” Chemical. Engineering Progress, vol. 48, pp. 89–94, 1952.
- C. Y. Wen and Y. H. Yu, “Mechanics of fluidization,” Chemical Engineering Progress Symposium, vol. 62, pp. 100–113, 1966.
- L. A. Giblaro, R. Difelice, and S. S. Waldram, “Generalized friction and drag coefficient correlations for flux particle interactions,” Chemical Engineering Science, vol. 40, pp. 1817–1823, 1985.
- R. Beetstra, M. A. Van Der Hoef, and J. A. M. Kuipers, “Drag force of intermediate reynolds number flow past mono- and bidisperse arrays of spheres,” AIChE Journal, vol. 53, no. 2, pp. 489–501, 2007.
- J. Van Der Schaaf, J. C. Schouten, and C. M. Van Den Bleek, “Origin, propagation and attenuation of pressure waves in gas-solid fluidized beds,” Powder Technology, vol. 95, no. 3, pp. 220–233, 1998.
- Z. He, W. Zhang, K. He, and B. Chen, “Modeling pressure fluctuations via correlation structure in a gas-solids fluidized bed,” AIChE Journal, vol. 43, no. 7, pp. 1914–1920, 1997.
- X. S. Lu and H. Z. Li, “Wavelet analysis of pressure fluctuation signals in a bubbling fluidized bed,” Chemical Engineering Journal, vol. 75, no. 2, pp. 113–119, 1999.
- J. Smagorinsky, “General circulation experiments with the primitive equation,” Journal of Chemical Physics, no. 91, pp. 99–164, 1963.
- S. Yuu, T. Umekage, and Y. Johno, “Numerical simulation of air and particle motions in bubbling fluidized bed of small particles,” Powder Technology, vol. 110, no. 1-2, pp. 158–168, 2000.
- T. McKeen and T. Pugsley, “Simulation and experimental validation of a freely bubbling bed of FCC catalyst,” Powder Technology, vol. 129, no. 1–3, pp. 139–152, 2003.
- J. Peltola, S. Kallio, and M. Honkanen P, “Saarenrinne. Image based measurement of particle phase Reynolds stresses in a laboratory scale circulating fluidized bed,” in Proceedings of the 7th International Conference on Multiphase Flow (ICMF '10), Tampa, Fla, USA, June2010.