Predictive Analysis of the Damage to Axial-Flow Pump ' s Impeller in Sandy Water

Hua HONG*, Zhang ZHI-ZHONG**, Liu XIAO-BING***, Zeng YONG-ZHONG****, Wang HUI-YAN***** *State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu 610065, Peoples R China, Key Laboratory of Fluid and Power Machinery (Xihua University), Ministry of Education, Chengdu 610039, Peoples R China, E-mail: huahong126@163.com **Guodian Dadu River Zhentouba Power Generation Co. LTD, Leshan, Peoples R China, E-mail: 653452406@qq.com ***Key Laboratory of Fluid and Power Machinery (Xihua University), Ministry of Education, Chengdu 610039, Peoples R China, E-mail: liuxb@mail.xhu.edu.cn ****Key Laboratory of Fluid and Power Machinery (Xihua University), Ministry of Education, Chengdu 610039, Peoples R China, E-mail: zyzzyzhome@163.com *****Key Laboratory of Fluid and Power Machinery (Xihua University), Ministry of Education, Chengdu 610039, Peoples R China, E-mail: 191758344@qq.com


Introduction
Axial flow pump is a low head and high flow pump, whose simple structure is available for urban water supply, agricultural irrigation and water diversion projects [1][2].In recent years, the ecological and environmental conditions are worsening, soil erosion is serious, and the sediment content in most rivers is higher and higher.At present, the designs of the axial flow pump impeller are mostly based on the premise of clean water medium, which leads to the serious damage of the impeller in the actual operation, reduces the operating efficiency of pumps and increases their energy consumption.This issue has attracted great attention from scholars, and has been studied theoretically and experimentally.Some results were obtained [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17].
Based on the Euler method, the 3D turbulence flow in the axial-flow pump impeller was numerically simulated by using the ANSYS-FLUENT, Eulerian multiphase flow model.The flow characteristics of solid-liquid two-phase flow were studied, and the results were compared with those in clear water.

Equations
The following assumptions have been made in this study: 1.The liquid phase (water) is incompressible.The solid phase (sand) is continuous.The physical properties of each phase are constants.
2. The solid phase consists of sand particles spherical in shape and uniform in size.
3. Neither the suspended matter nor the carrier liquid undergoes any phase change.
4. Interactions between particles, as well as between particles and the wall are neglected.

Basic equations of solid-liquid two-phase flow
The motion equations of solid-liquid two-phase flow in the Eulerian coordinate system are as follows [2]: Liquid phase continuity equation: Solid phase continuity equation: Liquid phase momentum equation: Solid phase momentum equation: Where: Ui and Vi are respectively the velocity components of the liquid and solid phase; ρ is the material density of phases; ν is the kinematic viscosity coefficient; P reprensents the gravity acceleration component; g is the component of gravity acceleration; and xi is the coordinate component.The equation of   , in which the subscripts L and S represent the liquid and solid phases, respectively, while i and k are the tensor coordinates.

Turbulence calculation model
The RNG k-ε model [2] was used in this study, which can be written as: Where: Gk is the turbulent kinetic energy caused by the average velocity gradient; Gb represents the turbulent kinetic energy caused by buoyancy influence; YM represents the influence of the total dissipation rate by the compressible turbulence fluctuation expansion; α k and αg are the reciprocals of effective turbulent Prandtl numbers of the turbulent kinetic energy k and dissipation rate ε.The calculation formula of turbulent viscosity coefficient is: where: , and Cν =100.And in Fluent, C1ε = 1.44,C2ε = 1.92,C3ε = 0.09 are all default constants.

Basic parameters of pump
The basic parameters of the pump are listed in Table 1.

Boundary conditions
The velocity inlet was used, and it was assumed that the inlet has uniform incoming flow, and the inlet velocity was vertical to the inlet boundary surface.The outlet was provided with flow outlet with the conditions of free development, that is to say, except the outlet pressure, the positive normal gradients of all flow variables were assumed as zero.The velocity on the solid wall can meet with the no-slip wall conditions, and the standard wall function was adopted for the near wall area.The average diameter of solid particles was 0.2 mm, and the density of solid phase material (sand in this study) was 2,719 kg/m 3 .The impact of gravity to the flow field during the calculation was considered and the direction of gravity was reverse to the normal of the pump outlet section.

Calculation conditions
The turbulent flow in an axial flow impeller was numerically simulated in the clear water (single-phase flow) and the sandy water (solid-liquid two-phase flow) conditions.According to the statistical data of the content of river sediment, three sediment concentrations were selected.The sand volume concentrations (fractions) of various conditions are listed in Table 2.

Calculation results and analysis
With the use of RNG k-ε turbulence model, SIM-PLEC algorithm and CFD software, the turbulent flow in an axial flow impeller was numerically simulated in the clear water (single-phase) and the sandy water (solid-liquid twophase) conditions.The distributions of solid concentration, velocity and pressure on the impeller of the axial flow pump were analyzed at the various conditions.

Static pressure distribution
Figs. 2-5 show the static pressure distributions on the leading side and the suction side of the blade in various conditions (including clear water and sandy water).In general, the pressure on the leading side is larger than that on the suction side.The pressure on the leading side of the blade is changing along the radius direction, and there is a small high pressure zone near the outlet flange and a small low pressure zone near the inlet flange.The pressure on the suction side is mainly changing in the circumferential direction, which reduces and then increases from the inlet to the outlet.There is a small high pressure zone near the outlet flange and a small low pressure zone near the inlet flange on the suction side.The minimum pressure on the blade is on the suction side near the inlet, which is the part most vulnerable to cavitation on the impeller.The minimum pressure in solid-liquid two-phase flow is lower than that in single-phase flow, which is shown in Fig. 6.The impeller is more vulnerable to cavitation in solid-liquid two-phase flow than that in single-phase flow.So the axial pump impeller will be damaged more easily and quickly in sandy water than in clear water.
With the increasing of the solid phase volume concentration in solid-liquid two-phase flow, the pressure gradually reduces.This indicates that in the solid-liquid twophase flow, due to the existence of the solid phase, the more the sediment content is, the more vulnerable of cavitation in the low pressure zone will be.This is consistent with the abrasion of impeller in practical projects.-9 show the volume concentration distributions of solid phase particles on the leading side and the suction side of blades in sandy water (solid-liquid two-phase flow) with various sediment concentrations.It shows that the distribution laws on the leading side and the suction side are different.In general, the solid phase concentration on the leading side of the blade is larger than that on the suction side.The solid phase concentration at the inlet edge is relatively higher on the leading side and the suction side.The abrasion on the part with higher solid phase concentration is serious, which is consistent with practical projects.There is a low concentration zone of solid phase on the leading side near the inlet and the flange, which is formed by the impact on the inlet wall as the pop-up of the solid phase particles with certain angle.The higher the solid phase volume concentration is, the higher the solid phase concentration on the blade surface will be, which is shown in Fig. 10.The higher the sediment concentration, the higher is the solid phase volume concentration near the blade surface, and the more serious of abrasion.
Under the combined effects of sand abrasion and cavitation, the extent of damage to the blades greatly intensified.This is consistent with the actual engineering situation.In the hydraulic design of this type of pump, some improvement measures should be taken to reduce the damage as much as possible.The figures show that the relative velocity near the blade surface is gradually increasing in radial direction from the hub to the flange.The relative velocity at the suction side is larger than that at the leading side.The relative velocity gradually decreases from the inlet to the outlet at the leading side and suction side.The relative velocity of the liquid phase in the solid-liquid two-phase flow is smaller than that in single-phase flow.This is caused by the solid phase which hinders the velocity of the liquid phase.Furthermore, with the increasing of sediment content, the velocity in the impeller reduces.

Conclusions
According to the simulation results, the pressure, velocity and the distributions of solid concentration on the impeller of an axial flow pump in different conditions are compared and analyzed.The following conclusions are obtained: 1.According to the simulation results, the axial pump impeller will be damaged more easily and quickly in sandy water than in clear water.
2. The area most prone to damage on the axial flow impeller is predicated, which is near the inlet edge on the leading side and the suction side of blades.Due to the minimum pressure on the blade near the inlet, the solid phase concentration at the inlet edge is relatively higher.Under the combined effects of sand abrasion and cavitation, the extent of damage to the blades is greatly intensified.
3. This study shows that the numerical simulation results are the same as the actual situation, and it has guiding significance for the wear design of the axial flow pump.In the hydraulic design of this type of pump, some improvement measures should be taken to reduce the damage as much as possible.

Acknowledgments
The

Fig. 1 3D
Fig. 1 3D model of the axial-flow pump impeller and grids: a -3D model of the axial flow impeller; b -flow passage meshing The entire flow passage consists of three parts, the tapered inlet stationary part, the impeller rotation part and the outlet stationary part.The 3D model of the axial-flow pump impeller is shown in the Fig. 1, a, and the direction of impeller rotation is forward along the Z-Axis.The entire

Fig. 6
Fig. 6 Max and minimum pressure in various conditions 4.2.Volume concentration distribution of solid phase particles Figs.7-9 show the volume concentration distributions of solid phase particles on the leading side and the suction side of blades in sandy water (solid-liquid two-phase flow) with various sediment concentrations.It shows that the distribution laws on the leading side and the suction side are different.In general, the solid phase concentration on the leading side of the blade is larger than that on the suction side.The solid phase concentration at the inlet edge is relatively higher on the leading side and the suction side.The abrasion on the part with higher solid phase concentration is serious, which is consistent with practical projects.There is a low concentration zone of solid phase on the leading side

Fig. 10
Fig. 10 Max and minimum Solid phase volume concentration in various conditions 4.3.Velocity distribution Figs.11-14 show the distributions of relative velocity of the liquid phase (water) near the surface of the axial-flow pump blades (hereinafter referred to as blade surface).The figures show that the relative velocity near the blade surface is gradually increasing in radial direction from the hub to the flange.The relative velocity at the suction side is larger than that at the leading side.The relative velocity gradually decreases from the inlet to the outlet at the leading side and suction side.The relative velocity of the liquid phase in the solid-liquid two-phase flow is smaller than that in single-phase flow.This is caused by the solid phase which hinders the velocity of the liquid phase.Furthermore, with the increasing of sediment content, the velocity in the impeller reduces.

Fig. 11
Fig. 11 Liquid phase relative velocity distribution near blades in condition 1 (clear water) : a -leading side; b -suction side work described in this paper is supported by the National Key Research and Development Program "Research and Application Demonstration of Complementary Combined Power Generation Technology between Distributed Photovoltaic and Cascade Small Hydropower" (Grant No.2018YFB0905200), the National Natural Science Foundation of China (Grant No. 51279172), the Fundamental Research Funds for the Central Universities (Grant No. 2012017yjsy169) and the Open Research Fund of the Key Laboratory of Fluid and Power Machinery, Ministry of Education (Xihua University) (Grant No. szjj2017-087).

Table 2
Statistics of CFD calculation conditions