A lack of thermohydrodynamic (THD) models adapted to discontinuous fluid film geometries can be noted in the lubrication-related literature. The finite differences method, largely used in THD analysis, cannot be easily applied to discontinuous flow domains. This article proposes two finite volumes based models, suitable for discontinuous domains. The first model, based on Reynolds generalized equation, lacks in accuracy near the discontinuity line, but gains in computing speed. The second model, based on Navier-Stokes equations, is more accurate but also significantly more demanding in terms of computing requirements. The two models are applied and compared on the simple geometry of a bi-dimensional Rayleigh step bearing. Different configurations are analyzed, allowing the validation of the Reynolds model for most of the cases encountered in fluid film lubrication.