Rotordynamic Analysis for a Turbo-Machine with Fluid-Solid Interaction and Rotation Effects

The rotordynamics and fluid dynamics of a turbo-machine considering the effect of fluid-solid interaction FSI are numerically investigated using finite element software ADINA. The iterative method is adopted in computation of coupled fields of displacement and fluid. What distinguishes the present study from previous ones is the use of ADINA’s rotational meshes and the FSI interface that separates the rotor surface from its surrounding fluid. The rotor’s center orbit and frequency response as well as the transient fluid dynamics are obtained with various axial flow speeds. By including real rotating motion of the rotor, this paper presents a better way to solve complicated rotordynamic problems of turbo-machines that are operated in FSI circumstances.


Introduction
Fluid solid interaction FSI analyses play very important roles in understanding rotordynamics of turbo-machines. It has been noticed that fluid and solid as parts of the rotor system influence each other and severe vibrations may be excited in some situations. Particularly, this is true for machines transferring heavy fluid, for example, pumps or other turbomachines. Historically, the FSI-induced vibration was studied empirically or experimentally. For instance, Childs 1, 2 carried out a series of experiments to investigate the FSI forces between the impeller and the shroud of a pump as well as the resulted axial vibrations. Marscher 3 determined the "wet" critical speeds of a multistage pump rotor. Numerical methods such as the finite element method FEM and coupled FEM with the boundary element method BEM have been widely applied to various FSI problems.  who is the developer of the ADINA software proposed application of FEM to solve the FSI problems for compressible and incompressible fluid flow by iterative method and direct method. Bermúdez et al. 7 demonstrated the plausibility and effectiveness of solving the FSI problems by using FEM as well. Kwon and Jo 8 modeled two 3D beams and external flow around them and used coupled lattice Boltzmann method LBM and FEM to analyze the characteristics of the flow field. Czygan and von Estorff 9, 10 and Soares et al. 11, 12 adopted both the FEM and BEM to solve the FSI problems. They used an iterative method which coupled of finite element and boundary element for the time domain to analyze fluid-solid systems. Michelin and Llewellyn Smith 13 introduced an unsteady point vortex method for two-dimensional coupled motion of a general sharp-edged solid body and a surrounding inviscid flow. The application of these methods makes integrated FSI analysis effective and computationally plausible for complex rotordynamic problems.
In this paper, the FSI problem is investigated for a turbo-machine to obtain its rotordynamic response by using finite element software ADINA with consideration of real rotating movement of the rotor surrounded by incompressible fluid flow. The fields of solid and fluid are modeled with ADINA's rotational meshes and solved iteratively for pressure and velocity of the fluid field as well as displacement of solid rotor. The first "wet" critical speed of the rotor is identified through a signal filtrating with different axial flow velocities. The natural frequency is changed from 15.63 Hz to 14.65 Hz at different Reynolds numbers.

Governing Equations
Consider a problem domain Ω consisting of a fluid part Ω f and a solid part Ω s . For the fluid domain Ω f we assume a flow of incompressible Newtonian fluid. The basic equations governing continuity and momentum are given by

2.1
where Einstein's summation convention is adopted here and throughout the rest of the paper. V υ 1 , υ 2 , υ 3 , υ i is the velocity vector with respect to Cartesian coordinates x i , ρ f the fluid density, D/Dt ∂/∂t V ·∇ is material derivative, and f fi are the body forces per unit mass. The stress tensor T ij for incompressible Newtonian fluids can be expressed as where p is pressure and μ f is dynamic viscosity. An empirical formula is used to compute Reynolds number, as follows Re Engineering   3 where D e is characteristic diameter denoted as

Mathematical Problems in
where A is flow area and χ is the length of fluid-solid interface attached to the flow area. In this model, the Reynolds number is approximately between 1.108 × 10 6 and 2.216 × 10 7 , which is larger than its empirical critical value 1100 for turbulent flows. Hence, the numerical simulation of the flow field is carried out by using one with the turbulence model. The viscosity coefficient for the turbulent flow is modified in ADINA, as where μ 0 is the laminar viscosity and μ t is its turbulence counterparts.
In the present paper, the K-ω model for high-Reynolds-number turbulences in incompressible flows is adopted. Based on the normal K-ε model, a parameter ω is introduced which is related to K and ε as where ε and K are dissipation rate and kinetic energy of turbulence, respectively. Then, μ t is computed by The governing equations for K and ω are referred to ADINA 14 .
The momentum equation for the solid domain Ω s is expressed as where u i is the displacement vector, σ ij denotes the Cauchy stress tensor, ρ s is the density of the solid material, and f si are body forces per unit mass acting on the solid. By assuming deformation is small, the strain tensor ε ij is related to the displacements u i by the kinematic relation On the fluid-solid interface the nonslip boundary conditions is applied for velocities and stresses on the interface where subscripts s and f denote the solid and fluid domains, respectively.

Computational Model and Method
The sketch of the analysis model and its finite element meshes are shown in Figure 1 unit: mm . The radiuses of shaft and impeller are 50 mm and 70 mm, respectively. The external diameter of the blade is 300 mm. The total lengths along the z-axis are 2002 mm for the rotor and 150 mm for the flow field, respectively. The external diameter of the flow field is 310 mm. The blades are straight and have a 3.814 • angle with respect to z-axis. 736 8-node-hexahedral elements and 1254 10-node-tetrahedral elements are used to mesh the solid field. For the fluid field, there are 4032 8-node-hexahedral elements. In ADINA a special kind of element called rotational meshes are used to simulate the rotation of the rotor surrounded with fluid. The material parameters are displayed in Table 1.
The rotation speed of the rotor is 600 rpm. The unbalanced mass eccentricity of the rotor is 2.33 μm which is applied to the model by loading it with a following force F see Figure 1 at the center point O. The magnitude of F is a constant 1 N, and its initial direction is chosen to be along the x axis in the Cartesian coordinate. To simulate the rotor in real rotating status, the direction of F is made to revolve synchronously with the rotor so that its loading relative to the rotation does not change. The displacement boundary condition is fixed end: U x U y U z θ x θ y 0. The initial velocity of v z is negative in the z-direction for all fluid domains. This velocity is kept unchanged on the cross-section at the inlet. The time step of FSI is 0.001 s and total computation lasts 2 s. In this case the sampling frequency for fast Fourier transform FFT analysis is 1000 Hz; the Fourier length for the FFT is 1024.
The displacement U x and U y of the rotor center O can be obtained from the computation results and are used in signal analysis to identify frequency components. The empirical mode decomposition EMD method 15 is used to purify displacement signals by numerically filtrating artificial noise. Then the filtrated signals are transformed by FFT and are analyzed to determine the orbital position point O. By cascading response of the solid part obtained from results with different flow velocities, the rotor's first critical speed upon different Reynolds numbers can be determined.

Dry Critical Speed
The rotor's first critical speed is analyzed by using the rotational meshes without considering of fluid. The results give the "dry" critical speed of the rotor that can be compared to the "wet" critical speed obtained with coupling fluid-solid interactions.
The displacements and the orbit of the axial center point O are illustrated in Figure 2. The displacement results of the rotor center are filtered using EMD to remove noise. Figures  2 a , 2 b and, 2 c are x-and y-directional displacements and the orbit of point O before the EMD filtering, respectively. In comparison, Figures 2 d , 2 e , and 2 f are those responses after the filtering. It is clearly seen that indeed the EMD effectively removes noise in the displacement signals while retain their dominating part in frequency domain.
Following the above filtering, the FFT analysis is carried out to identify the frequency components of the rotor response in Figure 3, where the "dry" natural frequencies are marked with their values over each peak.  Based on Figure 6 there is a short period of transient response when the rotor starts rotating along with a simultaneous axial flow see Figures 6 a , 6 b , and 6 c . Once the flow becomes stable, the transient response of the rotor is gone and the motion becomes periodic. Figures 6 d and 6 e present the FFT results for x-and y-displacements. There are clear peaks on FFTs in these two figures, suggesting that the first "wet" critical speed is 14.65 Hz when v z 5 m/s. Next, twelve different axial velocities are used to demonstrate the "wet" critical speed. The results are listed in Table 2 based on FFT analysis for x-and y-displacements.

Results of FSI Simulation
As it can be seen from

Conclusion
In this paper, the rotordynamic response of a turbo-machine is analyzed considering the FSI problem and real rotation movement of the rotor. ADINA's rotational elements are used iteratively to solve the responses of both the solid and fluid fields. The first "wet" critical As the increase of the fluid axial velocity, the "wet" first critical speed changes from 15.63 Hz to 14.65 Hz. Then, the peak gets lower and cannot be distinguished at higher Reynolds numbers. The present study shows that the first critical speed changes with the velocities of the flow. It means the structure has different "wet" first critical speeds in different Reynolds numbers with various flow conditions.