Next Article in Journal
Optimization of Design Variables of a Phase Change Material Storage Tank and Comparison of a 2D Implicit vs. 2D Explicit Model
Previous Article in Journal
De-Risking Wood-Based Bioenergy Development in Remote and Indigenous Communities in Canada
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental and Numerical Studies on the Influence of Blade Number in a Small Water Turbine

by
Grzegorz Peczkis
1,*,
Piotr Wiśniewski
1 and
Andriy Zahorulko
2
1
Department of Power Engineering and Turbomachinery, The Silesian University of Technology, 44-100 Gliwice, Poland
2
Volodymyr Martsynkovskyy Computational Mechanics Department, Sumy State University, 40007 Sumy, Ukraine
*
Author to whom correspondence should be addressed.
Energies 2021, 14(9), 2604; https://doi.org/10.3390/en14092604
Submission received: 23 March 2021 / Revised: 16 April 2021 / Accepted: 26 April 2021 / Published: 2 May 2021

Abstract

:
This paper demonstrates the procedure of blade adjustment in a Kaplan-type water turbine, based on calculations of the flow system. The geometrical adjustment of a twisted blade with varying chord length is described in the study. Computational fluid dynamics (CFD) analysis was used to characterise aerofoil and turbine performance. Furthermore, two turbines, with a different number of blades, were designed, manufactured, and tested experimentally. The numerical model results were then compared with the experimental data. The studies were carried out with different rotational velocities and different stator blade incidence angles. The paper shows a comparison of the turbine efficiencies that were assessed, using numerical and experimental methods, of a flow system with four- and five-bladed rotors. The numerical model results matched up well with those of the experimental study. The efficiency of the proposed turbines reached up to 72% and 84% for four-bladed and five-bladed designs, respectively. These efficiencies, when considered with the turbine’s simplicity, low production and maintenance costs, as well as their potential for harvesting energy from low energy flows, mean that Kaplan turbines provide a promising technology for processing renewable energy.

1. Introduction

In recent times, there has been a strong interest in replacing conventional power plants, powered by fossil fuels, with renewable energy sources [1]. Individually, many European countries have pushed development in their energy sectors towards nuclear, wind, and solar power. The introduction of biomass and weather-fuelled power generation to the energy grid has increased expectations of emission-free energy systems. Hydropower technologies can help to fulfil these expectations. In well-developed countries, hydropower technologies are already commonly used to harvest power from high-energy water sources; however, lower-energy water sources are mostly not [2]. Consequently, small water turbines are experiencing a renaissance [3,4,5].
The need for harvesting energy from relatively small water flows, which has been disregarded for years, creates challenges for engineers to design systems since ready-made algorithms are not available in the literature. New water turbine systems have to be designed since the operating conditions are very different than those used for large water flows. Many innovative, small turbines are now being developed to fill the previous void; one example of these is jet turbines. New prototypes are also being developed that combine new design concepts with tried and tested classical designs. Due to the ongoing improvement of numerical tools, it has become much easier to improve the efficiency of turbine designs. However, this has resulted in the temptation to continuously make small improvements, which, for relatively simple and cheap constructions, is not economically justified [6].
Results from Zhao et al. [7] have shown that small turbines present lots of potential for generating high-efficiency hydropower, whilst also showing increased flexibility over larger turbines during times of lower power demand.
Due to legislation that affects conventional power plants, the rise of environmental awareness, and the increase of fossil fuel cost, the price of energy continues to rise steadily. The dependence of energy cost on current weather conditions has become increasingly noticeable. Due to the restrictions in production and exceptionally favourable atmospheric conditions, Germany exceeded 55% of power production from renewable energy sources in 2020 [8]. However, in addition to a positive change for the environment, this is also a large challenge for the energy grid, which has to deal with lower demand for power. Small water turbines are a promising solution to increasing the flexibility of the power system so that it can deal with periods of lower demand. This is because small turbine units work with small water flows and can therefore be turned on at any time to meet the power demand caused by the change of energy production, which may be caused by weather-condition-dependent units. As in the case of other renewable energy sources, their effectiveness is going to depend on the scale of water stream utilisation. This reasoning prompted the authors to investigate small water turbine design, specifically to present the procedure of blade formation of water turbines for small power units.
The Kaplan-type water turbine presented in this paper is shown to harvest energy successfully in a way that does not emit harmful substances into the atmosphere. The development of small water turbines is important because of the future possibility of storing energy in pumped storage power plants using water, also in closed circuits. According to the authors, this is a real, structurally and operationally proven alternative to temporary power generation from fossil fuel sources. The presented turbine can generate electricity from the water flow within five minutes of activation and can be operated virtually. The presented, structurally novel 2 kW turbine is an example of the application of the computational algorithm quoted in this paper and the authors’ observations and described methods.

2. Blade

In the blade design process of the rotor, used in axial water turbines, aerofoils are usually selected from a database [9]. The sub-profiles are preselected based on the rotor application, for example, based on the working fluid, i.e., air or water, and the working parameters, i.e., Reynolds number. The sub-profiles are customarily catalogued based on the nomenclature adopted by research centres from all over the world [10], an example of which are the aerofoils that have been accepted by the National Advisory Committee for Aeronautics (NACA), which are described with several digits that distinctly specify the sub-profile geometry. The description of the characteristic elements of the sub-profile is shown in Figure 1. The research centres, which have manufactured and studied water turbines and pumps, recommend using the four-digit NACA sub-profiles. The adjustments of sub-profiles, based on the four- and five-digit NACA sub-profiles, are shown in this work. The modified sub-profiles are the results of analytical curves.
The aerofoil performance is customarily described with lift and drag coefficient curves, as a function of the angle of attack. The lift and drag coefficients are non-dimensional values that relate the lift and drag forces to the blade size and flow conditions around it. These coefficients are described as follows:
C l = L 1 2 ρ u 2 c
C d = D 1 2 ρ u 2 c
where L is the lift force, D is the drag force, ρ is the fluid density, u is the fluid velocity, and c is the chord length, which, for 2D applications, has a dimension of m2. It must be pointed out that the available databases of aerofoil characteristics are based on airflow, and due to the large difference in fluid viscosity, they usually cannot be used for water-based applications. However, the characteristics can be recomputed while maintaining a similarity number, named the Reynolds number, that is given by the following equation:
R e = u · c ν
where u is the free stream velocity, c denotes the chord length, and ν is the kinetic viscosity (equal to 1.46 · 10 5 m 2 s for air and to 1.07 · 10 6 m 2 s for water, if the temperature is equal to 20 °C). Due to the difference in viscosity, flow character changes, which means that the Reynolds number differs by an order of magnitude depending on whether the fluid used is air or water. Therefore, it is important to acknowledge this limitation in the aerofoils database.
The procedure of geometry change is shown, in which the length and thickness proportionality is not preserved. The geometrical adjustment is required to achieve the desired values of the chord length L and the thickness s. These optimal values are obtained from the flow system calculations. As a result, a modified sub-profile for a specific machine is created.
In the case of water turbines, one blade combines different aerofoils, whose weight centre is split radially. Subsequent aerofoils have different leaning angles β. Blade length is calculated from the empirical relations available in the literature [11,12,13,14,15]. The blade thickness is limited by material strength and the blade hub connection.
Blade geometry is transformed through subsequent aerofoils (see Figure 2), and the number of sub-profiles defines its twist and deformation rate. The number of sub-profiles must be chosen carefully to avoid blade face warping due to the imperfection of manufacturing processes [16,17].
In the design of water turbine blades, the first step is to obtain an optimal solution in regard to the water flow, followed by designing the blade hub connection [18]. The connection is defined by the constructional solution and the incidence angle regulation system. The size and construction of the mounting are determined by the material strength to withstand the forces exerted by the fluid. The mechanism of blade incidence angle change must give the possibility of the clearance cassation to prevent blade vibrations.

3. Sub-Profile Geometry Adjustment

The original aerofoil proposed for the water flow is shown in Figure 3. Its maximum thickness occurs 30% of the way along the chord length and is 8.25% of the chord length value. For further studies, the closed sub-profile has been changed into two splines, which create the upper and the lower part of the aerofoil. In Figure 3, points P00 and P16, namely, the leading and trailing points, are common for both curves.
The numerical analysis of the flow defines the blade chord length in subsequent sections. Based on the assumptions for flow calculation and the material strength, we obtain the blade thickness in the following subsections. Additionally, the vibration of excessively thin blades and water pollution (sand, mud, etc.) have to be taken into consideration. For clarity, the assumed sub-profile details have been provided in Table 1 and the coordinates of points along the aerofoil curves are shown in Table 2.
These values are the results of computations based on the Varchola algorithm [12]. If the blade size restrictions are known, the midspan sub-profile can be specified. Based on the maximum blade thickness of the original sub-profile and the desired thickness of the adjusted sub-profile, the ratio of thickness change, ks, is specified for each sub-profile as follows:
k s A = s N s A = 8.25 8 = 1.031
k s B = s N s B = 8.25 6 = 1.375
k s C = s N s C = 8.25 4 = 2.0625
where s N is the initial sub-profile thickness. The upper and lower coordinates in Y directions are scaled using the set of proportionality ratios for each sub-profile, as follows:
P A Y = P k s A
P B Y = P k s B
P C Y = P k s C
The identical procedure is applied to designate the chord length ratio, kL. The ratios are designated as follows:
k L A = L A L N = 103.74 100 = 1.0374
k L B = L B L N = 145.60 100 = 1.456
k L A = L C L N = 178.44 100 = 1.784
where L N is the initial chord length. The x coordinates of characteristic points along the upper and the lower curves are recomputed as follows:
P A X = P k L A
P B X = P k L B
P C X = P k L C
The coordinates of individual points along the sub-profile curves are shown in Table 2. Figure 4 shows the new sub-profiles that were obtained with the given procedure. In common points of upper and lower sub-profile, the round-off is carried out. The radius of the leading-edge round is equal on all sub-profiles, as well as the radius of the round on the trailing edge. The rounding-off process resulted in a slight reduction of the chord length. The round-off radius must be reasonably small to obtain a consistent blade, preventing unwanted deformations or edge chipping.
Another crucial aspect which must be considered during the blade design process is the design of an optimal blade hub mounting. To obtain an undisturbed flow, the presence of hub blade slots, sprues, or a protruding mounting is highly undesirable. In practice, we are looking for a compromise between the minimisation of flow disturbances and blade mounting machining costs.
An additional difficulty arises with the design of a consistent and clearly defined blade hub connection along the shaft. In the presented turbine, it was decided to connect the blades to the rotor hub using conical pins. Introducing additional holes in the rotor hub allowed for the possibility to adjust the blade incident angle, after the disassembly of the flow system.
The experimental studies were conducted for the obtained aerofoils and turbine geometries. The results for four- and five-bladed rotors with equal internal and external diameters were very promising. The experimental studies were conducted for both geometries with similar flow conditions.

4. Experiment

4.1. Test Rig

The turbine elements were manufactured in technical plants in Poland, whose many years of practical experience and mastery of manufacturing allowed for the production of turbomachines up to 12 meters in diameter. PA6 aluminium alloy was selected for the production of a prototype turbine. The turbine blades were made using multi-axis milling but were not subjected to surface hardening or forging since they did not need to operate in contaminated water. The plastic elements, such as the guide tube of the axial rotor, are made of polyamide. Zinc-coated screws with variable hardness, depending on the location of the screw, were used for the assembly of individual components. O-ring seals used in the turbine were made of NBR rubber since it was assumed that they would be operating at temperatures of up to 35 °C in non-contaminated water, containing no acids or bases. The bearing system was made of ceramic bearings with little frictional resistance. The sealing sleeves are individually made of BA type aluminium bronze.
Once manufactured, the water turbine was studied on a dedicated test rig at the Laboratory of Hydraulic Machines at the Silesian University of Technology. Figure 5 shows a schematic sketch of the test rig. The water, of temperature t, measured with a Pt100 thermocouple, is supplied by a one-stage pump (LFP 300 WM 270) with a varying rotational velocity, directly to a 1 m3 container, where its kinetic energy is reduced. The water flow rate is measured with an electromagnetic flowmeter. From the container, the water is directed to an open trough, before being directed into the rotor of the water turbine. The water that leaves the turbine is directed to a 90 m3 underground container. Table 3 shows the measurements from the flow meter. The minimum velocity of the supplied water, 0.5 ms−1, was achieved by reducing the diameter of the pipe downstream of the pump. The pipe sections, which were of reasonable length upstream and downstream of the flow meter, were kept to give reliable results. The water height at the inlet of the turbine’s stator was measured with a water gauge and pressure sensor.
The change in water height was achieved by changing the volume of the liquid in the underground container. As a result, the total water height, h, is the sum of the water height at the inlet to the turbine, hz, and the distance between the turbine and the underground container, hp.
h = h z + h p
The hydraulic power, Ph, i.e., the energy of water flowing through the turbine is defined as follows:
P h = g · ρ · Q ˙ · h
where g is the gravitational acceleration, ρ = f ( t ) is the density of water in the function of its temperature, and Q ˙ is the water volume flow rate.
The shaft power, Pw defined as follows:
P w = M · ω
where M is the torque, and ω is the rotational velocity. The torque and rotational velocity were measured with a strain gauge.
If the measuring range and class of the applied device are known, type 5 uncertainty can be applied to define the measurement accuracy. During the measurements, there were no significant changes in the measured values; therefore, it was assumed that the steady-state conditions were reached. The measurement of all quantities was simultaneous. Due to the interaction of the strain gauges with analogue-to-digital converters, the uncertainty of indications (C_2 = 0) that resulted from the class of transducers has been omitted—this is because they are measured with an accuracy that is an order of magnitude higher.

4.2. Measurement Uncertainty

Experimental research in the Laboratory of Fluid Flow Machines at the Department of Power Engineering and Turbomachinery has been ongoing for 60 years. During that time, the laboratory has always been equipped to carry out research for industry and academia. Measuring instruments that are used in the laboratory are systematically calibrated. Due to its size, the electromagnetic flowmeter used in the measurement (DM200) was calibrated outside of the laboratory. Calibration of the flow meter was performed using the volumetric method with a gravitational inflow of water. Due to the force required, the strain gauge was calibrated in the laboratory with a standard load over a length of 1 m. For the measurement of the rotational speed, the calibration was performed using a strobe lamp. The measurements that are recorded by the instruments are therefore validated by the manufacturer’s calibrations. The uncertainty of relevant quantities for turbine efficiency study is shown in Table 3.
The uncertainty determined further is related to the uncertainty of measuring instruments [19]. For the measuring scope used, the type B uncertainty is calculated for a rectangular distribution, as shown below:
u x = Δ x 3
Table 3. Relevant quantities measurement for turbine efficiency calculation.
Table 3. Relevant quantities measurement for turbine efficiency calculation.
Measured QuantityDesignationUnitMinimal Scope
Δ x
Measuring RangeStandard Uncertainty u x
Height h m 0.001 1 0.005774
Volume flow rate Q ˙ m 3 min 0.036 0–18 0.020785
Torque M Nm 0.02 20 0.011547
Rotational speed n RPM 2 6000 0.154701
Water temperature t °C0.2 50 ÷ 270 0.11547
For the demonstration case, the uncertainty of the water temperature measurement is neglected, and its constant density is assumed to equal to ρ = 998.62   k g 1 m 3 for t = 18 °C. The following equation includes the measured quantities needed to compute the combined standard uncertainty. The results of the following reasoning are shown in Table 4:
η T = M · π · n   30   · g · ρ · Q ˙ · h · 100 %
Partial derivatives of efficiency with respect to the subsequently measured quantities take the following form:
d η T d M = π · n   30   · g · ρ · Q ˙ · h · 100 %
d η T d n = π · M   30   · g · ρ · Q ˙ · h · 100 %
d η T d Q ˙ = M · π · n   30   · g · ρ · Q ˙ 2 · h · 100 %
d η T d h = M · π · n   30   · g · ρ · Q ˙ · h 2 · 100 %
Combined standard uncertainty is computed employing a general equation for x i = M , n , Q , h as follows:
u C ( η T ) = x i [ η T ( x i ) ( x i ) ] 2 · u x 2 ( x i )
By substituting the values, it is possible to obtain the efficiency of a turbine and the measurement uncertainty η T = 83.97 ( 1.57 ) % .

4.3. Experimental Study

The studies were performed for five-bladed rotors, as well as for an alternative solution with four blades (see Figure 6). The efficiencies of the turbines, based on the experimental studies, are presented in Figure 6. In both studies, the rotors worked in the same channel. The maximum efficiency and the working range within the efficiency range between η m a x and η m a x 10 % are similar. The system with four blades is more greatly affected by the stator blade incidence angle than the system with five blades. The five-blade system shows less sensitivity to the initial flow swirl. The sudden drop in the efficiency for the stator blade incidence angle, of around 60°, is induced by losses that are caused by the disturbance of the flow and the flow parameter drop at the inlet to the turbine. In the presented system, the maximum reasonable incidence angle is 55°. The high value of power that was obtained, 1.2 kW, and efficiency for a pico turbine, are the results of precise machining of the system’s elements and accurate blade machining, which resulted in smooth PA6 aluminium blades. In the proposed turbine, the water flowing through the bearing system provides cooling and lubrication. The construction of the turbine eliminates the necessity for seals. The power was measured by a tensometric manometer on the turbine’s hub, using the oil brake [20,21].

5. Numerical Studies

5.1. Blade Numerical Studies

To define the aerofoil characteristics, the numerical analysis of the sub-profiles obtained by the described method was performed. The numerical studies were carried out using commercial software, ANSYS Fluent (ver.2020R1, ANSYS Inc., Canonsburg, PA, USA). A mesh was created for each sub-profile. The domain size was chosen based on values from the literature, although it must be acknowledged that the domain size used in the literature was chosen for an airflow-based study [22,23,24]. A boundary layer was introduced to ensure that the y+ function value was approximately 1. The growth ratio of the following layers was equal to 1.05. The turbulence effect is modelled using the k–ω shear–stress transport (SST) model proposed by Mentner [25,26]. Figure 7 shows the domain size and the mesh boundary layer. The number of elements for sub-profiles A, B, and C is as follows: 460182, 466372, 472470.
Several simulations were conducted to calculate the lift and drag coefficients as a function of the incidence angle. The water was assumed to be incompressible and the energy equation was neglected. These assumptions are justified by the low flow velocities and negligible changes in the fluid pressure. The calculations were carried out for a range of incidence angles that are characteristic of a water turbine flow system—between 0° and 25°. The inflow velocity for the studies was assumed to be 3.5 ms−1, which is the maximal meridional velocity in the turbine, based on the empirical calculations. Figure 8 shows the drag and lift coefficient curves for each aerofoil. There is a sudden drop in the lift coefficient when the angle of attack equals 10°, which is induced by the boundary layer separation. This phenomenon is called a ‘stall’.

5.2. Flow System Numerical Studies

The flow system was modelled and assessed. In the presented case, the calculations were carried out for the whole system, which accounted for the inflow of the water, the stator stage with adjustable blades, and the outflow pipe. The outflow pipe was designed as a cylinder, in line with literature recommendations. In the literature, multi-stage outflow diffusers are usually recommended. In the described case, the drop mostly occurs in the outflow channel; therefore, to avoid the cavitation, the inflow velocity should not exceed 2 ms−1.
The numerical studies were performed in commercial software ANSYS CFX (ver.2020R1, ANSYS Inc., Canonsburg, PA, USA). The model consists of connected parts of the flow system to ensure reasonable mesh density and the number of elements. The computational domain and mesh are depicted in Figure 9; the mesh consisted of ~ 8 million elements. The steady-state condition was assumed; thus, no changes of flow character with time were considered. The boundary conditions were described using the mass flow inlet and pressure outlet, and the rotational velocity was given. The turbulence effect is modelled using the k–ω shear–stress-transport (SST) model proposed by Mentner. The boundary layer was introduced to ensure that the y+ function value equalled approximately 1, and the growth ratio of the following layers equalled 1.05. The flow is assumed to be isothermic, with a constant temperature of 25 °C.
Two series of computations were conducted for different inlet blade incidence angles, varying the volume flow in the range between 2.5 m3min−1 and 5.5 m3min−1 and varying rotational velocity within the range of 300 to 1000 RPM. The computed maximum efficiency of the turbine equalled 72.5% and 84% for the stator blade incidence angles of 0° and 45°, respectively. The computational results are in good agreement with experimental studies, which, for equivalent conditions, indicated efficiencies of 70% and 84%, for stator blade incidence angles of 0° and 45°, respectively. The velocity and pressure distribution in the flow system for the case study with maximum efficiency are shown in Figure 10, which shows no local pressure changes or boundary layer separation. Figure 11 presents the pressure distribution on the rotor’s blades. The pressure distribution, visible in Figure 11, confirms the assumed lack of local pressure changes on the surface of the blades. This confirms the correct design of the variability of the geometry of the flow profiles, which allows for a higher maximum value, for a given ratio of the speed of the flow system to the rotor efficiency. This was confirmed also by empirical research that was carried out in the laboratory.

6. Conclusions

This paper presents the experimental and numerical study of a small Kaplan turbine. Within the paper, an approach is used to tailor blade profiles, which are designed for airflow applications, for use with a specific, small, power water turbine. The simplified method of blade adjustment presented in this article allows for the design of turbine blade geometry, based on the airflow reference blade from available databases. The proposed approach was applied to the design of two turbines, differing in the number of blades. These were manufactured and studied both experimentally and numerically.
The results of the experimental study for the two turbines are presented. The influence of rotational velocity and stator blade incidence angle was studied. The efficiency of designed turbines reached up to 84%. The results obtained during the experimental studies confirmed the accuracy of the proposed method for turbine blades’ design.
In the numerical studies, the lift and drag coefficients were calculated as a function of the angle of attack, which was applied to the newly created sub-profiles of the turbine blade. For the A sub-profiles, stalling occurred at an angle of attack equal to ~10°. For sub-profiles B and C, stalling did not occur; however, due to boundary layer separation, the increase of lift beyond an angle of attack equal to ~10° was shown to result in significant drag.
Numerical studies were carried out for the whole flow system. The studies were performed with varying rotational velocity and volume flow rate. The numerical tool slightly underestimated the efficiency of the turbines, when compared to the experimental study. This can be caused by the use of assumptions within the numerical tool, e.g., isothermal flow, or may have been induced by the imperfect mesh, which may have resulted in gap leakage—which, for small turbines, has a significant influence on the efficiency. This issue needs to be studied further since it opens possibilities for the development of numerical models for water machinery computations.
The value of the small turbine does not reflect the time devoted to this study; however, it confirms the utility of the described method and its implementation in the design process. The presented turbines are promising, low-cost solutions for harvesting energy from relatively small water flows with low energy potential. Such a turbine has potential in numerous, so far undeveloped, small water flows in forested areas. Moreover, the presented results, especially due to the high efficiency of mechanical power generation, confirm the potential of small turbines in energy generation, especially when the added flexibility to the power grid is considered.

Author Contributions

Conceptualization, G.P. and P.W.; methodology, G.P.; software, P.W.; validation, G.P. and A.Z.; formal analysis, G.P.; investigation, P.W.; resources, G.P.; data curation, G.P.; writing—original draft preparation, P.W.; writing—review and editing, A.Z.; visualization, P.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Symbols
AoAangle of attack, °
cchord length, m
Ccoefficient, -
Ddrag, N
g gravity acceleration, m s 2
hwater height, m
kadjustment coefficient, -
Llength, m
Llift, N
Mtorque, N m
nrotational speed, RPM
Ppower, W
P, Rpoints along blade suction and pressure side, -
Q ˙ volume flow rate, m 3 s
ReReynolds number, -
sblade thickness, m
twater temperature, °C
ufree stream velocity m s
uuncertainty, -
βblade incident angle
η efficiency, -
ν kinetic viscosity, m 2 s
ρdensity,   k g m 3
Ω angular velocity, s 1
Subscripts
1inlet
2outlet
A, B, Csub-profile
bB type
ccombined
Ddrag
hhydraulic
Llift
Llength
maxmaximum
pturbine outlet
sthickness
Tturbine
wshaft
zturbine inlet

References

  1. Jedral, W. Res and Energy efficiency in context of challenges for electricity generation and utilization of Energy in Poland. Rynek Energii 2019, 1, 3–8. [Google Scholar]
  2. Sritram, P.; Suntivarakorn, R. Comperative Study of Small Hydropower Turbine Efficiency at Low Head Water. Energy Procedia 2017, 138, 646–650. [Google Scholar] [CrossRef]
  3. Wiliamson, S.J.; Stark, B.H.; Booker, J.D. Performance of a low-head pico-hydro Turgo turbine. Appl. Energy 2013, 102, 1114–1126. [Google Scholar] [CrossRef]
  4. Du, J.; Yang, H.; Shen, Z.; Chen, J. Micro hydro power generation from water supply system in high rise buildings using pump as turbines. Energy 2017, 137, 431–440. [Google Scholar] [CrossRef]
  5. Rentschler, M.; Maraongiu, J.; Neuhauser, M.; Parkinson, E. Overview of SPH-ALE applications for hydralic turbines in ANDRITZ hydro. J. Hydrodyn. 2018, 30, 114–121. [Google Scholar] [CrossRef]
  6. Lahimer, A.A.; Alghoul, M.A.; Sopian, K.; Amin, N.; Asim, N.; Fafhel, M.I. Research and development aspects of pico-hydro power. Renew. Sustain. Energy Rev. 2012, 16, 5861–5878. [Google Scholar] [CrossRef]
  7. Zhao, J.; Wang, L.; Liu, D.; Wang, J.; Zhao, Y.; Liu, T.; Wang, H. Dynamic model of Kaplan turbine regulating system suitable for power system analysis. Math. Probl. Eng. 2015, 2015, 1–12. [Google Scholar] [CrossRef] [Green Version]
  8. ISE Fraunhofer. 2020. Available online: https://www.ise.fraunhofer.de/de/presse-und-medien/presseinformationen/2020/nettostromerzeugung-im-ersten-halbjahr-2020-rekordanteil-erneuerbarer-energien.html (accessed on 1 January 2021).
  9. Ujwala, M.; Chowdary, C.K.; Srinivas Naik, L. Design and Analysis of Low Head, Light weight Kaplan Turbine Blade. Int. Refereed J. Eng. Sci. 2017, 6, 17–25. [Google Scholar]
  10. Bruschi, D.L.; Fernandes, J.C.; Falcao, A.F.; Bergmann, C.P. Analysis of the degradation in the Wells turbine blades of the Picooscillating-water-column wave energy plant. Renew. Sustain. Energy Rev. 2019, 115, 109368. [Google Scholar] [CrossRef]
  11. Soukala, J.; Pochylyb, F.; Varcholac, M.; Parygind, A.G.; Volkovd, A.V.; Khovanovd, G.P.; Naumov, A.V. Selection of Axial Hydraulic Turbines for Low Head Microhydropower Plants. Therm. Eng. 2015, 62, 862–867. [Google Scholar] [CrossRef]
  12. Sosilo, A.K.; Hadi, H.; Soehartanto, T. Design of Hydro Power by Using Turbines Kaplan on The Discharge Channel Paiton 1 and 2. In Proceedings of the E3S Web of Conferences, Yogyakarta, Indonesia, 1 November 2017. [Google Scholar]
  13. Chen, H.; Zhou, D.; Zheng, Y.; Jiang, S.; Yu, A.; Guo, Y. Load rejection Transient Process Simulation of a Kaplan Turbine Model by Co-Adjusting Guide Vanes and Runner Blades. Energies 2018, 11, 3354. [Google Scholar] [CrossRef] [Green Version]
  14. Verchola, M. Hydraulicky Navrh Odstredivych Cerpadiel; STU: Bratislava, Slovakia, 2016. [Google Scholar]
  15. Krzyżanowski, W. Turbiny Wodne Konstrukcja i Zasady Regulacji: Wydawnictwo Naukowo-Techniczne; WNT: Warszawa, Poland, 1971. [Google Scholar]
  16. Abeykoon, C.; Hantsch, T. Design and Analysis of a Kaplan Turbine Runner Wheel. In Proceedings of the 3rd World Congress on Mechanical, Chemical and Material Engineering, Roma, Italy, 8–10 June 2017. [Google Scholar]
  17. Singara, M.; Sammartano, V.; Arico, C.; Collura, A.; Tucciarelli, T. Cross-Flow turbine design for variable operating conditions. Procedia Eng. 2014, 70, 1539–1548. [Google Scholar] [CrossRef] [Green Version]
  18. Muis, A.; Sutkino, P.; Soewono, A.; Hartono, F. Design optimization of axial gydraulic turbine for very low head application. Energy Procedia 2015, 68, 263–270. [Google Scholar] [CrossRef] [Green Version]
  19. BIPM. JCGM 100:2008: Evaluation of Measurement Data—Guide to the Expression of Uncertainty in Measurement. 1995. Available online: https://ncc.nesdis.noaa.gov/documents/documentation/JCGM_100_2008_E.pdf (accessed on 30 April 2021).
  20. Dehkarqani, A.S.; Engstrom, F.; Aidanpaa, J.O.; Cervantes, M.J. Experimental Investigation of a 10 MW Prototype Kaplan Turbine during Start-Up Operation. Energies 2019, 12, 4582. [Google Scholar] [CrossRef] [Green Version]
  21. Chirag, T.; Agnalt, E.; Dahlhaug, O.G. Experimental investigation of a Francis turbine during exiqent ramping and transition into total load rejection. J. Hydraul. Eng. 2018, 144, 04018027. [Google Scholar]
  22. Piotr, W.; Francesco, B.; Zbigniew, B.; Alessandro, B. Numerical Analysis on the Effectiveness of Gurney Flaps as Power Augmentation Devices for Airfoils Subject to a Continuous Variation of the Angle of Attack by Use of Full and Surrogate Models. Energies 2020, 13, 1877. [Google Scholar]
  23. Janjua, A.B.; Khalil, M.S.; Saeed, M. Blade profile Optimization of Kaplan Turbine Using CFD Analysis. Mehran Univ. Res. J. Eng. Technol. 2013, 32, 559–574. [Google Scholar]
  24. Borowski, D.; Węgiel, M.; Ocłoń, P.; Węgiel, T. CFD model and experimental verification of water turbine integrated with electrical generatior. Energy 2019, 185, 875–883. [Google Scholar] [CrossRef]
  25. Ansys®. ANSYS CFX-Solver Theory Guide; Ansys: Canonsburg, PA, USA, 2020. [Google Scholar]
  26. Menter, F.R. Two-equation Eddy-Viscosity Turbulence Models for Engineering Applications. AIAA J. 1994, 32, 1598–1605. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Characteristic elements of the sub-profile.
Figure 1. Characteristic elements of the sub-profile.
Energies 14 02604 g001
Figure 2. Aerofoil geometry in three subsequent cross sections. Plane view.
Figure 2. Aerofoil geometry in three subsequent cross sections. Plane view.
Energies 14 02604 g002
Figure 3. Aerodynamic sub-profile divided into upper and lower splain and the coordinates of its characteristic points.
Figure 3. Aerodynamic sub-profile divided into upper and lower splain and the coordinates of its characteristic points.
Energies 14 02604 g003
Figure 4. Blade sub-profiles obtained during the adjustment procedure (a) and their transformation through the sub-profiles spaced according to the calculations (b).
Figure 4. Blade sub-profiles obtained during the adjustment procedure (a) and their transformation through the sub-profiles spaced according to the calculations (b).
Energies 14 02604 g004
Figure 5. Test rig. 1—Kaplan water turbine, 2—strain gauge, 3—auxiliary brake, 4—hydraulic brake, 5—supplying pump, 6—regulation valve, 7—electromagnetic flow-meter, 8—underground container (80 m3), 9—inflow trough (3.9 m), 10—turbine drain pump (1.4 ÷ 2.4 m), 11—expansion tank (1.3 m3).
Figure 5. Test rig. 1—Kaplan water turbine, 2—strain gauge, 3—auxiliary brake, 4—hydraulic brake, 5—supplying pump, 6—regulation valve, 7—electromagnetic flow-meter, 8—underground container (80 m3), 9—inflow trough (3.9 m), 10—turbine drain pump (1.4 ÷ 2.4 m), 11—expansion tank (1.3 m3).
Energies 14 02604 g005
Figure 6. Kaplan turbines: four- and five-bladed rotors and turbines’ efficiency for various stator blades incidence angle for four- and five-bladed rotors.
Figure 6. Kaplan turbines: four- and five-bladed rotors and turbines’ efficiency for various stator blades incidence angle for four- and five-bladed rotors.
Energies 14 02604 g006
Figure 7. Domain size (left) and mesh boundary layer (right).
Figure 7. Domain size (left) and mesh boundary layer (right).
Energies 14 02604 g007
Figure 8. Lift and drag coefficients in the function of incidence angle for blade A, B, and C cross sections (see Figure 4).
Figure 8. Lift and drag coefficients in the function of incidence angle for blade A, B, and C cross sections (see Figure 4).
Energies 14 02604 g008
Figure 9. The schematic view of the domain and the mesh.
Figure 9. The schematic view of the domain and the mesh.
Energies 14 02604 g009
Figure 10. Water pressure and velocity in the turbine based on the numerical studies.
Figure 10. Water pressure and velocity in the turbine based on the numerical studies.
Energies 14 02604 g010
Figure 11. The pressure distribution on the rotor’s blade surfaces, based on the numerical studies.
Figure 11. The pressure distribution on the rotor’s blade surfaces, based on the numerical studies.
Energies 14 02604 g011
Table 1. Details of blade sub-profiles.
Table 1. Details of blade sub-profiles.
PositionSub-ProfileLength, mmThickness, mm
Blade rootA103.748
Blade at 50% heightB145.606
Blade tipC178.444
Table 2. Coordinates of characteristic points along with the blade upper and lower curve for the sub-profile of thickness sA = 8 mm and length LA = 103.74.
Table 2. Coordinates of characteristic points along with the blade upper and lower curve for the sub-profile of thickness sA = 8 mm and length LA = 103.74.
Point NumberChord LengthUpper SurfaceLower SurfaceThickness
n L A , mm P A Y , mm P A X , mm s , mm
000.0001.2121.2120.000
011.2972.5210.4852.036
022.5943.3940.1943.200
035.1874.6550.0484.606
047.7815.6240.0005.624
0510.3746.3030.0006.303
0615.5617.3210.0487.273
0720.7497.9520.1457.806
0831.1238.2910.2918.000
0941.4978.0970.3887.709
1051.8727.5640.3887.176
1162.2466.5940.3396.255
1272.6205.3330.2425.091
1382.9943.9270.1453.782
1493.3692.1330.0482.085
1598.5561.1640.0191.144
16103.7430.0000.0000.000
Table 4. Example of values measured simultaneously for one measuring point and the calculated partial derivative.
Table 4. Example of values measured simultaneously for one measuring point and the calculated partial derivative.
Measured QuantityDesignationUnitMeasured ValuePartial Derivative
Height h m 2.133.365
Volume flow rate Q ˙ m 3 s 0.0053551338.9
Torque M Nm 0.1062267.497
Rotational speed n RPM 7210.009944
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Peczkis, G.; Wiśniewski, P.; Zahorulko, A. Experimental and Numerical Studies on the Influence of Blade Number in a Small Water Turbine. Energies 2021, 14, 2604. https://doi.org/10.3390/en14092604

AMA Style

Peczkis G, Wiśniewski P, Zahorulko A. Experimental and Numerical Studies on the Influence of Blade Number in a Small Water Turbine. Energies. 2021; 14(9):2604. https://doi.org/10.3390/en14092604

Chicago/Turabian Style

Peczkis, Grzegorz, Piotr Wiśniewski, and Andriy Zahorulko. 2021. "Experimental and Numerical Studies on the Influence of Blade Number in a Small Water Turbine" Energies 14, no. 9: 2604. https://doi.org/10.3390/en14092604

APA Style

Peczkis, G., Wiśniewski, P., & Zahorulko, A. (2021). Experimental and Numerical Studies on the Influence of Blade Number in a Small Water Turbine. Energies, 14(9), 2604. https://doi.org/10.3390/en14092604

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