Design Optimization of Francis Runners For Sediment Handling
Design Optimization of Francis Runners For Sediment Handling
Design Optimization of Francis Runners For Sediment Handling
Mette Eltvik
The Waterpower Laboratory
Norwegian University of Science
and Technology
Norway
Ole G. Dahlhaug
The Waterpower Laboratory
Norwegian University of Science
and Technology
Norway
Bhola Thapa
Turbine Testing Lab
Kathmandu University
P.O. Box 6250
Dhulikhel, Nepal
Kristine Gjster
The Waterpower Laboratory
Norwegian University of Science
and Technology
Norway
Abstract
Growth of energy demands is escalating rapidly. Crisis in conventional sources of energy is urging to expand
feasibility limits of the renewable energy resources such as hydropower. New development of hydropower projects
is shifting towards the unexplored regions of Asia and South America. These regions have their own specific
technical challenges and one of the major problems is sediment erosion of turbine components. Financial feasibility
of future hydropower developments across these basins would be largely influenced by technological advancements.
New innovations to prevent erosion of mechanical equipment exposed to sediments are important needs at present.
Erosion in hydro turbines is a complex phenomenon, which depends upon several parameters. Design of Francis
turbines is unique to each site and hence takes time and effort to produce the best design for specific conditions. This
makes design optimization of Francis turbines for erosive environment a challenging task. Several studies have been
conducted to estimate erosion behavior in turbines and other components. Recent advancements in computing tools
and software have added advantage to these studies. However, the design of Francis turbines to handle large amount
of sediments effectively has still not been established.
Kathmandu University, together with Norwegian University of Science and Technology, DynaVec and Nepal Hydro
& Electric as partners, has started a unique research project. The goal of the project is to establish a new philosophy
for design of Francis turbines effective to handle sediment problems. Jhimruk Hydroelectric Center in Nepal, which
represents a typical case of sediment erosion in hydro turbines across the Himalayas, has been taken as the reference
case. The research is funded by NORAD supported RenewableNepal Programme. This paper assimilates some of
the major findings of the ongoing project.
A new program Khoj has been developed to create and optimize the design of Francis runners. The program is also
featured to compare erosion in runner blades for different design cases. The final design can be exported to
Computational Fluid Dynamics (CFD) and Fluid Structure Interaction (FSI) for further analysis. Parametric survey
was carried out with this program to evaluate the relative effect of each design parameter on sediment erosion. The
result was compared with results obtained using CFD analysis to estimate effects of the design variables on
hydraulic performance. Several optimized designs were developed and analyzed to fulfill the desired condition of
erosion and efficiency.
This paper summarizes the standard procedures used for design of high head Francis runners. Use of computational
tools and methods for design optimization of Francis runners for sediment handling is also presented. Application of
design program Khoj to identify the effects of design parameters on sediment erosion in Francis runners is
explained. Results of CFD analysis of new runner designs, optimized for handling sediments at the reference power
plant, are also discussed.
Key Words: Francis runner, design optimization, CFD, erosion, RenewableNepal
ASIA-2012, March 26-27, Int. J. Hydropower & Dams, Chiang Mai, Thailand
1. Introduction
There is a huge potential of new hydropower developments across the basins of Himalaya and Alps of Andes.
Nepal, a small south Asian country, alone has more than 6,000 rivers with 42,000 MW of feasible hydropower
potential still to be harnessed [1]. However problems due to suspended sediment particles have been one of the
technical limiting factors for designers and developers in this region. Almost all the river and rivulets in this region
contains 60% -80% of minerals having hardness number above 6 in Mohs scale [2].
Sediment erosion of hydro turbines is a complex phenomenon, as it depends upon various parameters [3, 4]. Several
studies have been conducted to quantify the effects of the erosion in hydraulic components of power plants [5, 6, 7].
These studies conclude that sediment erosion removes the base material gradually. This leads to change in flow
pattern, losses in efficiency, vibrations and final breakdown of hydro turbine components. This has made sediment
erosion of hydro turbines a technical problem with huge economic losses.
Several techniques has been attempted to control the effects of sediment erosion in turbine components. This
includes from prevention of sedimentation in the catchment areas to tapping sediments at intakes and applying
preventative coatings on the turbine components exposed to high velocity water [3]. However, these and others
methods [8] to prevent the sediment erosion in Nepalese power plants have not shown successful results. Change in
turbine design philosophy so as to reduce relative velocity of water inside the runner without significant changes in
turbine dimensions has been identified as one of the new areas of research in the field of sediment erosion in Francis
turbine [9].
(a)
(b)
(c)
(d)
Fig. 1 Francis Turbines performing under basins of Himalaya and Andes: (a) Damage in Runner blade at Jhimruk
PP after one year of operation [9]. (b) Runner manufactured with new technology for Cahua PP [10]. (c, d)
Damage to guide vanes and cheek plates at Nathpa PP [11].
ASIA-2012, March 26-27, Int. J. Hydropower & Dams, Chiang Mai, Thailand
These calculations are based on hydraulic parameters like head and discharge, which are determined by the
topography and hydrology of the power plant site. Traditionally, velocity at the inlet and outlet of the runner,
triangles as shown in Fig. 3, are used in the design process.
B
C
D
U
W
N
Zp
Cm
Cu
hva
hb
hs
Notations
Runner height, [m]
Absolute velocity, [m/s]
Runner diameter, [m]
Peripheral velocity, [m/s]
Relative velocity, [m/s]
Synchronous speed, [m/s]
No of poles in generator, [-]
Meridian component of C, [m/s]
Tangential component of C,
[m/s]
Vapor pressure, [m]
Barometric pressure, [m]
Submergence, [m]
Angle between the relative and
absolute velocity, [degree]
Inlet Section
Outlet Section
The outlet diameter and speed is found by reorganizing the expression for flow rate and peripheral speed,
respectively. Cm2 is obtained from the known geometry in the velocity triangles. The number of poles, Z, in the
generator depends on the rotational speed and net frequency. With a grid frequency of 50 Hz, the number of poles is
obtained by using equation 4.
[m]
[rpm] - (2)
- (1)
[m/s]
- (3)
[-]
- (4)
The number of poles should be an integer. The value obtained from equation 4 must therefore be rounded off. With
the corrected number of poles, equation 4 is used to find the corrected synchronous speed; this value is again used to
calculate the corrected diameter at the outlet.
To avoid the cavitation at the out let of the runner, the required submergence for the runner is calculated using
equation 5. The NPSH, called as Net Positive Suction Head, is calculated using empirical equation 6 [13].
[m]
- (5)
[m]
- (6)
- (7)
At the inlet the values for diameter, D1, height of the inlet, B1, and inlet angle, 1 is calculated. In order to find these
values the Euler equation (equation 7) is used. By introducing reduced dimensionless values and assuming no
rotation at the outlet, the equation can be rewritten as equation 8. The turbine efficiency, t, is generally set to 0.96,
and U1 is chosen in the interval as given in equation 9. U1 can be calculated using equation 10.
- (8)
- (9)
[m/s] - (10)
[m]
- (11)
ASIA-2012, March 26-27, Int. J. Hydropower & Dams, Chiang Mai, Thailand
The inlet diameter can now be found by using equation 11. From the velocity triangles in Fig. 3 an expression for
the inlet angle can be derived as given in equation 12. In order to use the equation 12 the value for Cm1 is calculated
by the continuity equation given in equation 13. It is desirable with an acceleration of the water from the inlet to the
outlet, in order to avoid backflow through the runner. Therefore, the condition presented in equation 14 is used.
Cm1
tan 1
Cm1 A1 Cm2 A2 - (13)
Cm1 1.1 Cm 2 - (14)
[-] - (12)
[m] - (15)
U1 Cu1
The only value left to calculate is the inlet height. By combining equation (13) and equation (14) this value can be
obtained as presented in equation 15.
3.2. Shaping of Runner Blade Profile
When the main dimensions of the runner are known, the runner blades can be designed. The design procedure starts
by determining the shape of the blade in the axial view, then the radial view is established, and nally the runner
blade can be plotted in three dimensions [13-15].
Runner Axial view: At rst a streamline along the shroud or along the hub has to be dened. It is most common to
dene it along the shroud. Traditionally an elliptical or circular shape is chosen for the streamline. Based on the inlet
distribution and the denition of the rst streamline, the rest of the streamlines can be determined. After the rst
streamline is dened, the number of streamlines has to be chosen. The distribution of streamlines is determined by
the flow rate, Q, at the inlet of the runner, which initially is considered uniform. This gives a uniform distribution of
the streamlines between the hub and shroud at the inlet. Fig. 4 shows the axial view of a Francis runner after
removing upward curvature at the hub, cutting endpoints and redistributing the remaining points.
Runner Radial view: In order to simplify the design process of going from the axial view to the radial view, a GHplane is dened, shown in Fig. 5. G is the length of a streamline in the axial plane and H is the length of a streamline
in the radial plane. The values of G are calculated by using the values of R and Z from the axial view. Calculating
the values of H is more demanding, as they are dependent on the distribution of the blade angle, . The blade angle
is closely linked to the energy distribution along the blade. The energy distribution describes the transformation
from pressure energy to rotational energy along the blade. The relation between the energy distribution and the blade
angle is governed by equation 16.
Height Z [m]
Radius R [m]
Fig. 6 GH-plane
The blade angles at the inlet and outlet are known from the velocity triangles. The distribution in between has to be
determined. This can be done in two ways. Either choosing the U Cu distribution and then calculating the
distribution, or choosing distribution and then calculate U Cu distribution. Choosing the blade angle distribution
gives the designer full control of the design outcome. The blade angle distribution is also the most important factor
to shape the runner profile.
The values of H corresponding to each G can be obtained using equation 17. When this is performed for each
streamline, the G-H plane can be plotted as shown in Fig. 6. Now the radial view of the runner can be established
using equation 18, which is dened according to Fig. 7. When knowing all the coordinates for and R, the radial
view can be plotted as in Fig. 8. Combining the axial coordinates and the radial coordinates, the 3D shape of the
runner blade emerges as shown in Fig. 9.
(
[radians]
- (16)
- (17)
[radians]
- (18)
ASIA-2012, March 26-27, Int. J. Hydropower & Dams, Chiang Mai, Thailand
x
y
unit
Runner efficiency
m /s
2.35
%
96
4.2. Design Optimization Range and Methodology
The hydraulic design parameters are varied within
Table 2. Hydraulic design parameters
a defined range and its effects on erosion factor is
Value for
Range of
evaluated. Table 2 lists the range of variation of S.N.
Parameters
Symbol Unit Reference
optimization
the design parameters considered for this study.
design
1
[m3/s3]
2
3
4
4
5
Outlet diameter
Number of pole
pairs in generator
Reduced
peripheral velocity
at inlet
Acceleration of
flow through
runner
Height of runner
Blade angle
distribution
D2
0.54
0.4 - 0.75
ZP
3 - 12
U1
0.74
0.65 - 1
Acc
35
0 - 50
0.16
degrees
linear
0.05 - 0.4
4 different
nonlinear
(19)
Where n is the number of segment area (Ai) in the runner blade surface. Wi is the relative velocity of flow in each
segment area. The segment area is the area between the intersection of streamlines and stream points in the runner
blade surface.
Erosion Factor (Ef)
It is ratio of erosion tendency of each new design with respect to the reference design. Erosion factor is defined as
follows:
ASIA-2012, March 26-27, Int. J. Hydropower & Dams, Chiang Mai, Thailand
[-]
(20)
The erosion factor estimates a quantitative difference in sediment erosion of runner with the change in hydraulic
design alone. In this study the erosion factor is used as a means to compare the relative erosion in the optimized
designs of runner with respect to the reference design. The erosion factor for the reference design is 1.
4.3. CFD Analysis
To verify the reference design, a CFD simulation is carried out. Designs from Matlab are exported to Ansys CFX13. Simulations are done to evaluate the hydraulic performance and erosion on the runner blade surface. Exactly
same process has been repeated to all the design analysis to maintain the consistency. Table 3-6 presents the
parameters selected for the CFD analysis. Fig 10 shows mesh generated from Ansys Turbo Grid for the reference
design for the parameters given in table 3. Fig. 11 shows the computational domain generated by Ansys CFX-Pre for
the parameters given in table 4.
Table 3 Parameters for CFX-Turbo Grid
Paramater
Type Value
Grid Node Count
Fine
250000
Factor Ratio
2
Reynolds No
500000
Table 5 Parameters for CFX-Pre Sediment Data
Data
Value
Unit
Material
Quartz
Density
2.65
g/cm3
Diamter
0.1
Mm
Shape factor
1
Flow rate
0.07
kg/s
ASIA-2012, March 26-27, Int. J. Hydropower & Dams, Chiang Mai, Thailand
rate density on the pressure side of reference runner blade computed by Ansys CFX-Solver for the sediment data
presented in table 5. It shows that the erosion pattern spreading at the outlet section of the runner blade, which
resembles to the erosion pattern of reference runner as shown in Fig. 1, a.
CFD Analysis of Optimized Designs
Several optimized designs are developed by combining the changes anticipated by the design program to reduce
sediment erosion in runner blades. CFD analysis of the optimized designs having erosion factor below 1 has been
done to validate the reduction in erosion and also evaluate the effects on other hydraulic parameters. Effects of the
variation on other hydraulic parameters have also been computed. Table 8 presents corresponding erosion factors of
the two optimized designs calculated by the design program together with respective total runner efficiencies
computed by CFD.
Table 8. Optimized designs erosion factors and efficiencies
Fig. 18 to Fig. 20 shows the results of CFD analysis of optimized design 2. The erosion on the blade surface is less
than the reference runner but not lower than that on the optimized design 1 as predicted by the program.
Irregularities in the streamlines have further increased and also the rough transitions of the relative velocity. These
could have caused to further decrease in the total runner efficiency.
ASIA-2012, March 26-27, Int. J. Hydropower & Dams, Chiang Mai, Thailand
6. Conclusions
Problem of sediment erosion in turbine components at the hydropower plants across the basins of the Himalayas and
the Andes is a technical challenge. Several cases of severe erosion in Francis turbines operating in these regions
have been reported. Some attempts to modify manufacturing processes to allow full coating in runner blades has
shown an optimistic progress in controlling erosion in Francis turbine. However, new designs of Francis turbines for
erosive environments should be explored to make the future of hydropower more economic and sustainable even in
the sediment laden projects.
This study makes an attempt to optimize the existing hydraulic design of high head Francis turbines, with the
possibility to minimize erosion. The fundamental equations followed for designing a Francis runner has been
explored. A reference case has been established and the optimizing parameters with the range of optimization have
been identified. Effects of design parameters on sediment erosion in runner blade surface have been estimated by a
new design program named as Khoj. Several optimized designs of the runner have been developed and CFD
analysis of selected ones has been performed.
From the parametric survey it is found that the runner outlet diameter, peripheral velocity at inlet, and blade angle
distribution have the highest effect on sediment erosion of Francis runner. Optimized designs developed show the
considerable reduction of erosion in runner blades. Findings of this study will help for further exploration of the
research to establish a new design methodology for reduced erosion and increased efficiency of Francis runners.
References
[1]
[2]
[3]
[4]
[5]
[6]
[7]
[8]
[9]
[10]
[11]
[12]
[13]
[14]
[15]
D.P.Sangroula, Hydropower Development and its Sustainability with respect to Sedimentation in Nepal, Journal of the
Institute of Engineering, Vol. 7, No. 1, pp. 1-9.
B. Thapa, R. Shrestha, P. Dhakal, B. S. Thapa, Problems of Nepalese Hydropower Projects due to Suspended Sediments,
Journal of Aquatic Ecosystem Health and Management, pp 251-258, 2005.
B. Thapa, Sand Erosion in Hydraulic Machinery, Doctoral thesis at NTNU, 2004.
A.A. Hamed, W. Tabakoff, R.B. Rivir, K. Das, P. Arora, Turbine Blade Surface Deterioration by Erosion, Journal of
Turbomachinery, 2005, Vol. 127, pp 445-452.
G.F. Truscott, A literature survey on abrasive wear in hydraulic machinery, Wear, 1971.
M. Padhy, R.P. Saini, A Review on silt erosion in hydro turbines, Renewable and Sustainable Energy Reviews, pp.
19741987, 2008.
Eltvik M., Sediment Erosion in Francis Turbines, Masters thesis at NTNU, 2009.
J. Ruud, Sediment handling problems Jhimruk Hydroelectric Center, Masters thesis at NTNU, 2004.
B.S. Thapa, B. Thapa, O.G. Dahlhaug, Center of Excellence at Kathmandu University for R&D and Test Certification of
Hydraulic Turbine, Proceeding of International Conference on Hydraulic Efficiency Measurement, India, 2010.
H. Brekke, Discussion of Pelton turbine versus Francis turbines for high head turbines, IAHR, Colorado, 1978.
O.G. Dahlhaug, P.E. Skre, V. Mossing, A. Gutierrez, Sediment resistive Francis runner at Cahua Power Plant, The
International Journal of Hydropower & Dams, Issue 2, pp109-112, 2010.
H.K. Sharma, Power generation in sediment laden rivers, The International Journal of Hydropower & Dams, Issue 6, pp
112-116, 2010.
H. Brekke, Hydraulic Turbines Design, Erection and Operation, NTNU Publication, 2000.
M. Eltvik, G. Olimstad, E.C. Walseth, High Pressure Hydraulic Machinery, NTNU publication, 2009.
K. Gjosater, Hydraulic design of Francis Turbine Exposed to Sediment Erosion Masters thesis at NTNU, 2011.
The Authors
Biraj Singh Thapa is a MS by Research student in Mechanical Engineering Department at Kathmandu University (KU) with the
topic Hydraulic Design of Francis Turbine to Minimize Sediment Erosion. He is also a researcher for RenewableNepal Project
at KU to develop a new design philosophy for Francis turbine to minimize losses due to sediment erosion.
Corresponding address: [email protected]
Mette Eltvik has graduated from the Norwegian University of Science and Technology (NTNU) in 2009. Her field of
specialization is sediment erosion in Francis turbines. After one year work as a research assistant in the Waterpower Laboratory
at NTNU, she is continuing her research as a PhD candidate in the same field. She is also a supporting researcher for the
RenewableNepal Project.
Corresponding address: [email protected]
ASIA-2012, March 26-27, Int. J. Hydropower & Dams, Chiang Mai, Thailand
Kristine Gjosater has obtained her masters degree on 2011 from NTNU with the topic Hydraulic Design of Francis Turbine
Exposed to Sediment Erosion. Currently she is working as a research assistant in the Waterpower Laboratory at NTNU, and is
also a supporting researcher for the RenewableNepal Project.
Corresponding address: [email protected]
Ole G. Dahlhaug obtained his PhD in Mechanical Engineering at NTNU. From 1992 to 1998 he worked at SINTEF as Research
Scientist, with research and testing of pumps and turbines. Currently he is Professor at the Waterpower Laboratory at NTNU. He
has been actively working with research on sediment erosion of turbine components and efficiency measurements of hydropower
plants. He is supervisor of the RenewableNepal project.
Corresponding address: [email protected]
Bhola Thapa has done Bachelor and Master in Mechanical Engineering from India. He has done PhD from NTNU in the topics
Sand Erosion of Hydraulic Machinery in 2004. He has been working at Kathmandu University, School of Engineering since
1994. He is a Professor of Mechanical Engineering and he has been serving as Dean of School of Engineering since 2005. He has
been involved in research in the field of hydraulic turbines especially in the area design of erosion resistance, operation and
maintenance.
Corresponding address: [email protected]