Gas mixing in a tall narrow fluidized bed operated in the slugging fluidization regime is simulated with the aid of computational fluid dynamics. In the first part, a parametric study is conducted to investigate the influence of various parameters on the gas mixing. Among the parameters studied, the specularity coefficient for the partial-slip solid-phase wall boundary condition had the most significant effect on gas mixing. It was found that the solid-phase wall boundary condition needs to be specified with great care when gas mixing is modeled, with free slip, partial slip and no-slip wall boundary conditions giving substantial differences in the extent of gas back mixing. Axial and radial tracer concentration profiles for different operating conditions are generally in good agreement with experimental data from the literature. Detailed analyses of tracer back mixing are carried out in the second part. Two parameters, the tracer backflow fraction and overall gas backflow fraction, in addition to axial profiles of cross-sectional averaged tracer concentrations, are evaluated for different flow conditions. Qualitative trends are consistent with reported experimental findings. (C) 2010 American Institute of Chemical Engineers AIChE J, 56: 2280-2296, 2010