A method using reduced-order Kalman filters is developed to estimate the thin gas film axial force and moments in real-time for mechanical gas face seal systems in a flexibly mounted stator configuration. First-order Gauss-Markov stochastic models, combined with the stator motion equations, form the basis for the reduced-order Kalman filter estimators. Two schemes are presented to estimate axial force and moments based on stator motion measurements. In one scheme, the force and moments are directly estimated and, in another scheme, a set of Proper Orthogonal Decomposition (POD) weighting functions are estimated, from which the gas film force and moments are computed. Both estimators are shown to approximate the gas film axial force and moments successfully for different forcing functions over a wide range of compressibility numbers.
I. INTRODUCTIONNoncontacting mechanical gas face seals are often used in high speed, high performance rotating machinery, such as pumps, compressors and gas turbine engines. A schematic of a mechanical gas face seal is shown in Figure 1. Sealing is provided by forcing the gas to flow through a small gap between the stator and rotor mating faces. The seal dynamic behavior is characterized by the forces generated within the thin gas film and the flexible support. The thin gas film has nonlinear dynamic characteristics that affect the seal system behavior through the hydrostatic and hydrodynamic force and moments acting on the stator. Different methods have been developed to study gas film force effects. For example, the direct numerical method employs fully meshed discretizations of the gas film governing equations using finite element (FEM), finite volume (FVM) or finite difference (FDM) methods to couple the lubrication analysis with the stator equations of motion [1]. However, the direct numerical method is not conducive to parametric studies of seal geometry or changing operating conditions. Some analytical and semianalytical methods generate a linearized description of the gas film stiffness and damping properties and can be used to obtain the gas film effects. Examples include the small perturbation method [2-4], the step jump method [4], and the direct numerical frequency response method [5]. However, these analytical and semi-analytical methods are only applicable for small motions about a specific operating point, and the obtained force coefficients are valid only within that range [6]. When the system