Prediction of a Francis turbine prototype full load instability from investigations on the reduced scale model

The growing development of renewable energies combined with the process of privatization, lead to a change of economical energy market strategies. Instantaneous pricings of electricity as a function of demand or predictions, induces profitable peak productions which are mainly covered by hydroelectric power plants. Therefore, operators harness more hydroelectric facilities at full load operating conditions. However, the Francis Turbine features an axi-symmetric rope leaving the runner which may act under certain conditions as an internal energy source leading to instability. Undesired power and pressure fluctuations are induced which may limit the maximum available power output. BC Hydro experiences such constraints in a hydroelectric power plant consisting of four 435 MW Francis Turbine generating units, which is located in Canada's province of British Columbia. Under specific full load operating conditions, one unit experiences power and pressure fluctuations at 0.46 Hz. The aim of the paper is to present a methodology allowing prediction of this prototype's instability frequency from investigations on the reduced scale model. A new hydro acoustic vortex rope model has been developed in SIMSEN software, taking into account the energy dissipation due to the thermodynamic exchange between the gas and the surrounding liquid. A combination of measurements, CFD simulations and computation of eigenmodes of the reduced scale model installed on test rig, allows the accurate calibration of the vortex rope model parameters at the model scale. Then, transposition of parameters to the prototype according to similitude laws is applied and stability analysis of the power plant is performed. The eigenfrequency of 0.39 Hz related to the first eigenmode of the power plant is determined to be unstable. Predicted frequency of the full load power and pressure fluctuations at the unit unstable operating point is found to be in general agreement with the prototype measurements.

[1]  Jean Eustache Prenat,et al.  Surging of 140 MW Francis turbines at high load, analysis and solution , 1992 .

[2]  Yoshinobu Tsujimoto,et al.  A Theoretical Analysis of Rotating Cavitation in Inducers , 1993 .

[3]  G. G. Stokes "J." , 1890, The New Yale Book of Quotations.

[4]  Jason Foust,et al.  Overload Surge Investigation Using CFD Data , 2009 .

[5]  Jean Eustache Prenat,et al.  Comportement dynamique d'une turbine Francis à forte charge. Comparaisons modèle-prototype , 1988 .

[6]  C. Brennen,et al.  Dynamics of a Cavitating Propeller in a Water Tunnel , 2003 .

[7]  Christopher E. Brennen,et al.  The Dynamic Transfer Function for a Cavitating Inducer , 1976 .

[8]  F. Menter Two-equation eddy-viscosity turbulence models for engineering applications , 1994 .

[9]  Yoshinobu Tsujimoto,et al.  One-Dimensional Analysis of Full Load Draft Tube Surge , 2007 .

[10]  Christophe Nicolet,et al.  Overload Surge Event in a Pumped-Storage Power Plant , 2006 .

[11]  Christophe Nicolet,et al.  Influence of the vortex rope location of a Francis Turbine on the hydraulic system stability , 2008 .

[12]  Christophe Nicolet,et al.  Influence of the Francis turbine location under vortex rope excitation on the hydraulic system stability , 2009 .

[13]  Giuseppe Pezzinga,et al.  Second viscosity in transient cavitating pipe flows , 2003 .

[14]  Christophe Nicolet,et al.  Numerical Simulation of Nonlinear Self Oscillations of a Full Load Vortex Rope , 2009 .

[15]  Christopher E. Brennen Bubbly flow model for the dynamic characteristics of cavitating pumps , 1978 .

[16]  Mark E. Duttweiler,et al.  Surge instability on a cavitating propeller , 2001, Journal of Fluid Mechanics.