Assessment of Shaft Surface Structures on the Tribological Behavior of Journal Bearings by Physical and Virtual Simulation

Optimizing the surface topography of cast iron crankshafts offers the opportunity to use this material as an alternative to steel in high-performance combustion engines. In the past, this was not possible due to the higher wear on bearing shells and the higher friction losses in relation to forged steel shafts. In order to find an optimized shaft micro topography, the friction and wear behavior of steel and cast iron shafts with different surface treatments were compared to each other, using a combined physical (experimental) and a virtual (computational) simulation approach. The experiments were carried out with a rotary tribometer using a journal bearing test configuration with the possibility to test real-life bearing shells and shaft specimens, manufactured from real-life crankshafts. In the experiments, a polished steel shaft with low bearing wear was effective. The optimization of cast iron crankshafts by a novel surface treatment showed a significant reduction of bearing wear in relation to the classical surface finishing procedures of cast iron shafts. A computational simulation approach, considering the real-life micro topography by using the Navier–Stokes equations for the calculation of micro hydrodynamics, supports the assessment of fluid friction. The virtual simulation shows, in accordance to the experimental results, only a minor influence of the investigated shaft topographies on the fluid friction. Further optimization of shaft surfaces for journal bearing systems seems possible only by the usage of patterned micro topographies.


Introduction
Hydrodynamic journal bearings are widely used machine elements, and due to their advantages, they are now standard bearings in internal combustion engines. The components of hydrodynamic journal bearings include exchangeable bearing shells, a rotating shaft, and an intermediate lubricant. Optimization of each component can reduce mixed or fluid friction losses, subsequently increasing the efficiency of modern combustion engines. Furthermore, an unstable increase of wear, also known as seizure or scuffing, on the one hand and mild wear during standard start/stop operation on the other hand need to be taken into account. A holistic approach to evaluating friction, seizure resistance, and lifetime limitation by mild start-stop wear is important for the design of components especially for highly loaded journal bearings, such as crankshaft bearings in internal combustion engines.
Crankshafts for high-performance combustion engines are usually made from forged steel to guarantee lifetime in terms of fatigue and bearing wear [1]. Surface finishing processes generate a topography that results in good tribological performance. Crankshafts from cast iron characteristically show problematic metal burrs near the graphite regions, generating higher wear on the bearing shells and poorer tribological performance [2]. Apart from these disadvantages, cast iron crankshafts offer potential for lightweight design and cost-efficient production. Strength levels have been increased and come near to those of forged crankshafts [3].
Optimization of cast iron shafts in crankshaft bearing systems can be supported by the use of a tribological simulation technique. The term simulation refers to physical simulation (experiment) as well as virtual simulation (computational calculation) [4,5]. Both techniques try to mimic the real tribological conditions in the application.
Physical simulation is done by laboratory test models (also defined as physical simulation models) with a wide band of measurement possibilities and input parameter settings. Fundamental model test systems [6] compete with application-near journal bearing rigs [7,8]. A compromise in terms of result transfer and efficiency can be achieved by using an application-near tribometer test methodology [9]. Test strategies, which define the test program, are adapted for the respective goal of research [10]. Examples can be found in [11] for research on the seizure resistance of journal bearing materials in engine applications and in [12] for the assessment of lifetime due to mild wear during start/stop operation.
Using a virtual simulation requires a virtual simulation model containing the components of the tribosystems with specific material and geometry parameters coupled with the boundary conditions at operation (speed, load, and temperature) by mathematical relations. Estimation of friction and mild wear behavior are subjects of implemented virtual simulations of journal bearings [13,14]. Macro models, using elasto-hydrodynamic relations, can be expanded with coupled micro models to take the effects of micro hydrodynamics and solid contact into account [15,16].
A micro model, for consideration of micro hydrodynamics, calculates the flow through the rough lubricating gap and is coupled by flow factors with the macro model. First described by Patir & Cheng [17,18], the micro model uses an approach based on the Reynolds equation, calculating a flow that is disturbed or supported by the roughness topography in relation to ideal smooth surfaces. Exactly this relation between the flow of rough and ideal smooth lubricating gaps is known as the flow factor and can be implemented into the macroscopic Reynolds equation. Improvements to the original method of Patir & Cheng were made, for example, by the implementation of real surface topographies, and no virtual generated ones, and the consideration of cavitation effects [19,20]. Using the Reynolds equation is recommended as long as the errors due to the derivation from the Navier-Stokes equations are small. This is the case for the most macroscopic journal bearing models as well as for micro models with substantially smooth roughness profiles and would be the case for the most common bearing and shaft surfaces [21,22]. The trend in the research towards patterned micro topographies calls for updated simulation strategies using the Navier-Stokes equations for solving the micro model. Based on the method of de Kraker [23,24], authors presented, in a previous publication, a novel approach for modeling surface topography effects in hydrodynamic lubrication [25].
The combination of physical and virtual simulation, and thus experimental and computational results, shows high potential for increasing our basic understanding of tribosystems [26,27]. Therefore, this combination can be used for the optimization of critical journal bearing systems, for example, the crankshaft bearings of automotive combustion engines.
To make use of casted crankshafts in high-performance engines possible, different surface treatments are compared in this article with the goal of comparing and understanding wear and friction properties under application-near conditions. This task is carried out by a physical simulation on a rotational tribometer together with a virtual simulation, using a novel approach for modeling surface topography effects in hydrodynamic lubrication [25]. Besides finding an optimization potential for crankshaft surface topographies, this study allows for comparison and validation between physical and virtual simulation approaches.

Materials and Surface Structures
For an application-near estimation of the influence of shaft structures, the components for the simulation were directly taken out from the production line considering different materials with different surface treatments and subsequently with different roughness parameters, given in Table 1. Based on an optical surface measurement (400 × 800µm), according to ISO 25178 [28], we can show different maximum roughness heights Sz [29] for the collection of shaft variants. Furthermore, the reduced peak height Spk and the reduced valley depth Svk can give an impression of what the Abbott-curve and subsequently the distribution of peaks and valleys may look like. Microscopic images, which visualize the characteristic surface topographies, can be found in Section 5.3. For the five shaft variants, listed in Table 1, a physical simulation was carried out with a focus on wear and friction assessment. A virtual simulation for closer research of the hydrodynamic properties of the micro topography was established for Steel 1 and the most potential alternative variant, Cast Iron 3.
The FEDERAL-MOGUL IROX © is a high performance automotive crankshaft bearing shell, containing an aluminium-based lining material with a polymeric coating on top and an initial surface roughness of Sz = 6.8 µm. SEM (scanning electron microscopy) and EDX (Energy dispersive X-ray spectroscopy) analysis shows intermetallic hard phases of Al-Ni, Al-Cu, and Al-Fe-Si-Mn and soft phases of Sn in the lining material (see Figure 1). The top coating contains boron nitride, potentially as solid lubrication, and further fillers for advanced running performance (e.g., iron oxide and silicon carbide) embedded in the polymeric matrix [30].

The Physical Simulation Model
The model for physical simulation of a crankshaft bearing of an internal combustion engine is implemented on a rotational tribometer TE92 of Phoenix Tribology (see Figure 2). The changeable shaft specimen (turning with speed v) is symmetric loaded by two journal bearing shells with a 120°contact angle using pneumatic air bellows applying the normal force F N . The shaft specimens have a diameter of 50 mm and a width of 14.7 mm. The relative bearing clearance is manufactured to 1.6 ‰, which results in a bearing diameter of ≈ 50.08 mm. The system sits in an open oil reservoir with an 80 mL lubricant volume, which can be heated by electrical resistance heating. To guarantee reproducible hydrodynamic conditions, the heat control uses a thermo couple, measuring the oil temperature in the reservoir (accuracy < 1°K). The combination of free rotation of the bearing adapter, supported by a load cell on the main frame, enables a sensitive friction force measurement and subsequently a friction coefficient measurement with a relative error of ≈ 2.3% (details can be found in [9]). A pre-and post-test measurement of bearing shell weight and thickness is used for determining bearing wear. An exemplary procedure of thickness measurement of one bearing shell is shown in Figure 3. Reproducibility of the weight loss measurement was ±0.02 mg. For the thickness measurement, which uses a kind of outside micrometer, reproducibility was ±0.5 µm. The thickness loss measurement uses the maximum values per row averaged for one test (6 maximum values per test), named the averaged maximum wear height, to compare individual tests between each other. The gravimetrical measurement uses an average of the weight loss of both bearing shells per test.

The Physical Simulation Strategy
For a ranking of the wear height and wear volume, as well as for the evolution of mixed and fluid friction, a start/stop test strategy with exactly 20,040 start/stop cycles is used (see Figure 4). A running-in of the system is achieved at a constant speed of 0.26 m/s, followed by a heating phase with a still standing shaft. The start/stop test procedure is carried out at a 120°C oil temperature and a constant normal force to generate a pressure of 2.4 N/mm 2 on the projected area. We performed two parallel tests per shaft variant. For friction assessment, every 60th acceleration ramp is recorded by 1 kHz high speed data acquisition, resulting in 334 recorded Stribeck curves. An energy approach is used to compare the different shaft variants between each other. For the recorded Stribeck curves, frictional energy is separated into mixed friction and fluid friction (see Figure 5). Energy loss is then calculated for each friction regime by integrating the multiplication of coefficient of friction µ with the normal force F N and the relative velocity of the sliding partners v over time t:

The Virtual Simulation Model
The virtual simulation model is split into a macro model, containing the macroscopic geometry (shaft and bearing shells with specimen holders) and a micro model for considering the micro hydrodynamic behavior. Both are coupled via flow factors, so a change of the micro topography can affect the macroscopic bearing behavior. In the following, the most important information on the simulation models are given. Details can be found in [13] for the macro model and in [25] for the micro model. The EHD macro model is implemented in COMSOL Multiphysics by an FEM-Mesh with around 43,000 elements and uses the original geometry of the shaft specimen and the bearing shells mounted in the bearing shell holders. The Reynolds Equation (3) is solved on the refined mesh of the bearing shell surface, using the thin-film-flow module of COMSOL Multiphysics. The gap height is directly calculated from the elastic deformable geometry of the contact partners. Oil temperature T, dynamic viscosity η, and lubricant density ρ are held at constant values according to the experiment (T = 120 • C, η = 0.00548 Pa·s, ρ = 779.9 kg/m 3 ). Cavitation is considered by using the half-Sommerfeld boundary condition.
For the micro model, the CFD package OpenFOAM is used. On a finite-volume mesh of approximately 300,000 elements, which represents an area of the microscopic lubricating gap of 130 × 90 µm, the Navier-Stokes equations are solved. The model contains the real-life shaft surface topography as well as the real-life shell topography, separated by the nominal lubricant gap height.
According to the flow factor method, the model is completed by an inlet and outlet patch and by side faces with a symmetry boundary condition (no shear forces and no flow). A Newtonian viscosity model with a density-free viscosity of ν = 10 mm 2 /s is used.
Real-life micro topographies are obtained via surface scanning by a focus-variation microscope. Via MATLAB, the received 3D-point cloud is smoothed with a 5th order Savitzky-Golay filter and aligned using a linear fit of the surface. Output of this post-routine was an STL surface mesh, which can be used for the SnappyHexMesh utility of OpenFOAM for generation of the micro model mesh.

The Virtual Simulation Strategy
Details on the solvers and the used approximation schemes, as well as the used convergence criteria can be found in the correlating methodical publications ( [13] for the macro model and [25] for the micro model).
Patir & Cheng [17,18] presented an approach using the Reynolds equation for the flow calculation of the microscopic lubricating gap. Correlating the flow of the rough gap with a gap of the same nominal gap height but containing ideal smooth surfaces leads to flow factors. These factors, the pressure flow factor for a pressure-induced flow and the shear flow factor for a flow induced by the relative movement of the contact surfaces, are implemented in the macroscopic Reynolds equation, describing the overall bearing performance.
The calculation of the microscopic flow is a crucial part of this procedure. Patir & Cheng used the Reynolds equation for this purpose. However, the Reynolds equation is limited in its usage due to the derivation from the Navier-Stokes equations. A micro topography containing significant gradients in height, for example, patterned micro topographies, cannot be calculated properly by the Reynolds equation. The direct use of the Navier-Stokes equations is preferred.
Besides the different solution strategy of the micro model, implementation in the Reynolds equation of the macro model varies. The authors use an approach from de Kraker [23,24] for this coupling task that has been already used previously with virtual geometries [25]. The task in this article, to compute real-life surface topographies, is carried out with a similar procedure as published in [25], which is why only selected relations are written in this article.
Equation (3) depicts the macroscopic Reynolds equation used for the virtual simulation, whereby h * stands for the nominal lubricating gap height, η for the dynamic viscosity, p for the pressure, and U 1 for the velocity of the moving surface. The counter surface, in this specific case the bearing surface, is considered with U 2 = 0. The pressure flow factors χ x p /χ y p and the shear flow factor χ x s are a function of the nominal lubricating gap height h * and taken as a constant over the xand y -directions. Because of the major interest of this research in the fluid friction regime, a no pocket-squeeze effect, according to de Kraker [23], occurs, and the right side of Equation (3) can be set to zero. Furthermore, the flow correction factor χ sp is set to 1, due to the relatively mild flow conditions. Based on previous research [25], no major influence of χ sp is estimated.
The pressure flow factors χ x p and χ y p need to be separated into the individual direction of interest. For shear flow factor χ x s , only the direction of movement, in this case the x-direction, is considered. The flow factors define the relation between the flow through the rough lubricant gap Q p,rough /Q s,rough versus the flow through the ideal smooth lubricant gap Q p,smooth /Q s,smooth .
Flow rates of the ideal smooth surfaces in Equations (5) and (6) Q p,smooth and Q s,smooth can be calculated in an analytical way: In these equations, ∆p is the applied pressure gradient for the pressure flow calculation p Inlet − p Outlet , ∆x is the distance between inlet and outlet, and b is the width of the micro model. The flow rate of the rough lubricating gap is the output of the CFD simulation and can be determined by Equation (9). It calculates the flow through a surface A cross , which lies across the flow direction, using the velocity vectors u from the CFD calculation together with normal vectors n of A cross . Separation into pressure and shear-induced flow is done by the boundary conditions of the CFD simulation.
By solving the macroscopic Reynolds equation given in 3 on the bearing area A, it is possible to calculate the nominal gap height distribution h * by taking the force equilibrium between the resulting vector of the fluid film pressure F hyd and the external normal force vector F N into account. Because this article focuses on the hydrodynamic performance of different surface topographies, a solid contact model and the load carrying capacity of asperities are not considered. Assuming that only positive pressure can be load carrying (half-Sommerfeld boundary condition), the force equilibrium can be written as with the hydrodynamic force as the integrated fluid pressure times the normal vector of the surface A: Knowing pressure and height distribution of the lubricating gap, calculation of the friction force F Fluid f riction and thus the coefficient of friction COF Fluid is possible. Friction force acts against the movement direction of the shaft and is the integrated shear force over the bearing surface A.
A comparison between the fluid friction properties of the shaft variants used for the virtual simulation, as well as to the experimental data from the physical simulation is possible in that way.

Results of the Physical Simulation
All shaft variants vary in frictional behavior, especially static friction, and in terms of bearing lifetime under start/stop operation. For a qualitative comparison, accumulated frictional energy and bearing wear after 20,040 start/stop cycles were considered.

Friction Assessment
The continuous recording of friction data in the form of Stribeck curves allowed us to estimate both boundary and mixed friction, as well as the fluid friction regime. In Figure 6, the original Stribeck curves from a test with a Steel 1 shaft and a test with the Cast Iron 3 shaft are plotted, both taken after 20,040 start/stop cycles. In the fluid friction regime, both variants perform similarly. It should be recognized that COF is plotted on a logarithmic scale to visualize the fluid friction regime more clearly. Differences between these two curves can be detected in the form of a lower static friction for the steel shaft. Vibrations of the friction torque signal due to imperfections of alignment and the fast speed ramp have an impact on the COF curves in Figure 6. For a wider view of the impact of the shaft surface topography on the friction behavior, an energetic approach, described in Section 3.2, was used. The accumulated friction energies for mixed and fluid friction (averaged over parallel tests and normalized with the mixed friction losses of Steel 1) are given in Figure 7. Cast Iron 2 and Cast Iron 3 shafts show significant higher friction losses in the mixed friction regime compared to the steel shafts and the Cast Iron 1 shaft. In the fluid friction regime, no significant influence of the shaft surface topography is detectable.

Lifetime Assessment
Via measurement of weight and thickness loss of the bearing shells, a wear ranking of the different shaft variants could be generated. In Figure 8, the values of the averaged maximum wear height (see Section 3.1) together with the weight loss of the bearing shells are plotted for each shaft variant, averaged over the parallel tests and normalized with the wear results of the Steel 1 shaft. The results of the thickness losses demonstrate the good performance of the classical polished steel and cast iron shafts as well as the cast iron shaft with Novel Surface Treatment 2 in comparison to the high wear height values of the ground steel shaft and the Cast Iron 2 shaft. This ranking is also seen for the gravimetric wear results. Only for Cast Iron 1 and 2 is a changing wear relation observed, which can be attributed to slight geometrical differences in the worn volume. Another important observation is that no shaft variant was able to remove the top coating of the bearing material (coating thickness is around 8-12 µm) within 20,040 start/stop cycles, which visualizes the excellent sliding performance of optimized polymeric bearing materials.

Surface Analysis and Tribo Mechanisms
Friction and wear results indicate an influence of the shaft material and the finishing process on the tribological performance of journal bearings. In order to understand the mechanism behind the observed behavior, surface analysis was carried out. Figure 9 depicts a summary of pre-and post-test SEM images as well as hypotheses for the tribo mechanism taking place. All findings regarding topography changes were additionally confirmed by an Alicona InfiniteFocus 3D confocal microscope.
On the Steel 1 shaft, a fine structured surface with polishing grooves in the running direction is visible. After the test run with 20,040 start/stop cycles, no major changes were observed on the shaft surface, except the formation of tribo layers containing O, Ca, and Zn (Steel 1 in Figure 9 and EDX 1 in Table 2). The smooth surface results in low bearing wear and low frictional losses.
The Steel 2 shaft, which shows high bearing wear and moderate mixed friction, offers a surface topography with higher roughness. The initial as well as the surface after 20,040 start/stop cycles show grinding grooves in the running direction. Besides the formation of tribo layers, a slight smoothing of roughness peaks can be observed (Steel 2 in Figure 9 and EDX 2 in Table 2).
All cast iron shaft variants show for the initial state, as well as after testing, exposed graphite and differ only in the shape of the graphite-surrounding metal burrs. For the classical polished cast iron shaft (Cast Iron 1), the metal burrs were removed or pushed into the graphite during operation. A smooth surface with covering Ca tribolayers was generated. This explains the moderate bearing wear and the low mixed friction losses.
In contrast, the cast iron shaft with Novel Surface Treatment 1 (Cast Iron 2) shows sharp, compact edges around the graphite areas with no significant changes, due to the smoothing or flattening of these burrs during operation. Scratching of the uprising burrs generated poor performance for both aspects, the bearing wear and the frictional losses in the mixed friction regime.
Novel Surface Treatment 2 produces a surface topography given in Figure 9 under Cast Iron 3. The metal burrs were transferred into round elevations. The elevated structure contacts the counterpart and seems to be responsible for the high friction losses in the mixed friction regime. However, there is low bearing wear, which can be attributed to the round surface topography of the elevations. During tribological operation, the elevations are smoothed, and covering Ca-based tribolayers are formed.

Results of the Virtual Simulation
For the two shaft variants with the lowest wear (Steel 1 and Cast Iron 3), a virtual simulation focusing on the hydrodynamics of the micro topography was conducted. Both shaft surfaces were paired with the same bearing shell topography to build the lubricating gap model. For the comparison, pressure flow factors in and perpendicular to the running direction were calculated, plotted in Figure 10 as a function of nominal gap height. It should be recognized that first contact between roughness peaks in the micro model occur at a nominal gap height of 0.7 µm for Cast Iron 3 and 0.5 µm for Steel 1. Steel 1 shaft shows an increasing flow factor in the running direction and a decreasing flow factor perpendicular to the running direction. This can be interpreted as indicating that the flow through the rough gap is larger than that through an ideal gap with the same nominal gap height, in the case of the running direction. Flow perpendicular to the running direction is hindered by asperities and the flow appears to be less than that through an ideal smooth gap with the same nominal gap height.
For Cast Iron 3, the flow shows a decrease in both directions, in and perpendicular to the running direction. The graphite dimple in the center of the calculation area seems to generate isotropic micro hydrodynamic behavior. A hypothesis for the flow-increasing behavior in the range of 2 µm would be the primary influence of the flow-increasing area at the graphite dimple (flow increases with h 3 ). For gap heights below 2 µm, the decreasing gap height beside the dimple reduces the overall flow and generates a flow-decreasing behavior of the surface. Nevertheless, the flow factors of both surface variants, Steel 1 and Cast Iron 3, are only of small dimensions, so that no significant impact on the macroscopic hydrodynamic performance could be expected. The time-consuming calculation of shear flow factors was therefore neglected.
Using the elasto-hydrodynamic macro model of the bearing shell holder and the shaft specimen, a consideration of the received flow factors on the macroscopic friction performance could be conducted. Due to the fact that, for this research, the hydrodynamic properties of the micro topography are of interest, only the fluid friction regime is taken into account. Shear flow factors, as well as a temperature-dependent oil viscosity, are neglected for this basic investigation on the behavior of flow factors (see Section 4.2). Cavitation for the macro model is considered to be in a half-Sommerfeld condition. The results of the COF due to fluid friction force over the relative velocity are plotted in Figure 11 (Steel 1 (Virt.) and Cast Iron 3 (Virt.)) using the flow factors from Figure 10, together with the COF results of the experimental investigations ( Figure 11; Steel 1 (Phys.) and Cast Iron 3 (Phys.)). In addition, the difference in the calculated COF between the Steel 1 shaft and the Cast Iron 3 shaft COF Di f f = COF Steel 1 (Virt.) − COF Cast iron 3 (Virt.) is given (Figure 11; Difference COF Steel/Cast Iron). The calculation of COF f luid is performed down to a relative velocity of 0.6 m/s, where both surface variants are still separated by a lubricant gap. Figure 11. Coefficient of fluid friction force, lubrication gap height, and difference in the COF between the two shaft variants over the relative velocity.
There appears to be no significant difference between the experimental and calculated COF results. However, the simulation results match the experimental data quite well, at a speed above 0.6 m/s. COF Di f f shows a slightly lower COF for the steel shaft in the fluid friction regime. The virtual simulations were carried out up to a speed of 16 m/s, to show that both surfaces result in a similar fluid friction force, due to the increasing gap height and the decreasing influence of flow factors.
Details on the overlapping results of physical and virtual simulation can be found in Figure 12. For both surfaces variants, the results match within the spread of the experimental data.
The calculated gap heights (see Figure 13) indicate that the difference between Steel 1 and Cast Iron 3 is marginal. Comparing the two surfaces with an ideal smooth surface yields only minimal differences (a gap height of an elasto-hydrodynamic model with a smooth surface in Figure 13, "Smooth -elastic"). Regarding the neglecting of elastic deformations, as shown in Figure 13, "Smooth -rigid" plays a much greater role in the hydrodynamic performance. Especially in the area of higher relative speeds, this results in an underestimation of the gap height. This additional result should indicate the qualitative influence of surface topographies and flow factors on the results of a hydrodynamic calculation, which is marginal in relation to the consideration of elastic deformations.  Both the coefficient of friction based on fluid friction force and the lubricating gap height indicate only a minor impact of the surface topography. Moreover, with the inclusion of shear flow factors, no major impact of the micro hydrodynamics can be estimated, based on the already poor difference between the two shaft variants.

Conclusions and Outlook
The presented methodology of a physical simulation for the comparison of shaft topographies indicates significant differences in terms of bearing wear and friction losses. A novel approach for the virtual simulation of micro hydrodynamics shows similar results to those of the physical experiments and provides additional input on the understanding of tribological phenomena. While the possibilities of a combination of these two simulation approaches for research on tribosystems are shown, the main purpose of this article consists in the results of the optimization of cast iron shafts.
Using a smooth shaft surface topography, as can be found for the classical polished steel shaft, shows the best performance in terms of friction losses as well as for the mild wear of bearing shells. Surfaces with higher roughness, e.g., for the ground steel shaft, or inhomogeneous surface topographies, e.g., for the cast iron shafts, cannot achieve the tribological performance of a polished steel shaft. Optimizations of cast iron shafts by flattening the metal burrs around the graphite dimples can reduce the wear but leads to higher friction losses in the mixed friction regime. Virtual simulation could not confirm any potential major influence of the micro hydrodynamics. Comparing the effects of the micro hydrodynamics of the polished steel shaft with the novelly treated cast iron shaft reveals differences in the flow behavior, but overall any influence of micro hydrodynamics is negligible.
Following the fact that state-of-the-art finishing procedures generate surface topographies with subordinate micro hydrodynamic effects, only can a potential patterned micro topography of a bearing shaft show a major impact. Research in the field of patterned micro topographies, using physical and virtual simulation, could help to optimize journal bearings in terms of reducing bearing wear and friction losses.
For future research on this topic, it is planned to further improve both the physical and the virtual simulation methodology. The design process of future patterned micro topographies with a prospect for greater friction reduction can justify the implementation of more complex techniques in virtual simulations, such as shear flow factors or more complex cavitation models.