We present a verification study of three simulation techniques for fluid-particle flows, including an Euler-Lagrange approach (EL) inspired by Jackson's seminal work on fluidized particles, a quadrature-based moment method based on the anisotropic Gaussian closure (AG), and the traditional two-fluid model. We perform simulations of two problems: particles in frozen homogeneous isotropic turbulence (HIT) and cluster-induced turbulence (CIT). For verification, we evaluate various techniques for extracting statistics from EL and study the convergence properties of the three methods under grid refinement. The convergence is found to depend on the simulation method and on the problem, with CIT simulations posing fewer difficulties than HIT. Specifically, EL converges under refinement for both HIT and CIT, but statistics exhibit dependence on the postprocessing parameters. For CIT, AG produces similar results to EL. For HIT, converging both TFM and AG poses challenges. Overall, extracting converged, parameter-independent Eulerian statistics remains a challenge for all methods. V C 2017 American Institute of Chemical Engineers AIChE J, 63: [5396][5397][5398][5399][5400][5401][5402][5403][5404][5405][5406][5407][5408][5409][5410][5411][5412] 2017 In this section, the governing equations of the fluid and particle phases are briefly presented. The behavior of fluid phase Figure 4. St 5 5 particle fields at t 5 1 for TFM, AG, EL (left to right) and increasing resolution (top to bottom) forHIT case. a p shown for EE methods and particle locations for EL method.Statistics are captured using the various density estimation techniques shown in the legend. Fix and Var refer to kernel density estimation using fixed filter width and variable filter width, respectively. In the case of Fix, the number following refers to the filter width, in multiples of particle diameter, h/d p , and in the case of Var, the number refers to the number of particles sampled, N p .Spectra are shown for increasing grid resolution. For EL, Fix 10 was used for density estimation. Darker lines refer to higher resolutions. For TFM and AG, the black dashed lines are from the highest resolution EL simulation for comparison.Particle volume-fraction variance (left) and granular temperature (right). K refers to grid resolution, D, for EE methods and filter width, h, for EL.