The formulation presented below follows our earlier work on the soft-EHL problems [13,19] and on the mass-conserving cavitation model [20]. In the present paper, the formulation is further developed for the case of non-planar contact. In particular, the related issue of the choice of the domain on which the Reynolds equation is solved is discussed in detail. The corresponding effects are also studied in the numerical example of Section 4.1.