This work presents the development of an OpenFOAM solver aimed at correctly predicting dynamics of concentrated suspensions when subjected to non-uniform shear flows. The newly implemented solver is able to predict the behavior of a heterogeneous mixture whose characteristics depend on the solid particle local concentration. To simulate such behavior, the conservation equation expressing the time variation of the particle volume fraction has been implemented in OpenFOAM; this was achieved by modifying a pre-existing solver, pimpleFoam, which discretizes the Navier–Stokes system of equation through the PIMPLE algorithm. As a first step, the formulation of the momentum equation has been adapted to correctly solve cases with non-Newtonian fluids. Successively, the Krieger’s correlation has been used to model the viscosity variation in the domain to take in account heterogeneous particle distributions. Finally, the iterative cycle for the solution of the migration equation has been included within the time loop. The above-mentioned code has been successfully validated by comparing the numerical results with the measured data provided by experiments reported in literature.