Efficient solution method for the Reynolds equation with Herschel–Bulkley fluids
G. H.G. van der Meer (TU Delft - Mechatronic Systems Design)
Ron A.J. van Ostayen (TU Delft - Mechatronic Systems Design)
More Info
expand_more
Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.
Abstract
Thin film lubrication problems frequently involve the use of lubricants with non-Newtonian characteristics, and a relatively simple viscosity model that can describe several non-Newtonian fluids is the Herschel–Bulkley relation. This relation can model solid-like properties of a lubricant at low shear stress using a yield stress, while at higher shear stress values shear-thinning or thickening can be included. In literature, this viscosity model has been combined with various governing equations to solve the non-Newtonian thin film problem, resulting in models that range from full 3D CFD simulations, to 1D Reynolds equation based methods. However, something that all of these approaches have in common is that they are either computationally expensive, can only be used for 1D geometries, or use non-exact, regularised versions of the Herschel–Bulkley model for reasons of numerical stability. This paper therefore introduces a method for solving a thin film problem with a non-regularised Herschel–Bulkley lubricant using the 2D generalised Reynolds equation, and this approach is shown to be fast without compromising on accuracy. The increased speed will allow the model to be used more efficiently in complex simulations or design optimisation scenarios.
help_outline