Analysis of Radial Force and Vibration Energy in a Centrifugal Pump

Vibration is one of the main issues taken into consideration in the design and manufacture of the pump. (e radial force and vibration of the impeller induced by fluid in a centrifugal pump were investigated at different flow rates by numerical simulation. (e vibrations on the volute were measured by the experiment. (e variation trend of the radial displacements of the impeller is consistent with that of the radial forces, and the variation in the radial displacement lags that of the radial force. (e vibration energies on the impeller and the volute were analyzed based on root-mean-square (RMS) values in the frequency domain. (e distributions of energy loss in the pumps were calculated to determine the total entropy generation (TEG) and entropy generation rate (EGR). (e TEG values as calculated are in close accordance with hydraulic loss. (e vibration is a result of the poor flow patterns and consequently results in increased energy losses in the pump.


Introduction
Centrifugal pump is a liquid conveying machinery used in petrochemical, aerospace, and chemical industries. Vibration is one of the main issues taken into consideration in the design and manufacture of the pump, especially high-power high-speed centrifugal pumps [1,2]. Finite element analysis of fluid-structure interaction (FSI) has become an important method employed to investigate the vibration and deformation of fluid machinery. Many studies have been conducted on FSI in fluid machinery; most of the studies focused on turbines and compressors, while only few studies have investigated pumps [3][4][5][6]. Benra [7] carried out numerical and experimental studies on flow-induced vibrations of a single-blade pump impeller and reported good qualitative agreement between the FSI calculation and measurement. Katsutoshi et al. [8] evaluated the characteristics of a mixed-flow pump using FSI simulation and the predicted stresses were in good agreement with experimental values within −11 ± 6% prediction accuracy. Pei et al. [9] calculated the deformation and stress in the impeller of an axial-flow pump with two-way passage using FSI simulation and found that the deformation and stress significantly decreased as the flow rate increased. Schneideret et al. [10] numerically evaluated the deformation and stress in impellers of a multistage pump at different operating conditions using FSI simulations. e hydraulic exciting forces generated by unsteady flow are the main factors responsible for vibration and noise in the pump. Pei et al. [11] investigated the hydrodynamic force and flowinduced vibrations of a single-blade sewage water pump impeller at design and off-design flow rates using numerical and experimental methods. Although FSI problems have received significant research interest, only few researches have been carried out to experimentally validate complex FSI problems. In some cases, the deformations can be neglected or are difficult to be measured for moving or rotating structures [12]. Also, the energy loss caused by unsteady flow has a certain correlation with the vibration. Ren et al. [13] used entropy generation rate (EGR) to compare the energy loss in the multistage pump with and without the guide ring. Jia et al. [14] studied the internal instability characteristics of centrifugal pump based on entropy production and vibration energy and found that the internal flow loss has the same variation rule as the vibration energy.
In this study, the radial force and vibration of the impeller induced by fluid in a centrifugal pump were investigated using ANSYS 18.0 commercial software. e radial force and vibration displacement of the impeller were determined at different flow rates. e performance and vibration experiments were conducted on the centrifugal pump. e relationship between energy loss and vibration energy was analyzed based on the entropy generation method. ). e main geometrical parameters are as follows: the inlet diameter of the pump is 65 mm, the outlet diameter of the pump is 50 mm, the impeller diameter is 130 mm, and the number of blades is 5.

Geometrical Model and Numerical
A three-dimensional model of the entire flow field was established using Pro/E, including the inlet, impeller, volute, impeller side chamber, outlet, and annular seal at the impeller inlet. e computational domains for fluid and solid are shown in Figure 1 A grid number test of the computational flow domain was conducted to validate grid independence in the numerical simulation, as shown in Figure 3. e elements of the pump were 1891586, 3301394, 4402428, and 6381025. To ensure a comprehensive consideration of calculation accuracy and efficiency, the total number of elements used in the numerical simulation was 4402428. e elements for the inlet, impeller, volute, outlet, and impeller side room were 490930, 1391721, 1162718, 423291, and 933768, respectively. e grids near the solid wall were refined to satisfy the applicability of the turbulence model. e layers number in the clearance is ten for annular seal of impeller inlet. And it is fifteen for the chamber room of the impeller. e mean y + value of the impeller was 9.99915. For the solid region, there were 88,940 elements and 15,736 nodes, which were meshed with nonstructural grids. e temporal discretization accuracy depends on the time step size, which is of great significance to capture the pressure fluctuation. e influence of the time step was also investigated in this study. e impeller rotates at time steps of 1°, 1.5°, 2°, 3°, and 5°. Figure 4 shows the results of fast Fourier transformation (FFT) of the pressure fluctuation near the volute tongue for the five time steps. It can be observed that the main peak occurred at 245.6 Hz, corresponding to the blade passing frequency. When a smaller time step size was used, the peaks occurred at the same frequency with a minor difference in magnitude. e results indicate that impeller rotation at every 1.5°is sufficient.
In this study, the vibration of the impeller was investigated based on one-way coupling FSI simulations using CFX and Workbench in ANSYS 18.0. e periodically stable results in unsteady flow simulation were used as initial values for FSI simulation. For the fluid region, the SST (shear stress transport) k-ω turbulence model was used in numerical simulation for unsteady calculation and the transient rotor-stator model was applied to couple the interface between the rotating and static domains. e velocity inlet was set as the inlet boundary condition, while the opening was set as the outlet boundary condition. Because the rotating speed of the impeller was 2950 r/min, each time step

External Performance of Experiment and Simulation.
e external characteristics of the head and efficiency of the centrifugal pump obtained from the experiment are shown in Figure 5, and these were compared with those obtained from the unsteady numerical simulation. It can be observed in Figure 5 that the simulated head and efficiency of the centrifugal pump were slightly higher than the experimental values, but both have the same trend. e experimental and simulation values of the head and efficiency at the design flow rate were H e � 17.77 m and η e � 58.87% and H s � 18.40 m and η s � 63.17%, respectively. During simulation, hydraulic smooth surfaces were applied for the pump components. e measured efficiency includes mechanical losses, which were not considered in the simulation. is accounts for some of the observed differences between experimental and simulation results.

Static Pressure Distribution in Centrifugal Impeller.
e static pressure distributions in the middle section of the impeller are shown in Figure 6 at 0.6Q d , 0.8Q d , Q d , and 1.2Q d . From Figure 6, it can be seen that the static pressure in the impeller increases from the inlet to outlet. e pressure is larger on the pressure side than that on the suction side of the blade. e biggish pressure in the impeller is at outlet of the blade and it is relatively greater at 0.6Q d than that at 0.8Q d , Q d , and 1.2Q d . According to the law of performance curve, the head at 0.6 Q d is larger than that at 0.8Q d , Q d , and 1.2Q d . is may result in the biggish radial force at 0.6 Q d . e static pressure coefficient Cp (defined by (p-p in )/ (ρu 2 /2)) of monitoring points 1 and 2 near the outlet of the impeller marked in Figure 1

Radial Force and Vibration on Centrifugal Impeller.
Radial forces generated by unsteady flow are one of the important factors that affect operation stability, and these result in pump vibration. So the radial forces and vibrations on the impeller were investigated at 0.6Q d , 0.8Q d , Q d , and 1.2Q d .

Radial Force on the Impeller.
e time domain diagrams of the radial force in the X-direction, Y-direction, and resultant force on the impeller are shown in Figure 8 at 0.6Q d , 0.8Q d , Q d , and 1.2Q d .
As shown in Figure 8, the radial forces are periodic and vary with the time. e minimum radial force occurred at the design flow rate, while the radial forces increased, when the pump runs at off-design conditions. e radical force increased as the flow rate decreased. e radial forces and their fluctuations were relatively large at a flow rate of 0.6Q d .
eoretically, the velocity and pressure of fluid in the spiral volute around the impeller are uniform and axisymmetric, and the radial force on the impeller is zero at the design flow rate. In fact, the interaction between the volute tongue and impeller blade destroys the axisymmetric distribution of pressure along the circumference of the impeller. is results in the radial force on the impeller at the design flow rate. When the flow rate is less than the design flow rate, the fluid flow is equivalent to that in a diverging pipe relative to the design flow rate. e velocity of fluid in the volute becomes increasingly low from the tongue to the diffusion pipe; however, the pressure gradually increases. At a high flow rate, the opposite trend occurs relative to that of a low flow rate, which results in opposing directions of the radial force under the two conditions. As shown in Figures 8(a) and 8(b), the directions of the X-direction and Y-direction radial forces at flow rates of 0.6Q d and 0.8Q d are opposite to that at a flow rate of 1.2Q d for most time. As shown in Figure 8(c), the variation trend of the resultant force is similar at different flow rates, and five peaks and five troughs occurred in one cycle, which is the same as the number of blades. Owing to the interaction between the blade and tongue, the internal pressure fluctuates, which leads to variations in the radial forces.
e radial forces increased after the troughs and reached their peaks, when the blades sweep through the tongue. A possible explanation is that the pressure  fluctuation near the tongue increases as the distance between the tongue and blade decreases. e unbalance of the circumferential pressure leads to a higher radial force on the impeller. en, the radial forces rapidly decreased as the blades swept through the tongue. After passing the trough, pressures near the volute tongue gradually recovered in the pump and the radial forces reached their peak once again.  Figure 10. e amplitudes of the peak in the X-and Y-directions are mainly at the same, fourfold, sixfold, and ninefold shaft rotating frequency and same, twice, and thrice blade passing frequency. e amplitudes of the peak at Q d are less than that of other flow rates, when the frequency is less than triple blade passing frequency. ey are larger at blade passing frequency and frequency multiplication. e rotor-stator interaction between the blade and tongue has an important influence on vibration acceleration.

Vibration Experiment on the Volute and Vibration Energy.
Unsteady flow can cause severe energy loss in the centrifugal pump, where vibration energy as one of the main energy losses will severely affect the efficiency of the centrifugal pump.

Vibration Experiment on the Volute.
e pump is a magnetic force pump and vibration measurement cannot be conveniently carried out on the shaft. us, the vibration acceleration sensors indicated in Figure 1(b) were installed on the volute in the X-and Y-directions. e frequency range was 1-10 kHz and the mounted resonance frequency was 32 kHz. e vibration acceleration on the volute of the pump was measured at four different flow rates. e frequency domains of the vibration acceleration in the X-and Y-directions are shown in Figures 11(a) and 11(b), respectively. It can be observed in Figure 11 that there are peaks at 49.2 Hz, 147.5 Hz, 245 Hz, and 492.5 Hz in the X-and Y-directions, and these are the same and thrice the shaft rotating frequency and the same and twice the blade passing frequency, respectively. In general, the amplitude of vibration acceleration was relatively high at the shaft rotating frequency and at twice the blade passing frequency. In addition to vibrations induced by the fluid, there were also mechanical   vibrations during the experiment. us, the frequency domain of experimental results was more populated, because the pump is vertically installed under the platform and this is heavy and fixed on the base. Also the motor is mounted on the platform. e installation mode makes the vibration of volute very small. e most of dominant frequencies for vibration acceleration on the volute are at shaft frequency. But the most of dominant frequencies for vibration acceleration on the impeller are at blade passing frequency, which is same as the dominant one of pressure fluctuation for the monitoring point at the impeller outlet.

Vibration Energy Analyses.
To evaluate the varying trend of vibration energy in the frequency band, RMS method is applied to deal with discrete peaks at different frequencies, as defined in the following equation: where A n is the amplitude of the vibration acceleration at each frequency and n is the number of frequency points in each specific frequency band. Figure 12 presents the vibration energy of impeller and volute at different flow rates. e amplitude of the vibration acceleration on the impeller and volute mainly concentrated at the frequency range of 0-1000 Hz. From 0.6Q d to 1.0Q d , the vibration energy of impeller and volute decreases with the flow rate and reaches a larger value at 0.6Q d . And it is minimum at the design condition. And it increases at    Mathematical Problems in Engineering 1.2Q d with the increasing of flow rate. Under off-design flow rates, the inflow angle attack will change with the increase or decrease of the flow rates, which could cause vortex shedding from the surface of the blade. e unsteady flow has a significant effect on the pump vibration under off-design conditions.

Energy Loss of Centrifugal
Pump. e energy dissipation caused by the irreversible conversion of the available energy to unavailable energy is mainly manifested by the increase or decrease of entropy generation. e local entropy generation rate in the flow field is defined as follows [15]: Assuming an isothermal flow at a temperature T, u, v, and w are the time-averaged velocity components, μ is the eddy viscosity, and μ t is the dynamic viscosity. D is the contributions due to the time-averaged movement; D′ is the fluctuations of the velocity field. However, in the calculation of the Reynolds time-average method, the entropy generation caused by the fluctuations of the velocity field is not available. Kock [16] and Mathieu [17] proposed to correlate ε or ω in the turbulence model with the entropy generation by the fluctuations of the velocity field. In the calculation of the SST k-ω turbulence model, the local entropy generation due to the fluctuations of the velocity field can be expressed as k and ω are the turbulent kinetic energy and turbulent vortex frequency in the turbulence model, respectively. e constant β = 0.09. e total entropy generation (TEG) is the volume integration of the specific entropy generation rate, which can be expressed as Mathematical Problems in Engineering 9

Total Entropy Generation Rate Analyses.
e distributions of total entropy generation at different flow rates are shown in Figure 9, where Figure 13 e greater the total entropy generation value, the greater the energy loss and the lower the efficiency of the centrifugal pump. e change rule of vibration energy with flow rate on the impeller and the volute is consistent with that of total entropy generation. e vibration is a result of the poor flow conditions within the pump and it has an influence on the energy loss in centrifugal pump.
From Figure 13(a), it can be seen that the clearances of impeller side room result in greater energy loss and its proportion is about 40% at 1.0Q d . Also, the energy conversion of the volute leads to larger energy loss and its proportion is 37% at 1.2Q d . From Figure 13(b), it can be seen that the maximum TEG value of whole domain, volute domain, and impeller domain is 1.14 W/K, 0.35 W/K, and 0.29 W/K at 0.2Q d , respectively. is is due to a series of unsteady flow in the centrifugal pump at low flow rate conditions, especially near the volute tongue, resulting in large energy loss. e TEG value of the volute domain is higher than that of the impeller domain. From 0.2Q d to 1 e attack angle of inflow increases as the flow rate increases, also it decreases as the flow rate decreases. ey both cause impact loss to the surface of the blade.

Entropy Generation Rate Distribution.
e vibration is a result of the poor flow patterns and consequently high losses in the pump. e operating points of 0.2Qd and 1.0Qd are selected as objects for EGR and streamline analysis in the impeller. Figure 14 shows the EGR and streamline distribution in the impeller on different spanwise surfaces (Span 0.4, Span 0.6, and Span 0.8) at 0.2Qd and 1.0Qd. Span 0.4 is the surface close to the hub, and Span 0.8 is close to the shroud. At 0.2Q d , the energy loss of the impeller is mainly concentrated on the inlet and outlet of the blade. Because the inlet flow angle is less than the inlet install angle of blade, flow separation occurs at the front edge of the blade suction surface and gradually diffuses into the outlet of the impeller in the passage, resulting in a great separation loss. At the outlet of the blade, there is a great deal of energy loss generated by the backflow at the outlet of the impeller. As the flow rate increases, inflow attack angle gradually decreases. At 1.0Q d , the flow separation phenomenon and impact loss on the inlet of the blade become weakened, and the energy loss is drastically reduced.

Conclusions
Flow-induced impeller vibration and radial force in centrifugal pump were investigated through numerical analysis. And vibration energy and energy loss in the centrifugal pump were analyzed. e following conclusions can be drawn from the study: (1) e radial force and vibration displacement of the impeller were lowest at the design flow rate and increased at off-design flow rates. e variation trend of the radial displacements in the X-and Y-directions was consistent with that of the radial forces. e variation in the radial displacement lags that of the radial force.

Data Availability
No data were used to support this study.

Conflicts of Interest
e authors declare that they have no conflicts of interest.
Acknowledgments is work was supported by the Joint Project from the National Natural Science Foundation of China and Liaoning