A NONHOMOGENEOUS BULK FLOW MODEL FOR GAS IN LIQUID FLOW ANNULAR SEALS: AN EFFORT TO PRODUCE ENGINEERING RESULTS

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.

Author(s):  
Luis San Andrés ◽  
Adolfo Delgado

Oil seals in centrifugal compressors reduce leakage of the process gas into the support bearings and ambient. Under certain operating conditions of speed and pressure, oil seals lock, becoming a source of hydrodynamic instability due to excessively large cross coupled stiffness coefficients. It is a common practice to machine circumferential grooves, breaking the seal land, to isolate shear flow induced film pressures in contiguous lands, and hence reducing the seal cross coupled stiffnesses. Published tests results for oil seal rings shows that an inner land groove, shallow or deep, does not actually reduce the cross-stiffnesses as much as conventional models predict. In addition, the tested grooved oil seals evidenced large added mass coefficients while predictive models, based on classical lubrication theory, neglect fluid inertia effects. This paper introduces a bulk-flow model for groove oil seals operating eccentrically and its solution via the finite element (FE) method. The analysis relies on an effective groove depth, different from the physical depth, which delimits the upper boundary for the squeeze film flow. Predictions of rotordynamic force coefficients are compared to published experimental force coefficients for a smooth land seal and a seal with a single inner groove with depth equaling 15 times the land clearance. The test data represent operation at 10 krpm and 70 bar supply pressure, and four journal eccentricity ratios (e/c= 0, 0.3, 0.5, 0.7). Predictions from the current model agree with the test data for operation at the lowest eccentricities (e/c= 0.3) with discrepancies increasing at larger journal eccentricities. The new flow model is a significant improvement towards the accurate estimation of grooved seal cross-coupled stiffnesses and added mass coefficients; the latter was previously ignored or largely under predicted.


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.


1999 ◽  
Vol 121 (2) ◽  
pp. 363-369 ◽  
Author(s):  
Jiming Li ◽  
David Ransom ◽  
Luis San Andre´s ◽  
John Vance

Experiments and field applications have demonstrated that multiple-pocket gas damper seals effectively eliminate subsynchronous vibration and attenuate imbalance response at the critical speeds in turbomachinery. A one-control volume, turbulent bulk-flow model for the prediction of the seal leakage and rotordynamic force coefficients of centered multiple-pocket damper seals is hereby detailed. Comparisons of numerical predictions with experimental force coefficients for a four-pocket damper seal are presented. The bulk-flow model and experiments indicate the seal direct stiffness and damping force coefficients are insensitive to journal speed while the cross-coupled stiffnesses increase slightly. However, the current model overpredicts the direct damping coefficient and underpredicts the direct stiffness coefficient for increasing test pressure ratios. Computed results show that the force coefficients of multiple-pocket gas damper seals are also functions of the rotor excitation frequency.


2002 ◽  
Vol 129 (1) ◽  
pp. 195-204
Author(s):  
Luis San Andrés ◽  
Thomas Soulas ◽  
Patrice Fayolle

This 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 that novel model well predicts 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 shows evidence of 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, comprised of 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 ◽  
Jing Yang ◽  
Xueliang Lu

Subsea pumps and compressors must withstand multiphase flows whose gas volume fraction (GVF) or liquid volume fraction (LVF) varies over a wide range. Gas or liquid content as a dispersed phase in the primary stream affects the leakage, drag torque, and dynamic forced performance of secondary flow components, namely seals, thus affecting the process efficiency and mechanical reliability of pumping/compressing systems, in particular during transient events with sudden changes in gas (or liquid) content. This paper, complementing a parallel experimental program, presents a computational fluid dynamics (CFD) analysis to predict the leakage, drag power, and dynamic force coefficients of a smooth surface, uniform clearance annular seal supplied with air in oil mixture whose inlet GVF varies discretely from 0.0 to 0.9, i.e., from a pure liquid stream to a nearly all-gas content mixture. The test seal has uniform radial clearance Cr = 0.203 mm, diameter D = 127 mm, and length L = 0.36 D. The tests were conducted with an inlet pressure/exit pressure ratio equal to 2.5 and a rotor surface speed of 23.3 m/s (3.5 krpm), similar to conditions in a pump neck wear ring seal. The CFD two-phase flow model, first to be anchored to test data, uses an Euler–Euler formulation and delivers information on the precise evolution of the GVF and the gas and liquid streams' velocity fields. Recreating the test data, the CFD seal mass leakage and drag power decrease steadily as the GVF increases. A multiple-frequency shaft whirl orbit method aids in the calculation of seal reaction force components, and from which dynamic force coefficients, frequency-dependent, follow. For operation with a pure liquid, the CFD results and test data produce a constant cross-coupled stiffness, damping, and added mass coefficients, while also verifying predictive formulas typical of a laminar flow. The injection of air in the oil stream, small or large in gas volume, immediately produces force coefficients that are frequency-dependent; in particular the direct dynamic stiffness which hardens with excitation frequency. The effect is most remarkable for small GVFs, as low as 0.2. The seal test direct damping and cross-coupled dynamic stiffness continuously drop with an increase in GVF. CFD predictions, along with results from a bulk-flow model (BFM), reproduce the test force coefficients with great fidelity. Incidentally, early engineering practice points out to air injection as a remedy to cure persistent (self-excited) vibration problems in vertical pumps, submersible and large size hydraulic. Presently, the model predictions, supported by the test data, demonstrate that even a small content of gas in the liquid stream significantly raises the seal direct stiffness, thus displacing the system critical speed away to safety. The sound speed of a gas in liquid mixture is a small fraction of those speeds for either the pure oil or the gas, hence amplifying the fluid compressibility that produces the stiffness hardening. The CFD model and a dedicated test rig, predictions and test data complementing each other, enable engineered seals for extreme applications.


2002 ◽  
Vol 129 (1) ◽  
pp. 185-194 ◽  
Author(s):  
Thomas Soulas ◽  
Luis San Andres

A computational analysis for prediction of the static and dynamic forced performance of gas honeycomb seals at off-centered rotor conditions follows. The bulk-flow analysis, similar to the two-control volume flow model of Kleynhans and Childs (1997, “The Acoustic Influence of Cell Depth on the Rotordynamic Characteristics of Smooth-Rotor/Honeycomb-Stator Annular Gas Seals,” ASME J. Eng. Gas Turbines Power, 119, pp. 949–957), is brought without loss of generality into a single-control volume model, thus simplifying the computational process. The formulation accommodates the honeycomb effective cell depth, and existing software for annular pressure seals and is easily upgraded for damper seal analysis. An analytical perturbation method for derivation of zeroth- and first-order flow fields renders the seal equilibrium response and frequency-dependent dynamic force impedances, respectively. Numerical predictions for a centered straight-bore honeycomb gas seal shows good agreement with experimentally identified impedances, hence validating the model and confirming the paramount influence of excitation frequency on the rotordynamic force coefficients of honeycomb seals. The effect of rotor eccentricity on the static and dynamic forced response of a smooth annular seal and a honeycomb seal is evaluated for characteristic pressure differentials and rotor speeds. Leakage for the two seal types increases slightly as the rotor eccentricity increases. Rotor off-centering has a pronounced nonlinear effect on the predicted (and experimentally verified) dynamic force coefficients for smooth seals. However, in honeycomb gas seals, even large rotor center excursions do not sensibly affect the effective local film thickness, maintaining the flow azimuthal symmetry. The current model and predictions thus increase confidence in honeycomb seal design, operating performance, and reliability in actual applications.


Author(s):  
Thomas Soulas ◽  
Luis San Andres

A computational analysis for prediction of the static and dynamic forced performance of gas honeycomb seals at off-centered rotor conditions follows. The bulk-flow analysis, similar to the two-control volume flow model of Kleynhans and Childs [1], is brought without loss of generality into a single-control volume model, thus simplifying the computational process. The formulation accommodates the honeycomb effective cell depth, and existing software for annular pressure seals is easily upgraded for damper seal analysis. An analytical perturbation method for derivation of zeroth- and first-order flow fields renders the seal equilibrium response and frequency-dependent dynamic force impedances, respectively. Numerical predictions for a centered straight-bore honeycomb gas seal show good agreement with experimentally identified impedances, hence validating the model and confirming the paramount influence of excitation frequency on the rotordynamic force coefficients of honeycomb seals. The effect of rotor eccentricity on the static and dynamic forced response of a smooth annular seal and a honeycomb seal is evaluated for characteristic pressure differentials and rotor speeds. Leakage for the two seal types increases slightly as the rotor eccentricity increases. Rotor off-centering does have a pronounced non-linear effect on the predicted (and experimentally verified) dynamic force coefficients for smooth seals. However, in honeycomb gas seals, even large rotor center excursions do not sensibly affect the effective local film thickness, maintaining the flow azimuthal symmetry. The current model and predictions thus increase confidence in honeycomb seal design, operating performance and reliability in actual applications.


Author(s):  
Luis San Andre´s ◽  
Adolfo Delgado

Oil seals in centrifugal compressors reduce leakage of the process gas into the support bearings and ambient. Under certain operating conditions of speed and pressure, oil seals lock, becoming a source of hydrodynamic instability due to excessively large cross coupled stiffness coefficients. It is a common practice to machine circumferential grooves, breaking the seal land, to isolate shear flow induced film pressures in contiguous lands, and hence reducing the seal cross coupled stiffnesses. Published tests results for oil seal rings shows that an inner land groove, shallow or deep, does not actually reduce the cross-stiffnesses as much as conventional models predict. In addition, the tested grooved oil seals evidenced large added mass coefficients; while predictive models, based on classical lubrication theory, neglect fluid inertia effects. This paper introduces a bulk-flow model for groove oil seals operating eccentrically and its solution via the finite element method. The analysis relies on an effective groove depth, different from the physical depth, which delimits the upper boundary for the squeeze film flow. Predictions of rotordynamic force coefficients are compared to published experimental force coefficients for a smooth land seal and a seal with a single inner groove with depth equaling 15 times the land clearance. The test data represent operation at 10 krpm and 70 bar supply pressure, and four journal eccentricity ratios (e/c = 0, 0.3, 0.5, 0.7). Predictions from the current model agree with the test data for operation at the lowest eccentricities (e/c = 0.3); discrepancies increasing at larger journal eccentricities. The new flow model is a significant improvement towards the accurate estimation of grooved seal cross-coupled stiffnesses and added mass coefficients; the later previously ignored or largely under predicted.


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 ◽  
Jing Yang ◽  
Xueliang Lu

Subsea pumps and compressors must withstand multi-phase flows whose gas volume fraction (GVF) or liquid volume fraction (LVF) varies over a wide range. Gas or liquid content as a dispersed phase in the primary stream affects the leakage, drag torque, and dynamic forced performance of secondary flow components, namely seals, thus affecting the process efficiency and mechanical reliability of pumping/compressing systems, in particular during transient events with sudden changes in gas (or liquid) content. This paper, complementing a parallel experimental program, presents a computational fluid dynamics (CFD) analysis to predict the leakage, drag power and dynamic force coefficients of a smooth surface, uniform clearance annular seal supplied with an air in oil mixture whose inlet GVF varies discretely from 0.0 to 0.9, i.e., from a pure liquid stream to a nearly all gas content mixture. The test seal has uniform radial clearance Cr = 0.203 mm, diameter D = 127 mm, and length L = 0.36 D. The tests were conducted with an inlet pressure/exit pressure ratio equal to 2.5 and a rotor surface speed of 23.3 m/s (3.5 krpm), similar to conditions in a pump neck wear ring seal. The CFD two-phase flow model, first to be anchored to test data, uses an Euler-Euler formulation and delivers information on the precise evolution of the GVF and the gas and liquid streams’ velocity fields. Recreating the test data, the CFD seal mass leakage and drag power decrease steadily as the GVF increases. A multiple-frequency shaft whirl orbit method aids in the calculation of seal reaction force components, and from which dynamic force coefficients, frequency dependent, follow. For operation with a pure liquid, the CFD results and test data produce a constant cross-coupled stiffness, damping, and added mass coefficients, while also verifying predictive formulas typical of a laminar flow. The injection of air in the oil stream, small or large in gas volume, immediately produces force coefficients that are frequency dependent; in particular the direct dynamic stiffness which hardens with excitation frequency. The effect is most remarkable for small GVFs, as low as 0.2. The seal test direct damping and cross-coupled dynamic stiffness continuously drop with an increase in GVF. CFD predictions, along with results from a bulk-flow model (BFM), reproduce the test force coefficients with great fidelity. Incidentally, early engineering practice points out to air injection as a remedy to cure persistent (self-excited) vibration problems in vertical pumps, submersible and large size hydraulic. Presently, the model predictions, supported by the test data, demonstrate that even a small content of gas in the liquid stream significantly raises the seal direct stiffness, thus displacing the system critical speed away to safety. The sound speed of a gas in liquid mixture is a small fraction of those speeds for either the pure oil or the gas, hence amplifying the fluid compressibility that produces the stiffness hardening. The CFD model and a dedicated test rig, predictions and test data complementing each other, enable engineered seals for extreme applications.


Sign in / Sign up

Export Citation Format

Share Document