Effects of surface roughness on bearing performances are investigated numerically in this study, especially for the parallel thrust bearing. Although mating surfaces are parallel and separated by thin fluid film, the pressure distribution is formed due to asperities. Model surface is generated numerically with given autocorrelation function and some surface profile parameters. Then the average Reynolds equation is applied to predict the effects of surface roughness between hydrodynamic and mixed lubrication regimes. In this equation, flow factors are defined as correction terms to smooth out high frequency surface roughness. The correlation length is proposed to get the minimum load for the parallel thrust bearing for various sliding conditions.