## Abstract

It is more dangerous for the units sharing the same division pipeline system than those sharing individual system during the load rejection transient for a high water head pumped storage station. A new mathematical model for the load rejection transient in these two units sharing the same division pipeline was proposed using method of internal characteristic (MIC). The transient relationship between the geometric parameters and working parameters was established based on the dynamic equations of the hydraulic machinery. Meanwhile, the calculation of the load rejection transient was accomplished and compared with field experimental data from a pumped storage power station. The relative errors of the maximum rotational speed of the two units between the predicted results and the experimental data were −5.9% and −4.44%, respectively, and the relative errors of the maximum rising of hydraulic pressure were −6.7% and −6.98%, respectively. The results met with the design specification requirements, which indicated that the proposed new MIC could completely predict the load rejection transient in these two units with joint closing law of ball-valve and guide vane, and the corresponding calculated results were acceptable and reliable. The MIC is also applicable to the transient investigation of hydraulic systems with hydraulic machinery.

## NOMENCLATURE

*M*Moment, N·m

*ρ*Fluid mass density, kg/m

^{3}*W*Volume, m

^{3}*b*Infinitesimal thickness, m

*θ*Infinitesimal angle, °

*m*Infinitesimal axial meridian, m

*r*Radius of micro fluid to center of shaft, m

*v*_{u}Circumference component of absolute velocity, m/s

*t*Time, s

*v*_{m}Axial meridian velocity, m/s

*Q*Flow discharge, m

^{3}/s*v*_{u}_{1}Circumference component of absolute velocity at inlet, m/s

*v*_{u}_{2}Circumference component of absolute velocity at outlet, m

*r*_{1}Inlet radius of pump-turbine, m

*r*_{2}Outlet radius of pump-turbine, m

*M*_{H}Transient moment acting on runner, N·m

*Q*_{H}Transient flow discharge, m

^{3}/s*H*_{H}Turbine transient water head, m

*g*Gravitational acceleration, m/s

^{2}*η*Hydraulic efficiency

*ω*_{H}Transient angular speed, rad/s

*H*_{eH}Transient effective water head, m

*u*_{1}Circumference velocity at inlet, m/s

*u*_{2}Circumference velocity at outlet, m/s

*Γ*_{1}Velocity circulation at inlet, m

^{2}/s*Γ*_{2}Velocity circulation at outlet, m

^{2}/s*α*Angle between outlet of guide vane and circumferential direction, °

*b*_{0}Height of guide vane, m

*K*_{p}Cascade transparent coefficient

*l*Length of hydrofoil of cascade, m

*t*_{c}Pitch of hydrofoil of cascade, m

*r*_{a}Effective radius of runner cascade, m

*i*_{0}Coefficient of cascade in zero lift force

*β*_{2}Blade angle at outlet of meridian plane of runner, °

*F*Area, m

^{2}*q*Flow discharge of layer thickness, m

^{3}/s*v*_{r}Radial component of absolute velocity, m/s

*ω*Angular speed of unit, rad/s

*M*_{D}Driving moment, N·m

*M*_{R}Resisting moment, N·m

*J*Inertia moment of rotating parts, kg·m

^{2}*G*Weight of rotating parts, N

*D*Distance from center of mass of rotating parts to axle, m

*ω*_{0}Initial angular speed, rad/s

*K*Mockridge shaft moment coefficient

*M*_{M}Mockridge moment, N·m

*b*_{p}Blade width at outlet of pump, m

*D*_{p}Runner diameter at outlet of pump, m

*n*Rotational speed of unit, r/min

*Y*_{ul}Water level of upstream reservoir, m

*a*_{w}Wave velocity, m/s

*A*Cross-sectional area of pipe, m

^{2}*H*_{P}Transient water head in P point, m

*Q*_{p}Transient flow discharge in P point, m

*H*_{v}_{0}Initial pressure at inlet of ball-valve at full opening, m

*ξ*Pressure rising rate

*Q*_{cb}Static flow discharge of ball-valve, m

^{3}/s*Δx*Space increment, m

*f*Darcy–Weisbach friction factor

*Q*_{c}Static flow discharge, m

^{3}/s*y*Relative guide vane opening

*Q*_{0}Static flow discharge at full opening of guide vane, m

^{3}/s*a*Guide vane opening, °

*a*_{0}Full opening of guide vane, °

*k*Exponent

*H*_{zp}Transient water head of device, m

*H*_{z}_{0}Static water head of device, m

*D*_{1}Runner diameter of pump-turbine, m

*n*_{H}Transient rotational speed, r/min

*n*_{0}Initial rotational speed, r/min

*k*_{q}Rate of change of unit discharge with respect to unit speed

*H*_{0}Static water head at spiral case inlet, m

*σ*_{T}Characteristic coefficient of draft tube

*L*_{dr}Length of draft tube, m

*v*_{dr}Initial average velocity in draft tube, m/s

*t*_{s}Closing time of guide vane, s

*b*_{0}Height of guide vane, m

*K*_{1}Iterative coefficient

*c*_{1}Coefficient

*c*_{2}Coefficient

*K*_{2}Coefficient

*Q*_{T}Flow discharge at outlet of downstream surge tank, m

^{3}/s*Q*_{S}Flow discharge of downstream surge tank, m

^{3}/s*hs*Depth of downstream surge tank, m

*H*_{dS}Water head of downstream surge tank, m

*H*_{S}_{0}Water head at bottom of downstream surge tank, m

*f*_{S}Friction factor of downstream surge tank

*D*_{S}Diameter of downstream surge tank, m

*A*_{S}Cross-sectional area of downstream surge tank, m

^{2}*Y*_{dl}Water level of downstream reservoir, m

## INTRODUCTION

During load rejection transients, the sharp increment and the fluctuation of pipe pressure endanger the safety of the pumped storage power station (Calamak & Bozkus 2013), especially for two units sharing the same diversion pipeline (Yu *et al.* 2014). The ‘S’ and hump areas in full characteristic curves cause the unsteady dynamic process of pump-turbines with low specific speed under the transient process (Li *et al.* 2016; Zeng *et al.* 2017). Since the negative installation elevation of a draft tube is usually designed to prevent it from cavitation for high water head pumped storage power station, the pressure at the inlet of the draft tube will be much lower when load rejection occurs, which may cause reverse water hammer in the draft tube and result in huge damage of the device.

Developing a reasonable closing law of guide vanes is an effective method to reduce the pressure pulsation and the rotational speed increment (Fang & Koutnik 2012). Yao *et al.* (2013) investigated the key effect factor for the load rejection in pump-turbines by comparing five different guide closing laws. Vakil & Firoozabadi (2009) analyzed the effect of different multi-line closing laws, and found an optimum closing law for pressure pulsation reduction. The closing law based on ball-valve and guide vane was applied for the load rejection in a pumped storage power station (Zhang & Yang 2011). A numerical simulation of the closing law based on ball-valve and guide vane was adopted by Hou & Cheng (2005), where the results proved that the closing law can ensure the safety of the unit during load rejection. However, more tests were still required to figure out the effects of the closing law on pump-turbine units using joint adjustment of ball-valve and guide vane during load rejection. The one-dimensional method of characteristics (MOC) is widely used in analyzing the transient flow for pipeline systems due to its high efficiency (Freni *et al.* 2013; Acosta *et al.* 2015). Hwang *et al.* (2012) developed a novel particle method of characteristics to simulate unsteady pipe flows and the method was validated by solving some benchmark problems with significant transient effects in single pipe flows. Based on the MOC and the modified Suter transformations, Rezghi & Riasi (2016) simulated the simultaneous operation of two units at runaway.

The transition process of hydraulic machinery should consider unsteady flow characteristics in the hydraulic machinery while the traditional method employs the steady characteristics. The previous investigations on transient are applied with the method of external characteristics, which is based on the full characteristic curves of pump-turbines (Wylie & Streeter 1993; Rezghi & Riasi 2016). However, when the full characteristic curves obtained from a model turbine test in steady state are employed to calculate the transient state, this technique will cause an error. At the same time, the study on the closing law is insufficient, and its influence on dynamic performance in the transition should be further studied.

In this paper, the double units’ load rejection transient with the joint closing law of ball-valve and guide vane was calculated for a pumped storage power station. The numerical simulation of double load rejection was carried out with a new mathematical model which was based on method of internal characteristic (MIC) defining the boundary conditions such as pump-turbine, ball-valve and surge tank. In addition, the load rejection transient of a pumped storage power station was predicted and compared with field experimental data, in order to explore the characteristics of the flow discharge control mode.

## MIC FOR PUMP-TURBINE

The Euler equation is the energy conversion relation with the kinematics parameters at the inlet and outlet of the runner for turbine stability conditions. But the basic parameters such as load, flow discharge, and rotational speed always change frequently in the transition. Therefore, the relationship between the flow dynamics and kinematics parameters must be established for transient conditions. With the calculation MIC for the pump-turbine unit, the influences of the geometry parameters on the kinematics parameters of hydraulic machinery are analyzed, and the nonlinear equations for these parameters and the position adjustment are set up. As long as the law of the opening change of the adjusting element with time is given, the transient rules of dynamic parameters will be obtained using numerical calculation MIC for the transient solution of the nonlinear equations.

### General basic equations of hydraulic machinery

Under transient conditions, the flow in the runner of the Francis pump-turbine is shown in Figure 1 (Chang 2005).

*M*is the moment,

*ρ*is the fluid mass density,

*W*is the volume, d

*W*

*=*

*r*d

*θ*d

*m*d

*b*is the volume of micro fluid layer,

*b*,

*θ*, and

*m*are the infinitesimal thickness, angle and axial meridian, respectively,

*r*is the radius of the micro fluid to the center of the shaft,

*v*is the circumference component of the absolute velocity,

_{u}*v*is the time-dependent velocity moment, which is related to meridian streamline and polar coordinates (m,

_{u}r*θ*), which has no relationship to the polar coordinate for axisymmetric flow. The derivative of

*v*with respect to the time

_{u}r*t*, is obtained as: where

*v*=

_{m}*Q*/(2

*πr*d

*b*) is the axial meridian velocity, and

*Q*is the flow discharge.

*v*and d

_{m}*W*into the second of Equation (4), we get: where

*v*

_{u}_{1}and

*v*

_{u}_{2}are the circumference components of the absolute velocities at the inlet and outlet, respectively,

*r*

_{1}and

*r*

_{2}are the inlet and outlet radius of pump-turbine, respectively.

Equation (10) is the basic equation of water-turbine showing the conversion relationship between the flow energy and the mechanical energy. First, flow channel forces the fluid moment of momentum to change. Second, the fluid transmits the energy on the runner by the force on blades. Equation (10) is also called the Euler equation, building the relationship between the energy conversion and the basic flow kinematics parameters at the runner inlet and outlet under steady conditions. The Euler equation states that the energy conversion condition is the velocity moment at inlet and outlet of runner change.

The Euler equation is only applicable to the steady operations of the turbine. In the load rejection transient, the flow discharge and the rotational speed of the unit change dynamically and lead to non-zero value in the second item of Equation (8). At this time, the Euler equation is no longer applicable, especially for transient operations of the turbine, the operation parameters such as the flow velocity and the rotational speed of the unit change greatly with time because the second item is comparable to the first item and cannot be omitted. The transient moment and transient water head must be solved by Equations (6) and (8) which are called the generalized basic equations. Equations (9) and (10) are special cases of Equations (6) and (8), which are the basic expressions of transient moment and transient effective head for turbine, respectively.

The working state of the unit under pump condition is opposite to the state under turbine condition. As long as a minus sign is added to the right parts, Equations (6) and (8) will be the transient moment and the transient effective water head under pump condition. Equations (6) and (8) are also called the basic equations of hydraulic machinery, and are applied both under steady and unsteady conditions, especially under the transient for the unit.

*Γ*

_{1}and

*Γ*

_{2}are the velocity circulations at inlet and outlet,

*α*is the angle between the outlet of guide vane and the circumferential direction,

*b*

_{0}is the height of guide vane,

*K*= e

_{p}^{−(}

^{πl}^{/tc)}is the cascade transparent coefficient,

*l*and

*t*are the length and pitch of the hydrofoil of the cascade, respectively,

_{c}*r*is the effective radius of the runner cascade,

_{a}*i*

_{0}= −2π

*r*

_{2}cot

*β*

_{2}/

*F*

_{2}is the coefficient of the cascade in zero lift force,

*β*

_{2}is the blade angle at the outlet on the meridian plane of the runner,

*F*is the area,

*q*= 2π

*rv*is the flow discharge of layer thickness,

_{r}*v*=

_{r}*Q/F*is the radial component of the absolute velocity, and

*ω*is the angular speed of the unit.

### Motion equations of rotating parts

*M*and resisting moment

_{D}*M*, namely: where

_{R}*J*

*=*

*GD*

^{2}/(4

*g*) is the inertia moment of the rotating parts,

*G*is the weight of the rotating parts,

*D*is the distance from the center of mass of the rotating parts to the axle.

As shown in Figure 2, the load rejection process of the pump-turbine unit is composed of turbine, turbine-braking, and reverse pump modes. Then, the moment *M* is obtained by using Equation (16) and Figure 2.

#### Turbine mode

*M*is positive, and the motor-generator works in a generator mode when

_{D}*M*is negative, where

_{R}*M*only represents the mechanical friction resistance moment because the generator is not connective with the grid under load rejection condition. Meanwhile,

_{R}*M*is too small and can be ignored, so the angular speed is expressed as: where

_{R}*ω*

_{0}is the initial angular speed.

#### Turbine-braking mode

*K*and the ratio of the blade width

*b*to the runner diameter at outlet

_{p}*D*under zero flow discharge condition in 1943 (Stepanoff 1993), as shown by the dotted line 1 in Figure 3. Yu (1998) gave the relation curve between

_{p}*K*and

*b*/

_{p}*D*of Chinese pumps, which is shown by the solid line 2 in Figure 3, and also presented the fitting equation:

_{p}Under the same *b _{p}*/

*D*, the

_{p}*K*obtained by Yu (1998) is larger than that obtained by Mockridge. The main reason is that the braking torque obtained by Yu (1998) was drawn from the statistical methods, and thus it includes the mechanical torque. The Mockridge curve considered the secondary flow further by theoretical analysis, and thus is more appropriate for pump mode. In the transient process, the fluid flow in the pump-turbine runner is complex, and the secondary flow is serious. At the same time, the Mockridge curve obtained by theoretical analysis is more universal. Therefore, the Mockridge curve is employed during transient simulation.

#### Reverse pump mode

*M*on the runner acted by water becomes resistant and is replaced by −

_{D}*M*. The motor-generator works in a motor mode. Then, under the reverse pump condition, the angular speed is expressed as:

_{D}Although the proposed method in the present study is confined in a pipeline with a pump-turbine, it can easily be extended for simulation of transients in the hydraulic system with water turbine or pump as long as the flow direction and the signs of the torque and the head are changed accordingly.

The MIC expressions are based on the generalized equations of hydraulic machinery under unstable conditions, such as transient moment, transient water head, and transient axial force. The MIC can be used to calculate accurately the dynamic process of the turbine prototype device because it includes dynamic additional terms caused by the water inertia in the runner. The transition calculation can be completed without the full characteristic curve of the vane type hydraulic machinery. Thus, the MIC expressions can be used to determine the key structure parameters during the feasibility study and preliminary technological design stages for the hydropower stations, especially for large hydropower stations which have no complete characteristic curves. At the same time, the method of combining the MIC with the characteristic line theory of pressure transients in elastic piping systems is much simpler than the traditional method based on the complete characteristic curve of hydraulic turbines because the interpolation of the complete characteristic curve is necessary for traditional methods.

## NEW MODEL OF LOAD REJECTION BASED ON MIC

The internal characteristic theory describes the transient rules of the transient head, the transient rotational speed, and the transient flow discharge (Chang 2005), and then the related MIC together with the MOC defines the boundary condition of the pump-turbine system. Figure 4 is the schematic layout of two units sharing the same division pipeline, which consists of an upper reservoir, diversion tunnel, ball-valve, tailrace branches, downstream surge tank, tailrace tunnel, and downstream reservoir. The pipeline system consists of nine parts which are labeled *l*_{1} to *l*_{9}. Water flows from upstream reservoir (Node A), through the pipeline *l*1 to the Node B, and then it is divided into two branches which meet at Node G. After that, it flows through the pipeline *l*8, the surge tank (Node H) and the pipeline *l*9 into the down reservoir (Node I). One branch is composed of pipelines *l*2–*l*4, the valve (Node C) and the unit 1 (Node D), the other is composed of pipelines *l*5–*l*7, the valve (Node E) and the unit 2 (Node F). The boundary conditions are set as follows.

### Boundary condition of upstream reservoir

*Y*is the water level of the upstream reservoir, the subscript

_{ul}*i*is the section number,

*a*is the wave velocity,

_{w}*A*is the cross-sectional area of the pipe,

*H*and

_{P}*Q*are the transient water head and flow discharge in P point, respectively, where the first subscript is the pipeline number and the second subscript is the section number of the pipe.

_{p}### Boundary condition of branch pipes

### Boundary condition of ball-valve

*H*

_{v}_{0}is the initial pressure at the inlet of the ball-valve at full opening,

*ξ*is the pressure rising rate,

*Q*is the static flow discharge of the ball-valve,

_{cb}*Δx*is the space increment, and

*f*is the Darcy–Weisbach friction factor.

### Boundary condition of pump-turbine (MIC)

Generally, during the load rejection process, the pump-turbine can track all the three operating conditions. Calculation formulas of corresponding transient moment and transient flow discharge based on MIC for different operating conditions are presented (Node D and Node F in Figure 4).

#### Turbine mode

*Q*

_{c}*=*

*Q*

_{0}

*y*is the static flow discharge,

^{k}*Q*

_{0}is the static flow discharge at full opening of the guide vane,

*y*

*=*

*a*/

*a*

_{0}is the relative opening of the guide vane,

*a*is the opening of the guide-vane,

*a*

_{0}is the full opening of the guide vane,

*k*is the exponent which is dependent on the classification of turbine and the specific speed of turbine,

*H*is the transient water head of the device,

_{zp}*H*

_{z}_{0}is the static water head of the device,

*D*

_{1}is the runner diameter of pump-turbine,

*n*=

_{H}*ω*is the transient rotational speed,

_{H}r*n*

_{0}is the initial rotational speed,

*k*is the rate of change of unit discharge with respect to unit speed and is determined by specific speed,

_{q}*k*< 0, when pump-turbine is in a turbine mode,

_{q}*H*

_{0}is the static water head at the spiral case inlet,

*σ*

_{T}=

*L*

_{dr}

*v*

_{dr}/(

*gH*

_{z}_{0}

*t*

_{s}) is the characteristic coefficient of the draft tube,

*L*

_{dr}is the length of the draft tube,

*v*

_{dr}is the initial average velocity in the draft tube,

*t*

_{s}is the closing time of the guide vane, and

*b*

_{0}is the height of the guide vane.

#### Turbine-breaking mode

#### Reverse pump mode

### Boundary condition of parallel pipes

### Boundary condition of downstream surge tank

*Q*is the flow discharge at the outlet of the downstream surge tank,

_{T}*Q*is the flow discharge of the downstream surge tank,

_{S}*h*is the depth of the downstream surge tank,

_{S}*H*is the water head of the downstream surge tank,

_{dS}*HS*

_{0}is the water head at the bottom of the downstream surge tank,

*f*

_{S}is the friction factor of the downstream surge tank,

*D*

_{S}and

*A*

_{S}are the diameter and cross-sectional area of the downstream surge tank, respectively.

### Boundary condition of downstream reservoir

*Y*is the water level of the downstream reservoir.

_{dl}### Numerical method for solving MIC equations

The transient process is calculated as the unit operation conditions, and the ball-valve closing law and the guide vane closing law are given. The solving processes for the pump-turbine and the ball-valve using the MIC are briefly described as follows: (1) according to the closing laws of the guide vane and ball-valve, the initial steady angular speed of the unit is determined; (2) the steady flow inside the ball-valve is determined; (3) according to the movement of the guide vane, the steady flow inside the unit is calculated, and then the transient flow discharge is obtained; (4) the values of d*q*/d*t* and d*ω*/d*t* are calculated, where the linear interpolations of d*q*/dt and d*ω*/dt are adopted; (5) the transient torque is calculated; (6) the steady angular speed and transient angular speed are calculated, respectively; (7) the calculated angular speed is compared with the initial value. If the error does not meet the requirement, the calculated value will be taken as the initial value to calculate again until the value meets the requirement.

The MOC is employed to solve the equations along the characteristic line, where the central difference scheme is used for the flow discharge term and the pressure term and the trapezoidal integration method adopted for the rotational speed. The upstream reservoir is calculated by C− equation and the downstream reservoir is solved by C+ equation because the water level is constant for large reservoirs. Based on the above equations, the load rejection transition process is calculated by applying the MIC and MOC.

## CASE STUDIES AND ANALYSES

### Basic parameters of pump-turbine units

The numerical simulation of a simultaneous load rejection process of two units based on MIC was carried out in this paper. The field data of the load rejection transients, which were obtained from the pumped storage power station, were employed to validate the above-mentioned relevant models and methods. The numerical simulation of a simultaneous load rejection process of two units with nominal net head was compared with a field test. The schematic layout of the division pipeline is given in Figure 4. The basic characteristics of these two pump-turbine units are given in Table 1 and the parameters before load rejecting in Table 2. The flow discharge was regulated by the ball-valve and guide vane in the load rejection transient process. Figure 5 shows the closing laws of the guide vane and the ball-valve. During the first 11 seconds of the load rejection transient process, the guide vane was in the original position. From 11 to 26 seconds, the guide vane was closed linearly to zero opening from the initial opening. Finally, the guide vane was at zero opening. During the first 60 seconds of the load rejection transient process, the ball-valve was closed linearly to zero opening from the original position. Finally, the ball-valve was at zero opening.

Parameter | Value |
---|---|

Flow discharge (m^{3}/s) | 80.1 |

Nominal net head (m) | 308 |

Turbine output (MW) | 306.12 |

Rotational speed (r/min) | 333.3 |

Runner diameter (m) | 4.565 |

Guide vane number | 20 |

Guide vane height (m) | 0.5714 |

Parameter | Value |
---|---|

Flow discharge (m^{3}/s) | 80.1 |

Nominal net head (m) | 308 |

Turbine output (MW) | 306.12 |

Rotational speed (r/min) | 333.3 |

Runner diameter (m) | 4.565 |

Guide vane number | 20 |

Guide vane height (m) | 0.5714 |

Parameter | Value |
---|---|

Output (MW) | 300 |

Flow discharge (m^{3}/s) | 98.14 |

Water head (m) | 328 |

Rotational speed (r/min) | 333.3 |

Guide vane closing time (s) | 26 |

Guide vane delay time (s) | 11 |

Ball-valve closing time (s) | 60 |

Parameter | Value |
---|---|

Output (MW) | 300 |

Flow discharge (m^{3}/s) | 98.14 |

Water head (m) | 328 |

Rotational speed (r/min) | 333.3 |

Guide vane closing time (s) | 26 |

Guide vane delay time (s) | 11 |

Ball-valve closing time (s) | 60 |

### Measured and predicted data

Figures 6 and 7 show the water hammer data which were measured and predicted for unit 1 and unit 2 during the simultaneous load rejection process of the two units, respectively. At the beginning stage (0–10 s), the pressure at the spiral case inlet increased and that at the draft tube inlet decreased because of the rapid decrease of the discharge inside the runner, which was caused mainly by the centrifugal effect of the flowing fluid inside the runner of the pump-turbine during the rise of the rotational speed. When the change of the discharge reached the maximum, the pressure rose to the first extremum (point D) which lagged nearly 3 seconds behind the occurring time of the maximum rotational speed, because the centrifugal effect reached the maximum when the rotational speed reached extreme value while the discharge drop caused by the ball-valve was not large due to the slow change of flow area of the ball-valve at its large opening. Then, the pressure at the spiral case inlet dropped to the extremum (point E) and the pressure at the draft tube inlet rose to the maximum because the change of the discharge decreased, which was caused by the coupling effect of both the ball-valve and the rotational speed. After that, the pressure at the spiral case inlet increased again to the second maximum because the discharge decreased, which was caused by the guide vane closing. However, the first extremum was larger than the second extremum, which indicated that the centrifugal effect was obvious for the pump-turbine with low specific speed. Henceforth, the pressure fluctuations began to flatten out due to the resistance of the pipeline.

### Analyses of measured and predicted data

It can be seen from Figures 6 and 7 that the predicted pressure at the spiral case inlet was in good agreement with the experimental data at the beginning of the pressure rising process. The predicted maximum pressure at the spiral case inlet (point D) appeared nearly 1 second ahead of the occurring time of the measured maximum, and the predicted minimal pressure at the spiral case inlet (point E) appeared nearly 5 seconds ahead of the occurring time of minimum measured value, because the rotational speed rise caused the decrement of flow discharge and the pressure pulsation at the spiral case inlet changed with the discharge. The change trend of the predicted pressure at the draft tube inlet was similar to the measured one, but there was a little difference between the predicted and the measured extremums. The measured and predicted minimum relative pressures at the draft tube inlet of unit 1 were 25.6% and 18.9%, respectively, and the corresponding values of unit 2 were 25.86% and 18.88%, respectively. The relative errors of the maximum rise of hydraulic pressure at the spiral case inlet were −6.7% and −6.98%, respectively. The measured and predicted minimum absolute pressures at the draft tube inlet of unit 1 were 32.03 mH_{2}O and 39.59 mH_{2}O, respectively, and the corresponding values of unit 2 were 33.13 mH_{2}O and 41.74 mH_{2}O, respectively. The relative errors of the maximum rise of hydraulic pressure were 7.56% and 8.61%, respectively.

During the first 7 seconds of the load rejection process, the rotational speed of the unit had been rising because the driving moment of the fluid acting on the runner was greater than the resisting moment, and at this time, the pump-turbine unit worked in a turbine mode. When the driving moment equaled the resisting moment, the rotational speed reached the first extremum (point A) which appeared nearly 4 seconds ahead of the start time of the guide vane closing. Then, the unit operation moved into turbine-breaking mode and the rotational speed started to decrease to the extremum (point B), which was caused by the braking moment and the hydraulic moment due to the positive water hammer. Henceforth, the rotational speed gradually fluctuated to zero because of the joint actions of the above-mentioned two types of moments.

At the beginning of the rotational speed rising process, the predicted rotational speed was close to the experimental data. The maximum relative rotational speeds of the predicted and the predicted results of unit 1 were 40.6% and 34.7%, respectively, and the corresponding values of unit 2 were 39.1% and 34.66%, respectively. The relative errors of the maximum rotational speed of the two units between the predicted results and the experimental data were −5.9% and −4.44%, respectively. After the rotational speed reached the maximum, the predicted rotational speed decreased rapidly and the minimum of the predicted rotational speed (point B) appeared nearly 3 seconds ahead of the minimum of the measured rotational speed. The second extremum of the predicted rotational speed (point C) was bigger than the measured one in turbine-breaking mode because the Mockridge moment was possibly added on the runner earlier during the transient process. Thus, it is very important to investigate the suitable time to add the Mockridge moment when using the MIC in the future.

Figure 8 shows the field experimental curves for unit 1 and unit 2 during the simultaneous load rejection process of the two units. The change trend of water hammer data was similar for unit 1 and unit 2, because these two pump-turbine units had the same basic characteristics and similar hydraulic conditions which led to similar pressure pulsations and rotational speed rising. There was a little difference between the change curves of water hammer data including the phase and amplitude because the tailrace branch of unit 1 was much longer than unit 2.

In order to carry out clearer comparisons, the predicted and the measured extremums of two units during the simultaneous load rejection process of the two units are summarized in Table 3. The absolute errors of the maximum rotational speed and the maximum pressure at the spiral case inlet were both small, while the absolute error of the pressure at the draft tube inlet was relatively large. The maximum pressures at the spiral case inlet were far less than the maximum allowed value (500 mH_{2}O), and the minimum pressures at the draft tube inlet were larger than the minimum allowed value (20 mH_{2}O). The maximum of the rotational speed rise ratios was less than the maximum allowed value (45%). Thus it was concluded that the joint closing law of ball-valve and guide vane for two units could prevent the pressure pulsation and the rotational speed rise effectively.

Parameter | Unit | Measured | Predicted | Absolute error |
---|---|---|---|---|

Maximum rotational speed | ||||

Absolute value (r/min) | 1 | 468.66 | 449.01 | −19.65 |

2 | 463.59 | 448.87 | −14.72 | |

Rise ratio value (%) | 1 | 40.6 | 34.7 | −5.9 |

2 | 39.1 | 34.66 | −4.44 | |

Maximum pressure at spiral case inlet | ||||

Absolute value (mH_{2}O) | 1 | 464.35 | 439.54 | −24.81 |

2 | 465.32 | 439.52 | −25.8 | |

Rise ratio value (%) | 1 | 25.6 | 18.9 | −6.7 |

2 | 25.86 | 18.88 | −6.98 | |

Minimum pressure at draft tube inlet | ||||

Absolute value (mH_{2}O) | 1 | 32.03 | 39.59 | 7.56 |

2 | 33.13 | 41.74 | 8.61 |

Parameter | Unit | Measured | Predicted | Absolute error |
---|---|---|---|---|

Maximum rotational speed | ||||

Absolute value (r/min) | 1 | 468.66 | 449.01 | −19.65 |

2 | 463.59 | 448.87 | −14.72 | |

Rise ratio value (%) | 1 | 40.6 | 34.7 | −5.9 |

2 | 39.1 | 34.66 | −4.44 | |

Maximum pressure at spiral case inlet | ||||

Absolute value (mH_{2}O) | 1 | 464.35 | 439.54 | −24.81 |

2 | 465.32 | 439.52 | −25.8 | |

Rise ratio value (%) | 1 | 25.6 | 18.9 | −6.7 |

2 | 25.86 | 18.88 | −6.98 | |

Minimum pressure at draft tube inlet | ||||

Absolute value (mH_{2}O) | 1 | 32.03 | 39.59 | 7.56 |

2 | 33.13 | 41.74 | 8.61 |

Absolute error = predicted value − measured value; rise ratio = (predicted (or measured) value − initial value) / initial value.

The computational relative errors including the rotational speed and the pressures at the spiral case inlet for the two units were calculated from the time of initial load rejection to 26 seconds, as shown in Figure 9. All the errors were small at the beginning of load rejection and became a little larger in the later stage. Those were all within the allowable range. The main reasons were that the flow was extremely complex when the units were in the anti-pump operation condition and the Mockridge moment was probably not accurate enough. The computational errors of the pressures at the draft tube inlet for the two units were also calculated from the time of initial load rejection to 26 seconds, as shown in Figure 10. The errors were a little large, but most of them were within the allowable range. The main reasons were that the vortex rope in the draft tube affected the fluid flow severely and the flow was extremely unstable in the ‘S’ operation region for low specific speed pump-turbine, which caused large fluctuations in the measured data.

The average computational errors, including rotational speed, pressures at the spiral case inlet and at draft tube inlet for unit 1 were 4.5%, 4.9%, and 8.7 mH_{2}O, and the corresponding root mean square errors (RMSEs) were 5.3%, 5.9%, and 9.3 mH_{2}O. The average computational errors, including rotational speed, pressures at the spiral case inlet and at draft tube inlet for unit 2 were 3.1%, 4.7%, and 8.3 mH_{2}O, and the corresponding RMSEs were 3.9%, 5.9%, and 9.1 mH_{2}O. The comparisons showed that all the average errors were closer to the experimental values than the RMSE values.

## CONCLUSIONS

A new mathematical model based on MIC was proposed to simulate the transient process under the simultaneous load rejection of two units with joint closing law of ball-valve and guide vane in a pumped storage power station. Meanwhile, the hydraulic boundary conditions of pump-turbine, and other coupling boundary conditions of the ball-valve and the guide vanes were incorporated. The main conclusions were as follows.

In the first 10 seconds of the load rejection process, the pressure at the spiral case increased gradually to the maximum value. After a decreasing process of 10 seconds, the pressure began to increase and reached the second extreme. The trend of the pressure change at the draft tube inlet was almost opposite to the trend at the spiral case inlet. At the first stage (0–7 s) of the load rejection process, the rotational speed of the unit kept rising until it reached the first extremum. Then, the rotational speed decreased to the minimum in 13 seconds. After that, it started to increase to the second extremum and then gradually fluctuated to zero. The maximum rotational speed occurred about 5 seconds ahead of the moment at which the guide vane started to close, corresponding to nearly 3 seconds ahead of the moment when the maximum pressure occurred at the spiral case inlet.

At the spiral case inlet, the first pressure extremum subjected to the rotational speed increasing was greater than the second extremum, which was caused mainly by the guide vane closing process. It indicated that the centrifugal force of the fluid inside the runner had a great influence on the transient parameters in load rejection process of the pump-turbine with low specific speed.

The measured and predicted results confirmed that the maximum rotational speeds, the maximum pressures at the spiral case inlet, and the minimum pressures at the draft tube inlet all met the design specification requirements. During the simultaneous load rejection, the change trends of the transient parameters were similar for the two units which shared the same division pipeline. The technique for controlling flow discharge based on the joint closing law of ball-valve and guide vane for two units was acceptable.

The relative errors of the predicted maximum rotational speeds and the maximum pressures of the two units agreed well with the experimental data. The above analyses and comparisons indicated that the proposed models of MIC and numerical techniques were reliable in predicting the double load rejection process. The MIC can be used to analyze other transients of the hydraulic system with water turbine or pump as long as the signs of torque and head are changed properly.

## ACKNOWLEDGEMENT

This work was supported by the National Natural Science Foundation of China (Grant Nos 51779257, 51479196).

## REFERENCES

*.*

*,*

*.*