I am writing the code of FEM applied to lubrication problem in Matlab. I found little introduction of boundary conditions in the circumferential direction of finite width journal bearing in many reference literature. Without consideration of circumferential boundary condition, one pressure distribution can be solved then the negative pressure should be set zero which represent oil film rupture. But when I change the deviation angle 📷,the calculated pressure distribution also changed, and even appeared the positive pressure distribution in whole circumferential domain( 360 degree) which is apparently irrational. I think that the solution domain should be restricted in complete film domain rather than full circumferential domain ( 360 degree). The reason is that Reynolds equation is not applicable for oil film rupture (cavitation) domain. Therefore,
It is necessary to calculate the position of oil film rupture. Alternatively, full Somerfield condition or half Somerfield condition can be adopted to do this easily. But Reynolds condition is difficult to applied to FEM, Furthermore, which is even more so in the case stepped cavity journal bearing or herringbone groove bearing.
So, my question is:Is it necessary to consider the fluid film rupture in the circumferential direction not only consider axial side boundary condition? If it is necessary, how treat the Reynolds boundary?