Strojniski vestnik - Journal of Mechanical Engineering 56(2010)6, 377-384 UDC 629.4.027.31-272.2:531.2:004.8 Paper received: 10.03.2009 Paper accepted: 16.04.2010 Simulation of Transient Flow in Hydroelectric Power Plants Using Unsteady Friction Alireza Riasi* - Mehrdad Raisee - Ahmad Nourbakhsh University of Tehran, Department of Mechanical Engineering, Iran In this paper the influence of surge tanks on transient flow generated by changing turbine gate opening has been numerically studied. For this purpose, the governing equations of transient flow are solved using the method of characteristics and considering unsteady friction. The surge tank and the hydraulic turbine are considered as internal boundary conditions. Comparison of transient results for pressure oscillations with the available experimental data shows good agreement. The surge tanks lower down the pressure levels at the upstream side of turbine and stabilize the pressure at the downstream side. Also, surge tanks delay the time occurrence of the maximum and minimum pressure. Finally, the effect of wicket gate closing followed by wicket gate opening on the so called speed no load operation (SNLO) on the transient flow is investigated and transient conditions improved by this sequence are found. ©2010 Journal of Mechanical Engineering. All rights reserved. Key words: hydraulic turbine, hyperbolic PDE, surge tank, transient flow, unsteady friction 0 INTRODUCTION Unsteady flow analysis in water power stations is one of the most significant issues in the prediction of undesirable pressure variations in waterways as well as probable changes in rotor speed for the power plants safe operation. Hydraulic transient in a hydroelectric power plant is usually caused by various turbine operations such as start up, load acceptance and load rejection. Changing conditions in the power market have led to an increase in the demand of peak energy generation, short response time and fast frequency changes. Consequently, there is a need for accurate simulation of transient flow in hydroelectric power plants. Traditionally, transient flow in a water power station is simulated using well-known water hammer equations and considering the steady state skin friction. However, it was found that the use of steady friction approximation does not produce sufficient damping comparison to the measured data [1] to [3]. Therefore, unsteady friction models must be adopted to account for this damping deficit. Surge tanks are installed at the upstream or downstream sides of hydropower plant waterways. These are of the main components that equilibrate pressures and thus need a careful optimization. The main hydraulic requirements to be met by the designers of surge tanks include damping of the water hammer effect as well as accelerating and decelerating tunnel flows, improving the response capabilities of the mechanical equipment and sometimes used as a gate shaft at the same time. Wicket gate opening at speed no load operation (SNLO) is the wicket gate opening at which unit (turbine and generator) rotate at synchronous speed with zero output. At SNLO, the turbine output is equal to the turbo generator windage and friction losses at synchronous speed. When the load rejection phenomenon occurs in a hydropower plant, wicket gates are completely closed and in some conditions the wicket gates are reopened from complete closing to SNLO opening under the governor control. In this new situation the unit operates in the standby mode. In the present paper, this type of transient is called a wicket gate closing followed by an opening to SNLO. The main objectives of this paper are to simulate the transient flow in hydroelectric power stations using the method of characteristics and unsteady friction model and to investigate the effects of the surge tank and wicket gate closing followed by an opening to SNLO on transient conditions. To fulfill these aims, the governing equations including continuity and momentum equations are numerically solved along the water ways. The upstream and downstream reservoirs * University of Tehran, Department of Mechanical Engineering, Tehran, P.O. Box:11155-4563, Iran, ariasi@ut.ac.ir are considered as constant head boundary conditions. Also, surge tank and Francis hydraulic turbine in waterway are considered as internal boundary conditions. To verify the numerical method, the pressure oscillations at upstream and downstream of the turbine formed by the load rejection are compared with the measured data of a real test case. Finally, the influence of surge tanks and wicket gate closing followed by an opening to SNLO on transient conditions has been investigated. 1 MATHEMATICAL MODEL AND NUMERICAL PROCEDURE Unsteady closed conduit flow equations consist of one-dimensional, hyperbolic partial differential equations, in which the continuity and momentum equations are, [4] and [5]: ÔH ÔH_ Ôt Ôx ÔV +g ÔV=0 g Ôx —+VÔV + ÔH Ôt + Ôx g Ôx (1) pD Tw = 0 where H(x,t) is piezometric head and V(x,t) is average axial velocity. The wall shear stress is the sum of the steady and unsteady wall shear stresses as T — T + T (2) w ws wu in which the steady state wall shear stress is defined as — pf QQ T <2 8 (3) where Q is flow rate; A is cross section area of the conduit and f is Darcy-Weisbach friction factor. The above steady state friction formula (Darcy-Weisbach) fails to reproduce the pressure variation held beyond the first wave cycle [1] to [3]. To remedy this failure, the modified model of Brunone et al. [6] is implemented in the method of characteristics to solve the governing equations. Bergant et al. [7] and Pezzinga [8] found that the model of Brunone et al. [6] fails to correctly represent the unsteady friction for some transient event types and thus the model modifications are proposed. The modified model of Pezzinga [8] is expressed as: t = kPR(ÔQ + sign(QÔQ)aÔQ T 4A l Ôt + Sign(Q Ôx )a Ôx (4) where k is empirical coefficient and sign(x) is sign function. The term sign (Q-dQ/dx) is added to the original Brunone model by Pezzinga in order to make the model applicable for both upstream and downstream transient flows. For a power station with upstream and downstream waterways and also considerable length, this term is very important. This is because of the transient flow behaviour in upstream and downstream of the hydraulic turbine. The method of characteristics (MOC) is a common method for solving the water hammer equation. The convective term in the continuity and momentum equations can be neglected due to a small Mach number (V/a«1). With linear combination of water hammer equations, the following systems of compatibility equations for positive and negative characteristics lines are given [9]: C+ : dQ+gAdH +4At = 0 dt a dt pD w : dQ gAdH+4A T = 0 ' dt ' a dt pD w ' (5a) (5b) With reference to the grid shown in Fig.1, the discretization of Eqs. (5a) and (5b) along the positive and negative characteristics gives: (Qn+1 - Qi"-1) + C(Hn+1 - HI,) + 4AAt = 0 pD w Q+1 - Q+i) - C ( Hr1 - Hi+i) + 4 AAt pD < = 0 (6a) (6b) where At is time step; C is gA/a; superscript n is temporal location and subscript i is spatial location. t- n+1 At Ax=aAt Ax=aAt Fig. 1. x-t plane showing characteristics for Eq. (6) 1.1 Boundary Conditions Constant head boundary condition is applied for the upstream and the downstream reservoir. To simulate Francis hydraulic turbine as internal boundary condition, the relation between the net head and discharge has to be specified. Fluid flow through a Francis turbine depends upon various parameters such as net head, rotational speed of the unit and the wicket gate opening. Curves representing the relationship between these parameters are called turbine characteristics (Hill Chart). No data has been reported in the literature for turbine characteristics during transient conditions. Therefore, the steady state model test results (Hill Chart) are used to plot the experimental prototype turbine characteristics and these are assumed to be valid during the transient state as well. Also, very little data are available for small wicket gate openings and to cover the range, the characteristic curves are extrapolated. Furthermore, the positive and negative characteristics equations are valid for upstream and downstream of the turbine. The mentioned equations with the aid of continuity equating for flow before and after the turbine are solved simultaneously for the hydraulic turbine as an internal boundary condition [4], [5] and [10]. In contrast to reaction turbines, the impulse type turbine dynamic equations and water hammer equations are solved separately [11]. Also, the surge tank is an internal boundary condition which is fully described in the next section. 2 SURGE TANK One of the components to control the transients is the surge tank. This is an open chamber or a tank connected to the pipeline which reflects the pressure waves by supplying or storing excess liquid. As shown in Fig. 2, in this study a surge tank with a standpipe is considered for numerical simulation. The positive and negative characteristic equations are valid for nodes i and i+1, respectively Eqs. (6a) and (6b). The continuity equation is: Fig. 2. Notations for a surge tank with stand pipe where Qsp is the flow rate in the standpipe (flow in the upward direction is considered positive). The pressure losses between nodes i and i+1 can be neglected, thus Hn +1 _Hn +1 i i +1 ' (8) Since the length of the standpipe is usually short compared to the length of the pipe of the system, the liquid inside the standpipe may be considered as a lumped mass. Considering that the weight of the liquid in the standpipe is, pgAspLsp. PgAspfLsp Qsp \Qsp /(2 gA^L differential equation may be written for the acceleration of the flow in the standpipe [4] and the friction force is, ), the following L dQ sp Z-'SP gA dt = H+1 — fLP 2 gA2Ds -Qn Q n - Z" ZL-sp I s (9) where 2S is the height of the water surface in the tank above the datum. The elevation of the liquid surface in the tank above the datum line at the beginning and the end of the time step can be related to each other with the Eq. (10): +1 = z„ + _Ai_ (q„+1 s s 2A sp +q:p ). (10) Eqs. (6a, b) and (7) to (10) should be solved at the end of each time step for the six unknowns H?+1, qr1 = q« + q: (7) h:+ Qr, Q«1, QT and zj"+1. 3 CASE STUDY AND NUMERICAL RESULTS In the following, an application of the developed computer program for a prediction of the transient flow behaviour in Karun3 hydropower plant is presented. This power station is a 2000 MW (8 x 250 MW) Francis type hydropower plant, located in Khuzestan province on the Karun river in the South of Iran. Four surge tanks are located on the upstream water ways and eight surge tanks are located on the down stream water ways. The main characteristics of this hydropower plant are listed in Table 1. Table 1. The main characteristics of Karun 3 hydro power plant [12] Parameter Value Rated power [MW] 255.0 Rated head [m] 161.0 Rated flow [m3/s] 171.3 Rotational speed [rpm] 187.5 Head water level [m] 820.0 Tail water level [m] 660.0 Moment of inertia [ton.m2] 29500 Turbine elevation [m] 653.1 The profile of this power station is shown in Fig. 3 and its schematic for the four units is shown in Fig. 4. The average acoustic wave velocity is about 1100 [m/s]. Fig. 3. Power station profile [12] Fig. 4. Schematic of the waterways [12] Transient flow due to a sudden closing of the guide vanes (generator load rejection) has been considered in this case study for the two units fed with a main penstock (units 5 and 6). The wicket gates are closed in 20 seconds. The pressure oscillations were measured by pressure transducers attached to the spiral case and to the draft tube. In Fig. 5 pressure oscillations at the spiral case and the draft tube are compared with the measured data. The numerical results are in good agreement with the measured data. Moreover, the maximum and minimum pressure and also the pressure amplitude and the pressure phase of the transient wave are well computed by numerical method. Slight differences between the computed and the measured data are due to uncertainties in 200 ¥ 190 "3 S 150 Measured Computed 40 60 Time [s] the estimation of friction factor, waterway diameter, waterway length and especially the wave speed. The wave propagation speed varies along the concrete and steel lined tunnel and depends on the elastic properties of the fluid and pipe material, as well as on geometry and the pipe anchorage [13]. Therefore, the value of the wave speed is not precisely known due to errors in the estimation of the mentioned parameters. The path which is passed by Francis turbine during wicket gate closing is shown in the approximated characteristics curves (Hill chart) of this hydropower plant (Fig. 6). As shown in Figure 6, during wicket gate closing the discharge approaches zero whereas the net head slightly increases. [m I 10 & 1 I 0