Pressure Fluctuations in the S-Shaped Region of a Reversible Pump-Turbine

Numerical simulations were performed to investigate pressure fluctuations in the S-shaped region of a pump-turbine model. Analyses focused on pressure fluctuations in the draft tube and in the gap between the guide vanes and runner. Calculations were made under six different operating conditions with a constant guide vane opening, and the best efficiency point, runaway point, and low-discharge point in the turbine brake zone were determined. The simulated results were compared with experimental measurements. In the draft tube, a twin vortex rope was observed. In the gap between the guide vanes and runner, a low frequency component was captured at both the runaway and low-discharge points in the turbine brake zone, which rotated at 65% of the runner frequency. This low frequency component was induced by the rotating stall phenomenon. At the runaway point, a single stall cell was found in the gap between the guide vanes and runner, while at the low-discharge point, four stall cells were observed.


Introduction
Pumped hydro energy storage (PHES) is currently the only proven large-scale (>100 MW) energy storage technology.There are great benefits in using PHES in electricity generation systems.Its flexibility can provide upregulation and downregulation.Its quick start capabilities make it suitable for black starts, and for providing spinning and standing reserves.Interest in this technology has recently been renewed because of the increasing use of renewable energy such as wind-powered electricity generation.Such technology is weather dependent, and so is highly variable over time [1,2].
A key component of PHES stations is the pump-turbine.Pump-turbines usually use a single runner to perform the functions of both pump and turbine.In order to adapt to the load changes of the power system, a pump-turbine is required to frequently switch between turbine mode and pump mode.This leads to extended operations in off-design operating conditions, such as at start-up and during load rejection processes.In the associated transient processes of pump-turbines, the S-shaped characteristics of the performance curves in turbine mode usually induce instability problems.The most common problems are difficulties in synchronizing with power grids during turbine start-ups, and unstable performance during turbine load rejections [3,4].
In recent years, both experimental and modelling studies have investigated the unstable flow characteristics of pump-turbines operating in off-design conditions [5][6][7][8][9][10]. Husmatuchi et al. [5,6] experimentally investigated pressure fluctuations in the gap between the guide vanes and runner.In their research, a low frequency component was observed under runaway conditions and turbine brake conditions.This low frequency instability was induced by one stall cell rotating inside the runner.Widmer et al. [7] performed numerical simulations of flows in a model pump-turbine.It was concluded Energies 2017, 10, 96 2 of 13 that the unstable characteristics observed with low flow masses were induced by stationary vortex formation and rotating stall.Cavazzini et al. [8] carried out a numerical analysis of the load rejection process at a constant rate, and with a large guide vane opening of a pump-turbine in turbine mode.The flow analysis showed clearly the onset and development of unsteady phenomena progressively developing into an organized rotating stall during turbine brake operations.Both experimental and numerical studies have revealed complex flow characters in the S-shaped region of pump-turbines.Such vortexes may partially block the flow [7,9,10], or cause stall cells to rotate in the runner channels at 50%-70% of the runner rotation frequency [5,6,8].
As pump-turbines operate in or near the S-shaped region, vortex ropes are inevitably formed in the draft tubes [11].Ruprecht et al. [12] and Kirschner et al. [13] investigated unstable characteristics in the draft tube of a pump-turbine.They reported that pressure fluctuations are mainly induced by vortex ropes in the draft tube.Usually, a single-helical vortex rope is formed in the draft tube; therefore, most research has focused on regimes with a single helical vortex rope [11][12][13].However, for some hydraulic turbines, at small flow masses, the single helix is replaced by a double helix, and twin vortex ropes can be formed [11,14].One of the first detailed reports of twin vortex ropes in draft tubes was presented by Wahl [14], who was able to measure their precession frequency.
In this paper, numerical simulations are performed to study instability near, or in, the S-shaped region of a pump-turbine model operating under off-design conditions.The model runner was designed by Wang et al. [15][16][17] during the development of the Liyang Pumped Storage Power Station in Jiangsu Province, China.Model tests were conducted in a stand high-head test rig at the Harbin Institute of Large Electric Machinery of China.The analyses mainly focus on the pressure fluctuations and flow characteristics in the draft tube and in the gap between the guide vanes and runner.In the draft tube of the pump-turbine model, pressure fluctuations are mainly caused by a twin vortex rope.In the gap between the guide vanes and runner, a low frequency component induced by rotating stalls is captured both at runaway and low-discharge conditions.

Scaled Pump-Turbine Model
The scaled model pump-turbine runner used in this study is shown in Figure 1.The model was designed in order to conduct measurements on a standard test rig.The main parameters of the model runner are listed in Table 1.The scale of the model was one-tenth of the prototype turbines installed in aforementioned Liyang Pumped Storage Power Station.The rated specific speed of both model and prototype is n s = 145.80m • Kw.The specific parameters of the prototype pump-turbine are as follows.In turbine mode, the rated head is H r = 259 m and the rated output power is P r = 255 MW.The rotational speed is n r = 300 rpm for both turbine and pump modes.In pump mode, the maximum and minimum heads are H max = 298 m and H min = 239 m, respectively.
Energies 2017, 10, 96 2 of 13 was concluded that the unstable characteristics observed with low flow masses were induced by stationary vortex formation and rotating stall.Cavazzini et al. [8] carried out a numerical analysis of the load rejection process at a constant rate, and with a large guide vane opening of a pump-turbine in turbine mode.The flow analysis showed clearly the onset and development of unsteady phenomena progressively developing into an organized rotating stall during turbine brake operations.Both experimental and numerical studies have revealed complex flow characters in the S-shaped region of pump-turbines.Such vortexes may partially block the flow [7,9,10], or cause stall cells to rotate in the runner channels at 50%-70% of the runner rotation frequency [5,6,8].
As pump-turbines operate in or near the S-shaped region, vortex ropes are inevitably formed in the draft tubes [11].Ruprecht et al. [12] and Kirschner et al. [13] investigated unstable characteristics in the draft tube of a pump-turbine.They reported that pressure fluctuations are mainly induced by vortex ropes in the draft tube.Usually, a single-helical vortex rope is formed in the draft tube; therefore, most research has focused on regimes with a single helical vortex rope [11][12][13].However, for some hydraulic turbines, at small flow masses, the single helix is replaced by a double helix, and twin vortex ropes can be formed [11,14].One of the first detailed reports of twin vortex ropes in draft tubes was presented by Wahl [14], who was able to measure their precession frequency.
In this paper, numerical simulations are performed to study instability near, or in, the S-shaped region of a pump-turbine model operating under off-design conditions.The model runner was designed by Wang et al. [15][16][17] during the development of the Liyang Pumped Storage Power Station in Jiangsu Province, China.Model tests were conducted in a stand high-head test rig at the Harbin Institute of Large Electric Machinery of China.The analyses mainly focus on the pressure fluctuations and flow characteristics in the draft tube and in the gap between the guide vanes and runner.In the draft tube of the pump-turbine model, pressure fluctuations are mainly caused by a twin vortex rope.In the gap between the guide vanes and runner, a low frequency component induced by rotating stalls is captured both at runaway and low-discharge conditions.

Scaled Pump-Turbine Model
The scaled model pump-turbine runner used in this study is shown in Figure 1.The model was designed in order to conduct measurements on a standard test rig.The main parameters of the model runner are listed in Table 1.The scale of the model was one-tenth of the prototype turbines installed in aforementioned Liyang Pumped Storage Power Station.The rated specific speed of both model and prototype is    Model tests were conducted on a standard hydraulic machinery test rig at the Harbin Institute of Large Electric Machinery, China.Figure 2 gives a schematic diagram of the model test rig.The test rig had two test stations, A and B, of which Station A could be used to conduct the model tests for Francis turbines and reversible pump-turbines.The main specifications of the test rig and the pressure sensor (PCB112A22) used for the pressure fluctuation measurement are listed in Table 2.All measurements were conducted in accordance with International Electrotechnical Commission (IEC) Standard 60193 [18].
Energies 2017, 10, 96 3 of 13 Model tests were conducted on a standard hydraulic machinery test rig at the Harbin Institute of Large Electric Machinery, China.Figure 2 gives a schematic diagram of the model test rig.The test rig had two test stations, A and B, of which Station A could be used to conduct the model tests for Francis turbines and reversible pump-turbines.The main specifications of the test rig and the pressure sensor (PCB112A22) used for the pressure fluctuation measurement are listed in Table 2.All measurements were conducted in accordance with International Electrotechnical Commission (IEC) Standard 60193 [18].In Figure 3, the measured S-shaped characteristics in some smaller guide vane openings are given.Of these, the unit speed n 11 , unit discharge Q 11 , and unit moment M 11 are defined as follows: where n stands for the runner rotating speed, D 2 for the runner inlet diameter in turbine mode, H for the head in the model tests, Q for the discharge, and M for the moment.The four-quadrant test was conducted under a constant speed of n = 500 rpm.In Figure 3, only the first quadrant is illustrated, in order to clearly show the S-shaped characteristics.The normal operating range of the turbine mode is indicated by the two vertical lines.
Energies 2017, 10, 96 4 of 13 In Figure 3, the measured S-shaped characteristics in some smaller guide vane openings are given.Of these, the unit speed 11 n , unit discharge 11 Q , and unit moment 11 M are defined as follows: where n stands for the runner rotating speed, D2 for the runner inlet diameter in turbine mode, H for the head in the model tests, Q for the discharge, and M for the moment.The four-quadrant test was conducted under a constant speed of 500 rpm n = . In Figure 3, only the first quadrant is illustrated, in order to clearly show the S-shaped characteristics.The normal operating range of the turbine mode is indicated by the two vertical lines.

Numerical Solutions
Three-dimensional turbulent flow simulations were conducted for the full-passage pump-turbine.The extensively-used commercial code ANSYS CFX 15.0 (ANSYS Inc., Canonsburg, PA, USA) was used to conduct the numerical simulations.
Turbulence models are important factors in computational fluid dynamics (CFD).An advanced turbulence model, the detached eddy simulation (DES), was adopted to provide highly detailed simulations of flow patterns and unsteady phenomena.The shear stress transport ω k − model was used in the boundary layer, while the Smagorinsky-Lilly model was applied in detached regions [8].
The ANSYS ICEM (a powerful meshing software) and Turbo Grid were used for grid generation.The computational domain was divided into five parts for meshing (Figure 4), namely

Numerical Solutions
Three-dimensional turbulent flow simulations were conducted for the full-passage pump-turbine.The extensively-used commercial code ANSYS CFX 15.0 (ANSYS Inc., Canonsburg, PA, USA) was used to conduct the numerical simulations.
Turbulence models are important factors in computational fluid dynamics (CFD).An advanced turbulence model, the detached eddy simulation (DES), was adopted to provide highly detailed simulations of flow patterns and unsteady phenomena.The shear stress transport k − ω model was used in the boundary layer, while the Smagorinsky-Lilly model was applied in detached regions [8].
The ANSYS ICEM (a powerful meshing software) and Turbo Grid were used for grid generation.The computational domain was divided into five parts for meshing (Figure 4), namely the spiral casing, stay vanes, guide vanes, runner and draft tube.Structured mesh was created for all these parts except for the spiral casing tongue, where unstructured mesh was employed due to the irregular structure.
Energies 2017, 10, 96 5 of 13 the spiral casing, stay vanes, guide vanes, runner and draft tube.Structured mesh was created for all these parts except for the spiral casing tongue, where unstructured mesh was employed due to the irregular structure.At the inlet of the spiral casing, the mass flow rate was specified according to stochastic fluctuations of the velocities with a 5% free stream turbulent intensity.Static pressure was set at the outlet of the draft tube.The interfaces between stator-rotor blocks were set as the standard transient sliding interfaces.Walls were defined with no-slip wall boundary conditions.Runner blades, runner hub and shroud were fixed as rotating walls.The time-step size was set as 0.000139 s, such that one runner revolution was divided into 360 time steps.Second-order implicit time-stepping was adopted for time discretization.The convergence residuals for continuity and momentum equations were below 1.0 × 10 −5 .All calculations were conducted in a cluster computer with eight Intel 5645 2.4 GHz processors, 96 GB RAM, and 2 TB hard drives (Dell Inc., Round Rock, TX, USA).
As shown in Figure 3, at a guide vane opening of 17.0 mm A = , OP1 is the best efficient point (BEP), OP2-OP4 are the operating points in the normal operating range, while OP5 is the runaway point and OP6 is the low-discharge point in the turbine brake.Usually, the guide vane opening is given as an angle; however, its width was used in this study.The opening angle was o γ 18.0 = at A = 17.0 mm.Four meshes (Table 3) were used to test mesh independence.The calculated points were OP2-OP4 (Figure 4).The calculated efficiency errors for OP2-OP4 are given in Table 4. Here, the error was calculated as ( ) In the simulations, the head cover and stay ring were not included; therefore, some leakage losses and the mechanical losses were not incorporated.The calculated efficiency cal η is mainly the hydraulic efficiency, and is calculated by means of a time average after the calculated unsteady condition becomes stable.Simulated results show that mesh density has a weak influence on efficiency as the grid density is high (3.62 million).After considering the complex unsteadiness, Mesh III was chosen for the simulations.The validation of the numerical results in the next section also shows that the numerical solutions used can provide reliable results.Table 5 gives more information on Mesh III.There were a total of 4.5 million elements in the whole domain.Figure 3 shows the boundary mesh for different components.Ten element layers with normal stretching ratios were created from the walls for improving the velocity profile from the At the inlet of the spiral casing, the mass flow rate was specified according to stochastic fluctuations of the velocities with a 5% free stream turbulent intensity.Static pressure was set at the outlet of the draft tube.The interfaces between stator-rotor blocks were set as the standard transient sliding interfaces.Walls were defined with no-slip wall boundary conditions.Runner blades, runner hub and shroud were fixed as rotating walls.The time-step size was set as 0.000139 s, such that one runner revolution was divided into 360 time steps.Second-order implicit time-stepping was adopted for time discretization.The convergence residuals for continuity and momentum equations were below 1.0 × 10 −5 .All calculations were conducted in a cluster computer with eight Intel 5645 2.4 GHz processors, 96 GB RAM, and 2 TB hard drives (Dell Inc., Round Rock, TX, USA).
As shown in Figure 3, at a guide vane opening of A = 17.0 mm, OP1 is the best efficient point (BEP), OP2-OP4 are the operating points in the normal operating range, while OP5 is the runaway point and OP6 is the low-discharge point in the turbine brake.Usually, the guide vane opening is given as an angle; however, its width was used in this study.The opening angle was γ = 18.0 o at A = 17.0 mm.Four meshes (Table 3) were used to test mesh independence.The calculated points were OP2-OP4 (Figure 4).The calculated efficiency errors for OP2-OP4 are given in Table 4. Here, the error was calculated as Err = η cal − η exp /η num × 100.0%.In the simulations, the head cover and stay ring were not included; therefore, some leakage losses and the mechanical losses were not incorporated.The calculated efficiency η cal is mainly the hydraulic efficiency, and is calculated by means of a time average after the calculated unsteady condition becomes stable.Simulated results show that mesh density has a weak influence on efficiency as the grid density is high (3.62 million).After considering the complex unsteadiness, Mesh III was chosen for the simulations.The validation of the numerical results in the next section also shows that the numerical solutions used can provide reliable results.
Table 5 gives more information on Mesh III.There were a total of 4.5 million elements in the whole domain.Figure 3 shows the boundary mesh for different components.Ten element layers with normal stretching ratios were created from the walls for improving the velocity profile from the no-slip Energies 2017, 10, 96 6 of 13 boundaries.Near the wall, the value of no-dimensional distances was y + = 50 − 180.The mesh in the domain of the guide vanes, runner and diffuser section of the draft tube was finer than that in the other flow domain, in order to capture the complex flow characteristics in these components at higher resolution.

Validation of the Numerical Results
Numerical results were compared with experimental measurements (Table 6).In the numerical simulations, the discharge Q, and the runner rotating speed n, were specified with the same values as those in the measurements, while the head H, and the moment M, were time-averaged values within the simulations.The Reynolds number was Re = ωD 2 2 /ν = 2.82 × 10 6 , where ω is the runner rotation angular velocity, D 2 is the runner inlet diameter in turbine mode, and ν is the kinematic viscosity.
The agreement between numerical and experimental results was quite good in terms of n 11 , Q 11 , and M 11 , with an error less than 5.0% at OP1 and OP2-OP4.At OP5 and OP6, there was a significant increase in error (Table 5).As regards moment M, at OP5, its values should be zero, and at OP6, its value should be a small negative value.Minor numerical overestimation in the moment will induce a significant increase in errors.Moreover, mechanical friction losses were not considered in the numerical simulations.These losses are closely related to the runner rotation rate and they remain almost unchanged with the operating condition.Therefore, the corresponding reduction of hydraulic power at these two points will induce significant increases in percentage errors, as shown in Table 6.Even at these two points, the errors in n 11 and Q 11 are still smaller than 5.0%.This means that the hydraulic heads are still well predicted in spite of the simplification of mechanical friction losses.The comparisons demonstrate the capability of the numerical simulations to predict the hydraulic characteristics of the pump-turbine with sufficient accuracy.

Results and Discussion
Pressure monitors in the whole water passage are presented in Figure 5. Monitor SC1 was located on the spiral casing wall.In the gap between the guide vanes and runner, 20 monitors (GV1-GV20) were arranged evenly around the circumference.On the draft tube wall there were three monitors, DT1, DT2 and DT3.Pressure fluctuations on monitors GV5, GV15, and DT1-DT3 were registered in the model test [19,20].In a pump-turbine, pressure fluctuations in turbine mode are usually stronger than those in pump mode.Moreover, fluctuation amplitudes in the draft tube and vaneless gap between the guide vanes and runner are usually larger than those in the spiral casing and stay vane channels.

Results and Discussion
Pressure monitors in the whole water passage are presented in Figure 5. Monitor SC1 was located on the spiral casing wall.In the gap between the guide vanes and runner, 20 monitors (GV1-GV20) were arranged evenly around the circumference.On the draft tube wall there were three monitors, DT1, DT2 and DT3.Pressure fluctuations on monitors GV5, GV15, and DT1-DT3 were registered in the model test [19,20].In a pump-turbine, pressure fluctuations in turbine mode are usually stronger than those in pump mode.Moreover, fluctuation amplitudes in the draft tube and vaneless gap between the guide vanes and runner are usually larger than those in the spiral casing and stay vane channels.
where  p c is the pressure fluctuation coefficient, p is the instantaneous pressure, p is the time-averaged pressure, and E is the specific energy.For pressure fluctuation analyses, in frequency domains, the frequency spectrum was normalized as f / f n , where f is the frequency component and f n is the runner rotating frequency.In time domains, the pressure fluctuations are normalized by Equation (2).
where c p is the pressure fluctuation coefficient, p is the instantaneous pressure, p is the time-averaged pressure, and E is the specific energy.

Pressure Fluctuations in the Draft Tube
Three monitors, DT1-DT3, were located on the conical section in the draft tube.Table 7 shows the pressure fluctuations of monitor DT2 at six different operating conditions, along with the numerical and experimental results.The experimentally-measured pressure fluctuations are only listed in the table for OP2-OP4.This is because pressure fluctuations were not measured at operating points in the S-shaped region.Instead, pressure fluctuations were measured during the efficiency test [19,20].Therefore, only the operating points in the normal operating ranges are given in the table.Although the amplitudes have relative errors of a little more than 13.0% between numerical and experimental results, the first dominant frequency ( f 1 ) is in good agreement.The single helical vortex rope in the draft tube usually rotates at 20%-40% of the runner rotating speed [11].However, the dominant frequencies have higher values (Table 6).These higher frequencies may be induced by a twin vortex rope [14].The pressure fluctuations at monitor DT2 at OP3 are shown in Figure 6a and the frequency spectrum is shown in Figure 6b.A total of 32 runner revolutions were calculated after the unsteady simulation stabilized (Figure 6).The first dominant frequency measured during the experiment was 0.89 f n , compared with 0.86 f n in the simulation.Based on Table 7 and Figure 6, it can be seen that the pressure fluctuations in the draft tube are mainly induced by vortex ropes.

Pressure Fluctuations in the Draft Tube
Three monitors, DT1-DT3, were located on the conical section in the draft tube.Table 7 shows the pressure fluctuations of monitor DT2 at six different operating conditions, along with the numerical and experimental results.The experimentally-measured pressure fluctuations are only listed in the table for OP2-OP4.This is because pressure fluctuations were not measured at operating points in the S-shaped region.Instead, pressure fluctuations were measured during the efficiency test [19,20].Therefore, only the operating points in the normal operating ranges are given in the table.Although the amplitudes have relative errors of a little more than 13.0% between numerical and experimental results, the first dominant frequency ( 1f ) is in good agreement.The single helical vortex rope in the draft tube usually rotates at 20%-40% of the runner rotating speed [11].However, the dominant frequencies have higher values (Table 6).These higher frequencies may be induced by a twin vortex rope [14].The pressure fluctuations at monitor DT2 at OP3 are shown in Figure 6a and the frequency spectrum is shown in Figure 6b.A total of 32 runner revolutions were calculated after the unsteady simulation stabilized (Figure 6).The first dominant frequency measured during the experiment was 0.89 n f , compared with 0.86 n f in the simulation.Based on Table 7 and Figure 6, it can be seen that the pressure fluctuations in the draft tube are mainly induced by vortex ropes.Band-pass filtering was employed to extract the first dominant frequency component of 0.86 n f from the calculated pressure fluctuations of DT1, DT2 and DT3 at OP3.As shown in Figure 7, α stands for the phase difference of pressure between DT1 and DT2, and β represents the phase difference between DT2 and DT3.The α value was near zero, while β was not, which means that a twin vortex rope developed in the draft tube.If the same measurements were made at the other operating points, the twin vortex rope structure could also be detected.
Figure 8 shows the precession of a twin vortex rope at OP3 at a runner rotational period of T .The pressure contours at the planes where monitors DT1 and DT2 were located are also shown.The twin vortices are labelled A and B, and they rotate with the same angular velocity.This velocity Band-pass filtering was employed to extract the first dominant frequency component of 0.86 f n from the calculated pressure fluctuations of DT1, DT2 and DT3 at OP3.As shown in Figure 7, α stands for the phase difference of pressure between DT1 and DT2, and β represents the phase difference between DT2 and DT3.The α value was near zero, while β was not, which means that a twin vortex rope developed in the draft tube.If the same measurements were made at the other operating points, the twin vortex rope structure could also be detected.
alongside cavitation.This figure is simply used to verify the generation of twin vortex ropes because cavitation was not considered, and only a single flow field was simulated in the present study.It is well known that vortex ropes in draft tubes are usually induced by a high residual swirl in the input flow, which often forms a single helical vortex rope.Therefore, further research is required to understand the transition from a single to a double rope, and the influences of geometry on this process.Figure 8 shows the precession of a twin vortex rope at OP3 at a runner rotational period of T. The pressure contours at the planes where monitors DT1 and DT2 were located are also shown.The twin vortices are labelled A and B, and they rotate with the same angular velocity.This velocity corresponds to the dominant frequency of the twin vortex ropes, 0.86 f n .Figure 9 shows a photograph of a vortex rope observed during the experiment.The twin vortex rope is generated alongside cavitation.This figure is simply used to verify the generation of twin vortex ropes because cavitation was not considered, and only a single flow field was simulated in the present study.It is well known that vortex ropes in draft tubes are usually induced by a high residual swirl in the input flow, which often forms a single helical vortex rope.Therefore, further research is required to understand the transition from a single to a double rope, and the influences of geometry on this process.
Energies 2017, 10, 96 9 of 13 corresponds to the dominant frequency of the twin vortex ropes, 0.86 n f .Figure 9 shows a photograph of a vortex rope observed during the experiment.The twin vortex rope is generated alongside cavitation.This figure is simply used to verify the generation of twin vortex ropes because cavitation was not considered, and only a single flow field was simulated in the present study.It is well known that vortex ropes in draft tubes are usually induced by a high residual swirl in the input flow, which often forms a single helical vortex rope.Therefore, further research is required to understand the transition from a single to a double rope, and the influences of geometry on this process.

Pressure Fluctuations in the Vaneless Gap
Table 8 shows the amplitude and dominant frequency of pressure fluctuations on GV5.The experimentally measured pressure fluctuations are only listed in the table for the three operation points OP2-OP4 (as per Table 7).The guide vanes and runner blades interact to generate pressure fluctuations in the gap between them.The dominant frequencies are rotor-stator interaction frequency 7fn and its harmonics, for both simulated and experimental results.The pressure amplitudes in operating conditions of low mass flow rate at OP4-OP6 were much higher than those at OP1-OP3.These increased amplitudes are believed to be induced by serious flow separation at the inlet of the runner.Figure 10 shows the distribution of streamlines at the guide vanes and runner channels under operating conditions OP1, OP5 and OP6.At OP1, the streamlines were well distributed at both guide vanes and runner channels, and there was no obvious flow separation.At OP5, flow separation occurred at the inlet of the runner, leading to backflow cells and vortices in the runner channels.However, the streamlines at the guide vanes still remained uniform.At OP6, backflow cells and vortices entirely filled the runner channels and flow separation occurred at the guide vanes.

Pressure Fluctuations in the Vaneless Gap
Table 8 shows the amplitude and dominant frequency of pressure fluctuations on GV5.The experimentally measured pressure fluctuations are only listed in the table for the three operation points OP2-OP4 (as per Table 7).The guide vanes and runner blades interact to generate pressure fluctuations in the gap between them.The dominant frequencies are rotor-stator interaction frequency 7f n and its harmonics, for both simulated and experimental results.The pressure amplitudes in operating conditions of low mass flow rate at OP4-OP6 were much higher than those at OP1-OP3.These increased amplitudes are believed to be induced by serious flow separation at the inlet of the runner.Figure 10 shows the distribution of streamlines at the guide vanes and runner channels under operating conditions OP1, OP5 and OP6.At OP1, the streamlines were well distributed at both guide vanes and runner channels, and there was no obvious flow separation.At OP5, flow separation occurred at the inlet of the runner, leading to backflow cells and vortices in the runner channels.However, the streamlines at the guide vanes still remained uniform.At OP6, backflow cells and vortices entirely filled the runner channels and flow separation occurred at the guide vanes.

Pressure Fluctuations in the Vaneless Gap
Table 8 shows the amplitude and dominant frequency of pressure fluctuations on GV5.The experimentally measured pressure fluctuations are only listed in the table for the three operation points OP2-OP4 (as per Table 7).The guide vanes and runner blades interact to generate pressure fluctuations in the gap between them.The dominant frequencies are rotor-stator interaction frequency 7fn and its harmonics, for both simulated and experimental results.The pressure amplitudes in operating conditions of low mass flow rate at OP4-OP6 were much higher than those at OP1-OP3.These increased amplitudes are believed to be induced by serious flow separation at the inlet of the runner.Figure 10 shows the distribution of streamlines at the guide vanes and runner channels under operating conditions OP1, OP5 and OP6.At OP1, the streamlines were well distributed at both guide vanes and runner channels, and there was no obvious flow separation.At OP5, flow separation occurred at the inlet of the runner, leading to backflow cells and vortices in the runner channels.However, the streamlines at the guide vanes still remained uniform.At OP6, backflow cells and vortices entirely filled the runner channels and flow separation occurred at the guide vanes.Figure 11 presents the calculated pressure fluctuations and frequency spectrum at GV5 for operation conditions OP1, OP5 and OP6.For these three points, a low frequency component of 0.42fn was captured at OP1, and 0.65fn was well captured at OP5 as well as OP6, except for the blade passing frequency and its harmonics.The frequency component 0.42fn at OP1 was about half of the vortex rope frequency, according to Table 6.Consequently, the component 0.42fn was the rotation frequency of a single helical vortex rope.The frequency 0.65fn was also well captured at the other 19 monitors in the gap between the runner and guide vanes at OP5 and OP6 (Figure 12).However, as shown in Table 7, the frequency of the twin vortex rope is different at OP5 and OP6.Therefore, the frequency 0.65fn cannot be the rotation frequency of a single helical vortex rope at OP5 and OP6. Figure 12 shows the band-pass-filtered pressure signals of 0.65fn for monitors GV1 to GV20 at OP5 and OP6, respectively.It is mentioned by Brennen [21] that the rotating stall phenomenon may occur in one, or in a cascade of rotor or stator blades, operating at a high incidence angle.Usually, the stall appears on a few adjacent runner channels and stall cells propagate in the circumferential direction.In Figure 12, the propagation of stall cells is represented by parallel lines, and they have the same direction as the rotating direction of the runner.At OP5, there is one single stall cell in the gap between the runner and guide vanes, while at OP6, there are four stall cells.The band-pass-filtered pressure signals for the four circumferentially-distributed monitors, GV5, GV10, GV15 and GV20, are redrawn in Figure 13 for OP5 and OP6, respectively.At OP5, the pressure signals of GV5 and GV15 had a phase difference of π, which is the same as for monitors GV10 and Figure 11 presents the calculated pressure fluctuations and frequency spectrum at GV5 for operation conditions OP1, OP5 and OP6.For these three points, a low frequency component of 0.42f n was captured at OP1, and 0.65f n was well captured at OP5 as well as OP6, except for the blade passing frequency and its harmonics.The frequency component 0.42f n at OP1 was about half of the vortex rope frequency, according to Table 6.Consequently, the component 0.42f n was the rotation frequency of a single helical vortex rope.The frequency 0.65f n was also well captured at the other 19 monitors in the gap between the runner and guide vanes at OP5 and OP6 (Figure 12).However, as shown in Table 7, the frequency of the twin vortex rope is different at OP5 and OP6.Therefore, the frequency 0.65f n cannot be the rotation frequency of a single helical vortex rope at OP5 and OP6. Figure 11 presents the calculated pressure fluctuations and frequency spectrum at GV5 for operation conditions OP1, OP5 and OP6.For these three points, a low frequency component of 0.42fn was captured at OP1, and 0.65fn was well captured at OP5 as well as OP6, except for the blade passing frequency and its harmonics.The frequency component 0.42fn at OP1 was about half of the vortex rope frequency, according to Table 6.Consequently, the component 0.42fn was the rotation frequency of a single helical vortex rope.The frequency 0.65fn was also well captured at the other 19 monitors in the gap between the runner and guide vanes at OP5 and OP6 (Figure 12).However, as shown in Table 7, the frequency of the twin vortex rope is different at OP5 and OP6.Therefore, the frequency 0.65fn cannot be the rotation frequency of a single helical vortex rope at OP5 and OP6. Figure 12 shows the band-pass-filtered pressure signals of 0.65fn for monitors GV1 to GV20 at OP5 and OP6, respectively.It is mentioned by Brennen [21] that the rotating stall phenomenon may occur in one, or in a cascade of rotor or stator blades, operating at a high incidence angle.Usually, the stall appears on a few adjacent runner channels and stall cells propagate in the circumferential direction.In Figure 12, the propagation of stall cells is represented by parallel lines, and they have the same direction as the rotating direction of the runner.At OP5, there is one single stall cell in the gap between the runner and guide vanes, while at OP6, there are four stall cells.The band-pass-filtered pressure signals for the four circumferentially-distributed monitors, GV5, GV10, GV15 and GV20, are redrawn in Figure 13 for OP5 and OP6, respectively.At OP5, the pressure signals of GV5 and GV15 had a phase difference of π, which is the same as for monitors GV10 and Figure 12 shows the band-pass-filtered pressure signals of 0.65f n for monitors GV1 to GV20 at OP5 and OP6, respectively.It is mentioned by Brennen [21] that the rotating stall phenomenon may occur in one, or in a cascade of rotor or stator blades, operating at a high incidence angle.Usually, the stall appears on a few adjacent runner channels and stall cells propagate in the circumferential direction.In Figure 12, the propagation of stall cells is represented by parallel lines, and they have the same direction as the rotating direction of the runner.At OP5, there is one single stall cell in the gap between the runner and guide vanes, while at OP6, there are four stall cells.The band-pass-filtered pressure signals for the four circumferentially-distributed monitors, GV5, GV10, GV15 and GV20, are redrawn in Figure 13 for OP5 and OP6, respectively.At OP5, the pressure signals of GV5 and GV15 had a phase difference of π, which is the same as for monitors GV10 and GV20 because of a single stall cell.At OP6, as there were four stall cells all together, and there was no definite phase difference among the four monitors.
Energies 2017, 10, 96 12 of 13 GV20 because of a single stall cell.At OP6, as there were four stall cells all together, and there was no definite phase difference among the four monitors.

Conclusions
This study investigated flow in a reversible pump-turbine model.The unstable characteristics along a constant guide vane opening were simulated, and the numerical results were compared with experimental measurements.Analyses focused on the pressure fluctuations in the draft tube and in the gap between the runner and guide vanes.The main conclusions are summarized as: (1) In the draft tube, a twin vortex rope formed under all the operating conditions that were investigated.The twin vortex rope rotated at 0.75-1.19times the runner rotation speed, which was much faster than for single vortex ropes.(2) In the gap between the runner and guide vanes, pressure fluctuations were mainly caused by rotor-stator interactions and vortices due to flow separation.A low frequency component 0.65fn was well captured at both runaway and low-discharge points in the turbine brake zone, and was the effect of a rotating stall phenomenon.At the runaway point, a single stall cell was found in the gap, while at the low-discharge point, four stall cells were found.
parameters of the prototype pump-turbine are as follows.In turbine mode, the rated head is Hr = 259 m and the rated output power is Pr = 255 MW.The rotational speed is 300 rpm r n = for both turbine and pump modes.In pump mode, the maximum and minimum heads are max 298 m H = and min 239 m H = , respectively.

Figure 5 .
Figure 5. Pressure monitors in the whole passage.

For
pressure fluctuation analyses, in frequency domains, the frequency spectrum was normalized as n f f , where f is the frequency component and n f is the runner rotating frequency.In time domains, the pressure fluctuations are normalized by Equation (2).

Figure 5 .
Figure 5. Pressure monitors in the whole passage.

Figure 9 .
Figure 9. Photograph of a twin vortex rope taken during the experiment, with ropes 1 and 2 labelled.

Figure 9 .
Figure 9. Photograph of a twin vortex rope taken during the experiment, with ropes 1 and 2 labelled.

Table 1 .
Parameters of pump-turbine model in turbine mode.

Table 1 .
Parameters of pump-turbine model in turbine mode.

Table 2 .
Characteristics of test rig performance and pressure sensors.

Table 2 .
Characteristics of test rig performance and pressure sensors.

Table 5 .
Statistics and quality of the mesh.BEP: best efficient point.

Table 6 .
Investigated operating points and their values of n 11 , Q 11 , and M 11 .

Table 6 .
Investigated operating points and their values of n11, Q11, and M11.

Table 7 .
Amplitude and dominant frequency of pressure fluctuations at DT2.

Table 7 .
Amplitude and dominant frequency of pressure fluctuations at DT2.

Table 8 .
Amplitude and dominant frequency of pressure fluctuations at GV5.

Table 8 .
Amplitude and dominant frequency of pressure fluctuations at GV5.

Table 8 .
Amplitude and dominant frequency of pressure fluctuations at GV5.