A magnetorheological fluid is a type of smart fluid that can change its rheological properties when a magnetic field is applied to the fluid. The fluid is a suspension of magnetizable particles in a non-magnetic carrier fluid. When a magnetic field is applied to the fluid, the particles will form structures along the magnetic field lines. These structures resist flow, thereby increasing the viscosity of the fluid. In a non-uniform magnetic field that is generated by a permanent magnet, the particles separate from the fluid and aggregate on the surface of the magnet. This phenomenon is of use in hydrodynamic bearings, where textures are used to increase the load capacity of the bearing. Replacing the fluid in the bearing by a magnetorheological fluid and placing permanent magnets at the desired texture locations, results in particles separating from the fluid, aggregating on the permanent magnets and forming textures. These textures are of a self-healing nature, because any particle that is sheared off, returns into the fluid and is replaced by another particle. Currently, not much is known about the formation of these textures and the influence of their formation on the rheological behavior of the magnetorheological fluid. Therefore, several discrete element models are constructed, that can simulate a magnetorheological fluid in non-uniform field. These models use basic physical laws to determine the dynamics of the particles. A single core model is constructed to simulate the behavior of the magnetorheological fluid in small domains using two different methods for the magnetic interaction forces between the particles. A parallel model is made to simulate the magnetorheological fluid in larger domains. Finally, a model that is used to simulate particulate flows, is modified to research whether two-way coupling is required to determine the steady state behavior of the magnetorheological fluid. The models show that the particles do not separate from the fluid due to the attractive force of the magnet by itself, but rather by a combination of the attractive force and the deformation of the structures when a flow is applied. Furthermore, the rheological behavior of the fluid can still be approximated using the standard viscoplastic models. Finally, the modified particulate flow model showed that two-way coupling is not required to determine the steady state behavior of the magnetorheological fluid.