High-pressure hydraulic pumps utilizing external helical gears are fundamental components in power transmission systems. Their advantages, including simple construction, compact size, lightweight design, and high reliability, are complemented by a robust resistance to contamination, enabling operation in demanding environments. This combination makes them one of the most reliable and widely used hydraulic components. However, as operational pressure and rotational speed increase, significant challenges arise, such as cavitation, increased vibration and noise, and heightened internal leakage. These issues critically compromise the operational stability and service life of the gear pump. Consequently, a thorough investigation into the internal flow field characteristics of high-pressure pumps is essential for performance optimization.
Previous research has extensively employed computational fluid dynamics (CFD) to analyze gear pump behavior. Studies have utilized dynamic mesh techniques to simulate transient fluid behavior, analyzing forces on gears and the phenomenon of fluid trapping by varying boundary conditions like inlet pressure and rotational speed. Other investigations have focused on the hydraulic performance of low-speed pump impellers, examining the influence of blade geometry. However, a comprehensive analysis integrating the effects of critical design parameters like radial clearance and operational speed on the flow field, particularly for helical gear pumps, remains an area for deeper exploration. The flow pulsation coefficient is a key metric for evaluating pump performance. This study addresses gaps in existing models by investigating the influence of the helix angle, radial clearance, and rotational speed on both pressure pulsation and the flow pulsation coefficient through combined theoretical analysis and CFD simulation, providing valuable insights for designing high-performance helical gear pumps.
Theoretical Analysis of Flow Pulsation
The instantaneous flow rate is a critical parameter for understanding the dynamic behavior of a helical gear pump. The model can be derived by discretizing the helical gear into n segments along its axis, treating each segment as a spur gear pump. The total instantaneous flow rate is the sum of contributions from all segments.
The instantaneous flow rate for an external meshing high-pressure helical gear pump can be expressed as:
$$q_v = \sum_{i=0}^{n-1} \frac{\omega_1 \left[ R_a^2 – R^2 – f^2 \left( \theta + \frac{i b \tan \beta}{n R} \right)^2 \right] b}{n}$$
where \( \omega_1 \) is the angular velocity of the driving gear, \( R_a \) is the tip circle radius, \( R \) is the pitch circle radius, \( f \) is the distance from the meshing point to the pitch point, \( \theta \) is the rotation angle corresponding to the meshing point on the reference plane, \( b \) is the gear face width, \( \beta \) is the helix angle, and \( n \) is the number of axial segments.
The flow pulsation coefficient (\( \phi \)) quantifies the uniformity of the discharge flow. It is defined as the ratio of the difference between maximum and minimum instantaneous flow to the average flow rate:
$$\phi = \frac{q_{v_{\text{max}}} – q_{v_{\text{min}}}}{\bar{q}_v}$$
where \( \bar{q}_v = (q_{v_{\text{max}}} + q_{v_{\text{min}}}) / 2 \). Considering the gear meshing process where the contact point moves along the line of action, the parameter \( f \) varies within the range defined by the contact ratio \( \varepsilon \) and the base pitch \( r_j \): \( f \in (-\varepsilon r_j / 2, \varepsilon r_j / 2) \).
The minimum instantaneous flow occurs when a tooth pair either begins or ends meshing (\( f = \pm \varepsilon r_j / 2 \)):
$$q_{v_{\text{min}}} = \sum_{i=0}^{n-1} \frac{\omega_1 \left[ R_a^2 – R^2 – \frac{\varepsilon^2 r_j^2}{4} \left( \theta + \frac{i b \tan \beta}{n R} \right)^2 \right] b}{n}$$
The maximum instantaneous flow occurs when the meshing point is at the pitch point (\( f = 0 \)):
$$q_{v_{\text{max}}} = \sum_{i=0}^{n-1} \frac{\omega_1 \left( R_a^2 – R^2 \right) b}{n}$$
The theoretical flow rate, accounting for the volumetric displacement, is given by:
$$q_{v_t} = \sum_{i=0}^{n-1} \frac{\omega_1 \left[ R_a^2 – R^2 – \frac{K_v r_j^2}{12} \left( \theta + \frac{i b \tan \beta}{n R} \right)^2 \right] b}{n}$$
where \( K_v = 3\varepsilon^2 – 6\varepsilon + 4 \). Substituting these expressions into the definition of the pulsation coefficient yields a comprehensive formula that highlights the influencing parameters:
$$\phi = \frac{\varepsilon^{2} r_{j}^{2} \left( \theta + \frac{i b \tan \beta}{n R} \right)^{2}}{2 \left[ R_{a}^{2} – R^{2} – \frac{K_{v} r_{j}^{2}}{12} \left( \theta + \frac{i b \tan \beta}{n R} \right)^{2} \right]}$$
Here, \( r_j = m_n \pi \cos \alpha_n \), where \( m_n \) is the normal module and \( \alpha_n \) is the normal pressure angle. The tip radius is \( R_a = m_t z / 2 + h_a \), where \( m_t \) is the transverse module, \( z \) is the number of teeth, and \( h_a \) is the addendum.
Analysis of this formula reveals a direct relationship between the helix angle \( \beta \) and the flow pulsation coefficient \( \phi \). To avoid generating excessive axial thrust while maintaining the benefits of helical teeth, the helix angle is typically designed within a range of 8° to 15°. The calculated variation of \( \phi \) with \( \beta \) within this range is summarized in the table below.
| Helix Angle, β (degrees) | Flow Pulsation Coefficient, φ | Reduction Rate vs. 8° (%) |
|---|---|---|
| 8 | 0.67 | 0.0 |
| 9 | 0.66 | 1.5 |
| 10 | 0.64 | 4.5 |
| 11 | 0.62 | 7.5 |
| 12 | 0.60 | 10.4 |
| 13 | 0.57 | 14.9 |
| 14 | 0.55 | 17.9 |
| 15 | 0.52 | 22.4 |
The data clearly indicates that increasing the helix angle effectively reduces the flow pulsation coefficient, thereby improving the quality and smoothness of the outlet flow. On average, each degree increase in the helix angle within this practical range reduces φ by approximately 2%. This provides a clear theoretical guideline for designers: specifying a larger helix angle, within mechanical limits, is beneficial for mitigating flow pulsation in a helical gear pump.
Numerical Simulation Methodology
The study focuses on a specific external helical gear pump model. The primary geometrical and operational parameters are defined as follows:
| Parameter | Symbol | Value |
|---|---|---|
| Number of Teeth | z | 18 |
| Profile Shift Coefficient | X | -0.02636 |
| Normal Module | m_n | 2.383 mm |
| Helix Angle | β | 8° 35′ |
| Normal Pressure Angle | α_n | 20° |
| Center Distance | a | 43 mm |
| Fluid Density | ρ | 900 kg/m³ |
| Fluid Dynamic Viscosity | μ | 0.026 Pa·s |
| Rated Rotational Speed | n_rated | 2000 rpm |
| Nominal Radial Clearance | δ_rated | 1.0 mm |
To develop a comprehensive understanding of the internal flow field under various conditions, a 3×3 simulation matrix was designed, varying rotational speed and radial clearance around the nominal design point.
| Case ID | Rotational Speed, n (rpm) | Radial Clearance, δ (mm) | Turbulence Model |
|---|---|---|---|
| 1 | 1400 | 0.8 | k-ε |
| 2 | 2000 | 1.0 | |
| 3 | 2200 | 1.2 | |
| 4 | 1400 | 1.0 | k-ε |
| 5 (Nominal) | 2000 | 1.0 | |
| 6 | 2200 | 1.0 | |
| 7 | 1400 | 1.2 | k-ε |
| 8 | 2000 | 1.2 | |
| 9 | 2200 | 1.2 |
CFD Model and Governing Equations
The transient, three-dimensional, incompressible flow within the helical gear pump is governed by the Reynolds-Averaged Navier-Stokes (RANS) equations. The continuity and momentum equations are:
$$\nabla \cdot \vec{v} = 0$$
$$\frac{\partial \vec{v}}{\partial t} + (\vec{v} \cdot \nabla) \vec{v} = -\frac{1}{\rho} \nabla p + \nu \nabla^2 \vec{v}$$
where \( \vec{v} \) is the velocity vector, \( p \) is pressure, \( \rho \) is density, and \( \nu \) is the kinematic viscosity. The standard k-ε turbulence model was selected for its robustness and computational efficiency in simulating internal flows with recirculation and shear. To accurately capture the complex fluid motion driven by the rotating gears, a dynamic mesh approach with user-defined functions (UDFs) was employed to control the rigid body rotation of the gear domains. The computational domain included extended inlet and outlet sections to allow for fully developed flow and mitigate the influence of boundary conditions on the core pump chamber flow.
Mesh Generation and Boundary Conditions
Due to the complex geometry involving small clearances and the meshing region, an unstructured tetrahedral mesh was generated for the entire fluid domain. Particular attention was paid to the radial clearance and the inter-tooth gaps, where three layers of refined mesh were implemented to resolve the high-velocity gradient and leakage flows accurately. The total mesh count was approximately 10 million cells, ensuring grid-independent results for the parameters of interest. A rigorous time-step study was conducted, linking the time step size to the rotational speed to ensure accurate temporal resolution of the gear motion and flow phenomena.
The boundary conditions were set as follows: The inlet was defined as a pressure-inlet with a gauge pressure of 0.6 MPa. The outlet was defined as a pressure-outlet with a gauge pressure of 21 MPa, creating the high-pressure differential. All solid walls (gear teeth, pump casing) were treated as no-slip, adiabatic walls. The interfaces between the rotating gear zones and the stationary casing zone were handled using dynamic mesh sliding interfaces.
Results and Discussion
Pressure Distribution Characteristics
The pressure field within the helical gear pump is highly dynamic and non-uniform. Analysis of the nominal case (n=2000 rpm, δ=1.0 mm) reveals that the pressure in the meshing zone undergoes significant fluctuations during a rotation cycle. The pressure is lowest in the volume-expanding suction chamber and highest in the volume-contracting discharge chamber, with a steep gradient established across the meshing line that acts as a sealing barrier.
The effect of operational parameters on pressure was systematically investigated. For a fixed radial clearance, increasing the rotational speed leads to more pronounced pressure variations in the immediate vicinity of the meshing point. However, the pressure distribution on the gear teeth surfaces near the pump casing walls remains relatively less affected by speed changes, maintaining a more consistent profile. Conversely, at a fixed rotational speed, increasing the radial clearance results in a general reduction of the pressure magnitude within the pump chambers. This is attributed to the increased leakage path, which allows more fluid to bypass the sealing action of the meshing gears, thereby reducing the pressure build-up on the discharge side. This relationship can be conceptually summarized: the load borne by the gears and casing is inversely proportional to the radial clearance at a given speed.
Flow Pulsation Coefficient from CFD
The flow pulsation coefficient was extracted from the CFD results by monitoring the area-averaged velocity at the pump outlet over several complete rotation cycles. The maximum (\(v_{\text{max}}\)) and minimum (\(v_{\text{min}}\)) velocities were identified, and the pulsation coefficient was calculated using the relationship derived from the flow rate definition:
$$\phi = \frac{q_{v_{\text{max}}} – q_{v_{\text{min}}}}{\bar{q}_v} = \frac{2(v_{\text{max}} – v_{\text{min}})}{v_{\text{max}} + v_{\text{min}}}$$
For the nominal case, \(v_{\text{max}} = 3.54 \, \text{m/s}\) and \(v_{\text{min}} = 2.23 \, \text{m/s}\), yielding \( \phi \approx 0.12 \). The calculated pulsation coefficients for all cases are presented in the following analysis.
Effect of Rotational Speed: Holding the radial constant at δ=1.0 mm, the pulsation coefficient shows a distinct trend with speed. As the rotational speed increases from 1400 rpm to 2200 rpm, φ generally decreases. The most significant reduction occurs between 1700 rpm and 1800 rpm. For speeds exceeding 1900 rpm, the pulsation coefficient stabilizes, fluctuating around a low value of approximately 0.09. This indicates that operating the helical gear pump at higher speeds within its design range can effectively improve the smoothness of the output flow.
Effect of Radial Clearance: Holding the rotational speed constant, the influence of radial clearance on flow pulsation was analyzed. The results demonstrate that for a given speed, increasing the radial clearance leads to a reduction in the flow pulsation coefficient. This trend continues until the clearance reaches a threshold (around 1.0 mm for this model), beyond which the pulsation coefficient stabilizes and exhibits only minor fluctuations. This finding suggests that a larger radial clearance, within practical mechanical limits, can contribute to better flow quality by dampening the pulsations. However, this benefit must be carefully balanced against the associated increase in volumetric leakage.
Leakage Flow and Vortex Dynamics
The flow in the radial clearance is a critical aspect governing pump efficiency and pressure build-up. The pressure difference between the discharge (pressure side) and suction side of a gear tooth drives a high-velocity jet through the narrow clearance between the tooth tip and the pump housing. This leakage flow separates from the tooth surface as it enters the clearance, reattaches on the opposite side, and interacts with the main chamber flow to form complex vortex structures known as leakage vortices within the inter-tooth spaces on the suction side.
CFD results clearly captured these leakage vortices. Their strength and structure are sensitive to both rotational speed and radial clearance. For a constant radial clearance, an increase in rotational speed was found to reduce the intensity of the leakage vortex. While this is beneficial for reducing local losses, the concomitant larger pressure swings in the meshing zone at higher speeds elevate the risk of cavitation. For a constant rotational speed, increasing the radial clearance also leads to a reduction in leakage vortex intensity. This is because the wider gap offers less flow resistance, reducing the jet velocity and its shear interaction with the cavity flow. Analysis of the velocity vectors confirmed that the leakage velocity is highest at the sides of the tooth tip, not directly at the tip’s center. The trade-off is clear: while a larger clearance reduces vortex strength and potentially flow pulsation, it directly increases the volumetric leakage rate, reducing the pump’s volumetric efficiency. Furthermore, although a smaller clearance minimizes leakage, it increases viscous friction losses and the associated risk of wear and power loss.
Conclusion
This integrated study, combining theoretical modeling and high-fidelity CFD simulation, provides significant insights into the performance characteristics of external helical gear high-pressure pumps. The key conclusions are as follows:
1. Theoretical analysis establishes that the flow pulsation coefficient (φ) is inversely proportional to the helix angle (β). Increasing β from 8° to 15° can reduce φ by over 20%, with an average reduction of about 2% per degree. This provides a direct and effective design lever for improving flow smoothness in a helical gear pump.
2. CFD simulations validate the design approach and reveal the complex interplay between operational parameters. For the studied pump, operating at the nominal speed of 2000 rpm with a radial clearance of 1.0 mm results in a low and stable flow pulsation coefficient (φ ≈ 0.12) and manageable leakage vortex intensity. The results confirm that moderate increases in both rotational speed and radial clearance can synergistically improve outlet flow quality.
3. A detailed analysis of the internal flow field elucidates the transient pressure distribution and leakage mechanisms. Pressure fluctuates most significantly in the meshing region, while gear loads are inversely related to radial clearance. Leakage vortices formed in the tooth cavities are weakened by both higher speeds and larger clearances, although these changes come with trade-offs regarding cavitation risk and volumetric efficiency, respectively.
4. The employed CFD methodology serves as a powerful virtual prototyping tool. It not only allows for performance evaluation and validation of a helical gear pump design but also enables deep analysis of intricate flow details—such as leakage paths and vortex dynamics—that are difficult to assess experimentally. This capability is invaluable for guiding hydraulic design improvements, enhancing overall pump efficiency and reliability, and significantly shortening the product development cycle for advanced high-pressure helical gear pumps.