1. Introduction
The depletion of fossil fuels on a global scale and the processing of fosil fuels (heavy oils, coal, petroleum, etc.) releases environmentally harmful gas emissions (so called Greenhouse Gases which include but are not limite to Carbon dioxide – CO
2, Methane – CH
4, Nitrous Oxide – N
2O. and etc.) [
1] are the main motives for the search of new energy sources. These Greenhouse gases absorb the heat from the Earth’s surface and prevent its release into outer space. Most of the Greenhouse gases occurre naturally in the atmosphere, but human activities, especially during the Industrial Revolution, have contributed significantly to their accumulation. As a result, the Greenhouse effect raises the average temperatures and the climate of our planet is changing at alarming rates. This could cause extreme weather events, increased spreading of disieases and sea levels rising due to the ice caps melting.
In recent years, renewable energy a.k.a. Green energy gains popularity. The benefits of Green energy include inexhaustible supply, power independence, carbon-free energy production and an overall cleaner and healthier environment. Additionally, it helps to diversify the energy supply and reduces the dependence on imported fuels. The five major Green energy sources are: sun, wind, flowing water, biomass from plants, and geothermal energy. Solar power is among the cleanest energy sources and it has no harmful products produced in the process of electricity generation, but it’s still expensive compared to other methods of energy production because of it’s low power density and additionally the photovoltaic effect has an overall low efficiency when transforming the solar energy into electrical energy. Futhermore, solar power systems are not a constant source of energy and has a very variable energy production even during the daylight hours of the day. Wind energy is harnessed with wind turbines. They are installed in windy regions with an open farfield such as: offshore lines, hills, plains, and so on. Wind turbines can generate more electricity than solar panels. However, wind power energy generation is even more difficult to predict than solar, since wind often changes its speed and direction and in some cases it is absent and there is no generation at all.
When it comes to water, the energy generation is more consistant and predictable. Rivers normally have a steady flow, ensuring enough and constant supply of kinetic energy. Given that the density of water is almost а thousand times that of air, the same same power can be extracted from a water fluid flow with 1m/s of average velocity, compared to an air fluid flow with a 10m/s average velocity. Such winds are rare and can occur for a short period of time in certian locations, while stream velocity in most of the rivers vary from 0.5 to 1.1 m/s [
1,
2,
3]. Water kinetic energy is harnessed with so called hydrokinetic turbines, which work on the same principles as wind turbines. Both types of turbines are calculated according to the same theory, which allows the combining of the developed methods that increase their efficiency. The content of the article is focused on axial kinetic turbines, as the reports show that they have higher efficiency, better starting capabilities, and lower cut-in speed, compared to the vertical hydrokinetic turbines [
4,
5,
6,
7,
8,
9,
10,
11,
12].
Axial kinetic turbines are widely used in hydropower due to their simple, easy to maintain design and high efficiency. The design of the runner resembles that of a propeller turbine. It consists most often of two or three spanwise twisted blades with a variable cross section. The choice of generator is easy, since axial kinetic turbines works with high rotational speeds and low torque. Unlike reaction turbines, they are intended to generate electricity mainly from the kinetic energy of the stream. This makes them suitable for flows with no potential energy, such as ocean currents, rivers, tidal waters and man-made channels. The main advantage of the kinetic systems is that they don’t require installation of penstocks or powerhouses. In addition, they have almost negligible environmental impact.
According to the data from the reviewed literature, the global potential of the ocean currents ranges from 5000 GW (43.800 TWh/yr) to 450 GW (3932 TWh/yr [
13]). However, harvesting of the kinetic stream energy is beneficial only at sites with flow velocities in the range of 2-2.5 m/s and water depths between 5 m and 25 m [
13], as most of the hydrokinetic turbines generates maximum power at flow velocities at around 2.5 - 3m/s [
7,
8,
9,
10]. This requirments are met in limited number of sites in Russia, Australia, France, India, South Korea, New Zeland, Argentina, USA, China, Canada, and UK [
13]. On the other hand, most of the rivers have much less potential, as the monthly averaged speeds are in range 0.46 – 1.1 m/s. [
2,
14]. The influacing factors are the slope angle gradient, discharge, and the roughness of the channel. In these conditions, axial hydrokinetic turbines can operate with minimum power, as their cut-in speed is 0.6-0.7m/s [
7,
8,
9]. Some of the rivers, suitable for hydropower
Figure 1 through
Figure 3 present the main parameters of some of the commercial hydrokinetic turbines. The turbine runners of Smart Duofloat/Monofloat (
Figure 1) are attached inside a cylindrical casing that resembles a diffuser. The floating construction keeps the runner submerged at a certain depth. The generator is placed inside the rear fairing of the turbine. A debris protection is placed in front of the runner. Its purpose is to prevent fish from entering between the turbine blades. According to the power curve presented in
Figure 1 (c), the turbine can generate electrical power from 0.220 to 5 kW. The cut-in speed is 0.7m/s [
8] and the rated speed is 2.8m/s [
8]. Guinard Energys Axial kinetic turbine (
Figure 2) has a similar construction. The turbine is attached to a supporting frame, which keeps the runner submerged. The maximum output power is 3.5kW [
9] at a rated speed 3m/s [
9] (
Figure 2 (c)). The runner of the Smart Freestream turbine (
Figure 3) is attached inside a reactangular supporting frame. Its parameters are similar to the Duofloat/Monofloat turbine but in this case the maximum obtained power is at 3.1 m/s [
10]. According to the velocity data presented in
Table 1, the accumulated annual energy from one unit will be 1.928-43.829 kWh. This amount of electricity is enough to power up small residential area and support the global electricity grid.
2. Parameters Affecting Hydrokinetic Turbine Efficiency and Power Generation
The efficiency of kinetic turbines is determined by the power coefficient Cp, which represents the ratio between the power of the turbine Pef and the maximum theoretical power Pw that can be absorbed from the stream:
where P
w is a product of the power density P
s of the flow and the cross-section of the turbine runner S
r:
Power density of the air Pw is measured in W/m2 and depends on the water density ρwater and speed of the stream Cw:
The power of the turbine
Pef is product of the angular speed
ω and the torque
Mb:
Usally, the experiments are conducted on experimental benches. To reproduce the predicted power and efficiency one must use the laws of similarity of the flow thtough the model and the real kinetic turbine. This will limitate the influence of the scale effect on the performance curves of the real turbine. The main parameters that ensure the dynamic similarity are Reynolds number Re, Мах number Ma, and Strouhal number Sh, defined as [
31,
32,
33]:
In the above equations, Cw is the farfield stream velocity, L is a specific size referring to the blade chord length or the outer diameter of the runner, ν is the kinematic viscosity, ω is the angular speed of the runner, R is the outer radius of the turbine’s rotor, and as is the speed of sound. The obtained performance and power curves are presented in dimensionless form as , where λ refers to the so called tip speed reatio (TSR) of the turbine (Strouhal number from equation (7).
The average annual energy production of a particular kinetic turbine can be calculated for a specific location, following the International Electricity Commission (IEC) standards [
35]:
where
is the average stream velocity in the control section recorded every hour, 8766 is the number of hours in a Julian year,
is the efficiency of the mechanical transmission (0.92-0.96),
is the efficiency of the generator (0.87-0.9),
takes into account the heat losses in the wires (0.98-0.99),
is a coefficient that takes into account the probability that a given stream velocity will occur during the year. Before placing the hydrokinetic turbine in the river or open channel, it is necessary to evaluate the flow rate, solid content and their chemical composition. According to the studies, the most suitable channels are the ones in which strong turbulent currents develop and where the Reynolds number
reaches a value up to 10
5. The Froude number
determines the type of stream. At
(below critical flow), velocity of the waves is greater than the stream velocity (gravitational forces are greater than inertial forces). Hydrokinetic turbines are most often placed in sub critical streams to avoid the formation of large waves and eddies behind the runner. Absorbing part of the kinetic energy will reduce the stream velocity and will change the upper water level.
4. Calculation of Axial Kinetic Turbine
The most used theoretical model for calculating the power curves of hydrokinetic turbines are based on the actuator disk theory [
39] and the blade element theory (BEM) [
40,
41,
42]. The oscilating disk theory presents the turbine’s runner as rotating penetrant disk that absorbs part of the kinetic energy of the stream. It is assumed that:
– the outer diameters of the oscliating disk and the runner are identical;
– the flow is stationary, incompressible and continuous;
– viscous friction forces are neglected;
– the vorticity of the flow around the disk is neglected;
– the stream velocity near the disk is the same;
– the blade cascade is composed of an infinity number of infinitely thin blades.
Figure 7 shows the principle scheme of the flow (represented by stream tube) through the actuator disk. The control area is limited by the contour of the tube.
Assuming that the flow velocity across the disk remains constant, the velocities at sections 2 and 3 are equal:
wher ρ is the density of the flow, and Sdisk is the frontal area of the actuator disk.
According to the law of conservation of momentium in a closed system, the thrust force that acts on the disk is:
Since the flow is continuous along the entire length of the stream tube, Bernoulli’s equation can be used for each sections:
where p
0 is the pressure in section 1, p
d is the pressure in section 3, and p
u the pressure in section 2. The thrust force is expressed by the pressure drop:
Equating equations (11) and (12), give the pressure drop through the disk:
Substituting equation (15) in equation (16) expressing the thrust force as function of the disk area and the stream velocity:
By equating equations (10), (11), and (17) the stream velocity in the disk plane can be written as:
The axial induction coefficient
а expresses the drop in the stream velocity from section 1 to a section in the rotor plane. Its calculated by the formula:
The dependence of C
disk on C
w1 is derived from equations (18) and (19):
The relationship between velocities C
w1 and C
w2 is established by equating equations (18) and (20):
The velocity and pressure distribution along the length of the stream tube is shown in
Figure 8. A sudden drop in the dynamic pressure Δ
p is observed after the actuator disk, due to the absorbed kinetic energy. The pressure at the end sections can be assumed to be the same.
Theoretical power produced by the disk is calculated as equation (17) is multiplied by the in-plane velocity of the disk
Cdisk:
The theoretical power coefficient of the kinetic turbine is the ratio of the theoretical power of the disk to the maximum theoretical power it can absorb from the stream:
Equating equations (20), (21), (22), and (24), the theoretical power coefficient is represented as a function of the axial induction coefficient
a:
The maximum theoretical power coefficient is calculated by equaling to zero the first derivative of Eq (25):
and it follows that Cp th max = 16/27 = 0.5926 at a = 1/3.
This value shows that an ideal kinetic turbine can extract maximum around 2/3 of the kinetic energy of the stream, and is better known as the Betz’s law [
44]. Similarly, to equation (23), the thrust force coefficient is derived as:
Using equations (17), (20), and (26) the thrust force coefficient can be expressed as function of the axial induction coefficient:
Blade cascade of the kinetic turbine can be calculated with BEM theory. The blades are divided into separate segments of equal height. For each segment the critical parameters are calculated. These include the axial induction coefficient
a, tangential induction coefficient
b, the acting velocities and forces. The velocity triangles for one blade segment are shown in
Figure 9.
Cw is the farfield velocity in front of the runner,
a and
b are the induction coefficients, and
u is the blade linear velocity for this cylindical section. The angle
β∞ between the mean relative velocity
w∞ and the blade linear velocity
u is sum of the angle of attack
i and the pitch angle of the blade
β:
From the velocity triangle in figure 9, the dependence of
β∞ can be deduced to:
where
Va and
Vt are the induced axial and tangential velocities on the blade.
Figure 10 shows the forces that act in a given section of the working blade. The drag force vector
Fx is collinear with the mean relative velocity vector
w∞, which is rotated at angle
β∞ relative to the rotor plane. The lift force vector
Fy is perpendicular to the drag force vector
Fx. The aerodynamic coefficients of an isolated air foil,
cx and
cy, are used to calculate these forces:
where
is the chord length of the blade, and
ρ is the water density.
Mean relative velocity is calculated by the following equation:
The vector of the normal force
FN is projection of the resultant force
FR =
along the axis of the runner. The normal force is used to determine the axial load on the blade cascade at given cylindrical section. The tangential vector
FT that creates the torque
Mb is a projection of the resultant force
FR on the tangent of the circle with radius
Ri, which cuts the blade at the corresponding cylindrical section.
FN and
FT are calculated by the equations below:
The runner torque
Mb is obtained with integration of the elementary tangential forfes applied to the individual blade segments (
Figure 11):
where
is the coefficient of the tangential force, and z is the number of the blades in the blade cascade.
The mechanical power of the runner for a given operating regime is obtained by multiplying the angular velocity ω of the runner by its torque:
5. An Overview of Published Articles
This section summarizes some of the articles published within the field. The first article is presented by Marianela et. al [
45]. The experiments are conducted at the Energy and Environment Laboratory of the University of Brasilia [
46]. The runner, composed of 4 blades, has an outer diameter
D1 = 220 mm. The axis of the hub is placed 234 mm above the floor of the channel. The control chamber has dimensions of 1350x1650x2000 mm. The operating regime is regulated by an electric generator. The torque values and rotational speed are monitored by integrated torque sensor. A comparison is made with a numerical experiment using finite element methods. The mathematical model was created with Ansys CFX. The computational mesh is composed of 557 216 cells. The surrounding space is rectangular with dimensions 15
D1x5
D1x
D1. A k-omega turbulent model is used. The boundary conditions are set as a velocity inlet, a pressure outlet and free slip walls. A structured mesh is used along the contour of the turbine’s runner with minimum dimensionless distance between the contour and the cell nodes
y+ = 1.6, which is within the permissible limits of the used turbulence model. Outside the runner region is used a tetrahedral mesh due to its better adaptability to the geometry of the surrounding space. Numerical studies were performed with the full runner geometry and with the actuator line model (ALM). The results are presented in
Figure 12. The physical experiment shows a maximum value of the power coefficient
Cp = 0.37 with tip speed ratio TSR
λ = 1.5, while the numerical experiments shows
Cp = 0.39 at TSR
λ = 1.7 (Ansys CFX) and
Cp = 0.43 at TSR λ = 2.15 (ALM method).
Next article belongs to Vimal Patel et al. [
47]. They investigated an axial hydrokinetic turbine placed inside a hydrodynamic channel with dimensions 1000x900 mm. The runner has an outer diameter
D1 = 360 mm and its composed of 5 blades. The torque is measured with a dynamometer and a rope. The rotational speed is measered with a chronometer. Experiments are conducted by submerging the runner at six different depth
Hs = 160, 175, 250, 286, 320 and 360 mm. The stream velocity is around 0.7m/s. According to the results (
Figure 13), the highest value of
Cp = 0.16 is achieved at
λ = 0.75 and
H = 320 mm.
Pournazeri et al. [
48] describes in their work an optimization process for the calculating of an axial kinetic turbine with an outer diameter
D1 = 110 mm [
3]. The blades are using Göttingen 804 profile. They are calculated with two main theories - BEM [
40,
41] and Momentum theory[
42] with 949 variants of the chord length along the blade span and 572 variants of the twist angle. The objective functions are the maximum values of the power coefficient Cp and the thrust force coefficient C
T. An optimal version of the working blade (
Figure 14 (b)) was obtained, with maximum
Cp = 0.3. The studies are carried out in a hydrodynamic channel with stream velocity of 0.45 m/s. In order to validate the calculated geometric parameters, additional studies were made with an offset of the pitch angle of the blades within ± 2 deg. Maximum values of
Cp = 0.28 and
CT = 0.58 were achieved (
Figure 14 (b)).
Paul Danca et al. [
49] conduct an experimental study of an axial hydrokinetic turbine with outer diameter
D1 = 200 mm, positioned in a hydrodynamic channel with dimensions 1.55D1x1.5D1. The runner has 3 blades with profile NACA4400 (
Figure 15 (a)) [
49]. It is offset by 0.8x
D1 form the floor of the channel. The experiments are carried out at stream velocity of 1.1 m/s at 5 different operating rotational velocities (
n = 70, 100, 120, 160 and 210 min
-1). The stream velocity is measured with a Pitot Prandtl tube and a differential manometer. Maximum power of
6 W was achieved at a rotational speed of 120 min
-1 (
Figure 15 (b)).
D. M. Bucur et al. [
50] studied two runners of an axial kinetic turbine H
T1 and H
T2 with outer diameter
D1 = 200 mm with 3 profiled blades which differ in the pitch angle along the span. At H
T1, the pitch angle changes uniformly from 62.63 deg to 29 deg, while at H
T2 it changes from 57.75 deg to 17.07 deg. The experiments are carried out in a hydrodynamic channel at stream velocity
Cw = 0.6, 0.7, 0.8, 0.9 and 1 m/s. The water level in the channel is about 275 mm and the dimensions of the control area are 375x300x1015 mm. According to the results of the study (
Figure 16), maximum
Cp = 0.22 is achieved with runner H
T1 at stream velocity
Cw = 1.1 m/s, and with H
T2 Cp = 0.30. The efficiency of the runner increases with the smaller pitch angles and the large difference in chord lengths in the hub and the peripheral sections.
C. Patel et al. [
51] studied the influence of the radius curvature
Rtu of the blade’s tip of an axial hydrokinetic turbine with an outer diameter
D1 = 86 mm. According to the authors, edge vortices are forming around the leading edge at the tip of the blades, which create additional hydraulic losses thus reducing the efficiency of the turbine. The rounding of the inlet edge at the tip stabilize the flow between the blades and increases the efficiency of the turbine. Five different runners with values of
Rtu = 2, 4, 6, 8 and 10 mm were manufactured. The results are presented in
Figure 17. The lowest power coefficient Cp = 0.35 is obtained without rounding of the leading edge (
Rtu = 0mm), the highest
Cp = 0.82 is achived with
Rtu = 2 mm. At
Rtu = 4 mm C
p decrease with 0.01. Further increasing of R
tu to 10 mm decrease
Cp to 0.68.
Rareş-Andrei Chihaia et al. [
52] investigate an axial hydrokinetic turbine with casing, placed in a hydrodynamic channel with dimensions of the control section 375x300x1015 mm. The experiments are conducted individually with confusor (converged shroud) and diffusor (divergent shroud), placed at both ends of the cylindrical casing. The runner of the turbine has an outer diameter
D1 = 200 mm and is composed of 4 blades with variable cross section along the span. The hub diameter is 40 mm. The cylindrical part of the casing has a length
L = 70 mm and an internal diameter 206 mm. The conical part has a length
L = 47 mm with diameters of the end sections of 206 and 250 mm. The casing can be translated along the axis of the runner with special mechanism. The different operating modes are controlled with the help of an electromagnetic brake. The torque of the runner is transmitted to a vertical shaft via a clutch and bevel gears. It is measured by a torque sensor. Experiments are performed at constant stream velocity of 0.9 m/s. The power and the power coefficient are obtained for 5 different positions of the runner along the casing (
Figure 18. Location of the different sections (a) and (b) ; picture of the runner (c) [
52]). These positions are: inlet section of the cylindrical part (section 1-1), middle section of the cylindrical part (section 2-2), outlet section of the cylindrical part (section 3-3), middle section of the conical part (section 4-4) and the outlet section of the conical part (section 5-5). According to the results (
Figure 19), the highest power coefficient 0.94 is achieved with the help of a diffuser, placed after the runner with the leading edge of the blades overlapping section 3-3. In this case, the shroud is oriented
with its cylindrical part upstream. The lowest value of the power coefficient 0.68 is recorded orienting the conical part of the casing against the flow. The highest value of power coefficient 0.7 is reported in sections 4-4 and 5-5. As the leading edge of the blades approaches the conical part, the power coefficient drops to 0.34. The experimental data in
Figure 25 show that placing a diffuser after the outlet section of the hydrokinetic turbine significantly increases its efficiency.
The runner generates maximum power when the leading edge is placed in the middle of the cylindrical part of the casing. The maximum Cp drops about 0.2 when shroud is reversed (i.e. swap the inlet and the outlet). The best results are achieved in sections 4-4 and 5-5, where the leading edge of the blades falls in the middle and in the end of the cylindrical part of the casing.
Dharmesh D. et al. [
53] investigated an axial hydrokinetic turbine with an outer diameter
D1 = 360 mm and hub diameter 104 mm. The experiments are conducted in an open hydrodynamic channel with a width of 3600 mm, water height of 800 mm at stream velocity of 0.5 m/s [
7]. The turbine runner is composed of 5 blades. The loading system is composed of a rope, dynamometer and screw gear that adjust the force of the rope on the shaft. According to the results of the experimental studies (Table
2), maximum power coefficient
Cp = 0.178 is achieved at rotational speed 53.61 min
-1.
Nedelcul et al. [
54] present a study of an axial hydrokinetic turbine placed in front of a 155 mm long conical casing. The diameters of the cross section varies from 142 to 170 mm. The outer diameter of the runner is
D1 = 137 mm and of the hub - 80 mm. The experiments are conducted at stream velocities of C
w = 0.68, 0.78, 0.9 and 1 m/s in a hydrodynamic channel with dimensions 375x300x1015 mm. The flow rate is measured by PIV method with monitoring of the trajectory and velocity of the particles. The load on the shaft is created by an electronic brake, positioned above the water level in the channel. The torque of the runner is transmitted to the loading system by bevel gear. Particle velocity measurements are carried out using Dantec system composed of 200 mJ Litron laser and FlowSenseEO_4M 32 CDD camera with resolution of 2072x2072 pixels. Five hundred pictures are recorded of each operating regime. The time interval between laser beams depends on the stream velocity. In this case it is 700 2000 µs. The laser beam is parallel to the flow and the CCD camera is placed perpendicular to the measurement plane. The results from the conducted experiments are presented in
Figure 20 in tabular and graphical form.
Optimal operating modes are marked with a red rectangle. According to the data, the highest power of 1.4 W is obtained at stream velocity Cw = 1 m/s and rotational speed n = 190 min-1. It should be noted that the calculated power coefficient reaches values greater than 1. That data is not presented in the publication which casts doubt on the accuracy of the conducted physical experiment. One possible reason is the outside diameter of the runner, which is smaller than the allowable (350 mm) for model studies of turbine runners.
Antonio C et al. [
55] describe an analytical mathematical model for the calculation of the blade cascade of an axial kinetic turbine. Three different runners with 2, 3 and 4 blades were manufactured, optimized with the presented methology. The experimental studies are carried out in an open aerodynamic tunnel with dimensions of the control section 1200x1200 mm. The numerical results are presented in
Figure 21 (a), and the optimized blade geometry at fixed number of blades and tip speed ratio, in
Figure 21 (b). The dependence of
Cp on
λ resembles a logarithmic function with a maximum value of 0.55 for an infinite number of blades, neglecting flow losses. With consideration of blade cascade hydraulic losses, the maximum power coefficient drops to 0.53. The difference in
Cp.increase with decreasing of TSR - λ. According to the data, as the number of blades
Nb and
λ increases, the chord length along the span of the blades decreases (
Figure 21 (b)), and the shape sharpens towards the tip section. The results of the physical experiment are shown in
Figure 22 (a). The experimental operating curves (
Figure 22 (a)) approach the theoretical at
λ = 1-1.75 with a maximum difference in
Cp of 0.02 at the peak point. Maximum power coefficient of 0.4 is achieved with a four bladed runner, 0.35 with a three blades runner and 0.3 with a two bladed runner.
Kyozuka Y. et al.[
56] present a numerical optimization of a three bladed axial hydrokinetic turbine with diffuser casing. The runner of the model turbine has an outer diameter
D1 = 0.3 m. It is composed of 3 profiled NACA4616 blades. The shape of the diffuser casing is optimized with a numerical design of experiment (DOE) with 64 control points based on artificial neural network developed at the von Karman Institute for Fluid Dynamics[
57]. The objective function is the power coefficient
Cp. The numerical model is composed of 450000 cells (70000 tetrahedral and 380000 hexacore -
Figure 23 (a)). A k-omega turbulence model is used. The dimensionless height of the first layer cells along the inner contour of the casing is y
+ = 20. The inlet section of the surrounding space is located at a distance of 2.7x
D1 from the cylindrical part of the diffuser casing, and the outlet section at 3.3x
D1 from the conical part of the diffuser casing. The boundary conditions are: Velocity inlet with velocity set at 1m/s and Pressure outlet with a pressure set at 0 Pa. The different operating modes of the turbine are calculated by the Frozen Rotor method with the technique of sliding mesh. The velocity and pressure distribution in the casing and the torque of the turbine runner are calculated by solving the Naiver Stokes differential equations. The results of the numerical optimization are presented in figure 38 in tabular and graphical form. The shape of the runner’s fairing and diffuser casing are determined with 7 control points.
The physical experiments are conducted in two hydrodynamic channels with dimensions 6x1.5x1.3 m and 65x5x7 m in the laboratory of Applied Mechanics at RIAM Kyushu University [
58]. The diameter of the model runner is
D1 = 0.3 m, and the optimal external diameter of the conical part of the casing is 1.64x
D1 with an inclination angle of 64 deg. Adhering the optimization data, the inlet part of the casing has the contour of a diffusor with slope angle 15.9 deg. The operating curves of the model turbine with an without a diffuser casing are obtained for several operating modes at stream flow
Cw = 0.5, 0.7 and 0.9 m/s. The results are presented in
Figure 23. The maximum power coefficient is 0.2 at tip speed ratio
λ = 2, and those with a casing have power coefficient 1.2 at λ = 3.
Ohya O et al. [
59] present optimization of the casing of a horizontal axis wind turbine (HAWT) using an innovative diffuser that increases the efficiency of the turbine by 2.5 times. With slight modifications, the design can be used in the field of water kinetic turbines. The physical experiments are conducted in wind tunnel with dimensions 3.6x2x15 m at the institute of Applied Mechanics, Kyushu University in Japan [
58]. The studied shapes of the casing have a rectangular cross section. These shapes are a confusor, a straight cylindrical pipe and a diffuser. According to the data (
Figure 24), the confusor creates a sharp drop in the velocity in the inlet section. The flow gradually recovers its initial value along the confusor until it reaches the exit section. A sharp peak in the static pressure is reported near the inlet section. The reason for this is the large frontal resistance created by the narrowest cross section. It determines the amount of flow that can pass through the casing while the rest of it tends to bypass the confusor. The straight pipeline casing does not significantly affect the flow parameters. The parameters of the stream along the length of the diffusor casing change opposite to those in the confusor. An increase in the dimensionless velocity in the inlet section of 1.8 is observed, which gradually decrease to 0.6 in the outlet section. The pressure coefficient in the inlet section reaches negative value of -1.8, which is an indicator of the creation of vacuum. Approaching the exit section, the pressure reaches the environmental value. The authors claim that the reason for the local increase in velocity at the inlet of the diffuser casing is the vortex wake that growing towards the exit section.
Figure 25 (a) shows the results from the experimental study of the dependence of the inlet velocity on the dimensionless length of the casing
L/
D, with different modifications of the geometry. The velocity increases along with the L/D ratio. According to the article, placing a flange on the outer edge of the diffuser near its exit section increases the velocity even more, as a result of the larger pulling force. The researchers claim that the reason is the formation of Kármán vortices around the outlet section of the diffuser, which create additional underpressure (
Figure 25 (b)). On theother hand, the excessive length of the diffuser prevents its incorporation in hydrokinetic turbines. From a technical and economical point of view, it is more advantageous to use a short diffuser with a flange. The researchers found that there are no significant differences in the parameters of the flow through a casing with rectangular and circular cross section.
Figure 26 shows the influence of the relative height of the flange
h/
D on the flow parameters through the turbine. According to the results, the optimal flange length is
h/
D = 0.25. The minimum value of the dimensionless velocity is 1.4 and the maximum value 1.75. The peak point of the pressure coefficient C
P is around x/L = 0.25.
The performance characteristics of the kinetic turbine with and without diffuser casing are shown in dimensionless form in
Figure 27. The maximum value of the power coefficient without casing is about 0.28 at tip speed ratio
λ = 2.4. The diffuser casing increases the maximum power coefficient up to 1.4 (nearly 5 times).