Research Article  Open Access
Filippo Cangioli, Paolo Pennacchi, Leonardo Nettis, Lorenzo Ciuchicchi, "Design and Analysis of CFD Experiments for the Development of BulkFlow Model for Staggered Labyrinth Seal", International Journal of Rotating Machinery, vol. 2018, Article ID 9357249, 16 pages, 2018. https://doi.org/10.1155/2018/9357249
Design and Analysis of CFD Experiments for the Development of BulkFlow Model for Staggered Labyrinth Seal
Abstract
Nowadays, bulkflow models are the most timeefficient approaches to estimate the rotor dynamic coefficients of labyrinth seals. Dealing with the onecontrol volume bulkflow model developed by Iwatsubo and improved by Childs, the “leakage correlation” allows the leakage massflow rate to be estimated, which directly affects the calculation of the rotor dynamic coefficients. This paper aims at filling the lack of the numerical modelling for staggered labyrinth seals: a onecontrol volume bulkflow model has been developed and, furthermore, a new leakage correlation has been defined using CFD analysis. Design and analysis of computer experiments have been performed to investigate the leakage massflow rate, static pressure, circumferential velocity, and temperature distribution along the seal cavities. Four design factors have been chosen, which are the geometry, pressure drop, inlet preswirl, and rotor peripheral speed. Finally, dynamic forces, estimated by the bulkflow model, are compared with experimental measurements available in the literature.
1. Introduction
Labyrinth seals and their sealing principles are commonplace in turbomachinery and they can arise in various configurations [1]. The most popular are the straightthrough, the staggered, the slanted, and the stepped labyrinth seals. By their nature, labyrinth seals are noncontact seals that can harmfully rub against stationary parts [2]. Their working principle consists in reducing the leakage massflow rate by dissipating the flow kinetic energy via sequential cavities arranged to impose a tortuous path to the fluid (see the scheme of a staggered labyrinth seal in Figure 1). The speed and pressure at which they operate are bounded by their structural design.
The clearance is defined by aerothermomechanical conditions that stave off the contact against the shroud under radial and axial excursions. In straightthrough labyrinth seal, the angle at which the flow approaches to the teeth is usually 90°; the coefficients of discharge (), peculiar of labyrinth seals, emphasize the effectiveness of the sealing element. The discharge coefficient is defined as the ratio between the actual massflow rate and that in ideal flow conditions. The relation between the sharpness of the tooth and the ability to restrict the flow has been given by Mahler [3] and more recently explored by means of computational fluid dynamics (CFD) in [4].
Labyrinth seals are effective in reducing the flow but have a strong effect on the dynamic behaviour and often lead to dynamic instabilities. These problems have been addressed by several researchers starting with Thomas [5] and Alford [6]. They recognized that dynamic forces can lead to instabilities. Childs [7] addressed the root cause to the swirl velocity and introduced the swirl brake at the seal inlet to mitigate the circumferential velocity of the fluid.
The current trends in the oil and gas market lead the manufacturers to maximize the turbomachinery efficiency. This target pushes the rotor dynamic design to be challenging. Therefore, the accuracy of rotor dynamic analyses has become of primary importance to avoid undesired instability phenomenon. Because labyrinth seals are one of the major sources of destabilizing effects, it is straightforward that the accurate prediction of labyrinth seal dynamic coefficients is crucial.
The bulkflow model is the most used approach for estimating the seal dynamic coefficients thanks to the low computational effort compared to CFD calculations. It was developed by Iwatsubo [8] and improved by Childs and Scharrer [9] in the 80s. The bulkflow model is a control volume approach: the continuity and circumferential momentum equations are solved for each cavity, allowing the circumferential fluid dynamics within the seal to be investigated. For the axial fluid dynamics, empirical correlations are required, aiming at estimating the leakage and the pressure distribution. No turbulence model is employed in the bulkflow model; therefore also the effect of the turbulence is demanded to correlations that estimate the wall shear stresses. The leakage strictly influences the swirl velocity profile along the seal and consequently the resulting dynamic coefficients.
Dealing with staggered labyrinth seals, few models allow the dynamic coefficients to be estimated, despite the fact that this seal configuration is widely used in turbomachinery, especially in steam turbines. The advantage of this configuration, with respect to straightthrough labyrinth seal, is the reduction of the leakage massflow rate given equal clearance.
Wang et al. [10] use the straightthrough bulkflow model for the staggered labyrinth seal; Childs [7] considered the kinetic energy carryover equal to the unity for all the cavities. Kwanka and Ortinger [11] compared the numerical results of the bulkflow theory with experimental results performed on a staggered labyrinth seal. They revealed the strong connection between the circumferential flow and the axial one. The numerical results were not accurate compared to the experimental ones. The leakage correlation used by Kwanka et al. did not reproduce well the “zigzag” behaviour of the discharge coefficients, because they considered an average value. Li et al. [12] investigated numerically and experimentally the effects of the pressure ratio and rotational speed on the leakage flow in a staggered labyrinth seal. In their paper, the authors do not report any correlation for the prediction of the leakage flow.
This paper aims at filling the lack of analytical models for staggered labyrinth seals. Firstly, the authors performed a design and analysis of computer experiments (DACE) using steadystate CFD analysis to develop a new empirical correlation for the prediction of staggered labyrinth seal leakage. The CFD analysis investigates several seal geometries and various operating conditions, which are pressure drop, inlet preswirl, and rotor peripheral speed. To capture the nonlinear dependency of these factors on the flow conditions, three levels each are chosen, except two for the preswirl. The correlation obtained is finally verified numerically, comparing the results of the bulkflow model with those of a CFD analysis in terms of leakage massflow rate, pressure distributions, and Mach number.
In conclusion, the rotor dynamic coefficients prediction has been compared with the experimental measurements performed by Kwanka and Ortinger [11], highlighting the more accurate estimation of the crossforce with respect to other bulkflow models.
2. BulkFlow SteadyState Solution
The estimation of the steady leakage, pressure, and circumferential velocity distribution represents the socalled zerothorder solution in the bulkflow model, as described in [13]. For a fixed pressure drop and assuming, at a first stage, an isenthalpic process, the leakage and the pressure distribution in the seal cavities are determined by solving iteratively the continuity equation, until the massflow rate is equal for each control volume (see Figure 2). Then, the enthalpy and the circumferential velocity in each cavity are determined by solving the coupled circumferential momentum and energy equation. Multivariate NewtonRaphson algorithm is employed to find the solution. At the following step, the leakage is calculated one more time, considering the enthalpy variation previously estimated. The iteration continues until the leakage mass flowrate converges.
The continuity, circumferential momentum, and energy zerothorder equations for each cavity are stated as follows:
Continuity Equation
Circumferential Momentum Equation
Energy Equationwhere is the steady leakage massflow rate, and are the steady circumferential velocity and enthalpy for the th cavity. and are the shear stresses on the rotor and stator parts, whereas and are the lengths of the rotor and stator parts, respectively.
The leakage correlation () used in the model is given by the Neumann formula, as described in [14]:where and are, respectively, the discharge coefficient and kinetic energy carryover coefficient, is the clearance, and and are the pressure and density in the th cavity.
The shear stresses both on the rotor and on the stator ( and ) have been evaluated by implementing the SwameeJain [15] explicit formula for the estimation of the Darcy friction factor. Compared with the Blasius correlation, which is used in several models [9, 16–19], the SwameeJain correlation is more accurate in the calculation of the friction factor, because it takes into account the roughness of the surface and it is validated also for high Reynolds numbers. A comprehensive investigation of the correlation has been performed by Kiijarvi in [20]. The rotor and stator friction factors are defined aswhere is the circumferential hydraulic diameter, is the kinematic viscosity, is the absolute roughness of the rotor/stator surface, and is the rotational speed of the rotor. The Reynolds numbers are defined for the rotor and stator parts, respectively, as follows:
The shear stresses generated by the fluid/structure interaction are evaluated considering both the axial and circumferential velocity, as described in [19] by Wyssmann et al., as follows:where is the axial velocity of the fluid in the th cavity, which is estimated as a function of the calculated massflow rate as
3. Preliminary Results
To investigate the accuracy of the leakage correlation used for straightthrough labyrinth seals in the case of staggered seals, the zerothorder solutions are compared with the CFD results. Steadystate CFD analysis is very accurate in the prediction of the leakage, static pressure, and circumferential velocity within the seal; therefore CFD has been considered as the reference for the zerothorder solution.
A detailed description of the CFD analysis is shown in the next chapter of the paper. The seal geometry is shown in Figure 3, and the geometrical parameters are reported in Table 1, using the nomenclature shown in Figure 4. The working fluid is the steam.

The leakage can be calculated in two manners: considering the kinetic energy carryover coefficient (Wang’s model), as suggested by Wang et al. [10], or assuming that the kinetic energy carryover coefficient is equal to the unity for all the constrictions (Childs’ model), as suggested by Childs in [7]. Both models use the empirical Neumann’s correlation (see (4)) for the leakage [14] and Chaplygin’s formula (see (9)) for the discharge coefficient. The kinetic energy carryover coefficient is calculated using the semiempirical correlation developed by Hodkinson in [21].
The results in terms of leakage massflow rate are shown in Figure 5, while the swirl velocity and the pressure distribution are shown in Figure 6. The massflow rate is overestimated for both bulkflow models: the relative error is, respectively, equal to +75% and +272% for Childs’ model and Wang’s model. The swirl ratios along the seal cavities are strongly influenced by the mass flow. Even if the asymptotic value of the swirl ratio for Childs’ model is equal to that of the CFD, the values in the first cavities are far from the CFD ones. The predicted static pressure distribution is very similar to Childs’ model.
Eventually, Childs’ model is the more accurate onecontrol volume bulkflow model in the literature for staggered labyrinth seals; however, the results in terms of leakage, pressure, and circumferential velocity distribution are not accurate enough compared to CFD results. The zerothorder solution strongly influences the prediction of the seal dynamic coefficients; hence, the authors decide to develop a new empirical correlation for staggered labyrinth seal.
4. Design and Analysis of Computer Experiments
Design and analysis of CFD experiments allow several staggered seal geometries and operating conditions to be investigated for the development of a new leakage correlation. A priori, it is not possible to know which are the geometric parameters that influence the leakage flow rate, as well as the operating conditions.
Three seal geometries are chosen according to the OEMs design criteria, while additional two are selected in order to extend the domain of the investigation in the case of possible axial shifting of the rotor. Axial excursion of the rotor commonly happens during the startup condition or even during the variation of the process parameters (i.e., the output power density). The validation of the correlation also for offdesign seal operations allows extending the seal dynamic coefficients predictability. The geometric parameters that define the seal geometry are those reported in Table 1. Three different rotor diameters have been investigated (220 mm, 500 mm, and 800 mm); the teeth are 20, while the rotor steps are 10. The additional two seal geometries are equal to those with the rotor diameter equal to 500 mm (geometry #4) and 800 mm (geometry #5) but with a different length of the rotor step. The three main geometries are shown in Figure 7; the steam flows from right to left.
Because CFD analysis is very timeconsuming, a factorial DACE is performed, even if not all the combinations of the parameters have been investigated. The list of the CFD analysis performed is shown in Table 2.

5. CFD Analysis and Results
The domain discretization has been carried out using Ansys™ [22]. The meshing procedure consists of creating, at a first stage, the 2D computational domain by means of quadrilateral element. After that, the 2D domain is rotated circumferentially to obtain the 3D structured mesh, composed by hexahedra elements. The 3D domain corresponds to a fivedegree sector along the circumferential direction. For instance, in the case of the seal geometry with the diameter equal to 220 mm, one element per degree has been used for the mesh rotation, while for the other diameters, two elements per degree are necessary to avoid excessive stretched elements in the circumferential direction.
It is worth mentioning that the discretization has been refined close to the wall to achieve a low Reynolds approach in the boundary layer (see Figure 8). Details concerning the discretization have been summarized in Table 3.

Among the several approaches to solve numerically the NavierStokes equations, the most employed in the turbomachinery field is the Reynolds Averaged NavierStokes (RANS) due to the good compromise between accuracy and computational time required. The IAWPS steam tables have been employed to properly model the fluid properties [23].
Turbulence has been modelled using the  model, firstly introduced by Menter [24] and widely tested in turbomachinery application in the last decade. This turbulence model is associated with the automatic near wall treatment, where a smooth transition between the low Reynolds number approach and the standard wall functions is adopted. However, due to the low value of + achieved, using the grids employed in this paper, this treatment corresponds to the low Reynolds model; hence, the viscous sublayer is not modelled but directly solved. This approach is surely preferable due to the final purpose of the present analysis, where the exact amount of flow rate for a given pressure drop must be assessed through the clearances.
The boundary condition scheme is hereafter introduced in Figure 9. The same set of conditions are imposed to all the test cases. The value of inlet total pressure is changed at each new iteration to target the desired pressure drop; moreover, the circumferential velocity component is adjusted to obtain the desired preswirl ratio.
The results of some simulations are shown in Figures 10 and 11. Static pressure and temperature contours are reported in Figure 10. Dimensionless values are obtained using the downstream value equaling the unity. Furthermore, streamlines in the meridional plane are illustrated in Figure 11, focusing on the inlet and outlet locations. The flow structure before the first clearance is different from the one occurring upstream the consecutive teeth.
(a)
(b)
(a)
(b)
6. Discharge Coefficients Correlation
The objective of the paper is the development of a new leakage correlation to accurately predict the leakage massflow rate and the pressure distribution in the seal cavities in steadystate condition (centred rotor). In a nozzle or other constrictions, the discharge coefficient is defined as the ratio of the actual massflow rate at the discharge end of the constriction to that in which an ideal working fluid expands with the same initial and exit conditions. Starting from the definition of the discharge coefficient, it is possible to define the leakage formula, which is equal to the product between the discharge coefficient and the ideal massflow rate. The ideal massflow rate in a constriction is derived from the Bernoulli equation [25], assuming that the outlet area is much smaller than the inlet one. The coefficient of discharge is defined as
Since the massflow rate, pressure, and density can be evaluated by the CFD, the discharge coefficient for each cavity can be evaluated. Using the CFD data of all the tests investigated by the DACE, the idea is to develop a new correlation based on the geometric parameters and on the operating conditions.
The authors have observed that the discharge coefficients have different values depending on the constriction type (long teeth or short teeth coupled with the rotor step); however, the values along the seal orifices almost remain constant, for both the constriction types, expect for the first tooth (see Figure 12). Since the discharge coefficient depends on the angle at which the fluid approaches to the tooth, the different discharge coefficient is likely due to the fluid, which approaches to the long teeth and to the short ones in two different manners. For instance, as can be observed in Figure 12, for test 1 shown in Table 2, the discharge coefficient in correspondence of the short teeth is about 0.3, whereas for long teeth it is practically constant and equal to 0.48. For the first tooth, the discharge coefficient is about 0.6.
After several investigations, the authors have observed that two geometrical parameters mainly influence the discharge coefficient. These parameters are listed in Table 4. The parameter is related to the ratio between the height of the constriction and the height of the seal. The other parameters are related to the axial position and length of the rotor step. In addition, one parameter regarding the operating conditions has been considered. From a primary screening of the CFD results, the authors observed that the only operating condition parameter that significantly affects the leakage is the pressure drop. The inlet preswirl and the rotor peripherical speed have practically no influence on the mass flow and therefore on the pressure distribution.

The authors assume that the crossdependence terms between the geometrical parameters and operating conditions are negligible. Thus, at a first stage, the discharge coefficients have been fitted with a polynomial as a function of the geometrical parameters for both the strips teeth and the short strips teeth combined with the rotor steps. At a second stage, the residual between the discharge coefficients evaluated using the CFD data and the discharge coefficients calculated with the polynomial in (12)(13) have been fitted using the operating parameter as shown in (14)(15). The superscripts and subscripts “” indicate the long strips (long teeth), whereas the superscripts and subscripts “s” indicate the short jstrips combined with the rotor steps.
For confidentiality reasons, the values of the coefficients and , shown in (12)–(15), are not reported in the paper.
Figures 13 and 14 show the polynomial surface that fits the discharge coefficients evaluated with the CFD analysis for both short teeth (Figure 13) and long teeth (Figure 14). The blue points represent the average value of the discharge coefficients, in correspondence of the rotor step, for each CFD analysis. The CFD points are well fitted, since the adjusted [26] is about 0.97 for the short teeth and 0.98 for the long ones, whereas the root mean square error (RMSE) [26] is equal to 0.01406 and 0.01399, respectively. These statistical parameters highlight the goodness of fit of the data.
Furthermore, Figures 15 and 16 show the fitting of the residual by the operating parameter (). The blue points represent the average value of the residual, in correspondence of the rotor step, for each CFD analysis. The adjusted is about 0.9689 and the RMSE is 0.01406.
The discharge coefficient in correspondence of the first tooth has been considered constant for all the tests and it seems to be uncorrelated from the geometric parameters and from the pressure drop.
6.1. Prediction of the Leakage and Pressure Distribution
The new discharge coefficient correlation has been implemented in the bulkflow model. The bulkflow leakage massflow estimations compared to the CFD leakage are shown in Figure 17. The bulkflow leakage is in the range of ±6% of relative error (the blue region defines the ±10% of relative error).
Moreover, the pressure distribution predicted with the bulkflow model is very similar to that evaluated with the CFD. In Figure 18(a), the pressure shows a “zigzag” behaviour along the seal, which is due to the nonconstant discharge coefficients (see Figure 18(b)).
(a)
(b)
6.2. Prediction of the Circumferential Velocity
In this section, the authors have compared the fluidwall shear stresses estimated by the bulkflow model with those computed using the CFD. For the sake of brevity, only the results related to test case 1 are reported (see Figure 19) considering both cases with and without the axial velocity component in the bulkflow model. The accuracy of the results is still long valid for the other tests.
The CFD shear stresses are averaged along the length of each cavity to be compared with the bulkflow results. In Figure 20, the circumferential velocity profile is shown. Typically, in the 1CV bulkflow model (see bulkflow models in [9, 10, 18]), the shear stresses are evaluated considering only the circumferential velocity, whereas the authors have considered also the contribution of the axial velocity (see (7)) as suggested in [19].
6.3. Numerical Verification of the New Correlation
The new correlation has been numerically verified using the results of the staggered seal shown in Figure 3, which has about twice the teeth of the DACE seal geometries and moreover the operating conditions are far from those investigated in the DACE. The relative error between the leakage evaluated with the bulkflow model and the one with the CFD is equal to 10% (see Figure 21). Moreover, the pressure distribution is accurately reproduced, and the sonic flow condition is reached as in the CFD analysis (see Figure 23). The results of Childs’ model and Wang’s model are reported in Figures 21 and 22.
7. Prediction of Staggered Labyrinth Seal Rotor Dynamic Coefficients
7.1. BulkFlow FirstOrder Mathematical Treatment
The bulkflow model, developed by the authors, is based on the onecontrol volume bulkflow model described in [13, 27–29] and here customized for staggered labyrinth seal. Each cavity of the seal is modelled by onecontrol volume. The model includes the steam properties evaluated with the IAPWS tables [23]. The model is based on the perturbation analysis, as described in [7].
The fluid dynamics within the seal are governed by the continuity equation, circumferential momentum equation, and energy equation (only in the zerothorder problem). The first two governing equations are stated as follows:(i)Continuity equation [9]:(ii)Circumferential momentum equation [9]:
7.2. Perturbation Analysis
The rotor whirls within the seal in the and directions based on the equations of motion are as follows:where is the radius of the orbit.
By using the complex notation, the perturbation of the clearance can be written as follows:
The solution of the firstorder equation represented by the perturbed terms of pressure and velocity can be given by the same mathematical expression of the perturbed clearance due to the linearity of the system as follows:where , and represent the amplitude of the harmonic solution. Four linearized algebraic equations can be obtained for each cavity as follows:
The matrices in (21) are defined in the appendix. Lastly, the matrices of each cavity can be assembled in a band matrix, where is the number of cavities. The linear system has unknowns, that is, , and , for each th cavity. Then, the complete seal linear system can be obtained as follows:
7.3. Dynamic Coefficients
By considering the contribution of the shear stresses on the rotor surface, the dynamic force acting on the rotor [18] is given by
The dynamic coefficients () of the seal are calculated by means of the following equations:where and are complex numbers that correspond to the dynamic force acting on the rotor and are produced by the forward and the backward orbit, respectively, which can be given as follows:
7.4. Numerical Results
A preliminary investigation on the dynamic behaviour of staggered labyrinth seals has been performed in the paper. The numerical results are compared with the experimental measurements performed by Kwanka and Ortinger in [11]. The radial force generated by a fourteencavity staggered labyrinth seal was obtained by integrating the pressure measured along the circumferential direction. The maximum circumferential velocity at the seal inlet is equal to 22 m/s. The seal geometry and the operating conditions are reported in Tables 5 and 6.


Figure 24 shows the measured crossforce as a function of the pressure drop, considering the rotor eccentricity equal to 0.4 mm. The geometry of the seal and the operating condition are far from those used for the development of the leakage correlation, but the prediction obtained by the authors’ model is more accurate than the numerical results reported in the paper by Kwanka and Ortinger.
8. Conclusion
The paper focuses on the numerical modelling of staggered labyrinth seals to predict the rotor dynamic coefficients. By considering the onecontrol volume bulkflow model, a new leakage correlation has been developed to estimate accurately the zerothorder solution: leakage massflow rate, static pressure distribution, circumferential velocity distribution, and temperatures.
Design and analysis of CFD experiments have been used to develop the leakage correlation using four design factors: the geometry, the pressure drop, the inlet preswirl, and the rotor peripherical speed. Only the geometry and the pressure drop have been confirmed to be effective on the estimation of the leakage and static pressure.
The procedure used by the authors for the development of the correlation relies on the definition of the discharge coefficient. The discharge coefficients have different values in correspondence of the long teeth with respect to the short ones. Polynomial functions have been used to calculate the discharge coefficients as a function of the seal geometry parameters and of the pressure drop. Once these coefficients have been evaluated, the static pressure in the cavities is estimated till the leakage converges.
Finally, a detailed description of the onecontrol volume bulkflow model for staggered labyrinth seals has been reported in the paper. The model allows the rotor dynamic coefficients to be estimated. Numerical results are compared to experimental measurements by Kwanka and Ortinger. The results of the model developed by the authors are more accurate than the other numerical models.
Appendix
Definitions of the system matrices and vector are as follows:
Definitions of the firstorder continuity and circumferential momentum equation coefficients are as follows:
Nomenclature
:  Dimensionless length of the rotor and stator of the th cavity 
:  Unsteady and steady crosssectional area of the th cavity 
:  Discharge coefficient under the th tooth 
Direct damping coefficient of the seal  
Crosscoupled damping coefficient of the seal  
:  Unsteady and steady hydraulic diameter of the th cavity 
:  Absolute roughness of the rotor and stator surface 
Steps width  
:  Lateral forces acting on the rotor 
:  Friction factors on the rotor and stator surfaces 
Steptocasing radial distance  
Steps height  
:  Unsteady and steady enthalpy of the th cavity 
:  Unsteady and perturbed clearance 
Direct stiffness coefficient of the seal  
Crosscoupled stiffness coefficient of the seal  
Length of the th cavity  
:  Unsteady and steady mass flow rate in the th orifice 
Number of long teeth  
Number of rotor steps  
Unsteady and steady pressure in the th cavity  
:  Seal inlet pressure 
:  Seal outlet pressure 
:  Polynomial coefficients for the discharge coefficients in correspondence of the short teeth combined with the rotor step 
:  Polynomial coefficients for the discharge coefficients in correspondence of the long teeth 
:  Radius of the circular orbit of the rotor 
:  Rotor radius 
:  Rotor radius in the tooth location and rotor base radius 
:  Clearance of the th cavity 
Long toothtostep axial distance  
:  Axial velocity of the th cavity 
:  Unsteady and steady tangential velocity in the th cavity 
Rotor displacement in the lateral directions  
:  Geometrical parameters for the strip 
:  Geometrical parameters for the short teeth and rotor steps 
Tooth width  
Perturbation parameter  
:  Unsteady and steady density in the th cavity 
:  Unsteady and steady rotor shear stress in the th cavity 
:  Unsteady and steady stator shear stress in the th cavity 
:  Kinetic energy carryover coefficient in the th cavity 
:  Kinetic viscosity in the th cavity 
Whirling speed of the orbit of the rotor  
Rotational speed of the rotor  
Imaginary part. 
CFD:  Computational fluid dynamics 
DACE:  Design and analysis of computer experiments 
OEMs:  Original equipment manufacturers 
SST:  Shear stress transport. 
Conflicts of Interest
The authors declare that they have no conflicts of interest.
References
 R. E. Chupp, R. C. Hendricks, S. B. Lattime, and B. M. Steinetz, “Sealing in turbomachinery,” Journal of Propulsion and Power, vol. 22, no. 2, pp. 313–349, 2006. View at: Publisher Site  Google Scholar
 F. Cangioli, P. Pennacchi, A. Vania, S. Chatterton, and P. V. Dang, “Analysis of the Dynamic Behavior of Two HighPressure Turbines for the Possible Detection of Rub Symptoms,” in Proceedings of ASME Turbo Expo 2016: Turbomachinery Technical Conference and Exposition, Seoul, South Korea. View at: Publisher Site  Google Scholar
 F. Mahler, “Advanced Seal Technology,” Pratt & Whitney Aircraft Report PWA−4372, 1972. View at: Google Scholar
 H. A. Tanvir, Evaluation of steam turbines triangular tooth on stator labyrinth seal [Master, thesis], Texas A & M University, 2012.
 H. J. Thomas, “Unstable oscillations of turbine rotors due to steam leakage in the sealing glands and the buckets,” Bulletin Scientifique A J M, vol. 71, no. 11, pp. 223–236, 1958. View at: Google Scholar
 J. S. Alford, “Protection of labyrinth seals from flexural vibration,” Journal of Engineering for Gas Turbines and Power, vol. 86, no. 2, pp. 141–147, 1964. View at: Publisher Site  Google Scholar
 D. Childs, Turbomachinery Rotordynamics  Phenomena, Modeling, and Analysis, John Wiley & Sons, New York, NY, USA, 1993.
 T. Iwatsubo, “Evaluation of instability forces of labyrinth seals in turbines or compressors,” in NASA CP 2133, Proceedings of a Workshop at Texas A & M University,Rotordynamic Instability Problems in High Performance Turbomachinery, pp. 205–222, 1980. View at: Google Scholar
 D. W. Childs and J. K. Scharrer, “An iwatsubobased solution for labyrinth seals: Comparison to experimental results,” Journal of Engineering for Gas Turbines and Power, vol. 108, no. 2, pp. 325–331, 1986. View at: Publisher Site  Google Scholar
 W. Wang, Y. Liu, and P. Jiang, “Numerical investigation on influence of real gas properties on nonlinear behavior of labyrinth sealrotor system,” Applied Mathematics and Computation, vol. 263, pp. 12–24, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 K. Kwanka and W. Ortinger, “Rotordynamic Coefficients of Long Staggered Labyrinth Gas Seals,” International Journal of Rotating Machinery, vol. 1, no. 34, pp. 285–291, 1995. View at: Publisher Site  Google Scholar
 Z. Li, J. Li, X. Yan, and Z. Feng, “Effects of pressure ratio and rotational speed on leakage flow and cavity pressure in the staggered labyrinth seal,” Journal of Engineering for Gas Turbines and Power, vol. 133, no. 11, Article ID 114503, 2011. View at: Publisher Site  Google Scholar
 F. Cangioli, P. Pennacchi, G. Vannini, and L. Ciuchicchi, “Effect of energy equation in one controlvolume bulkflow model for the prediction of labyrinth seal dynamic coefficients,” Mechanical Systems and Signal Processing, vol. 98, pp. 594–612, 2018. View at: Publisher Site  Google Scholar
 K. Neumann, “Zur Frage der Verwendung von Durchblicktungen im Dampgturbinebau,” Maschinentechnik, vol. 13, no. 4, 1964. View at: Google Scholar
 P. K. Swamee and A. K. Jain, “Explicit equations for pipeflow problems,” Journal of the Hydraulics Division, vol. 102, no. 5, pp. 657–664, 1976. View at: Google Scholar
 D. W. Childs and J. K. Scharrer, “Theory Versus Experiment for the Rotordynamic Coefficient of Labyrinth Gas Seals: Part II—A Comparison to Experiment,” Journal of Vibration, Acoustics, Stress, and Reliability in Design, vol. 110, no. 3, p. 281, 1988. View at: Publisher Site  Google Scholar
 D. Eser and Y. Dereli, “Comparisons of rotordynamic coefficients in stepped labyrinth seals by using ColebrookWhite friction factor model,” Meccanica, vol. 42, no. 2, pp. 177–186, 2007. View at: Publisher Site  Google Scholar
 D. Eser and J. Y. Kazakia, “Air flow in cavities of labyrinth seals,” International Journal of Engineering Science, vol. 33, no. 15, pp. 2309–2326, 1995. View at: Publisher Site  Google Scholar
 H. R. Wyssmann, T. C. Pham, and R. J. Jenny, “Prediction of stiffness and damping coefficients for centrifugal compressor labyrinth seals,” Journal of Engineering for Gas Turbines and Power, vol. 106, no. 4, pp. 920–926, 1984. View at: Publisher Site  Google Scholar
 J. Kiijarvi, “Darcy friction factor formulae in turbulent pipe flow,” Lunowa Fluid Mechanics Paper, Article ID 110727, 2011, no. 110727. View at: Google Scholar
 B. Hodkinson, “Estimation of the Leakage through a Labyrinth Gland,” Proceedings of the Institution of Mechanical Engineers, vol. 141, no. 1, pp. 283–288, 2016. View at: Publisher Site  Google Scholar
 Ansys.com, “Ansys,” http://www.ansys.com/. View at: Google Scholar
 W. Wagner, “The IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam,” Journal of Engineering for Gas Turbines and Power, vol. 122, no. 1, pp. 150–184, 2000. View at: Google Scholar
 F. R. Menter, “Twoequation eddyviscosity turbulence models for engineering applications,” AIAA Journal, vol. 32, no. 8, pp. 1598–1605, 1994. View at: Publisher Site  Google Scholar
 F. White, Fluid Mechanics, McGraw Hill, 7th edition, 2009.
 J. Rice, “Mathematical statistics and data analysis,” Mathematical statistics and data analysis, 2006. View at: Google Scholar
 F. Cangioli, S. Chatterton, P. Pennacchi, L. Nettis, and L. Ciuchicchi, “Thermoelasto bulkflow model for labyrinth seals in steam turbines,” Tribology International, vol. 119, pp. 359–371, 2018. View at: Publisher Site  Google Scholar
 F. Cangioli, P. Pennacchi, G. Vannini et al., “On the Thermodynamic Process in the BulkFlow Model for the Estimation of the Dynamic Coefficients of Labyrinth Seals,” Journal of Engineering for Gas Turbines and Power, vol. 140, no. 3, Article ID 032502, 2018. View at: Publisher Site  Google Scholar
 F. Cangioli, P. Pennacchi, G. Riboni et al., “Sensitivity analysis of the onecontrol volume bulkflow model for a 14 teethonstator straightthrough labyrinth seal,” in Proceedings of the ASME Turbo Expo 2017: Turbomachinery Technical Conference and Exposition, GT 2017, Charlotte, NC, USA, June 2017. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Filippo Cangioli et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.