This article proposes a rotordynamic identification method using the accelerating sweep frequency excitation method. The computational fluid dynamics transient solution combined with the moving grid method is utilized to obtain the transient flow field of the seal excited by the whirling rotor with an accelerating frequency. Rotordynamic coefficients at swept frequencies are obtained by analyzing the transient response force acting on the rotor. Rotordynamic coefficients of three published experimental seals including a labyrinth seal, a fully partitioned pocket damper seal and a honeycomb seal are identified to validate the proposed method. The results show that the predicted rotordynamic coefficients are all well agreement with the trend of experimental data. Compared with the existing numerical models based on transient solutions, the central processing unit consumption of the proposed method is substantially reduced by 98% when achieving the same frequency resolution. In addition, the impact of the exciting acceleration on the identification accuracy is also illustrated.

The identification method using infinitesimal theory is proposed to predict rotordynamic coefficients of annular gas seals. The transient solution combined with moving grid method was unitized to obtain the fluid reaction force at a specific position under different whirling frequencies. The infinitesimal method is then applied to obtain the rotordynamic coefficients, which agrees well with published experimental results for both labyrinth seals and eccentric smooth annular seals. Particularly, the stability parameter of the effective damping coefficient can be solved precisely. Results show that the whirling frequency has little influence on direct damping coefficient, effective damping coefficient, and cross-coupled stiffness coefficient for the labyrinth seal. And the effective damping coefficients decrease as the eccentricity ratio increases. A higher eccentricity ratio tends to destabilize the seal system, especially at a low whirling frequency. Results also show that the fluid velocity in the maximum clearance in the seal leakage path is less than that in the minimum clearance. The inertial effect dominates the flow field. Then it results in higher pressure appearing in maximum clearances. The pressure difference aggravates the eccentricity of rotor and results in static instabilities of the seal system.

Nonuniform mesh is beneficial to reduce computational cost and improve the resolution of the interest area. In the paper, a cell-based adaptive mesh refinement (AMR) method was developed for bearing cavitation simulation. The bearing mesh can be optimized by local refinement and coarsening, allowing for a reasonable solution with special purpose. The AMR algorithm was constructed based on a quadtree data structure with a Z-order filling curve managing cells. The hybrids of interpolation schemes on hanging nodes were applied. A cell matching method was used to handle periodic boundary conditions. The difference schemes at the nonuniform mesh for the universal Reynolds equation were derived. Ausas’ cavitation algorithm was integrated into the AMR algorithm. The Richardson extrapolation method was employed as an a posteriori error estimation to guide the areas where they need to be refined. The cases of a journal bearing and a thrust bearing were studied. The results showed that the AMR method provided nearly the same accuracy results compared with the uniform mesh, while the number of mesh was reduced to 50–60% of the number of the uniform mesh. The computational efficiency was effectively improved. The AMR method is suggested to be a potential tool for bearing cavitation simulation.