The clustering characteristic of purely hydrodynamically interacting particles suspended in pressure-driven flow in a circular cylinder is studied using direct numerical simulation based on the solution of the lattice-Boltzmann equation. We find a universal scaling relation for the cluster size distribution in the subcritical regime for all of the cases considered in this study. This scaling relation is independent of particle shape and concentration.