Research Papers: Fundamental Issues and Canonical Flows

Measurement of the Self-Oscillating Vortex Rope Dynamics for Hydroacoustic Stability Analysis

[+] Author and Article Information
Andres Müller

EPFL Laboratory for Hydraulic Machines
Lausanne 1007, Switzerland
e-mail: andres.mueller@epfl.ch

Keita Yamamoto

EPFL Laboratory for Hydraulic Machines
Lausanne 1007, Switzerland
e-mail: keita.yamamoto@epfl.ch

Sébastien Alligné

Power Vision Engineering LLC,
Ecublens 1024, Switzerland

Koichi Yonezawa

Assistant Professor
Graduate School of Engineering Science,
Osaka University,
Toyonaka, Osaka 560-8531, Japan

Yoshinobu Tsujimoto

Professor Emeritus
Graduate School of Engineering Science,
Osaka University,
Toyonaka, Osaka 560-8531, Japan

François Avellan

EPFL Laboratory for Hydraulic Machines,
Lausanne 1007, Switzerland

1Corresponding author.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received December 2, 2014; final manuscript received October 8, 2015; published online November 4, 2015. Assoc. Editor: Olivier Coutier-Delgosha.

J. Fluids Eng 138(2), 021206 (Nov 04, 2015) (8 pages) Paper No: FE-14-1721; doi: 10.1115/1.4031778 History: Received December 02, 2014; Revised October 08, 2015

Flow instabilities in hydraulic machines often feature oscillating cavitation volumes, which locally introduce compliance and mass flow gain effects. These unsteady characteristics play a crucial role in one-dimensional stability models and can be determined through the definition of transfer functions for the state variables, where the cavitation volume is commonly estimated from the discharge difference between two points located upstream and downstream of the cavity. This approach is demonstrated on a test rig with a microturbine, featuring a self-oscillating vortex rope in its conical draft tube. The fluctuating discharges at the turbine inlet and the draft tube outlet are determined with the pressure–time method using differential pressure transducers. The cavitation volume is then calculated by integrating the corresponding discharge difference over time. In order to validate the results, an alternative volume approximation method is presented, based on the image processing of a high-speed flow visualization. In this procedure, the edges of the vortex rope are detected to calculate the local cross section areas of the cavity. It is shown that the cavitation volumes obtained by the two methods are in good agreement. Thus, the fluctuating part of the cavitation volume oscillation can be accurately estimated by integrating the difference between the volumetric upstream and downstream discharges. Finally, the volume and discharge fluctuations from the pressure–time method are averaged over one mean period of the pressure oscillation. This enables an analysis of the key physical flow parameters’ behavior over one characteristic period of the instability and a discussion of its sustaining mechanisms.

Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.


Susan-Resiga, R. , Ciocan, G. , Anton, I. , and Avellan, F. , 2006, “ Analysis of the Swirling Flow Downstream a Francis Turbine Runner,” ASME J. Fluids Eng., 128(1), pp. 177–189. [CrossRef]
Rheingans, W. J. , 1940, “ Power Swings in Hydroelectric Power Plants,” Trans. ASME, 62, pp. 171–184.
Dörfler, P. , 1985, “ Francis Turbine Surge Prediction and Prevention,” Waterpower 85, ASCE, pp. 952–961.
Jacob, T. , Prenat, J. , and Maria, D. , 1988, “ Dynamic Behavior at High Load of a Francis Water Turbine. Model/Prototype Comparison,” Houille Blanche, 3–4, pp. 293–300. [CrossRef]
Paynter, H. M. , 1953, “ Electrical Analogies and Electronic Computers: A Symposium: Surge and Water Hammer Problems,” Trans. Am. Soc. Civil Eng., 118(1), pp. 962–989.
Nicolet, C. , 2007, “ Hydroacoustic\ing and Numerical Simulation of Unsteady Operation of Hydroelectric Systems,” Ph.D. thesis, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland.
Dörfler, P. , 2009, “ Evaluating 1D Models for Vortex-Induced Pulsation in Francis Turbines,” 3rd Meeting of the IAHR Workgroup on Cavitation and Dynamic Problems in Hydraulic Machinery and Systems, Brno, Czech Republic, pp. 14–16.
Alligné, S. , 2011, “ Forced and Self Oscillations of Hydraulic Systems Induced by Cavitation Vortex Rope of Francis Turbines,” Ph.D. thesis, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland.
Landry, C. , Nicolet, C. , Bergant, A. , Müller, A. , and Avellan, F. , 2012, “ Modeling of Unsteady Friction and Viscoelastic Damping in Piping Systems,” IOP Conf. Ser. Earth Environ. Sci., 15(5), p. 052030.
Alligne, S. , Nicolet, C. , Tsujimoto, Y. , and Avellan, F. , 2014, “ Cavitation Surge Modelling in Francis Turbine Draft Tube,” J. Hydraul. Res., 52(3), pp. 399–411. [CrossRef]
Landry, C. , Favrel, A. , Müller, A. , Nicolet, C. , and Avellan, F. , 2015, “ Local Wave Speed and Bulk Viscosity in Francis Turbines at Part Load Operation,” J. Hydraul. Res. (in press).
Koutnik, J. , Nicolet, C. , Schohl, G. A. , and Avellan, F. , 2006, “ Overload Surge Event in a Pumped Storage Power Plant,” 23th IAHR Symposium on Hydraulic Machinery and Systems, Yokohama, Japan.
Chen, C. , Nicolet, C. , Yonezawa, K. , Farhat, M. , Avellan, F. , and Tsujimoto, Y. , 2008, “ One-Dimensional Analysis of Full Load Draft Tube Surge,” ASME J. Fluids Eng., 130(4), p. 041106. [CrossRef]
Brennen, C. , and Acosta, A. , 1976, “ Dynamic Transfer Function for a Cavitating Inducer,” ASME J. Fluids Eng., 98 Ser 1(2), pp. 182–191. [CrossRef]
Braisted, D. , and Brennen, C. , 1980, “ Auto-Oscillation of Cavitating Inducers,” Polyphase Flow and Transport Technology, ASME, pp. 157–166.
Tsujimoto, Y. , Kamijo, K. , and Yoshida, Y. , 1993, “ A Theoretical Analysis of Rotating Cavitation in Inducers,” ASME J. Fluids Eng., 115(1), pp. 135–141. [CrossRef]
Koutnik, J. , and Pulpitel, L. , 1996, “ Modeling of the Francis Turbine Full-Load Surge,” Modeling, Testing and Monitoring for Hydro Power Plants, Lausanne, Switzerland.
Philibert, R. , and Couston, M. , 1998, “ Francis Turbines at Part Load, Matrix Simulating the Gaseous Rope,” 19th IAHR Symposium on Hydraulic Machinery and Cavitation, Singapore, Sept. 9–11.
Yonezawa, K. , Konishi, D. , Miyagawa, K. , Avellan, F. , Dörfler, P. , and Tsujimoto, Y. , 2012, “ Cavitation Surge in a Small Model Test Facility Simulating a Hydraulic Power Plant,” Int. J. Fluid Mach. Syst., 5(4), pp. 152–160. [CrossRef]
Dörfler, P. , and Ruchonnet, N. , 2012, “ A Statistical Method for Draft Tube Pressure Pulsation Analysis,” Earth Environ. Sci., 15(6), p. 062002.
Gibson, N. R. , 1923, The Gibson Method and Apparatus for Measuring the Flow of Water in Closed Conduits, ASME, pp. 343–392.
Kashima, A. , Lee, P. J. , Ghidaoui, M. S. , and Davidson, M. , 2013, “ Experimental Verification of the Kinetic Differential Pressure Method for Flow Measurements,” J. Hydraul. Res., 51(6), pp. 634–644. [CrossRef]
Yamamoto, K. , Müller, A. , Ashida, T. , Yonezawa, K. , Avellan, F. , and Tsujimoto, Y. , “ Experimental Method for the Evaluation of the Dynamic Transfer Matrix Using Pressure Transducers,” J. Hydraul. Res., 53(4), pp. 466–477. [CrossRef]
Müller, A. , Dreyer, M. , Andreini, N. , and Avellan, F. , 2013, “ Draft Tube Discharge Fluctuation During Self-Sustained Pressure Surge: Fluorescent Particle Image Velocimetry in Two-Phase Flow,” Exp. Fluids, 54(4), pp. 1–11. [CrossRef]
Stoer, J. , Bulirsch, R. , Bartels, R. , Gautschi, W. , and Witzgall, C. , 2002, “ Introduction to Numerical Analysis,” Introduction to Numerical Analysis, (Texts in Applied Mathematics), Springer, New York.
Otsu, N. , 1979, “ Threshold Selection Method From Gray-Level Histograms,” IEEE Trans. Syst. Man Cybern., 9(1), pp. 62–66. [CrossRef]
Bendat, J. S. , and Piersol, A. G. , 2010, Random Data: Analysis and Measurement Procedures, 4th ed., Wiley, Hoboken, NJ.
Müller, A. , Favrel, A. , Landry, C. , Yamamoto, K. , and Avellan, F. , 2014, “ On the Physical Mechanisms Governing Self-Excited Pressure Surge in Francis Turbines,” Earth Environ. Sci., 22(3), p. 032034.
Müller, A. , 2014, “ Physical Mechanisms Governing Self-Excited Pressure Oscillations in Francis Turbines,” Ph.D. thesis, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland.
Alligné, S. , Nicolet, C. , Ruchonnet, N. , Hasmatuchi, V. , Maruzewski, P. , and Avellan, F. , 2009, “ Numerical Simulation of Nonlinear Self Oscillations of a Full Load Vortex Rope,” Vol. 2, 3rd Meeting of the IAHR Workgroup on Cavitation and Dynamic Problems in Hydraulic Machinery and Systems, Brno, Czech Republic, pp. 325–338.


Grahic Jump Location
Fig. 1

Closed-loop test configuration for the simulation of a hydraulic powerplant with microturbine and conical diffuser

Grahic Jump Location
Fig. 2

Microturbine runner (top) and volute geometry (bottom)

Grahic Jump Location
Fig. 3

Close view of the draft tube cone with an inlet and outlet diameter of D1¯ = 31mm and cone Dcone = 57 mm. The straight part measures L1 = 50.8 mm and the conical part has a length of L2 = 186 mm with a half-angle of 4 deg.

Grahic Jump Location
Fig. 4

Draft tube flow visualization showing one period T of the vortex rope oscillation: (a) t = 0, (b) t = 0.11 · T, (c) t = 0.22 · T, (d) t = 0.33 · T, (e) t = 0.44 · T, (f) t = 0.55  ·  T, (g) t = 0.66 · T, (h) t = 0.77 · T, (i) t = 0.88 · T, and (j) t = T

Grahic Jump Location
Fig. 5

Amplitude (solid line) and phase angle (dashed line) of the cross-spectral power density between the pressure sensors p3 and p6

Grahic Jump Location
Fig. 6

Control volume in the draft tube containing the volume Vc

Grahic Jump Location
Fig. 7

Draft tube cone in cavitation-free conditions with sections I, II, III, IV, and V for image processing

Grahic Jump Location
Fig. 8

Subtraction of a reference image from a video frame (top) and masked binary conversion (bottom) for the edge detection in sections I, III, and V

Grahic Jump Location
Fig. 9

Pixel value curve along a vertical line for the edge detection in sections II and IV. The dashed horizontal lines represent the locations of the estimated cavity edges according to the local minima of the pixel value curve.

Grahic Jump Location
Fig. 10 Top

low-pass filtered cp from p2 with instant phase (dashed line) and bottom: mean phase average with standard deviation (vertical error bars)

Grahic Jump Location
Fig. 11

Nondimensional volume Vc* estimated from flow visualization (point markers) with low-pass filtered signal (solid line) as a function of the number of runner revolutions

Grahic Jump Location
Fig. 12

Nondimensional vortex rope volume Vc* from flow visualization, wall-pressure coefficient from p2, and nondimensional discharge fluctuations Q̃1∗, and Q̃2∗ as a function of the number of runner revolutions

Grahic Jump Location
Fig. 13

Nondimensional, centered Vc* from flow visualization (point markers) and from discharge measurements (dashed line) as a function of the number of runner revolutions

Grahic Jump Location
Fig. 14

Mean pressure phase-averaged vortex rope volume (top) and discharge fluctuations (bottom) with standard deviations (vertical error bars)

Grahic Jump Location
Fig. 15

Vortex rope volume as a function of the upstream and downstream discharge fluctuations during one mean period of the pressure oscillation




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In