In the field of aviation fuel systems, the helical gear pump stands out due to its smooth operation, reduced noise, and high load-bearing capacity. However, the complex geometry of helical gears makes it challenging to analyze the instantaneous geometric flow characteristics. In this article, I will explore a calculation method for the instantaneous geometric flow rate of an external meshing involute helical gear pump used in aviation fuel applications. The goal is to clarify how the helical structure influences flow dynamics, leveraging parametric modeling and mathematical formulations. The method involves discretizing the pump’s sealing discharge chamber into axial slices, modeling boundary curves, and applying Green’s theorem for area calculation. By summing volumes across slices, the total discharge chamber volume variation with gear rotation is obtained, leading to instantaneous flow characteristics. This approach provides insights into design parameters like helical angle and tooth width, impacting flow pulsation and average output.
The helical gear pump operates on the principle of varying sealed volumes to facilitate fluid transfer. Unlike spur gears, helical gears introduce a phase lag along the axial direction due to their twist, complicating the密封容腔 geometry. Traditional methods, such as the swept-area approach or energy conservation, have yielded mixed results on the effect of helical gears on flow pulsation. Here, I propose a method that treats the discharge chamber as a stack of thin axial slices. Each slice’s cross-section is a single-connected domain bounded by piecewise smooth curves derived from gear tooth profiles. Using parameterized equations for involute curves, transition curves, and circles, the boundary is mathematically defined. Green’s formula converts area integration into a line integral, simplifying computation. By updating the gear rotation angle, the volume change per time step gives the instantaneous geometric flow rate. This method is validated against CAD models and computational fluid dynamics (CFD) simulations, ensuring accuracy.
The axial slicing principle is central to this method. Consider a helical gear pump with left-handed driving and right-handed driven gears. At a reference starting point, one pair of teeth is in full-width meshing, isolating the discharge chamber. As the gear rotates, subsequent teeth engage gradually along the axis due to the helical angle. The discharge chamber is divided into n axial slices, each with a thickness of B/n, where B is the tooth width. For each slice, the cross-sectional area is computed via Green’s formula. The total volume is the sum of slice volumes. The helical nature causes a phase lag between slices, calculated as: $$ \theta_z = \frac{h \tan \beta}{r} $$ where h is the axial distance from the reference slice, β is the helical angle, and r is the pitch circle radius. This lag affects the boundary curves in each slice, requiring coordinate transformations.
To model the boundary curves, parametric equations for tooth profiles are essential. For an involute helical gear, the tooth profile on any axial slice can be derived from the base profile adjusted by the phase lag. The involute curve in a coordinate system attached to the gear is given by: $$ x_k = r_b \sin \mu_k – r_b \mu_k \cos \mu_k, $$ $$ y_k = r_b \cos \mu_k + r_b \mu_k \sin \mu_k, $$ where r_b is the base circle radius and μ_k is the roll angle, related to the pressure angle α_k by μ_k = \tan α_k. A parameter t_1 (0 ≤ t_1 ≤ 1) is introduced to simplify calculations: α_k = t_1 α_{km}, with α_{km} being the maximum pressure angle at the tip circle. For a helical gear, the coordinates transform to account for axial position. If (x_1, y_1) are coordinates in the front-face system, then for a slice at axial distance h, the coordinates (x, y) in a fixed global system are: $$ \begin{bmatrix} x \\ y \\ z \end{bmatrix} = \begin{bmatrix} \cos \theta_z & -\sin \theta_z & 0 \\ \sin \theta_z & \cos \theta_z & 0 \\ 0 & 0 & 1 \end{bmatrix} \begin{bmatrix} x_1 \\ y_1 \\ 0 \end{bmatrix} + \begin{bmatrix} 0 \\ r_w \\ h \end{bmatrix}, $$ where r_w is the pitch circle radius. Similar transformations apply to the driven gear.
The root transition curve, often generated by a hob cutter, is modeled with additional parameters. For a double-arc hob, the transition curve coordinates (x_g, y_g) are: $$ x_g = r \sin \phi – \left( e – \frac{\zeta m}{\sin \gamma} + \rho \right) \cos(\gamma – \phi), $$ $$ y_g = r \cos \phi – \left( e – \frac{\zeta m}{\sin \gamma} + \rho \right) \sin(\gamma – \phi), $$ where r is the pitch radius, φ and e are tool-related parameters, ρ is the tool arc radius, ζ is the modification coefficient, m is the module, and γ varies from the transverse pressure angle to 90°. These coordinates are then rotated to align with the gear’s orientation.
The tip and root circles are simple circles: for the driving gear, the tip circle equation is x_1^2 + (y_1 – r_w)^2 = r_a^2, and the root circle is x_1^2 + (y_1 – r_w)^2 = r_f^2, where r_a and r_f are the tip and root radii. For the driven gear, similar equations apply with adjusted signs. The meshing line, which is the common tangent to the base circles, has the equation: $$ y = \frac{\sqrt{a^2 – (2r_b)^2}}{2r_b} x, $$ where a is the center distance. Meshing points, crucial for boundary definition, are found by solving for intersections between the meshing line and tooth profiles. For a slice at axial distance h, the meshing point coordinates depend on the gear rotation angle θ and the phase lag θ_z. The distance f from the pitch point to the meshing point is: $$ f = r_b (\omega t – \theta_z), $$ with ω as the angular velocity and t as time. The parameters t_1 and t_2 for the driving and driven gears are computed using inverse trigonometric functions based on distances from gear centers.
The condition to avoid oil churning and ensure reliable separation of suction and discharge chambers in a helical gear pump is: $$ \frac{z(\tan \alpha_{at} – \tan \alpha’)}{\pi} – \frac{B \sin \beta}{\pi m_n} \geq 1, $$ where z is the number of teeth, α_{at} is the tip circle pressure angle, α’ is the operating pressure angle, m_n is the normal module, and B is the tooth width. When equality holds, the pump has no churning or leakage; otherwise, relief grooves are needed. The angle for double-tooth meshing is: $$ \beta_c = (\varepsilon_\alpha – 1) \frac{2\pi}{z} = 2(\tan \alpha_{at} – \tan \alpha’) – \frac{2\pi}{z}, $$ where ε_α is the transverse contact ratio. Relief grooves are typically placed at the axial mid-plane where the churning volume is minimal, corresponding to a gear rotation angle of: $$ \theta_f = \frac{\beta_c + \beta_0}{2}, $$ with β_0 as the lag angle between front and rear faces.
The calculation process involves dividing the gear rotation period into intervals based on meshing states. For a rotation angle θ less than β_0, only one tooth pair is in full-width meshing; the discharge chamber boundary includes the meshing point of that pair. When β_0 ≤ θ < θ_f, two tooth pairs are in full-width meshing, creating a churning zone connected to the discharge chamber via relief grooves. For θ_f ≤ θ ≤ 2π/z, the churning zone connects to the suction chamber, and the discharge chamber boundary shifts to the later meshing pair. This segmentation ensures accurate volume computation.
Using Green’s formula, the area A_i of the discharge chamber cross-section for the i-th slice is: $$ A_i = \frac{1}{2} \oint_L (x \, dy – y \, dx), $$ where L is the closed boundary curve. The total volume V is: $$ V = \frac{B}{n} \sum_{i=1}^n A_i. $$ The instantaneous geometric flow rate q_V is then: $$ q_V = \frac{V(t + \Delta t) – V(t)}{\Delta t}. $$ By iterating over small time steps, the flow rate curve over a gear rotation cycle is generated. This method accounts for the helical gear’s axial variations, providing a precise tool for design analysis.
To validate the method, I applied it to a case study with parameters typical of aviation fuel helical gear pumps. The gear specifications are summarized in Table 1.
| Parameter | Symbol | Value |
|---|---|---|
| Number of Teeth | z | 10 |
| Normal Module | m_n | 3 mm |
| Helical Angle | β | 11° |
| Tooth Width | B | 12.5 mm |
| Normal Pressure Angle | α_n | 28° |
| Center Distance | a | 30 mm |
| Modification Coefficient | ζ | -0.1098 |
| Tip Radius | r_a | 16.5 mm |
| Root Radius | r_f | 13.2 mm |
The helical gear pump was sliced into n = 20 layers for computation. At a reference rotation angle of 0°, the discharge chamber volume calculated was 2539.971 mm³, matching CAD measurements within 0.0069% error. The volume variation with gear rotation is plotted, showing a discontinuity near 9° due to churning volume changes. The instantaneous flow rate curve was derived and compared with existing formulas. For instance, one common formula for helical gear pump flow is: $$ q_V = \sum_{i=1}^n \frac{B \omega}{n} \left[ r_a^2 – r_w^2 – r_b^2 \left( \theta – \frac{i B \tan \beta}{n r} \right) \right]. $$ Another simplified version is: $$ q_V = B \omega (r_a^2 – r_w^2 – u_0^2) – \frac{\omega \cos^2 \alpha_n \tan^2 \beta}{12} B^3, $$ where u_0 is the distance from the pitch point to the meshing point at the mid-plane. My method yielded an average flow rate of 17.4602 L/min, closely aligning with the simplified formula’s result of 17.4701 L/min, confirming consistency.
The influence of the helical angle on flow characteristics was investigated by varying β while keeping tooth width constant. Table 2 summarizes the results for pumps with relief grooves, showing that as the helical angle increases, the average flow rate slightly decreases and the flow pulsation rate marginally increases. This is attributed to the phase lag reducing the effective discharge volume per rotation. For pumps without relief grooves, increasing the helical angle reduces the churning angle interval and volume, thereby lowering pulsation. This highlights the dual role of helical gears in design trade-offs.
| Helical Angle β (°) | Average Flow Rate q_{Vm} (L/min) | Flow Pulsation Rate δ (%) |
|---|---|---|
| 3 | 17.5354 | 19.307 |
| 5 | 17.5274 | 19.314 |
| 7 | 17.5143 | 19.323 |
| 9 | 17.4961 | 19.336 |
| 11 | 17.4602 | 19.359 |
The tooth width also plays a significant role. With a fixed helical angle of 11°, varying B affects the flow output. As shown in Table 3, for pumps with relief grooves, increasing tooth width raises the average flow rate but slightly increases pulsation. Without relief grooves, larger tooth widths reduce pulsation by expanding the discharge interval. This underscores the importance of holistic design when optimizing helical gear pumps for minimal flow fluctuation.
| Tooth Width B (mm) | Relief Grooves | Average Flow Rate q_{Vm} (L/min) | Flow Pulsation Rate δ (%) |
|---|---|---|---|
| 5 | Yes | 7.0125 | 19.310 |
| 5 | No | 6.9809 | 25.247 |
| 9 | Yes | 12.6039 | 19.329 |
| 9 | No | 12.5778 | 22.053 |
| 12.5 | Yes | 17.4602 | 19.359 |
| 12.5 | No | 17.3463 | 29.849 |
To further validate the method, computational fluid dynamics (CFD) simulations were conducted for selected helical gear configurations. The CFD model accounted for fluid properties of aviation fuel RP-3, with an outlet pressure of 2 MPa and consideration of cavitation using the Singhal model. The mesh was refined to ensure independence, with gear side clearances set to 10 µm to minimize leakage effects. The instantaneous flow curves from CFD showed good agreement with the proposed method, though with slightly lower pulsation amplitudes due to fluid dynamics and leakage. For example, for a helical gear pump with β = 11° and B = 12.5 mm, the CFD average flow was 16.2417 L/min, about 1.18% lower than the geometric calculation, but the pulsation frequency and trend matched well.
The CFD results reinforced the parametric trends: increasing helical angle reduces average flow and increases pulsation for pumps with relief grooves, while increasing tooth width boosts flow but can alter pulsation depending on relief design. This alignment between geometric and CFD analyses confirms the robustness of the slicing method for helical gear pump design. The helical gear’s ability to smooth engagement is evident, but its impact on flow must be balanced with other factors like efficiency and noise.
In conclusion, the axial slicing method combined with Green’s theorem offers a precise way to compute instantaneous geometric flow rates for helical gear pumps. The helical gear structure introduces complexities that are effectively captured through phase-lagged boundary modeling. Key findings include that for pumps with relief grooves, larger helical angles slightly reduce average flow and increase pulsation, while larger tooth widths raise flow but may heighten pulsation. For pumps without relief grooves, both increasing helical angle and tooth width can reduce pulsation by mitigating churning. These insights aid in optimizing helical gear pump designs for aviation fuel systems, where low pulsation and reliability are critical. Future work could extend this method to variable operating conditions or incorporate more advanced fluid-structure interactions for even greater accuracy in predicting helical gear pump performance.