Combined experimental and multiphase computational fluid dynamics analysis of surface textured journal bearings in mixed regime of lubrication

This paper investigates the effect of surface texturing in a partial pad journal bearing through a series of controlled experiments at operating conditions, promoting mixed or boundary regimes of lubrication. Improvements to load carrying capacity are observed under certain operating conditions. A comprehensive computational finite volume multiphase fluid dynamics analysis, including vapour transport equation and modified finite‐size cavity Rayleigh‐Plesset model, is used to study the effect of indented surface textures in the microscale contact domain and within the individual textures themselves. The results show improved conditions with a textured journal through promotion of micro‐hydrodynamic effect, delaying the effect of lubricant rupture, thus extending the effective load bearing region. A very good agreement is obtained between measurements and predictions.


| INTRODUCTION
Surface texturing is progressively used as an approach to improve load carrying capacity of various tribological conjunctions, as well as mitigating frictional losses. 1,2 Dimples or other shallow cavities have been used to act as reservoirs of lubricant particularly in instances that lead to adverse operating conditions (eg, high load, low sliding speed, or poor lubricant availability). Such conditions inhibit the formation of a coherent film of lubricant, for example, during piston reversals. Research into surface texturing has been pursued vigorously because of the high share of piston frictional losses in all forms of crank-slider mechanisms (including in internal combustion engines). Some degree of success has been reported in various studies in palliation of friction and improvements to the output power [2][3][4][5] Furthermore, the underlying microhydrodynamic mechanism, causing pressure perturbations and enhanced lubricant film thickness, has been demonstrated by a number of detailed analytical and numerical analyses. [6][7][8] The same effects have also been noted for other tribological applications, particularly in the case of seals and thrust bearings, where the initial surface texturing work was conceived. [9][10][11] The unfavourable operating conditions are not confined to pistons, seals, or thrust bearings but to all contact conjunctions that are subjected to stop-start, reciprocating, or inlet reversing conditions. These conditions result in mixed or boundary regimes of lubrication because of a diminishing lubricant entrainment into the contact. Use of surface coatings 5 or other forms of overlay 12 and lubricant rheology, particularly lubricant additives such as friction modifiers, have been commonplace in mitigating excessive boundary friction and thus the ensuing wear. In the case of engine conjunctions, the widely varying nature of transient regime of lubrication in different load bearing engine conjunctions means that an overall system solution to lubricant-surface combination cannot be realised. Thus, surface texturing as a localised approach may be regarded as a promising alternative.
Journal bearings are the most commonly used supporting element for rotating machinery, subjected to radial loading. In internal combustion engines, they are used as the main supporting elements for both the crankshaft and the camshaft. They are also used as the load transmitting elements for the connecting rod-crankshaft and the piston assemblies. Altogether, several journal bearings are used in any engine configuration, and they account for 25% of all the frictional losses of an engine, or equivalent to nearly 5% of the fuel consumption. In particular, the elliptic bore big-end bearings are subject to significant variations in the applied generated cyclic loading under a wide range of engine operating speeds. 12 These conditions are further exacerbated owing to the introduction of new technologies such as cylinder deactivation and intermittent stop-start engine operation in crawling traffic. While the growing take-up of these technologies heralds improvements in fuel efficiency and reduced emissions, recent studies have shown deteriorating frictional losses and other adverse tribological conditions, such as cavitation and transient thermal imbalance. 13,14 As the lubricant used is necessarily the same for all the engine conjunctions, there is only limiting opportunities for its rheological optimisation to suit all engine conjunctions. Thus, surface texturing and coatings are viewed as practical alternatives for improved contact conditions. The salient practical implications of multilayered overlays are adequately addressed in some studies. 12,15 This paper focusses on the effect of surface texturing applied to the journal surface of a partial pad bearing. Lu and Khonsari 16 investigated the effect of dimple features in journal bearings, using several different texture dimensions. They showed that there was little accrued benefit under pure hydrodynamic regime of lubrication. The improvements were noted when running conditions pertained to a mixed regime of lubrication. Morris et al 17 also showed the same for a partial pad journal bearing, operating at low journal speeds, as did Galda et al. 6 Morris et al 17 observed that at higher journal speeds, resulting in hydrodynamic regime of lubrication, not only no significant benefit was noted with a textured journal but also increased side leakage, and thus oil loss was evident. This would be an issue of concern in terms of oil consumption, one which was also noted by Henry et al 11 for the case of thrust bearings. Therefore, a fundamental study of textured bearings using detailed numerical analysis is regarded as very timely. For the case of journal bearings, dimple textures have been the texture form of choice. [18][19][20] Therefore, dimple textures are used in the current study through use of indentation process.
When dealing with tribological problems, using Reynolds equation, where the roughness of the moving surface cannot be ignored or it incorporates special surface features, a transient form of Reynolds equation should be used. The differences emerge as a nonsmooth moving surface is shown to produce different effects compared to a nonsmooth fixed surface as shown in various studies. [21][22][23][24][25][26] Owing to this issue, the current study uses the full Navier-Stokes equations for 2-phase flow, including the effect of cavitation, which does not suffer from the usual assumptions made in deriving Reynolds equation. Therefore, inherently, the approach used includes both the spatial and temporal terms, such as the squeeze film effect and fluid inertia, as well as the usual viscous shear. In this approach, the problem is defined without the need for readjusting the coordinates or modifying any of the terms in the original set of equations. This certainly is not the case with Reynolds equation. It is noteworthy that a study of pure squeeze of an impacting sphere on an oily plate performed by Al-Samieh and Rahnejat 27 with the coordinate system attached to the top layer of the lubricant to study the rapid squeeze caving phenomenon showed direct correspondence with the results obtained when the coordinate system was fixed to the stationary surface, while the moving surface fell under an equivalent gravitational effect (ie, squeeze effect). This is in fact known as the equivalence principle as an allusion to a relativistic framework. The results obtained by Al-Samieh and Rahnejat 27 agreed with the experimental work of Safa and Gohar. 28 Nevertheless, the use of Navier-Stokes approach in the current analysis alleviates the aforementioned usual assumptions of Reynolds approach.
In the current study, a multiphase flow through the journal-pad conjunction is expected as the lubricant film ruptures in the vicinity of the contact exit, preceding cavitation. It is also expected that with pressure perturbation occurring over the introduced microbearing areas (textures), the fall in generated pressure may induce localised vapour cavities in the scale of minutiae. To ensure continuity of multiphase flow, it is necessary to use a mass and momentum continuity analysis, using a combined solution of Navier-Stokes and vapour transport equations. This is the approach undertaken in the current analysis, which is in line with that described by Shahmohamadi et al, 14 who provided detailed experimental verification. The work by Shahmohamadi et al 14 is extended to the scale of minutiae and applied to the texture features on the surface of the journal, an approach not hitherto reported in literature.

| Experimental rig
The experimental rig shown in Figure 1 comprises a turned aluminium alloy journal, which is partially submerged in a reservoir of lubricant (sump). The journal is driven by an electric motor and a power pack at low speed (0-250 rpm), resulting in low speed of lubricant entrainment, thus promoting a mixed regime of lubrication. These conditions enable the study of enhanced lubrication when dimple textures are subsequently introduced on the journal surface. A 180°toughened acrylic bearing pad (bushing arc) conforms closely to the journal, forming a partial pad journal bearing arrangement. The pad is cut from an open cylinder. Weights are hung from the bearing pad to load the conjunction. A pointer is attached to the centre of the pad arc so that the attitude angle can be directly measured. The attitude angle is formed between the line of journal and the bearing pad centres and the vertically applied load. This angle is a direct measure of generated friction torque. The angle is read from a finely graduated scale attached to the test rig structure. Therefore, the attitude angle is used to infer the journal eccentricity and thus the generated friction. The properties of the lubricant selected for the experiment are listed in Table 1.
The lubricant used has the viscosity characteristics similar to those of formulated engine oils at start-up conditions ( Table 1). The tests were performed at ambient temperature of 21°C ± 2°C.
The topography of the journal surface was measured using an infinite focus white light interferometer with a nominal vertical resolution of ±1 nm and a lateral sensitivity of 0.175 μm. The topography of the pad is measured using a stylus owing to its translucence. The stylus has a 90°conical angle and a tip radius of 2 μm. The geometrical forms of the journal and the pad were measured using a coordinate measurement machine with a measurement repeatability of ±1 μm. These are shown in Figure 2. The dimensions, out-of-roundness, and surface roughness of the journal and pad are provided in Table 2.

| Surface texturing
Two alternative journals are used. One is untextured, while the other has a textured surface with indented features (Figure 3). A survey of similar previous work drew on the reported experiences of others with regard to the density  distribution, depth, and diameter of indentations. 6,16,17 The contact conditions in the experimental study reside in the mixed regime of lubrication. A significant body of research has shown that under these conditions, low texture density produces the greatest reduction in friction. 18,[29][30][31][32] Therefore, a low texture density is used for the journal surface.
The chosen values for the surface textures are shown in Table 3. The sizes of texture features and particularly their depth were also chosen in such a way that the size of the machining scars and marks on the shaft could be neglected in comparison. These values represent a compromise between the convenient production of surface textures using the readily available techniques such as indentation process and those suggested in open literature. A Vicker indenter was used to produce the surface textures and ensure a good repeatability for the indentation depth. A mesh of indented textures was fabricated onto the journal surface. Owing to the nature of the indention process, there was small variation in the size of dimples. However, this variation is quite small and can be neglected. Table 3 shows the average/typical size of the texture features, which is also used to create virtual textures in the numerical model.
Surface roughness measurement of the textured surface indicated that few protrusions or burrs were created during the surface texturing process. To remove any burrs, the journal was subjected to a running-in test procedure comprising 10 minutes' running period with a 2.2-kg load and at 1.5 rev s −1 . A separate pad was used during the running-in procedure, different to that used in the main tests. The surfaces then were examined using optical microscopy. It was observed that any burrs were completely removed during the running-in period.

| Experimental procedure
A repeatable test procedure was adopted to determine the effect of surface textures upon the tribological   performance of the partial pad journal bearing. The test rig was run at a given applied load for a range of journal speeds (36-180 rpm). The rotational speed of the journal was measured using a digital tachometer. The low speed of revolution results in an insufficient lubricant entrainment into the contact conjunction, thus promoting a mixed regime of lubrication. Under such conditions, the use of a textured journal is expected to promote microhydrodynamic effect, thus enhancing the lubricant film thickness and as a result palliating the effect of boundary generated friction. It should be noted that surface texturing is expected to reduce frictional losses under such circumstances. Therefore, the study reported here is instructive for such operating conditions and should not be regarded as optimal in a generic sense. A steady speed was maintained for a set interval of time prior to each measurement. This was performed to avoid the complications arising from an initial squeeze film effect. After 4 measurements were taken at each chosen journal speed, the load was altered, and the same range of speed tests was undertaken. Owing to short duration of each test (of the order of a minute), the bulk lubricant temperature remains largely unaffected. Preliminary tests showed that the system reaches a stable state within a minute of testing. A short break was made between successive tests to ensure that the temperature of the lubricant sump was maintained at the recorded ambient temperature of 21°C ± 2°C. The bath of lubricant was maintained at the constant level of 20-mm immersion depth of the journal surface so that there would be a sufficient volume of lubricant to be entrained into the contact for an ideal fully flooded inlet.
At the beginning of each test, the bearing pad and the journal were separated and cleaned thoroughly so that any debris produced by a preceding test was removed. After subsequent reassembly, the rig was recalibrated, and a weight was hung from the pad ensuring, the orthogonal loading condition.
The attitude angle was measured using the pointer and the graduated scale. As the attitude angle varied slightly at a given applied load and speed, the average value of measured attitude angle is plotted under each test condition.
Each testing condition was repeated 4 times to ensure the repeatability of the presented outcomes.
Higher attitude angles are limited owing to the physical constraint of the loading mechanism, while the rotational speeds are limited by the splash of excess of lubricant from the oil sump.

| Experimental results
The experimental results for the untextured journal at various applied loads and speeds are shown in Figure 4.
The results show that as the load decreases, the attitude angle increases as the ratio between the generated frictional torque and applied load alters.
The attitude angle can be related to the friction coefficient as μ ∝ tan φ. The results shown in Figure 5 compare the textured and untextured results at different applied loads and journal speeds. At a load of 0.8 kg ( Figure 5A), the boundary interactions are the least for all the tests performed. Under this condition, there is little difference between the coefficient of friction for the 2 cases (textured and untextured journals). As the load is increased in Figure 5B and 5C, the beneficial effect of surface textures becomes evident on the account of a reducing attitude angle at the higher journal speeds. At low speed conditions, the textured and untextured cases have very similar attitude angles, indicating that the textures have little benefit under the prescribed conditions. This is because at the very low sliding speeds, the contact is particularly starved of lubricant. These findings are in line with those of Lu and Khonsari 16 and Galda et al. 6 The results in Figure 5D and 5E show the 2 clearest examples of benefits of surface texturing, with significant reductions in the recorded attitude angles (ie, proportional to the generated frictional torque).
Overall, the experimental results indicate that the introduced surface textures are beneficial under certain operating conditions. When the load is quite low, little difference is noted between the textured and untextured cases. When comparatively moderate loads are applied, the surface textures produce a beneficial effect. At higher loads and medium-to-low speeds, little advantage is accrued through use of a textured journal.
The benefit of surface texturing appears to be closely related to the prevailing regime of lubrication as it would be expected. All the experiments run in the mixed regime of lubrication. However, the extent of direct asperity interactions changes under different conditions. A numerical analysis of the experimental tests is required to further understand the salient underlying phenomena.

| Navier-Stokes equations
The governing equations for conservation of mass and momenta for each phase in the multiphase flow through the contact are 33 where the indices l and v denote the liquid and vapour phases, respectively. D/Dt is the covariant derivative operator, ρ is the lubricant density, p is the pressure, and τ ij is the viscous shear stress tensor. The body force field vector is F ! . This is typically negligible for lubrication analysis, but left in the model for sake of generality. Furthermore, is the velocity vector with u, v, and w being the components of the velocity in the 3 orthogonal Cartesian directions (x, y, and z). The viscous shear stress tensor is where η is the effective lubricant dynamic viscosity and Δ ij is the Kronecker delta.

| Cavitation model
To incorporate the cavitation effect, a model proposed by Singhal et al 34 was used, in which the growth of cavities is confined to a finite-size conjunction. This is particularly important in the small scale of indented textures. The transport equation for the vapour mass fraction is expressed as where Γ, f v , and V * v denote diffusion coefficient, vapour mass fraction, and the vapour velocity vector, respectively. The vapour generation and condensation rates were accounted for by using source terms, R e and R c , respectively. Following generalised Rayleigh-Plesset equations, Singhal et al 34 proposed the following expressions to take into account such phase change rates: for p>p sat : (6) In the relationships above, S represents the surface tension of the lubricant, while V ch stands for the characteristic velocity, which can be considered as the local relative velocity between the liquid and vapour phases. It is assumed that no gasses are dissolved in the lubricant. Then, the saturation pressure, p sat , can be considered to be equal to the bubble pressure in the cavitation zone at a given temperature. Finally, the empirical constants, C e and C c , specified by Singhal et al, 34 are 0.02 and 0.01, respectively.

| Asperity contact model for mixed and/or boundary regimes of lubrication
Direct interaction of solid surfaces occurs with a reducing journal speed. The interaction of asperities on the counterface surfaces result in mixed or boundary regimes of lubrication. Therefore, a portion of applied load is carried by the interacting asperities, while the remainder is supported by hydrodynamic reaction, thus where the hydrodynamic reaction force is obtained as The load supported by asperities is a function of material properties and topographical properties of the counterfaces. It is important to determine the nature of topography of the counterface surfaces through measurement as reported in Section 2.1. Figure 6 shows that the measured topography of both the journal and the bearing pad conforms very well at 2.7σ from the mean roughness height as assumed by Greenwood and Tripp 35 and as for the majority of engineering surfaces. The area of interaction is shown in Figure 6A.
Then, the convoluted asperity height distribution or the probability density function of the 2 surfaces is described as The Greenwood and Tripp 35 statistical functions are defined as where s is the intermediate integration variable. For asperity load carrying capacity of the counterfaces, n = 5/2 and for asperity contact area n = 2. Therefore, the component of the load supported by the asperity tips in the direction opposite to the direction of FIGURE 6 Measured probability distribution of asperity heights on the counterface surfaces applied load can be expressed as 35 where E′ is the composite effective modulus of elasticity of the contacting pairs, ξ is the unit area density of asperities, and β is the mean radius of curvature of tip of asperities. The statistical function F 5/2 (λ) is a function of the ratio of the film thickness to the combined standard deviation of the roughness of the contacting surfaces; ie, Stribeck oil film parameter λ = h/σ. For ease of computations, a polynomial-fit curve is used to approximate this statistical function, as follows 36 : It is shown that the roughness parameter ξβσ is reasonably constant with a value in the range of 0.03 to 0.05 for engineering surfaces. [35][36][37][38] Furthermore, the ratio σ/β is normally in the range 10 −4 to 10 −2 . 27,30 The composite roughness parameter is obtained through according to Gohar and Rahnejat. 38 In the current study, the measured roughness parameter, ξβσ, and average asperity slope, σ/β, provided mean values close to the upper limits of those typical ranges; ie, ξβσ = 0.05 and σ/β = 0.001.

| Conjunctional friction
In the mixed regime of lubrication in journal bearings, the total frictional torque is obtained as where the viscous torque, M vis , is obtained through integration of the generated shear stress in the lubricantsurface interface as The frictional torque due to boundary friction can be described as follows: This model for friction (hence, torque) accounts for the direct asperity contact along with the non-Newtonian shear of pockets of lubricant entrapped between the asperities. The latter is accounted for by considering the limiting Eyring shear stress 39 for the lubricant used, which in this case is τ 0 = 2 MPa. Furthermore, the parameter ς, which is defined as the pressure coefficient for boundary shear strength of asperities on the softer counterface, is measured using an atomic force microscope in lateral force mode. To measure this parameter for the purpose of the current study, a similar procedure to that described by Buenviaje et al 40 and Styles et al 41 was followed.
The cumulative area of asperity tips, A a , is found as 35 in which F 2 (λ) is a statistical function of Stribeck oil film ratio and a surface with Gaussian asperity distribution is given by 36 where λ cr ≈ 3 is selected as the critical film ratio. Below this critical value, some degree of direct solid boundary interactions would be expected.

| BOUNDARY CONDITIONS AND SOLUTION PROCEDURE
Equations 1, 2, and 4 are solved through the use of finite volume method. For this purpose, a model is constructed using the ANSYS DesignModeler, which is subsequently spatially discretised using tetrahedral size-controlled elements. The entire flow domain is meshed using approximately 1 million tetrahedral cells. Each texture feature contains 260 finite volumes to ensure that the flow is accurately modelled. The current model assumes isothermal conditions. 42 A semi-implicit method for pressurelinked equations was used to couple the velocity field with pressure. In addition, a second-order upwind scheme is used for the discretisation of momenta and for minimising the discretisation errors.
With a constant applied load, the bearing eccentricity ratio depends on the pressure equilibrium over the journal surface. To compute eccentricity, the dynamic mesh method is used. 42 An atmospheric pressure of 101 kPa is assumed at the outlet boundary, implying its exposure to the atmosphere, although Shen et al 43 have shown that cavitation pressures depend on machine running conditions and may be as low as 30 kPa for certain conditions. To take into account the cavitation effect in the tribological contacts, cavitation models and implementation algorithms such those proposed by Jakobsson and Floberg 44 and Olsson 45 or the Elrod cavitation algorithm 46 are commonly used. However, in the current study, the cavitation phenomenon is implemented through conjugate solution of the transport equation and the general Navier-Stokes equations. This removes the need for imposing any particular boundary condition(s) to specify either lubricant film rupture or reformation boundaries. This is a significant fundamental improvement upon the traditional cavitation implementation methods commonly used in the study of tribological conjunctions. Similar advantages can be obtained from alternative methods such as complementary-based models described by Woloszynski et al 47 and Giacopini et al, 48 who provided a 1D mass conserving method, based on a complementary algorithm and Reynolds equation. The algorithm is similar to the Elrod-Adams cavitation algorithm.
The bearing pad is stationary while the journal rotates. Both the solid boundaries are impermeable with lubricant slip assumed at these interfaces. The inlet and outlet pressures are set at the ambient pressure. Finally, the lateral boundaries of the bearing are open to atmosphere.
The analysis is advantageous as it allows for a more accurate description of the cavitation phenomenon. In addition, it allows for the variations of lubricant flow properties across the film. However, flow of lubricant through rough surfaces can also be significant. The average flow model developed by Patir and Cheng, 49 based on the original Reynolds equation, addresses this issue through introduction of flow factors. However, such a model has not been developed for a full set of Navier-Stokes equations. Therefore, it should be noted that the current analysis does not incorporate the roughness effect on the hydrodynamic flow of the lubricant. Figures 7 and 8 show the predicted attitude angle for textured and untextured journal surfaces, as well as those measured under the same conditions and for 2 different applied loads. Note that an attitude angle of zero denotes the journal in static equilibrium (stationary) and subjected only to a normal vertical applied load. In all cases, the attitude angle is quite small, indicating a high eccentricity ratio, corresponding to mixed or boundary regimes of lubrication. This is due to the intentional low journal speeds chosen to promote such conditions, where surface texturing is likely to improve lubricant retention and mitigate undue boundary interactions. Good agreement is noted between the predictions and measurements. The numerical results show consistently slightly lower attitude angle for the textured journal, indicating lower contributions due to asperity interactions. The experimental results show the same trend with an increasing applied load.

| Pressure distribution
With the validated numerical predictions, it is instructive to seek the underlying mechanism for the improved performance of the textured journal relative to the untextured case. Figure 9 shows the central pressure profiles of the generated hydrodynamic contact pressure distributions for the textured and untextured cases. These are coincident, as would be expected, with pressures generated in the first part of the arc of contact (0 < θ ≤ θ c , where θ c is the lubricant film rupture point) with cavitation thereafter. There are pressure perturbations at the positions of the indented surface textures as indicated in the figure. Morris et al 8 noted the same for the case of a thin strip with a contacting face of a typical piston compression ring traversing a series of laser textured chevrons on a flat plate. They showed that the textured features retained a volume of lubricant, which caused pressure perturbation due to micro-hydrodynamic effect, enhancing the load carrying capacity of the contact, thus reducing the measured frictional torque. The results here indicate that the same phenomenon is responsible for improved conditions in the case of the textured journal.

| State of flow
The state of flow through the contact and the extent of lubricated domain are good indications of the role of surface texturing. Figure 10 shows the predictions for multiphase flow through the contact, represented by the volume fraction of vapour content. A zero volume fraction represents regions of full liquid lubricant film (single-phase lubricant flow). It can be seen that the area of single-phase flow is extended further when a textured journal is used ( Figure 10A), when compared with the untextured case ( Figure 10B). The full lubricant film rupture boundary is at the onset of a finite vapour volume fraction (onset of vapour cavities). Therefore, the extent of full lubricant film region has increased with the use of texture features, which induce a micro-hydrodynamic effect through pressure perturbations ( Figure 9) and thus enhance the contact load carrying capacity. Cavitation occurs at 2 distinct physical scales. At the larger scale (overall contact), there is the usual cavitation region beyond the lubricant film rupture boundary. In this region, and extending to the outlet of the contact, multiphase flow occurs with varying levels of vapour content. Cavitation also occurs at the local scale of an individual texture feature within the overall cavitation region. The multiphase flow dynamics within these textures is quite complex as can be seen in the various insets of Figure 10A. In all cases, there is some degree of swirl and counterflow, which yield a varying multiphase-layered content. This complex flow pattern gives rise to localised pressure perturbations even within the cavitation region, which is discernible in Figure 9 after the lubricant rupture boundary, prior to θ = π/2 (contact outlet). The layered nature of the flow in texture features indicates pressure and density variations into the depth of film in all the regions of contact (ie, ∂p/∂z ≠ 0), which necessitates the combined solution of Navier-Stokes and vapour transport equations as highlighted in this paper, instead of the usual solutions using Reynolds equation. Figure 11 shows the variations in load shared by asperities at different operational speeds and under various applied loads. As expected, at lower journal speeds, the asperity load share increases and reaches 11% of the total applied load, whereas at higher journal speeds, this contribution is merely 2% of the total load. This effect can be attributed to the build-up of hydrodynamic effect at higher speeds, generating a thicker lubricant film. The effect of the applied load is also shown in Figure 11. At higher loads, the asperities contribution is higher at all speeds particularly at lower ones. The effect of increasing the journal speed is to reduce the asperity load share, through formation a thicker film. In addition, the contribution of the asperities on carrying load seems to be largely insensitive to the application of the textures. The effect of texturing can best be observed at higher loads, where the mixed regime of lubrication is encountered, in line with expectation.
It is shown that surface textures delay the onset of lubricant film rupture, thus extending the effective load carrying region of the conjunction. Furthermore, the textures cause micro-hydrodynamic pressure perturbation, which improve the local as well as overall load carrying capacity of the contact. Remarkably, their effect is reminiscent of the action of rough cartilage in natural endoarticular natural joints. 38 The flow of fluid within each individual texture is quite complex, comprising entrainment, reverse, and counterflows, but overall demonstrates a layered nature, clearly with density variations heralding some form of discretisation. This would obviously require further in-depth analysis, particularly with very shallow features to ascertain when the continuum nature of such microflows may be breached.

| CONCLUSIONS
Effect of texturing of journal surface in a partial pad journal bearing is investigated through a series of experiments at operating conditions, which promote mixed or boundary regimes of lubrication. Improvements to load carrying capacity are observed under certain conditions. A comprehensive computational finite volume multiphase fluid dynamics analysis, including vapour transport equation and modified finite-size cavity Rayleigh-Plesset void fraction model, is used to study the effect of indented surface textures on the microscale of the contact domain as well as the scale of the individual textures themselves. More than 1 million tetrahedral elements are used in the analysis.
The presented analysis is isothermal. In most relevant applications, generated temperatures can lead to more complex thermal cavitating flows, which in turn may need the inclusion of compressible liberated gas behaviour. These issues constitute the future direction of research into the optimisation of texture geometries.
It is also shown that any gain in tribological performance is marginal and confined to certain envisaged operating conditions such as stop-start and under conditions that promote mixed regime of lubrication. However, this marginal gain, if attained, over a number of bearing conjunctions in, for example, a typical engine and over a large number of vehicles, can have significant accrued economic and environmental impacts. More comprehensive analysis would be required, including texture geometry, density, and distribution as well as a broader range of operating conditions prior to practical applications under FIGURE 11 Fraction of applied load supported by the asperities at various journal speeds the broad engine operating conditions. The work of Tala-Ighil et al 18 demonstrates the complex behaviour of surface texturing in the hydrodynamic regime of lubrication. In hydrodynamic conditions, increased texture area can lead to side leakage and oil loss.