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:
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 m
2. 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:
where
is the free stream velocity,
denotes the chord length, and
is the kinetic viscosity (equal to
for air and to
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:
where
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:
The identical procedure is applied to designate the chord length ratio,
kL. The ratios are designated as follows:
where
is the initial chord length. The x coordinates of characteristic points along the upper and the lower curves are recomputed as follows:
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 m
3 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 m
3 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.
The hydraulic power,
Ph, i.e., the energy of water flowing through the turbine is defined as follows:
where
g is the gravitational acceleration,
is the density of water in the function of its temperature, and
is the water volume flow rate.
The shaft power,
Pw defined as follows:
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:
Table 3.
Relevant quantities measurement for turbine efficiency calculation.
Table 3.
Relevant quantities measurement for turbine efficiency calculation.
Measured Quantity | Designation | Unit | Minimal Scope
| Measuring Range | Standard Uncertainty |
---|
Height | | | | | |
Volume flow rate | | | | 0–18
| |
Torque | | | | | 0.011547 |
Rotational speed | | | 2 | | |
Water temperature | | °C | 0.2 | | |
For the demonstration case, the uncertainty of the water temperature measurement is neglected, and its constant density is assumed to equal to
for
°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:
Partial derivatives of efficiency with respect to the subsequently measured quantities take the following form:
Combined standard uncertainty is computed employing a general equation for
as follows:
By substituting the values, it is possible to obtain the efficiency of a turbine and the measurement uncertainty .
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
and
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
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 m
3min
−1 and 5.5 m
3min
−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.