This work reports on an evaluation of the computational fluid dynamics–discrete element method (CFD-DEM) numerical approach to study the behavior of polymer-filled suspensions in a parallel-plate rheometer. For this purpose, an open-source CFD-DEM solver is used to model the behavior of such suspensions considering different particle volume fractions and different types of fluid rheology. We first validate the numerical approach for the single-phase flow of the continuum phase (fluid phase) by comparing the fluid’s azimuthal velocity and shear stress components obtained from the open-source solver against the analytical expressions given in cylindrical coordinates. In addition, we compare the numerical torque given by the numerical procedure with analytical expressions obtained for Newtonian and power law fluids. For both cases, there is a remarkable agreement between the numerical and analytical results. Subsequently, we investigated the effects of the particle volume fraction on the rheology of the suspension. The numerical results agree well with the experimentally measured ones and show a yield stress phenomenon with the increase of the particle volume fraction.