A spouted bed is simulated in three dimensions by a discrete element method (DEM) in a cylindrical coordinate system. The numerical scheme is based on a second order finite difference method in space and a second order Adams-Bashforth method for time advancement. Gas-particle interaction is assumed to obey the Ergun equation (for void fraction less than 0.8) and its corrected model by Wen and Yu (for void fraction greater than 0.8). The spouted bed vessel is a flat-bottomed cylinder 0.5 m in height and 0.15 m in diameter. The gas inlet diameter is 0.02 m. Three hundred thousand monosized spheres of diameter 2.4 mm are used in the simulation. The typical characteristics of spouted beds, such as spout, annulus and fountain, are reproduced. Particle velocity profiles show good agreement with experimental results and self-similarity of the radial distribution of axial particle velocities is reported. Gas flow patterns are also studied and the effect a vortex ring fixed at the bottom of the vessel is investigated. The simulation is validated through comparisons with results reported in the literature. (C) 2004 Elsevier Ltd. All rights reserved.