Numerical Investigation of the FSI Characteristics in a Tubular Pump ShuoWang ,

Flow condition was simulated in a shaft tubular pump by using the Shear-Stress Transport (SST) k-ω turbulence model with high quality structured grids in design condition. Corresponding structural vibration characteristics were then analyzed based on twoway coupled Fluid-Structure Interaction (FSI) method. Fluid results showed that flow in the outlet flow passage was a combination of the axial flow and circumferential rotationmotion. Time and frequency domain analysis of pressure pulsation of typical measure points indicated that larger pulsation amplitudes appeared in the tip of the blades and the main vibration source was the pressure pulsation induced by rotation of the blades.The fluid pulsation amplitudes decreased gradually along the flow direction, which can be ascribed to the function of fixed guide vane. Structural analysis of the blades in both pressure and suction side indicated that significant stress concentration was formed at the blade and hub connection near the leading edge. Maximum effective stress of the blades varied periodically, so prevention measures of the fatigue of blades should be taken. This research can provide important reference for the design of the tubular pump.


Introduction
Rapid development of the industry and agriculture puts great demand on the water conservancy facilities, which gives rise to more pumping stations.Because of its less hydraulic loss, high efficiency, and compact structure, convenient installation, and maintenance, tubular pump, a kind of axial flow pump, presents wide popularity and development potential worldwide.Nevertheless, it is a bit difficult for the design and performance prediction of the tubular pump.Experimental model testing is one of the solutions for performance prediction but it is costly and time-consuming [1].By contrast, numerical simulation is a powerful tool to provide information of the fluid flow behavior accurately, thus helping the scholars to obtain a thorough performance evaluation of a specific design [2][3][4].
Numerous CFD simulations were performed to analyze the flow features of the pump and then corresponding optimization design was provided.Jafarzadeh et al. [3] presented a general three-dimensional simulation of a highspeed centrifugal pump to predict velocity and pressure field.The optimum pumping operation was explored in terms of turbulence models and the number of blades.Shojaeefard et al. [5] investigated the effects of the blade outlet angle and passage width on the centrifugal pump performance during the pumping of water and oil, and they found pump head and efficiency can be increased with a proper modification of the original geometry.Zhu et al. [6] applied a multiobjective optimization design system to the design of a middle-highhead pump-turbine runner and concluded that the choice of blade loading and the meridional channel shape is crucial to efficiency and operation stability.All their researches have beneficial enlightenment, but without taking the coupled effect of the fluid and solid into account, their concern is mainly the fluid feature.
Analysis of the stress distribution characteristics of a centrifugal pump was conducted in [7,8].Results showed that the maximum stress on the blade appeared on the pressure side near the hub, and the maximum static stress increased with the decreasing of the flow rate.However, like previously mentioned researches, they focused on the centrifugal pump, and enough importance was not attached to the axial flow pump.
Shi and Wang [9], Tang et al. [10], and Zhang et al. [11] calculated the stress and deformation of the axial impeller and found similar law of stress and deformation distribution.But guide vanes and tip clearance were omitted for simplification in their structure simulation, which could not reflect the actual characteristics of the pump.There is tremendous need for the detailed study of the axial flow pump based on twoway coupled Fluid-Structure Interaction approach [12].
On account of the scarcities above, this research performed simulation of the actual flow in a shaft tubular pump to acquire the general features of the fluid flow and vibration characteristics of the blades in design condition.Based on two-way coupled FSI approach, a three-dimensional unified model for a shaft tubular pump was established, including the whole flow passage, the structure of the impeller, and the fixed guide vane.Analysis of flow pattern was performed using the SST - turbulence model to obtain the general law in the flow passage.Time and frequency domain analysis of the fluid pressure pulsation was conducted to explore its vibration characteristics.For the structure, features of effective stress, strain, deformation, and velocity distribution in both pressure and suction sides of the blades were analyzed and the property of maximum dynamic stress was evaluated.Note that the whole calculation was performed with the commercial software ADINA.

Basic Theory and Method
2.1.Basic Equations.For the fluid flow analysis, the continuity equation and Reynolds-Averaged Navier-Stokes equation for an incompressible flow have been used in the following form [13]: where   is an average velocity component, and       is the turbulent stress.In this research Shear-Stress Transport (SST) - turbulence model proposed by Menter is used to acquire the features of unsteady flow in CFD calculation.SST - turbulence model can yield more accurate and reliable results for simulating flows with adverse pressure gradient [14], flows around complex geometry [15], transonic flows [16], and so forth.Equation for SST - turbulence model is expressed as follows [17][18][19]: where  and  are turbulent kinetic and turbulence dissipation rate, respectively,   is a productive term of the turbulent kinetic,   is a productive term of the turbulence dissipation rate, Γ  and Γ  are the effective diffusion coefficients of  and , respectively,   and   are the dissipation terms of  and , respectively,   is the cross-diffusion term, and   and   are source terms.The fundamental conditions applied to the FSI [20] are the kinematic condition (or displacement compatibility) and the dynamic condition (or traction equilibrium), and corresponding equations are expressed in (3) and (4), respectively.
where   and   are, respectively, the fluid and solid displacements, and   and   are, respectively, the fluid and solid stresses.The underlining denotes that the values are defined on the fluid-structure interfaces only.
When no-slip or slip condition is applied separately, the fluid velocity condition that resulted from the kinematic condition is expressed in ( 5) and ( 6), respectively.
Fourier analysis is employed in vibration analysis of the pressure pulsation in fluid domain and dynamic stress in solid domain.Given a time history (), defined from  0 to  1 , the Fourier series corresponding to () is expressed as follows: where 7) could be rewritten as (8) due to sine and cosine transforms: where Conceptually, the function () is considered to contain the frequencies 0,  1 , 2 1 , . ... The constants   and Φ  are the amplitude and phase angle of that portion of () which oscillates at frequency  1 .compatibility, continuity, and completeness conditions along the nonconforming interfaces.Wall boundary condition is imposed on the outmost layer of the flow body with noslip velocity assumed.For the interface between fluid and structure, FSI boundary condition is assigned to make the mutual transmission of velocity and pressure.

Numerical
The commercial software ADINA is adopted to study the performance of the pump under transient state in ADINA CFD module and dynamic-implicit analysis in ADINA structure module.The entire fluid domain of the pump is formed by combining the components with an interface between impeller and guide vanes.Total calculation time is set as 10 seconds to obtain reliable simulation results.There are 2000 steps calculated and the time step is 0.005 seconds.the water flows into the impeller.Although most circulation generated by the work of impeller is rectified and recovered by the fixed guide vane, there is still certain velocity circulation and reflux in the outlet flow passage, which can also be ascribed to the diffuse type of the outlet flow passage.The streamlines deflect in the outermost section and have friction with the pump casing.In total, the design of the pump flow channel is reasonable, so the stability and safety of operation can be guaranteed.

Analysis of Fluid Pressure
Pulsation.Vibration analysis of pressure pulsation enhances both hydraulic and dynamic performance and assures safety and reliability of the pumping station.As shown in Figure 5, four measure points are chosen along the radius direction to conduct vibration response analysis in time and frequency domain, the locations of which are from blade tip to hub, respectively.The measure sections where the measure points are placed are before the impeller, between impeller and the guide vanes, and behind the waterguide cone.Time history curves and spectrum analysis of pressure pulsation are presented in Figures 6-8 and first main frequency and amplitude are shown in Table 1.Since fluctuation of the nodal pressure remains periodically steady at about  = 6s, the regularity of pressure pulsation is analyzed within the time range of  = 6s to  = 10s.It is clear that the pressure pulsation of the measure points at the same cross section has the similar trend.For the section before impeller, the pressure magnitude varies between 30000 Pa and 44000 Pa and increases with less amplitude from tip to hub with a dominant frequency of 5.25 Hz, shown in Figure 6  measure points between impeller and guide vanes varies between 40000 Pa and 50000 Pa and has similar vibration features, which illustrates that the pulsation is mainly caused by blades rotation.For the section after water-guide cone, the pressure magnitude varies between 47000 Pa and 51000 Pa and decreases with more amplitude from tip to hub with a dominant low frequency of 1.50 Hz, illustrated in Figure 8(b).
When analyzing these measure points together, we find that the pulsation amplitudes decrease gradually along the flow direction; that is, the pressure pulsation amplitude before the impeller is larger than that between impeller and the guide vanes, and both amplitudes of them are larger than that behind the water-guide cone.Due to the     constraint function of the guide vanes, the pressure tends to be stable and the pulsation amplitude is greatly reduced.Great attention should be paid to the inlet of the impeller chamber in future hydraulic design as the disturbance of the blades has already existed before the water flowed into the impeller.

Effective Stress and Strain Analysis of the Blades.
For pump design, the major mission is to analyze the structural response to the flow induced excitation and diagnose the faults early.Figure 9 shows the effective stress distribution of the blades in both pressure and suction side at  = 10 s in dynamic design operation.The maximum effective stress of blades is 23.59 MPa, appearing at the blade root near the leading edge.Because of the huge flow impact generated in the leading edge of the blades, the stress concentration is formed and it is prone to crack and fracture.The effective stress distribution basically declines from the joint of the blade root to the trailing edge and the tip, where the stress is much lower.Figure 10 is the corresponding strain at  = 10 s.The distribution of the strain is coincident with that of the stress.11 is the deformation distribution of the blades at  = 10 s in dynamic design operation.The deformation magnitude augments with the increase of the radius in both pressure and suction side.On account of the larger centrifugal force, thinner thickness, and insufficient rigidity of the tip, the maximum deformation appears at the blade tip, which is consistent with the result of the maximum effective stress distribution.

Velocity Analysis of the Blades.
Velocity distribution of the blades at  = 10 s in dynamic design operation is shown in Figure 12.We can see that velocity increases gradually with the increase of radius.The maximum velocity magnitude appears at the blade tip, which is coincident with the manual calculation result; that is, V =  ⋅  = 10.99 × 1.625 = 17.86 m/s.

Dynamic Stress
Analysis of the Blades.Figure 13 is the time history curve and Fourier frequency analysis of maximum effective stress of the blades and Table 2 is corresponding magnitude of the main frequency and amplitude.It is clear that the maximum effective stress changes periodically as time changes, with magnitude varying between 22 and 25 MPa.Frequency analysis illustrated in Figure 13(b) and Table 2 shows that the first main frequency of the maximum effective stress is 1.75 Hz, and the corresponding amplitude is 0.43 MPa, which illustrates the periodically varying stress is mainly influenced by the rotational frequency of impeller.The static stress can not make the blades damaged because the value of maximum effective stress is much less than the limit stress of material.However, practice shows that the alternating stress could cause abrupt fracture of the structure with no significant plastic deformation before occurrence of fracture, even if the stress magnitude is much lower than the yield limit.Therefore, prevention measures of the fatigue of impeller should be taken due to the period variation of the effective stress.

Conclusion
In present investigation, numerical simulation of flow condition in a shaft tubular pump was performed and the vibration characteristics of the blades in design condition were analyzed.The numerical method was based on the SST - turbulence model and the two-way coupled FSI approach under transient state.The following conclusions can be made by analysis of the results: (1) The flow pattern is steady in the inlet flow passage in design condition.Due to the incomplete recovery of the flow circulation by the fixed guide vane and the diffusion shape of the outlet flow passage, vortex occurs in the center and the streamlines deflect in the outer section of the outlet flow passage, which illustrates that the flow in outlet flow passage is a combination of the axial flow and circumferential rotation motion.In total, the design of the pump flow channel is rational and the stability and safety of operation can be guaranteed.
(2) Time and frequency domain analysis shows that the fluid pressure magnitude increases with less amplitude from blade tip to hub before impeller and between impeller and guide vanes with a dominant frequency of 5.25 Hz.The fluid pressure magnitude decreased with more amplitude from tip to hub after water-guide cone, and its dominant frequency is 1.50 Hz.Larger pulsation amplitudes appear in the tip of the blades, which indicates that optimizing the airfoil shape and structural parameters of the blade tip is an important way to improve the pressure pulsation condition in pump operation.The dominant frequency shows the main vibration source in fluid flow is the pressure pulsation induced by blade rotation.The pulsation amplitude decreases along the flow direction  with the maximum before the impeller, which illustrates that the fixed guide vane has a function of restricting pressure pulsation and greatly improving the flow pattern in the pump.
(3) Analysis of the stress, strain, deformation, and velocity distribution of the blades is performed based on FSI.The result demonstrates that the maximum effective stress occurs between blade and hub near the leading edge, which requires the designers to attach importance to the strength check of the root.The stress at trailing edge and the tip is comparatively smaller, which is coincident with the result of strain.The tip has maximum deformation, which is prone to friction damage of the blade tip and pump casing because of the tiny clearance between them.Maximum velocity occurs at the tip, which is coincident with the manual calculation results.Although the value of maximum effective stress is much less than the yield limit, prevention measures of the fatigue of impeller should be taken due to the period variation of the effective stress.

Figure 3 :
Figure 3: Three-dimension distribution of velocity magnitude and streamlines.

Figure 4 :
Figure 4: Distribution of velocity magnitude and streamlines in typical sections.
(b), which equals the product of blades number 3 and rotation frequency 1.75 Hz.Pressure magnitude of the

Figure 5 :
Figure 5: Measure point distribution in the flow passage of the pump.

Figure 9 :
Figure 9: Effective stress distribution of the blades.

Figure 12 :
Figure 12: Velocity distribution of the blades.

Figure 13 :
Figure 13: Time and frequency domain graph of maximum effective stress of the blades.

Table 1 :
First main frequency and amplitude of measure points.

Table 2 :
Main frequency and amplitude of maximum effective stress of the blades.