The operational temperature of cylindrical gear transmissions, especially under demanding conditions, is a critical factor influencing their service life, efficiency, and reliability. In aerospace applications, such as helicopter transmission systems, gears are frequently subjected to high speeds and heavy loads. A common design strategy to enhance load-carrying capacity and reduce stress is the use of High Contact Ratio (HCR) spur cylindrical gears, defined as those with a contact ratio greater than 2.0. Unlike standard spur cylindrical gears, HCR gears maintain at least two pairs of teeth in contact simultaneously during operation, leading to a longer total contact line and a reduced average load per unit length. However, the power losses generated during meshing, primarily in the form of frictional heat, must be effectively managed to prevent failures like scuffing or thermal fatigue. Jet oil lubrication is the predominant method for cooling and lubricating such high-performance cylindrical gears. Therefore, a comprehensive understanding of the temperature field and convective heat transfer characteristics of HCR spur cylindrical gears under jet lubrication is of paramount engineering significance.
This analysis focuses on investigating the steady-state temperature field of HCR spur cylindrical gear pairs using a Computational Fluid Dynamics (CFD) based approach. The methodology involves a conjugate heat transfer simulation that couples the fluid flow of the lubricant-air mixture with the solid conduction within the cylindrical gears. Key aspects such as the distribution of temperature and convective heat transfer coefficient (HTC) on the gear tooth surfaces are obtained. Furthermore, a parametric study is conducted to evaluate the influence of lubrication parameters, operational conditions, and cylindrical gear design parameters on the thermal behavior. Finally, experimental measurements from a gear test rig are presented to validate the trends observed in the numerical simulations.
1. Theoretical Foundation for CFD Analysis
The flow within a gearbox, involving a complex interaction between air, oil droplets, and the rotating solid components, is governed by the fundamental laws of fluid dynamics and heat transfer. The following sections outline the core theoretical models employed.
1.1 Governing Equations
The fluid flow is described by the Reynolds-Averaged Navier-Stokes (RANS) equations, which encompass conservation of mass, momentum, and energy. For a fluid domain, these are expressed as:
Mass Conservation (Continuity):
$$ \frac{\partial \rho_f}{\partial t} + \nabla \cdot (\rho_f \vec{u}) = 0 $$
Momentum Conservation:
$$ \frac{\partial (\rho_f u_i)}{\partial t} + \nabla \cdot (\rho_f u_i \vec{u}) = \nabla \cdot (\mu \, \nabla u_i) – \frac{\partial p}{\partial x_i} + S_i $$
Energy Conservation:
$$ \frac{\partial (\rho_f T)}{\partial t} + \nabla \cdot (\rho_f \vec{u} T) = \nabla \cdot \left( \frac{k_f}{C_p} \nabla T \right) + S_T $$
where \( \rho_f \) is the fluid density, \( t \) is time, \( \vec{u} \) is the velocity vector, \( u_i \) are its components, \( \mu \) is the dynamic viscosity, \( p \) is pressure, \( S_i \) is a momentum source term, \( T \) is temperature, \( k_f \) is thermal conductivity, \( C_p \) is specific heat capacity, and \( S_T \) is the viscous dissipation term.
1.2 Multiphase Flow Modeling: The VOF Method
To capture the interface between the lubricating oil and the air within the gearbox, the Volume of Fluid (VOF) method is employed. This model is particularly suitable for tracking the transient motion of distinct, immiscible fluids. In each computational cell, the volume fractions of oil (\( \alpha_{oil} \)) and air (\( \alpha_{air} \)) sum to unity:
$$ \alpha_{oil} + \alpha_{air} = 1 $$
The properties in each cell, such as density and viscosity, are calculated as weighted averages based on these volume fractions.
1.3 Heat Generation Model for Cylindrical Gears
The primary heat source in the system is the frictional power loss from the meshing cylindrical gear teeth. A comprehensive model based on the Anderson and Loewenthal method is used to estimate the total power loss \( Q \). This includes rolling, sliding, and windage losses. The heat generation on the tooth flank surfaces is then implemented as a volumetric heat source in the solid gear domains adjacent to the meshing interfaces.
The average sliding power loss \( P_s \) and rolling power loss \( P_r \) are critical components. They can be expressed in a generalized form that highlights their dependence on operational parameters:
$$ P_s \propto f \cdot F_n \cdot \bar{V}_s $$
$$ P_r \propto \bar{V}_t \cdot \bar{h} \cdot b \cdot e_p $$
where \( f \) is the coefficient of friction, \( F_n \) is the average normal load, \( \bar{V}_s \) is the average sliding velocity, \( \bar{V}_t \) is the average rolling velocity, \( \bar{h} \) is the lubricant film thickness, \( b \) is the face width, and \( e_p \) is the contact ratio. The total heat flux \( Q \) is partitioned between the driving and driven cylindrical gears based on their thermal properties and tangential velocities at the meshing point.
2. CFD Simulation Model Setup
2.1 Geometry and Parameters
The analysis focuses on a pair of spur cylindrical gears with high contact ratio. The key geometric and operational parameters are summarized in the table below.
| Parameter | Driving Gear (Large) | Driven Gear (Small) |
|---|---|---|
| Normal Module, \( m_n \) (mm) | 3.25 | 3.25 |
| Number of Teeth, \( z \) | 32 | 25 |
| Face Width, \( b \) (mm) | 16.0 | 16.5 |
| Pressure Angle, \( \alpha \) (°) | 20 | 20 |
| Addendum Coefficient, \( h_a^* \) | 1.32 | 1.32 |
| Contact Ratio, \( e_p \) | 2.2 (High Contact Ratio) | |
| Rotational Speed, \( n \) (rpm) | 1500 | 1920 |
| Load Level | 9 (Heavy Load) | |
| Material Property | Value (20CrMnMoA Steel) |
|---|---|
| Thermal Conductivity, \( k_s \) (W/m·K) | 46 |
| Specific Heat Capacity, \( C_{p,s} \) (J/kg·K) | 470 |
| Density, \( \rho_s \) (kg/m³) | 7850 |
| Lubricant Property (Shell 555) | Value |
|---|---|
| Density at 15.6°C, \( \rho_{oil} \) (kg/m³) | 993 |
| Kinematic Viscosity at 37.8°C, \( \nu_{oil} \) (mm²/s) | 29 |
| Kinematic Viscosity at 98.9°C, \( \nu_{oil} \) (mm²/s) | 5.4 |
2.2 Computational Domain and Meshing
A 3D model of the gearbox encompassing the two spur cylindrical gears, the housing, the oil jet nozzle, and the outlet is created. To efficiently model the rotation of the gears without resorting to transient sliding meshes for steady-state analysis, the Multiple Reference Frame (MRF) method is adopted. In this approach, cylindrical fluid domains surrounding each cylindrical gear are defined as rotating reference frames with their respective speeds. The interaction between these rotating zones and the stationary fluid domain (housing) is handled through interface boundaries. The mesh is primarily composed of unstructured tetrahedral cells, with local refinement applied in the gear meshing zone, near the jet impingement area, and on the tooth surfaces to accurately resolve the flow and thermal boundary layers.
2.3 Boundary Conditions and Solution Strategy
The volumetric heat source calculated from the power loss model is applied to a thin layer of elements on the active flanks of both spur cylindrical gears. The boundary conditions for the fluid domain are listed below:
| Boundary | Type | Specification |
|---|---|---|
| Oil Jet Inlet | Velocity Inlet | Velocity = 40 m/s, Oil Temperature = 60°C |
| Gearbox Outlet | Pressure Outlet | Gauge Pressure = 0 Pa (Atmospheric) |
| Gear Surfaces | Coupled Wall | Conjugate Heat Transfer |
| Gearbox Housing Walls | Wall | Natural Convection: HTC = 5 W/m²·K, Ambient Temp. = 26.85°C |
| Rotating Fluid Zones | MRF Zones | Speeds: 1500 rpm (Driving), 1920 rpm (Driven) |
The simulation solves the continuity, momentum, energy, and volume fraction equations. A pressure-based solver with the SIMPLE algorithm for pressure-velocity coupling is used. Second-order discretization schemes are employed for higher accuracy. The convergence is monitored through the residuals of all solved equations.
3. Numerical Results and Parametric Analysis
3.1 Temperature and Convective Heat Transfer Distribution
The steady-state temperature distribution on the surfaces of the HCR spur cylindrical gears reveals distinct patterns. For both gears, the highest temperatures are located on the tooth flanks, specifically in the region between the pitch line and the tooth tip. This area experiences significant sliding friction, leading to greater heat generation. The temperature distribution is symmetric about the mid-plane of the face width, with temperatures being highest at the center and gradually decreasing towards the gear sides due to better heat dissipation at the exposed ends.
The convective heat transfer coefficient (HTC) distribution provides insight into the cooling effectiveness. The HTC is generally higher on the driven cylindrical gear (smaller, faster) compared to the driving gear. On each gear, the HTC increases with radius on the gear side faces, reaching a maximum near the tip circle where the peripheral speed is highest. The mesh zone, where oil is entrained and squeezed between the teeth, also shows locally elevated HTC values.
3.2 Influence of Operating and Design Parameters
A controlled parametric study was conducted by varying one parameter at a time while keeping others at their baseline values. The effects on the maximum tooth flank temperature of the driving spur cylindrical gear are summarized conceptually below. It is important to note that the driven gear follows qualitatively similar trends.
| Parameter | Trend of Maximum Gear Temperature | Primary Physical Reason |
|---|---|---|
| Oil Inlet Temperature | Increases linearly | Higher initial energy input to the system. |
| Oil Jet Flow Rate | Decreases, then asymptotes | Improved cooling with more oil; diminishing returns and increased churning loss at very high flow. |
| Gear Rotational Speed | Increases, then slightly decreases at very high speed | Frictional heat increases with speed; at very high speeds, reduced contact time and altered flow may enhance cooling. |
| Transmitted Load (Torque) | Increases | Increased normal force leads to higher frictional power loss. |
| Gear Face Width | Decreases | Increased heat dissipation surface area outweighs the increase in total frictional heat. |
| Pressure Angle | Decreases | Larger pressure angle can reduce sliding friction and relative sliding velocity. |
| Contact Ratio (via Addendum Coeff.) | Increases | More teeth in contact increases the total frictional heat generation area. |
The relationship between load \( F \), speed \( N \), and the resulting maximum temperature rise \( \Delta T_{max} \) can be conceptually framed. The frictional heat generation has a strong dependence on load and sliding velocity, which is itself a function of gear geometry and speed. A highly simplified representation focusing on the sliding loss component could be:
$$ \Delta T_{max} \propto \frac{P_{heat}}{h_{conv} \cdot A} \approx \frac{\mu \cdot F_n \cdot V_s}{h_{conv} \cdot A} $$
where \( \mu \) is friction coefficient, \( F_n \) is normal load, \( V_s \) is sliding velocity, \( h_{conv} \) is an average convective heat transfer coefficient, and \( A \) is the effective cooling area. This illustrates why both increasing load and speed (which increases \( V_s \)) tend to raise temperature, while parameters that improve \( h_{conv} \) or \( A \) (like higher oil flow or face width) help lower it.
4. Experimental Validation and Discussion
To validate the trends observed in the CFD analysis, experimental tests were performed on a CL-100 type closed-loop power recirculating gear test rig. The test spur cylindrical gear pair had geometry matching the simulation model. Thermocouples were embedded near the tooth flank surface of both gears to measure bulk temperature. The gearbox was lubricated via a jet oriented at the mesh entry zone, with oil pre-heated to the desired inlet temperature. Tests were run under various load and oil temperature conditions until a steady-state temperature was reached.
The experimental results confirmed the core findings from the simulation. Firstly, a clear correlation between increased transmitted load and higher measured gear temperature was established, aligning with the CFD predictions. Secondly, comparative tests between a standard contact ratio (approx. 1.7) spur cylindrical gear pair and the HCR (approx. 2.2) pair were conducted. The data demonstrated that under identical operating conditions, the HCR gears consistently operated at a higher temperature. This difference became more pronounced as the load increased. This experimental evidence supports the simulation finding that the increased number of simultaneous contact points in an HCR spur cylindrical gear, while beneficial for stress reduction, leads to greater total frictional heat generation that must be accounted for in the thermal design.
The quantitative comparison between simulated and measured absolute temperatures showed acceptable agreement, with the simulation capturing the correct trends and order of magnitude. Discrepancies can be attributed to simplifications in the heat generation model, uncertainties in material properties, and the challenge of exactly replicating boundary conditions like external housing convection in the simulation.
5. Conclusion
This integrated CFD and experimental study provides detailed insights into the thermal behavior of high contact ratio spur cylindrical gears under jet lubrication. The CFD model, utilizing the VOF method for multiphase flow and the MRF approach for rotation, successfully predicted the characteristic temperature and convective heat transfer coefficient distribution on the gear teeth. The parametric analysis systematically quantified the influence of key variables: gear temperature increases with higher oil inlet temperature, rotational speed, transmitted load, and contact ratio; it decreases with larger oil flow rate (up to a point), increased face width, and larger pressure angle. The experimental validation confirmed the critical trend of rising temperature with load and, importantly, provided empirical evidence that HCR spur cylindrical gears run hotter than their standard contact ratio counterparts, a phenomenon that intensifies under heavier loads. These findings underscore the necessity of incorporating detailed thermal analysis early in the design process of high-performance spur cylindrical gear transmissions, particularly when employing high contact ratio designs for their mechanical advantages. Effective thermal management through optimized lubrication parameters and gear geometry is essential to fully leverage the benefits of HCR spur cylindrical gears in demanding applications.