A Novel Bulk-Flow Model for Improved Predictions of Force Coefficients in Grooved Oil Seals Operating Eccentrically

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 ◽  
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.


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 ◽  
Yingkun Li

Tilting pad journal bearings (TPJBs) supporting high-performance turbomachinery rotors have undergone steady design improvements to satisfy ever stringent operating conditions that include large specific loads, due to smaller footprints, and high surface speeds that promote flow turbulence and hence larger drag power losses. Simultaneously, predictive models continuously evolve to include minute details on bearing geometry, pads and pivots' configurations, oil delivery systems, etc. In general, predicted TPJB rotordynamic force coefficients correlate well with experimental data for operation with small to moderately large unit loads (1.7 MPa). Experiments also demonstrate bearing dynamic stiffnesses are frequency dependent, best fitted with a stiffness-mass like model whereas damping coefficients are adequately represented as of viscous type. However, for operation with large specific loads (>1.7 MPa), poor correlation of predictions to measured force coefficients is common. Recently, an experimental effort (Gaines, J., 2014, “Examining the Impact of Pad Flexibility on the Rotordynamic Coefficients of Rocker-Pivot-Pad Tiling-Pad Journal Bearings,” M.S. thesis, Mechanical Engineering, Texas A&M University, College Station, TX) produced test data for three TPJB sets, each having three pads of unequal thickness, to quantify the effect of pad flexibility on the bearings' force coefficients, in particular damping, over a range of load and rotational speed conditions. This paper introduces a fluid film flow model accounting for both pivot and pad flexibility to predict the bearing journal eccentricity, drag power loss, lubricant temperature rise, and force coefficients of typical TPJBs. A finite element (FE) pad structural model including the Babbitt layer is coupled to the thin film flow model to determine the mechanical deformation of the pad surface. Predictions correlate favorably with test data, also demonstrating that pad flexibility produces a reduction of up to 34% in damping for the bearing with the thinnest pads relative to that with the thickest pads. A parametric study follows to quantify the influence of pad thickness on the rotordynamic force coefficients of a sample TPJB with three pads of increasing preload, r¯p  = 0, 0.25 (baseline) and 0.5. The bearing pads are either rigid or flexible by varying their thickness. For design considerations, dimensionless static and dynamic characteristics of the bearings are presented versus the Sommerfeld number (S). Pad flexibility shows a more pronounced effect on the journal eccentricity and the force coefficients of a TPJB with null pad preload than for the bearings with larger pad preloads (0.25 and 0.5), in particular for operation with a small load or at a high surface speed (S > 0.8).


Author(s):  
Luis San Andrés ◽  
Yingkun Li

Tilting pad journal bearings (TPJBs) supporting high performance turbomachinery rotors have undergone steady design improvements to satisfy ever stringent operating conditions that include large specific loads due to smaller footprints, and high surface speeds that promote flow turbulence and thus larger drag power losses. Simultaneously, predictive models continuously evolve to include minute details on bearing geometry, pads and pivots’ configurations, oil delivery systems, etc. In general, predicted TPJB rotordynamic force coefficients correlate well with experimental data for operation with small to moderately large unit loads (1.7 MPa). Experiments also demonstrate bearing dynamic stiffnesses are frequency dependent, best fitted with a stiffness-mass like model whereas damping coefficients are adequately represented as of viscous type. However, for operation with large specific loads (> 1.7 MPa), poor correlation of predictions to measured force coefficients is common. Recently, an experimental effort [1] produced test data for three TPJB sets, each having three pads of unequal thickness, to quantify the effect of pad flexibility on the bearings’ force coefficients, in particular damping, over a range of load and rotational speed conditions. This paper introduces a fluid film flow model accounting for both pivot and pad flexibility to predict the bearing journal eccentricity, drag power loss, lubricant temperature rise and force coefficients of typical TPJBs. A finite element pad structural model including the Babbitt layer is coupled to the thin film flow model to determine the mechanical deformation of the pad surface. Predictions correlate favorably with test data, also demonstrating that pad flexibility produces a reduction of up to 34% in damping for the bearing with the thinnest pads relative to that with the thickest pads. A parametric study follows to quantify the influence of pad thickness on the rotordynamic force coefficients of a sample TPJB with three pads of increasing preload, rp = 0, 0.25 (baseline) and 0.5. The bearing pads are either rigid or flexible by varying their thickness. For design considerations, dimensionless static and dynamic characteristics of the bearings are presented versus the Sommerfeld number (S). Pad flexibility shows a more pronounced effect on the journal eccentricity and the force coefficients of a TPJB with null pad preload than for the bearings with larger pad preloads (0.25 and 0.5), in particular for operation with a small load or at a high surface speed (S>0.8).


Author(s):  
Luis San Andrés ◽  
Bonjin Koo

Abstract Practice and experiments with squeeze film dampers (SFDs) sealed with piston rings (PRs) show the lubricant exits through the PR slit, i.e., the gap made by the PR abutted ends when installed, forced as a jet during the portion of a rotor whirl cycle generating a positive squeeze film pressure. In the other portion of a whirl cycle, a subambient dynamic pressure ingests air into the film that mixes with the lubricant to produce a bubbly mixture. To reduce persistent air ingestion, commercial air breathing engines utilizing PRSFDs demand of a sufficiently large lubricant supply pressure (Ps), and hence a larger flow rate that is proportional to the journal squeeze velocity (vs = amplitude r × frequency of motion ω). The stringent requirement clearly limits the applicability and usefulness of SFDs. This paper presents a computational physics model for a sealed-end SFD operating with a mixture and delivers predictions benchmarked against profuse laboratory test data. The model implements a Reynolds equation adapted for a homogeneous bubbly mixture, includes temporal fluid inertia effects, and uses physics-based inlet and outlet lubricant conditions through feed holes and PR slit, respectively. In the experiments for model validation, a SFD damper, 127 mm in diameter D, film land length L = 25.4 mm (L/D = 0.2), and radial clearance c = 0.371 mm, is supplied with an air in ISO VG2 oil bubbly mixture of known gas volume fraction (GVF), zero (pure oil) to 50% in steps of 10%. The mixture supply pressure varies from Ps = 2.06 bar-g (30 psig) to 6.20 bar-g (90 psig). Located in grooves at the top and bottom of the journal, a PR and an O-ring (OR) seal the film land. The OR does not allow any oil leakage or air ingestion; hence, the supplied mixture discharges through the PR slit into a vessel submerged within a large volume of lubricant. Dynamic load tests with a single frequency ω, varying from 10 Hz to 60 Hz, produce circular centered orbits (CCO) with amplitude r = 0.2c. The measurements record the exerted forces and journal motions and an analysis delivers force coefficients, damping and inertia, representative of the exerted frequency range. The model predicts the pressure field and evolution of the GVF within the film land and, in a simulated process replicating the experimental procedure, delivers representative force coefficients. For all Ps conditions, both predictions and tests show the SFD added mass coefficients significantly decrease as the inlet GVF (βs) increases. The experimentally derived damping coefficients do not show a significant change, except for tests with the largest concentration of air (βs = 0.5). The predicted damping differs by 10% with the test derived coefficient which does not readily decrease as the inlet GVF (βs) increases. The added mass coefficients, test and predicted, decrease with βs, both being impervious to the magnitude of supply pressure. The test PRSFD shows a quadrature stiffness due to the sliding friction between the PR being pushed against the journal. An increase in supply pressure exacerbates this unique stiffness that may impair the action of the squeeze film to dissipate mechanical energy. The comprehensive test results, first of their kind, demonstrate that accurate modeling of SFDs operating with air ingestion remains difficult as the flow process and the paths of its major components (air and liquid) are rather complex.


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

Aircraft engine rotors are particularly sensitive to rotor imbalance and sudden maneuver loads, since they are always supported on rolling element bearings with little damping. Most engines incorporate squeeze film dampers (SFDs) as means to dissipate mechanical energy from rotor vibrations and to ensure system stability. The paper quantifies experimentally the forced performance of a SFD comprising two parallel film lands separated by a deep central groove. Tests are conducted on two open ends SFDs, both with diameter D = 127 mm and nominal radial clearance c = 0.127 mm. One damper has film lands with length L = 12.7 mm (short length), while the other has 25.4 mm land lengths. The central groove has width L and depth 3/4 L. A light viscosity lubricant flows into the central groove via three orifices, 120 deg apart and then through the film lands to finally exit to ambient. In operation, a static loader pulls the bearing to various eccentric positions and electromagnetic shakers excite the test system with periodic loads to generate whirl orbits of specific amplitudes. A frequency domain method identifies the SFD damping and inertia force coefficients. The long damper generates six times more damping and about three times more added mass than the short length damper. The damping coefficients are sensitive to the static eccentricity (up to ∼ 0.5 c), while showing lesser dependency on the amplitude of whirl motion (up to 0.2 c). On the other hand, inertia coefficients increase mildly with static eccentricity and decrease as the amplitude of whirl motion increases. Cross-coupled force coefficients are insignificant for all imposed operating conditions on either damper. Large dynamic pressures recorded in the central groove demonstrate the groove does not isolate the adjacent squeeze film lands, but contributes to the amplification of the film lands’ reaction forces. Predictions from a novel SFD model that includes flow interactions in the central groove and feed orifices agree well with the test force coefficients for both dampers. The test data and predictions advance current knowledge and demonstrate that SFD-forced performance is tied to the lubricant feed arrangement.


Author(s):  
Je´roˆme Gehannin ◽  
Mihai Arghir ◽  
Olivier Bonneau

The present work presents a bulk flow based numerical analysis of squeeze film dampers (SFD) provided with circumferential grooves, feeding orifices and piston ring sealing devices. This means that the bulk flow model for inertia dominated flow regimes (high Re) was adapted for taking into account thin film discontinuities, mass sources located in the thin film zone (for modelling feeding orifices) or on the boundaries (piston ring seals with open slots). In the present work the rotor is whirling on a centered whirl and the bulk flow equation are solved by following a time step integration.


Author(s):  
Filippo Cangioli ◽  
Steven Chatterton ◽  
Paolo Pennacchi ◽  
Leonardo Nettis ◽  
Lorenzo Ciuchicchi ◽  
...  

Over the last few decades, the increasing demand on efficiency and performance for steam turbines has resulted in OEMs operating machines near critical conditions of their structural and thermal capabilities. Consequently, a more accurate estimation of the dynamic behavior of the machine has become mandatory as well as the stability assessment. Steam turbines are subjected to high temperatures, high pressures and centrifugal forces that could change the nominal geometry, especially the clearance profile in correspondence of the sealing components, occasionally generating a convergent or divergent annulus. In this paper, a new thermo-elasto bulk-flow model for labyrinth seals has been introduced. The model includes the bulk-flow model for estimating the dynamic coefficients, heat transfer model for evaluating the temperature distribution in the rotating and stationary parts and structural-mechanics model for calculating the radial growth. By considering a staggered labyrinth seal installed in the balancing drum of a steam turbine, different inlet pre-swirl ratios, as well as the stability of the seal are investigated in this paper. The model can be extremely useful for the dynamic characterisation of a wide class of labyrinth seals considering the effect of the surrounding environment on the rotordynamic coefficient prediction.


Author(s):  
Filippo Cangioli ◽  
Paolo Pennacchi ◽  
Giacomo Riboni ◽  
Giuseppe Vannini ◽  
Lorenzo Ciuchicchi ◽  
...  

Since the 80s, academic research in the rotordynamics field has developed mathematical treatment for the prediction of the dynamic coefficients of sealing components. Dealing with the straight-through labyrinth seal, Iwatsubo [1], at a first stage, and Childs [2], later on, have developed the one-control volume bulk flow model. The model allows evaluating the surrounding fluid forces acting on the rotor, analyzing the fluid dynamics within the seal: the continuity, circumferential momentum and energy equations are solved for each cavity. To consider axial fluid dynamics, correlations, aiming to estimate the leakage and the pressure distribution, are required. Several correlations have been proposed in the literature for the estimation of the leakage, of the kinetic energy carry-over coefficient (KE), of the discharge coefficient and of the friction factor. After decades of research in the field of seal dynamics, the bulk-flow model has been confirmed as the most popular code in the industries, however, it is not clear which is the best set of correlations for the prediction of seal dynamic coefficients. This paper allows identifying the most accurate combination of correlations to be implemented in the bulk-flow model. The correlations are related to: the leakage formula, the flow coefficient, the KE and the friction factor. Investigating the results of several models (32 models), which consider different sets of correlations, in comparison to the experimental data (performed by General Electric Oil & Gas), it is possible to observe the dependence, of the model correlations, on the operating conditions. The experimental results, performed using a 14 teeth-on-stator labyrinth seal, investigate several operating conditions of pressure drop.


Sign in / Sign up

Export Citation Format

Share Document