Particle-laden flows in a vertical channel were simulated using an Eulerian–Eulerian, Anisotropic Gaussian (EE-AG) model. Two sets of cases varying the overall mass loading were done using particle sizes corresponding to either a large or small Stokes number. Primary and turbulent statistics were extracted from these results and compared with counterparts collected from Eulerian–Lagrangian (EL) simulations. The statistics collected from the small Stokes number particle cases correspond well between the two models, with the EE-AG model replicating the transition observed using the EL model from shear-induced turbulence to relaminarization to cluster-induced turbulence as the mass loading increased. The EE-AG model was able to capture the behavior of the EL simulations only at the largest particle concentrations using the large Stokes particles. This is due to the limitations involved with employing a particle-phase Eulerian model (as opposed to a Lagrangian representation) for a spatially intermittent system that has a low particle number concentration.