A Comprehensive Numerical Study of a Wedge-Shaped Textured Convergent Oil Film Gap

: The modification of surface geometries to reduce friction is an omnipresent topic of research. In nature, different low-friction surfaces, such as fish skins, exist. To transfer this knowledge to technical applications, for example, to journal or plain bearings, many numerical and experimental studies of textured surfaces have been performed. In this work, the influence of the geometric parameters (texture length l , width b , angle α and start position x start ) of a wedge-shaped texture on three different convergent oil film gaps was analyzed in full-film lubrication and compared with untextured oil film gaps. With the aid of a CFD (computational fluid dynamics) model, a comprehensive variation study was conducted


Introduction
To achieve climate targets, it is necessary to save energy and resources.One opportunity is a reduction in friction.Holmberg et al. [1] reported that 20% of all energy consumption is used to overcome friction.So, the reduction in friction is an omnipresent topic of research.In particular, nature demonstrates perfectly shaped surfaces, such as shark and fish skins [2] and bird wings [3], for drag reduction.Many researchers have tried to transfer these bio-inspired surfaces to technical applications and add additional indentations, called texture, to the surface, which should improve the tribological performance.The limiting factor for the geometry of the texture is the manufacturing process.Hence, simplified geometries, such as dimples, pockets, chevrons and so on, are used.Different surface texture manufacturing techniques exist, and according to [4], these methods can be divided into the following categories:
The advantages, disadvantages and applications of these methods are discussed in [4].The quality of the generated surface is strongly influenced by the variable parameters of the manufacturing technique.In [5], the influence of the surface-texturing laser's pulse duration on the surface topography is depicted.Too long a pulse duration leads to microcracks, high-heat-affected zones and quite irregular surface profiles [5].Another efficient manufacturing process is milling, and this process can be further optimized for texture manufacturing, for instance, by implementing spindle speed modulation, which leads to low costs for texture generation [6].
In [7,8], different textured surfaces for a parallel sliding contact were studied, and a significant drag reduction could be detected.Nevertheless, many technical applications are not parallel.Hence, Bei et al. and Rosenkranz et al. [9,10] examined the tribological performance of a pocket-, dot-and line-cross-like-shaped, textured convergent oil film gap and showed that the geometry of the oil film gap had a significant influence on the load-carrying capacity.Sharma et al. [11] performed a numerical study on a partially textured journal bearing and reported that texturing can enhance the performance for low and average journal eccentricities.The application of this knowledge to a convergent oil film gap demonstrated that partial surface texturing significantly enhances the performance of parallel sliding contacts, which is described in [7, 8,12] and, for a convergent oil film gap with a low pitch angle α 0 , in [10].
Different mechanisms underlying the operating principles of textured surfaces are discussed in the literature.Grützmacher et al. [13] published five different hypotheses about the performance enhancement of textured surfaces:

•
Increase the hydrodynamic pressure; • Supply the surfaces with additional lubricant; • Store lubricant; • Trap wear particles;

•
Reduce the real contact area.
Besides these hypotheses, Hamilton et al. [14] reported that additional hydrodynamic pressure is a consequence of the existence of cavitation.
However, Babu et al. and Hu et al. [15,16] performed a detailed numerical and experimental analysis of different texture shapes and found that the geometrical parameters of every texture, such as the area density and texture height, have a strong influence on the enhancement of a textured surface.Besides the geometry itself, the position of the texture affects the tribological performance.According to Rosenkranz et al. [10], textures that are located at the inlet lead to a friction reduction when the pitch angle α 0 of a convergent oil film gap is low.Furthermore, Yagi et al. [17] found that an opened dimple on the inlet side increased the load-carrying capacity more than an enclosed dimple at the inlet position.Similar results can be observed for a journal bearing.The load-carrying capacity can be increased when the texture is located in a certain region in the convergent region, whereas textures that are placed in the region with the maximum pressure can only impair the performance [18].
As mentioned at the beginning, the surfaces of animals can inspire the development of surface texture.In [2] and Figure 1, a fish skin is depicted.The skin itself consists of different scales.The height profile of such a scale is also visualized in Figure 1.A simplification of the measured height profile, presented in the blue graph, leads to two wedges, presented in the dot-dashed turquoise graph.These simplified profiles served as inspiration for the development of a wedge-shaped texture.
The current paper focuses on a comprehensive study of such a texture on a convergent oil film gap.In [19,20], a wedge-shaped texture, among others, was examined under mixed lubrication.In contrast to van Liem et al. [19], who focused on the size and distribution of textures in mixed lubrication, Maier et al. [20] concentrated on the wear simulation of wedge-shaped surfaces.In both publications, specific geometries were chosen and analyzed.This leads to the question of the influence that the geometric parameters of a wedge-shaped texture have and which geometry will best enhance the performance.As in many technical applications, full-film lubrication is present, and this was not determined for wedge-shape textures by [19,20].In this study, such a lubrication regime will be analyzed.
determined for wedge-shape textures by [19,20].In this study, such a lubrication regime will be analyzed.

Simulation Model
This section focuses on the methodical approach of the simulation model and its validation.

Mathematical Formulation of a Multiphase Flow
Numerical methods have the advantage, compared to experimental methods, of achieving time-and cost-efficient results.More than half of the research methods published in the field of surface texturing are simulation models [21].In 77 % of cases, the Reynolds differential equation is solved [21].Dobrica et al. [22] published a paper discussing the validity of the Reynolds equation in detail.This differential equation is valid for low Reynolds numbers , high dimple aspect ratios  (the ratio between texture length  and texture height ℎ) and small variations in lubricant thickness.The Reynolds number  is the product of density , velocity  and untextured oil film gap height ℎ divided by the kinematic viscosity  [22].In [23], the influence of inertia effects in textured lubricated sliding contacts was discussed noting that these effects become more pronounced at high Reynolds numbers.For this study, a maximum Reynolds number of about  10 is expected, so the influence of inertia effects should be small.This means that the multiphase Stokes flow, which is presented in Equation ( 1), is applicable to this tribological system.
In this equation,  ⃗ represents the velocity,  the pressure field,  ⃗ represents the body force,  the kinematic viscosity,  the second viscosity,  the identity matrix and  the density of the fluid.In the case of multiphase flow, ,  and  are scalar fields.In addition to this equation, the continuity equation (cf.Equation (2)) is needed for timeindependent density .In [24], the mathematical inequality for the viscosity coefficients  and  is presented.Both coefficients were set to the same values using this inequality and an existing simulation model, which is used for validation (cf.Section 2.2).For numerical implementation, the software package Netgen/NGsolve, based on a 3D finite element method (FEM), is used [25].

Simulation Model
This section focuses on the methodical approach of the simulation model and its validation.

Mathematical Formulation of a Multiphase Flow
Numerical methods have the advantage, compared to experimental methods, of achieving time-and cost-efficient results.More than half of the research methods published in the field of surface texturing are simulation models [21].In 77 % of cases, the Reynolds differential equation is solved [21].Dobrica et al. [22] published a paper discussing the validity of the Reynolds equation in detail.This differential equation is valid for low Reynolds numbers Re, high dimple aspect ratios λ (the ratio between texture length l and texture height h) and small variations in lubricant thickness.The Reynolds number Re is the product of density ρ, velocity U and untextured oil film gap height h 0 divided by the kinematic viscosity ν [22].In [23], the influence of inertia effects in textured lubricated sliding contacts was discussed noting that these effects become more pronounced at high Reynolds numbers.For this study, a maximum Reynolds number of about Re ≈ 10 is expected, so the influence of inertia effects should be small.This means that the multiphase Stokes flow, which is presented in Equation (1), is applicable to this tribological system.
In this equation, → u represents the velocity, p the pressure field, → f represents the body force, ν the kinematic viscosity, ζ the second viscosity, I the identity matrix and ρ the density of the fluid.In the case of multiphase flow, ρ, ν and ζ are scalar fields.In addition to this equation, the continuity equation (cf.Equation (2)) is needed for time-independent density ρ.In [24], the mathematical inequality for the viscosity coefficients ν and ζ is presented.Both coefficients were set to the same values using this inequality and an existing simulation model, which is used for validation (cf.Section 2.2).For numerical implementation, the software package Netgen/NGsolve, based on a 3D finite element method (FEM), is used [25].
The multiphase Stokes flow (cf.Equation ( 1)) only needs to be solved if cavitation occurs; otherwise, the incompressible Stokes flow presented in Equation (3) can be solved numerically.Cavitation is determined by finding the minimum pressure and comparing it with the vapor pressure p v .If the minimum pressure is less than p v , cavitation occurs.
In [26], the Merkle [27] cavitation model, among others, for a journal bearing was compared with the measured results in [28] and yielded favorable outcomes.This cavitation algorithm is also used in this study and presented in Equations ( 4)- (9).In this application, the stationary state is of interest, so the time derivative, which is generally a part of the transport equation (Equation (4)), is not required.
The source term, located on the right hand side of the transport equation (cf.Equation ( 4)), describes the vaporization and condensation process.α v represents the vapor volume fraction and α l represents the liquid volume fraction.The sum of both fractions must be 1, and the vapor and liquid volume fractions describe the volume of vapor V v or liquid V l , relative to the total volume (cf.Equation (7)). . .
The vaporization .m − and the condensation term .m + (cf.Equations ( 5) and ( 6)) depend on the condensation constant Cc, the vaporization constant Cv, vapor density ρ v , liquid density ρ l , input velocity U ∞ , reference time t ∞ and vapor pressure p v .
In addition to the previously discussed equations, it is essential to incorporate Equation (8) to compute the density field ρ and Equation (9) to characterize the kinematic viscosity field (ν) when dealing with a multiphase flow.

Comparison of Different Simulation Models
To compare the self-programmed simulation model implemented in this study with an established software tool, OpenFOAM-8 [29] and its interPhaseChangeFoam solver were used.Since mesh generation for a wedge-shaped texture in OpenFOAM is quite complex, the geometry of a block-shaped texture on a convergent oil film gap, as depicted in Figure 2, was created.The bottom surface of the geometry is moving with a velocity U = 10 m/s, while the top surface and the texture surface are static.A Dirichlet boundary condition of p abs = 1.013 bar was applied to the lateral surfaces.Additional boundary conditions for the liquid volume fraction were required.The lateral surface has a pressure of 1.013 bar, resulting in a liquid volume fraction α l of 1 at these surfaces.The fluid parameters and parameters for the Merkle cavitation algorithm are presented in Section 2.3.A grid study was conducted for both models, and the results in the equilibrium state were analyzed.Figure 3 depicts the pressure distribution for different y-coordinate positions for both simulation models.Despite the differences in mathematical approach (finite volume method (FVM) in OpenFOAM and finite element method (FEM) in NGSolve) and element types (hexahedrons in OpenFOAM and tetrahedrons in NGSolve), similar results can be achieved.The comparison of the lifting and drag forces, calculated according to Section 2.3, showed a difference of about 6%.These small deviations between the different simulation models can be explained by different numerical methods and potentially by the fact that the NGSolve model does not include the inertia effects of the multiphase flow.Nevertheless, the use of this simulation model is valid because this variational study requires a large number of simulations (more than 50,000) to understand the various influencing geometric parameters.Neglecting the inertia effects significantly reduces computational times, as the nonlinear term of the Navier-Stokes equation is disregarded.As shown in the Results section, the improvement in tribological performance is much greater than the difference between both simulation models.
Lubricants 2024, 12, 121 5 of 21 volume method (FVM) in OpenFOAM and finite element method (FEM) in NGSolve) and element types (hexahedrons in OpenFOAM and tetrahedrons in NGSolve), similar results can be achieved.The comparison of the lifting and drag forces, calculated according to Section 2.3, showed a difference of about 6 %.These small deviations between the different simulation models can be explained by different numerical methods and potentially by the fact that the NGSolve model does not include the inertia effects of the multiphase flow.Nevertheless, the use of this simulation model is valid because this variational study requires a large number of simulations (more than 50,000) to understand the various influencing geometric parameters.Neglecting the inertia effects significantly reduces computational times, as the nonlinear term of the Navier-Stokes equation is disregarded.As shown in the Results section, the improvement in tribological performance is much greater than the difference between both simulation models.

Simulation Methodology for Wedge-Shaped Textured Surfaces
In Figure 4, the oil film with a global length l 0 , width b 0 and pitch angle α 0 is visible.It was assumed that the surface was smooth, and the roughness of a surface was not included in this simulation model, unlike [30,31], who used flow factors to describe the real surface topography.Within this study, the performance of the texture geometry parameters themselves should be analyzed in detail and the effects of surface roughness should be neglected.
position was varied in a defined interval  = 0.1, 0.2, … .0.8 mm So, all in all, 9000 simulations for geometry I, 20,000 for geometry II an etry III were conducted.In Figure 5, the boundary conditions of an unscaled model are d sliding velocity of  = 10 m/s is defined on the bottom surface.This for technical applications in hydrodynamic lubrication and can also be disc test rig, as can be observed in [32].The lateral surface of the textu Three different global geometries were analyzed in detail, and the parameters themselves are presented in Table 1.Two criteria were important for the selection of these parameters.On the one hand, the tribological performance should be enhanced, and on the other hand, the computation time has to be quite short, so that it is possible to conduct a huge number of simulations in an adequate time.Therefore, a low α 0 and a small global geometry were selected (cf.[10]).The parameters α 0 and h 0 were kept constant to vary as few parameters as possible and to focus this study on the influencing parameters of the texture parameters themselves.The geometrical parameters of the texture, length l, width b angle α and start position x Start were varied in a defined interval, with the requirement that the texture does not overlap the oil film gap.The texture angle α was varied between 0.2 • and 2.2 • .The step size of the geometry parameters between each simulation was 0.1 mm or 0.2 mm for l and b, depending on the width b 0 and length l 0 of the oil film gap.A step size of 0.2 • for α and 0.1 mm for x Start was selected.For only a few diagrams, some additional simulations, where a smaller step size was chosen, were added.The texture start position was varied in a defined interval x Start = [0.1, 0.2, . . . 0.8] mm for all geometries.So, all in all, 9000 simulations for geometry I, 20,000 for geometry II and 20,000 for geometry III were conducted.
In Figure 5, the boundary conditions of an unscaled model are depicted.A constant sliding velocity of U = 10 m/s is defined on the bottom surface.This velocity is typical for technical applications in hydrodynamic lubrication and can also be tested on a pin-on-disc test rig, as can be observed in [32].The lateral surface of the texture, the top surface of the oil film gap and the texture itself are stationary.On the lateral surface of the oil film gap, the ambient pressure p 0 is applied.Additionally, a liquid volume fraction of α l = 1 is defined on this surface.The boundary condition of zero gradient is applied to all the other surfaces that are not mentioned.Table 2 shows the oil and cavitation parameters used for the simulation model, as presented in [33,34].The geometry was discretized using second order 3D tetrahedron-shaped elements.To achieve sufficient precision, simulations with different maximum mesh sizes were conducted.These analyses were used to identify an adequate mesh size where only small changes, less than 1%, could be detected in the results.Based on these results, a maximum mesh size of 0.025 mm was selected for discretizing the geometries.A direct solver called the unsymmetric multi-frontal package UMFPACK was used to solve the partial differential equations.In postprocessing, the lifting force  ⃗ , which is equivalent to the load that can be absorbed by the oil film gap, can be calculated by integrating the pressure field on the bottom surface (cf.Equation ( 10)).To calculate the drag force  ⃗ , the shear stresses  are integrated over the bottom surface.To validate the model, an analytical solution according to [35] for an untextured convergent oil film gap was used to calculate the lifting force  ⃗ and the drag force  ⃗ .The comparison showed similar results.
⃗ .=  ( , , ) •  ⃗ (10) The geometry was discretized using second order 3D tetrahedron-shaped elements.To achieve sufficient precision, simulations with different maximum mesh sizes were conducted.These analyses were used to identify an adequate mesh size where only small changes, less than 1%, could be detected in the results.Based on these results, a maximum mesh size of 0.025 mm was selected for discretizing the geometries.A direct solver called the unsymmetric multi-frontal package UMFPACK was used to solve the partial differential equations.
In postprocessing, the lifting force → F li f t , which is equivalent to the load that can be absorbed by the oil film gap, can be calculated by integrating the pressure field on the bottom surface (cf.Equation ( 10)).To calculate the drag force → F drag , the shear stresses τ x are integrated over the bottom surface.To validate the model, an analytical solution according to [35] for an untextured convergent oil film gap was used to calculate the lifting force F drag -it is necessary to combine these two physical quantities for comparison with an untextured convergent oil film gap.This leads to the performance enhancement ratio (PER) according to Sharma et al. [11] (cf. Equation ( 11)).The two parameters are combined into one to simplify the visualization.If PER is greater than one, the performance is enhanced; otherwise, it is impaired.However, it should be noted that a reduction in friction-→ F drag and lifting force → F li f t can lead to a performance enhancement ratio PER greater than one.

Results
Overall, numerous simulations were conducted to provide an overview of the results.Scatter plots are quite useful for this purpose, while standard plots are commonly used to investigate the influence of different parameters.This section will also focus on a detailed analysis of the velocity and pressure fields of selected texture geometries.

Analysis of Variation Parameters Depending on PER
Figure 6a-h depict the results of geometry I for different start values x start .Each point in each plot represents one simulation for a specific texture geometry with the parameters l, b, α and x start .The color of each point indicates the performance enhancement ratio (PER).The lower limit of the color bar for Figures 6-8 was set as one, allowing the visualization to focus on geometries that improve tribological performance, by increasing lifting force → F li f t , reducing friction or both.The upper limit for these colored plots was determined by finding the maximum PER.For geometry I, it is 2.65, which is achievable at a texture length l = 1.8 mm, width b = 1.7 mm, angle α = 1.2 • and at a start position x start = 0.2 mm.Figures 7 and 8 present the same information as Figure 6, but for geometry II and III.The maximum improvement for geometry II can be achieved with a texture length l = 2.8 mm, width b = 2.6 mm, angle α = 0.8 • and a start position x start = 0.2 mm.For geometry III, the maximal PER is 1.67 at a texture length l = 3.7 mm, width b = 3.4 mm, angle α = 0.6 • and a start position x start = 0.3 mm.All geometries were simulated with the same start positions, and a selection of these results is visible in Figures 6-8.The simulation effort for geometry II and III is higher, compared to geometry I.The results show that the best performance is achieved with low texture start positions x start , and consequently for geometry II's start positions x start = 0.1 mm and x start = 0.2 mm a smaller step size of 0.1 mm for texture length l and width b was defined (cf. Figure 7a) in contrast to Figure 7b-d   Although the maximum PER for all geometries is different, similar trends can be observed.The best performance for every geometry can be achieved, with textures where the sum of texture length  and start position  is equivalent to the global length  , resulting in an open configuration at the inlet side.An investigation of the texture width  shows that tribological performance can be enhanced with increasing texture width  In Table 3, the performance enhancement ratio (PER) for the best operating textures for each geometry is shown.Additionally, the ratio, describing the improvement in lifting PER li f t and drag force PER drag , is illustrated.It is evident for all geometries that the lifting force increases in comparison to the untextured convergent oil film gap.Furthermore, the drag force of these textured geometries decreases compared to the untextured geometry.Additionally, it is noticeable that as the global length l 0 and width b 0 increase, PER li f t in particular decreases.Although the maximum PER for all geometries is different, similar trends can be observed.The best performance for every geometry can be achieved, with textures where the sum of texture length l and start position x start is equivalent to the global length l 0 , resulting in an open configuration at the inlet side.An investigation of the texture width b shows that tribological performance can be enhanced with increasing texture width b to a local maximum, beyond which the PER decreases.The influence of the texture angle α is also depicted in Figures 6-8.As the global length and width increase, the optimum texture angle α decreases and the range for performance enhancement becomes smaller.
The detailed influence of the different parameters is illustrated in Although the absolute values of the PER are different, the same trends for geometry I, II and III can be observed.As already mentioned, an open texture against the movement direction improves the performance the most.An increasing texture length  leads to a rising PER, as shown in Figures 9a-11a.For small texture lengths , a significant deterioration can be noticed.In Figure 9a, it is visible that for texture lengths  > 1.3  the performance is enhanced.This corresponds to a percentage length / of 65 %.Similar      Although the absolute values of the PER are different, the same trends for geometry I, II and III can be observed.As already mentioned, an open texture against the movement direction improves the performance the most.An increasing texture length l leads to a rising PER, as shown in Figures 9a, 10a and 11a.For small texture lengths l, a significant deterioration can be noticed.In Figure 9a, it is visible that for texture lengths l > 1.3 mm the performance is enhanced.This corresponds to a percentage length l/l 0 of 65%.Similar analogies can be observed for geometry II and III.For l > 2.1 mm (geometry II) and l > 2.6 mm (geometry III), the PER is greater than one, which corresponds to a percentage length l/l 0 of 70% and 65%.
The PER depending on texture width b is presented in Figures 9b, 10b and 11b.For all geometries, the curve progression is similar and the maximal PER can be achieved at a related width b/b 0 of 85%, 86.7% and 85.0% for geometry I, II and III.An analysis of the influence of α showed that with increasing global geometry, the maximum enhancement is located at a lower texture angle α (cf.Figures 9c, 10c and 11c).For geometry I, the best improvement can be achieved with a texture angle α = 1.2 • , for geometry II α = 0.8 • and for geometry III α = 0.6 • .
In Figures 9d, 10d and 11d, the PER as a function of the start position x start is presented.In these figures, the texture length l of each simulation point is different, because it was assumed that the texture is open at the inlet side.So, the texture length l was calculated by subtracting x start from l 0 .An optimum relative start position x start /l 0 can be reached at 10.0%, 6.7% and 7.5% for geometry I, II and III.
In summary, as can be seen from the figures discussed, small variations in texture geometry parameters have a significant effect on tribological performance.
For additional analysis, the pressure and velocity field of the best performing texture geometries were investigated.

Pressure and Velocity Field Analysis
In Figures 12-14, the pressure distribution and the calculated volume flow Q for an untextured and best performing textured oil film gap are presented.The small variations that violate volume continuity can be attributed to rounding and numerical errors.The qualitative characteristic of the pressure field for all textured and untextured geometries is similar.Compared to the untextured geometry, the maximum pressure of the textured geometry is located closer to the outlet.Additionally, the texture increases the pressure, as shown in Figures 12b, 13b and 14b.In addition to the pressure field, the volume flow Q between the textured and untextured geometries differs, especially the volume flow Q perpendicular to the direction of movement is significantly larger, which is visually represented by the width of the arrows.This means that 14.9%, 16.5% and 17.4% of the inlet flow are flowing out there for the textured geometry I, II and III.In contrast, only 6.9%, 9.4% and 11.6% of the inlet flow exit on the faces perpendicular to the movement direction for the untextured geometries I, II and III.

Pressure and Velocity Field Analysis
In Figures 12-14, the pressure distribution and the calculated volume flow  for an untextured and best performing textured oil film gap are presented.The small variations that violate volume continuity can be attributed to rounding and numerical errors.The qualitative characteristic of the pressure field for all textured and untextured geometries is similar.Compared to the untextured geometry, the maximum pressure of the textured geometry is located closer to the outlet.Additionally, the texture increases the pressure, as shown in Figures 12b-14b.In addition to the pressure field, the volume flow  between the textured and untextured geometries differs, especially the volume flow  perpendicular to the direction of movement is significantly larger, which is visually represented by the width of the arrows.This means that 14.9 %, 16.5 % and 17.4 % of the inlet flow are flowing out there for the textured geometry I, II and III.In contrast, only 6.9 %, 9.4 % and 11.6 % of the inlet flow exit on the faces perpendicular to the movement direction for the untextured geometries I, II and III.In Figures 15-17, the velocity field of the untextured geometries I, II, III and the best performing textures of these convergent oil film gaps are visible.The texture depicted in these pictures had the highest performance enhancement ratio of all conducted simulations.Since the real height of an oil film gap in comparison to its horizontal and vertical size is small, the z-coordinates of these figures were scaled by a factor of eight.The size of the flow vectors is constant, and the magnitude of the velocity is represented by the color.It is apparent in all slices B-B that the velocity magnitude is larger on the right side.Since all slices were created parallel to the bottom surface with a displacement of 10 µm, the velocity increases towards the right side where the oil film gap is higher.This can be explained by the simple shear case, where the velocity is linear over the oil film gap height and the flow due to pressure is not present.An increasing height consequently leads to a higher velocity for the same vertical position.This can also be observed in the outer regions of the textured geometry, where there is a sudden change in velocity between the crossover of the textured and untextured region.Compared to the untextured geometries, the magnitude of the velocity of the textured geometries is significantly higher, as seen in Figures 15b-17b, due to the increased oil film gap height caused by the texture.Additionally, the y-coordinates of the textured flow vectors are higher compared to the untextured geometry, causing an increase in volume flow  on the faces normal to the directional movement (cf.Figures 12b-14b).Slice A-A of Figures 15a-17a shows an expected velocity contribution where the fluid on the bottom surface is moving with , while the top surface is stagnant.On the right side of the textured geometry slices (cf.Figures 15b-17b), a backstreaming of the oil with a low magnitude can be observed.As the global geometry In Figures 15-17, the velocity field of the untextured geometries I, II, III and the best performing textures of these convergent oil film gaps are visible.The texture depicted in these pictures had the highest performance enhancement ratio of all conducted simulations.Since the real height of an oil film gap in comparison to its horizontal and vertical size is small, the z-coordinates of these figures were scaled by a factor of eight.The size of the flow vectors is constant, and the magnitude of the velocity is represented by the color.It is apparent in all slices B-B that the velocity magnitude is larger on the right side.Since all slices were created parallel to the bottom surface with a displacement of 10 µm, the velocity increases towards the right side where the oil film gap is higher.This can be explained by the simple shear case, where the velocity is linear over the oil film gap height and the flow due to pressure is not present.An increasing height consequently leads to a higher velocity for the same vertical position.This can also be observed in the outer regions of the textured geometry, where there is a sudden change in velocity between the crossover of the textured and untextured region.Compared to the untextured geometries, the magnitude of the velocity of the textured geometries is significantly higher, as seen in Figures 15b, 16b and 17b, due to the increased oil film gap height caused by the texture.Additionally, the y-coordinates of the textured flow vectors are higher compared to the untextured geometry, causing an increase in volume flow Q on the faces normal to the directional movement (cf.Figures 12b, 13b and 14b).Slice A-A of Figures 15a, 16a and 17a shows an expected velocity contribution where the fluid on the bottom surface is moving with U, while the top surface is stagnant.On the right side of the textured geometry slices (cf.Figures 15b, 16b and 17b), a backstreaming of the oil with a low magnitude can be observed.As the global geometry increases, the length of the backstreaming region decreases.

Discussion
The principal mechanisms of textured surfaces, such as the increase in hydrodynamical pressure and the supply of surfaces with additional lubricant as published by Grützmacher et al. [13], were also observed in this study.The increase in pressure and volume flow  can be seen in Figures 12-14.The volume flow  also increases, when the texture is not open against the movement of direction, probably due to an increasing velocity at the inlet.This phenomenon is illustrated in Figure 18, which shows an untextured convergent oil film gap in the center, surrounded by different textures.It is qualitatively presented that the volume flow at the inlet for all four textures increases compared to the untextured oil film gap.For an open texture, the relative inlet surface area  , is higher compared to closed textures or untextured lubrication gaps.The average velocity magnitude  , on this surface is lower than the average magnitude  , of a smooth oil film cap, resulting in an increased volume flow at the inlet  due to the larger surface area  , .These observations for an opened texture are presented in Figures 13 and 14.The three principal mechanisms of surface textures, trap wear particles, reduce the real contact area and store lubricant, according to Grützmacher et al. [13], and so cannot be described with this simulation model.
The hypothesis that cavitation is the reason for additional pressure in some cases [14], or that cavitation can increase the load-carrying-capacity [36], is disproved for the best performing texture geometries, as the minimum pressure for these textures and many others was not lower than the vapor pressure  .In Figure 18, the pressure distribution for different textures is qualitatively presented and it can be observed that for an at-theinlet opened texture, the minimal pressure is not lower than the ambient pressure.Additionally, it is evident that in some cases, cavitation can negatively impact performance by decreasing the lifting force  ⃗ , when the pressure falls below the ambient pressure.However, it is important to note that textures where cavitation occurs and enhance tribological performance do exist.For the texture, which is analyzed in detail in this study, the highest performing enhancement ratio (PER) can be achieved and cavitation does not occur.
In Figures 12-14, it is depicted that textures increase the volume flow  , , perpendicular to the movement direction.Figure 18 similarly presents this qualitatively,

Discussion
The principal mechanisms of textured surfaces, such as the increase in hydrodynamical pressure and the supply of surfaces with additional lubricant as published by Grützmacher et al. [13], were also observed in this study.The increase in pressure and volume flow Q can be seen in Figures 12-14.The volume flow Q also increases, when the texture is not open against the movement of direction, probably due to an increasing velocity at the inlet.This phenomenon is illustrated in Figure 18, which shows an untextured convergent oil film gap in the center, surrounded by different textures.It is qualitatively presented that the volume flow at the inlet for all four textures increases compared to the untextured oil film gap.For an open texture, the relative inlet surface area A in,rel is higher compared to closed textures or untextured lubrication gaps.The average velocity magnitude v in,rel on this surface is lower than the average magnitude v in,rel of a smooth oil film cap, resulting in an increased volume flow at the inlet Q in due to the larger surface area A in,rel .These observations for an opened texture are presented in Figures 13 and 14.The three principal mechanisms of surface textures, trap wear particles, reduce the real contact area and store lubricant, according to Grützmacher et al. [13], and so cannot be described with this simulation model.
The hypothesis that cavitation is the reason for additional pressure in some cases [14], or that cavitation can increase the load-carrying-capacity [36], is disproved for the best performing texture geometries, as the minimum pressure for these textures and many others was not lower than the vapor pressure p v .In Figure 18, the pressure distribution for different textures is qualitatively presented and it can be observed that for an at-the-inlet opened texture, the minimal pressure is not lower than the ambient pressure.Additionally, it is evident that in some cases, cavitation can negatively impact performance by decreasing the lifting force → F li f t , when the pressure falls below the ambient pressure.However, it is important to note that textures where cavitation occurs and enhance tribological performance do exist.For the texture, which is analyzed in detail in this study, the highest performing enhancement ratio (PER) can be achieved and cavitation does not occur.deviation of 0.1  for the medium class and 0.05  for the fine class is acceptable [38].For angle dimension deviations for nominal lengths up to 10 °, a deviation of 1 ° is permissible [38].Comparing these tolerances with Figures 9-11 shows that the texture angle  and the start position  would change significantly with these tolerances, while the texture length  and width  would not be as strongly influenced.

Conclusions
Within this work, a comprehensive numerical study of a wedge-shaped texture on three different convergent oil film gaps was conducted and the influence of the geometry parameters was analyzed in detail.The following conclusion for the already mentioned global geometries can be summarized:

•
The performance of a wedge-shaped texture is strongly influenced by the geometry parameters.An increasing texture length  improves the tribological behavior and the maximal PER can be achieved with an open texture at the inlet.This means that the texture increases the area of the face where the fluid advects.In Figures 12-14, it is depicted that textures increase the volume flow Q out, f ,b perpendicular to the movement direction.Figure 18 similarly presents this qualitatively, showing that the average magnitude of the velocity on the front and back surface v f ,b,y can be increased, compared to the untextured lubricant gap.This results in an increasing volume flow Q out, f ,b perpendicular to the movement direction.However, it should be noted that this is not a universally valid mechanism, because an infinitely expanded untextured oil film gap is also providing a higher lifting force, in contrast to a finitely expanded one [35].Therefore, the additional inflow of the textured oil film gaps potentially improves performance, which is one operational principle according to Grützmacher et al. [13].
As can be seen in Figures 9a, 10a and 11a and Figure 18, an open texture at the inlet noticeably increases the PER.Similarly, in [17], a parallel sliding surface with multiple dimples exhibits a similar behavior.This open texture may be providing additional lubrication, thereby enhancing the tribological performance.
In Figures 9a, 10a and 11a, it can be observed that for texture lengths l less than l < 1.3 mm (geometry I), l < 2.1 mm (geometry II) and l < 2.6 mm (geometry III), the performance enhancement ratio is less than one.When compared to the pressure field of the untextured geometries, cf.Figures 12a, 13a and 14a, it is evident that the area of maximum pressure is located within these lengths.This suggests that textures located only in the area of maximum pressure negatively impact performance (cf. .This phenomenon is also illustrated in the top left corner of Figure 18, showing the pressure distribution of an oil film gap, where the texture is located in the region of the maximal pressure of a smooth lubrication gap.Due to cavitation, the lifting force → F li f t is reduced, despite the higher maximal pressure compared to a smooth wedge.In summary, it can be concluded that textures located in the region of maximal pressure have a negative impact on tribological performance.Liang et al. [18] also published similar results for partially textured journal bearings. Bei et al. [9] published that a maximal load-carrying capacity for a convergent oil film gap and pocket-shaped textures could be achieved at a dimensionless texture width b/b 0 = 70%, which is not very different from the values found in this study.Here, the optimal related width b/b 0 is between 85% and 87% (cf. Figure 18).In [9], the texture is defined through the parameters length, width and height.In this work, the texture height h, (also referred to as texture depth in the literature) is not constant over the texture length l.This led to the conclusion that for this variation study, a texture angle α was used instead of a texture height h.
According to Bei et al. [9], the load-carrying capacity first decreases, then increases and finally decreases over the texture width and depth.Within this work, it was observed that there is an initial increase followed by a decrease in performance, depending on texture width b and angle α.
In [37], the influence of square-shaped dimple depths is discussed in detail and the best performance could be achieved with a dimple depth of h = 30 µm.In [37], nearly the same density ρ and doubled dynamic viscosity η as in this study were used.To compare this value with the results of this variation study, trigonometric calculations are necessary to calculate the maximal height of the texture.For geometry I, the highest PER could be acquired for a maximal texture height of h = 33 µm, for geometry II h = 31.8µm and for geometry III h = 29.1 µm, as visualized in Figure 18.These values are in the same range as those mentioned by Ullah et al. [37].Another interesting observation regarding the texture height h could be made in this study.The absolute oil film gap height at the inlet, which is the sum of the maximal texture height and the maximal height of the untextured convergent oil film gap, is nearly equivalent for all geometries (51.6 µm for geometry I, 53.2 µm for geometry II and 54.5 µm for geometry III).The conclusion of Ullah et al. [37] that a decreasing texture height h improves performance cannot be observed at all in this study.It can be shown that with increasing texture height, the PER first increases and then decreases.One possible reason why this behavior was detected in this study is that in comparison to Ullah et al. [37], simulations with a very small texture angle α, which corresponds to small heights, were also conducted.
In general, it can be observed that the geometry parameters of the texture strongly influence the performance of a lubricated sliding oil film gap.A comparison with the general tolerances outlined in ISO 2768 reveals that for nominal lengths less than 3 mm, a deviation of ±0.1 mm for the medium class and ±0.05 mm for the fine class is acceptable [38].For angle dimension deviations for nominal lengths up to 10 • , a deviation of ±1 • is permissible [38].Comparing these tolerances with  shows that the texture angle α and the start position x start would change significantly with these tolerances, while the texture length l and width b would not be as strongly influenced.

Conclusions
Within this work, a comprehensive numerical study of a wedge-shaped texture on three different convergent oil film gaps was conducted and the influence of the geometry parameters was analyzed in detail.The following conclusion for the already mentioned global geometries can be summarized: • The performance of a wedge-shaped texture is strongly influenced by the geometry parameters.An increasing texture length l improves the tribological behavior and the maximal PER can be achieved with an open texture at the inlet.This means that the texture increases the area of the face where the fluid advects.

•
Depending on the dimensions of the oil film gap, a relative texture width b/b 0 of 85% to 87% leads to the best performance.• With increasing dimensions of the oil film gap, the optimum texture angle α decreases.
The best performance can be achieved if the oil film gap at the inlet side has an absolute height of about 52 − 55 µm, more or less independent of the geometry of the oil film gap.• A relative start position of the texture x start /l 0 , referring to the length of the oil film gap l 0 , between 6.7% and 10% enhances the performance in the best way.• The texture increases the fluid velocity perpendicular to the movement direction, leading to an increasing volume flow perpendicular to the movement direction.• The texture also induces additional pressure.
• Textures that are only located in the area of maximal pressure deteriorate the tribologi- cal performance.

Figure 5 .
Figure 5. Boundary conditions for velocity  ⃗, pressure  and liquid volume fraction  of simulation model, as well as lifting force  ⃗ and drag force  ⃗ (unscaled).

Figure 5 .
Figure 5. Boundary conditions for velocity → u , pressure p and liquid volume fraction α l of simulation

→FF
li f t and the drag force → F drag .The comparison showed similar results.li f t and drag force → F drag can be calculated for different textured geometries.Since different geometries result in different lifting-→ F li f t and friction forces →

Figure 8 .
Figure 8. Overview of all simulations of geometry III for (a) x start = 0.1 mm, (b) x start = 0.3 mm, (c) x start = 0.5 mm and (d) x start = 0.7 mm.

Figures 9 - 21 Figure 9 .
Figure 9. Performance enhancement ratio (PER) for geometry I as a function of texture (a) length , (b) width , (c) angle  and (d) start position  .

Figure 9 .
Figure 9. Performance enhancement ratio (PER) for geometry I as a function of texture (a) length l, (b) width b, (c) angle α and (d) start position x start .

Figure 10 .
Figure 10.Performance enhancement ratio PER for geometry II as a function of texture (a) length , (b) width , (c) angle  and (d) start position  .

Figure 11 .
Figure 11.Performance enhancement ratio PER for geometry III as a function of texture (a) length , (b) width , (c) angle  and (d) start position  .

Figure 10 . 21 Figure 10 .
Figure 10.Performance enhancement ratio PER for geometry II as a function of texture (a) length l, (b) width b, (c) angle α and (d) start position x start .

Figure 11 .
Figure 11.Performance enhancement ratio PER for geometry III as a function of texture (a) length , (b) width , (c) angle  and (d) start position  .

Figure 11 .
Figure 11.Performance enhancement ratio PER for geometry III as a function of texture (a) length l, (b) width b, (c) angle α and (d) start position x start .

Figure 18 .
Figure 18.Qualitative visualization of influencing parameters in wedge-shaped textures.

Figure 18 .
Figure 18.Qualitative visualization of influencing parameters in wedge-shaped textures.

Table 2 .
Fluid and cavitation model parameters.
Parameters Density oil liquid ρ l 857 kg/m 3 Density oil vapor ρ v 0.13 kg/m 3 Dynamic viscosity liquid η l 22.4 mPas Dynamic viscosity vapor η

Table 3 .
Maximal performance enhancement ratios of best operating geometries.