The efficiency and reliability of gear transmissions, particularly cylindrical gears operating under high-speed and heavy-load conditions, are profoundly dependent on effective lubrication. The lubrication system must accomplish two primary objectives: reducing friction and wear at the contact interfaces and dissipating the substantial heat generated during operation. Oil injection lubrication is the predominant method employed for such demanding applications. However, the efficacy of this method is critically influenced by complex aerodynamic interactions. The high-speed rotation of the gear pair generates a circumambient air flow field, which can create a “gas barrier” or “air shield” effect, impeding the oil jet from reaching the critical meshing zone. This interference can lead to insufficient lubrication, increased localized temperatures, and accelerated wear. For conventional spur or helical cylindrical gears, significant research has been conducted to understand this phenomenon and optimize parameters like nozzle position, injection angle, and speed. Yet, the lubrication challenge becomes even more intricate for gears with complex tooth geometries.
The Variable Hyperbolic Circular-Arc-Tooth-Trace (VH-CATT) cylindrical gear represents one such advanced geometry. It distinguishes itself from standard cylindrical gears through its tooth trace, which is a circular arc, and its transverse profile, which consists of a family of variable hyperbolas enveloped from a base involute. This unique design offers superior load-bearing capacity, improved meshing characteristics, and eliminates axial thrust forces, making it highly attractive for advanced transmission systems. However, this very complexity in tooth surface geometry fundamentally alters the surrounding airflow patterns during operation. The curved tooth trace and variable profile disrupt the otherwise more predictable flow fields generated by straight-tooth cylindrical gears. Consequently, the established rules for optimizing oil injection for standard gears may not be directly applicable. There is a pressing need to investigate how the special aerodynamic behavior induced by the VH-CATT geometry interacts with the injected oil stream and to establish new guidelines for its lubrication system design. The primary objective of this work is to employ Computational Fluid Dynamics (CFD) to dissect the complex gas-liquid two-phase flow during the oil injection lubrication process for VH-CATT cylindrical gears. We aim to identify the interference mechanisms between the high-speed air flow and the oil jet, determine the optimal injection parameters—especially the injection angle—and systematically analyze the influence of various operational and geometric parameters on the resulting lubrication performance of the tooth surface.
The analysis of VH-CATT cylindrical gears begins with their mathematical definition, which originates from a specific machining process using a dual-cutter head. The coordinate systems for the gear generation are established, where one frame is attached to the cutter head and another to the gear blank. Through the principles of gearing and coordinate transformation, the tooth surface position vector, \(\mathbf{r}_i^{(d)}\), and the unit normal vector, \(\mathbf{n}_i^{(d)}\), for the pinion (i=p) and gear (i=g) in the gear-fixed coordinate system can be derived as:
$$
\begin{cases}
\mathbf{r}_i^{(d)}(u_i, \theta_i, \phi_i) = \mathbf{M}_{di}(\phi_i) \cdot \mathbf{M}_{i1}(\theta_i) \cdot \mathbf{r}_i^{(1)}(u_i, \theta_i) \\
\mathbf{n}_i^{(d)}(\theta_i, \phi_i) = \mathbf{L}_{di}(\phi_i) \cdot \mathbf{L}_{i1}(\theta_i) \cdot \mathbf{n}_i^{(1)}(\theta_i)
\end{cases}
$$
Here, \(u_i\) and \(\theta_i\) are the surface parameters of the generating tool, \(\phi_i\) is the rotation angle of the gear during generation, \(\mathbf{M}\) and \(\mathbf{L}\) represent coordinate transformation matrices for position and unit normal vectors respectively, and the superscripts denote the coordinate system. Solving these equations within the parameter ranges yields the point cloud data for the concave and convex tooth surfaces, which can then be used to construct a precise three-dimensional solid model essential for CFD analysis.
The contact characteristics of VH-CATT cylindrical gears are also distinct. They engage in point contact under no load, which expands into an elliptical contact area under load due to elastic deformation. The size and orientation of this contact ellipse are governed by the principal curvatures and directions of the mating tooth surfaces at the theoretical contact point. For a convex pinion surface (\(\Sigma_p\)) and a concave gear surface (\(\Sigma_g\)), with principal curvatures \(K_1^p, K_2^p\) and \(K_1^g, K_2^g\), and corresponding direction vectors \(\mathbf{e}_1^p, \mathbf{e}_2^p\) and \(\mathbf{e}_1^g, \mathbf{e}_2^g\), the contact ellipse is defined by:
$$ A x^2 + B y^2 = \delta $$
where \(\delta\) is the deformation at the contact point, and the coefficients are:
$$ A = \frac{1}{4}\left[ (K_{\Sigma}^p – K_{\Sigma}^g) – \sqrt{(g_p – g_g)^2 + 4 g_p g_g \cos^2 \sigma} \right] $$
$$ B = \frac{1}{4}\left[ (K_{\Sigma}^p – K_{\Sigma}^g) + \sqrt{(g_p – g_g)^2 + 4 g_p g_g \cos^2 \sigma} \right] $$
with \(K_{\Sigma}^i = K_1^i + K_2^i\), \(g_i = K_1^i – K_2^i\) for \(i = p, g\), and \(\sigma\) being the angle between \(\mathbf{e}_1^p\) and \(\mathbf{e}_1^g\). The semi-major axis \(a\), semi-minor axis \(b\), and the orientation angle \(\alpha_p\) of the ellipse relative to \(\mathbf{e}_1^p\) are:
$$ a = \sqrt{\delta / A}, \quad b = \sqrt{\delta / B}, \quad \tan(2\alpha_p) = \frac{2g_p g_g \sin(2\sigma)}{g_p^2 – g_g^2 + 2g_p g_g \cos(2\sigma)} $$
Understanding this contact ellipse is crucial for defining the target “inspection surface” or “detection zone” in our CFD model to quantitatively assess lubrication performance where it matters most—within the loaded contact path.
To simulate the oil injection lubrication process, a robust mathematical framework for multiphase flow is required. The Volume of Fluid (VOF) model within the Eulerian multiphase framework is well-suited for this task, as it tracks the volume fraction of each immiscible fluid (air and oil) throughout the computational domain while solving a single set of momentum equations. Let \(\alpha_r\) represent the volume fraction of the \(r^{th}\) fluid phase, with the constraint:
$$ \sum_{r=1}^{N} \alpha_r = 1 $$
where \(N=2\). The mixture density \(\rho\) and viscosity \(\mu\) are calculated as weighted averages based on these volume fractions. The governing equations are the continuity and momentum equations for the mixture:
$$ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{U}) = 0 $$
$$ \frac{\partial (\rho \mathbf{U})}{\partial t} + \nabla \cdot (\rho \mathbf{U} \mathbf{U}) = -\nabla p + \nabla \cdot [\mu (\nabla \mathbf{U} + \nabla \mathbf{U}^T)] + \rho \mathbf{g} + \mathbf{F} $$
Here, \(\mathbf{U}\) is the velocity vector, \(p\) is pressure, \(\mathbf{g}\) is gravity, and \(\mathbf{F}\) represents other external forces. To account for the turbulent nature of the flow at high rotational speeds, the standard \(k\)-\(\epsilon\) turbulence model is employed. It introduces two transport equations for the turbulent kinetic energy \(k\) and its dissipation rate \(\epsilon\):
$$ \frac{\partial (\rho k)}{\partial t} + \frac{\partial (\rho k u_i)}{\partial x_i} = \frac{\partial}{\partial x_j} \left[ \left( \mu + \frac{\mu_t}{\sigma_k} \right) \frac{\partial k}{\partial x_j} \right] + G_k – \rho \epsilon $$
$$ \frac{\partial (\rho \epsilon)}{\partial t} + \frac{\partial (\rho \epsilon u_i)}{\partial x_i} = \frac{\partial}{\partial x_j} \left[ \left( \mu + \frac{\mu_t}{\sigma_\epsilon} \right) \frac{\partial \epsilon}{\partial x_j} \right] + C_{1\epsilon} \frac{\epsilon}{k} G_k – C_{2\epsilon} \rho \frac{\epsilon^2}{k} $$
where \(\mu_t = \rho C_\mu k^2 / \epsilon\) is the turbulent viscosity, and \(G_k\) represents the generation of turbulent kinetic energy due to mean velocity gradients. The model constants \((C_\mu, \sigma_k, \sigma_\epsilon, C_{1\epsilon}, C_{2\epsilon})\) have standard values.
Capturing the motion of the rotating cylindrical gears is achieved through dynamic meshing techniques and User-Defined Functions (UDFs). The spring-based smoothing method is used to update the mesh in deforming zones. The displacement of a node \(\Delta x_i\) is calculated based on forces from all connected edges, analogous to a network of springs:
$$ \Delta x_i = \frac{\sum_{j}^{n_i} k_{ij} \Delta x_j}{\sum_{j}^{n_i} k_{ij}} $$
where \(k_{ij}\) is the spring constant between nodes \(i\) and \(j\), and \(n_i\) is the number of neighboring nodes. The node position is updated iteratively: \(x_i^{n+1} = x_i^n + \Delta x_i^{converged}\). To prevent excessive mesh distortion in the narrow gap of the meshing zone, a small offset (e.g., 0.5°) is applied to expand the initial clearance between the teeth in the 3D fluid domain model.
The specific case study focuses on the lubrication of a VH-CATT cylindrical gear pair with oil injected at the meshing-in side. Key geometric and operating parameters are summarized in the following table:
| Parameter | Pinion | Gear |
|---|---|---|
| Number of Teeth, \(z\) | 21 | 29 |
| Module, \(m\) (mm) | 4 | 4 |
| Pressure Angle, \(\alpha_0\) (°) | 20 | 20 |
| Rotational Speed, \(n\) (rpm) | 6000 | 4345 |
| Center Distance, \(a\) (mm) | 100 | |
| Face Width, \(B\) (mm) | 40 | 40 |
Other fixed parameters include a nozzle diameter \(D\) of 5 mm, lubricant density \(\rho_{oil}\) of 960 kg/m³, and a baseline injection velocity \(v\) of 60 m/s. The key variables for optimization are the injection angle \(\alpha\) (defined as the angle between the nozzle axis and the common tangent to the pitch circles) and the injection height \(d\) (the distance from the nozzle to the gear center line).
The first critical step is to understand the aerodynamic interference and find the optimal injection angle. A steady-state, single-phase (air-only) simulation of the gears rotating at their operational speeds reveals the complex streamlines around the VH-CATT cylindrical gears. The curved tooth traces create distinct vortices and flow patterns not seen in simpler gears. By analyzing this flow field, particularly in the region above the meshing zone, a trajectory of weakest airflow resistance can be identified. Conceptually, injecting oil along this path would minimize the deflection and breakup of the oil jet by the opposing air flow. For the specific gear geometry studied, this optimal trajectory was found to form an angle of approximately 10.73° relative to the pitch line tangent, pointing towards the driving pinion.
To validate this finding, a series of transient, multiphase (oil-air) simulations were conducted with different injection angles: -5°, 0°, 5°, and the proposed 10.73°. The performance metrics were the oil volume fraction distribution on a defined inspection surface covering the meshing path and the resulting oil pressure in the meshing zone. A higher and more concentrated oil volume fraction indicates better lubricant delivery, while a larger pressure difference (between positive and negative pressure zones) near the meshing point indicates a stronger “pumping” or entrainment effect, drawing oil into the contact. The results strongly validated the theoretical prediction. The configuration with \(\alpha = 10.73^\circ\) showed the most concentrated oil distribution on the tooth flank and yielded the highest average oil volume fraction in the inspection zone during the meshing cycle. The pressure statistics further confirmed this, showing the largest positive pressure, the largest negative pressure (in absolute value), and consequently, the greatest pressure differential for the 10.73° angle, signifying the strongest capability to draw lubricant into the contact interface of these cylindrical gears.
| Injection Angle, \(\alpha\) (°) | Avg. Max Oil Volume Fraction | Max Positive Pressure (MPa) | Max Negative Pressure (MPa) | Pressure Differential (MPa) |
|---|---|---|---|---|
| -5.0 | 0.42 | 0.65 | 0.32 | 0.97 |
| 0.0 | 0.58 | 0.78 | 0.45 | 1.23 |
| 5.0 | 0.68 | 0.85 | 0.52 | 1.37 |
| 10.73 | 0.75 | 0.96 | 0.61 | 1.57 |
With the optimal angle established, the influence of other key operational parameters was systematically investigated by holding \(\alpha\) at 10.73°. The effect of injection velocity \(v\) was studied from 30 m/s to 80 m/s. Unsurprisingly, higher injection momentum helped the oil jet penetrate the air barrier more effectively. The average oil volume fraction on the inspection surface increased significantly as \(v\) rose from 30 to 60 m/s. However, beyond 60 m/s, the rate of improvement diminished, suggesting a point of diminishing returns where increasing pump power yields minimal additional lubrication benefit. The pressure differential in the meshing zone followed a similar trend, increasing steadily with injection velocity, which promotes better oil entrainment.
The injection height \(d\) is another crucial design parameter. Simulations varied \(d\) from 30 mm to 55 mm. A shorter distance (e.g., 30-40 mm) resulted in a more focused oil impingement pattern and higher localized oil volume fractions and pressures. As the nozzle was moved farther away, the oil jet had more time to disperse and lose momentum due to air drag before reaching the tooth flank, leading to a more spread-out but thinner oil film and a noticeable reduction in the meshing zone pressure differential. This indicates that, within practical limits related to gearbox geometry and avoiding physical interference, a closer nozzle placement is beneficial for lubricating these cylindrical gears.
Finally, the impact of the gear rotational speed \(n_1\) (pinion speed) was analyzed from 3000 rpm to 8000 rpm. This parameter directly controls the intensity of the “gas barrier” effect. At lower speeds (3000-4000 rpm), the air flow is less disruptive, allowing for good oil deposition. As speed increased, the centrifugal force acting on oil droplets and the velocity of the opposing air stream both rose substantially. This led to a pronounced decrease in the average oil volume fraction adhering to the tooth surface during meshing. Furthermore, the high-speed rotating air field altered the pressure distribution, reducing the beneficial pressure differential in the meshing zone. This inverse relationship between rotational speed and lubrication effectiveness highlights the heightened cooling and lubrication demands for very high-speed VH-CATT cylindrical gear applications.
| Parameter Varied | Trend on Lubrication Performance | Key Observation |
|---|---|---|
| Injection Angle, \(\alpha\) | Performance peaks at a specific angle. | Optimal angle was 10.73° towards the pinion. |
| Injection Velocity, \(v\) | Improves up to a limit. | Significant gain up to ~60 m/s; diminishing returns thereafter. |
| Injection Height, \(d\) | Decreases with increasing distance. | Shorter distances (30-40 mm) yield more focused and effective lubrication. |
| Gear Speed, \(n_1\) | Decreases with increasing speed. | Higher speeds severely reduce oil film coverage and entrainment pressure. |
In conclusion, this comprehensive CFD-based investigation into the oil injection lubrication of Variable Hyperbolic Circular-Arc-Tooth-Trace cylindrical gears has yielded critical insights for system design. The unique tooth geometry of these cylindrical gears creates a specific aerodynamic environment that necessitates a tailored lubrication strategy. The most significant finding is the existence of an optimal injection angle (10.73° towards the driving pinion for the studied case), which aligns the oil jet with the path of least airflow resistance, maximizing lubricant delivery to the meshing zone. Furthermore, the analysis provides quantitative guidance: lubrication effectiveness is enhanced by employing a sufficiently high injection velocity (e.g., ~60 m/s), placing the nozzle as close as practically feasible to the gears, and accounting for the significantly increased lubrication challenge at higher operational speeds. These findings establish a foundational understanding and provide practical design guidelines for developing efficient lubrication and cooling systems for high-performance VH-CATT cylindrical gear transmissions, ensuring their reliability and longevity in demanding applications.