An Euler–Euler anisotropic Gaussian approach (EE‐AG) for simulating gas–particle flows, in which particle velocities are assumed to follow a multivariate anisotropic Gaussian distribution, is used to perform mesoscale simulations of homogeneous cluster‐induced turbulence (CIT). A three‐dimensional Gauss–Hermite quadrature formulation is used to calculate the kinetic flux for 10 velocity moments in a finite‐volume framework. The particle‐phase volume‐fraction and momentum equations are coupled with the Eulerian solver for the gas phase. This approach is implemented in an open‐source CFD package, OpenFOAM, and detailed simulation results are compared with previous Euler–Lagrange simulations in a domain size study of CIT. The results demonstrate that the proposed EE‐AG methodology is able to produce comparable results to EL simulations, and this moment‐based methodology can be used to perform accurate mesoscale simulations of dilute gas–particle flows. © 2017 American Institute of Chemical Engineers AIChE J, 63: 2630–2643, 2017