A Computational Fluid Dynamics Modified Bulk Flow Analysis for Circumferentially Shallow Grooved Liquid Seals

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.

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.


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.


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):  
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.


2013 ◽  
Vol 136 (1) ◽  
Author(s):  
C. I. Papadopoulos ◽  
L. Kaiktsis ◽  
M. Fillon

The paper presents a detailed computational study of flow patterns and performance indices in a dimpled parallel thrust bearing. The bearing consists of eight pads; the stator surface of each pad is partially textured with rectangular dimples, aiming at maximizing the load carrying capacity. The bearing tribological performance is characterized by means of computational fluid dynamics (CFD) simulations, based on the numerical solution of the Navier–Stokes and energy equations for incompressible flow. Realistic boundary conditions are implemented. The effects of operating conditions and texture design are studied for the case of isothermal flow. First, for a reference texture pattern, the effects of varying operating conditions, in particular minimum film thickness (thrust load), rotational speed and feeding oil pressure are investigated. Next, the effects of varying texture geometry characteristics, in particular texture zone circumferential/radial extent, dimple depth, and texture density on the bearing performance indices (load carrying capacity, friction torque, and friction coefficient) are studied, for a representative operating point. For the reference texture design, the effects of varying operating conditions are further investigated, by also taking into account thermal effects. In particular, adiabatic conditions and conjugate heat transfer at the bearing pad are considered. The results of the present study indicate that parallel thrust bearings textured by proper rectangular dimples are characterized by substantial load carrying capacity levels. Thermal effects may significantly reduce load capacity, especially in the range of high speeds and high loads. Based on the present results, favorable texture designs can be assessed.


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):  
Koya Yamada ◽  
Atsushi Ikemoto ◽  
Tsuyoshi Inoue ◽  
Masaharu Uchiumi

Rotor-dynamic fluid force (RD fluid force) of turbomachinery is one of the causes of the shaft vibration problem. Bulk flow theory is the method for analyzing this RD fluid force, and it has been widely used in the design stage of machine. The conventional bulk flow theory has been carried out under the assumption of concentric circular shaft's orbit with a small amplitude. However, actual rotating machinery's operating condition often does not hold this assumption, for example, existence of static load on the machinery causes static eccentricity. In particular, when such a static eccentricity is significant, the nonlinearity of RD fluid force may increase and become non-negligible. Therefore, conventional bulk flow theory is not applicable for the analysis of the RD fluid force in such a situation. In this paper, the RD fluid force of the annular plain seal in the case of circular whirling orbit with static eccentricity is investigated. The case with both the significant static eccentricity and the moderate whirling amplitude is considered, and the perturbation analysis of the bulk-flow theory is extended to investigate the RD fluid force in such cases. In this analysis, the assumption of the perturbation solution is extended to both static terms and whirling terms up to the third order. Then, the additional terms are caused by the coupling of these terms through nonlinearity, and these three kinds of terms are considered in the extended perturbation analysis of the bulk flow theory. As a result, a set of nonlinear analytical equations of the extended perturbation analysis of the bulk flow theory, for the case with both the significant static eccentricity and the moderate whirling amplitude, is deduced. The RD fluid force for such cases is analyzed, and the occurrence of constant component, backward synchronous component, and super-harmonic components in the RD fluid force is observed in addition to the forward synchronous component. The representation of RD fluid force coefficients (RD coefficients) are modified for the case with significant static eccentricity, and the variation of RD fluid force coefficients for the magnitude of static eccentricity is analyzed. These analytical results of RD fluid force and its RD coefficients are compared with the numerical results using finite difference analysis and experimental results. As a result, the validity of the extended perturbation analysis of the bulk-flow theory for the case with both the significant static eccentricity and the moderate whirling amplitude is confirmed.


Author(s):  
Riccardo Da Soghe ◽  
Cosimo Bianchini ◽  
Antonio Andreini ◽  
Lorenzo Mazzei ◽  
Giovanni Riccio ◽  
...  

Combustor liner of present gas turbine engines is subjected to high thermal loads as it surrounds high temperature combustion reactants and is hence facing the related radiative load. This generally produces high thermal stress levels on the liner, strongly limiting its life expectations and making it one of the most critical components of the entire engine. The reliable prediction of such thermal loads is hence a crucial aspect to increase the flame tube life span and to ensure safe operations. The present study aims at investigating the aerothermal behavior of a GE Dry Low NOx (DLN1) class flame tube and in particular at evaluating working metal temperatures of the liner in relation to the flow and heat transfer state inside and outside the combustion chamber. Three different operating conditions have been accounted for (i.e., lean–lean partial load, premixed full load, and primary load) to determine the amount of heat transfer from the gas to the liner by means of computational fluid dynamics (CFD). The numerical predictions have been compared to experimental measurements of metal temperature showing a good agreement between CFD and experiments.


Sign in / Sign up

Export Citation Format

Share Document