Abstract
Dynamic property of a gas face seal is about an internal ability of the massspringdamper system consisting of the nonlinear gas film, the flexible support, and the component mass or moment of inertia to handle ambient disturbances. Stability and tracking properties are important dynamic components, the current research for which mainly focuses on highspeed and highpressure conditions in turbomachinery applications. In this paper, perturbation theory is used to obtain the linearized properties of the gas film in a spiral groove gas face seal, where the BoltzmannReynolds model with FukuiKaneko approximation is involved to account for slip flow under lowparameter conditions. Transmissibility of the step jump method is introduced into the perturbation method, and a critical ambientdisturbance amplitude is further proposed to cover the defect of transmissibility in characterizing the tracking property. The effects of slip flow on the stability and tracking properties of gas face seals are explored. The results show that slip flow induces an extra “relaxation region” in the critical transverse moment of inertia, providing a softer demand for angular stability. With regard to transmissibility, slip flow will induce an advantageous effect on transmitting the rotor motion to the stator. However, after considering the space required for tracking motions, the results of critical amplitude indicate that slip flow will induce an extra “rigorous region”, leading gas face seals to require more conservative design and more rigorous manufacture and assembly for noncontacting tracking.
1. Introduction
A spiral groove gas face seal (Fig. 1) is a type of noncontacting face seals that prevents sealed gas from escaping from one compartment to another. A well designed gas face seal should have an appropriate dynamic property to head off excessive leakage and face contact [1, 2]. As depicted in Fig. 1, the dynamic system of a gas face seal is a typical massspringdamper system which consists of the nonlinear gas film stiffness and damping, the flexible support, and the mass and the transverse moment of inertia of the stator. Thus, the dynamic property of a gas face seal is about a systematic internal ability to handle ambient disturbances, which typically include axial endplay and angular wobble. The axial endplay of the rotor is caused by the axial endplay of the shaft, and the angular wobble of the mated rings is a result of manufacturing and assembly tolerances such as the runouts of both rings and the imperfection in the stator’s flexible support. When subjected to these disturbances, the stator may become completely unstable. Even in a stable state, an improper tracking property can still either induce excessive leakage or local face contact [3, 4]. It is therefore imperative to study the stability and tracking properties of gas face seals which have played important dynamic roles in other applications [5, 6].
Dynamic research on gas face seals has been performed using both direct numerical methods [711] and semianalytical methods [1217]. In such research, the stability and tracking properties of gas face seals act as important focal points. Green and Barnsby [8] used the statespaceform method developed by Miller and Green [9] to gauge the influence of sealing parameters on the stability threshold of coned gas face seals. The tracking property of the stator was also discussed by a parametric investigation [10]. Zirkelback and San Andres obtained the axial stiffness and damping coefficients of the gas film in a spiral groove gas face seal using a perturbation method [12], and subsequently presented a corresponding stability analysis [13]. Ruan [14] added angular motions to the dynamics to further develop the perturbation method, and presented a simple stability analysis. Meanwhile, Miller and Green developed a step jump method [15, 16] to provide a stability and tracking analysis of spiral groove gas face seals [17]. Because gas face seals are usually used in highspeed and highpressure turbomachinery, all of the above works have been performed under highparameter conditions.
Fig. 1Schematic of a typical spiral groove face seal and equivalent massspringdamper representation
With the continuing development of gas face seals, the applications of these seals are expanding to devices such as mixers and agitators, where the speed and pressure represent comparatively lowparameter conditions [18]. Such lowparameter conditions can even exist in highparameter turbomachinery in scenarios such as debugging and engine warm up. In such applications, the continuumflow assumption of the gas film breaks down and the rarefaction effect becomes significant. First [19], second [20], and 1.5order [21] slip flow models have been developed to modify the classical continuumflow Reynolds equation, accounting for rarefaction effects by allowing molecular slipping to occur at the gaswall interface. Fukui and Kaneko developed a BoltzmannReynolds model for an arbitrary Knudsen number ${K}_{n}$ based on the linearized Boltzmann equation and provided a database of values in terms of the inverse Knudsen number $D$ [22]. Subsequently, they published a set of polynomials for the curve fitting of the above database [23]. These slip flow models have been used in the analysis of gas face seals. Pecht and Netzel [24] and Ruan [18] gauged the effect of slip flow on the steadystate performance of gas face seals including liftoff speed, leakage rate, frequencyindependent gas film stiffness and load carrying capacity. Ruan [25] also considered slip flow in a study on the transient performance of spiral groove gas face seals during the startup and shutdown operations. However, few previous works have researched the stability and tracking properties of gas face seals under lowparameter conditions.
The objective of the present study is to provide an investigation of the effect of slip flow on the stability and tracking of spiral groove gas face seals under lowparameter conditions. In addition, the approach used to analyze stability and tracking properties can still be improved. In direct numerical methods [711], the stability and tracking properties are analyzed numerically, but this approach is inefficient for intensive parametric analyses. In the perturbation method [1214], it is convenient to obtain expressions for the critical stability condition because of the existence of linearized gas film coefficients in the time domain. Yet, the tracking property is still characterized by axial relative displacement or angular relative tilt angle. In the step jump method [17], the transmissibility is presented in the form of an expression that can be used to characterize the tracking property. However, the step response, which acts as the constitutive model of gas film, is easy to analyze in the frequency domain [16] but difficult to analyze in the time domain. Therefore, a balanced approach is proposed that introduces the transmissibility into the perturbation method. Furthermore, because the transmissibility only reflects the tracking ability while ignores the space required for tracking motions, critical ambientdisturbance amplitude, above which the seal will achieve face contact, is provided based on the transmissibility to cover the defect of transmissibility in characterizing the tracking property.
2. Model
2.1. Gas film lubrication and sealing dynamics
A schematic diagram of a typical gas face seal is shown in Fig. 1. A dynamic schematic of the mated rings and the rotor geometry are depicted in Fig. 2. A total of $N$ spiral grooves are processed on the rotor face at a depth of ${\delta}_{g}$. $\alpha $ is the spiral angle and varies from 0 to 180°. The groove width ${W}_{g}$ relative to the land width ${W}_{l}$ is measured by $\lambda ={W}_{g}/({W}_{g}+{W}_{l})$, and the groove length relative to the dam length is measured by $\beta =({r}_{g}{r}_{i})/({r}_{o}{r}_{i})$. Here, ${r}_{o}$, ${r}_{i}$ and ${r}_{g}$ are the dimensionless outer radius, inner radius and boundary of the groove region and dam region, respectively, where ${r}_{o}={r}_{o}^{*}/{r}_{i}^{*}$, ${r}_{i}=1$, and ${r}_{g}={r}_{g}^{*}/{r}_{i}^{*}$. Here, superscript $*$ represents the dimensional form of a variable in the present study.
The gas pressure distribution ${p}^{*}$ in the gas face seal is governed by the lubrication equation. Assuming that the gas flow is ideal and isothermal, the compressible Reynolds equation is then given by [811, 1417, 25]:
where $\mu $ is the gas viscosity, ${r}^{*}$ is the radius, ${\mathbf{i}}_{\theta}$ is the unit vector in the $\theta $ direction, ${h}^{*}$ is the gas film thickness, $\omega $ is the shaft speed, and ${t}^{*}$ is the time. The dimensionless form of Eq. (1) and its corresponding boundary conditions are given by:
where $r={r}^{*}/{r}_{i}^{*}\text{,}$$h={h}^{*}/{\delta}_{g}\text{,}$$p={p}^{*}/{p}_{i}^{*}\text{,}$$t=\omega {t}^{*}$ and the compressibility number $\mathrm{\Lambda}=6\mu {r}_{i}^{*2}\omega /{p}_{i}^{*}{\delta}_{g}^{2}$. Also, ${p}_{i}^{*}$ and ${p}_{o}^{*}$ are the pressures at the inner and outer radii. Here, a flow factor ${Q}_{p}$ is incorporated into Eq. (2a) to modify the compressible Reynolds equation to account for the slip flow [18, 22, 23]:
where ${R}_{gas}$ is the gas constant, and ${T}_{a}$ is the ambient temperature.
The dimensionless rotor tilt angle ${\mathbf{\gamma}}_{r}$, the stator tilt angle ${\mathbf{\gamma}}_{s}$ and the relative tilt angle ${\mathbf{\gamma}}_{rel}$ can be decomposed into components along the $X$ and $Y$ axes:
where ${\mathbf{\gamma}}_{r}={\mathbf{\gamma}}_{r}^{*}{r}_{i}^{*}/{\delta}_{g}$, ${\mathbf{\gamma}}_{s}={\mathbf{\gamma}}_{s}^{*}{r}_{i}^{*}/{\delta}_{g}$ and ${\mathbf{\gamma}}_{rel}={\mathrm{\gamma}}_{rel}^{*}{r}_{i}^{*}/{\delta}_{g}$. The gas film thickness $h$ at any point ($x$, $y$) on the stator face is given by:
where ${c}_{0}$ is the dimensionless central film thickness at equilibrium, $z$ is the dimensionless axial displacement of the stator from ${c}_{0}$, ${z}_{r}$ is the dimensionless axial endplay of the rotor, and $\u23291\u232a$ indicates a dimensionless groove depth in the groove region. Here, ${c}_{0}={c}_{0}^{*}/{\delta}_{g}$, $z={z}^{*}/{\delta}_{g}$ and ${z}_{r}={z}_{r}^{*}/{\delta}_{g}$. By considering the stiffness and damping properties of the springs and the secondary seal, the dimensionless dynamic equations for the axial and angular modes are given by [9, 14, 17, 25]:
In Eq. (6), $m$ is the dimensionless mass, and $I$ is the dimensionless transverse moment of inertia of the stator, where $m={m}^{*}{\delta}_{g}{\omega}^{2}/{p}_{i}^{*}{r}_{i}^{*2}$ and$I={I}^{*}{\delta}_{g}{\omega}^{2}/{p}_{i}^{*}{r}_{i}^{*4}$. ${k}_{sz}$ is the dimensionless axial stiffness of the springs and ${c}_{sz}$ is the dimensionless axial damping of the flexible support, where ${k}_{sz}={k}_{sz}^{*}{\delta}_{g}/{p}_{i}^{*}{r}_{i}^{*2}$ and ${c}_{sz}={c}_{sz}^{*}{\delta}_{g}\omega /{p}_{i}^{*}{r}_{i}^{*2}$. According to the work of Green and Etsion [26], the dimensionless angular stiffness ${k}_{s\gamma}$ and the damping ${c}_{s\gamma}$ are given by:
assuming that the support forces act at the outer radius. ${F}_{closing}$ is the total load of the fluid pressure and the spring support load at equilibrium, which is applied to the rear of the stator to balance the gas film load ${F}_{g}$ caused by the fluid effect. ${M}_{gx}$ and ${M}_{gy}$ are the gas film moments about the $X$ and $Y$ axes. ${M}_{si}$ is a moment caused by the dimensionless static stator tilt ${\gamma}_{si}$, which is arbitrarily assumed to be oriented about the $X$ axis. The constant moment is given by ${M}_{si}={k}_{s\gamma}{\gamma}_{si}$. ${F}_{g}$, ${M}_{gx}$ and ${M}_{gy}$ can all be obtained by integrating $p$ over the dimensionless sealing area ${A}_{n}$, as follows:
Fig. 2Dynamic schematic of a spiral groove gas face seal and the geometry of the rotor
2.2. Perturbation method
Malanoski and Pan [27] performed a pioneering work in studying the axial gas film coefficients in a spiral groove gas thrust bearing by a perturbation method. Zirkelback and San Andres [12] used the perturbation method to analyze a spiral groove gas face seal. Ruan [14] added angular modes to further develop the perturbation method. In the present study, the perturbation method described in Ref. [14] is adopted. Because our dimensionless approach differs from that of Ref. [14], and because the flow factor ${Q}_{p}$ is involved, the whole process should be deduced again. If small oscillatory motions ($\u2206z$,$\u2206{\gamma}_{x}$,$\u2206{\gamma}_{y}$) with excitation frequency ${\omega}_{e}$ occur near the equilibrium ${c}_{0}$, Eq. (5) can then be rewritten as:
where, the excitation frequency ratio $\mathrm{\Omega}={\omega}_{e}/\omega $ is introduced because of the dimensionless process for ${e}^{i{\omega}_{e}{t}^{*}}$. These small motions cause a small pressure perturbation $\u2206p$ about the equilibrium pressure ${p}_{0}$. Using a firstorder Taylor series expansion, the new pressure can be expressed as:
After substituting Eqs. (9, 10) into Eq. (3a) and only remaining terms of order O(1), O($\u2206z$), O($\u2206{\gamma}_{x}$), O($\u2206{\gamma}_{y}$), O($\u2206\dot{z}$), O($\u2206\dot{{\gamma}_{x}}$), O($\u2206\dot{{\gamma}_{y}}$), the following seven equations can be produced:
$=\mathrm{\Lambda}\frac{\partial ({c}_{0}{p}_{z}+{p}_{0})}{\partial \theta}2\mathrm{\Lambda}{\mathrm{\Omega}}^{2}{c}_{0}\stackrel{}{{p}_{z}},$
$=\mathrm{\Lambda}\frac{\partial ({c}_{0}{p}_{x}+y{p}_{0})}{\partial \theta}2\mathrm{\Lambda}{\mathrm{\Omega}}^{2}{c}_{0}\stackrel{}{{p}_{x}},$
$=\mathrm{\Lambda}\frac{\partial ({c}_{0}{p}_{y}x{p}_{0})}{\partial \theta}2\mathrm{\Lambda}{\mathrm{\Omega}}^{2}{c}_{0}\stackrel{}{{p}_{y}},$
The corresponding boundary conditions are given by Eq. (2b) and:
These seven equations can be classified into four groups. Eq. (11) governs the steadystate pressure distribution ${p}_{0}$ at equilibrium ${c}_{0}$. Eqs. (1214) are three pairs of perturbation Reynolds equations that govern three pairs of variables (${p}_{z}$, $\stackrel{}{{p}_{z}}$), (${p}_{x}$, $\stackrel{}{{p}_{x}}$) and (${p}_{y}$, $\stackrel{}{{p}_{y}}$), respectively. For each pair, an iterative process is required to produce stable variables.
Based on Eq. (8), the generalized force in response to the perturbation can be expressed as:
$+{\int}_{{A}_{n}}^{}\left\{\begin{array}{ccc}\stackrel{}{{p}_{z}}& \stackrel{}{{p}_{x}}& \stackrel{}{{p}_{y}}\\ y\stackrel{}{{p}_{z}}& y\stackrel{}{{p}_{x}}& y\stackrel{}{{p}_{y}}\\ x\stackrel{}{{p}_{z}}& x\stackrel{}{{p}_{x}}& x\stackrel{}{{p}_{y}}\end{array}\right\}dA\left\{\begin{array}{c}\u2206\dot{z}\\ \u2206\dot{{\gamma}_{x}}\\ \u2206\dot{{\gamma}_{y}}\end{array}\right\}.$
Therefore, the stiffness and damping matrices are:
$\left[\begin{array}{ccc}{c}_{zz}& {c}_{zx}& {c}_{zy}\\ {c}_{xz}& {c}_{xx}& {c}_{xy}\\ {c}_{yz}& {c}_{yx}& {c}_{yy}\end{array}\right]={\int}_{{A}_{n}}^{}\left\{\begin{array}{ccc}\stackrel{}{{p}_{z}}& \stackrel{}{{p}_{x}}& \stackrel{}{{p}_{y}}\\ y\stackrel{}{{p}_{z}}& y\stackrel{}{{p}_{x}}& y\stackrel{}{{p}_{y}}\\ x\stackrel{}{{p}_{z}}& x\stackrel{}{{p}_{x}}& x\stackrel{}{{p}_{y}}\end{array}\right\}dA.$
Because the works of Ruan [14] and Miller and Green [16, 17] indicate that the crosscoupling coefficients can be ignored and that only three terms are unique, the reduced matrices are given by:
2.3. Stability and tracking
2.3.1. Stability
The dimensionless form of the linearized dynamics used for the stability research is given by:
After substituting a small axial displacement $z={z}_{\epsilon}{e}^{{s}_{e}t/\omega}$ into the equation of axial motion, the characteristic equation can then be obtained as:
where the complex number ${s}_{e}={\lambda}_{e}+{\omega}_{e}i$. For the critical condition of axial stability, the real part ${\lambda}_{e}$ is zero and Eq. (19) can then be reduced to:
Here, subscript $cr$ represents the critical condition in the present study. To satisfy Eq. (20), both the real and imaginary parts should be zero. Therefore, the critical excitation frequency ratio and the critical mass for axial stability should satisfy:
Here, an iterative process is inevitable to calculate ${\mathrm{\Omega}}_{cr}$. It is because before ${\mathrm{\Omega}}_{cr}$ can be solved by Eq. (21a), ${c}_{zz}$ is required. However, ${c}_{zz}$ is related to ${\mathrm{\Omega}}_{cr}$ in the perturbation method.
The same method is also used to analyze the angular stability. After substituting small angular displacements ${\gamma}_{x}={\gamma}_{x\epsilon}{e}^{{s}_{e}t/\omega}$ and ${\gamma}_{y}={\gamma}_{y\epsilon}{e}^{{s}_{e}t/\omega}$ into the equations of angular motions, the dimensionless angular dynamics are transformed into:
Because the determinant of the coefficient should be zero, the characteristic equation is then obtained as:
${A}_{2}=2\left(\frac{I}{{\omega}^{2}}\right)\left({k}_{xx}+{k}_{s\gamma}\right)+\frac{{\left({c}_{xx}+{c}_{s\gamma}\right)}^{2}+{c}_{xy}^{2}}{{\omega}^{2}},$
${A}_{3}=2\frac{\left({c}_{xx}+{c}_{s\gamma}\right)\left({k}_{xx}+{k}_{s\gamma}\right)+{k}_{xy}{c}_{xy}}{\omega},{A}_{4}={({k}_{xx}+{k}_{s\gamma})}^{2}+{k}_{xy}^{2}.$
By considering the real part of ${s}_{e}$ to be zero under the critical condition, the critical excitation frequency ratio and the critical transverse moment of the inertia for angular stability can be obtained from Eq. (23), as follows:
where an iterative process is also necessary to calculate of ${\mathrm{\Omega}}_{cr}$.
2.3.2. Tracking
When the stator satisfies the stability demand, an appropriate tracking property is expected to avoid both excessive leakage and local face contact. Because the typical ambient disturbances are the axial endplay and the angular wobble of the rotor, where ${M}_{si}$ in Eq. (6) is ignored, ${\omega}_{e}$ is equal to $\omega $ (i.e., $\mathrm{\Omega}=$ 1) in the tracking study. The dimensionless sealing dynamics for the axial endplay and the angular wobble can then be expressed as:
A Laplace transform approach [14, 17] is then used to solve Eq. (25). The axial and angular motions in the Laplace domain are given by:
${b}_{1}=\left({c}_{xx}s+{k}_{xx}\right){\mathrm{\Gamma}}_{rx}\left(s\right)+\left({c}_{xy}s+{k}_{xy}\right){\mathrm{\Gamma}}_{ry}\left(s\right){c}_{xx}{\gamma}_{rx}\left(0\right)$
${c}_{xy}{\gamma}_{ry}\left(0\right)+I\dot{{\gamma}_{x}}\left(0\right)+\left(Is+{c}_{xx}+{c}_{s\gamma}\right){\gamma}_{x}\left(0\right)+{c}_{xy}{\gamma}_{y}\left(0\right),$
${b}_{2}=\left({c}_{xy}s+{k}_{xy}\right){\mathrm{\Gamma}}_{rx}\left(s\right)+\left({c}_{xx}s+{k}_{xx}\right){\mathrm{\Gamma}}_{ry}\left(s\right)+{c}_{xy}{\gamma}_{rx}\left(0\right)$
${c}_{xy}{\gamma}_{ry}\left(0\right)+I\dot{{\gamma}_{y}}\left(0\right){c}_{xy}{\gamma}_{x}\left(0\right)+\left(Is+{c}_{xx}+{c}_{s\gamma}\right){\gamma}_{y}\left(0\right),$
where $s$ is the Laplace variable. $Z\left(s\right)$ and ${Z}_{r}\left(s\right)$ are the Laplace transforms of $z\left(t\right)$ and ${z}_{r}\left(t\right)$, respectively. ${\mathrm{\Gamma}}_{x}\left(s\right),$${\mathrm{\Gamma}}_{y}\left(s\right),$${\mathrm{\Gamma}}_{rx}\left(s\right)\mathrm{}$and ${\mathrm{\Gamma}}_{ry}\left(s\right)$ are the Laplace transforms of ${\gamma}_{x}\left(t\right)$, ${\gamma}_{y}\left(t\right)$, ${\gamma}_{rx}\left(t\right)\mathrm{}$and ${\gamma}_{ry}\left(t\right)$, respectively. The variables with zero in parentheses are the initial values at $t=$ 0. If the axial relative displacement or the angular relative tile angle is used to characterize the tracking property [14], the inverse Laplace transform of Eq. (26) can be used to obtain $z\left(t\right)$, ${\gamma}_{x}\left(t\right)$ and ${\gamma}_{y}\left(t\right)$ numerically or analytically. Here, the transmissibility in the step jump method [17] is introduced into the perturbation method to derive expressions for the characterization of the tracking property.
In the axial mode, the dimensional axial endplay and the corresponding response are assumed to be ${z}_{r}^{*}\left({t}^{*}\right)={A}_{rz}^{*}{e}^{i\omega {t}^{*}}$ and ${z}^{*}\left({t}^{*}\right)={A}_{z}^{*}{e}^{i(\omega {t}^{*}{\mathrm{\Psi}}_{z})}\text{,}$ and their dimensionless forms are ${z}_{r}\left(t\right)={A}_{rz}{e}^{it}$ and $z\left(t\right)={A}_{z}{e}^{i(t{\mathrm{\Psi}}_{z})}$ respectively, where ${\mathrm{\Psi}}_{z}$ is the axial phase difference, ${A}_{rz}={A}_{rz}^{*}/{\delta}_{g}$ and ${A}_{z}={A}_{z}^{*}/{\delta}_{g}$. By substituting ${z}_{r}\left(t\right)$ and $z\left(t\right)$ into the axial motion of Eq. (25), one can obtain:
where $\stackrel{}{z}$ is the axial transmissibility, and ${\stackrel{}{z}}_{rel}$ is the axial relative transmissibility that is demanded to be less than 0.3 in a proper tracking design [28].
In the angular modes, the dimensional angular wobble and the corresponding response are assumed to be ${\mathbf{\gamma}}_{r}^{*}\left({t}^{*}\right)={A}_{r\gamma}^{*}{e}^{i\omega {t}^{*}}$ and ${\mathbf{\gamma}}_{s}^{*}\left({t}^{*}\right)={A}_{s\gamma}^{*}{e}^{i(\omega {t}^{*}{\mathrm{\Psi}}_{\gamma})}={\stackrel{}{\gamma}}_{s}{A}_{r\gamma}^{*}{e}^{i(\omega {t}^{*}{\mathrm{\Psi}}_{\gamma})}$, and their dimensionless forms are ${\mathbf{\gamma}}_{r}\left(t\right)={A}_{r\gamma}{e}^{it}$ and ${\mathbf{\gamma}}_{s}\left(t\right)={\stackrel{}{\gamma}}_{s}{A}_{r\gamma}{e}^{i(t{\mathrm{\Psi}}_{\gamma})}$ respectively, where ${\stackrel{}{\gamma}}_{s}$ is the angular transmissibility, ${\Psi}_{\gamma}$ is the angular phase difference, and ${A}_{r\gamma}={A}_{r\gamma}^{*}{r}_{i}^{*}/{\delta}_{g}$. A relation similar to Eq. (27) is more difficult to be obtained in the angular modes. The steadystate response to the harmonic angular wobble takes the form of harmonic functions, as follows:
${\gamma}_{y}\left(t\right)={A}_{r\gamma}{A}_{y,c}\mathrm{cos}t+{A}_{r\gamma}{A}_{y,s}\mathrm{sin}t,$
where ${A}_{x,c}$, ${A}_{x,s}$, ${A}_{y,c}$ and ${A}_{y,s}$ are the angular amplitude coefficients. By transforming Eq. (28) into the Laplace domain and then substituting the results into Eq. (26b26c), these angular amplitude coefficients are gathered as:
${A}_{x,s}=\frac{1}{\mathrm{\Delta}}\left({a}_{22}\left({c}_{xy}s+{k}_{xy}\right){a}_{12}({c}_{xx}s+{k}_{xx})\right),$
${A}_{y,c}=\frac{1}{\mathrm{\Delta}}\left({a}_{21}\left({c}_{xx}s+{k}_{xx}\right){a}_{11}({c}_{xy}s+{k}_{xy})\right),$
${A}_{y,s}=\frac{1}{\mathrm{\Delta}}\left({a}_{21}\left({c}_{xy}s+{k}_{xy}\right)+{a}_{11}({c}_{xx}s+{k}_{xx})\right).$
The real part of ${\gamma}_{x}$ corresponds to the tilt angle about the $X$ axis, and the imaginary part of ${\gamma}_{x}$ corresponds to the tilt angle about the $Y$ axis. Likewise, the real part of ${\gamma}_{y}$ corresponds to the tilt angle about the $Y$ axis, and the imaginary part of ${\gamma}_{y}$ corresponds to the negative tilt angle about the $X$ axis. Therefore, ${\gamma}_{x}$ and ${\gamma}_{y}$ satisfy the following:
If ${\gamma}_{x}$ and ${\gamma}_{y}$ are expressed as:
the angular transmissibility ${\stackrel{}{\gamma}}_{x}$ and the phase difference ${\mathrm{\Psi}}_{rx}$ in the $X$ mode can be given by:
Because the transverse moments of inertia and the support properties are symmetrical about the $X$ and $Y$ axes,${\stackrel{}{\gamma}}_{x}={\stackrel{}{\gamma}}_{y}={\stackrel{}{\gamma}}_{s}$ and ${\mathrm{\Psi}}_{\gamma x}={\mathrm{\Psi}}_{\gamma y}={\mathrm{\Psi}}_{\gamma}$. The angular relative tilt angle ${\mathbf{\gamma}}_{rel}$ and the relative transmissibility${\stackrel{}{\gamma}}_{rel}$ can then be given by:
The transmissibility can be used to characterize the tracking property. However, it is a ratio to only convey the information about the amplitude rather than the datum line. It is possible to obtain the tracking ability in detail from the transmissibility but it is impossible to see how much space is reserved for the tracking motions. A feasible approach to observe the tracking property including both “amplitude” and “datum line” is to investigate how much the amplitude of axial endplay and angular wobble can be allowed before face contact occurs. Therefore, the critical ambientdisturbance amplitude is provided based on the transmissibility to cover the defect of transmissibility in characterizing the tracking property. Because the linearized coefficients of the gas film are based on small perturbations from the parallel face position, the following results may be not extremely accurate when subject to large deviations. Yet, the validity of this method can be acceptable before a more elaborate analysis becomes available [26]. Even if a more elaborate analysis can be applied, the time cost should be considered. Therefore, the critical condition for face contact is given by:
where the surface texture is not considered. The equation then can be reduced to an expression that contains the transmissibility:
3. Results and discussion
Table 1 lists the parameters of a spiral groove gas face seal. The subsequent simulations under the slip flow and nonslip flow conditions are both performed based on the data in this table. Note that the value of ${c}_{sz}^{*}$ is not measured but selected from Refs. [14, 25] (1 kN∙s/m). It is not a significant factor because the present study focuses on the simulations and is not related to any experimental work. Compressibility number $\Lambda $ ranging from 2.5233 to 5.8878 (i.e., shaft speed $\omega $ ranges from 300 to 700 r/min) is selected as the representative lowparameter conditions in the simulation. The calculation is performed by using a finite element method. Numerical accuracy and stability are important for a complex dynamic system. A very fine mesh would result in excessive computation time. Therefore, a balance between accuracy and economical computing must be maintained. After the experiments, it is determined that a 3160finite element mesh is adopted. The calculation is carried out on a 2.3 GHz workstation.
Table 1Parameters of a spiral groove gas face seal
Parameter  Value 
Outer radius ${r}_{o}^{*}$  61.7 mm 
Inner radius ${r}_{i}^{*}$  51.6 mm 
Balance radius ${r}_{g}^{*}$  53.1 mm 
Spiral angle $\alpha $  15º 
Number of grooves $N$  12 
Groove to land width ratio $\lambda $  0.5 
Groove to dam length ratio$\beta $  0.6 
Groove depth ${\delta}_{g}$  6 μm 
Stator mass ${m}^{*}$  0.13 kg 
Transverse moment of inertia ${I}^{*}$  2.5×10^{}^{4} kg∙m^{2} 
Pressure at outer radius ${p}_{o}^{*}$  0.2 MPa 
Pressure at inner radius ${p}_{i}^{*}$  0.1 MPa 
Support axial stiffness ${k}_{sz}^{*}$  16.4 kN/m 
Support axial damping ${c}_{sz}^{*}$  1 kN∙s/m 
Gas constant ${R}_{gas}$  300∙K 
Ambient temperature ${T}_{a}$  287 J/(kg∙K) 
3.1. Fundamental sealing performance
The effect of slip flow on the fundamental sealperformance variables, including the central film thickness at equilibrium, the liftoff speed and the stiffness and damping of the gas film, are analyzed to form the foundation of the following discussion on stability and tracking. Fig. 3 shows the central film thicknesses at equilibrium at different shaft speeds. It can be seen that the central film thickness at equilibrium increases as the shaft speed increases. It is because the closing force ${F}_{closing}$ is constant in these cases and the hydrodynamic effect is enhanced by the increased shaft speed. A greater central film thickness at equilibrium helps to weaken the hydrodynamic effect to balance the constant closing force. A significant deviation between the two flow approaches is apparent. The central film thickness at equilibrium is smaller under the slip flow condition than under the nonslip flow condition. The relative error decreases as the shaft speed increases because the slip flow is diminishing, which can be visualized from the variation of the minimum inverse Knudsen number ${D}_{min}$ from 11.985 at 300 r/min to 22.767 at 700 r/min. A decreased central film thickness at equilibrium caused by the slip flow means that the slip flow will induce a reduced load carrying capacity because of the loss of viscous pumping due to the slippage at the solid boundary. Fig. 4 displays the liftoff speeds at different liftoff thicknesses of 0.6, 0.7, 0.8, 0.9 and 1.0 μm. The liftoff speed increases with the increased liftoff thickness because of the decreased hydrodynamic effect. The liftoff speed is higher under the slip flow condition, and the relative error decreases with the increased liftoff thickness. A higher liftoff speed under the slip flow condition means that if a gas face seal were to start up from rest, its corresponding liftoff would be delayed, and the face of the seal would suffer more severe wear. Therefore, the slip flow should be noted during the startup and shutdown operations accounting for the service life of seals.
Fig. 3Comparison of central film thicknesses at equilibrium c0 at different shaft speeds ω from 300 to 700 r/min (Λ from 2.5233 to 5.8878)
Fig. 4Comparison of liftoff speeds at different liftoff film thicknesses from 0.6 to 1.0 μm (dimensionless values from 0.1000 to 0.1667)
Fig. 5Comparison of frequency dependent stiffness and damping coefficients of gas film. a) Λ varies from 2.5233 to 5.8878 (ω from 300 to 700 r/min) where Ω= 1; b) Ω varies from 0.5 to 50 where Λ= 2.5233 (ω= 300 r/min)
a)
b)
Fig. 5 shows the linearized stiffness and damping coefficients of the gas film by using the perturbation method. In Fig. 5(a), $\mathrm{\Lambda}$ is varied from 2.5233 to 5.8878 (i.e., $\omega $ ranges from 300 to 700 r/min) while $\mathrm{\Omega}$ is maintained at 1, thus showing the linearized gas film coefficients as a function of $\mathrm{\Lambda}$ (i.e., $\omega $). In Fig. 5(b), $\mathrm{\Lambda}$ is maintained at 2.5233 (i.e., $\omega =$ 300 r/min) while $\mathrm{\Omega}$ is varied from 0.5 to 50, thus showing the linearized gas film coefficients as a function of $\mathrm{\Omega}$. For each figure, the upper component corresponds to the stiffness, while the lower corresponds to the damping. In contrast with the frequencyindependent stiffness of Ref. [18], the stiffness presented here is frequency dependent. For the stiffness, the direct terms, i.e., ${k}_{zz}$ and ${k}_{xx}$, increase with the increased $\mathrm{\Lambda}$ or $\mathrm{\Omega}$, but the cross terms, i.e., ${k}_{xy}$, decrease with the increased $\mathrm{\Lambda}$ or $\mathrm{\Omega}$. The tendency of the stiffness is opposite to that of the damping. The tendency of linearized gas film coefficients is in line with Ref. [14]. With respect to the effect of the slip flow, it can be drawn that the linearized stiffness and damping coefficients of the gas film are both larger under the slip flow condition because of the smaller central film thickness at equilibrium under the slip flow condition, as shown in Fig. 3. These linearized gas film coefficients serve as the basis of the research with regard to stability and tracking.
3.2. Stability
Fig. 6 displays a stability example where the stator has an axial initial deviation (i.e., ${z}^{*}\left(0\right)=$ 10 μm) and an angular deviation (i.e., ${\gamma}_{x}^{*}$(0)$=$ 10 μrad) when $\omega =$ 300 r/min (i.e., $\mathrm{\Lambda}\mathrm{}=$ 2.5233). Note that in the case, the slip flow is considered and the face contact is ignored. It can be seen that the oscillations of the stator in the axial and angular modes both rapidly relax to stable state because of the nonlinear gas film and the flexible support. As the sealsystem parameters have been properly designed here, the stator does not exhibit underdamping vibrations.
Fig. 6Stability example where the stator has axial and angular initial deviations (i.e., z*(0)= 10 μm, γx*(0)= 10 μrad) when ω= 300 r/min (i.e., Λ= 2.5233)
Fig. 7Comparison of critical transverse moment of inertia at different shaft speeds ω from 300 to 700 r/min (Λ from 2.5233 to 5.8878), and “relaxation region” offers a softer demand of angular stability when considering the slip flow
In the axial mode, it is clear that the axial damping coefficient ${c}_{zz}$ shown in Fig. 5 never becomes negative over the range of $\mathrm{\Lambda}$ from 2.5233 to 5.8878 or over the excessively wide range of $\mathrm{\Omega}$ from 0.5 to 50. Therefore, the stability of the gas face seal in the axial mode is not a significant factor to be focused on. However, the gas film coefficients shown in Fig. 5 can satisfy Eq. (24), and thus the angular stability of the gas face seal should be analyzed. Fig. 7 displays the dimensionless critical transverse moment of inertia ${I}_{cr}$ as $\omega $ increases from 300 to 700 r/min. It can be seen that the critical transverse moment of inertia is greater under the slip flow condition, and that the relative error decreases as $\omega $ increases. Here, the deviation is marked in yellow and is called “relaxation region”. When the ratio of the real transverse moment of inertia to the critical transverse moment of inertia is larger than 1, the face seal is known to exceed its critical speed, indicating that an instability appears in the angular modes. Therefore, for a specific safety threshold (i.e., less than 1), the engineers are provided with an extra relaxation region for the stator design (in terms of the mass, the transverse moment of inertia and the geometry) to satisfy the demand for angular stability when considering the slip flow under lowparameter conditions. This advantageous effect under the slip flow condition is caused by the relative large linearized stiffness and damping coefficients of the gas film under the slip flow condition, which are related to the reduced central film thickness at equilibrium under the slip flow condition. Therefore, these successional phenomena essentially occur because of the loss of the viscous pumping due to the slippage at the solid boundary.
3.3. Tracking
Fig. 8 displays a tracking example where the rotor has a periodic vibration with ${A}_{rz}^{*}=$ 100 μm, ${A}_{r\gamma}^{*}=$ 100 μrad and $\omega =$ 300 r/min (i.e., $\mathrm{\Lambda}\mathrm{}=$ 2.5233). In the case, the slip flow is considered. The stator tracks the rotor in the axial and angular modes. In the close look of Fig. 8(a), the vibration amplitude of the stator is reduced, and a phase difference occurs between the mated rings. These phenomena also exist in the angular modes in Fig. 8(b).
Fig. 8Tracking example where Arz*= 100 μm, Arγ*= 100 μrad, ω= 300 r/min (i.e., Λ = 2.5233) and the slip flow is considered: a) axial motion; b) angular motions
a)
b)
Fig. 9 shows the axial and angular relative transmissibility at different shaft speeds $\omega $ ranging from 300 to 700 r/min (i.e., $\mathrm{\Lambda}$ ranges from 2.5233 to 5.8878). The relative transmissibility increases with an increase in shaft speed, indicating that the tracking property worsens as $\omega $ increases. Moreover, it is obvious that the relative transmissibility is smaller under the slip flow condition. A reasonable explanation is that the slippage phenomenon will induce relative large linearized gas film coefficients which help to transmit the rotor motion to the stator. As the red curves show, this acceleration decreases as the slip flow diminishes. However, this advantageous effect of the slip flow on the tracking property should be doubted because the transmissibility is a ratio that used to only reflect the message about the amplitude rather than about the datum line. Namely, the improved tracking property under the slip flow condition is caused by the relative large linearized gas film coefficients under the slip flow condition, which are related to the decreased central film thickness at equilibrium under the slip flow condition. Although the tracking ability is better under the slip flow condition, the space for tracking motions has been compressed because of the decreased central film thickness at equilibrium. Therefore, an index including both “amplitude” and “datum line” is required. The critical ambientdisturbance amplitude, above which the seal will have face contact, are provided based on the transmissibility to cover the defect of transmissibility in characterizing the tracking property.
In engineering projects, the axial endplay and the angular wobble occur simultaneously. In Fig. 10, single axialendplay and single angularwobble situations are discussed respectively. The critical amplitude is smaller under the slip flow condition. It means that if an engineer ignored the slip flow and then selected a value in the “rigorous region” marked in yellow as a design criterion, face contact would occur. Therefore, the slip flow makes gas face seals require more conservative design and more rigorous manufacture and assembly for a noncontacting tracking. This disadvantageous conclusion drawn from the critical amplitude is in contrast with the advantageous conclusion obtained by the transmissibility, because the former considers the tracking property more comprehensively by adding the space for tracking motions. In addition, with the increase of $\omega $, ${A}_{rzcr}$ and ${A}_{r\gamma cr}$ vary nonmonotonously. It is because the critical amplitude is a function of the central film thickness at equilibrium (Fig. 3) and the relative transmissibility (Fig. 9).
Fig. 9Comparison of relative transmissibility at different shaft speeds ω from 300 to 700 r/min (Λ from 2.5233 to 5.8878)
Fig. 10Comparison of critical amplitudes of axial endplay and angular wobble for face contact at different shaft speeds ω from 300 to 700 r/min (Λ from 2.5233 to 5.8878), and “rigorous region” requires more rigorous working conditions when considering the slip flow
4. Conclusions
The dynamics of a spiral groove gas face seal are established and the BoltzmannReynolds model with FukuiKaneko approximation is used to account for slip flow. Transmissibility in the step jump method is introduced into the perturbation method to characterize the tracking property. Based on the transmissibility, the critical amplitudes of axial endplay and angular wobble, above which the seal will have face contact, are provided to cover the defect of transmissibility in characterizing the tracking property.
The effect of slip flow on the stability and tracking properties of a spiral groove gas face seal under lowparameter conditions is researched. The results show that slip flow will induce a decreased central film thickness at equilibrium, a reduced load carrying capacity and a higher liftoff speed (i.e., a longer startup operation and more severe face wear). The decreased central film thickness induces to relative large linearized gas film characteristics that affect the stability and tracking. Unlike the axial stability, the angular stability should be focused on. Slip flow will induce an extra “relaxation region” in the critical transverse moment of inertia to provide a softer demand for angular stability. With regard to the transmissibility, slip flow will induce an advantageous effect on transmitting the rotor motion to the stator. However, after introducing the space for tracking motions, the results of the critical amplitudes of the axial endplay and the angular wobble show that slip flow will induce an extra “rigorous region” to make gas face seals require more conservative design and more rigorous manufacture and assembly for a noncontacting tracking. The proposition of “relaxation region” and “rigorous region” demonstrates the necessity of considering slipflow under the lowparameter applications of gas face seals.
References

Huang W., Lin Y., Liu Y., Liu X., Gao Z., Wang Y. Face rubimpact monitoring of a dry gas seal using acoustic emission. Tribology Letters, Vol. 52, Issue 2, 2013, p. 253259.

Huang W., Lin Y., Gao Z., Fan W., Suo S. An acoustic emission study on the starting and stopping processes of a dry gas seal for pumps. Tribology Letters, Vol. 49, Issue 2, 2013, p. 379384.

Hu S., Huang W., Liu X., Wang Y. A fluctuation transient response analysis of spiral groove gas face seals. STLE Tribology Frontiers Conference, Chicago, USA, 2014, p. 41.

Hu S., Huang W., Liu X., Wang Y. Numerical model for dynamic behavior of spiral groove gas face seals using fractal contact model. 70th STLE Annual Meeting, Dallas, USA, 2015, p. 150.

Sayed M., Hamed Y. S. Stability analysis and response of nonlinear rotorseal system. Journal of Vibroengineering, Vol. 16, Issue 8, 2014, p. 41524170.

Soltani M., Keshmiri M., Misra A. K. Dynamic analysis and trajectory tracking of a tethered space robot. ACTA Astronautica, Vol. 128, 2016, p. 335342.

Shapiro W., Colsher R. Steadystate and dynamic analysis of a jet engine, gas lubricated shaft seal. ASLE Transactions, Vol. 17, Issue 3, 1974, p. 190200.

Green I., Barnsby R. M. A simultaneous numerical solution for the lubrication and dynamic stability of noncontacting gas face seals. Journal of Tribology, Vol. 123, Issue 2, 2000, p. 388394.

Miller B. A., Green I. Numerical formulation for the dynamic analysis of spiralgrooved gas face seal. Journal of Tribology, Vol. 123, Issue 2, 2000, p. 395403.

Green I., Barnsby R. M. A parametric analysis of the transient forced response of noncontacting conedface gas seals. Journal of Tribology, Vol. 124, Issue 1, 2002, p. 151157.

Green I. A transient dynamic analysis of mechanical seals including asperity contact and face deformation. Tribology Transactions, Vol. 45, Issue 3, 2002, p. 284293.

Zirkelback N., Andres L. S. Effect of frequency excitation on force coefficients of spiral groove gas seals. Journal of Tribology, Vol. 121, Issue 4, 1999, p. 853863.

Zirkelback N. Parametric study of spiral groove gas face seals. Tribology Transactions, Vol. 43, Issue 2, 2000, p. 337343.

Ruan B. A semianalytical solution to the dynamic tracking of noncontacting gas face seals. Journal of Tribology, Vol. 124, Issue 1, 2002, p. 196202.

Elrod H. G., McCave J. T., Chu T. Y. Determination of gasbearing stability by response to a stepjump. Journal of Lubrication Technology, Vol. 89, Issue 4, 1967, p. 493498.

Miller B. A., Green I. Numerical techniques for computing rotordynamic properties of mechanical gas face seals. Journal of Tribology, Vol. 124, Issue 4, 2002, p. 755761.

Miller B. A., Green I. Semianalytical dynamic analysis of spiralgrooved mechanical gas face seals. Journal of Tribology, Vol. 125, Issue 2, 2003, p. 403413.

Ruan B. Finite element analysis of the spiral groove gas face seal at the slow speed and the low pressure conditionsslip flow consideration. Tribology Transactions, Vol. 43, Issue 3, 2000, p. 411418.

Burgdorfer A. The Influence of the Molecular Mean Free Path on the Performance of Hydrodynamic Gas Lubricated Bearings. AECU3771, IA20492, Franklin Institute, Philadelphia, 1958.

Hsia Y. T., Domoto G. A. An experimental investigation of molecular rarefaction effects in gas lubricated bearings at ultralow clearances. Journal of Lubrication Technology, Vol. 105, Issue 1, 1983, p. 120130.

Mitsuya Y. Modified Reynolds equation for ultrathin film gas lubrication using 1.5order slipflow model and considering surface accommodation coefficient. Journal of Tribology, Vol. 115, Issue 2, 1993, p. 289294.

Fukui S., Kaneko R. Analysis of ultrathin gas film lubrication based on linearized Boltzmann equation: first reportderivation of a generalized lubrication equation including thermal creep flow. Journal of Tribology, Vol. 110, Issue 2, 1988, p. 253261.

Fukui S., Kaneko R. A database for interpolation of Poiseuille flow rates for high Knudsen number lubrication problems. Journal of Tribology, Vol. 112, Issue 1, 1990, p. 7883.

Pecht G. G., Netzel J. P. Design and application of noncontacting gas lubricated seals for slow speed services. Lubrication Engineering, Vol. 55, Issue 7, 1999, p. 2025.

Ruan B. Numerical modeling of dynamic sealing behaviors of spiral groove gas face seals. Journal of Tribology, Vol. 124, Issue 1, 2002, p. 186195.

Green I., Etsion I. Stability threshold and steadystate response of noncontacting conedface seals. ASLE Transaction, Vol. 28, Issue 4, 1985, p. 449460.

Malanoski S. B., Pan C. H. T. The static and dynamic characteristics of the spiralgrooved thrust bearing. Journal of Basic Engineering, Vol. 87, Issue 3, 1965, p. 547558.

Wilcock D. F, Bjerklie J., Cheng H. Design of floated shoe close clearance seals for supersonic jet engine compressors. Journal of Lubrication Technology, Vol. 90, Issue 2, 1968, p. 500509.
Cited by
About this article
This work was supported by the National Key Basic Research Program of China (973) (Grant No. 2012CB026003), the National Science and Technology Major Project (Grant No. ZX06901), and the National Science and Technology Support Plan Projects (Grant No. 2015BAA08B02).