We present the results of direct numerical simulations (DNS) of turbulent flows seeded with millions of passive inertial particles. The maximum Reynolds number is Re_λ∼ 200. We consider particles much heavier than the carrier flow in the limit when the Stokes drag force dominates their dynamical evolution. We discuss both the transient and the stationary regimes. In the transient regime, we study the growth of inhomogeneities in the particle spatial distribution driven by the preferential concentration out of intense vortex filaments. In the stationary regime, we study the acceleration fluctuations as a function of the Stokes number in the range St ∈ [0.16:3.3]. We also compare our results with those of pure fluid tracers (St = 0) and we find a critical behavior of inertia for small Stokes values. Starting from the pure monodisperse statistics we also characterize polydisperse suspensions with a given mean Stokes, .