Pump Grooved Seals: A Computational Fluid Dynamics Approach to Improve Bulk-Flow Model Predictions

2019 ◽  
Vol 141 (10) ◽  
Author(s):  
Tingcheng Wu ◽  
Luis San Andrés

In multiple stage centrifugal pumps, balance pistons, often comprising a grooved annular seal, equilibrate the full pressure rise across the pump. Grooves in the stator break the evolution of fluid swirl and increase mechanical energy dissipation; hence, a grooved seal offers a lesser leakage and lower cross-coupled stiffness than a similar size uniform clearance seal. To date, bulk-flow modelbulk-flow models (BFMs) expediently predict leakage and rotor dynamic force coefficients of grooved seals; however, they lack accuracy for any other geometry besides rectangular. Note that scalloped and triangular (serrated) groove seals are not uncommon. In these cases, computational fluid dynamics (CFD) models seals of complex shape to produce leakage and force coefficients. Alas, CFD is not yet ready for routine engineer practice. Hence, an intermediate procedure presently takes an accurate two-dimensional (2D) CFD model of a smaller flow region, namely a single groove and adjacent land, to produce stator and rotor surface wall friction factors, expressed as functions of the Reynolds numbers, for integration into an existing BFM and ready prediction of seal leakage and force coefficients. The selected groove-land section is well within the seal length and far away from the effects of the inlet condition. The analysis takes three water lubricated seals with distinct groove shapes: rectangular, scalloped, and triangular. Each seal, with length/diameter L/D = 0.4, has 44 grooves of shallow depth dg ∼ clearance Cr and operates at a rotor speed equal to 5,588 rpm (78 m/s surface speed) and with a pressure drop of 14.9 MPa. The method validity is asserted when 2D (single groove-land) and three-dimensional (3D) (whole seal) predictions for pressure and velocity fields are compared against each other. The CFD predictions, 2D and 3D, show that the triangular groove seal has the largest leakage, 41% greater than the rectangular groove seal does, albeit producing the smallest cross-coupled stiffnesses and whirl frequency ratio (WFR). On the other hand, the triangular groove seal has the largest direct stiffness and damping coefficients. The scalloped groove seal shows similar rotordynamic force coefficients as the rectangular groove seal but leaks 13% more. For the three seal groove types, the modified BFM predicts leakage that is less than 6% away from that delivered by CFD, whereas the seal stiffnesses (both direct and cross-coupled) differ by 13%, the direct damping coefficients by 18%, and the added mass coefficients are within 30%. The procedure introduced extends the applicability of a BFM to predict the dynamic performance of grooved seals with distinctive shapes.

Author(s):  
Tingcheng Wu ◽  
Luis San Andrés

Abstract In multiple stage centrifugal pumps, balance pistons, often comprising a grooved annular seal, equilibrate the full pressure rise across the pump. Grooves in the stator break the evolution of fluid swirl and increase mechanical energy dissipation; hence, a grooved seal offers a lesser leakage and lower cross-coupled stiffness than a similar size uniform clearance seal. To date bulk-flow models (BFMs) expediently predict leakage and rotor dynamic force coefficients of grooved seals; however, they lack accuracy for any other geometry besides rectangular. Note scalloped and triangular (serrated) groove seals are not uncommon. In these cases, computational fluid dynamics (CFD) models seals of complex shape to produce leakage and force coefficients. Alas CFD is not yet ready for routine engineer practice. Hence, an intermediate procedure presently takes an accurate two-dimensional (2D) CFD model of a smaller flow region, namely a single groove and adjacent land, to produce stator and rotor surface wall friction factors, expressed as functions of the Reynolds numbers, for integration into an existing BFM and ready prediction of seal leakage and force coefficients. The selected groove-land section is well within the seal length and far away from the effects of the inlet condition. The analysis takes three water lubricated seals with distinct groove shapes: rectangular, scalloped and triangular. Each seal, with length/diameter L/D = 0.4, has 44 grooves of shallow depth dg ∼ clearance Cr, and operates at a rotor speed equal to 5,588 rpm (78 m/s surface speed) and with a pressure drop of 14.9 MPa. The method validity is asserted when 2D (single groove-land) and 3D (whole seal) predictions for pressure and velocity fields are compared against each other. The CFD predictions, 2D and 3D, show the triangular groove seal has the largest leakage, 41% greater than the rectangular groove seal does, albeit producing the smallest cross-coupled stiffnesses and whirl frequency ratio. On the other hand, the triangular groove seal has the largest direct stiffness and damping coefficients. The scalloped groove seal shows similar rotordynamic force coefficients as the rectangular groove seal but leaks 13% more. For the three seal groove types, the modified BFM predicts leakage that is less than 6% away from that delivered by CFD, whereas the seal stiffnesses (both direct and cross-coupled) differ by 13%, the direct damping coefficients by 18%, and the added mass coefficients are within 30%. The procedure introduced extends the applicability of a BFM to predict the dynamic performance of grooved seals with distinctive shapes.


Author(s):  
Luis San Andrés ◽  
Tingcheng Wu ◽  
Hideaki Maeda ◽  
Ono Tomoki

In straight-through centrifugal pumps, a grooved seal acts as a balance piston to equilibrate the full pressure rise across the pump. As the groove pattern breaks the development of fluid swirl, this seal type offers lesser leakage and lower cross-coupled stiffnesses than a similar size and clearance annular seal. Bulk-flow models predict expediently the static and dynamic force characteristics of annular seals; however they lack accuracy for grooved seals. Computational fluid dynamics (CFD) methods give more accurate results, but are not computationally efficient. This paper presents a modified bulk-flow model to predict the rotordynamic force coefficients of shallow depth circumferentially grooved liquid seals with an accuracy comparable to a CFD solution but with a simulation time of bulk-flow analyses. The procedure utilizes the results of CFD to evaluate the bulk flow velocity field and the friction factors for a 73 grooves annular seal (depth/clearance dg/ Cr = 0.98 and length/diameter L/D = 0.9) operating under various sets of axial pressure drop and rotor speed. In a groove, the flow divides into a jet through the film land and a strong recirculation zone. The penetration angle (α), specifying the streamline separation in the groove cavity, is a function of the operating conditions; an increase in rotor speed or a lower pressure difference increases α. This angle plays a prominent role to evaluate the stator friction factor and has a marked influence on the seal direct stiffness. In the bulk-flow code the friction factor model (f = nRem) is modified with the CFD extracted penetration angle (α) to account for the flow separation in the groove cavity. The flow rate predicted by the modified bulk-flow code shows good agreement with a measured result (6% difference). A perturbation of the flow field is performed on the bulk-flow equations to evaluate the reaction forces on the rotor surface. Compared to the rotordynamic force coefficients derived from the CFD results, the modified bulk-flow code predicts rotordynamic force coefficients within 10%, except that the cross-coupled damping coefficient is over-predicted up to 14%. An example test seal with a few grooves (L/D = 0.5, dg/Cr = 2.5) serves to further validate the predictions of the modified bulk-flow model. Compared to the original bulk-flow analysis, the current method shows a significant improvement in the predicted rotordynamic force coefficients, the direct stiffness and damping coefficients in particular.


Author(s):  
Luis San Andrés ◽  
Tingcheng Wu ◽  
Hideaki Maeda ◽  
Ono Tomoki

In straight-through centrifugal pumps, a grooved seal acts as a balance piston to equilibrate the full pressure rise across the pump. As the groove pattern breaks the development of fluid swirl, this seal type offers lesser leakage and lower cross-coupled stiffnesses than a similar size and clearance annular seal. Bulk-flow models (BFMs) predict expediently the static and dynamic force characteristics of annular seals; however they lack accuracy for grooved seals. Computational fluid dynamics (CFD) methods give more accurate results, but are not computationally efficient. This paper presents a modified BFM to predict the rotordynamic force coefficients of shallow depth, circumferentially grooved liquid seals with an accuracy comparable to a CFD solution but with a simulation time of bulk-flow analyses. The procedure utilizes the results of CFD to evaluate the bulk flow velocity field and the friction factors for a 73 grooves annular seal (depth/clearance dg/Cr = 0.98 and length/diameter L/D = 0.9) operating under various sets of axial pressure drop and rotor speed. In a groove, the flow divides into a jet through the film land and a strong recirculation zone. The penetration angle (α), specifying the streamline separation in the groove cavity, is a function of the operating conditions; an increase in rotor speed or a lower pressure difference increases α. This angle plays a prominent role to evaluate the stator friction factor and has a marked influence on the seal direct stiffness. In the bulk-flow code, the friction factor model (f = nRem) is modified with the CFD extracted penetration angle (α) to account for the flow separation in the groove cavity. The flow rate predicted by the modified bulk-flow code shows good agreement with the measured result (6% difference). A perturbation of the flow field is performed on the bulk-flow equations to evaluate the reaction forces on the rotor surface. Compared to the rotordynamic force coefficients derived from the CFD results, the modified bulk-flow code predicts rotordynamic force coefficients within 10%, except that the cross-coupled damping coefficient is over-predicted up to 14%. An example test seal with a few grooves (L/D = 0.5, dg/Cr = 2.5) serves to further validate the predictions of the modified BFM. Compared to the original bulk-flow analysis, the current method shows a significant improvement in the predicted rotordynamic force coefficients, the direct stiffness and damping coefficients, in particular.


2021 ◽  
pp. 1-31
Author(s):  
Xueliang Lu ◽  
Luis San Andres ◽  
Jing Yang

Abstract Seals in multiple phase rotordynamic pumps must operate without compromising system efficiency and stability. Both field operation and laboratory experiments show that seals supplied with a gas in liquid mixture (bubbly flow) can produce rotordynamic instability and excessive rotor vibrations. This paper advances a nonhomogeneous bulk flow model (NHBFM) for the prediction of the leakage and dynamic force coefficients of uniform clearance annular seals lubricated with gas in liquid mixtures. Compared to a homogeneous BFM (HBFM), the current model includes diffusion coefficients in the momentum transport equations and a field equation for the transport of the gas volume fraction (GVF). Published experimental leakage and dynamic force coefficients for two seals supplied with an air in oil mixture whose GVF varies from 0 (pure liquid) to 20% serve to validate the novel model as well as to benchmark it against predictions from a HBFM. The first seal withstands a large pressure drop (~ 38 bar) and the shaft speed equals 7.5 krpm. The second seal restricts a small pressure drop (1.6 bar) as the shaft turns at 3.5 krpm. The first seal is typical as a balance piston whereas the second seal is found as a neck-ring seal in an impeller. For the high pressure seal and inlet GVF = 0.1, the flow is mostly homogeneous as the maximum diffusion velocity at the seal exit plane is just ~0.1% of the liquid flow velocity. Thus, both the NHBFM and HBFM predict similar flow fields, leakage (mass flow rate) and drag torque. The difference between the predicted leakage and measurement is less than 5%. The NHBFM direct stiffness (K) agrees with the experimental results and reduces faster with inlet GVF than the HBFM K. Both direct damping (C) and cross-coupled stiffness (k) increase with inlet GVF < 0.1.Compared to the test data, the two models generally under predict C and k by ~ 25%. Both models deliver a whirl frequency ratio (fw) ~ 0.3 for the pure liquid seal, hence closely matching the test data. fw raises to ~0.35 as the GVF approaches 0.1. For the low pressure seal the flow is laminar, the experimental results and both NHBFM and HBFM predict a null direct stiffness (K). At an inlet GVF = 0.2, the NHBFM predicted added mass (M) is ~30 % below the experimental result while the HBFM predicts a null M. C and k predicted by both models are within the uncertainty of the experimental results. For operation with either a pure liquid or a mixture (GVF = 0.2), both models deliver fw = 0.5 and equal to the experimental finding. The comparisons of predictions against experimental data demonstrate the NHBFM offers a marked improvement, in particular for the direct stiffness (K). The predictions reveal the fluid flow maintains the homogeneous character known at the inlet condition.


2019 ◽  
Vol 141 (6) ◽  
Author(s):  
Jongin Yang ◽  
Alan Palazzolo

Part II presents a novel approach for predicting dynamic coefficients for a tilting pad journal bearing (TPJB) using computational fluid dynamics (CFD) and finite element method (FEM), including fully coupled elastic deflection, heat transfer, and fluid dynamics. Part I presented a similarly novel, high fidelity approach for TPJB static response prediction which is a prerequisite for the dynamic characteristic determination. The static response establishes the equilibrium operating point values for eccentricity, attitude angle, deflections, temperatures, pressures, etc. The stiffness and damping coefficients are obtained by perturbing the pad and journal motions about this operating point to determine changes in forces and moments. The stiffness and damping coefficients are presented in “synchronously reduced form” as required by American Petroleum Institute (API) vibration standards. Similar to Part I, an advanced three-dimensional thermal—Reynolds equation code validates the CFD code for the special case when flow Between Pad (BP) regions is ignored, and the CFD and Reynolds pad boundary conditions are made identical. The results show excellent agreement for this validation case. Similar to the static response case, the dynamic characteristics from the Reynolds model show large discrepancies compared with the CFD results, depending on the Reynolds mixing coefficient (MC). The discrepancies are a concern given the key role that stiffness and damping coefficients serve instability and response predictions in rotordynamics software. The uncertainty of the MC and its significant influence on static and dynamic response predictions emphasizes a need to utilize the CFD approach for TPJB simulation in critical machines.


2020 ◽  
Vol 142 (9) ◽  
Author(s):  
Zhigang Li ◽  
Zhi Fang ◽  
Jun Li

Abstract Liquid annular seals with parallelly grooved stator or rotor are used as replacements for smooth plain seals in centrifugal pumps to reduce leakage and break up contaminants within the working fluid. Parallelly grooved liquid annular seals have advantages of less leakage and smaller possibility of abrasion when the seal rotor–stator rubs in comparison to smooth plain seals. This paper deals with the static and rotordynamic characteristics of parallelly grooved liquid annular seals, which are limited in the literature. Numerical results of leakage flow rates, drag powers, and rotordynamic force coefficients were presented and compared for a grooved-stator/smooth-rotor (GS-SR) liquid annular seal and a smooth-stator/grooved-rotor (SS-GR) liquid annular seal, utilizing a modified transient computational fluid dynamics-based perturbation approach based on the multiple-frequency elliptical-orbit rotor whirling model. Both liquid annular seals have identical seal axial length, rotor diameter, sealing clearance, groove number, and geometry. The present transient computational fluid dynamics-based perturbation method was adequately validated based on the published experiment data of leakage flow rates and frequency-independent rotordynamic force coefficients for the GS-SR and SS-GR liquid annular seals at various pressure drops with differential inlet preswirl ratios. Simulations were performed at three pressure drops (4.14 bar, 6.21 bar, and 8.27 bar), three rotational speeds (2 krpm, 4 krpm, and 6 krpm) and three inlet preswirl ratios (0, 0.5, and 1.0), applying a wide rotor whirling frequency range up to 200 Hz, to analyze and compare the influences of operation conditions on the static and rotordynamic characteristics for both the GS-SR and SS-GR liquid annular seals. Results show that the present two liquid annular seals possess similar sealing capability, and the SS-GR seal produces a slightly larger (∼2–10%) drag power loss than the GS-SR seal. For small rotor whirling motion around a centered position, both seals have the identical direct force coefficients and the equal-magnitude opposite-sign cross-coupling force coefficients in the orthogonal directions x and y. For all operation conditions, both the GS-SR and SS-GR liquid annular seals possess negative direct stiffness K and positive direct damping C. The GS-SR seal produces purely positive Ceff throughout the whirling frequency range for all operation conditions, while Ceff for the SS-GR seal shows a significant decrease and transitions to negative value at the crossover frequency fco with increasing rotational speed and inlet preswirl. From a rotordynamic viewpoint, the GS-SR liquid annular seal is a better seal concept for pumps.


Author(s):  
Cori Watson ◽  
Houston Wood

Helical groove seals are non-contacting annular seals used in pumping machinery to increase the efficiency and, in the case of the balance drum, to manage the axial force on the thrust bearing. Prior work has shown that optimization of helical groove seals can reduce the leakage by two thirds given a desired pressure differential or, conversely, can significantly increase the pressure differential across the helical groove seal given a flow rate. This study evaluates the dependency of the optimal helical groove seal design on the inlet preswirl, which is the ratio of the inlet circumferential velocity to the rotor surface speed. To accomplish this goal, second stage optimization from the previously optimized helical groove seal with grooves on the stator and water as the working fluid were conducted at a series of preswirls ranging from −1 to 1. Optimization is performed using ANSYS CFX, a commercial computational fluid dynamics software and mesh independence is confirmed for the baseline case. For each preswirl case, design of experiments for the design parameters of groove width, groove depth, groove spacing, and number of grooves was performed using a Kennard-Stone Algorithm. The optimized solution is interpolated from the simulations run by using multi-factor quadratic regression from the 30 simulations in each optimization and the interpolated solution is simulated for comparison. In addition to evaluating the optimized solution’s dependency on preswirl, the viability of using swirl breaks or swirl promoting inlet passages to improve the overall efficiency of the seal is discussed. Finally, the power loss performance is evaluated for each of the seal designs simulated so that potential trade-offs can be evaluated. Overall, the results show that increasing preswirl can increase the efficiency of the helical groove seal both by improving power loss and by improving leakage.


Author(s):  
Luis San Andre´s ◽  
Thomas Soulas ◽  
Florence Challier ◽  
Patrice Fayolle

The paper introduces a bulk-flow model for prediction of the static and dynamic force coefficients of angled injection Lomakin bearings. The analysis accounts for the flow interaction between the injection orifices, the supply circumferential groove, and the thin film lands. A one control-volume model in the groove is coupled to a bulk-flow model within the film lands of the bearing. Bernoulli-type relationships provide closure at the flow interfaces. Flow turbulence is accounted for with shear stress parameters and Moody’s friction factors. The flow equations are solved numerically using a robust computational method. Comparisons between predictions and experimental results for a tangential-against-rotation injection water Lomakin bearing show the novel model predicts well the leakage and direct stiffness and damping coefficients. Computed cross-coupled stiffness coefficients follow the experimental trends for increasing rotor speeds and supply pressures, but quantitative agreement remains poor. A parameter investigation evidences the effects of the groove and land geometries on the Lomakin bearing flowrate and force coefficients. The orifice injection angle does not influence the bearing static performance, although it largely affects its stability characteristics through the evolution of the cross-coupled stiffnesses. The predictions confirm the promising stabilizing effect of the tangential-against-rotation injection configuration. Two design parameters, comprising the feed orifices area and groove geometry, define the static and dynamic performance of Lomakin bearing. The analysis also shows that the film land clearance and length have a larger impact on the Lomakin bearing rotordynamic behavior than its groove depth and length.


Author(s):  
Luis San Andrés ◽  
Joshua Norsworthy

High speed rotors supported on bump-type foil bearings (BFBs) often suffer from large subsynchronous whirl motions. Mechanically preloading BFBs through shimming is a common, low cost practice that shows improvements in rotordynamic stability. However, there is an absence of empirical information related to the force coefficients (structural and rotordynamic) of shimmed BFBs. This paper details a concerted study toward assessing the effect of shimming on a first generation BFB (L = 38.1 mm and D = 36.5 mm). Three metal shims, 120 deg apart, are glued to the inner surface of the bearing cartridge and facing the underside of the bump foil strip. The shim sets are of identical thickness, either 30 μm or 50 μm. In static load tests, a bearing with shims shows a (nonlinear) structural stiffness larger than for the bearing without shims. Torque measurements during shaft acceleration also demonstrate a shimmed BFB has a larger friction coefficient. For a static load of 14.3 kPa, dynamic loads with a frequency sweep from 250 Hz to 450 Hz are exerted on the BFB, without and with shims, to estimate its rotordynamic force coefficients while operating at ∼50 krpm (833 Hz). Similar measurements are conducted without shaft rotation. Results are presented for the original BFB (without shims) and the two shimmed BFB configurations. The direct stiffnesses of the BFB, shimmed or not, increase with excitation frequency, thus evidencing a mild hardening effect. The BFB stiffness and damping coefficients decrease slightly for operation with rotor speed as opposed to the coefficients when the shaft is stationary. For frequencies above 300 Hz, the direct damping coefficients of the BFB with 50 μm thick shims are ∼30% larger than the coefficients of the original bearing. The bearing structural loss factor, a measure of its ability to dissipate mechanical energy, is derived from the direct stiffness and damping coefficients. The BFB with 50 μm thick shims has a 25% larger loss factor—average from test data collected at 300 Hz to 400 Hz—than the original BFB. Further measurements of rotor motions while the shaft accelerates to ∼50 krpm demonstrate the shimmed BFB (thickest shim set) effectively removes subsynchronous whirl motions amplitudes that were conspicuous when operating with the original bearing.


Sign in / Sign up

Export Citation Format

Share Document