Liquid-phase migration in highly concentrated suspensions undergoing constant-force squeeze flow is modeled numerically by taking into account the time and position dependence of the rheological properties due to changes in the volume fraction of solids. This is done by coupling the equation of motion for a non-Newtonian material that behaves approximately as a Bingham plastic with a continuity equation that includes diffusive flux. The developed model was first tested with experimental data and then used to study the effect of various parameters on liquid-phase migration.