The pursuit of higher power density and efficiency in aerospace, marine, and heavy-duty industrial transmissions has consistently driven the development of advanced gear systems. Among these, herringbone gears stand out due to their superior characteristics, including compact structure, smooth operation, high load-carrying capacity, and inherent axial force balance. In high-speed transmission applications, power loss becomes a critical performance metric, directly influencing overall system efficiency, thermal management, and energy consumption. Total power loss in a gearbox comprises load-dependent meshing losses and load-independent losses, primarily churning (oil drag) and windage losses. Windage loss is defined as the power dissipated due to the resistance exerted by the surrounding fluid—air or an air-oil mixture—on the rotating gears. When the pitch line velocity exceeds approximately 50 m/s, the agitation of the surrounding fluid intensifies significantly, and windage loss can become the dominant contributor to the total no-load or even part-load power loss. Therefore, a detailed understanding and accurate prediction of windage losses in high-speed herringbone gear pairs are essential for optimal design and efficiency enhancement.
This study focuses on the transient aerodynamic analysis of high-speed herringbone gear pairs operating in air. A three-dimensional computational fluid dynamics (CFD) model is developed to simulate the complex, turbulent flow field around the meshing gears. The analysis aims to elucidate the underlying mechanisms of windage loss, quantify the power dissipation, and systematically investigate the influence of key operational and geometrical parameters, namely rotational speed, direction of rotation, and helix angle. The findings provide valuable insights for the design and optimization of efficient herringbone gear transmissions for high-speed applications.
1. Computational Methodology and Numerical Modeling
The accurate simulation of windage losses requires a robust numerical model that can capture the transient interaction between the rotating gear geometries and the surrounding fluid. The methodology encompasses geometric modeling, domain discretization, the application of appropriate physical models, and the use of dynamic mesh techniques to handle gear motion.
1.1 Geometric Model and Computational Domain
The analysis is performed on a herringbone gear pair with standard design parameters. The basic geometrical specifications are summarized in Table 1. The pinion (driver) rotates at a nominal speed of 9,375 rpm, and the wheel (driven) at 7,500 rpm, resulting in a pitch line velocity of approximately 60.7 m/s, placing it firmly in the high-speed regime where windage effects are pronounced.
| Parameter | Wheel | Pinion |
|---|---|---|
| Number of Teeth, \(z\) | 35 | 28 |
| Rotational Speed, \(n\) (rpm) | 7,500 | 9,375 |
| Normal Module, \(m_n\) (mm) | 4 | 4 |
| Normal Pressure Angle, \(\alpha_n\) (°) | 20 | 20 |
| Helix Angle, \(\beta\) (°) | 25 | 25 |
| Face Width (per helix), \(B\) (mm) | 50 | 50 |
| Central Gap Width, \(b\) (mm) | 20 | 20 |
| Center Distance, \(a\) (mm) | 139 | 139 |
A hexahedral enclosure (control volume) is created to represent the fluid domain surrounding the herringbone gear pair. The dimensions of this domain are proportionally defined based on gear geometry to minimize boundary effects: Length \(L = 3d_1\), Width \(W = 1.2(2B+b)\), and Height \(H = 2d_1\), where \(d_1\) is the pitch diameter of the wheel. The meshing zone of the herringbone gear pair is positioned at the center of the domain. Supporting structures like shafts and bearings are omitted to isolate the windage effect of the gears themselves.
A critical challenge in transient CFD simulation of meshing gears is the extremely small clearance at the contact region, which can lead to negative cell volumes during mesh deformation. To circumvent this, the “tooth surface modification” method is employed. This technique involves slightly thinning the tooth profiles to create a small but consistent working clearance while strictly maintaining the original center distance, thus preserving the realistic meshing configuration without introducing geometric interference or significantly altering the tooth surface area—a factor crucial for accurate windage prediction.
1.2 Fluid Domain Discretization and Dynamic Mesh Technique
The fluid domain is discretized using an unstructured tetrahedral mesh, chosen for its adaptability to the complex gear geometry and compatibility with dynamic mesh remeshing algorithms. Surface meshes on the gear teeth are refined to resolve the high-velocity gradients and pressure fields near the rotating surfaces. The volume mesh is generated with a growth rate of 1.2, ensuring a minimum of 2-3 cell layers exist in the narrowed meshing gap. The final mesh consists of approximately 2.62 million surface cells and 2.8 million volume cells.
To simulate the rotation of the herringbone gears, Fluent’s dynamic mesh technology is employed. The enclosure walls remain stationary, while both herringbone gears are prescribed rigid-body rotation via User-Defined Functions (UDFs). The mesh update is handled by a combination of smoothing and local remeshing methods. The smoothing method adjusts the position of interior nodes based on the boundary displacement, while the remeshing method locally reconstructs cells that become overly distorted or violate size criteria, thus maintaining mesh quality throughout the simulation.
1.3 Governing Equations and Turbulence Modeling
The flow is assumed to be transient, three-dimensional, and turbulent. The fluid (air) is treated as incompressible, as the maximum Mach number (\(Ma = v/c\)) based on the pitch line velocity is well below 0.3. The governing equations are the Unsteady Reynolds-Averaged Navier-Stokes (URANS) equations:
Continuity Equation:
$$ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \vec{u}) = 0 $$
Momentum Equation:
$$ \frac{\partial (\rho \vec{u})}{\partial t} + \nabla \cdot (\rho \vec{u} \vec{u}) = -\nabla p + \nabla \cdot (\mu_{eff} \nabla \vec{u}) $$
where \(\rho\) is the fluid density, \(\vec{u}\) is the velocity vector, \(t\) is time, \(p\) is pressure, and \(\mu_{eff}\) is the effective viscosity (sum of molecular and turbulent viscosity).
For turbulence closure, the RNG \(k\)-\(\epsilon\) model is selected. This two-equation model is derived using Re-Normalization Group theory, which provides an analytically-derived differential formula for effective viscosity and accounts for the effects of smaller scales of motion. It offers improved accuracy for flows involving high strain rates, streamline curvature, and rapid transients—characteristics present in the flow around rotating herringbone gears. The transport equations for turbulent kinetic energy \(k\) and its dissipation rate \(\epsilon\) are:
$$ \frac{\partial (\rho k)}{\partial t} + \frac{\partial (\rho k u_i)}{\partial x_i} = \frac{\partial}{\partial x_j} \left( \alpha_k \mu_{eff} \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( \alpha_\epsilon \mu_{eff} \frac{\partial \epsilon}{\partial x_j} \right) + C_{1\epsilon} \frac{\epsilon}{k} G_k – C_{2\epsilon} \rho \frac{\epsilon^2}{k} – R_\epsilon $$
Here, \(G_k\) represents the generation of turbulent kinetic energy due to mean velocity gradients. The term \(R_\epsilon\) is a special addition of the RNG model, which improves accuracy for rapidly strained flows. The model constants are: \(C_{1\epsilon}=1.42\), \(C_{2\epsilon}=1.68\), \(\alpha_k = \alpha_\epsilon ≈ 1.393\).
A scalable wall function approach is used to model the near-wall flow, providing consistent results regardless of the near-wall mesh refinement (\(y^+\)), which is advantageous for complex deforming meshes where \(y^+\) values can vary.
2. Mechanism of Windage Loss in Herringbone Gears
Windage loss originates from the torque required to overcome the resistance offered by the fluid to the rotating gear. This resisting torque has two primary components: pressure drag torque and viscous friction torque.
As a herringbone gear rotates, it imparts momentum to the surrounding air. In an absolute frame of reference, the air near the gear circumference moves in the same direction as the gear’s tangential velocity, but at a lower speed. In a frame rotating with the gear, the air appears to flow opposite to the direction of rotation. This relative airflow separates at the tooth tip. One part impacts the leading face (windward side) of the downstream tooth, creating a region of high pressure, while the other part flows over the tip. The tooth’s trailing face (leeward side) typically experiences a region of lower pressure, often with recirculating flow in the tooth space. The pressure difference between the windward and leeward faces generates a net pressure torque that opposes the gear’s motion.
Simultaneously, viscous shear stresses act tangentially on all wetted surfaces—the tooth flanks, the two end faces, and the surfaces of the central gap. The component of this viscous force opposing the motion contributes to the viscous friction torque. For a herringbone gear, the total windage torque \(T_w\) can be decomposed based on the contributing surfaces:
$$ T_w = T_{wt} + T_{wf} + T_{wg} $$
where \(T_{wt}\), \(T_{wf}\), and \(T_{wg}\) are the torques on the tooth flanks, gear end faces, and central gap faces, respectively.
The instantaneous windage power loss \(P_w\) for a single gear is then:
$$ P_w = T_w \cdot \omega = (T_{wt} + T_{wf} + T_{wg}) \cdot \frac{2\pi n}{60} $$
For a herringbone gear pair, the total windage power loss \(P_{total}\) is the sum of the losses from both the wheel and the pinion:
$$ P_{total} = P_{w, wheel} + P_{w, pinion} $$
Simulation results for the baseline case confirm that the torque on the tooth flanks is orders of magnitude larger than the torque on the end faces and central gap. This indicates that pressure drag on the teeth is the dominant mechanism for windage loss in herringbone gears, accounting for over 96% of the total loss. The contribution from viscous friction on the relatively smooth end and gap surfaces is minor. The total windage power for the baseline herringbone gear pair stabilizes after an initial transient period, with an average value of approximately 1.33 kW. In a transmission system delivering 400 kW, this windage loss alone could constitute a significant portion of the total power loss budget, highlighting its importance in high-speed design.
3. Parametric Analysis of Influencing Factors
A systematic analysis was conducted to evaluate the impact of rotational speed, direction of rotation, and helix angle on the windage performance of the herringbone gear pair. The parameter variations are listed in Table 2.
| Case Description | Speed, \(n\) (rpm) | Helix Angle, \(\beta\) (°) | Normal Module, \(m_n\) (mm) |
|---|---|---|---|
| Baseline | 7,500 / 9,375 | 25 | 4.00 |
| Reverse Rotation | -7,500 / -9,375 | 25 | 4.00 |
| Speed Variation | 2,500 / 3,125 | 25 | 4.00 |
| 5,000 / 6,250 | 25 | 4.00 | |
| 10,000 / 12,500 | 25 | 4.00 | |
| 12,500 / 15,625 | 25 | 4.00 | |
| 15,000 / 18,750 | 25 | 4.00 | |
| Helix Angle Variation | 7,500 / 9,375 | 10 | 4.35 |
| 7,500 / 9,375 | 15 | 4.26 | |
| 7,500 / 9,375 | 20 | 4.15 | |
| 7,500 / 9,375 | 30 | 3.82 | |
| 7,500 / 9,375 | 35 | 3.62 | |
| 7,500 / 9,375 | 40 | 3.38 |
Note: For helix angle variations, the normal module was adjusted to maintain a constant center distance and transverse module, ensuring a fair comparison by keeping the gear volume and domain size relatively constant.
3.1 Effect of Rotational Speed
Rotational speed is the most significant factor affecting windage loss. The relationship between average windage power and speed for the herringbone gear pair is shown in Figure 1 (conceptual data). A power-law curve fit reveals a nearly cubic relationship:
$$ P_{total} \approx 2.774 \times 10^{-9} \cdot n^{3.0139} $$
This aligns well with established empirical formulas for spur gears, which typically suggest windage power is proportional to the cube of speed (\(P \propto n^3\)). The physical explanation is that the dynamic pressure (\(\frac{1}{2}\rho v^2\)) and the rate of kinetic energy imparted to the fluid both scale with the cube of the tangential velocity. As speed increases, the flow field becomes more turbulent and unsteady, leading to greater fluctuation in the instantaneous windage power signal. The proportion of tooth flank loss to total loss remains consistently high (above 96%) and increases slightly with speed, further underscoring the dominance of tooth-related pressure drag.
3.2 Effect of Direction of Rotation
Reversing the direction of rotation of the herringbone gear pair has a profound effect on the flow field structure but a negligible effect on the magnitude of windage power loss. In the baseline (forward) rotation, the central gap acts as the primary air intake region; air is drawn axially inwards from the gap and expelled radially and axially from both outer ends of the gears, creating two symmetric vortex structures above the gear pair. In reverse rotation, this pattern is inverted: air is drawn in from both outer ends and expelled through the central gap. This phenomenon is due to the axial pumping action induced by the helical teeth, analogous to a fan. For a given helix hand, reversing rotation reverses the direction of the axial pressure gradient and thus the flow direction.
Consequently, the location of the maximum pressure difference on a tooth flank shifts. In forward rotation, the highest drag occurs on teeth near the intake region (central gap), where the airflow is strongest. In reverse rotation, the highest drag shifts to the regions near the end faces, which become the intakes. However, the integrated effect over all teeth—the total windage torque—remains statistically unchanged, as the flow physics are symmetric with respect to rotation direction. The windage power curves for forward and reverse rotation are virtually identical.
3.3 Effect of Helix Angle
The helix angle is a critical design parameter for herringbone gears, influencing load distribution, noise, and axial force balance. Its effect on windage loss is significant and stems from two main factors: the change in the projected frontal area of the teeth and the modification of the axial flow component. Increasing the helix angle reduces the normal module (to maintain center distance with fixed teeth numbers), which slightly reduces tooth height. More importantly, it increases the tooth angle relative to the axis of rotation, effectively reducing the area projected perpendicular to the direction of the primary circumferential flow. This reduction in “frontal area” leads to a decrease in pressure drag.
Simulation results confirm that increasing the helix angle reduces the total windage power. The reduction is more pronounced at higher helix angles. For instance, increasing the helix angle from 25° to 40° leads to a more significant reduction in loss compared to an increase from 10° to 25°. This trend highlights a potential efficiency benefit in using larger helix angles for high-speed herringbone gears, although this must be balanced against other design constraints such as bending strength and manufacturing considerations.
4. Transient Torque Characteristics on a Single Tooth
Analyzing the instantaneous windage torque on a single, monitored tooth provides deeper insight into the interaction between gear meshing and the flow field. The torque on a single tooth of the wheel or pinion exhibits distinct periodic behavior, synchronized with its rotation period.
When the tooth is far from the meshing zone, its torque remains relatively stable, corresponding to its interaction with a largely uniform circumferential flow. However, as the tooth approaches and passes through the meshing region, its torque undergoes a characteristic fluctuation sequence: Increase → Decrease to a negative value → Second increase → Decrease back to a stable level.
This pattern is directly linked to the localized pressure field in the meshing zone. As the gear pair rotates, the meshing teeth compress the air in the ingoing (approach) side, creating a localized high-pressure region. Conversely, on the outgoing (recess) side, the separating teeth create a low-pressure region. As a monitored tooth enters the meshing zone:
- Its leading face encounters the high-pressure region, causing a sharp increase in positive (resisting) torque.
- As it moves through the mesh, its leading face may enter the low-pressure zone while its trailing face is still in the high-pressure zone, causing the net pressure torque to reverse direction, resulting in a brief negative torque value.
- Subsequently, as the tooth exits, its trailing face experiences the low-pressure region, leading to another increase in positive torque.
- Finally, away from the meshing zone, the torque stabilizes again.
The amplitude of this torque fluctuation is strongly influenced by operational and geometric parameters:
- Speed: Higher rotational speeds increase the magnitude of the pressure zones in the mesh, leading to larger fluctuation amplitudes in single-tooth torque.
- Helix Angle: Larger helix angles tend to reduce the fluctuation amplitude. This is likely because the more gradual engagement of helical teeth smears the pressure concentration over a longer contact line and time, reducing the peak pressure intensities compared to a spur-like engagement.
It is important to note that while individual tooth torques fluctuate significantly, the total torque on the entire gear is the sum of the torques on all its teeth. Since the phases of these fluctuations are distributed around the gear’s circumference, they average out, resulting in the relatively stable total torque signal observed for the complete herringbone gear.
5. Conclusion
This study presents a comprehensive transient CFD analysis of windage losses in high-speed herringbone gear pairs. Utilizing the RNG \(k\)-\(\epsilon\) turbulence model coupled with dynamic mesh techniques provides a detailed view of the complex, three-dimensional, and unsteady airflow generated by rotating herringbone gears. The key conclusions are as follows:
- Dominant Loss Mechanism: Windage loss in herringbone gears is predominantly caused by pressure drag on the tooth flanks, which accounts for over 96% of the total loss. Viscous friction on the end faces and central gap surfaces contributes minimally.
- Speed Dependency: The total windage power loss of a herringbone gear pair exhibits a strong correlation with rotational speed, following an approximately cubic relationship (\(P \propto n^{3.01}\)). This aligns with established empirical knowledge for spur gears and underscores the dramatic increase in aerodynamic loss with speed.
- Effect of Rotation Direction: Reversing the direction of rotation fundamentally alters the flow pattern around the herringbone gears, swapping the primary intake and exhaust regions. However, this change in flow topology does not measurably affect the total magnitude of windage power loss; the system exhibits symmetry with respect to rotation direction.
- Beneficial Role of Helix Angle: Increasing the helix angle is an effective geometrical strategy for reducing windage loss in herringbone gears. The reduction becomes more significant at larger helix angles (e.g., above 25°), primarily due to the decrease in the projected frontal area of the teeth against the circumferential flow.
- Single-Tooth Torque Dynamics: The instantaneous windage torque on a single tooth undergoes a characteristic cyclic fluctuation as it passes through the meshing zone, driven by the localized high and low-pressure fields generated by tooth engagement and disengagement. The amplitude of these fluctuations increases with speed and decreases with increasing helix angle.
The findings from this analysis provide valuable guidelines for the design and optimization of herringbone gear transmissions intended for high-speed operation. Designers should pay particular attention to minimizing tooth-related drag, possibly through aerodynamic shaping or shroud design, and can consider the beneficial effect of a larger helix angle within the constraints of strength and dynamics. The methodologies and results established here form a foundation for further research into more complex scenarios, such as windage loss in oil-mist environments or in multi-stage herringbone gear systems.