“…The numerical calculation of the Reynolds lubrication equation together with the Jakobsson-Floberg-Olsson (JFO) complementary boundary conditions for cavitation [87,88,115,160] is not straightforwardly accomplished in such domains, often discretized by unstructured meshes; such complementary conditions aim to ensure the whole mass conservation of the lubricant flow, which are in turn accounted for into the solution by means of the p − θ Elrod-Adams cavitation model [77,78] that defines a modified diffusion-convection Reynolds equation. This solution on irregular grids traditionally uses Finite Element Method (FEM) schemes, which impose difficulties on the discretization process and on satisfying the fluid flow conservation on the moving cavitation boundaries, as well as in providing reliable convergence robustness [30,32,35,82,100,128,158,162]. More recently, a novel FEM formulation that restates the lubrication modelling as a linear complementary problem (LCP) was proposed by [27,93], and an efficient algorithm for the fluid pressure calculation was established by [217].…”