Semi-analytical solution of the Reynolds equation considering cavitation

Abstract

The numerical effort of rotordynamic simulations with hydrodynamic bearings is often dominated by the solution of the Reynolds equation. Because of the nonlinear fluid–structure interaction, this equation needs to be solved in every time step. Although computationally efficient analytical approximations and look-up table techniques exist, these approaches often do not provide the necessary modeling depth and accuracy. Thus, in many cases, numerical methods are preferred despite their computational cost.In recent studies, a semi-analytical solution of the Reynolds equation has been developed based on the scaled boundary finite element method (SBFEM) with the objective of improving the numerical efficiency. The developed approach assumed Gümbel conditions, which means that cavitation was handled in a highly simplified manner. In this work, the SBFEM solution is for the first time combined with a more complex, nonlinear cavitation model based on the Elrod algorithm. Because of the semi-analytical technique, this requires a simplification of the switch function that defines the locations of the pressure and cavitation zones. However, as demonstrated in this paper, the resulting bearing forces are still in good agreement with a standard numerical reference solution based on the finite volume method (FVM). Since the SBFEM model uses only a one-dimensional discretization, its high-order formulation is very straightforward, even in case of a non-equidistant grid with varying shape functions that take the smooth and non-smooth regions of the solution into account. As a result, the numerical effort is reduced significantly in comparison to the standard FVM.