Being able to fully model the flow dynamics of grease, including phase separation, will be highly valuable in the design of lubricated machine elements such as rolling element bearings. Complete models will also be a valuable tool in the process of providing tailor-made greases for different applications. An understanding of the grease flow dynamics enables prediction of grease distribution for optimum lubrication and for the migration of wear- and contaminant particles. In this paper the potential of combined analytical modelling, flow visualizations, and numerical modelling in grease flow dynamics is presented. Specifically, the relation between the rheology of the grease and its impact on the flow motion is of interest in combination with validation of the numerical models in simplified geometries. The numerical models then enable simulations in more complex geometries of particular interest for the grease and bearing industry. It is shown that grease flow is heavily influenced by its non-Newtonian properties and the shear rates in the contact, resulting in distinct regions of yielded and un-yielded grease. Further, the numerical models are shown to match well with experiments and analytical models, enabling numerical models on more complicated geometries.