A gas-solid two-fluid model with the second-order moment method is presented to close the set of equations applied to fluidization. With the kinetic theory of granular flow, transport equations for the velocity moments are derived for the particle phase. Closure equations for the third-order moments of velocity and for the fluid-particle velocity correlation are presented. The former is based on a modified model with the contribution of the increase of the binary collision probability, and the latter uses an algebraic model proposed by Koch and Sangani [1999. Particle pressure and marginal stability limits for a homogeneous monodisperse gas-fluidized bed: kinetic theory and numerical simulations. journal of Fluid Mechanics 400, 229-263]. Boundary conditions for the set of equations describing flow of particles proposed by Strumendo and Canu [2002. Method of moments for the dilute granular flow of inelastic spheres. Physical Review E 66, 041304/1-041304/20] are modified with the consideration of the momentum exchange by collisions between the wall and particles. Flow behavior of gas and particles is performed by means of gas-solid two-fluid model with the second-order moment model of particles in the bubbling fluidized bed. The distributions of velocity and moments of particles are predicted in the bubbling fluidized bed. Predictions are compared with experimental data measured by Muller et al. [2008. Granular temperature: comparison of magnetic resonance measurements with discrete element model simulations. Powder Technology 184, 241-253] and Yuu et al. [2000. Numerical simulation of air and particle motions in bubbling fluidized bed of small particles. Powder Technology 110, 158-168]. in the bubbling fluidized beds. The simulated second-order moment in the vertical direction is 1.1-2.5 [Muller, C.R., Holland, D.J., Sedeman, A.J., Scott, S.A., Dennis, J.S., Gladden. L.F.. 2008. Granular temperature: comparison of magnetic resonance measurements with discrete element model simulations. Powder Technology 184, 241-253] and 1.1-4.0 [Yuu, S., Umekage, T., Johno, Y., 2000. Numerical simulation of air and particle motions in bubbling fluidized bed of small particles. Powder Technology 110, 158-168] times larger than that in the lateral direction because of higher velocity fluctuations for particles in the bubble fluidized bed. The bubblelike Reynolds normal stresses per unit bulk density used by Gidaspow et al. [2004. Hydrodynamics of fluidization using kinetic theory: an emerging paradigm 2002 Flour-Daniel lecture. Powder Technology 148, 123-141.] are computed from the simulated hydrodynamic velocities. The predictions are in agreement with experimental second-order moments measured by Muller et al. [2008. Granular temperature: comparison of magnetic resonance measurements with discrete element model simulations. Powder Technology 184, 241-253] and fluctuating velocity of particles measured by Yuu et al. [2000. Numerical simulation of air and particle motions in bubbling fluidized bed of small particles. Powder Technology 110, 158-168]. (C) 2009 Elsevier Ltd. All rights reserved.