To develop a predictive, experimentally verified computational fluid dynamic model for gas-liquid-solid flow, a 3-D transient computer code for the coupled Navier-Stokes equations for each phase was developed. The principal input into the model is the viscosity of the particulate phase, which was determined from a measurement of the random kinetic energy of the 800-micron glass beads and a Brookfield viscometer. The computed time-averaged particle velocities and concentrations agree with PIV measurements of velocities and concentrations, obtained using a combination of gamma-ray and X-ray densitometry, in a slurry bubble column, operated in the bubbly-coalesced fluidization regime with a continuous flow of water. Both the experiment and simulation show a downflow of particles in the center of the column, and an upflow near the walls, and a nearly uniform particle concentration. Normal and shear Reynolds stresses were constructed from the computed instantaneous particle velocities. The PIV measurement and the simulation produced similar nearly flat horizontal profiles of turbulent kinetic energy of particles.