The relentless pursuit of higher power density and efficiency in advanced aero-engines imposes increasingly severe demands on the transmission systems, particularly the high-speed gears operating within them. These components are pushed towards extreme conditions of rotational speed and load. Among various gear types, the helical gear is often favored for such demanding applications due to its inherent advantages: smooth and quiet operation, high contact ratio, and superior load distribution which reduces stress on individual teeth, making it well-suited for high-speed, high-load scenarios. However, under these extreme operating regimes, the lubrication condition at the tooth contact interface becomes the critical factor determining reliability, efficiency, and service life. The lubricant film separating the meshing teeth of a helical gear is subjected to pressures in the gigapascal range, significant shear, and substantial temperature rises. Traditionally, lubrication analysis for such contacts is based on the assumption of perfectly smooth surfaces—an idealization that rarely holds true in practice. In reality, manufacturing processes impart a certain degree of surface roughness to the gear teeth. When the characteristic dimensions of this surface roughness become comparable to the thickness of the lubricating film, its influence on the lubrication characteristics can no longer be neglected. The presence of roughness asperities disrupts the otherwise smooth pressure and film thickness distributions predicted by classical elastohydrodynamic lubrication (EHL) theory, leading to localized pressure spikes, elevated flash temperatures, and potential areas of severely diminished film thickness. This can significantly increase the risk of surface distress, wear, and ultimately, gear failure. Therefore, a comprehensive understanding of how surface roughness, in conjunction with key operational and material parameters, affects the lubrication performance of high-speed helical gears is paramount for their optimal design and reliable operation.
This work presents a detailed investigation into the thermo-elastohydrodynamic lubrication (TEHL) characteristics of a high-speed helical gear pair used in an aero-engine transmission system. A finite-line-contact TEHL model is employed, which is the most appropriate representation for the contact condition along the line of action in a helical gear mesh. Crucially, this model incorporates the influence of surface roughness on the meshing teeth. The analysis systematically explores the effects of roughness parameters—namely wavelength and amplitude—as well as the classical dimensionless parameters governing EHL: the material parameter (G), the speed parameter (U), and the load parameter (W). To quantify the destabilizing effect of roughness on the lubricant film, the concept of the standard deviation of fluctuation for key film characteristic parameters (pressure, thickness, and temperature rise) is introduced and analyzed. The findings provide critical insights into the interplay between surface topography and lubrication performance under high-speed conditions, offering valuable guidance for the design and manufacturing of more robust and efficient helical gear transmissions.
Geometrical Analysis and Contact Mechanics of Helical Gears
The contact geometry of a helical gear pair is fundamentally different from that of a spur gear. During meshing, contact begins as a point, which then extends into a line as the teeth engage more fully. This line of contact sweeps across the tooth flank, and its length varies during the meshing cycle. This progressive engagement characteristic is a key contributor to the smooth operation of helical gears. For the purpose of lubrication analysis, the complex three-dimensional contact of helical gear teeth can be effectively simplified. At any given instant during meshing, a pair of contacting teeth can be equivalently modeled as two opposed conical rollers in contact. The contact line at that instant corresponds to a common generator of these two cones.
This conical model can be further simplified for the purpose of calculating the geometrical gap between the surfaces. At any point along the contact line, the local radii of curvature perpendicular to the contact line are denoted as $R_1$ and $R_2$ for the two gear teeth, respectively. These are related to the radii in the plane containing the gear axis by the base helix angle $\beta_b$: $R_1 = r_1 / \cos\beta_b$ and $R_2 = r_2 / \cos\beta_b$. The equivalent radius of curvature at that point, used in Hertzian contact theory, is then given by the standard formula for two cylinders in contact:
$$ r_K = \frac{R_1 \cdot R_2}{R_1 + R_2} $$
The undeformed geometrical gap, $h_g$, between the two surfaces can be approximated using a parabolic expression, valid because the contact width is much smaller than the radii of curvature:
$$ h_g = h_0 + \frac{x^2}{2 r_K} $$
Here, $h_0$ is a constant reference gap and $x$ is the coordinate along the direction of rolling/sliding (i.e., along the contact line’s length direction when projected onto the plane of action). This geometrical term forms the foundation upon which the elastic deformations and the lubricant film profile are calculated in the full EHL model. The analysis of the helical gear contact thus reduces to solving a finite-line-contact EHL problem, where the contact length is finite in the $y$-direction (across the face width) and the equivalent curvature $r_K$ may vary along it. For numerical simplicity, a representative section or an average/equivalent curvature is often used for a specific meshing position, such as the pitch point analyzed in this study.
Mathematical Model for Finite-Line-Contact TEHL with Surface Roughness
The core of the analysis is a set of coupled equations that describe the physics of lubrication, elasticity, rheology, and heat transfer. The governing equations are presented here in their dimensional forms before being normalized for numerical solution.
Reynolds Equation
The generalized steady-state Reynolds equation governing lubricant flow in a thin film for a finite line contact, accounting for time-varying effects (squeeze film) in a transient sense, is given by:
$$ \frac{\partial}{\partial x}\left( \frac{\rho h^3}{\eta} \frac{\partial p}{\partial x} \right) + \frac{\partial}{\partial y}\left( \frac{\rho h^3}{\eta} \frac{\partial p}{\partial y} \right) = 12 \left( u_e \frac{\partial (\rho h)}{\partial x} + \frac{\partial (\rho h)}{\partial t} \right) $$
where $p$ is the film pressure, $h$ is the film thickness, $\rho$ is the lubricant density, $\eta$ is the lubricant viscosity, $u_e = (u_1 + u_2)/2$ is the entrainment velocity in the $x$-direction, and $t$ is time. The boundary conditions for pressure are $p=0$ at the inlet, outlet, and sides of the computational domain.
Film Thickness Equation
The film thickness equation incorporates the geometrical gap, the elastic deformation of the contacting surfaces under pressure, and the combined surface roughness of the two helical gear teeth:
$$ h(x,y,t) = h_0 + \frac{x^2}{2 r_K} + \frac{2}{\pi E’} \iint_\Omega \frac{p(x’,y’)}{\sqrt{(x-x’)^2 + (y-y’)^2}} dx’ dy’ – s_{12}(x,y,t) $$
Here, $E’$ is the effective elastic modulus $\left( \frac{1}{E’} = \frac{1-\nu_1^2}{E_1} + \frac{1-\nu_2^2}{E_2} \right)$, and the integral represents the elastic deformation calculated using the half-space Boussinesq approximation. The term $s_{12}(x,y,t)$ is the composite surface roughness. To model a realistic yet tractable roughness profile, a sinusoidal representation with end relief is adopted:
$$
s_{12}(x,y,t) =
\begin{cases}
A_{xy} – A_{xy} \cos\left[ \frac{2\pi}{\lambda_x} (x – u_e t) \right], & \text{for } y < -\frac{l}{2} + l_x \text{ or } y > \frac{l}{2} – l_x \\
A_{xy} – A_{xy} \cos\left( \frac{2\pi x}{\lambda_x} \right) \cos\left[ \frac{2\pi}{\lambda_y} (y + l_x) \right], & \text{for } -\frac{l}{2} + l_x \leq y \leq \frac{l}{2} – l_x
\end{cases}
$$
where $A_{xy}$ is the composite roughness amplitude, $\lambda_x$ and $\lambda_y$ are the roughness wavelengths in the entrainment ($x$) and axial ($y$) directions, respectively, $l$ is the total contact length, and $l_x$ is the length of end relief at both ends of the contact line. The time dependence $t$ accounts for the movement of roughness textures through the contact zone.
Rheological Equations: Viscosity and Density
The dramatic changes in lubricant properties with pressure and temperature are captured by empirical relations. The Roelands equation is used for viscosity:
$$ \eta(p, T) = \eta_0 \exp\left\{ (\ln \eta_0 + 9.67) \left[ \left(1 + 5.1 \times 10^{-9} p \right)^Z \left( \frac{T – 138}{T_0 – 138} \right)^{-S} – 1 \right] \right\} $$
where $\eta_0$ is the viscosity at ambient pressure and temperature $T_0$, and $Z$ and $S$ are parameters derived from the pressure-viscosity coefficient $\alpha$ and temperature-viscosity coefficient $\beta_T$, respectively. The density variation is modeled using the Dowson-Higginson relationship:
$$ \rho(p, T) = \rho_0 \left[ 1 + \frac{0.6 \times 10^{-9} p}{1 + 1.7 \times 10^{-9} p} – 0.00065 (T – T_0) \right] $$
Energy and Heat Transfer Equations
The temperature field within the lubricant film is governed by the energy equation, which balances convection, conduction, viscous dissipation, and compressive heating:
$$ \rho c_p \left( u \frac{\partial T}{\partial x} + v \frac{\partial T}{\partial y} + w \frac{\partial T}{\partial z} \right) = k \frac{\partial^2 T}{\partial z^2} – \frac{T}{\rho} \frac{\partial \rho}{\partial T} \left( u \frac{\partial p}{\partial x} + v \frac{\partial p}{\partial y} \right) + \eta \left[ \left( \frac{\partial u}{\partial z} \right)^2 + \left( \frac{\partial v}{\partial z} \right)^2 \right] $$
The velocities $u$, $v$, and $w$ across the film are obtained from the integrated momentum equations. The heat generated in the film is conducted into the contacting solids (the helical gear teeth). Assuming semi-infinite bodies and a moving heat source, the surface temperatures $T_1$ and $T_2$ are given by:
$$
\begin{aligned}
T_1(x,y) &= T_{1,0} + \frac{1}{\sqrt{\pi \rho_1 c_{s1} k_{s1} u_1}} \int_{-\infty}^{x} \frac{k \frac{\partial T}{\partial z}\big|_{z=0}}{\sqrt{x – s}} ds \\
T_2(x,y) &= T_{2,0} + \frac{1}{\sqrt{\pi \rho_2 c_{s2} k_{s2} |u_2|}} \int_{-\infty}^{x} \frac{k \frac{\partial T}{\partial z}\big|_{z=h}}{\sqrt{x – s}} ds
\end{aligned}
$$
Load Balance Equation
The integrated pressure over the contact area must balance the applied external load per unit length $w$:
$$ \iint_\Omega p(x,y) \, dx \, dy – w = 0 $$
Numerical Solution and Dimensionless Parameters
The system of equations is solved numerically using an iterative procedure involving successive over-relaxation for the Reynolds equation, the fast Fourier transform (FFT) for efficient calculation of elastic deformations, and a finite difference method for the energy equation. Prior to solution, the equations are normalized using characteristic parameters from the Hertzian dry contact solution at a reference point (e.g., one end of the contact line). The key scaling parameters are the Hertzian half-width $b$ and maximum Hertzian pressure $p_H$. The dimensionless coordinates, pressure, film thickness, and time are defined as $X=x/b$, $Y=y/b$, $P=p/p_H$, $H=hR_d/b^2$, and $\tau = u_{e,d} t / b$, where $R_d$ and $u_{e,d}$ are the equivalent radius and entrainment speed at the reference point.
The overall behavior of the EHL contact is often discussed in terms of three dimensionless parameters:
$$
\begin{aligned}
G &= \alpha E’ \quad &\text{(Material Parameter)} \\
U &= \frac{\eta_0 u_{e,d}}{E’ R_d} \quad &\text{(Speed Parameter)} \\
W &= \frac{w}{E’ R_d L} \quad &\text{(Load Parameter)}
\end{aligned}
$$
These parameters will be used to generalize the discussion of results concerning the helical gear lubrication.
Results and Discussion: Influence of Roughness and Operational Parameters
The analysis focuses on a specific high-speed helical gear pair operating at a design-point condition representative of an aero-engine application. The base parameters for the gear geometry, material properties, and lubricant (4106# aviation oil) are summarized in the tables below.
| Helical Gear Geometrical Parameters | |
|---|---|
| Number of Teeth (Driver/Driven) | 40 / 57 |
| Normal Module ($m_n$) | 2.0 mm |
| Normal Pressure Angle ($\alpha_n$) | 20° |
| Helix Angle ($\beta$) | 18° |
| Face Width ($B$) | 30 mm |
| End Relief Length ($l_x$) | 1.5 mm |
| Material Properties (at 120°C) | |
|---|---|
| Material | 16Cr3NiWMoVNbE Steel |
| Thermal Conductivity ($k_s$) | 29.82 W/(m·K) |
| Specific Heat ($c_{s}$) | 726.3 J/(kg·K) |
| Density ($\rho_s$) | 7980 kg/m³ |
| Elastic Modulus ($E$) | 209 GPa |
| Lubricant Properties (4106# at 120°C) | |
|---|---|
| Thermal Conductivity ($k$) | 0.143 W/(m·K) |
| Specific Heat ($c_p$) | 2108 J/(kg·K) |
| Density ($\rho_0$) | 896.9 kg/m³ |
| Dynamic Viscosity ($\eta_0$) | 3.191 mPa·s |
| Pressure-Viscosity Coeff. ($\alpha$) | 13.6 GPa⁻¹ |
| Temperature-Viscosity Coeff. ($\beta_T$) | 0.0169 K⁻¹ |
The design-point operating condition for the helical gear is a driver speed of 24,000 rpm, transmitting 300 kW of power, with an oil supply temperature of 120°C. The reference roughness parameters are an amplitude $A_{xy} = 0.05 \mu m$, wavelengths $\lambda_x = 22.5 \mu m$ and $\lambda_y = 2.0 mm$.
General Effect of Surface Roughness on Helical Gear Lubrication
A comparison between the smooth surface solution and the rough surface solution at the pitch point of the helical gear mesh reveals a profound impact. For the smooth case, the pressure, temperature rise, and film thickness distributions exhibit the classic EHL features: a central flat pressure plateau with steep sides, a corresponding central film thickness constriction, and a central peak in temperature rise due to viscous shear. When surface roughness is introduced, these profiles become highly perturbed. The pressure and temperature rise profiles show pronounced oscillations, following the sinusoidal pattern of the underlying roughness, with local pressure spikes exceeding the smooth-surface maximum pressure by up to 50%. More critically, the film thickness profile also oscillates, creating localized valleys where the film thickness is significantly lower than the central nominal value. These valleys, particularly near the contact exit and at the ends of the contact line (where side leakage further thins the film), represent zones of high risk for metal-to-metal contact and the onset of mixed or boundary lubrication. This confirms that when the roughness amplitude is on the same order as the nominal central film thickness (typically sub-micron in high-speed gears), its influence is dominant and cannot be ignored in the design and analysis of helical gear contacts.
Quantifying Roughness Impact: The Standard Deviation of Fluctuation
To move beyond qualitative observation and quantitatively assess the degree of disturbance caused by roughness, the standard deviation of the fluctuation of key parameters is calculated. For any film parameter $f$ (e.g., pressure $p$, temperature rise $\Delta T$, or film thickness $h$), its fluctuation $\delta f$ at a given point relative to the smooth solution $f_{\text{smooth}}$ is $\delta f = f_{\text{rough}} – f_{\text{smooth}}$. The standard deviation $\sigma_f$ of this fluctuation over the domain quantifies the overall instability or variability introduced by the rough surface to that parameter. A larger $\sigma_f$ indicates a greater destabilizing influence of roughness on the lubricant film for that particular condition.
Effect of Roughness Wavelength
The influence of the spatial distribution of roughness asperities is investigated by varying the wavelengths $\lambda_x$ and $\lambda_y$ independently.
Entrainment Direction Wavelength ($\lambda_x$): Increasing $\lambda_x$ while keeping amplitude constant leads to a gradual smoothing of the oscillations in pressure, temperature, and film thickness profiles. The fluctuations become less frequent and their peak magnitudes decrease. This is intuitive: a longer wavelength means fewer asperities per unit length in the direction of lubricant flow, resulting in a more continuous and less perturbed film. The standard deviations decrease significantly with increasing $\lambda_x$, as shown in the table below for the $x$-direction centerline ($y=0$).
| $\lambda_x$ ($\mu m$) | $\sigma_{P}$ (GPa) | $\sigma_{\Delta T}$ (°C) | $\sigma_{H}$ ($\mu m$) |
|---|---|---|---|
| 15.0 | 0.0682 | 3.37 | 0.0374 |
| 22.5 | 0.0578 | 2.66 | 0.0361 |
| 30.0 | 0.0428 | 2.37 | 0.0352 |
Axial Direction Wavelength ($\lambda_y$): Similarly, increasing the wavelength $\lambda_y$ across the face width of the helical gear tooth also dampens fluctuations. It also weakens the pronounced side leakage or “necking” effect at the ends of the finite contact line, leading to a more uniform pressure and film thickness distribution across the face width. The corresponding standard deviations along the $y$-direction centerline ($x=0$) also show a clear decreasing trend.
| $\lambda_y$ (mm) | $\sigma_{P}$ (GPa) | $\sigma_{\Delta T}$ (°C) | $\sigma_{H}$ ($\mu m$) |
|---|---|---|---|
| 1.5 | 0.0574 | 2.31 | 0.0341 |
| 2.0 | 0.0370 | 1.84 | 0.0270 |
| 2.5 | 0.0361 | 1.69 | 0.0262 |
Comparing the effects, roughness in the entrainment direction ($x$) generally causes more severe disturbance to the pressure and temperature fields than roughness in the axial direction ($y$). This is because the primary flow of lubricant is in the entrainment direction, and asperities aligned perpendicular to this flow (i.e., with a short $\lambda_x$) create more significant obstacles and flow restrictions.
Effect of Roughness Amplitude
The amplitude of the composite surface roughness $A_{xy}$ is the most direct measure of surface finish quality. As expected, increasing the roughness amplitude dramatically amplifies the fluctuations in all lubrication parameters. The pressure spikes become sharper, the temperature fluctuations more extreme, and the film thickness valleys deeper. The quantitative impact is severe, as evidenced by the rapid growth in the standard deviation of fluctuations for both $x$ and $y$ profiles.
| $A_{xy}$ ($\mu m$) | $\sigma_{P_x}$ (GPa) | $\sigma_{P_y}$ (GPa) | $\sigma_{\Delta T_x}$ (°C) | $\sigma_{\Delta T_y}$ (°C) | $\sigma_{H_x}$ ($\mu m$) | $\sigma_{H_y}$ ($\mu m$) |
|---|---|---|---|---|---|---|
| 0.025 | 0.0337 | 0.0483 | 1.59 | 2.60 | 0.0173 | 0.0258 |
| 0.05 | 0.0625 | 0.0735 | 3.12 | 4.36 | 0.0348 | 0.0370 |
| 0.10 | 0.1040 | 0.1180 | 6.71 | 8.09 | 0.0700 | 0.0670 |
The percentage increases in $\sigma$ are substantial, underscoring that controlling the peak roughness amplitude is paramount for ensuring stable lubrication in high-speed helical gears.
Influence of Dimensionless Parameters G, U, and W
Beyond the roughness parameters themselves, the operational regime defined by G, U, and W fundamentally alters how roughness affects the helical gear contact.
Material Parameter (G): An increase in G, which primarily reflects an increase in the pressure-viscosity coefficient $\alpha$ of the lubricant, leads to a stiffer, more viscous film under pressure. While this generally increases the central film thickness and maximum pressure, it also intensifies the pressure and temperature fluctuations caused by a given roughness profile. The lubricant’s strong piezoviscous response amplifies the local pressure perturbations at asperities. The standard deviations $\sigma_P$ and $\sigma_{\Delta T}$ increase sharply with G, while $\sigma_H$ remains relatively constant, indicating that the film thickness “follows” the roughness contour more rigidly under high-G conditions.
| $G$ | $\sigma_{P}$ (GPa) | $\sigma_{\Delta T}$ (°C) | $\sigma_{H}$ ($\mu m$) |
|---|---|---|---|
| 2000 | 0.0195 | 1.19 | 0.0354 |
| 2800 | 0.0469 | 3.56 | 0.0351 |
| 3600 | 0.0893 | 6.58 | 0.0349 |
Speed Parameter (U): Increasing U, equivalent to increasing entrainment speed or lubricant viscosity, is highly beneficial for EHL film formation. It results in a much thicker lubricant film and higher temperatures due to increased shear. Interestingly, a thicker film tends to “smother” the roughness to some extent. The pressure fluctuations become less severe relative to the increased central pressure, and their standard deviation $\sigma_P$ actually decreases with increasing U. Although the absolute temperature rise increases, the relative fluctuation $\sigma_{\Delta T}$ also decreases, suggesting a more thermally stable film. The film thickness fluctuation $\sigma_H$ decreases slightly as the film becomes thicker relative to the fixed roughness amplitude.
| $U / U^*$ | $\sigma_{P}$ (GPa) | $\sigma_{\Delta T}$ (°C) | $\sigma_{H}$ ($\mu m$) |
|---|---|---|---|
| 0.6 | 0.0517 | 3.39 | 0.0349 |
| 1.0 | 0.0352 | 2.76 | 0.0340 |
| 1.4 | 0.0253 | 2.14 | 0.0331 |
Load Parameter (W): An increase in W, meaning higher transmitted load on the helical gear teeth, has a detrimental effect on roughness sensitivity. The contact pressure rises significantly, and the film thickness decreases slightly. The reduced film thickness brings the surfaces closer together, effectively amplifying the relative height of the roughness asperities. Consequently, the pressure and temperature fluctuations become much more violent. The standard deviations $\sigma_P$ and $\sigma_{\Delta T}$ exhibit a very strong increase with W, highlighting that heavily loaded helical gear contacts are particularly vulnerable to the destabilizing effects of surface roughness.
| $W / W^*$ | $\sigma_{P}$ (GPa) | $\sigma_{\Delta T}$ (°C) | $\sigma_{H}$ ($\mu m$) |
|---|---|---|---|
| 0.5 | 0.0083 | 0.66 | 0.0383 |
| 1.0 | 0.0240 | 2.02 | 0.0352 |
| 1.5 | 0.0480 | 4.17 | 0.0354 |
Conclusion
This investigation into the thermo-elastohydrodynamic lubrication of high-speed helical gears, incorporating the critical effect of surface roughness, yields several important conclusions for the design and operation of such components in aero-engines:
- Significance of Roughness: When the amplitude of gear tooth surface roughness is comparable to the nominal EHL film thickness (sub-micron level), its influence is dominant and cannot be neglected. It causes severe oscillations in pressure and temperature, and creates localized valleys in film thickness, especially near the contact exit and edges, significantly increasing the risk of lubrication failure.
- Roughness Parameter Effects: The destabilizing effect is most strongly governed by the roughness amplitude. Increasing amplitude linearly exacerbates fluctuations in all film parameters. Reducing the spatial frequency of asperities (increasing wavelength) in either the entrainment or axial direction alleviates these fluctuations. Roughness oriented across the lubricant flow direction (short $\lambda_x$) is particularly detrimental.
- Dimensionless Parameter Interplay: The operational regime defined by G, U, and W modulates the system’s sensitivity to roughness.
- High Material Parameter (G) conditions, while promoting thicker films, intensify pressure and temperature fluctuations caused by roughness.
- High Speed Parameter (U) conditions are beneficial, as increased entrainment generates thicker films that partially suppress roughness-induced pressure fluctuations and lead to more stable thermal conditions.
- High Load Parameter (W) conditions are detrimental, as the reduced film thickness magnifies the relative effect of roughness, leading to dramatically increased pressure and temperature fluctuations.
Practical Implications for Helical Gear Design: For high-speed helical gears operating under severe conditions, particular attention must be paid to surface finish. The primary goal should be to minimize the peak-to-valley roughness amplitude ($A_{xy}$). Furthermore, manufacturing processes should aim to avoid a high spatial frequency of sharp asperities (i.e., strive for a longer effective wavelength). While increasing speed (U) is beneficial, it is often dictated by engine performance. Increasing the equivalent radius of curvature (e.g., through gear geometry optimization) to reduce the load parameter (W) can be an effective strategy to mitigate roughness sensitivity. Although selecting a lubricant with a higher pressure-viscosity coefficient (increasing G) boosts nominal film thickness, it comes at the cost of increased pressure fluctuation severity and must be carefully considered. In summary, a holistic approach combining superior surface finishing with optimal gear geometry and operational parameter selection is essential for ensuring reliable elastohydrodynamic lubrication and long service life for high-speed helical gears in aviation applications.