Next Article in Journal
Clionaterpene, a New Cadinene Sesquiterpene from the Marine Sponge Cliona sp.
Next Article in Special Issue
CFD Simulation of Centrifugal Pump with Different Impeller Blade Trailing Edges
Previous Article in Journal
Effects of Anticyclonic Eddies on the Unique Tropical Storm Deliwe (2014) in the Mozambique Channel
Previous Article in Special Issue
Effect of Flow Rate on Regular Patterns of Pressure Load Distribution on Helico-Axial Pump Impeller Blade Surface
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analyzing the Inlet Gas Void Fraction on the Flow Characteristics for a Multiphase Pump Based on Cavitation Cases

1
Key Laboratory of Fluid and Power Machinery, Xihua University, Ministry of Education, Chengdu 610039, China
2
Key Laboratory of Fluid Machinery and Engineering, Xihua University, Chengdu 610039, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2023, 11(1), 130; https://doi.org/10.3390/jmse11010130
Submission received: 4 December 2022 / Revised: 15 December 2022 / Accepted: 16 December 2022 / Published: 6 January 2023
(This article belongs to the Special Issue CFD Simulation of Floating Offshore Structures)

Abstract

:
Inlet gas void fraction (IGVF) affects the cavitation evolution in a multiphase pump and easily results in a drop of the head and efficiency when cavitation is more serious. In this paper, a numerical method was performed to qualitatively and quantitatively analyze the effect of the inlet gas void fraction on the pressure and velocity characteristics of the multiphase pump at different cavitation stages. The results show that with the increase of IGVF and the development of cavitation, the pressure in the impeller flow passage is reduced, and the pressure corresponding to the cavitation region drops sharply to the saturated vapor pressure. With the decrease of the cavitation coefficient, and due to the expulsion effect of the cavitation bubbles, the relative velocity in the cavitation region becomes larger. Because of the large pressure gradient at the end of the cavitation bubbles, the kinetic energy of the fluid is insufficient to overcome the effect of the inverse pressure gradient, resulting in a backflow vortex. Investigations on cavitation evolution in the multiphase pump at different IGVFs are of great significance for improving its performance.

1. Introduction

With the continuous expansion of the global energy demand, the exploitation of oil and gas in human society is increasing, and multiphase transportation technology is also emerging [1,2,3]. As the core piece of equipment of the multiphase transportation system, the multiphase pump has the advantages of large flow, small volume and insensitivity to solid particles [4,5]. Therefore, many scientific research institutions, universities and oil companies have carried out plenty of investigations on it [6,7,8,9]. However, the inlet gas void fraction (IGVF) always changes during the operation of the pump, which gives rise to the variation of the flow field. In some cases, a low-pressure region appears, which easily leads to cavitation in the multiphase pump [10,11,12,13]. Cavitation seriously affects the operation of the multiphase pump, and it is of great engineering significance to investigate the effect of an IGVF on the flow characteristics in a multiphase pump in a cavitation case.
At present, the investigations on the flow characteristics of multiphase pumps mainly focus on the internal flow behaviors and structural optimization design, without considering the effect of cavitation. Yu Zhiyi et al. [14] employed numerical and experimental methods to investigate the gas distribution in the multiphase pump and the effect of the IGVF on the head under the gas-liquid two-phase condition. Zhang Jinya et al. [15] took water and air as the medium and conducted steady and unsteady numerical simulations of a multistage multiphase pump. When the IGVF was less than 10% and greater than 90%, the gas was regarded as an incompressible fluid; at other IGVFs, the gas was regarded as a compressible fluid. Yang Xiao qiang et al. [16] conducted an experiment on the external characteristics of the twin-screw multiphase pump under different pressure differences and the IGVFs. The relationship between output power and IGVF and the inlet pressure difference were studied after analyzing the experimental data. Ma Xijin et al. [17] carried out a three-dimensional numerical simulation of the flow characteristics of a multiphase pump at different IGVFs, and investigated the effect of impeller blade number on multiphase pump performance. The results showed that the appropriate increase in blade number increased the pump head. Yang, X. et al. [18] adopted N32 oil and air as the medium to conduct theoretical and experimental investigations on the multiphase pump at high IGVF. The experimental data were obtained by changing the IGVF and inlet pressure. The leakage in the multiphase pump at a given pressure was increased with the increase of the IGVF, and the shaft power wasted via the high-pressure reflux reduced the pump’s efficiency. Zhang, J.Y. et al. [19,20] used a multi-stage multiphase pump as the research object and adopted PIV and high-speed photography to capture the gas-liquid two-phase flow at the pump inlet. As the IGVF increased, the gas-liquid flow patterns were divided into isolated bubble flow, bubbly flow, gas pocket flow and segregated gas flow. Shi, Y. et al. [21] compared the experimental data of the three-stage multiphase pump with the numerical results under different turbulence models, and uncovered the effect of the turbulence model, wall roughness, bubble size and interphase resistance model on the multiphase pump’s performance. The simulation results with the SST turbulence model were more consistent with the experimental case. Kim, J.H. et al. [22] took the multiphase pump efficiency as the objective function for optimization. By defining the diffusion-area ratio, the straight-blade length ratio, and the length ratio between the trailing edge (TE) of the impeller blade and the leading edge (LE) of the diffuser blade, a radial basis function neural network was used to optimize the multiphase pump blade. The best efficiency point of the multiphase pump was increased by 9.75%, and the multiphase pump efficiency was also improved at the large flow rate. By establishing a more appropriate and reliable numerical analysis method, Suh, J.W. et al. [23] revealed the effect of the IGVF on the flow characteristics of the multiphase pump. Zhang, W.W. et al. [24] investigated the pressure fluctuation characteristics of the multiphase pump under the gas-liquid two-phase case, and disclosed the relationship between the IGVF and the pressure fluctuation. The low IGVF had an inhibitory effect on the pressure pulsation caused by interference at a certain degree. When the IGVF reached a certain level, the amplitude of the pressure fluctuation became larger due to the enhancement of the interphase effect. In terms of the thermal effect of cavitation, Ge Mingming, et al. [25,26,27,28] studied the influence law of temperature on cavitation dynamics in a Venturi tunnel with experimental methods, divided the law of cavitation evolution, and studied the influence of temperature on the cavitation structure. The research deepened scholars’ understanding of the cavitation law and has a very important significance. In terms of optimal design, Shi, G. et al. [29] used orthogonal optimization and the CFD method to optimize the gas-liquid conveying capacity of a multiphase pump. When the IGVF was 15%, the optimized pump head and efficiency were increased by 2.81% and 5.6% respectively. By considering the pump head and efficiency, Wang C et al. [30] established an optimization numerical model and used an artificial intelligence algorithm to optimize, so the pump performance had been improved. Nguyen, V et al. [31] proposed some design methods for centrifugal pumps that had great effects on the improvement of pump performance. Similarly, Zhang, W.W. et al. [32], Kim, H. et al. [33], Liu, M. et al. [34] and Peng, C. et al. [35] adopted different optimization methods to optimize the performance of the multiphase pump and improve the working capacity of the pump.
To sum up, there are currently few investigations on the flow characteristics of multiphase pumps in cavitation case. Meanwhile, the effect of IGVF on the flow characteristics of multiphase pumps under cavitation is seldom reported. Therefore, under different cavitation stages, the effect of IGVF on the flow characteristics of multiphase pumps was qualitatively and quantitatively investigated, which had important scientific significance in improving the flow stability of multiphase pumps.

2. Computational Model and Method of Multiphase Pump

2.1. Computational Model

An impeller was selected to investigate the cavitation characteristics of a multiphase pump in this paper, and the computational domain was composed of inlet pipe, impeller pipe, and outlet pipe. To achieve full flow at the impeller inlet and outlet, the inlet and outlet were respectively extended to 2 and 6 times the axial length of the impeller, as shown in Figure 1. The design parameters of the pump are as follows: design flow rate Q = 100 m3·h−1, speed n = 3000 rpm, impeller blade number Z = 4, impeller blade wrap angle 179.6°, impeller inlet hub ratio 0.79, impeller outlet hub ratio 0.74, and impeller diameter 230.5 mm.

2.2. Numerical Simulation Theory

2.2.1. Turbulence Model

The multiphase pump blade has a large curvature and high speed, and flow separation occurs in the multiphase media case. Therefore, the SST k-ω turbulence model was selected in this paper.
The k-epsilon turbulence model has weak results in predicting strongly separated flows, flows containing a large curvature and flows with a strong pressure gradient. The SST k-ω turbulence model combines the advantages of the k-ω model and the k-ε model, which better deals with the viscous flow in the near-wall region and the turbulent flow in the far field. The expressions of the turbulent kinetic energy k and the turbulent pulsation frequency w of the SST k-ω model are as follows:
ρ k t + ρ v j k x j = x j u + u t σ k k x j + u t v j x i v j x i + v i x j β * ρ k w
ρ w t + ρ v j w x j = x j u + u t σ w w x j + α w k u t v j x i v j x i + v i x j β ρ k w 2 + 2 1 F 1 ρ 1 w σ w k x j w x j
where ρ is the density, β , β * and σ w are the empirical coefficient, and F1 is the blending function.

2.2.2. Multiphase Flow Model

In this paper, the range of IGVF in the multiphase pump was 0~20%, and the volume fraction of bubbles varied greatly and its distribution range was wide. The cavitation flow was more complicated, and the interface between the vapor and liquid phases was not clearly defined. Considering the computational cost, the mixture model was selected to perform the numerical simulation of the multiphase flow in the multiphase pump.
The mixture model is used commonly for multiphase flow in engineering, which allows two phases to intersect each other. The volume fraction of the two phases in a control body can be any value between 0 and 1, and there is a speed slip between the two phases. The expressions of the continuity equation, momentum equation, energy equation, the volume fraction equation of the second phase and the relative velocity in the mixture model are as follows:
t ρ m + ρ m V m = 0
t ρ m V m + ρ m V m V m = p + μ m V m + V m T + ρ m g + F + k = 1 n α k ρ k V d r , k V d r , k
t k = 1 n α k ρ k E k + k = 1 n α k V k ρ k E k + ρ = k e f f T + S E
t α p ρ p + α p ρ p V m = α p ρ p V d r , p
μ m = k = 1 n α k μ k
k = 1 n α k = 1
ρ m = k = 1 n α k ρ k
V m = k = 1 n α k ρ k V k ρ m
Drift velocity of the k phase is
V d r , k = V k V m
where ρ m and ρ k are the density of the mixture and k phase; V d r , k , V k and V m are the drift velocity of the k phase, the velocity of the k phase and the average velocity of the mass; α k is the volume fraction of the k phase; p is the pressure; μ m and μ k are the dynamic viscosity of the mixture phase and k phase; g is the acceleration of gravity; F is body force; k e f f is the effective thermal conductivity; E k is the energy of the k phase; T is the temperature; and SE includes the contribution of all other volumetric heat sources.
The mixture model is a heterogeneous model that uses slip velocity to allow for a different velocity between phases. The slip velocity is also the relative velocity, which refers to the velocity of the secondary phase (p) with respect to the primary phase (q)
V q p = V p V q
The slip velocity is expressed as follows:
V q p = τ p f d r a g ρ p ρ m ρ p a
where a is the acceleration of the secondary phase (gas or vapor in cavitation); τ p is the relaxation time of the p phase; and f d r a g is the drag function.
a = g ( V m · ) V m V m t
τ p = ρ p d p 2 18 μ q
where d p is the p phase particle diameter.

2.2.3. Cavitation Model

The simulation results based on the Singhal cavitation model have a large deviation from the experimental data and over-predict the range of cavitation. The Schnerr–Sauer cavitation model predicts that the evolution period of the cavitation flow is smaller than the experimental value, and cannot accurately predict the evolution and development of cavitation. The Zwart–Gerber–Belamri cavitation model accurately simulates the quasi-periodical and evolution process of cavitation. Therefore, the Zwart–Gerber–Belamri cavitation model was used to simulate the cavitation flow in the multiphase pump. The interphase transmission rate of this model is as follows:
when P P v
R e = F v a p 3 α n u c 1 α v ρ v R B 2 P v P 3 ρ l
when P > P v
R c = F c o n d 3 α v ρ v R B 2 P P v 3 ρ l
In the formula, RB is the bubble radius, 10−6 m; αnuc is the volume fraction of the vapor nucleus position, 5 × 10−4; Fvap and Fcond represent the vapor evaporation and condensation coefficients, and the values are 50 and 0.01, respectively.

2.3. Mesh and Independent Verification

The structural hexahedral mesh was adopted on the single-flow passage of the pump impeller, and ICEM CFD software was performed to rotate and copy the single-flow passage mesh into a full-flow passage. Meanwhile, the ICEM software was also used to arrange structural hexahedral mesh on the inlet and outlet extensions. The advantage of structural hexahedral mesh is that the mesh nodes can be adjusted to optimize the local mesh, and the flow details can be clearly presented. The meshes of the inlet pipe, outlet pipe and impeller are shown in Figure 2, and the mesh quality is shown in Table 1.
To improve the accuracy and efficiency of the numerical simulation, the mesh independence was verified. A total of 7 sets of meshes on the multiphase pump were performed and simulated at water case and design flow rate. The appropriate mesh number was selected for the final numerical simulation. The mesh independent verification is shown in Figure 3.
From Figure 3, as the mesh number increases, the head gradually decreases and stabilizes. When the mesh number is greater than 3.38 million, the head changes only by 0.31%, which is less than 0.5%. Therefore, 3.38 million mesh was more appropriate in the final simulation. The meshes of the inlet pipe, impeller and outlet pipe were 500,000, 2,170,000 and 900,000, respectively.

2.4. Boundary Condition Settings

ANSYS Fluent software was used to simulate the steady cavitation flow in the multiphase pump at the design case. The turbulence model, multiphase flow model and cavitation model were employed as described above. The inlet and outlet boundaries were set to pressure inlet and mass flow, respectively. The IGVF was set to 0 at the inlet, and the saturated vapor pressure of water at 25 °C was 3170 Pa. Cavitation occurs in the pump by gradually reducing the pressure at the multiphase pump inlet. The interfaces between the impeller and the inlet pipe, the impeller and outlet pipe are set as the interface. The impeller hub, blade and shroud are set as relatively no-slip walls, and the other walls are absolutely no-slip walls. The convergence accuracy is set to 10−5.

3. Experimental Rig and Numerical Verification

3.1. Experimental Rig

The experimental rig of the multiphase pump included the motor, multiphase pump, gas-liquid mixing tank, lubrication system, cooling system, control system, water supply system, gas supply system, test system, pipeline and valves, etc. The experimental system for the multiphase pump is shown in Figure 4.
In the experiment, the inlet and outlet pressure gauges were used to measure the inlet and outlet pressure of the multiphase pump. A torque meter was used to measure speed, power and torque. Repeated tests were carried out for each operating point to prevent errors in the test process.

3.2. Numerical Verification

Figure 5 shows the comparison between the CFD and the experiment, and the head and efficiency in the CFD agreed well with the experimental one. Meanwhile, the relative errors of the head and efficiency at the optimal point were 4.1% and 4.1%, respectively. The maximum relative errors of the head and the efficiency were both no more than 5%, and this shows that the simulation is reliable.

4. Results and Discussion

4.1. Cavitation Characteristic Prediction

The head coefficient ψ and cavitation coefficient σ are shown as follows:
ψ = P o u t P i n 0.5 ρ U 2
σ = P i n P v ρ U 2 / 2
where Pin and Pout denote the multiphase pump inlet and outlet pressure, respectively. U is the circumferential velocity of the impeller hub, expressed as U = π D n / 60 . Pv is the saturated vapor pressure of water at 25 °C, 3170 Pa.
Figure 6 shows that the head coefficient ψ of the multiphase pump varies with the cavitation coefficient σ. For Figure 6, under IGVF = 0, 0.1 and 0.2, with the decrease of the cavitation coefficient σ, the head coefficient ψ remains unchanged at first, then steadily decreases, and finally decreases sharply.
Under IGVF = 0, when the cavitation coefficient σ is greater than 0.86, the head coefficient ψ remains unchanged, indicating that there is no cavitation, or the cavitation in the multiphase pump is relatively weak and does not affect the head. When the cavitation coefficient σ is between 0.106~0.86, the head coefficient ψ gradually decreases. When the cavitation coefficient σ is 0.28, the head coefficient ψ is reduced by 3% compared to the no-cavitation case, which is called critical cavitation in engineering. When the cavitation coefficient σ is less than 0.106, the head coefficient ψ decreases sharply. When σ decreases to 0.077, the head coefficient decreases by 7.68%, which indicates that the degree of cavitation in the multiphase pump is more serious. When the cavitation coefficient σ is 0.051, the head reduction is exceeded by 20%, and the cavitation reaches the fractured cavitation.
In the case of IGVF = 0.1, the cavitation coefficient σ is greater than 0.86, and the head coefficient ψ remains unchanged. When the cavitation coefficient σ gradually decreases between 0.0769 and 0.86, the head coefficient ψ decreases steadily, and the head coefficient ψ corresponding to the critical cavitation is 0.24. When the cavitation coefficient σ is less than 0.0769, the head coefficient ψ drops sharply. When the head coefficient ψ decreases by 7.3%, the cavitation coefficient σ is 0.057. The cavitation coefficient σ corresponding to a decrease of more than 20% in the head coefficient ψ is 0.033.
Under IGVF = 0.2, the head coefficient ψ remains unchanged in the range of cavitation coefficient σ greater than 0.86. When the cavitation coefficient σ gradually decreases in the range of 0.107~0.86, the head coefficient ψ of the multiphase pump is reduced steadily, and the critical cavitation coefficient σ is 0.208. When the cavitation coefficient σ gradually decreases in the range of less than 0.107, the head coefficient ψ decreases greatly. When the cavitation coefficient σ is 0.048, the head coefficient ψ decreases by 7.4%; when the cavitation coefficient σ is 0.029, the head coefficient ψ decreases by more than 20%.

4.2. Effect of IGVF on Pressure Characteristics at Cavitation Case

To analyze the effect of the cavitation development on the pressure characteristics in the multiphase pump impeller, the pressure was analyzed under the different IGVFs and cavitation states. Figure 7 shows that the pressure in the impeller passage at 0.1, 0.5, and 0.9 span IGVF = 0, 0.1, and 0.2 under different cavitation stages.
Under different operating points, because the impeller blade works on the fluid, the pressure in the flow passage gradually increases from the impeller inlet to the outlet; as the span increases, the pressure gradually rises. In the critical cavitation case, with the increase of the IGVF, the pressure in the impeller flow passage at the same span gradually decreases; that is, the increase of the IGVF reduces the pressurization performance of the multiphase pump. When the cavitation develops to the second stage, the pressure in the impeller flow passage becomes smaller than that during critical cavitation, and the pressure drop corresponding to the cavitation region is larger. When the cavitation develops to the third stage, the degree of cavitation is very serious, and cavitation appears at the suction side (SS) and pressure side (PS). The pressure in the cavitation region is very small. In the IGVF = 0 case, the cavitation occupies the entire blade SS, and the pressure in this region is almost 0.
To explore the effect of cavitation development on the pressure on the blade surface, the pressure on the 0.5 blade span was analyzed. Figure 8 displays the pressure along the streamwise at the 0.5 blade span under different cavitation stages. The pressure on the PS is higher, and the pressure on the SS is lower. The area enclosed by the pressure between the PS and SS can reflect the blade load.
In Figure 8a, when the cavitation is weak, the pressure on the PS increases steadily from the inlet to the outlet along the streamwise. However, the pressure on the SS first maintains 80 kPa in the streamwise of 0~0.3, then increases, and finally remains at a steady increasing trend. When the cavitation coefficient is reduced to 0.077, the pressure on the blade SS is basically unchanged from the critical cavitation. Meanwhile, the pressure on the blade SS in the streamwise of 0~0.32 is about 3170 Pa. When the cavitation coefficient σ is reduced to 0.051, the cavitation on the blade PS extends to the streamwise of 0.05~0.62, and the SS is completely occupied by cavitation. Then, the blade load in the streamwise of 0.05~0.62 becomes 0, and the blade load from the streamwise of 0.62 to the blade outlet becomes larger.
From Figure 8b,c, with the increase of IGVF, in the critical cavitation the pressure on the blade PS at the same location gradually decreases, and the blade load gradually decreases. When the cavitation is in the second stage—that is, σ = 0.057 and 0.048—since cavitation occurs in the streamwise of 0~0.35 on the blade SS, the pressure on the blade SS is reduced to about the saturated vapor pressure. When the cavitation coefficient σ is reduced to 0.033 and 0.029, the cavitation appears in the streamwise of 0.05~0.38 and 0~0.7 on the blade PS and SS, and the pressure on the blade in this range is reduced to about 3170 Pa. The blade load in the cavitation intersects between the PS and SS (in the streamwise of 0.05~0.38) decreases to 0. However, the blade load in the streamwise of 0.38~0.7 gradually increases. In addition, in the second and third stages of cavitation at different IGVFs, the pressure gradient corresponding to the end of the cavitation on the blade surface is relatively large.
Based on the above analysis, the increase of the IGVF makes the pressure in the impeller flow passage decrease, and the blade load decreases. The pressure corresponding to the bubble region is reduced to about the saturated vapor pressure.

4.3. Effect of IGVF on Velocity Characteristics at Cavitation Case

To analyze the effect of cavitation evolution on the flow pattern in the impeller, the velocity in the impeller passage at different spans is analyzed below. Figure 9 displays the relative velocity in the impeller at different spans with IGVF of 0, 0.1, and 0.2 under different cavitation states.
In Figure 9, in the critical cavitation at different IGVFs, the kinetic energy and pressure energy of the fluid are mutually converted during the fluid flows from the impeller inlet to the outlet, resulting in a gradual decrease in velocity. The relative velocity also gradually decreases from the hub to the shroud. When the cavitation develops to the second stage, the relative velocity in the cavitation region becomes larger, and as the span increases, the relative velocity gradually decreases. When the cavitation develops to the third stage, the cavitation in the impeller is more serious, and the cavitation range increases. The relative velocity in the cavitation region increases more so than in the previous two cavitation stages. Compared to the first two cavitation stages, as the blade span increases, the relative velocity in the cavitation region gradually increases; this indicates that the flow state in the impeller changes greatly, and the velocity distribution from the end of the cavitation to the outlet is more turbulent when the cavitation is severe. In addition, there is a low-velocity region at the end of the cavitation zone. From the pressure distribution, this is because the pressure gradient at the end of the cavitation region is large, which results in a larger velocity reduction.
To quantitatively analyze the effect of the development of cavitation on the velocity distribution under the different IGVFs, the velocity distribution at 0.5 span along the streamline direction is analyzed below. Figure 10 shows the velocity distribution at 0.5 span along the streamwise with IGVF of 0, 0.1 and 0.2 under different cavitation stages.
From Figure 10, at different IGVFs, when cavitation develops to the second stage, compared to the critical cavitation, the relative velocity gradient at blade SS becomes larger in the streamwise of 0~0.1. From the previous distribution of bubbles, in this cavitation stage, the bubbles are attached to the blade SS, which thickens the boundary layer around the inlet of the blade SS and results in a large change in the velocity gradient. However, in the streamwise of 0.3 to 0.4, the velocity gradient changes greatly compared to the critical cavitation, which is caused by the larger pressure gradient at the end of the bubbles on the blade surface. When the cavitation develops to the third stage, at IGVF = 0 the cavitation is very serious, and the cavitation occupies the entire blade SS and most of the impeller passage; this leads to an increase in the velocity on the blade SS from the streamwise of 0.3 to the outlet. For IGVF = 0.1 and 0.2, because the cavitation phenomenon extends along the blade SS to the streamwise of 0.7, the pressure gradient in the streamwise of 0.6 to 0.8 becomes larger, resulting in a larger change in the velocity gradient. Under different IGVFs, the large velocity gradient on the PS is also caused by the large pressure gradient at the end of the bubbles. Due to the squeezing effect of the bubbles, the relative velocity outside of the bubble area increases.
To observe the flow pattern in the impeller passage more intuitively, the velocity streamlines at 0.5 span are used to represent the fluid behaviors. Figure 11 is the velocity streamlines at 0.5 span with IGVF = 0, 0.1 and 0.2 under different cavitation states.
From Figure 11, in the critical cavitation, the velocity streamlines in the flow passage are relatively smooth at different IGVFs. When the cavitation develops to the second stage, due to the large pressure gradient at the end of the bubbles, a backflow vortex appears at the end of the bubbles on the blade SS. When the cavitation develops to the third stage, at IGVF = 0 the backflow vortex appears at the end of the bubbles on the blade PS due to the large reverse pressure gradient. However, a backflow vortex appears at the end of the bubbles on the PS and SS with IGVF = 0.1 and 0.2.

5. Conclusions

The effect of cavitation flow in the multiphase pump was investigated through the analysis of the pressure, velocity and turbulent kinetic energy characteristics. The main conclusions are as follows:
(1)
The increase of IGVF and the development of cavitation reduce the pressure in the impeller flow passage, and the pressure corresponding to the bubbles drops sharply to about the saturated vapor pressure. When the bubbles extend to the blade PS, the load near the cavitation intersection area is reduced to 0. The load on the end of the bubbles between the PS and SS increases, and the pressure gradient at the end of the bubbles is very large.
(2)
As the cavitation coefficient decreases, the relative velocity near the cavitation becomes larger due to the squeezing effect of the bubbles. Due to the large pressure gradient at the end of the bubbles, the kinetic energy of the fluid is not enough to overcome the effect of the reverse pressure gradient, resulting in a backflow vortex.
(3)
As the blade span increases, the relative velocity in the cavitation region gradually increases. Because the pressure gradient at the end of the cavitation region is large, a low-velocity region occurs at the end of the cavitation zone. To further avoid the cavitation phenomenon, a convex structure can be set on the blade to improve the flow state in the pump.

Author Contributions

Conceptualization, W.L. and H.W.; methodology, W.L.; software, W.L., H.W. and G.S.; writing—original draft preparation, H.W.; writing—review and editing, W.L. and H.W.; supervision, S.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Open Research Fund Program of the State Key Laboratory of Hydroscience and Engineering (sklhse-2021-E-03, sklhse-2022-KY-06); the Key scientific research fund of Xihua University of China (Z1510417); the Central Leading Place Scientific and Technological Development Funds for Surface Project (2021ZYD0038); the National Key Research and Development Program (2018YFB0905200); and The National Natural Science Foundation of China (52279088).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Liu, M.; Tan, L.; Xu, Y.; Cao, S. Optimization design method of multi-stage multiphase pump based on Oseen vortex. J. Pet. Sci. Eng. 2020, 184, 106532. [Google Scholar] [CrossRef]
  2. Liu, M.; Tan, L.; Cao, S. Influence of viscosity on energy performance and flow field of a multiphase pump. Renew. Energy 2020, 162, 1151–1160. [Google Scholar] [CrossRef]
  3. Sano, T.; Wakai, T.; Reclari, M.; Xu, Y.; Cao, S. Investigation of internal flow pattern of a multiphase axial pump. IOP Conf. Ser. Earth Environ. Sci. 2019, 240, 062058. [Google Scholar] [CrossRef]
  4. Shi, G.; Wang, S.; Xiao, Y.; Liu, Z.; Li, H.; Liu, X. Effect of cavitation on energy conversion characteristics of a multiphase pump. Renew. Energy 2021, 177, 1308–1320. [Google Scholar] [CrossRef]
  5. Wang, C.; Zhang, Y.; Zhang, J.; Zhu, J. Flow pattern recognition inside a rotodynamic multiphase pump via developed entropy production diagnostic model. J. Pet. Sci. Eng. 2020, 194, 107467. [Google Scholar] [CrossRef]
  6. Zhang, Y.; Zhu, J.; Zhu, H.; Cai, S.; Wang, J. 3D Blade Hydraulic Design Method of the Rotodynamic Multiphase Pump Impeller and Performance Research. Adv. Mech. Eng. 2014, 6, 803972. [Google Scholar] [CrossRef]
  7. Zhang, J.; Cai, S.; Zhu, H.; Zhang, Y. Experimental investigation of the flow at the entrance of a rotodynamic multiphase pump by visualization. J. Pet. Sci. Eng. 2015, 126, 254–261. [Google Scholar] [CrossRef]
  8. Xu, Y.; Cao, S.; Sano, T.; Wakai, T.; Reclari, M. Experimental Investigation on Transient Pressure Characteristics in a Helico-Axial Multiphase Pump. Energies 2019, 12, 461. [Google Scholar] [CrossRef] [Green Version]
  9. Han, W.; Li, X.; Su, Y.; Su, M.; Li, R.; Zhao, Y. Effect of Thickness Ratio Coefficient on the Mixture Transportation Characteristics of Helical–Axial Multiphase Pumps. Appl. Sci. 2020, 10, 345. [Google Scholar] [CrossRef] [Green Version]
  10. Liu, X.; Hu, Q.; Wang, H.; Jiang, Q.; Shi, G. Characteristics of unsteady excitation induced by cavitation in axial-flow oil–gas multiphase pumps. Adv. Mech. Eng. 2018, 10, 168781401877126. [Google Scholar] [CrossRef]
  11. Huan, Y.-Y.; Liu, Y.-Y.; Li, X.-J.; Zhu, Z.-C.; Qu, J.-T.; Zhe, L.; Han, A.-D. Experimental and numerical investigations of cavitation evolution in a high-speed centrifugal pump with inducer. J. Hydrodyn. 2021, 33, 140–149. [Google Scholar] [CrossRef]
  12. Lomakin, V.; Bibik, O. Numerical prediction of the gas content effect on the cavitation characteristics of the pump using the simplified Rayleigh-Plesset equation. IOP Conf. Ser. Mater. Sci. Eng. 2019, 492, 012037. [Google Scholar] [CrossRef]
  13. Li, J.; Liu, L.-J.; Feng, Z.-P. Two-dimensional analysis of cavitating flows in a centrifugal pump using a single-phase Reynolds averaged Navier—Stokes solver and cavitation model. Proc. Inst. Mech. Eng. Part A J. Power Energy 2006, 220, 783–791. [Google Scholar] [CrossRef]
  14. Yu, Z.; Zhang, Q.; Huang, R.; Cao, S. Numerical simulation of unsteady flow in multiphase rotodynamic pumps. J. Drain. Irrig. Mach. Eng. 2013, 31, 284–288. [Google Scholar]
  15. Zhang, J.; Cai, S.; Zhu, H.; Yang, K.; Qiang, R. Numerical Investigation of Compressible Flow in a Three-stage Helico-axial Multiphase Pump. Trans. Chin. Soc. Agric. Mach. 2014, 45, 89–95. [Google Scholar]
  16. Yang, X.; Xia, Y.; Jin, L.; Cao, F. Experimental Study on Pumping Behavior of Twin-Screw Multiphase Pump. J. Xi’an Jiaotong Univ. 2013, 47, 30–35. [Google Scholar]
  17. Ma, X.; Zhao, J.; Shao, L. Study on the Performance of Oil-gas Pump by Varying the Blade Number. Fluid Mach. 2009, 37, 24–27. [Google Scholar]
  18. Yang, X.; Hu, C.C.; Qu, Z.C. Theoretical and experimental study of a synchronal rotary multiphase pump at very high inlet gas volume fractions. Appl. Therm. Eng. 2017, 110, 710–719. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, J.; Cai, S.; Li, Y. Visualization study of gas–liquid two-phase flow patterns inside a three-stage rotodynamic multiphase pump. Exp. Therm. Fluid Sci. 2016, 70, 125–138. [Google Scholar] [CrossRef]
  20. Zhang, J.; Cai, S.; Zhu, H. Experimental Study of Gas-Liquid Two-Phase Flow Pattern in a Helico-Axial Multiphase Pump by Visualization. J. Eng. Thermophys. 2015, 36, 1937–1941. [Google Scholar]
  21. Shi, Y.; Zhu, H.; Zhang, J. Experiment and numerical study of a new generation three-stage multiphase pump. J. Pet. Sci. Eng. 2018, 169, 471–484. [Google Scholar] [CrossRef]
  22. Kim, J.H.; Kim, K.Y. Analysis and Optimization of a Vaned Diffuser in a Mixed Flow Pump to Improve Hydrodynamic Performance. J. Fluids Eng. 2012, 134, 071104–071113. [Google Scholar] [CrossRef]
  23. Suh, J.W.; Kim, J.W.; Choi, Y.S.; Kim, J.H.; Joo, W.G.; Lee, K.Y. Development of numerical Eulerian-Eulerian models for simulating multiphase pumps. J. Pet. Sci. Eng. 2018, 162, 588–601. [Google Scholar] [CrossRef]
  24. Zhang, W.; Yu, Z.; Zhu, B. Numerical Study of Pressure Fluctuation in a Gas- Liquid Two-Phase Mixed-Flow Pump. Energies 2017, 10, 634. [Google Scholar] [CrossRef] [Green Version]
  25. Ge, M.; Zhang, G.; Petkovek, M.; Long, K.; Coutier-Delgosha, O. Intensity and regimes changing of hydrodynamic cavitation considering temperature effects. J. Clean. Prod. 2022, 338, 130470. [Google Scholar] [CrossRef]
  26. Ge, M.; Petkovek, M.; Zhang, G.; Jacobs, D.; Coutier-Delgosha, O. Cavitation dynamics and thermodynamic effects at elevated temperatures in a small Venturi channel. Int. J. Heat Mass Transf. 2021, 170, 120970. [Google Scholar] [CrossRef]
  27. Ge, M.; Sun, C.; Zhang, G.; Coutier-Delgosha, O.; Fan, D. Combined suppression effects on hydrodynamic cavitation performance in Venturi-type reactor for process intensification. Ultrason. Sonochem. 2022, 86, 106035. [Google Scholar] [CrossRef]
  28. Ge, M.; Manikkam, P.; Ghossein, J.; Subramanian, R.; Coutier-Delgosha, O.; Zhang, G. Dynamic mode decomposition to classify cavitating flow regimes induced by thermodynamic effects. Energy 2022, 254, 124426. [Google Scholar] [CrossRef]
  29. Shi, G.; Li, H.; Liu, X.; Wang, B. Transport performance improvement of a multiphase pump for gas–liquid mixture based on the orthogonal test method. Processes 2021, 9, 1402. [Google Scholar] [CrossRef]
  30. Wang, C.; Yang, F.; Nguyen, V.; Vo, N. CFD analysis and optimum design for a centrifugal pump using an effectively artificial intelligent algorithm. Micromachines 2022, 13, 1208. [Google Scholar] [CrossRef]
  31. Nguyen, V.T.T.; Vo, T.M.N. Centrifugal Pump Design: An Optimization. Eurasia Proc. Sci. Technol. Eng. Math. 2022, 17, 136–151. [Google Scholar] [CrossRef]
  32. Zhang, W.; Zhu, B.; Wang, Z.; Wang, F. Optimization design for an impeller of the multiphase rotodynamic pump handling gas-liquid two-phase flow. Proc. Inst. Mech. Eng. Part A J. Power Energy 2022, 236, 09576509221098223. [Google Scholar] [CrossRef]
  33. Kim, J.H.; Lee, H.C.; Yoon, J.Y.; Lee, K.; Choi, Y. Multi objective optimization of a multiphase pump for offshore plants. In Proceedings of the ASME 2014 4th Joint US-European Fluids Engineering Division Summer Meeting, Chicago, IL, USA, 3–7 August 2014; American Society of Mechanical Engineers: New York, NY, USA, 2014; p. 46223. [Google Scholar]
  34. Liu, M.; Tan, L.; Cao, S. Design method of controllable blade angle and orthogonal optimization of pressure rise for a multiphase pump. Energies 2018, 11, 1048. [Google Scholar] [CrossRef] [Green Version]
  35. Peng, C.; Zhang, X.; Chen, Y.; Gong, Y.; Li, H.; Huang, S. A Method for the Integrated Optimal Design of Multiphase Pump Based on the Sparse Grid Model. Processes 2022, 10, 1317. [Google Scholar] [CrossRef]
Figure 1. Three-dimension model of multiphase pump.
Figure 1. Three-dimension model of multiphase pump.
Jmse 11 00130 g001
Figure 2. Computational mesh (a) inlet pipe; (b) outlet pipe; and (c) impeller.
Figure 2. Computational mesh (a) inlet pipe; (b) outlet pipe; and (c) impeller.
Jmse 11 00130 g002
Figure 3. Mesh independent verification.
Figure 3. Mesh independent verification.
Jmse 11 00130 g003
Figure 4. Experimental rig of the multiphase pump. (a) Schematic diagram. (b) Impeller and diffuser.
Figure 4. Experimental rig of the multiphase pump. (a) Schematic diagram. (b) Impeller and diffuser.
Jmse 11 00130 g004
Figure 5. Hydraulic characteristics between CFD and experiment.
Figure 5. Hydraulic characteristics between CFD and experiment.
Jmse 11 00130 g005
Figure 6. Cavitation characteristics.
Figure 6. Cavitation characteristics.
Jmse 11 00130 g006
Figure 7. Pressure in impeller flow passage at IGVF = 0, 0.1, 0.2. (a) Critical cavitation; (b) severe cavitation; (c) fracture cavitation.
Figure 7. Pressure in impeller flow passage at IGVF = 0, 0.1, 0.2. (a) Critical cavitation; (b) severe cavitation; (c) fracture cavitation.
Jmse 11 00130 g007
Figure 8. Pressure on 0.5 blade span at IGVF = 0, 0.1, 0.2. (a) IGVF = 0; (b) IGVF = 0.1; (c) IGVF = 0.2.
Figure 8. Pressure on 0.5 blade span at IGVF = 0, 0.1, 0.2. (a) IGVF = 0; (b) IGVF = 0.1; (c) IGVF = 0.2.
Jmse 11 00130 g008
Figure 9. Velocity in impeller under IGVF = 0, 0.1 and 0.2. (a) Critical cavitation; (b) severe cavitation; and (c) fracture cavitation.
Figure 9. Velocity in impeller under IGVF = 0, 0.1 and 0.2. (a) Critical cavitation; (b) severe cavitation; and (c) fracture cavitation.
Jmse 11 00130 g009
Figure 10. Velocity distribution at 0.5 span under IGVF = 0, 0.1, 0.2. (a) IGVF = 0; (b) IGVF = 0.1; (c) IGVF = 0.2.
Figure 10. Velocity distribution at 0.5 span under IGVF = 0, 0.1, 0.2. (a) IGVF = 0; (b) IGVF = 0.1; (c) IGVF = 0.2.
Jmse 11 00130 g010aJmse 11 00130 g010b
Figure 11. Velocity streamlines at 0.5 span under IGVF = 0, 0.1 and 0.2. (a) Critical cavitation; (b) severe cavitation; and (c) fracture cavitation.
Figure 11. Velocity streamlines at 0.5 span under IGVF = 0, 0.1 and 0.2. (a) Critical cavitation; (b) severe cavitation; and (c) fracture cavitation.
Jmse 11 00130 g011aJmse 11 00130 g011b
Table 1. Mesh quality.
Table 1. Mesh quality.
Computational DomainMax AngleMax WarpSkewAspect RatioQuality
Inlet pipe90.0377~108.2840~0.4160.80~0.990.00091~0.990.95
Impeller90.002~121.610~2.460.68~10.0018~0.990.48
Outlet pipe90.001~119.160~0.270.688~10.0006~0.990.87
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lv, W.; Wen, H.; Shi, G.; Wang, S. Analyzing the Inlet Gas Void Fraction on the Flow Characteristics for a Multiphase Pump Based on Cavitation Cases. J. Mar. Sci. Eng. 2023, 11, 130. https://doi.org/10.3390/jmse11010130

AMA Style

Lv W, Wen H, Shi G, Wang S. Analyzing the Inlet Gas Void Fraction on the Flow Characteristics for a Multiphase Pump Based on Cavitation Cases. Journal of Marine Science and Engineering. 2023; 11(1):130. https://doi.org/10.3390/jmse11010130

Chicago/Turabian Style

Lv, Wenjuan, Haigang Wen, Guangtai Shi, and Shan Wang. 2023. "Analyzing the Inlet Gas Void Fraction on the Flow Characteristics for a Multiphase Pump Based on Cavitation Cases" Journal of Marine Science and Engineering 11, no. 1: 130. https://doi.org/10.3390/jmse11010130

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop