Next Article in Journal
Suitability of Remediated PFAS-Affected Soil in Cement Pastes and Mortars
Next Article in Special Issue
Effect of Residents’ Involvement with Small Hydropower Projects on Environmental Awareness
Previous Article in Journal
Alternative Metrics for Assessing the Social Impact of Tourism Research
Previous Article in Special Issue
Probabilistic Reliability Enhancement Strategies of Hydro Dominant Power Systems under Energy Uncertainty
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Study of a Francis Turbine over Wide Operating Range: Some Practical Aspects of Verification

Waterpower Laboratory, NTNU—Norwegian University of Science and Technology, 7491 Trondheim, Norway
*
Author to whom correspondence should be addressed.
Current address: Design Engineer, Rainpower Norge AS, Norway.
Sustainability 2020, 12(10), 4301; https://doi.org/10.3390/su12104301
Submission received: 7 January 2020 / Revised: 2 May 2020 / Accepted: 22 May 2020 / Published: 25 May 2020

Abstract

:
Hydropower plays an essential role in maintaining energy flexibility. Modern designs focus on sustainability and robustness using different numerical tools. Automatic optimization of the turbines is widely used, including low, mini and micro head turbines. The numerical techniques are not always foolproof in the absence of experimental data, and hence accurate verification is a key component of automatic optimization processes. This work aims to investigate the newly designed Francis runner for flexible operation. Unsteady simulations at 80 operating points of the turbine were conducted. The numerical model consisted of 16 million nodes of hexahedral mesh. A SAS-SST (scale adaptive simulation-shear stress transport) model was enabled for resolving/modeling the turbulent flow. The selected time-step size was equivalent to one-degree angular rotation of the runner. Global parameters, such as efficiency, torque, head and flow rate were considered for proper verification and validation. (1) A complete hill diagram of the turbine was prepared and verified with the reference case. (2) The relative error in hydraulic efficiency was computed and the over trend was studied. This allowed us to investigate the consistency of the numerical model under extreme operating conditions, far away from the best efficiency point. (3) Unsteady fluctuations of runner output torque were studied to identify unstable regions and magnitude of torque oscillations.

1. Introduction

Flexibility of power generation is an important requirement for future turbine designs. The turbines need to operate outside the guaranteed region with more start–stop cycles and high-ramping rates. However, power generation at off-design load brings specific challenges, such as high-amplitude pressure pulsations, vortex breakdown, cavitation etc. Present day variable-speed operation of a hydro turbine is considered as an alternative option to meet fluctuating energy demand, as it allows a high-ramping rate. New design approaches are focusing on design optimization for the variable-speed operation including high, medium, low, mini and micro head turbines [1,2,3]. Furthermore, mini and micro hydro stations are good options for smart grid configurations. This will enable local flexibility of energy need and sustainability. One of the challenges for the very low head turbines is the lower efficiency. However, modern design approach and automatic optimization with proper verification will certainly be useful for improving the efficiency of very low head turbines.
Modern algorithms and tools allow designers to perform optimization with reduced effort. Numerical methods have played significant roles in design and development of hydro turbines over the last 50 years [4]. Presently, almost all turbine designs and the optimization are highly dependent on the numerical outcome. High quality results and the accurate prediction of physical phenomena are always a valuable input for the turbine designers. During the past few years, several verification approaches have been developed to quantify the numerical errors [5,6,7]. However, these approaches have certain limitations and require a customized solution for hydro turbines. The present work is the continuation of previous efforts initiated to study the practical aspects of accurate verification and validation of turbine simulations. During 2009–2014, the focus was the Francis-99 test case [8,9]. The reviewed parameters were passage modeling, component modeling (runner and draft tube separately), influence of boundary conditions, discretization techniques, turbulence models and time step sizes. Later (2014–2016), the work extended to other hydro turbines, and the data of several turbines were compiled [10]. The analysis underlined that the hybrid models with newly developed techniques of passage modeling (blade row) are able to predict turbine flow with reasonable accuracy. The gain in simulation speed is high; however, the accuracy related to flow periodicity, variable torque from each blade channel and flow asymmetry at runner outlet is expected to worsen, due to the assumption of periodical boundary type. This may be a compromise solution for the industry standard modeling of hydro turbines. The work also covered the fluid structure interactions and related numerical errors in hydro turbines [11,12]. Recently [13,14,15], a comprehensive study (including the compressibility effect) was conducted to develop a systematic verification and validation guide for the turbines. The conclusion was that more comprehensive study is needed to determine the precise trend and fit (systematic/random behavior) of numerical errors. In this paper, we focus on the complete hill diagram of a Francis turbine that allows credible quantification of numerical errors at several operating points and helps to investigate the systematic or random behavior of errors.

2. Test Case

A newly designed and optimized Francis runner was used for the verification and validation study, hereafter referred to as a HydroFlex runner. The specific speed ( N Q E —see IEC 60193:1999-11, sub clause 1.3.3.12.11) of the turbine is 0.027. The runner is designed for flexible operations. Other components of the turbine, except for the runner, are the same as the Francis-99 runner [9]. The modeled turbine includes spiral casing with 14 stay vanes, 28 guide vanes, a runner with 17 blades and the elbow type draft tube. Detailed description of the test facility is presented in a previous publication [16]. Table 1 describes the boundary conditions and the numerical parameters enabled for the simulations. Figure 1 shows the geometry of the runner and the hexahedral mesh. The complete turbine includes around 16 million nodes (around four node points per cubic millimeter volume). Unsteady simulations at operating points (a total of 80 flow rate and speed combinations) were conducted.

3. Verification and Validation

The HydroFlex runner has geometrical parameters identical to the Francis-99 runner, e.g., the runner inlet and outlet diameters and the profiles of the crown and band seals are same. The other geometrical dimensions, such as blade, crown and band profiles, are based on optimization. Figure 2 shows the iso-efficiency hill diagram of the Francis-99 runner. The iso-efficiency hill diagram of a turbine represents the overall performance of the turbine for a wide range of operating conditions. Figure 3 shows the iso-efficiency hill diagram of the newly designed HydroFlex runner. The hill diagrams include guide vanes angular positions, and hydraulic efficiency on the axis of speed factor, n E D (indirect representation of turbine rotational speed) and discharge factor, Q E D (indirect representation of flow rate through the runner). The hill diagrams do not include runaway and no-load discharge characteristics. The hill diagram of the Francis-99 runner was prepared using available data of model acceptance tests in the laboratory. The best efficiency point is located at n E D = 0.18 ,   Q E D = 0.153 and α = 100 % , and the hydraulic efficiency is 93.52 ± 0.16 % . The operating range for the Francis-99 runner is n E D = 0.15 0.22 and Q E D = 0.05 0.22 . The experimental data of the model acceptance test of the Francis-99 runner in the laboratory were used to prescribe the boundary conditions for the HydroFlex runner. The numerical model of the HydroFlex runner predicted the maximum efficiency of 95.3%. It is important to note—all simulations are unsteady, and the data selected for the analysis were averaged over one complete rotation ( φ = 360 or s = 1 ) of the runner.
Since the blade profile of the HydroFlex runner is somewhat different from that of the Francis-99 runner, a direct validation can be misleading. The main purpose of using the Francis-99 runner is for relative comparison, as no experimental data of the HydroFlex runner are available. A relative (normalized) comparison with the maximum efficiency could be a good alternative. The HydroFlex runner was designed by keeping the best efficiency point at the same location (head, discharge and speed) as the Francis-99 runner. Additional peaks, apart from the best efficiency point, are believed to be numerical uncertainty. Figure 4 shows the relative deviation in hydraulic efficiency for both runners. The efficiency at all operating points of the Francis-99 runner was normalized by the maximum efficiency of 93.5%. The efficiency of the HydroFlex runner was normalized by the maximum efficiency of 95.3%. The speed factor, n E D , of 0.18 represents the synchronous speed of the turbine. Model test data of the Francis-99 (Figure 4a) shows a deviation of 25%, which corresponds to deep part load conditions, i.e., a 40 % guide vane opening at the higher n E D side. In the case of the HydroFlex runner (Figure 4b), the trend of efficiency across the different n E D values is similar to the Francis-99 runner. Figure 4c shows the exact difference in efficiency values for the corresponding n E D and Q E D . On the left side ( n E D = 0.15 0.17 ) , the numerical model shows positive errors in the efficiency, whereas on the right side ( n E D = 0.18 0.22 ), it shows negative errors in the efficiency. The numerical model seems to struggle with predicting the realistic flow phenomena at higher n E D values, although the simulations were unsteady in nature.
In this work, pressure-based boundary conditions were prescribed for all simulations; therefore, head values were predicted accurately, at 30±0.1 m. However, the reduced numerical error and the minimized fluctuations in runner torque are very important as they reflect the global performance of the runner. Figure 5 shows the standard deviation of torque for the HydroFlex runner. The torque is normalized by the corresponding value at the maximum efficiency point. The standard deviation σ corresponds to samples acquired during one complete rotation ( φ = 360 ) of the runner. The deviation is high for some operating conditions, especially the 40%, 50%, 130% and 140% guide vane positions and n E D = 0.19 0.22 . Unsteady torque fluctuations at some of these operating points are presented in Figure 6. The scale on the y-axis is different so as to preserve clarity in the high frequency fluctuations. The pattern of torque fluctuations across different operating conditions is quite interesting. Two operating points at 40% guide vane position, two operating points at 70% guide vane position and one operating point ( n E D = 0.18 ) at each 90%, 100%, 120% and 140% guide vane positions are shown in the figure. Two operating points at 40%, which guide vane opening, show quite different behaviors. At n E D = 0.18 , the amplitudes ( f 3 f n ) of torque fluctuations are small; however, at n E D = 0.20 , the amplitudes ( f = 0.5 f n ) are predominantly high. While comparing the results at the same speed factor ( n E D = 0.18 ), but different load/guide vane positions, the results are even more surprising, especially the frequency variation. Stochastic fluctuations are predominant at high load conditions and the 120% and 140% guide vane positions. At the location of maximum efficiency (Figure 6f), the amplitudes of torque fluctuations are moderate, and the frequency corresponds to the runner blades and the rotor–stator interactions. While comparing the signature of torque fluctuations, low frequency oscillations ( f = 0.5 0.7 Hz) were recorded at almost all operating points. Presently, it is unclear what causes such low frequency oscillations. This could be valuable to investigate in the future with longer simulation time, 20 – 25 revolutions of the runner, which would allow a large enough number of samples for spectral analysis and for examining the signature of low frequency fluctuations.

4. Conclusions

The present work focused on the numerical study of a newly designed Francis runner with a specific emphasis on verification. For many situations, especially automatic optimization, experimental data are not available to validate the numerical results. In such cases, relative verification is important with the reference case. For this work, we considered the Francis-99 runner as a reference case and optimized the new runner, HydroFlex, for flexible operations. The best efficiency point of the HydroFlex runner was maintained at the same location as the Francis-99 runner. The final optimized geometry of the HydroFlex runner was considered for the hill diagram simulations. The overall trend of efficiency variation across the n E D and Q E D values was found to be identical. This indicated that the numerical model performs consistently over the turbine operating range. However, the efficiency is a combination of several quantities, such as head, flow rate and torque, which can be misleading due to arithmetic cancellation of errors in these quantities. It is also important to examine the individual quantity, mainly torque, and its unsteady behavior over the hill diagram points. This will allow us to identify unstable regions. Only low frequency torque fluctuations exhibited at the 40% guide vane opening and n E D = 0.20 (see Figure 6b) were almost 100 times that of the best efficiency point. On the other hand, torque fluctuations at all other points were of high frequency. When examining the characteristics of torque closely, a quite interesting signature was observed along the synchronous speed line of the turbine ( n E D = 0.18 ). It would be interesting to see a similar signature from the experiments in the future and validate the numerical results. Overall, investigations of unsteady torque fluctuations, in addition to efficiency, were found to be useful during the runner optimization process and the identification of the stable regions of turbine operation or the need for further optimization.
It is important to note that the simulations were unsteady, and that the computational domains carry 16 million nodes. For industry-scale simulations, such simulations may be expensive and time consuming. However, it will be worth to check a few operating conditions, especially along the synchronous speed, with unsteady simulations and scale-resolving models. This will help to identify the torque characteristics. The authors are also of the opinion that the appropriate method for one particular application depends on the expectations of the designer and available computational resources. The authors tend to favor hybrid models evolving continuously from the RANS (Reynolds-averaged Navier–Stokes) to the LES (large eddy simulation) mode. These models allow us to characterize the inter-blade vortices in the blade channels and vortex breakdown in the draft tube. Perhaps, the time has come to use a fully customized approach—a combination of steady RANS, unsteady RANS, SAS-SST, SBES (stress blended eddy simulation) and RANS-LES for selected points of the hill diagram.
Supplementary data of flow parameters for both runners are provided in the Appendix A.

Author Contributions

Optimization, I.I.; conceptualization, C.T.; methodology, C.T.; software, C.T.; validation, C.T.; formal analysis, C.T.; investigation, C.T.; resources, C.T.; data curation, C.T.; writing—original draft preparation, C.T.; writing—review, I.I., O.G.D. and C.T.; visualization, C.T.; project administration, C.T.; funding acquisition, O.G.D. All authors have read and agreed to the published version of the manuscript.

Funding

This project has received funding from the European Union’s Horizon 2020 “Secure, Clean and Efficient Energy” program, H2020-LCE-07-2016-2017, under grant agreement no 764011. Project: Increasing the value of hydropower through increased flexibility — HydroFlex (www.h2020hydroflex.eu). The computational resources are used under the Notur/Norstore project (number nn9504k)—Numerical investigations of a Francis turbine.

Conflicts of Interest

The authors declare no conflict of interest and the funders had no role in the design of the study in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
LESLarge eddy simulation
RANSReynolds-averaged Navier–Stokes equation
SASScale adaptive simulation
SBESStress blended eddy simulation
SSTShear stress transport
DRunner reference diameter (m), D = 0.349 m
e ^ r η Relative error in efficiency with respect to maximum efficiency
fFrequency (Hz)
gGravity (m s−2), g = 9.821465 m s−2
HHead (m)
NNumber of data points
N Q E Specific speed, N Q E = n E D Q E D
N s Specific speed, N s = N P / H 5 / 4 (rpm, kW, m)
n E D Speed factor
PPower (W)
pPressure (Pa)
QFlow rate (m3 s−1)
Q E D Discharge factor
sRunner pitch
TTorque (Nm)
tTime (s)
vVelocity (m s−1)
α Guide vane opening position (%)
η Efficiency
σ Standard deviation

Appendix A

Equations used to compute numerical errors and hill diagram and the selected data.
H = p a b s 1 p a b s 2 ρ g + v 1 2 v 2 2 2 g
P 1 = ρ g H Q
P 2 = 2 π n T 60
η = P 2 P 1
n E D = n D g H
Q E D = Q D 2 g H
e ^ r η = η max η η max
e ^ r η η = e ^ r η HydroFlex e ^ r η Francis 99
T + = T ˜ ( t ) T ¯
σ T * = T i T ¯ N
Table A1. Flow parameters correspond to the Francis-99 runner based on model acceptance tests.
Table A1. Flow parameters correspond to the Francis-99 runner based on model acceptance tests.
n E D Q E D n (rpm)Q (m3 s−1) p 1 (Pa) p 2 (Pa)H (m)T (Nm) η (%) α (°)
0.1800.063531.500.133355,02560,63329.99603.4086.1203.95
0.1800.080531.400.168355,82561,52630.04792.589.2805.01
0.1800.095531.400.199356,32562,52330.03955.590.9706.02
0.1800.110531.400.230353,72560,32230.061126.1092.2906.99
0.1800.125531.500.261353,62561,31030.021288.9093.2308.00
0.1800.139531.400.291353,72562,21730.011436.9093.4709.01
0.1800.153531.500.319350,72560,14330.001578.4093.52010.02
0.1800.166531.600.347352,02562,75630.021710.6093.16010.99
0.1800.179531.600.374349,82560,74030.091843.8092.97012.00
0.1800.191531.700.399348,92561,14530.041952.3092.40013.05
0.1800.202531.700.423348,32561,54930.032054.2091.76014.02
Table A2. Flow parameters correspond to the HydroFlex runner based on numerical simulations.
Table A2. Flow parameters correspond to the HydroFlex runner based on numerical simulations.
n E D Q E D n (rpm)Q (m3 s−1) p 1 (Pa) p 2 (Pa)H (m)T (Nm) η (%) α (°)
0.1800.057531.500.118353,84060,57929.98549.8788.1203.95
0.1800.073531.400.153354,39061,47229.99722.0490.6385.01
0.1800.087531.400.181354,40062,50829.93879.2792.1376.02
0.1800.100531.400.209351,12060,39029.861020.4092.6626.99
0.1800.113531.500.237349,95061,53629.981247.7092.6258.00
0.1800.130531.400.272348,80062,65829.911353.4092.0959.01
0.1800.148531.500.310371,85059,70330.011530.0095.30910.02
0.1800.156531.600.326350,45062,28829.891641.9095.29510.99
0.1800.169531.600.352348,00060,14529.941771.2095.14812.00
0.1810.180531.700.374346,09060,75629.751838.1093.90613.05
0.1810.191531.700.397345,39061,15029.721930.7093.00814.02

References

  1. Iliev, I.; Trivedi, C.; Dahlhaug, O.G. Variable-speed operation of Francis turbines: A review of the perspectives and challenges. Renew. Sustain. Energy Rev. 2019, 103, 109–121. [Google Scholar] [CrossRef]
  2. Quaranta, E. Optimal rotational speed of Kaplan and Francis turbines with focus on low-head hydropower applications and dataset collection. J. Hydraul. Eng. 2019, 145, 1–15. [Google Scholar] [CrossRef]
  3. Quaranta, E.; Muller, G. Optimization of undershot water wheels in very low and variable flow rate applications. J. Hydraul. Res. 2020, 58, 1–5. [Google Scholar] [CrossRef]
  4. Keck, H.; Sick, M. Thirty years of numerical flow simulation in hydraulic turbomachines. Acta Mech. 2008, 201, 211–229. [Google Scholar] [CrossRef]
  5. Celik, I.B.; Ghia, U.; Roache, P.J.; Freitas, C.J. Procedure for estimation and reporting of uncertainty due to discretization in CFD applications. J. Fluids Eng. 2008, 130, 078001. [Google Scholar] [CrossRef] [Green Version]
  6. Roache, P. Perspective: Validation—What does it mean? J. Fluids Eng. 2009, 131, 034503. [Google Scholar] [CrossRef]
  7. Roy, C.; Oberkampf, W. A comprehensive framework for verification, validation, and uncertainty quantification in scientific computing. Comput. Method Appl. Mech. Eng. 2011, 200, 2131–2144. [Google Scholar] [CrossRef]
  8. Cervantes, M.J.; Trivedi, C.; Dahlhaug, O.G.; Nielsen, T.K. Francis-99 Workshop 1: Steady operation of Francis turbines. J. Phys. Conf. Ser. 2015, 579, 011001. [Google Scholar] [CrossRef] [Green Version]
  9. Trivedi, C.; Cervantes, M.J.; Dahlhaug, O.G. Experimental and numerical studies of a high-head Francis turbine: A review of the Francis-99 test case. Energies 2016, 9, 74. [Google Scholar] [CrossRef] [Green Version]
  10. Trivedi, C.; Cervantes, M.J.; Dahlhaug, O.G. Numerical techniques applied to hydraulic turbines: A perspective review. Appl. Mech. Rev. 2016, 68, 010802. [Google Scholar] [CrossRef]
  11. Trivedi, C.; Cervantes, M.J. Fluid-structure interactions in Francis turbines: A perspective review. Renew. Sustain. Energy Rev. 2017, 68, 87–101. [Google Scholar] [CrossRef]
  12. Trivedi, C. A review on fluid structure interaction in hydraulic turbines: A focus on hydrodynamic damping. Eng. Fail. Anal. 2017, 77, 1–22. [Google Scholar] [CrossRef]
  13. Trivedi, C. Investigations of compressible turbulent flow in a high-head Francis turbine. J. Fluids Eng. 2018, 140, 011101. [Google Scholar] [CrossRef]
  14. Trivedi, C. A systematic validation of a Francis turbine under design and off-design loads. J. Verif. Valid. Uncertain. Quantif. 2019, 4, 011003. [Google Scholar] [CrossRef]
  15. Trivedi, C.; Dahlhaug, O.G. A comprehensive review of verification and validation techniques applied to hydraulic turbines. Int. J. Fluid Mach. Syst. 2019, 12, 345–367. [Google Scholar] [CrossRef]
  16. Trivedi, C.; Agnalt, E.; Dahlhaug, O.G. Experimental study of a Francis turbine under variable-speed and discharge conditions. Renew. Energy 2018, 119, 447–458. [Google Scholar] [CrossRef]
Figure 1. Francis runner and hexahedral mesh used for numerical study.
Figure 1. Francis runner and hexahedral mesh used for numerical study.
Sustainability 12 04301 g001
Figure 2. Iso-efficiency hill diagrams of the Francis-99 runner (model acceptance test).
Figure 2. Iso-efficiency hill diagrams of the Francis-99 runner (model acceptance test).
Sustainability 12 04301 g002
Figure 3. Iso-efficiency hill diagrams of newly designed HydroFlex runner (the present numerical simulations).
Figure 3. Iso-efficiency hill diagrams of newly designed HydroFlex runner (the present numerical simulations).
Sustainability 12 04301 g003
Figure 4. Relative deviation in hydraulic efficiency for all operating points of a performance hill diagram: (a) Francis-99 runner—model acceptance tests; (b) HydroFlex runner—numerical simulations; and (c) difference between the Francis-99 and HydroFlex runners.
Figure 4. Relative deviation in hydraulic efficiency for all operating points of a performance hill diagram: (a) Francis-99 runner—model acceptance tests; (b) HydroFlex runner—numerical simulations; and (c) difference between the Francis-99 and HydroFlex runners.
Sustainability 12 04301 g004
Figure 5. Normalized standard deviation in torque, HydroFlex runner.
Figure 5. Normalized standard deviation in torque, HydroFlex runner.
Sustainability 12 04301 g005
Figure 6. Signature of torque fluctuations at selected operating points across the hill diagram of the HydroFlex runner. The scale of the y-axis is kept different to maintain clarity in the high frequency fluctuations and the amplitude level.
Figure 6. Signature of torque fluctuations at selected operating points across the hill diagram of the HydroFlex runner. The scale of the y-axis is kept different to maintain clarity in the high frequency fluctuations and the amplitude level.
Sustainability 12 04301 g006
Table 1. Boundary conditions and the numerical parameters enabled for the simulations. SAS-SST: scale adaptive simulation-shear stress transport.
Table 1. Boundary conditions and the numerical parameters enabled for the simulations. SAS-SST: scale adaptive simulation-shear stress transport.
ParametersDescription
Meshspiral casing—3.56 million, guide vanes—4.75 million,
runner—5.63 million, draft tube—2.9 million. 0.1 < y + < 30
Boundary typesTotal pressure inlet and static pressure outlet
Turbulence modelSAS-SST
Advection schemeHigh Resolution
Time marching schemeSecond order backward Euler
TimeTime step: 1° of runner rotation. Total time: three rotations

Share and Cite

MDPI and ACS Style

Trivedi, C.; Iliev, I.; Dahlhaug, O.G. Numerical Study of a Francis Turbine over Wide Operating Range: Some Practical Aspects of Verification. Sustainability 2020, 12, 4301. https://doi.org/10.3390/su12104301

AMA Style

Trivedi C, Iliev I, Dahlhaug OG. Numerical Study of a Francis Turbine over Wide Operating Range: Some Practical Aspects of Verification. Sustainability. 2020; 12(10):4301. https://doi.org/10.3390/su12104301

Chicago/Turabian Style

Trivedi, Chirag, Igor Iliev, and Ole Gunnar Dahlhaug. 2020. "Numerical Study of a Francis Turbine over Wide Operating Range: Some Practical Aspects of Verification" Sustainability 12, no. 10: 4301. https://doi.org/10.3390/su12104301

APA Style

Trivedi, C., Iliev, I., & Dahlhaug, O. G. (2020). Numerical Study of a Francis Turbine over Wide Operating Range: Some Practical Aspects of Verification. Sustainability, 12(10), 4301. https://doi.org/10.3390/su12104301

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop