Next Article in Journal
Experimental Research on Cavitation Evolution and Movement Characteristics of the Projectile during Vertical Launching
Previous Article in Journal
Application of Feature Point Matching Technology to Identify Images of Free-Swimming Tuna Schools in a Purse Seine Fishery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Internal Flow Characteristics of a Startup Pump Turbine at the Lowest Head under No-Load Conditions

1
Department of Energy and Power Engineering, Tsinghua University, Beijing 100084, China
2
Anhui Jinzhai Pumped Storage Power Co., Ltd., Lu’an 237333, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2021, 9(12), 1360; https://doi.org/10.3390/jmse9121360
Submission received: 15 October 2021 / Revised: 22 November 2021 / Accepted: 22 November 2021 / Published: 1 December 2021
(This article belongs to the Section Ocean Engineering)

Abstract

:
The instability of the no-load working condition of the pump turbine directly affects the grid connection of the unit, and will cause vibration and damage to the components of the unit in severe cases. In this paper, a three-dimensional full flow numerical model including the runner gap and the pressure-balance pipe was established. The method SST k-ω model was used to predict the internal flow characteristics of the pump turbine. The pressure pulsation of the runner under different operating conditions during the no-load process was compared. Because the rotation speed, flow rate, and guide vane opening of the unit change in a small range during the no-load process, the pressure pulsation characteristics of the runner are basically the same. Therefore, a working condition was selected to analyze the transient characteristics of the flow field, and it was found that there was a high-speed ring in the vaneless zone, and a stable channel vortex was generated in the runner flow passage. Analyzing the axial water thrust of each part of the runner, it was found that the axial water thrust of the runner gap was much larger than the axial water thrust of the runner blades, and it changed with time periodically. It was affected by rotor stator interaction. The main frequency was expressed as a multiple of the number of guide vanes, that is, vanes passing frequency, 22fn. During the entire no-load process, the axial water thrust of the runner changed slowly with time and fluctuated slightly.

1. Introduction

With the increasing demand of renewable energy, one of the mature technologies for large-scale energy storage is pumped-storage power stations [1,2]. The vigorous development of pumped-storage power plants can provide sufficient storage for energies generated by both the wind and the sun. As the core part of pumped-storage power station, pump turbines can achieve peak shaving and valley filling of the power grid through the conversion between pump mode and turbine mode to optimize grid load [3]. Frequent startup and shutdown and operating conditions’ conversion are the basic requirements for pump turbines [4,5]. During startup, the guide vane is gradually opened from the closed state, and the unit experiences the S characteristics area before the unit is connected to the grid [6].
A typical full characteristic curve of a pump turbine may have two S characteristic areas. One is the direction of rotation of the turbine and the other is the direction of rotation of the pump. The S characteristic of the turbine rotation direction is more obvious than the pump rotation direction, which brings a series of difficulties to the power generation mode [7,8]. In the S area, the same speed may correspond to 3 different flow rates, including the negative flow in the reverse pump mode. When the unit is operating in the S area, especially under no-load conditions, the flow is unstable [9,10]. Staubli [11] pointed out that the instability of S characteristic area is related to the slope of the unit characteristic curve under no-load conditions. Gentner [12] found that the blocking effect of the vortex flow in the vaneless zone under no-load opening is an important source of the instability of the S area. The study of vortex flow in the vaneless zone has also been reported by many scholars [13,14,15]. The misaligned guide vanes (MGV) have a controlling effect on the stability of the S area [16]. However, for actual operation, the three-dimensional flow field model not only includes the spiral case, state vane, guide vane, runner, and draft tube, but also needs to consider the influence of runner gap and pressure-balance pipe. As shown in Figure 1, the space between the top cover and the upper crown of the runner is the upper crown gap, and the space between the bottom cover and the lower band of the runner is the lower band gap. There are two pressure-balance pipes of the unit, which are located between the main shaft seal and the upper crown gap seal ring, and are connected to the draft tube through the top cover of the unit.
As mentioned above, the unstable flow in the unit under no-load conditions is related to the instability of S area. In addition, it complicates the force on the runner. Under startup no-load conditions, the axial force and radial force of the unit will change greatly [17]. Li [18,19] conducted that the high-speed vortex in the vaneless zone at a no-load condition caused a significant increase in the radial force of the runner. Meanwhile, the pump turbine is installed vertically, and the changing characteristics of the axial water thrust directly affect the safe and stable operation of the unit. Excessive axial water thrust causes accidents such as unit lifting and tile burning [20,21]. Runner gap and pressure-balance pipe are important factors that affect the axial water thrust of the unit. The instantaneous characteristics of the axial water thrust of the runner under no-load conditions have not been studied in depth. In particular, the influence of the upper crown gap and the lower band gap on the axial water thrust of the runner needs further discussion.
In this paper, a numerical simulation was conducted to analyze the internal flow characteristics of the pump turbine at the lowest head under no-load condition. Firstly, a pump turbine model including the upper crown gap, the lower band gap, and the pressure-balance pipe was established. The mesh division was performed, and the accuracy of the numerical method was verified by empirical results. Then, the pressure pulsation of the runner at different working conditions during the no-load process is compared. Finally, a working condition was selected to show that the vortex characteristics of the vaneless zone, and analyzed the instantaneous characteristics of the axial water thrust of each part of the runner.

2. Methods

The Navier–Stokes equations of fluid momentum and the flow continuity equation were used to analyze the complicated flow field in the pump turbine.
ρ t + x j ρ u j = 0
t ρ u i + x j ρ u i u j = p x i + x j u u i t ρ u i u j ¯ + S M
where ρ is water density, p is the pressure, u is the flow velocity, μ is the dynamic viscosity, S M is the external momentum source term, and ρ u i u j ¯ is the Reynolds stresses. The relationship between the Reynolds stresses and the mean velocity gradients is established by the eddy viscosity assumption.
ρ u i u j ¯ = μ t u i x j + u j x i 2 3 ρ k δ i j
where δ i j is the Kronecker delta (when i = j , δ i j = 1; when i j , δ i j = 0), and μ t is the eddy viscosity and can be solved by the turbulent kinetic energy, k and the turbulent frequency, ω .
μ t = ρ k ω
The starting point of the present formulations is the k-ω model developed by Wilcox as follows:
ρ k t + x j ρ u j k = x j μ + μ t σ k 1 k x j + P k β ρ k ω
ρ ω t + x j ρ u j ω = x j μ + μ t σ ω 1 ω x j + α 1 ω k P k β 1 ρ ω 2
Since the Wilcox model is very sensitive to freestream conditions and the inlet value has a significant influence on the calculation results. Menter [22] developed the SST k-ω model (Shear Stress Transport k-ω), which used the k-ω model near the surface and the k-ε model in the outside region. The k-ω model was thereby multiplied by a blending function F1 and the k-ε model by a function (1 − F1).
ρ k t + x j ρ u j k = x j μ + μ t σ k 3 k x j + P k β ρ k ω
ρ ω t + x j ρ u j ω = x j μ + μ t σ ω 3 ω x j + 1 F 1 2 ρ 1 σ ω 2 ω k x j ω x j + α 3 ω k P k β 3 ρ ω 2
The functions in the model can be solved by the linear relationship in Equation (7).
Φ 3 = F 1 Φ 1 + 1 F 1 Φ 2
where β = 0.09 , α 1 = 5 / 9 , β 1 = 0.075 , σ k 1 = 2 , σ ω 1 = 2 , α 2 = 0.44 , β 2 = 0.0828 , σ k 2 = 1 , σ ω 2 = 1 / 0.856 .
The k-ω based SST model accounts for the transport of the turbulent shear stress by limiting eddy viscosity shown in Equations (8)–(10). It gives highly accurate predictions of the onset and the amount of flow separation under adverse pressure gradients [23,24]. Therefore, the SST k-ω is used in present simulation.
ν t = μ t / ρ
ν t = a 1 k max a 1 ω , S F 2
F 2 = tanh a r g 2 2
a r g 2 = max 2 k β ω y , 500 ν y 2 ω
where a 1 = 5 / 9 , S is an invariant measure of the strain rate, y is the distance to the nearest wall, and ν is the kinematic viscosity.
The three-dimensional model of the pump turbine of a power station is shown in Figure 1, including the spiral case, stay vane, guide vane, runner, upper crown seal gap, lower band seal gap, pressure-balance pipe, and draft tube. The pump turbine studied in this paper featured 22 stay vanes, 22 guide vanes, and 13 runner blades.
The simulation was carried out in ANSYS CFX with the inlet set to total pressure and the outlet set to static pressure. The other walls were set as no slip walls. The runner and the gap between runners were set as a rotating domain, and the other regions are static domains. The transient motion interfaces of the rotating domain and the static domain were set to transient rotor stator. The turbulence model selected the SST k-ω model. The time step was 5 × 10 4 . In the solution of the incompressible flow field by RANS (Reynolds-averaged Navier-Stokes equations), SIMPLEC (SIMPLE Consistent) was used to realize the coupling of the velocity field and the pressure field, and the convergence residual was 1 × 10 4 .
The pump turbine adopted a combination of structured mesh and unstructured mesh. The structured hexahedral mesh was used for the runner and the seal gap. Moreover, mesh refinement was performed on the blade surface and the gap labyrinth ring. The mesh result was shown in Figure 2.
The efficiency of the unit under different water heads was shown in Figure 3. The numerical simulation results of the pump turbine basically agree with the curve of the test results. Compared with the test results, this paper considered the energy loss of the runner gap, so the simulation result was slightly smaller.

3. Results and Discussion

The relative rotation speed, relative flow rate, and relative guide vane opening during the startup process when the unit at the lowest water head are shown in Figure 4. From t = 0 s to t = 18 s, when the opening of the guide vane increases the flow rate increases, resulted in the increase of the energy acting on the runner, then the speed rose rapidly, approaching the rated speed. The guide vane slightly closed up to prevent the speed from rising too fast and exceeding the rated speed. From t = 18 s to t = 28 s, the opening of the guide vane decreased and the flow rate decreased. At the same time, the water was still working on the runner, the speed of the runner was also increasing. After the unit was connected to the power grid, it entered the no-load operation area.
Figure 5 shows the changes of the flow rate, the rotation speed, and the guide vane opening when the relative torque was 0. The relative flow rate varied from 0.135 to 0,128, the speed varied from 1.0 to 0.996, and the guide vane opening varied from 0.175 to 0.16. The pressure pulsation under three working conditions A, B, and C were selected for comparison. The operating parameters of the three working conditions are shown in Table 1.
Figure 6 shows the distribution of monitoring points in the runner, where rn1-rn3 are located on the suction surface of the blade, rn4-rn6 are located on the pressure surface of the blade, and rn7 in the runner’s flow channel.
As shown in Figure 7, rn1 is located at the inlet of the blade suction surface. Under the influence of the rotor–stator interaction by the 22 guide vanes, the main frequency of rn1 was 22fn (fn is the rotational frequency of the runner). rn2 is located in the middle of the suction surface of the blade, away from the inlet side, and the influence of the rotor-stator interaction was weakened. The pressure pulsation amplitude of 22fn decreased, and the pressure pulsation amplitude of 1.0fn increased. rn3 was located at the outlet of the suction surface of the blade, and the dominant frequency was 1.0fn. The dominant frequency at the pressure surfaces of rn4, rn5, and rn6 were the same as that of rn1, rn2, and rn3. The monitoring points nearby the inlet of the blade was affected by the rotor–stator interaction. rn7 is located in the middle of the flow channel, and still had the frequency of the rotor–stator interaction of 22fn.
As shown in Equation (14), the relative axial water thrust coefficient, Fy*, is defined. The calculation of Fy is based on the results of CFD (Computational Fluid Dynamics) calculations and extracts the axial water thrust of the runner.
F y * = F y m g
where m is the weight of the rotating parts of the pump turbine unit. g is the acceleration due to gravity, 9.8 kg/N.
Figure 8 is an exploded view of the axial water thrust of the runner, including the axial water thrust of the outer surface of the upper crown of the runner (UCO-Fy), the axial water thrust of the inner surface of the upper crown of the runner (UCI-Fy), the total axial water thrust of the upper crown of the runner(UC-Fy), the axial water thrust of the outer surface of the lower band of the runner (LBO-Fy), the axial water thrust of the inner surface of the lower band of the runner (LBI-Fy), the total axial water thrust of the lower band of the runner(LB-Fy), the axial water thrust of the pressure and suction surfaces of the runner blades (BFy), and the axial water thrust of the runner (RN-Fy). Figure 9 shows the axial water thrust of the runner under conditions A, B, and C. The results show that the axial water thrust of the runner was smaller than the gravity of the runner. In the no-load process, it changed steadily with time, showing a trend of small increase and gradual decrease.
During the no-load process, the values of the flow rate, the rotation speed, and the guide vane opening of the unit are in small ranges. Therefore, under the different working conditions, A, B and C, the main frequency of pressure pulsation at each monitoring point in the runner was basically the same. The amplitude of the pressure pulsation had some differences, but the pulsation magnitude was the same. Moreover, the axial water thrust of the runner did not change much during the no-load process. Therefore, the unsteady calculation of a certain operating point can be selected to reflect the instantaneous flow characteristics of the unit in the no-load process. The following is the unsteady flow field analysis of the pump turbine under A working condition.
Figure 10 shows the velocity vector and streamline distribution of the guide vane-runner area. It can be seen from the figure that in the no-load process, the guide vane had a small opening, a high rotation speed, and a low flow rate. Therefore, the speed difference between the pressure surface and the suction surface of the guide vane was relatively high, and a high-speed ring was formed in the vaneless area. A channel vortex was generated near the inlet of the runner, and the shape of the vortex was stable in one period.
The axial water thrust of each part shown in Figure 8 was analyzed, and the contribution of the gap flow to the axial water thrust of the runner was clarified. In Figure 11, the axial water thrust of the blade (BFy) is small, about 1/10 of the gravity of the rotating part of the unit, and its direction is the opposite of the gravity (+Y). Moreover, the BFy changes periodically, and its main frequency is the rotor–stator interaction frequency, 22fn.
Figure 12 shows the axial water thrust on the inner surface of the upper crown of the runner (UCI-Fy). The results show that the UCI-Fy was very large, with a small fluctuation range over time, which was basically about 6.5 times the gravity of the rotating part of the unit (6.5 mg), and the direction was at the +Y axis. It changed with time periodically and was affected by rotor–stator interaction. The dominant frequency was 22fn. Figure 13 shows the axial water thrust on the outer surface of the upper crown of the runner (UCO-Fy). The UCO-Fy periodically changes along the -Y axis, and the absolute value was greater than that of the inner surface of the upper crown, about 6.7 mg. The main frequency was 22fn. The total axial water thrust on the upper crown of the runner was downward (-Y), about 0.24 mg, and its vibration is shown in Figure 14.
Figure 15 shows the axial water thrust on the inner surface of the lower band of the runner (LBI-Fy). The LBI-Fy changed along the -Y axis, about 5.6 mg. The main frequency of the periodic change was 22fn. In Figure 16, it can be seen that the axial water thrust on the outer surface of the lower band of the runner (LBO-Fy) is about 6.45 mg at the positive direction of the Y axis, and greater than LBI-Fy. Moreover, LBO-Fy did not change periodically, but decreased slightly over time.
In Figure 17, the total axial water thrust of the lower band of the runner is along the positive direction of the Y axis, about 0.82 mg. Combining the axial water thrust of the blades in Figure 11 and the axial water thrust of the upper crown in Figure 14, the total water thrust of the runner (RN-Fy) in Figure 18 is about 0.6 mg. Affected by LB-Fy, it decreased slightly over time.

4. Conclusions

In this paper, the CFD method considering the gap flow and pressure-balance pipe flow of the runner was applied to analyze the internal flow of the pump turbine at the lowest head under no-load condition.
The SST k-ω model was used to predict the three-dimensional internal flow characteristics of pump turbine. The numerical simulation of the efficiency of the unit under different heads was basically consistent with the test results. Considering the energy loss of gap flow, the simulation results of the pump turbine were slightly smaller.
Three working conditions (A, B, C) of no-load process were selected for analysis and comparison. The flow rate, rotation speed, and guide vane opening of the unit under no-load conditions all change within a small range. It can be seen that the main frequency of the monitoring points in the runner under each working condition was basically the same, and the amplitude of pressure pulsation was little different, which was mainly affected by rotor–stator interaction. Among them, the inlet of the runner was close to the guide vane, which was mainly affected by the 22 guide vanes, and there was an obvious 22fn pressure pulsation. The outlet of the runner was far away from the guide vane, which was represented by the rotation frequency of the runner, 1.0fn.
A working condition was selected for the instantaneous internal flow analysis. The speed difference between the pressure surface and the suction surface of the guide vane was relatively large, forming a high-speed ring in the vaneless area. A stable channel vortex was generated in the flow passage near the inlet of the runner. Compared with the axial water thrust of the runner blades BFy, the axial water thrust of the upper crown of the runner (UCI-Fy, UCO-Fy) and the axial water thrust of the lower band of the runner (LBI-Fy, LBO-Fy) were very large, about 560 times of BFy. The BFy, UCI-Fy, UCO-Fy, and LBI-Fy changed periodically. Affected by rotor–stator interaction, it was 22fn. The total axial water thrust of the runner RN-Fy was in the opposite direction of gravity, about 0.6 mg. During the entire no-load process, the axial water thrust of the runner changed slowly with time.

Author Contributions

Conceptualization, W.W.; methodology, X.W.; software, M.N.; investigation, W.W.; writing—original draft preparation, W.W.; writing—review and editing, Z.W.; supervision, Z.W.; project administration, C.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by National Natural Science Foundation of China (No.: 51876099).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Thanks for the test data provided by relevant units and the support of National Natural Science Foundation of China (No.: 51876099).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Colmenar-Santos, A.; Linares-Mena, A.R.; Fernandez Velazquez, J.; Borge-Diez, D. Energy-efficient three-phase bidirectional converter for grid-connected storage applications. Energy Convers. Manag. 2016, 127, 599–611. [Google Scholar] [CrossRef]
  2. Mei, Z. Pumped Storage Power Generation Technology; Machinery Industry Press: Beijing, China, 2000. [Google Scholar]
  3. Tagare, D.M. Electricity Power Generation: The Changing Dimensions; Wiley-IEEE Press: Hoboken, NJ, USA, 2011. [Google Scholar]
  4. Zhang, Y.; Zhang, Y.; Wu, Y. A review of rotating stall in reversible pump turbine. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2017, 231, 1181–1204. [Google Scholar] [CrossRef] [Green Version]
  5. Zhai, L.M.; Wang, Z.W.; Luo, Y.Y.; Li, Z. TEHD analysis of a bidirectional thrust bearing in a pumped storage unit. J. Ind. Lubr. Tribol. 2016, 68, 315–324. [Google Scholar] [CrossRef]
  6. Xia, L.S.; Cheng, Y.G.; You, J.F.; Zhang, X.X.; Yang, J.D.; Qian, Z.D. Mechanism of the S-Shaped Characteristics and the Runaway Instability of Pump turbines. J. Fluids Eng. 2017, 139, 1–14. [Google Scholar] [CrossRef]
  7. Meier, L. Hydro pump Storage Machines: The importance of the product today for Power Control and other Ancillary Services. In Proceedings of the HydroVision International 2010: A Global Gathering, Charlotte, NC, USA, 27–30 July 2010. [Google Scholar]
  8. Demin, L.; Xiaobing, L. Numerical optimization on the “S” characteristics of model pump turbine. In Proceedings of the 16th International Symposium on Transport Phenomena and Dynamics of Rotating Machinery, Honolulu, HW, USA, 10–15 April 2016. hal-01894412. [Google Scholar]
  9. Zuo, Z.; Liu, S. Flow-induced instabilities in pump turbines in China. Engineering 2017, 3, 504–511. [Google Scholar] [CrossRef]
  10. Cavazzini, G.; Covi, A.; Pavesi, G.; Ardizzon, G. Analysis of the unstable behavior of a pump turbine in turbine mode: Fluid-dynamical and spectral characterization of the S-shape characteristic. J. Fluids Eng. 2016, 138, 021105. [Google Scholar] [CrossRef]
  11. Staubli, T.; Widmer, C.; Tresch, T.; Sallaberger, M. Starting Pump-Turbines with Unstable Characteristics. In Proceedings of the Hydro 2010, Lisbon, Portugal, 27–29 September 2010. [Google Scholar]
  12. Gentner, C.; Sallaberger, M.; Widmer, C.; Bobach, B.J.; Jaberg, H.; Schiffer, J.; Guggenberger, M. Comprehensive experimental and numerical analysis of instability phenomena in pump turbines. IOP Conf. Ser. Earth Environ. Sci. 2014, 22, 032046. [Google Scholar] [CrossRef] [Green Version]
  13. Widmer, C.; Staubli, T.; Ledergerber, N. Unstable Characteristics and Rotating Stall in Turbine Brake Operation of Pump turbine. J. Fluids Eng. 2011, 133, 041101. [Google Scholar] [CrossRef]
  14. Hosseinimanesh, H.; Devals, C.; Nennemann, B.; Reggio, M.; Guibault, F. A numerical study of Francis turbine operation at no-load condition. J. Fluids Eng. 2017, 139, 011104. [Google Scholar] [CrossRef]
  15. Fu, X.; Li, D.; Wang, H.; Zhang, G.; Li, Z.; Wei, X. Analysis of transient flow in a pump turbine during the load rejection process. J. Mech. Sci. Technol. 2018, 32, 2069–2078. [Google Scholar] [CrossRef]
  16. Sun, H.; Xiao, R.; Liu, W.; Wang, F. Analysis of S characteristics and pressure pulsations in a pump turbine with misaligned guide vanes. J. Fluids Eng. 2013, 135, 511011. [Google Scholar] [CrossRef]
  17. Zhang, Y.; Chen, T.; Li, J.; Yu, J. Experimental study of load variations on pressure pulsations in a prototype reversible pump turbine in generating mode. J. Fluids Eng. 2017, 139, 074501. [Google Scholar] [CrossRef]
  18. Li, Z.; Bi, H.; Wang, Z.; Yao, Z. Three-dimensional simulation of unsteady flows in a pump turbine during start-up transient up to speed no-load condition in generating mode. Proc. Inst. Mech. Eng. Part A J. Power Energy 2016, 230, 570–585. [Google Scholar] [CrossRef]
  19. Li, Z.; Bi, H.; Karney, B.; Wang, Z.; Yao, Z. Three-dimensional transient simulation of a pump turbine during normal turbine shutdown. J. Hydraul. Res. 2017, 55, 1–18. [Google Scholar] [CrossRef]
  20. Mao, Z.; Tao, R.; Chen, F.; Wang, Z. Investigation of the Starting-Up Axial Hydraulic Force and Structure Characteristics of Pump Turbine in Pump Mode. J. Mar. Sci. Eng. 2021, 9, 158. [Google Scholar] [CrossRef]
  21. Sun, H.; Lv, Y.; Ni, J.; Jiang, X.; Wang, Z. Effect of seal locations of pump-turbine on axial hydraulic trust. J. Mar. Sci. Eng. 2021, 9, 623. [Google Scholar] [CrossRef]
  22. Menter, F.R. Two-equation eddy-viscosity turbulence models for engineering applications. AIAA J. 1994, 32, 1598–1605. [Google Scholar] [CrossRef] [Green Version]
  23. Cademartori, C.; Cravero, C.; Marini, M.; Marsano, D. CFD Simulation of the Slot Jet Impingement Heat Transfer Process and Application to a Temperature Control System for Galvanizing Line of Metal Band. J. Appl. Sci. 2021, 11, 1149. [Google Scholar] [CrossRef]
  24. Huang, H.; Sun, T.; Zhang, G.; Li, D.; Wei, H. Evaluation of a developed SST k-ω turbulence model for the prediction of turbulent slot jet impingement heat transfer. Int. J. Heat Mass Transf. 2019, 139, 700–712. [Google Scholar] [CrossRef]
Figure 1. Three-dimensional model of pump turbine.
Figure 1. Three-dimensional model of pump turbine.
Jmse 09 01360 g001
Figure 2. Parts of the mesh in the runner and the upper crown seal gap. (a) Refinement of the blade surface mesh (b) Refinement of the gap labyrinth ring mesh.
Figure 2. Parts of the mesh in the runner and the upper crown seal gap. (a) Refinement of the blade surface mesh (b) Refinement of the gap labyrinth ring mesh.
Jmse 09 01360 g002
Figure 3. Comparison of numerical simulation and experimental results.
Figure 3. Comparison of numerical simulation and experimental results.
Jmse 09 01360 g003
Figure 4. The relative flow rate, relative speed, relative guide vane opening, and relative torque changes of the unit under the lowest head start-up.
Figure 4. The relative flow rate, relative speed, relative guide vane opening, and relative torque changes of the unit under the lowest head start-up.
Jmse 09 01360 g004
Figure 5. The change of flow rate, rotation speed, and guide vane opening during the no-load process of the lowest head. (a) Relative flow rate (b) Relative rotation speed (c) Relative guide vane opening.
Figure 5. The change of flow rate, rotation speed, and guide vane opening during the no-load process of the lowest head. (a) Relative flow rate (b) Relative rotation speed (c) Relative guide vane opening.
Jmse 09 01360 g005
Figure 6. Runner monitoring point distribution map.
Figure 6. Runner monitoring point distribution map.
Jmse 09 01360 g006
Figure 7. The pressure pulsation of the monitoring point in the runner under different working conditions. (a) rn1 (b) rn2 (c) rn3 (d) rn4 (e) rn5 (f) rn6 (g) rn7.
Figure 7. The pressure pulsation of the monitoring point in the runner under different working conditions. (a) rn1 (b) rn2 (c) rn3 (d) rn4 (e) rn5 (f) rn6 (g) rn7.
Jmse 09 01360 g007aJmse 09 01360 g007b
Figure 8. The exploded view of the axial water thrust of the runner.
Figure 8. The exploded view of the axial water thrust of the runner.
Jmse 09 01360 g008
Figure 9. The axial water thrust change of the runner during no-load process.
Figure 9. The axial water thrust change of the runner during no-load process.
Jmse 09 01360 g009
Figure 10. The internal flow characteristics of guide vanes and runners in one period. (a) 0 T (b) 0.2 T (c) 0.4 T (d) 0.6 T (e) 0.8 T (f) 1.0 T.
Figure 10. The internal flow characteristics of guide vanes and runners in one period. (a) 0 T (b) 0.2 T (c) 0.4 T (d) 0.6 T (e) 0.8 T (f) 1.0 T.
Jmse 09 01360 g010aJmse 09 01360 g010b
Figure 11. The axial water thrust change of the runner blade. (a) Time domain diagram of BFy. (b) Frequency domain diagram of BFy.
Figure 11. The axial water thrust change of the runner blade. (a) Time domain diagram of BFy. (b) Frequency domain diagram of BFy.
Jmse 09 01360 g011
Figure 12. The axial water thrust of the runner upper crown inner surface. (a) Time domain diagram of UCI-Fy. (b) Frequency domain diagram of UCI-Fy.
Figure 12. The axial water thrust of the runner upper crown inner surface. (a) Time domain diagram of UCI-Fy. (b) Frequency domain diagram of UCI-Fy.
Jmse 09 01360 g012
Figure 13. The axial water thrust of the runner upper crown outer surface. (a) Time domain diagram of UCO-Fy. (b) Frequency domain diagram of UCO-Fy.
Figure 13. The axial water thrust of the runner upper crown outer surface. (a) Time domain diagram of UCO-Fy. (b) Frequency domain diagram of UCO-Fy.
Jmse 09 01360 g013
Figure 14. The axial water thrust of the runner upper crown.
Figure 14. The axial water thrust of the runner upper crown.
Jmse 09 01360 g014
Figure 15. The axial water thrust of the runner lower band inner surface. (a) Time domain diagram of LBI-Fy. (b) Frequency domain diagram of LBI-Fy.
Figure 15. The axial water thrust of the runner lower band inner surface. (a) Time domain diagram of LBI-Fy. (b) Frequency domain diagram of LBI-Fy.
Jmse 09 01360 g015
Figure 16. The axial water thrust of the runner lower band outer surface.
Figure 16. The axial water thrust of the runner lower band outer surface.
Jmse 09 01360 g016
Figure 17. The axial water thrust of the runner lower band.
Figure 17. The axial water thrust of the runner lower band.
Jmse 09 01360 g017
Figure 18. The total axial water thrust of the runner.
Figure 18. The total axial water thrust of the runner.
Jmse 09 01360 g018
Table 1. The operating parameters of the three working conditions.
Table 1. The operating parameters of the three working conditions.
ConditionsABC
Parameters
Relative flow rate0.1340.130.128
Relative rotation speed0.9990.9960.995
Relative guide vane opening0.1750.1650.16
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, W.; Wang, X.; Wang, Z.; Ni, M.; Yang, C. Analysis of Internal Flow Characteristics of a Startup Pump Turbine at the Lowest Head under No-Load Conditions. J. Mar. Sci. Eng. 2021, 9, 1360. https://doi.org/10.3390/jmse9121360

AMA Style

Wang W, Wang X, Wang Z, Ni M, Yang C. Analysis of Internal Flow Characteristics of a Startup Pump Turbine at the Lowest Head under No-Load Conditions. Journal of Marine Science and Engineering. 2021; 9(12):1360. https://doi.org/10.3390/jmse9121360

Chicago/Turabian Style

Wang, Wei, Xi Wang, Zhengwei Wang, Mabing Ni, and Chunan Yang. 2021. "Analysis of Internal Flow Characteristics of a Startup Pump Turbine at the Lowest Head under No-Load Conditions" Journal of Marine Science and Engineering 9, no. 12: 1360. https://doi.org/10.3390/jmse9121360

APA Style

Wang, W., Wang, X., Wang, Z., Ni, M., & Yang, C. (2021). Analysis of Internal Flow Characteristics of a Startup Pump Turbine at the Lowest Head under No-Load Conditions. Journal of Marine Science and Engineering, 9(12), 1360. https://doi.org/10.3390/jmse9121360

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