This paper presents a numerical study of flow through static random assemblies of monodisperse, spherical particles. A lattice Boltzmann approach based on a two relaxation time collision operator is used to obtain reliable predictions of the particle drag by direct numerical simulation. From these predictions a closure law F(Re-p, phi) of the drag force relationship to the bed density phi and the particle Reynolds number Re phi is derived. The present study includes densities phi ranging from 0.01 to 0.35 with Re phi ranging up to 300, that is compiled into a single drag correlation valid for the whole range. The correlation has a more compact expression compared to others previously reported in literature. At low particle densities, the new correlation is close to the widely-used Wen & Yu - correlation. Recently, there has been reported a discrepancy between results obtained using different numerical methods, namely the comprehensive lattice Boltzmann study of Beetstra et al. (2007) and the predictions based on an immersed boundary - pseudo-spectral Navier-Stokes approach (Tenneti et al., 2011). The present study excludes significant finite resolution effects, which have been suspected to cause the reported deviations, but does not coincide exactly with either of the previous studies. This indicates the need for yet more accurate simulation methods in the future. (C) 2014 Elsevier Ltd. All rights reserved.