The real area of contacts, average film thickness, mean real pressure, and mean hydrodynamic pressure are investigated numerically in this study, especially for the parallel thrust bearing. Model surface is generated numerically with given autocorrelation function and some surface profile parameters. Then the average Reynolds equation contained flow factors and contact factor is applied to predict the effects of surface roughness in mixed lubrication regimes. In this equation, flow factors are defined as correction terms to smooth out high frequency surface roughness and contact factor is introduced to relieve from obtaining the average film thickness. Therefore the computation time to obtain barh h can be reduced.