Performance evaluation of rough thrust pad bearing under thermo-elastohydrodynamic lubrication using an improved iterative method

. The asperities present on interacting surfaces of a bearing in ﬂ uence the ﬁ lm formation when the oil ﬁ lm becomes thinner and thinner. The aim of this article is to study the effect of stochastic roughness on bearing performance under thermo-piezoviscous and elastic condition using an average ﬂ ow model. To investigate the present operating conditions, progressive mesh densi ﬁ cation method as a fast and simple algorithm has been applied. The results obtained indicate that transverse roughness generates higher pressure compared to other orientational roughness at various ﬁ lm thicknesses. Maximum pressure, mass ﬂ ow rate and load capacity are larger in transverse roughness compared to other orientational roughness for all values of hydrodynamic roughness parameters. A large sensitivity in load capacity for transverse orientation compared to longitudinal at higher ﬁ lm thickness and small ﬁ lm thickness ratios are witnessed. The frictional coef ﬁ cient in longitudinal orientation is large compared to transverse orientation at all values of ﬁ lm thickness ratio and step ratios. Materials with low elastic modulus undergo large deformation, resulting in generation of two sharp pressure peaks. These results may possess good acceptability to practical applications for studying the effect of surface roughness under thermo-elastohydrodynamic lubrication condition.


Introduction
More attention has been given to lubricant's flow under thin film lubrication in recent times. For large operating efficiency, lubricants with low viscosity are used under high load and high velocity. Under this lubrication regime, operating conditions become much more severe when surface roughness or asperities present on interacting surfaces influence the film formation. The effect of these asperities on performance parameters of the bearing cannot be ignored with reduced film thickness.
Compared to counterformal contacts, the film formation is more in conformal contacts that appear in thrust bearings, journal bearings, and mechanical seals. Previous studies suggest that at high load and high sliding velocity, surface deformation and piezoviscous effect appear under hydrodynamic lubrication. The theory of hydrodynamic lubrication is believed to be valid under small load and small sliding velocity. With decrease in film thickness, the deflection of bounding surfaces becomes comparable to film thicknesses that are caused due to high generated pressure. Carl [1] investigated that small surface deformation caused by a nominal pressure of 7 MPa or less cannot be neglected. Hemingway [2] suggested that the surface deformation under small pressure of 2 MPa plays a significant role in film formation. Kawabata et al. [3] analyzed the elastohydrodynamic lubrication condition in thrust pad bearing with circular bumps and different inlet conditions. Many researches [4][5][6][7] gave more importance to thermal expansion compared to elastic deformation in conformal contacts. Cameron [8] studied the effect of generated heat on viscosity of the lubricant under conformal contacts. The combined effect of elastic deformation and thermopiezoviscous effect on lubricant in conformal contact has been given less importance till date compared to counterformal contact.
Lundberg [9] investigated that microscopic scale roughness gets imparted on engineering surfaces during finishing operations like grit blasting, lapping, and grinding. Lubrication is affected by surface roughness in various ways. Bakolas [10] suggested that microscopic roughness disturbs the inlet pressure built up along with coefficient of friction and stress field at subsurface. Two different models were utilized to evaluate the effect of roughness on lubricant film. The first model was based on stochastic method. Tzeng and Sabiel [11] introduced this model by treating the surface roughness as a random variable having characteristics of a probability distribution function. Based on stochastic model, high amount of work was done on one-dimensional roughness pattern [12][13][14][15][16]. An average flow model was derived by Patir & Cheng [17][18][19], where coefficients called flow factors were used to modify the Reynolds equation to study the effect of two-dimensional roughness. This average model was utilized to study the roughness effect on average film thickness and pressure distribution. This model was called as PC model for average flow and is widely accepted till date [20][21][22][23][24][25][26]. These works utilized the flow factor formula suggested by Patir and Cheng, which used isotropic Gaussian surface and neglected the effect of roughness directional pattern, kurtosis, and skewness. Qiu and Raeymaekers [27] studied the effect of stochastic roughness present inside textured bearing surface on performance parameters of a parallel slider bearing. Wang et al. [28] utilized the PC model to study the effect of roughness parameters along with roughness direction, skewness, and kurtosis on performance parameters of a thrust pad bearing lubricated with water.
The second model was based on deterministic approach in which roughness values are directly implemented in film thickness equation. The film thickness and pressure distribution are simulated using this deterministic model to study its effect on performance parameters of the bearing. This deterministic model was implemented in lubrication problem having a well-defined geometry [29][30][31][32][33][34][35][36].
All these studies suggest that the effect of elastic deformation, thermo-piezoviscous effect of lubricant, and surface roughness (stochastic) in conformal contacts has been given less importance till date compared to counterformal contacts under thin film lubrication. Although, these conditions affect the performance parameters of a bearing, no literature reports the combined effect of all these conditions on performance parameters under thin film lubrication.
The present article deals with the thermo-elastohydrodynamic simulation of a one-dimensional Rayleigh step bearing with rough (stochastic) interacting surfaces using an improved method of direct iteration. The effect of directional orientation of surface roughness and hydrodynamic roughness parameter on load capacity, frictional coefficient, total mass flow rate, and maximum temperature has been studied. The results obtained under different orientations have been compared with those obtained under smooth conditions to showcase a detailed view of roughness effect on bearing performance. Figure 1 represents the illustrative sketch of a onedimensional Rayleigh step bearing with stochastic roughness on both the surfaces with same directional pattern (g 1 = g 2 = g) and same standard deviation of roughness heights (s 1 = s 2 )with combined standard deviation s ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi

Mathematical modeling
The contact area consists of a stationary stepped pad of width l and a flat runner moving with velocity U s . The directional advancement of runner is from left to right.

Governing equations
The modified form of one-dimensional average Reynolds equation considering stochastic roughness [37,38] to calculate hydrodynamic pressure is given as follows: where h is the nominal film thickness, h avg is the average film gap between two surfaces, h is the lubricant's viscosity, r is the lubricant's density, p is the generated pressure in the lubricant film, f x is the pressure flow factor, and U s is the moving surface's sliding velocity. The average film gap between two surfaces is given by where normal probability density function f (d) for Gaussian distribution is defined as By introducing the dimensionless parameters in equation (1), the dimensionless form of modified average Reynolds equation is as follows: The discretized form (using FDM) of the dimensionless equation (4) is given by The boundary conditions satisfying equation (5) are P = 0 for X = 0 and X = 1.

Film thickness equation
The film thickness equation considering elastic deformation of the interacting surfaces (assuming them semiinfinite bodies) is given by where d represents the equivalent elastic deformation of the interacting surfaces. The equivalent elastic deformation of the interacting surfaces at any point x due to pressure at point s is governed by the following equation: where E 0 is the equivalent elastic modulus,E 1 is the Young's moduli of bearing pad, and E 2 is the Young's moduli of moving surface. # 1 and # 2 are the Poission's ratios of the respective materials. The dimensionless form of the equations (6) and (7) are The discretized form of equation (9) can be written as where The discretized (dimensionless) form of average film gap is where s Ã ¼ s s is the error function and is defined as The discretized (dimensionless) form of pressure flow factor [17] is ; for g > 1: The constants r and c are function of roughness pattern g (as shown in Fig. 2) and have been defined in Table 1. The pressure flow factor (' x ) is the correction factor for pressure-induced flow (Poiseuille flow) in rough condition. It compares the pressure-induced flow in smooth surface with rough surface. Its dependency on hydrodynamic roughness parameter L ¼ H is shown in Figure 3.

Lubricant's density and viscosity equations
The discretized form of dimensionless density-pressuretemperature relationship with linear thermal correction factor [39] used in the present study is expressed as where T 0 is the lubricant's temperature at inlet, b is the thermal expansivity coefficient of the lubricant, and T i * is the dimensionless temperature in lubricant film at node i.
The discretized form of dimensionless viscosity-temperature-pressure relationship [40] used in the present study is expressed as where z and a are thermal viscosity coefficient and pressure viscosity coefficient of the lubricant, respectively.

Energy equation
Energy equation for one-dimensional Rayleigh step bearing with stochastic roughness is given by Integrating the energy equation (18) across the average film thickness by assuming that (a) h, T, and p do not change across the film thickness, (b) convective heat transfer across film thickness, conductive heat along film, and effect of temperature on r is negligible [41], the one-dimensional energy equation for temperature distribution is where ∂p ∂x is the flow rate per width and C p is the lubricant's specific heat.
Using the dimensionless parameter, The dimensionless form of temperature distribution equation (19) is where The discretized form of dimensionless energy equation (20) is

Mass flow rate
Total mass flow rate per unit width of the lubricant is where,q c , q p , and q are the Couette, Poiseuille, and total mass flow rate, respectively. Using the dimensionless parameters, the dimensionless form total mass flow rate equation (22) is The discretized form of dimensionless mass flow rate equation (24) is

Load capacity
Integrating the pressure distribution over the entire contact area, the bearing's load capacity per unit length is given by The dimensionless form of load capacity is where w Ã ¼ wh 2 2 6h 0 U s l 2 is the dimensionless parameter of the load capacity.
The discretized form of dimensionless load capacity equation (27) is

Frictional coefficient
Integrating the shear stress developed over the entire contact area, the frictional force on moving surface is given by Using the dimensionless parameter, the dimensionless form of friction equation (29) is where F ¼ fh 2 h 0 U s l is the dimensionless parameter. The dimensionless coefficient of friction is The above set of discretized dimensionless equations are solved iteratively using an improved method for direct iteration (progressive mesh densification (PMD) method) [42] for the present analysis. The steps involved in computation are in Figure 4. Pressure P and temperature T * at inlet and outlet are set at zero and one respectively as the boundary conditions for temperature and pressure distribution.

Results and discussion
To study the effect of directional orientation of stochastic roughness on performance parameters of Rayleigh step bearing, the minimum film thickness(h 2 ), film thickness ratios (k), and step ratios (L 1 ) are changed by keeping other parameters fixed. This is similar to changing of load by keeping other parameters fixed. For the present study, equivalent elastic modulus (E ' = 230.76GPa) has been taken, which serves as a consolidation of two steel (E = 210GPa, # = 0.3) bodies. The other essential parameters with fixed values are catalogued in Table 2. Figure 5 compares the result obtained using present algorithm and that of Yagi and Sugimura [43] under elastohydrodynamic lubrication (EHL) condition. The results obtained using the present algorithm is in good agreement with the published one. This shows that present algorithm can serve as a simple and fast algorithm to study EHL condition. Using the present algorithm, effect of directional orientation on performance parameters of the Rayleigh step bearing has been studied under rough thermo-EHL condition in subsequent sections.

Effect of directional orientation of roughness
To study the effect of directional orientation of surface roughness on bearing performance parameters, a modified average Reynolds equation is derived on the basis of probabilistic properties and distribution of roughness height. A flow factor method was developed by Patir and Cheng [17] to study the effect of different directional orientations of surface roughness. Being a powerful method, it has wide acceptance till date in the field of thin film lubrication. This method uses a modified average Reynolds equation (1) implicating flow factor to inspect the response of fluid over the rough surface. The ratio of the minimum film thickness and standard deviation (SD) of roughness height (say hydrodynamic roughness parameter (L)) is used to express the flow factor (Eq. (14)). L characterizes the surface irregularities. A higher value of L represents smaller heights of asperities resulting into almost smoother surfaces. From average film gap and flow factor definition, f x and H avg tend to 1 and H, respectively, when s * tends to 0. In such a case, modified average Reynolds equation reduces to Reynolds equation of smooth surface. Generally, the values at L ≥ 6 is equivalent to values obtained considering smooth surfaces. For the present analysis (thermo-piezoviscous and elastic condition), L = 3 has been taken.
For different values of minimum film thickness (h 2 ), the pressure generation is maximum (0.75 MPa at 1000 nm and 4 MPa at 100 nm) in case of bearing with transverse directional orientation of surface roughness (as shown in Fig. 6). Minimum pressure generation (0.42 MPa at 1000 nm and 3.85 MPa at 100 nm) is witnessed in the bearing with longitudinal orientation of roughness. The trend followed by pressure distribution under different directional orientation is as follows: This trend depends on the pressure flow factor (also called correction factor) (f x ) multiplied with the value of  Fig. 3). Due to this, longitudinal roughness offers less resistance to pressure-induced flow, allowing only a very small amount of side flow. The pressureinduced flow is slightly larger in longitudinal roughness compared to smooth surface because the average gap in valleys are slightly larger than the nominal film thickness (as shown in Fig. 2). Under transverse roughness, the side flow increases and pressure-induced flow decreases, thus offering more resistance to the pressure-induced flow and results into larger pressure generation. For isotropic roughness, the side flow and pressure-induced flow are of same order. Due to this, the flow has to pass around the roughness orientation resulting into higher resistance to the pressure-induced flow compared to smooth surfaces but lesser resistance to the pressure-induced flow compared to transverse roughness.
As h 2 decreases from 1000 to 100 nm, the pressure generation in bearing increases for all orientations of roughness. But the difference between the maximum pressure of transverse and longitudinal roughness decreases from 0.3 to 0.15 MPa with decrease in h 2 from 1000 to 100 nm (as shown in Fig. 6). Reason behind such type of change in maximum pressure is that the pressureinduced flow rate depends upon the film thickness and viscosity. In the value of f x rh 3 h of the pressure-induced flow term, film thickness raises to power equal to 3. Due to reduction in h and increase in viscosity due to pressure, the pressure-induced flow rate decreases, resulting into high pressure generation at smaller film thickness. At this high pressure, the surfaces deform and result into generation of a large no-pressure area at inlet at very small film thickness (100 nm). Asperities under different types of roughness orientations also deform under the influence of this high pressure and become almost close to smooth surface (as shown by film thickness distribution in Fig. 6). Due to asperity deformation, the difference between maximum pressure under different roughness orientations at very small film thickness (100 nm) becomes small compared to that at higher film thickness (1000 nm). This is because the difference between pressure-induced flows under different roughness orientations becomes very small at small film thickness (100 nm) (as shown in Fig. 8).
Temperature reduces the viscosity of the lubricant resulting into an increase in the pressure-induced flow rate. Thus, the pressure generated under thermo-piezoviscous and elastic condition for all roughness orientations is slightly lower compared to piezoviscous and elastic condition.    Figure 7 showcases the temperature distribution at different h 2 for various directional orientations of surface roughness. The temperature distribution with longitudinal and transverse roughness is maximum and minimum, respectively, at all h 2 . This is because the temperature distribution depends on mass flow rate, average film thickness, and difference between pressures of two adjacent nodes (as shown in Eq. (21)). The combined effect of all these parameters causes more viscous shear heating in lubricant, thus, resulting into high temperature under longitudinal orientation compared to transverse. The difference between outlet temperatures of different orientations are negligible for all h 2 . The difference between outlet temperature of transverse and longitudinal roughness increases from 0.0125°C at 1000 nm to 0.25°C at 100 nm. This increase in difference is more at lower value of h 2 (100 nm) because of more pressurized shear heating and large difference in total mass flow rate of transverse and longitudinal roughness (Eq. (21)). Due to existence of nonpressure area at smaller value of film thickness, there exists the parabolic nature of temperature distribution at inlet. Figure 8 shows the flow rates of different roughness orientations at h 2 = 100 nm and are compared to flow rates of smooth surfaces. A constant value of the pressureinduced flow rate and couette flow rate exist in 40% of the contact area from the inlet side for all roughness orientations. This is due to existence of no-pressure zone in this area, which makes the film thickness almost flat (as shown in Fig. 6). The flow rates take the shape of the film thickness, which consists of a flat zone (0 X 0.4), an increasing zone (0.4 X 0.7), a step zone (0.7 X 0.72) , and a slightly diverging and converging zone (0.72 X 1). The difference between pressure-induced flow rates of different orientations is very small. Longitudinal roughness has slightly larger value of pressure-induced flow rate compared to other orientations. Transverse roughness has slightly lower value of pressure-induced flow rate compared to other orientations, thus slightly increasing the generated pressure under this orientation compared to others. The difference between Couette flow rates of different orientations is very small. Transverse roughness has slightly higher value of Couette flow rate compared to other orientations. Longitudinal roughness has slightly lower value of this flow rate compared to others. This is because of the reason that Couette flow rate is proportional to the area of triangle with height h and base U s . This area is less in longitudinal roughness compared to transverse, thus resulting in lower value of Couette flow rate in longitudinal roughness compared to others. Figure 9 shows the effect of L on various performance parameters of the bearing. Due to large resistance to pressure-induced flow in transverse roughness, the maximum pressure and load capacity at all values of L is large compared to other orientations (as shown in Figs. 9a and 9b). The difference between load capacities and maximum pressure of different orientation is large at small value of L. This difference decreases with increase in L. At very large value of L (L ≥ 6), this difference becomes almost negligible and equals to that of smooth surfaces. Total mass flow rate of transverse roughness is large compared to other orientations (as shown in Fig. 9c). This is because of the fact that it is combination of pressure-induced flow rate and Couette flow rate. Pressure-induced flow rate is small and coquette flow rate is large in transverse roughness compared to others. The combined effect of these two flow rates increases the total mass flow rate in transverse roughness. The difference between total mass flow rates of different orientations is large at low values of L. This difference decreases with increase in L and becomes almost negligible at very high value of L. The maximum temperature of film at outlet is generated in longitudinal roughness for all values of hydrodynamic roughness parameter compared to other orientations (as shown in Fig. 9d). This is because of the fact that fluid under longitudinal roughness undergoes large viscous shear heating. The difference between outlet temperature for all orientations is maximum at lower value of L. This difference decreases with increase in L and becomes negligibly small at very high values of L. Figure 10 shows the variation of dimensionless load capacity with respect to film thickness ratios for transverse (Fig. 10a) and longitudinal (Fig. 10b) orientations of surface roughness. It is evident from the graph that the load capacity increases with increase in film thickness. The load capacity for different film thickness increases up to a certain value of k. With further increase in k, the load capacity decreases and becomes almost constant at very high value of k. The difference between load capacity for different h 2 is large at small value of k, which decreases with increase in k. This type of behavior in load capacity is witnessed in both type of roughness orientations. With increase in film thickness, the position of maximum load capacity shifts toward low value of film thickness ratios. In case of transverse orientation, the value of k shifts from 5 to 2.1 when the film thickness increases from 100 to 1000 nm, whereas in case of longitudinal roughness, the value of k shifts from 5.2 to 2.3 with increase in film thickness from 100 to 1000 nm.

Load capacity
The load capacity in transverse roughness (as shown in Fig. 10a) is large compared to that in longitudinal (as shown in Fig. 10b) for all values of h 2 . At smaller value of k, the sensitivity in load capacity is large in case of transverse orientation compared to longitudinal. Due to this, with increase in k, the change in load capacity is small in case of longitudinal roughness compared to transverse for all values of h 2 . This small change tries to slightly flatten the load capacity curve at smaller value of k in case of longitudinal roughness for all h 2 . Figure 11 showcases the effect of step ratios (L 1 ) on load capacity for different value of h 2 under transverse (as shown in Fig. 11a) and longitudinal (as shown in Fig. 11b) orientations. In case of transverse roughness (as shown in Fig. 11a), with increase in film thickness, the value of load capacity increases for all step ratios with some exceptions. At step ratio of around 0.2, the load capacity at film thickness of 800 nm are more compared to other h 2 . However, at very small step ratio (around 0.1), the load capacity at film thickness of 600 nm is large compared to that at other h 2 . With increase in h 2 , the position of maximum load capacity shifts toward large value of step ratios.
In longitudinal roughness (as shown in Fig. 11b), the load capacity increases with increase in film thickness. These values of load capacities are less compared to that under transverse roughness due to increase in the pressure-induced flow (as mentioned in Sect. 3.1), and have some exceptions at lower value of step ratios. At step ratio of around 0.3 and less, the load capacity at film thickness of 600 nm are large compared to that at other h 2 . At very small step ratio (around 0.2 and less), the load capacity at film thickness of 400 nm is large compared to that at 800 and 1000 nm. At larger value of step ratio, the load capacity at film thickness of 800 nm is large compared to that at 1000 nm. The position of maximum load capacity follows the trend of transverse roughness. Figure 12 shows the variation of dimensionless frictional coefficient with respect to k at different minimum film thickness for both transverse (as shown in Fig. 12a) and longitudinal (as shown in Fig. 12b) roughness patterns. It is evident from the graph that at lower value of k, the frictional coefficient is maximum for smaller film thickness. This is due to large decrement in load capacity at smaller film thickness and low film thickness ratios (as shown in Figs. 10a and 10b). With increasing film thickness, the value of frictional coefficient decreases. At higher value of film thickness ratio (k > 8), the frictional coefficient becomes more for thick film (1000 nm) compared to thin film (100 nm). The difference between frictional coefficient for various film thickness becomes small at higher value of k. This trend of frictional coefficient variation with respect to film thickness ratios is same for both transverse and longitudinal roughness. The position of minimum frictional coefficient shifts toward lower value of k with increasing film thickness. In case of transverse orientation, the value of k shifts from 5.5 to 1.8 with increase in h 2 from 100 to 1000 nm for minimum frictional coefficient, whereas the value of k shifts from 6 to 2 in case of longitudinal orientation with increase in film thickness from 100 to 1000 nm for minimum frictional coefficient. The frictional coefficient in longitudinal roughness is slightly larger in comparison to transverse roughness due to decrement in load capacity under longitudinal roughness (as shown in Fig. 10). Figure 13 showcases the effect of step ratios on frictional coefficient. With increasing h 2 , the frictional coefficient decreases up to certain film thickness and then increases for transverse (as shown in Fig. 13a) and longitudinal roughness (as shown in Fig. 13b). This is due to change in load capacity at various step ratios (as shown in Fig. 11). At lower value of L 1 (L 1 < 0.3), the difference between frictional coefficient for different h 2 is large. The minimum value of frictional coefficient exists at lower step ratio (at around 0.3). At very low value of step ratio, the coefficient of friction for 100 and 200 nm is small compared to other film thicknesses for both orientations. This is due to reduction in frictional force at this step ratio. With further increment in step ratio, the difference between frictional coefficient of different h 2 decreases. There exists an exception at 100 and 200 nm where this difference remains almost constant. The value of frictional coefficient becomes almost constant under both orientations for all other h 2 with increasing step ratio (0.3 L 1 0.7). But, at 100 and 200 nm, the value of frictional coefficient constantly increases with increase in step ratio. At very high value of step ratio(L 1 > 0.7), the difference between frictional coefficients of various h 2 further decreases. The value of coefficient of friction at high step ratio increases for all h 2 . This trend of frictional coefficient under both orientations is due to combined effect of load capacity (as shown in Fig. 11) and frictional force on moving surface. The value of frictional coefficient under transverse roughness is small  compared to that under longitudinal roughness. This is due to reduction in load capacity under longitudinal roughness compared to transverse roughness (as shown in Fig. 11). Figure 14 compares the effect of different material parameters on pressure distribution for both transverse (Fig. 14a) and longitudinal (Fig. 14b) orientations. It is evident from the graph that materials with low modulus of elasticity undergo more surface deformation even at small pressure compared to material with high elastic modulus.

Effect of material parameters
With increase in elastic modulus, the area under no pressure zone decreases. At low elastic modulus, two prominent pressure peaks get generated. In between these two pressure peaks, there exists a small area where pressure is very low (almost zero). This low pressure area decreases with increase in elastic modulus. Also, the generated pressure under longitudinal orientation of roughness for all elastic modulus is less compared to transverse orientation of roughness due to reason mentioned in Section 3.2. The no-pressure zone under longitudinal orientation is slightly small compared to that under transverse orientation due to less surface deformation and more pressure-induced flow rate under longitudinal orientation.

Conclusion
This study investigates the effect of stochastic roughness on performance parameters of Rayleigh step bearing under thin film lubrication. The following conclusions have been drawn from the generated results: -The nature and directional orientation of the stochastic roughness play a significant role in pressure generation. Among all orientations, transverse roughness generates more pressure compared to others at various film thicknesses. The difference between generated pressures under different types of orientations goes on decreasing with decrease in film thickness. -The temperature distribution under longitudinal orientation is slightly more compared to other orientations at different film thicknesses. This is due to slightly increased shear heating under this orientation resulting from combined effect of mass flow rate, average film thickness, and difference between pressures of two adjacent nodes. The temperature distribution between different roughness cases is almost identical, with insignificant difference for all values of film thickness. -Maximum pressure, mass flow rate, and load capacity is large in transverse roughness compared to other for all values of hydrodynamic roughness parameters. With increasing hydrodynamic roughness parameter, these values decrease and become almost constant and equal to that of smooth surface. At very high hydrodynamic roughness parameters, these values become almost equal to that of smooth surface. -Maximum temperature at outlet is large in longitudinal roughness compared to others for all hydrodynamic roughness parameters. At very high hydrodynamic roughness parameter, the maximum temperature decreases and becomes almost equal to that under smooth surface.
-Load capacity under transverse roughness at all values of film thickness ratios and step ratios for different film thickness is large in comparison to longitudinal roughness. There exists large sensitivity in load capacity for transverse roughness compared to longitudinal roughness at higher film thickness at different film thickness ratios. At very high film thickness ratios, load capacity for different film thickness becomes almost constant for both (longitudinal and transverse) orientations. With increasing film thickness, the value of maximum load capacity shifts toward low film thickness ratios and high step ratios. -Coefficient of friction in longitudinal roughness is large compared to transverse at all values of film thickness ratios and step ratios for different film thickness. At high film thickness ratios and step ratios, the difference between coefficient of friction for different film thickness becomes less compared to that at low film thickness ratios and step ratios. With increasing film thickness, the position for minimum coefficient of friction shifts toward low value of film thickness ratios and high value of step ratios. -Materials with low elastic modulus undergo more surface deformation even at small pressure compared to material with high elastic modulus, and two prominent pressure peaks get generated. With increase in elastic modulus, the area under no pressure zone decreases. The current investigation suggests that directional orientation of surface roughness influences the performance parameters of the bearing greatly under thin film lubrication. However, a clear understanding of phenomenon related to rough thin film lubrication needs further investigation. The main reason is that the present analysis is based on the assumption of same roughness orientations and same standard deviation of roughness heights on interacting surfaces. However, these surfaces may have same roughness orientations but different standard deviations of roughness height, this condition needs further