  Magnetic Resonance Imaging (MRI) experiments have been performed in
conjunction with Direct Numerical Simulations (DNS) to study neutrally buoyant
particle-laden pipe flows. The flows are characterized by the suspension liquid
Reynolds number ($Re_s$), based on the bulk liquid velocity and suspension
viscosity obtained from Eilers' correlation, the bulk solid volume fraction
($\phi_b$) and the particle-to-pipe diameter ratio ($d/D$). Six different cases
have been studied, each with a unique combination of $Re_s$ and $\phi$, while
$d/D$ is kept constant at 0.058. These cases ensure that the comparison is
performed across different flow regimes, each exhibiting characteristic
behavior. In general, an excellent agreement is found between experiment and
simulation for the average liquid velocity and solid volume fraction profiles.
This study presents, for the first time, accurate and quantitative velocity and
volume fraction profiles of semi-dilute up to dense suspension flows using both
experimental and numerical methods. The discrepancy between the experiments and
simulations can be explained by various reasons, including a difference in
particle size distribution, uncertainty in experimental parameters used as an
input for the DNS, slight variations in particle roughness and frictional
collisions in the simulations. Eventually, three different flow regimes are
identified. For low bulk solid volume fractions a drag increase is observed.
For moderate $\phi_b$ the drag is found to decrease, due to particle
accumulation at the pipe centre. For high volume fractions the drag is found to
decrease further. For solid volume fractions of 0.4 a drag reduction higher
than 25% is found. This drag reduction is linked to the strong viscosity
gradient in the radial direction, where the relatively low viscosity near the
pipe wall acts as a lubrication layer between the pipe wall and the dense core.
