In this study, we seek to investigate the rheological behavior of phosphate suspension with the aim of performing numerical simulations of concentrated slurry flows using the open source framework OpenFoam. Adequate experimental protocols were presented to identify the main characteristics of these complex fluids and to establish relevant rheological models with a view to simulate numerical flow in a cylindrical pipeline. The solver nonNewtonianIcoFoam is used to compute slurry flows through an axisymmetric pipe. For validation, the numerical results are compared to semi-analytical solutions for pipe flow of non-Newtonian fluids. It is shown that the axial velocity distribution and pressure loss computations using the incompressible non-Newtonian solver agree well with the semi-analytical results.