D5.4 H SHP: Andbook of Innovative Technologies To Promote
D5.4 H SHP: Andbook of Innovative Technologies To Promote
is the projection of force F on the tangent drawn to the path of motion of the blade axis.
The summary torque includes the general component of the resistance force F
h
. The torque
moment generated by the turbine consists of the torques generated by each separate blade.
Currently only one blade will not generate positive moment (it will generate a negative moment
the resistance one). Thus, the torque generated by the proposed turbine will be essentially bigger
35
than the torque produced by the existing turbines for the same geometrical (blades dimensions)
and kinematical parameters of water. The proposed micro hydro power plant allows the
transformation of the water flow kinetic energy into mechanical or electrical energy with an
increased utilization coefficient of water energy.
In the floating micro hydro power plant (figure 28) [10] an additional centrifugal pump 2 is
mounted on the resistance structure 1. It is connected cinematically to the multi-blade rotor spindle
3 by belt transmissions 4 and 5. The electric generator 6 is connected to the multi-blade rotor
spindle by belt transmissions 4 and 7. As well, the resistance structure 1 is connected to the shore
by the metal truss 8 and supporting cables with cross-ties 9.
Figure 28 Floating micro hydropower plant with electric generator and hydraulic pump.
In the floating micro hydro power plant (figure 29, a) [11] the rotor 1 contains an odd number of
blades 2 that are fitted with the possibility of rotation on vertical axes O'O' (Figure 29,b) mounted
on the extreme end of each horizontal bar 3. On frame 4, in the front part (through which water
flows pass) a rigid bar 5 is installed on which, in front of the hydro turbine relative to the water flow
direction, a sensor 6 is fixed that determines the water flow direction and connects to the rotation
gear 7. The water flow moves in the direction of vector V
0
(figure 29, b). Angle is the entering
angle of the blades formed by the hydrodynamic surface string and the working lines of the water
flow vector V
0
V
0
. The angle depends on the form of the hydrodynamic surface and on the position
in the plane surface. By changing the water flow direction due to the change of water discharge
and river bed, the water currents will divert by angle Ao modifying the entering angle . To meet the
angle of attack, optimal in terms of conversion, it is necessary to correct the position of all blades
by angle o. When changing the water flow direction, the positioning of all blades 2 is corrected
simultaneously by angle Ao using the rotation mechanism 7:
1,2
= Ao. (3)
36
Figure 29 Floating micro hydropower plant with influence compensation of water currents flow
direction change.
In the floating micro hydro power plant (figure 30) [8,13] a technical solution is proposed
ensuring the transverse stability of platform 1 of the floating micro hydro power plant that is
mounted on floating bodies 2 and 3, placed on the same side (shore side of the rotor spindle 4).
Due to the fact that the rotor 4 blades 5 are hollow, the hydrostatic Archimedes force of the blades
5 fulfills the role of the floating bodies (figure 30,a, b).
The analysis of the application points motion path of the Archimedes force F
a
(point N in fig. 30, c)
has shown that the distance from this point to the plane that crosses the rotor spindle 4 (O
1
-O
1
,
figure 38, b) will differentiate depending on the positioning angle of the rotor. Thus, these
distances, for the blades that are placed in the upper semiplane defined by axis O1O1-OO differ
from the distances of those blades placed in the lower semiplane. The migration of the points of
application of the Archimedes force causes the pitching moment:
M
r
=M
Eas
-M
Ead
, (4)
where M
as
is the summary moment developed by the Archimedes forces that react on the blades
currently located in the upper semiplane;
M
ad
is the summary moment developed by the Archimedes forces that react on the blades
currently located in the lower semiplane.
The summary moments developed by the Archimedes forces that react on the blades are
determined by the relations:
M
as
= F
ai
l
si
and M
ad
= F
ai
l
di
(5)
where F
ai
is the Archimedes forces that react on the blades 5 currently located in the upper
semiplane;
l
si
is the distance from the point of application of the Archimedes force that reacts on the
blades 5 currently located in the upper semiplane;
37
Figure 30 Micro hydropower plant with increased transverse stability.
38
l
di
is the distance from the point of application of Archimedes force that reacts on the blades
5 currently located in the lower semiplane.
Distances l
si
and l
di
are calculated from the relation:
2 2 2
2 cos( )
M M
l R c Rc o = + + + , (6)
Where R is the rotor radius 4;
C
M
is the distance from the point of application of the Archimedes force and the blade fixing
point to the turbine rotor;
is the angle formed by the blade chord and the water flow direction;
is the angle formed by the rotor lever and the perpendicular direction on the watercourse.
To compensate the pitching moment M
r,
the rotor spindle 10 is settled in plane
'
1
'
1
O O at distance
e compared to the longitudinal axial plane of the floating bodies O
1
-O
1
. Distance e is calculated
from the relation:
1
n
i
i
y
e
n
=
=
, (7)
where n is the number of rotor blades, and y
i
is the distance from the centre of application of
Archimedes force on blade i up to the longitudinal axial plane (figure 30, c). For each distance, y
i
is
calculated by the relation:
360
cos sin( ( 1) )
o
i M
y c R i
n
o = + +
, (8)
where R is the rotor radius;
C
M
is the distance from the point of application of the Archimedes force and the blade fixing
point to the turbine impeller, O
i
N
i
in figure 30,c;
n is the number of rotor blades.
So, the distance e is calculated by the relation:
cos
M
e c o = , (9)
where is the angle formed by the blade chord and the water flow direction.
Conclusion: To ensure the floating stability of the micro hydro power plants the rotor is mounted
on the main structure with displacement e against the water stream. Thus, the micro hydro power
plants designed to be anchored on the left bank cannot be anchored on the right bank.
5.1.3. Design of the hydrodynamic rotor
5.1.3.1. Theoretical justification of the hydrodynamic profile selection of
the blade in normal section
Let consider the symmetrical profile of the blade placed in a fluid stream that moves uniformly at
velocity V
(figure 31). In the fixing point O' of the symmetrical blade with lever OO let consider
two coordinate systems, that is: the system O'xy with axis O'y oriented in the direction of the
velocity vector V
, and axis O'x - normal for this direction; and the system O'xy with axis O'y
39
oriented to the lever direction O'O, and axis O'x - normal for this direction. Point A corresponds to
the rear edge, and point B corresponds to the entering edge. The entering angle o is the angle
between the chord AB of the profile and the direction of the velocity vector V
, and the positioning
angle is the angle formed by the velocity vector direction and lever O'O
Figure 31 Hydrodynamic profile blade.
The components of the hydrodynamic force F
is flow velocity, S
p
=ch (c is the length of chord AB, and h is the blade
height) represents the area of the blade lateral surface, and C
L
and C
D
are hydrodynamic
dimensionless coefficients, called the lift coefficient and drag coefficient. The hydrodynamic
coefficients C
L
and C
D
are functions of the entering angle , Reynolds number Re and the
hydrodynamic shape of the blade profile. The components of the hydrodynamic force in the
coordinate system O'xy are
sin cos ,
cos sin .
x L D
y L D
F F F
F F F
'
'
= +
= +
(12)
The torque moment of the rotor spindle OO developed by blade i is
,
,
r i x
T F OO
'
' = (13)
and the summary torque moment developed by blades is
1
,
Npal
r ri
i
T T
E
=
=
(14)
where Npal is the number of rotor blades.
40
Generally, the hydrodynamic force has no point of application in the origin of the blade axes
system O so as it produces a resulting moment. The produced moment is determined by
comparing it to a certain point of reference. The point situated at distance of the chord from the
entering edge B will be considered as point of reference. The moment, also called the pitching
moment, is calculated according to formula
2
1
,
2
M p
M C V cS
=
(15)
where C
M
is the profile number of turns.
5.1.3.2. Determination of the hydrodynamic coefficients CL and CM.
Plane potential (cyclic) motion
The profile chord is considered unitary for simplicity. Initially, the fluid is considered incompressible
and non-viscous, and its motion plane and cyclic. In the case of an incompressible fluid in plane
motion the velocity components
( ) , V u v =
u
t
u =
}
(20)
In the relations (19, 20) s represents the measured distance of profile C, and (r, ) are the polar
coordinates of point P'(x,y) reported to the point on the contour corresponding to distance s (figure
32).
41
Figure 31 Fluid cyclic motion around profile . C
Therefore the potential in point P'(x,y) is given by formula:
( ) ( )
( ') cos sin ln( ) .
2 2
C C
q s s
P V x V y r ds ds
o o u
t t
u = + +
} }
(21)
To calculate the cyclic motion potential the collocation method is used, namely: the boundary of
profile C is approximated by a closed polygon
1
,
N
j
j
C E
=
~
sides E
j
having their points (vertex) P
j
and P
j+1
placed on C. The numbering of points starts from the
rear edge on the lower side in the direction of the entering edge, passing further to the upper side
(Fig. 33). It is considered that the intensity of vortexes (s) distributed on profile C is constant at
the boundary having value , and the intensity of sources q(s) distributed on the profile is constant
at each boundary element E
j
having value q
j
, where j=1,, N. Specifying the above, equation (21)
becomes:
1
cos sin ln( ) ,
2 2
j
N
j
j
E
q
V x V y r ds
o o u
t t
=
| |
u = + +
|
\ .
}
(22)
The unknown being and , 1, , .
j
q j N =
Figure 32 Digitization of profile . C
Let consider the boundary element E
j
with points P
j
and P
j+1
(Fig. 34). The normal and tangent unit
vectors of the element E
j
are given in formulas:
( sin ,cos ),
(cos ,sin ),
j j j
j j j
n u u
t u u
=
=
(23)
where
1
sin ,
j j
j
j
y y
l
u
+
=
1
j
cos .
j j
j
x x
l
u
+
=
42
Figure 33 Boundary element .
j
E
The unknown and q
j
, where j=1,, N from the relation (22), are determined from the boundary
conditions and Kutta condition. In the case of non-viscous fluid, the boundary condition is the
sliding condition at the profile boundary that is watertight and rigid, that, in the particular case of
plane and potential motion of the incompressible fluid, is written as follows:
0, V n =
(24)
where n
is the normal of the profile. It is necessary to satisfy the condition (24) in the points of
collocation. Points
( )
, j
j j
M x y the centers of sides E
j
, are selected as points of collocation:
1 1
, , 1, , .
2 2
j j j j
j
j
x x y y
x y j N
+ +
+ +
= = =
Velocity components in the point of collocation M
j
are written by:
( , ),
v ( , ).
j
j j
j
j j
u u x y
v x y
=
=
Thus, condition (21) delivers N algebraic relations:
sin cos 0, 1, ,
i i i i
u v i N u u + = = (25)
that are used to determine those N+1 unknown and q
j
, where j=1,, N. Kutta condition will
deliver the final relation, namely:
1
,
N
E E
V V t t =
(26)
where t
is the tangent versor of the boundary element. In our notations, condition (26) takes the
form:
1 1 1 1
cos sin cos sin .
N N N N
u v u v u u u u + = + (27)
Velocity components in point
i
M are determined by the contributions of velocities induced by the
distribution of sources and vortexes on each boundary element E
j
:
43
1 1
1 1
cos ,
sin ,
N N
s v
i j ij ij
j j
N N
s v
i j ij ij
j j
u V q u u
v V q v v
o
o
= =
= =
= + +
= + +
(28)
where , , ,
s s v v
ij ij ij ij
u v u v are so-called induction (influence) coefficients. For instance,
s
ij
u represents the
component of velocity direction x in point M
i
, induced by the unitary intensity source distribution
from the element
j
E
.
The induction coefficients can be calculated in the following way:
, 1
, 1
, 1
, 1
1
ln cos sin ,
2 2
1
ln sin cos ,
2 2
1
cos ln sin ,
2 2
1
sin ln cos ,
2 2
i j ij s
ij j j
ij
i j ij s
ij j j
ij
ij i j v
ij j j
ij
ij i j s
ij j j
ij
r
u
r
r
v
r
r
u
r
r
v
r
|
u u
t t
|
u u
t t
|
u u
t t
|
u u
t t
+
+
+
+
| |
= |
|
\ .
| |
= + |
|
\ .
| |
= |
|
\ .
| |
= + |
|
\ .
(29)
where
ij
is the angle formed by sides
j i
P M and
1 i j
M P
+
, for ij, and
ij
= , i,j=1,,N, and r
ij
is the
distance between points M
i
and P
j
. Let substitute expressions (28) and (29) in the boundary
conditions (25) and in Kutta condition (27) to obtain the linear system N+1 of equations with N+1
unknowns: and q
j
, where j=1,,N :
, 1
1
1, 1, 1 1
1
, 1, ,
,
N
ij j i N i
j
N
N j j N N N
j
A q A b i N
A q A b
+
=
+ + + +
=
+ = =
+ =
(30)
where coefficients A
ij
and b
i
, i, j = 1,, N+1 are calculated by formulas:
( ) ( )
, 1
1 1
sin ln cos , , 1, , ,
2 2
i j
ij ij ij ij
ij
r
A i j N
r
|
t t
+
| |
= A + A = |
|
\ .
( ) ( )
, 1
, 1
1
1
cos ln sin , 1, , ,
2
N
i j
i N ij ij ij
j ij
r
A i N
r
|
t
+
+
=
( | |
= A A = ( |
|
(
\ .
( ) ( )
( ) ( )
1, 1 1, ,
1, 1 , 1
1
1, ,
1
sin sin
2
cos ln cos ln ,
N j j j Nj N j
j N j
j Nj
j N j
A
r r
r r
| |
t
+
+ +
= A + A
( | | | |
A A ( | |
| |
(
\ . \ .
( ) ( )
( ) ( )
1, 1 , 1
1, 1 1
1 1, ,
1 1
1
sin ln sin ln
2
cos cos ,
N
j N j
N N j Nj
j j N j
j j Nj Nj
r r
A
r r t
| |
+ +
+ +
=
| | | |
= A + A | |
| |
\ . \ .
(
+ A + A
(
44
( )
1 1
sin , 1, , ,
cos( ) sin( ),
i i
N N
b V i N
b V V
u o
u o u o
+
= =
=
and .
ij i j
u u A =
The linear system (30) will give the searched values: and q
j
, where j=1,, N, that will help further
to calculate the tangential components of velocity in the points of collocation M
i ,
i = 1,, N. Let
remind that the normal component of velocity in the points of collocation is null. The below relation
gives the tangential component:
cos sin .
i i i i i
u u v
t
u u = +
Let substitute the relation (28) in the above relation to obtain:
1 1 1 1
cos cos sin sin .
N N N N
s v s v
i j ij ij i j ij ij i
j j j i
u V q u u V q v v
t
o u o u
= = = =
| | | |
= + + + + +
| |
\ . \ .
Consequently, the following relations will be obtained for the tangential components of velocity:
( ) ( ) ( )
( ) ( )
, 1
1
, 1
1
cos sin cos ln
2
sin ln cos .
2
N
i j
i
i i ij ij ij
j ij
N
i j
ij ij ij
j ij
r
q
u V
r
r
r
t
u o |
t
|
t
+
=
+
=
( | |
= + A A ( |
|
(
\ .
( | |
+ A A ( |
|
(
\ .
(31)
Bernoulli equation
2 2
1 1
2 2
p V p V
+ = +
implies that
2 2
1 1
.
2 2
p p V V
= +
Thus, the local coefficient of pressure can be rewritten as follows:
2
2
2
1 .
1
2
p
p p V
C
V
V
= =
(32)
Accordingly, the local pressure coefficient on the discretized contour profile can be calculated from
the relation
2
,
1 ,
i
p i
u
C
V
t
| |
=
|
\ .
(33)
where components
i
u
t
are supplied by formula (31).
The hydrodynamic forces that react on the boundary element E
j
are obtained from the relations
relations:
( )
( )
, 1
, 1
,
,
xj p j j j
yj p j j j
f C y y
f C x x
+
+
=
=
(34)
and the pitching moments reported to the point of reference
( )
( )
, ,0 ,
4
ref ref
c
x y = are calculated by
formula:
1 1
,
.
2 2 4
j j j j
m j xj yj
y y x x
c
c f f
+ +
| | | |
= +
| |
\ . \ .
(35)
45
The total force is the sum of contributions of each boundary element:
1
1
,
,
N
x xj
j
N
y yj
j
F f
F f
=
=
=
=
(36)
and the lift coefficient and the moment coefficient are calculated as follows:
sin cos ,
L x y
C F F o o = + (37)
,
1
.
N
M m j
j
C c
=
=
(38)
5.1.3.3. Selection of the optimal hydrodynamic profile of blades
The optimization of the hydrodynamic blade turbine performance demands blade optimal
hydrodynamic profile. The numerical calculation methods, previously described, are used to
calculate the coefficients
, L ref
C
and
, D ref
C
for the symmetrical profiles from the library of NACA
aerodynamic profiles with a chord length c
ref
= 1 m. It should be remarked that the calculation
method converges for the entering angles that do not exceed 20 25
o o
dependent on the
selected profile and the corresponding Reynolds number (Re = 1300000). For the entering angles
exceeding this critical value, the rates corresponding to a flat (plane) profile are considered. Some
of the considered profiles are shown in figure 35: NACA 0012, 0016, 63018 and 67015. Figure 36
shows the hydrodynamic lift
, L ref
C and drag
, D ref
C coefficients depending on the entering angle.
Taking into account the data from Fig. 36, the NACA 0016 hydrodynamic profile is being selected
as the reference profile. Subsequently, this profile will be optimized in order to increase the turbine
performance.
0.2 0.4 0.6 0.8 1
-0.3
-0.2
-0.1
0
0.1
0
0.2
0.3
Profil NACA 0012
X (m)
Y
(
m
)
0.2 0.4 0.6 0.8 1
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
Profil NACA 0016
X (m)
Y
(
m
)
0.2 0.4 0.6 0.8
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
Profil NACA 63018
X (m)
Y
(
m
)
0.2 0.4 0.6 0.8
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
Profil NACA 67015
X (m)
Y
(
m
)
Figure 34 Symmetric hydrodynamic profiles: NACA 0012, 0016, 63018 and 67015.
46
15 30 45 60 75 90
0
0.5
1
1.5
2
Profil NACA 67015
Unghiul de atac, (Deg)
C
L
,
C
D
15 30 45 60 75 90
0
0.5
1
1.5
2
Profil: NACA 63018
Unghiul de atac, (Deg)
C
L
,
C
D
15 30 45 60 75 90
0
0.5
1
1.5
2
Profil: NACA 0016
Unghiul de atac, (Deg)
C
L
,
C
D
0 15 30 45 60 75 90
0
0.5
1
1.5
2
Profil: NACA 0012
Unghiul de atac, (Deg)
C
L
,
C
D
C
L
C
D
Figure 35 Hydrodynamic lift
L
C and drag
D
C coefficients dependant on the entering angle for
NACA 0012, 0016, 63018 and 67015 profiles.
5.1.3.4. The torque moment and the forces applied on the multi-blade
hydrodynamic rotor
The hydrodynamic coefficients for the NACA 0016 reference profile with chord length, for instance,
1,3 . c m = , are calculated below. The coefficients corresponding to the profile with the chord length 1.3 m
are calculated from the relations:
47
,
2
,
,
1,3 ,
(1,3) ,
1,3 .
L L ref
M M ref
D D ref
C C
C C
C C
=
=
=
(39)
The values of the lift and drag coefficients
dependant on the entering angle are
shown in figure 37. Taking into account
these values, the angle 18
o
o = is selected
as the working entering angle.
The blade changes its entering angle
during its motion depending on the position
(figure 38). Thus, in sector I the entering
angle (angle formed by the blade and water
flow) is 18; in sector II the entering angle
shifts from 18 up to -18, but the blade
does not contribute to the total moment
developed at the rotor shaft. In this sector,
extended up to approximately 60, the
blade is carried freely by the water flow and its re-positioning takes place at an angle of -18 at the
end of sector III. The entering angle is -18 in sector III. In sectors IV-VI the hydrodynamic effect is
minimal and the blade has to be re-positioned from angle -18 to angle 18. In order to use the
kinetic energy in the sectors IV-VI it is proposed to re-position the blade from -18 to 90 in sector
IV; in sector V the blade remains under an angle of 90, and in sector VI the entering angle returns
to 18. Knowing the values of the hydrodynamic coefficients
L
C
and
D
C
, the lift force
L
F
and drag
force
D
F
are calculated by the formulas (10) and (11), and the formula (12) supplies the
hydrodynamic force that reacts on the blade (figure 39).
0 10 20 30 40 50 60 70 80 90
0
0.25
0.5
0.75
1
1.25
1.5
1.75
2
2.25
2.5
Unghiul de atac, (Deg)
C
L
,
C
D
Coeficientii C
L
si C
D
in functie de unghiul de atac. Profil: NACA 0016
Lungimea palei1.3 m
Numarul Reynolds1300000
C
L
C
D
Figure 36 - Hydrodynamic lift
L
C
and drag
D
C
coefficients dependant on the entering angle for
NACA 0016 profile.
0 30 60 90 120 150 180 210 240 270 300 330 360
4000
3500
3000
2500
2000
1500
1000
500
0
500
1000
1500
2000
2500
3000
3500
4000
Fortele care actioneaza pe pala , Profil:NACA 0016
Unghiul de pozitionare, (Deg)
F
o
r
t
e
l
e
,
(
N
)
Profil: NACA 0016
Viteza fluxului de apa = 1 m/s
Unghiul de atac = 18 Deg
Raza rotorului = 2 m
Numarul palelor = 5
Inaltimea palei = 1.4 m
Lungimea palei = 1.3 m
Modulul rezultantei
Componenta tangentiala
Componenta normala
Figure 37 Blade position and working
areas.
Figure 38 Module, tangential component and normal
component of the hydrodynamic force of a rotor
blade depending on the angle of positioning.
48
The module of the hydrodynamic force F
m
)
Profil: NACA 0016
Viteza fluxului de apa = 1 m/s
Unghiul de atac = 18 Deg
Raza rotorului = 2 m
Numarul palelor = 5
Inaltimea palei = 1.4 m
Lungimea palei = 1.3 m
0 30 60 90 120 150 180 210 240 270 300 330 360
0
2000
4000
6000
8000
10000
12000
14000
Momentul total la arborele rotorului functie de unghiul de pozitionare
Unghiul de pozitionare (Deg)
M
o
m
e
n
t
,
(
N
m
)
Profil: NACA 0016
Viteza fluxului de apa = 1 m/s
Unghiul de atac = 18 Deg
Raza rotorului = 2 m
Numarul palelor = 5
Inaltimea palei = 1.4 m
Lungimea palei = 1.3 m
Figure 39 Moment
, r i
T
developed by the rotor
blade depending on the angle of positioning.
Figure 40 Total moment
r
T
E
developed by 5
blades at rotor shaft depending on the angle of
positioning.
0 30 60 90 120 150 180 210 240 270 300 330 360
0
0.5
1
1.5
2
2.5
3
3.5
4
4.5
x 10
4
Unghiul de pozitionare, (Deg)
M
o
m
e
n
t
,
(
N
m
)
Momentul total la diferite viteze
1 m/s
1.3 m/s
1.6 m/s
Figure 41 Total moment
r
T
E
at rotor shaft
depending on the angle of positioning for
various velocities of the water flow
Figure 42 Number of turns
, M ref
C depending on
the entering angle for NACA 0016 profile
49
Taking into account the fact that the hydrodynamic force is not applied in the blade fixed coordinate
system O' (figure 44) this force produces a moment of torsion called the pitching moment. This
moment is determined given the point of reference. Point P will be considered as the point of
reference situated at distance of the chord from the entering edge B (figure 44). For the working
values of the entering angle = 18 it is obtained C
M,ref
= -0.026. Thus, from the relation (40)
results that C
M
= 0.0439. The moment of torsion compared to the point P is
2
1
39,92 ,
2
M p
M C V cS N m
= =
(40)
where V
= 1 m/s, c = 1.3 m and H = 1.4 m. In the system of coordinates O x y ' '' '' , the components
of the hydrodynamic forces are delivered by the relation (12). Applying the values F
L
and F
D
obtained previously we have:
1601,2 ,
413,8 .
x
y
F N
F N
''
''
=
=
(41)
Then
0,0249 25 .
x
O P M F m mm
''
' = = ~ (42)
In order to ensure the stability of the blade motion, the fixing point W should be selected in the
interval
25 mm O W H ' s s
,
where
min max
H H H s s . Values
H
min
and H
max
are taken under
the condition that the frictional
force, appearing in the
kinematical couples of the
orientation mechanism, must
be minimal.
To determine the optimal working entering angle it is necessary to calculate the value of the
moment developed by one blade and the total moment for several values of the entering angle,
namely: 15 , 17 , 18 , 20 ,
o o o o
o = (figure 45-46). In this context the entering angle for the blade with
hydrodynamic profile NACA 0016 is 18.
o
o =
0 30 60 90 120 150 180 210 240 270 300 330
0
1000
2000
3000
4000
5000
Momentul dezvoltat de o pala la diferite unghiuri de atac
Unghiul de pozitionare, (Deg)
M
o
m
e
n
t
,
(
N
m
)
15 Deg
17 Deg
18 Deg
20 Deg
0 30 60 90 120 150 180 210 240 270 300 330 360
0.5
0.6
0.7
0.8
0.9
1
1.1
1.2
1.3
1.4
1.5
1.6
1.7
x 10
4
Momentul total la diferite unghiuri de atac
Unghiul de pozitionare, (Deg)
M
o
m
e
n
t
,
(
N
m
)
15 Deg
17 Deg
18 Deg
20 Deg
Figure 44 Moment developed by the blade
, r i
T
depending on the positioning angle for various
values of the entering angle 15 , 17 , 18 , 20.
o o o o
o =
Figure 45 Total moment
r
T
E
depending on the
positioning angle for various values of the
entering angle 15 , 17 , 18 , 20.
o o o o
o =
Figure 43 Location of the blade fixing point.
50
Also, the performance of 3-, 4- and 5-blades rotor was analysed. The total moment developed by
the rotor shaft was calculated and the results are presented in figure 47.
0 30 60 90 120 150 180 210 240 270 300 330 360
0.5
0.6
0.7
0.8
0.9
1
1.1
1.2
1.3
1.4
1.5
1.6
1.7
x 10
4
Momentul total la diferite configuratii
Unghiul de pozitionare, (Deg)
M
o
m
e
n
t
,
(
N
m
)
3 pale
4 pale
5 pale
Figure 46 Total moment
r
T
E
developed at the 3-, 4- and 5-blade rotor shaft depending on the
positioning angle.
5.1.3.5. Optimisation of NACA 0016 hydrodynamic profile
In order to maximize the moment of torsion produced by the micro hydro power plant rotor, the
optimization of the hydrodynamic profile will be considered. The moment of torsion depends on the
lift and drag hydrodynamic forces given by formulas (10) and (11). The hydrodynamic forces
through the hydrodynamic coefficients depend on the entering angle , Re number and the shape
of the hydrodynamic profile. The hydrodynamic shape of the profile was selected from the NACA
library of 4 and 5 figures having as parameters (with account of the profile symmetry) only the
maximal thickness. The entering angle constitutes the second parameter. The optimization aims at
maximizing the lift force and, at the same time, does not allow the pitching moment and the
resistance force to take very big values. The following issue of optimization should be considered:
Maximize ( , )
L L
C C u o =
with constraints imposed to the coefficients
D
C and
M
C , (43)
where is the maximum thickness and is the entering angle. The values of the inferior and
superior borders are determined, as follows: the negative maximum value for the pitching
coefficient will correspond to the solution for the entering angle 0. The maximum value for the
resistance coefficient will correspond to the solution for the entering angle = 18. Also, restrictions
have been added to the optimization parameters 10% 20% u s s and 0 20
o o
o s s . To find the
optimal values of function
1
( , , )
n
f f x x = an iterative method is used:
As long as the demanded accuracy is not reached the solution will be,
( )
i i i
B s f x = V ,
1 i i i i
x x s o
+
= + , (44)
51
where
i
o are the multipliers and
i
B are the definite positive approximations of the Hessian
function f . The partial derivation of function f related to the component i is approximated with
the help of the finite difference formulas:
( ) ( )
( )
2
i i
i
f x he f x he f
x
x h
+ c
=
c
, (45)
where
i
e is the basis vector.
The optimization is done by the MATLAB optimization soft: Sequential quadratic
programming algorithm with a line search and a BFGS Hessian update. The quadratic sub-tasks
are solved by modified projection method. The gradient of function ( , )
L L
C C u o = is calculated by the
finite difference formulas with the constant pitch
4
10 h
= . NACA 0016 profile was considered as
the initial profile (figure 48). The result of optimization is presented in Fig. 49. The results of the
carried out research were published in [1 - 5].
Figure 47 NACA 0016 hydrodynamic rack profile
standard.
Figure 48 NACA 0016 hydrodynamic rack
profile standard and the optimised profile.
In order to optimize the hydrodynamic
profiles of the blades a prototyping 5-axis
machine has been purchased (figure 50).
For the manufacturing of the prototypes
of hydrodynamic blades in the framework
of the current project several special
modules are at purchasing stage. The
purchased equipment will allow the
prototyping, manufacturing and testing of
the profiles sugested by modelling and
computer simulations.
Figure 49 Blades prototyping 5-axis machine
52
5.1.3.6. Floating stability of the micro hydropower plant
The micro hydro power plant is posted in the river water flow. The position of blades compared to
the water level is ensured by the Archimedes forces that react on the floating blades. The blade
cavity generates the Archimedes force determined by
the relation
,
A
F Vg = (46),
where is the water density, V is the interior volume of
the blade and g is the gravitational acceleration. The
analysis of the path of motion (motion trajectory) of the
points of application of the Archimedes force F
A
(points
N
i
, i = 1, 2, 3, figure 51) has shown that the distance
from these points to the rotor axis O will oscillate
depending on the positioning angle . Thus, for the
blades located in the superior semi-plane defined by the
straight line OO these distances are different from the
respective distances of the blades located in the inferior
semi-plane.
This fact leads to the appearance of the pitching
moment with respect to the axis of longitudinal
symmetry of the floating bodies:
, ,
,
r S I
M M M
E E
= (47)
where
,S
M
E
is the total moment developed by the Archimedes forces that react on the blades,
currently located in the superior semi-plane, and
,I
M
E
is the total moment developed by the
Archimedes forces that react on the blades currently located in the inferior semi-plane.
The total moments developed by the Archimedes forces that react on the blades currently located
in the superior and the inferior semi-plane are determined by the relation
, , ,
,
S A i A i
M F D
E
=
(48)
where
, A i
F are the Archimedes forces that react on the blades;
, A i
D are the distances from the
point of application of the Archimedes forces to the rotor spindle, and the summarization is done for
all blades located in the superior semi-plane. Similarly,
, , ,
.
I A i A i
M F D
E
=
(49)
Distances
, A i
D are calculated by the formula:
2 2 2
,
2 cos( ),
A i A A i
D R c Rc o = + + + (50)
where R is the rotor radius;
A
c is the distance between the point of application of the Archimedes
forces and the point of blade fixing to the rotor lever; o is the angle formed by the blade chord
Figure 50 Floating stability analysis.
53
AB and the direction of the water flowing, and
i
is the angle formed by the rotor lever and the
direction OO' .
To compensate the pitching moment
r
M it is proposed to locate the rotor spindle in the plane
shifted at distance e compared to the plane of longitudinal symmetry of the floating bodies.
Distance e is calculated by the relation:
1
,
pal
N
i
i
pal
y
e
N
=
=
(51)
where
pal
N is the number of rotor blades and
i
y
is the distance from the central point of
application of the Archimedes force on the blade i till the plane of longitudinal symmetry (figure
51). For each blade, distance
i
y
is calculated by the relation:
360
cos sin( ( 1) ).
o
i A
y c R i
n
o = + +
(52)
Let introduce (52) into (52) and obtain:
cos
A
e c o = (53)
The point of application of the Archimedes force on each blade is the centre of gravity (mass point)
of the applied hydrodynamic profile, in our case NACA 0016 profile. The central point of application
of the Archimedes force system that reacts on a number
pal
N
of the submersible blades will
describe a migration trajectory generated by the rotor revolution. The migration trajectory
generated by a complete revolution of 3- and 5-blade rotor represents closed curves described in
figure 52 (a) and figure 52 (b). One point on the closed curve represents the position of the central
point of application of the Archimedes force system corresponding to a concrete angular position of
the rotor.
To identify the technical solution ensuring the floating stability of the micro hydro power plant it is
0.03 0.01 0.01 0.03 0.05 0.07 0.09 0.11 0.13
0.22
0.24
0.26
0.28
0.3
0.32
0.34
0.36
0.38
X, (m)
Y
,
(
m
)
Traiectoria de migrare a punctului c
A
0.01 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07
0.29
0.3
0.31
0.32
0.33
0.34
0.35
0.36
0.37
X, (m)
Y
,
(
m
)
Traiectoria de migrare a punctului c
A
necessary to estimate
the values of the
distance between the
central point of
application of the
Archimedes force
system and the
longitudinal symmetry
axis of the floating
Figure 51 Migration trajectory of the central point of application of the
Archimedes forces for the 3-blade (a) and 5-blade rotor (b).
bodies Figure 53 (a, b).
shows
54
0 30 60 90 120 150 180 210 240 270 300 330 36
0.24
0.26
0.28
0.3
0.32
0.34
0.36
Unghiul de pozitionare, (Deg)
e
,
(
m
)
Distanta e
0.3296
0 30 60 90 120 150 180 210 240 270 300 330
0.29
0.3
0.31
0.32
0.33
0.34
0.35
0.36
Unghiul de pozitionare, (Deg)
e
,
(
m
)
Distanta e
0.3296
the distance e
depending on the
positioning angle of
the 3-blade (a) and 5-
blade (b) rotor. It has
been stated that in the
case of 3-blade rotor,
distance e takes values
comprising
0,238
min
e m = and
max
0,363 e m = .
Figure 52 Dependence of distance e of the central point of application of
the Archimedes forces on the positioning angle of the 3-blade rotor (a)
and of 5-blade rotor (b).
And in the case of 5-blade rotor: 0,289
min
e m = , and
max
0,363 e m = . Let calculate the mean value
of distance e depending on the positioning angle : for 3- and 5-blade rotors the same mean
distance 0,33
med
e m = is obtained.
Conclusions:
1. To ensure the floating stability of the micro hydro power plant it is necessary that the spindle of
3- and 5-blade rotor shifts from the axis of longitudinal symmetry of the floating bodies at rate
0,33
med
e m = in the direction opposite to the water flow.
2. Micro hydro power plants anchored on the left bank differ from those anchored on the right bank
by the spatial truss constructions and, in particular, by the constructional elements of the
hydrodynamic rotor shifted at rate 0,33 .
med
e m =
5.1.3.7. Turbulence and stability of the hydrodynamic rotor
Figure 54 shows the field of fluid velocities around NACA 0016 profile at entering angle 18 and the
Reynolds number calculated from the relation [2]:
Re ,
cV cV
q v
= =
(54)
where the fluid density is
3
998,4 kg m = at 20
o
C,
the kinematical viscosity
is
6 2
1,012 10 m s v
= ,
and the length of the
profile chord is
1,3 c m = .
For the fluid flow rate
1 , 1,5 , 2 V m s m s m s
=
the following values of the
Figure 53 Velocity field around NACA 0016 profile at the entering
angle18
o
.
55
Reynolds number Re = 1284600, 1798400 and 2312300 are obtained. The transition and
separation of the boundary layer on the inferior and superior surfaces of the blade profile is shown
in figure 55. Points T.U. and T.L. stand for the points of separation from the laminar flow to the
turbulent flow on the inferior C
inf
and superior C
sup
surfaces of the blade (figure 52). Respectively,
points S.U. and S.L. represent the points of separation on the inferior and superior surface. In all
cases it was observed that transition from the laminar flow to the turbulent flow takes place in the
proximity of the point of stagnation, and the separation of flow from the profile surface is foreseen
at an approximate 40-50% distance of the blade chord length. The transition from the laminar flow
to the turbulent flow as well as the separation of the turbulent boundary layer will take place in the
proximity of the rear edge on the inferior surface.
5.1.3.8. Design and manufacturing of the hydrodynamic rotor
The hydrodynamic rotor has been designed in the Autodesk MotionInventor software (figure 56 - 3
blades rotor, and figure 57 - 5 blades rotor). The
hydrodynamic rotor is the main working
element of the micro hydropower plant, which
converts kinetic energy of the water flow and
transmit it via the kinematical linkage to the
production units of electrical (generator 1)
energy or mechanical (hydraulic pump 2)
energy (figure 58).
From the point of view of its design, the
rotor comprises the main shaft 1 (figures 56,
57), the casing with radial bars 2, on which
ends the hydrodynamic profile blades 3 are
mounted with the help of node 4. The main
shaft 1 and the casing with the bars 2 are
mounted removable. The hydrodynamic rotor
is a spatial structure, strained complex with
the bending and twisting moments. The
casing with radial bars are made of aluminium
alloy plates with calculated dimensions able to
ensure design positioning (calculated) of
blades with minimal deviations (the deflexion
of blade axles till 5 mm, the angle of twist of
the radial bars 1
o
).
a. b.
Figure 54 Point of separation for the flow velocities 1 m/s (a) and 2 m/s (b).
Figure 55 3-blade hydrodynamic rotor.
1 2 4
3
Figure 56 5-blade hydrodynamic rotor.
2 1
4
3
56
The technical solutions adopted in the final
construction of the hydrodynamic rotor have
resulted into computer simulated research, by
utilising ANSYS CFX 5.7, sub programmes in the
package of mathematic modelling MathCAD, etc.,
by applying possible loading in real conditions.
The node 4 allows the changing of blade 3
positioning towards their axis of rotation, aiming
at the insurance of the optimal pitching moment
(the technical solution contains Know-How
elements). Figure 59 shows 3-blade (a) and,
respectively, 5-blade (b) rotors, the diameter of
blade location is D = 4 m. The shape and
dimensions of the hydrodynamic profile of blades
have been justified within the carried out research
[1,2].
Figure 57 Multiblade rotor connected
kinematically to the electric (generator 1) energy
or mechanical hydraulic pump 2) energy
production units.
Pilot station destination. The micro hydro power plant is a complex system, which includes the
main working element the hydrodynamic rotor, destined for the conversion of a part of the river
water kinetic energy potential into useful energy, and the units with distinct functions, such as the
multiplier, generator, hydraulic pump etc. The results of research on the construction and functional
parameters of the mentioned units, shown separately from each other regardless of their functional
interdependence as an integral system, are presented in [1,2].
To carry out research on the functional parameters, in field conditions, of all units participating in
the process of conversion of flowing water kinetic energy into useful energy as an integral technical
system, the pilot plant design was developed in MotionInventor AutoDesk design environment,
Figure 58 Rotors with 3- (a) and 5-blades (b) with hydrodynamic profile, manufactured in the
laboratory of the Centre for Renewable Energy Conversion Systems Design, TUM.
using the developed conceptual diagrams [1,2] and the results of theoretical research findings, of
computer modelling and simulation of hydrodynamic profile blade interaction with the fluid (water).
The pilot station is intended for research and verification of the technical solutions adopted at each
stage of unit development, in real conditions, and, if necessary, changes will be made at the design
stage of the industrial prototypes of micro hydro power plants.
Pilot station positioning. The pilot station is installed on the river Prut, v. Stoieneti, district
Cantemir, under the following positioning conditions:
- Rotor axis is perpendicular to the water mirror;
57
- Longitudinal axis of spatial housing is perpendicular to the water flow velocity vector;
- Level of blades submersion meets the project rate (h = 1.4 m). Perpendicularity of the
hydrodynamic rotor shaft to the water mirror is ensured by changing the buoyancy of four floating
bodies, and the perpendicularity of the longitudinal axis of the spatial housing to the water flow
stream vector is ensured through supporting cables secured to the anchoring rods. The
submersion level of blades with hydrodynamic profile h = 1.4 m is maintained by changing the
buoyancy of the four equal floating bodies.
Structural and functional parameters of hydrodynamic rotor, planetary multiplier, generator and low
speed centrifugal pumps, determined separately for each working body, needs to be checked by
experimental investigations in real conditions of their operation in an integral technical system.
In this context, the time table of experimental research on the pilot station in field conditions
includes:
1. study of diversity of water flow speed cadastre within the boundaries of rotor effective
section (the width of the rotor blades and level of submersion) and assessment of the water flow
energy potential;
2. study of the influence of force factors on the stability of the hydrodynamic profile blades
positioning (angle of attack o) and the kinetostatic analysis of the blades positioning mechanism;
3. study of energy and kinetic conversion efficiency parameters (for the electric generator
terminals and the input shaft of hydraulic pumps);
4. study of kinematic parameters of hydrodynamic rotor and mechanical losses in the
kinematic chain (linkage);
5. setting the influence of structural and functional parameters of hydrodynamic rotor on the
hydrodynamic effects and water turbulence flow mode under field conditions;
6. study of functional characteristics of electrical generator and centrifugal pumps.
The experimental research of the integral technical system - hydrodynamic rotor coupled
kinematically with component units of micro hydro power plants aims at increasing efficiency of
water flow kinetic energy conversion into useful energy by identifying and, where necessary,
introducing in the technical documentation of partial structural changes and, sometimes, of
conceptual and technical solutions adopted previously.
When developing industrial prototypes of micro hydropower plants for river water kinetic energy
conversion the following criteria and requirements were taken into account:
- Exclusion of dam construction and of the negative impact on the environment, implicitly;
- Lowest cost;
- Simplicity of construction and operation;
- Increased reliability at dynamic overload in operating conditions;
- Resistant composite materials, including conditions of high humidity;
- Automatic adjustment of the micro hydropower plant platform position in conditions of water
level changing.
The technical solutions adopted in the process of micro hydropower plants design result from
theoretical and experimental research presented in [1,2]. To justify the constructive and functional
parameters, additional numerical modelling and simulations were performed, using ANSYS
software CFX5.7, and the sub-software developed by the authors for MathCAD, Autodesk
MotionInventor, etc., namely, simulation of the fluid blade interaction and floating stability,
hydrodynamic optimization of the blade profile in order to increase kinetic energy conversion
efficiency of river water at its different velocities using 3-, 4- and 5-blades rotor. The efficient
operation of micro hydropower plants by individual customers for particular destination depends on
their constructive configuration choices and on the functional characteristics of component units
58
participating in the conversion of flowing water kinetic energy into useful energy.
To meet the objectives and consumer demands for micro hydro power plants, and also to increase
the efficiency of conversion of flowing water kinetic potential in the certain river area, the authors
have developed the following structural and functional concepts based on modular assembly:
1. micro hydropower plant with
hydrodynamic rotor for river water kinetic energy
conversion into mechanical energy for water
pumping (MHCF D4x1,5 M);
2. micro hydropower plant with
hydrodynamic rotor for river water kinetic energy
conversion into electric and mechanical energy
(MHCF D4x1,5 ME);
3. micro hydropower plant with
hydrodynamic rotor for river water kinetic energy
conversion into mechanical energy at low
speeds (MHCF D4x1,5 ME);
4. micro hydropower plant with
hydrodynamic rotor for river water kinetic energy
conversion into electric energy (MHCF D4x1,5
E).
The micro hydropower plants designed
modularly, allow change of destination and of
functional characteristics by replacing some units
with others (generator, pump, blades with other
hydrodynamic profile, 3- and 5- blade rotor).
The micro hydropower plants possess similar resistance in structure construction calculated for
strength and rigidity to dynamic applications. Buoyancy and maintainance of hydropower plant
rotor axis perpendicularity at variable river water level is ensured by patented technical solutions
[8-14]. Continuous orientation mechanism of the blades at the constant entering angle relative to
the direction of fluid stream contains Know-how elements and is not described here. The main
working body, which depends mainly on the amount of kinetic energy converted into useful energy,
is the blade with hydrodynamic NACA 0016 profile, developed on the basis of carried research
[1,2].
Two types of 3- and 5-blade rotors were designed for the described micro hydropower plants.
Installed capacity of micro hydropower plants with diameter D = 4 m, submerged height of blades h
= 1,4 m and the length of blade chord l = 1,3 m at water flow velocity V = 1...2 m/s can be in the
boundaries P = 2...19 kW (see [figure 60]).
5.2. Industrial prototypes of micro hydropower plant with
hydrodynamic rotor
5.2.1. Pilot station of micro hydropower plant with hydrodynamic rotor
for river water kinetic energy conversion into mechanical energy (MHCF
D4x1,5 M)
This model of hydro power plant is designed to convert river water kinetic energy into mechanical
energy used to pump water into irrigation and sewerage systems, and supply industrial water etc.
with the flow rate Q = 40 m
3
/h at the height pumping H = (10 15) m.
Figure 59 Generated power at rotor shaft.
59
Static description of the microhydro power plant. The blades 1 (figure 61) are connected to the
hydrodynamic rotor 2 by roller friction bearings to ensure their orientation under a certain entering
angle o.
Figure 60 Micro hydropower plant with hydrodynamic rotor for river kinetic energy conversion into
mechanical energy for water pumping (flow rate Q = 40m
3
/h, pumping height H =10...15 m)
(MHCF D4x1,5 M)
1.1. Blade with hydrodynamic 0016
NACA profile; 2 3-blade rotor; 3
planetary multiplier with multiplying
ratio i=112; 4 belt drive with
multiplying ratio i = 1,9; 5 - permanent
magnet generator (characteristics
see p. 5.4); 6. centrifugal pump
PSS4010/50 (characteristics pump
flow rate Q=40m
3
/h at pumping height
(10...15)m; 7 plastic pontoons; 8
guide; 9 spatial housing.
60
Functioning principle. The river flowing water with the energy potential dependent on the flow velocity
drives the hydrodynamic profile blades 1 (figure 62), oriented continuously by the entering The hydrodynamic
rotor 2 is mounted on the input shaft of the planetary multiplier 3 through an auxiliary shaft, which is fixed on
the bearings. The belt pulleys of the transmission 4 are mounted on the output shaft of the planetary
multiplier - the big one, and the small one - on the input shaft of the centrifugal pump 5. The hydrodynamic
rotor 2 and blades 1, the multiplier 3, the centrifugal pump 5 and guides 6 are mounted on the spatial
housing 7, installed on the
pontoons 8.
The micro hydropower
rotor 2 comprises three
blades oriented at an
entering angle o, which is
dependent on the water
flow velocity. In the areas
of blades 1 location,
inefficient from the point of
view of river kinetic energy
conversion, under
hydrodynamic forces the
blades 1 are repositioned
at an angle of 90
o
to the
currents of water or are
carried by the water
unhampered to the angle
o = 0. Thus, the respective
positioning of blades allows
the increase of water
kinetic energy rate
converted into useful
energy. As result, the water
currents transmit a part of
their kinetic energy to the
blades 1, stressing them
under the hydrodynamic forces and
reporting rotational motion with angular
frequency e
1
and torque T
1
to the rotor 2.
The summary torque T
1
, developed by the
hydrodynamic forces and applied to the 3-
blade rotor shaft at water flow velocities
1.3, 1.6 and 1.8 m/s and at the entering
angle of blades o = 18
o
, is presented in
figure 63.
For rotor diameter D = 4 m, the submerged
height of blades h=1,4 m and length of
blade chord l = 1,3 m, the torque is:
T
1
= 11938 Nm for water flow velocity V =
1,3 m/s; T
1
= 18084 Nm for V = 1,3 m/s;
T
1
= 22887 Nm for V = 1,8 m/s. The
calculations of kinematics and lifting
capacity of all constructive elements as well
as of all functional and energy parameters
Figure 61 Kinematics of micro hydropower plant MHCF D4x1,5 M.
0 30 60 90 120 150 180 210 240 270 300 330 360
0
0.5
1
1.5
2
2.5
3
3.5
4
x 10
4
Unghiul de pozitionare, (Deg)
T
,
(
n
m
)
Momentul sumar T la diferite viteze, Profil:NACA 0016
v
0
=1.3 m/s
v
0
=1.6 m/s
v
0
=1.8 m/s
22887 N m
18084 N m
11938 N m
valoare medie
Figure 62 Torque T
1
at the hydrodynamic rotor shaft with
NACA 0016 profile blades.
61
of micro hydropower plant have been carried out for the torque value T
1
= 18084 Nm.
Rotor 2, rigidly coupled by means of auxiliary shaft with the input shaft of the multiplier 3, transmits
rotational motion to the last with angular frequency e
1
and torque T
1
. The multiplier reproduces the
rotor 2 revolutions up to
1 1
2 1
30
n i (min ),
e
t
=
where i
1
represents the multiplying ratio of the multiplier
(i
1
=112). Rotational motion at angular frequency
1 2
2
n
( s )
30
t
e
=
is transmitted from the multiplier
input shaft via a transmission belt 4 of the centrifugal pump input shaft with multiplying ratio i
1
=
2,25. As result, the input shaft of the centrifugal pump swivels with angular frequency e
3
= e
1
i
1
i
2
(s
-1
) and is stressed at torque:
1 1 2 r
3
1 2
T
T ,( Nm)
i i
q q q
=
,
where: q
1
is the multiplier mechanical efficiency (q
1
= 0,9);
q
2
- is belt transmission mechanical efficiency (q
1
= 0,95);
q
r
- mechanical efficiency of hydrodynamic rotor bearings (q
1
= 0,99).
According to the experimental research presented in [1,2] the mechanical efficiency of centrifugal
pump is q
1
= 0,72 at rated speed frequency
1 3
3
30
n 500 min .
e
t
= =
The mechanical efficiency of the micro hydropower plant with hydrodynamic rotor for river water
kinetic energy conversion into mechanical energy, with account of all mechanical losses in the
linkage (at the hydraulic pump shaft) is:
1 2 3 r
0,9 0,95 0,99 0,846. q q q q q
= = =
Accordingly, the micro hydropower plant (MHCF D4x1,5 M) ensures the transformation into useful
energy of 84,6% of the kinetic energy potential of the flowing water transmitted to the
hydrodynamic rotor.
62
Figure 63 Industrial prototype of the microhydropower station for the river kinetic energy
conversion into electrical and mechanical energies (diameter of rotor d = 4m, submersed height of
the blade h = 1,4m, length of blade l =1,3m) (MHCF D4x1,5 ME).
According with obtained results the industrial prototype of the microhydropower station for the river
kinetic energy conversion into electrical and mechanical energies (diameter of rotor d = 4m,
submersed height of the blade h = 1,4m, length of blade l =1,3m) (MHCF D4x1,5 ME) was
produced (figure 64). Now is installed on the river Prut, v. Stoieneti, Cantemir for testing in real
conditions (figure 65).
Figure 64 Industrial prototype of the microhydropower station for the river kinetic energy
conversion into mechanical energy installed on the river Prut, v. Stoieneti, Cantemir.
5.2.2. Micro hydropower plant with hydrodynamic rotor for river water
kinetic energy conversion into electrical and mechanical energy (MHCF
D4x1,5ME)
The micro hydropower plant MHCF D4x1,5 ME for river water kinetic energy conversion into
electrical and mechanical energy (figure 66) is polyfunctional and can be utilised for electrical
lighting of streets, heating, water pumping in drip irrigation systems, and also for draining
agricultural lands adjacent to rivers.
Rigging the NACA 0016 profile blades 1 in the hydrodynamic rotor 2 and its attachment to the
multiplier input shaft 3 is done similar to micro hydropower plant MHCF D4x1, 5 M. Kinematic and
construction peculiarities of MHCF D4x1,5 ME are as follows: rotational movement of the
63
hydrodynamic rotor 2 (Figure 67) with an angular frequency (velocity) e
1,
by means of multiplier 3
and belt transmission 4 with effective multiplying ratio i = 212.8, is multiplied to the operating
angular frequency of the permanent magnet low speed generator 5:
e
3
=e
1
i
1
(s
-1
).
Torque T
3
, applied to rotor 5, is:
1 1 2 r
3
T
T ,( Nm)
i
q q q
= ,
where: q
1
is multiplier mechanical efficiency (q
1
= 0,9);
q
2
is belt transmission mechanical efficiency (q
1
= 0,95);
q
r
mechanical efficiency of hydrodynamic rotor bearings (q
1
= 0,99).
i effective multiplying ratio equal to the multiplying ratios product of the planetary multiplier
and belt transmission
64
Figure 65 Micro hydropower plant with hydrodynamic rotor for river kinetic energy conversion into
electrical and mechanical energy (rotor diameter D = 4 m, submerged height of blade h = 1,4 m,
length of blade chord l = 1,3 m) (MHCF D4x1,5 ME)
The electric energy produced by the permanent magnet generator 5 (figure 67) can be utilised to satisfy the
energy needs of the private consumers and, as well, to supply the centrifugal pump 6 (model CH 400) with
electrical energy in order to pump water in drip irrigation systems or for drainage of agricultural land adjacent
to river (with relocation of the centrifugal pump 6). In the case of electric energy production, with account of
mechanical losses both in the micro hydropower plant linkage and in the permanent magnet generator, the
efficiency of energy utilisation at generators terminals is,
1 2 r g
0,9 0,95 0,99 0,87 0,736, q q q q q
= = =
And in the case of water pumping (at the centrifugal pump shaft) the efficiency is:
1 2 r g me
0,9 0,95 0,99 0,87 0,91 0,67, q q q q q q
= = =
1. Blade with hydrodynamic NACA
0016 profile; 2 3-blade rotor; 3
planetary multiplier with multiplying
ratio i=112; 4 belt drive with
multiplying ratio i = 1,9; 5 - permanent
magnet generator (characteristics
see p. 5.4); 6. centrifugal pump CH
400 (characteristics pump flow rate
Q=(20-40)m
3
/h at pumping height
15...32m); 7 plastic pontoons; 8
guide; 10 space housing.
65
where: q
g
is the
generator efficiency;
q
me
is the
efficiency of the
hydraulic pump
electromotor.
Thus, the micro
hydropower plant
MHCF D4x1,5 ME
ensures the
transformation into
useful energy of
73,6% and 67% of the
energy potential of
flowing water, picked
up by the
hydrodynamic rotor, in
producing electrical
energy and,
respectively, in water
pumping.
5.2.3. Micro hydropower plant with hydrodynamic rotor for river water
kinetic energy conversion into electrical and mechanical energy at
small speeds (MHCF D4x1,5ME)
Micro hydropower plant MHCF D4x1,5 ME (figure 68) is designed to convert river water kinetic
energy into electrical and mechanical energy, by utilising low speed permanent magnet generator
5 (n = 375 min
-1
) and three-stage low speed centrifugal pump (PSS 40-10/50 (n = 500 min
-1
) 7
designed, in particular, for the micro hydropower plant and manufactured at Hidrotehnica S.A.,
Chiinu. Research results and functional characteristics of the low speed pump are presented in
[1,2].
Kinematics and functional principle of the micro hydropower plant are analogic to the micro
hydropower plant presented above (figure 67). Constructive peculiarities of this micro hydropower
plant refer, in particular, to the driving mechanism unit of the centrifugal pump (fig. 69) and to the
supply of the pump low speed electromotor 2 from the permanent magnet low speed generator 5
(figure 68). This design configuration can be used both to meet the needs of irrigation by pumping
water at relatively low heights (10 ... 15) m (e.g. over the river dam) and to perform the draining
works of the land adjacent to the river. When the micro hydropower plant is used for draining
works, the centrifugal pump driving mechanism (Fig. 70) is relocated from the spatial housing of
the micro hydropower plant to the floating platform in the flooded area of agricultural land adjacent
to the river.
Water current
Figure 66 Kinematics of micro hydropower plant MHCF D4x1,5 ME.
66
Figure 67 Micro hydropower plant with hydrodynamic rotor for river kinetic energy conversion into
electrical and mechanical energy used for water pumping (rotor diameter D = 4 m, submerged height
of blade h = 1,4 m, length of blade chord l = 1,3 m).
Three-stage centrifugal pump PSS 40-10/50 1 (fig. 69) is coupled to the electromotor 2 through the
toroidal coupling 3 and the housing 4 for the transmission of the reactive torque. Taking into
account the mechanical losses both in the linkage and in the electric generator, the kinetic energy
efficiency transmitted by the water flow to the hydrodynamic rotor at the permanent magnet
generator terminals makes up:
1 2 r g
0,9 0,95 0,99 0,87 0,736, q q q q q
= = =
1. Blade with hydrodynamic NACA
0016 profile; 2 3-blade rotor; 3
planetary multiplier with multiplying
ratio i=112; 4 belt drive with
multiplying ratio i = 1,9; 5 -
permanent magnet generator
(characteristics see p. 5.4); 6.
asynchronous electromotor; 7
centrifugal pump PSS-40-10/50
(characteristics Q and H see
p.5.4.3); 8 plastic pontoons; 9
guide; 10 spatial housing.
67
and at the inlet shaft of the hydraulic pump PSS4010/50 it is:
Figure 68 Unit of three-stage hydraulic pump
driving mechanism PSS 40-10/50.
Figure 69 Unit of low speed electric generator
driving mechanism (MCHF D4x1,5E).
1 2 r g m.e.
0,9 0,95 0,99 0,87 0,81 0,596. q q q q q q
= = =
The micro hydropower plant (MHCF D4x1,5 ME) with the given configuration ensures the
transformation of 73,6% of river water energy potential into useful energy and only 59,6% - at
water pumping. Relatively small efficiency at water pumping is due to a more reduced efficiency of
the low speed electromotor.
5.2.4. Micro hydropower plant with hydrodynamic rotor for river water
kinetic energy conversion into electrical energy (MHCF D4x1,5E)
The micro hydro power plant (figure 54) is designed to convert river water kinetic energy into
electrical energy only. The construction peculiarities are as follows: the hydraulic rotor comprises
five blades 1; the permanent magnet generator 3 (figure 53) is assembled coaxially with a
68
planetary multiplier 1 through toroidal coupling 2 and housing 4 for taking over the reactive moment
of torsion.
Figure 70 Micro hydro power plant with hydrodynamic rotor for river water kinetic energy
conversion into electrical energy (5-blade rotor diameter D = 4 m, submerged height of blade h = 1,4
m, length of blade chord l = 1,3 m).
1. Blade with hydrodynamic
NACA 0016 profile; 2 5 blade
rotor; 3 planetary multiplier with
multiplying ratio i=112; 4
permanent magnet generator
(characteristics see p. 5.3); 5
plastic pontoons; 6 guide; 7
spatial housing.
69
Dependence of summary torque T
1
applied to
the 5-blade rotor shaft depending on water flow
velocity (V = 1,3...1,8) m/s is presented in
figure 72. Kinematic and lifting capacity
calculations of all structural elements, including
functional parameters and technical
characteristics of micro hydropower plants have
been carried out for the torque value T
1
= 19893
Nm, corresponding to the water flow velocity V
= 1,3 m/s (maximum velocity specified for Prut,
Dniester and Raut rivers).
Efficiency of the kinetic energy transmitted by
the water flow to the hydraulic rotor can be
considered (at the permanent magnet generator
terminals):
1 r g
0,9 0,99 0,87 0,775 q q q q
= = = .
In conclusion, we state that micro hydropower plant MHCF D4x1, 5E ensures the transformation of
77.5% of the flowing water potential energy into useful electrical energy transmitted to the
hydrodynamic rotor.
According with obtained research results the industrial prototype of the microhydropower station for
the river kinetic energy conversion into electrical and mechanical energies (diameter of rotor d =
4m, submersed height of the blade h = 1,4m, length of blade l =1,3m) (MHCF D4x1,5 E) was
produced (figure 73).
Figure 72 Industrial prototype of the microhydropower station for the river kinetic energy
conversion into electrical energy (diameter of rotor d = 4m, submersed height of the blade h = 1,4m,
length of blade l =1,3m) (MHCF D4x1,5 E)
0 30 60 90 120 150 180 210 240 270 300 330 360
0
0.5
1
1.5
2
2.5
3
3.5
4
4.5
5
x 10
4
Unghiul de pozitionare, (Deg)
T
,
(
n
m
)
Momentul sumar T la diferite viteze, Profil:NACA 0016
v
0
=1.3 m/s
v
0
=1.6 m/s
v
0
=1.8 m/s
19893 N m
38137 N m
30133 N m
Valoare medie
Figure 71 Torque T1 at the shaft of 5-blade
hydrodynamic rotor with NACA 0016 profile
70
6. Summary and Conclusions
The paper presents a modern monitoring system of the hydraulic, mechanical and electrical
parameters related to a SHP and also a new concept with regard to the possibility of catching the
kinetic energy of a water stream.
The monitoring system of the SHP parameters has been implemented to a SHP on Arges river,
here being presented the equipment of the monitoring system, their arrangement within the power
house, the connection between the equipment, and the information processing and presentation
method. The paper presents also the fact that several small hydropower plants, having
implemented this monitoring system, can be interconnected and managed from a dispatcher
center. The most important issue is that the different failures, that can occur, can be analyzed and
interpreted accurately, especially the electrical failures and that cannot be other way interpreted.
As innovative technology for the catchment of the water kinetic energy it is presented a kinetic
turbine, completely new, of last generation, fully designed and manufactured by the technical
University of Moldavia and tested in situ in Prut river (located at the border between Romania and
Moldavia).
In the Annex there are presented new technologies developed for the catchment of the water
kinetic energy.
71
Annex
FREE FLOW TURBINE VERDANT POWER
-a three-blade horizontal-axis turbine designed to capture energy from both river and tidal currents
-the turbines are installed and operate fully under water, invisible from the shore
-spun slowly and steadily by underwater currents, the turbines rotor drives a gearbox, which in turn
drives a grid-connected generator
-the gearbox and generator are encased in a waterproof streamlined nacelle mounted on a pylon
-the pylon assembly has internal yaw bearings allowing it to pivot the turbine with the direction of
the rivers currents
-the pylon is bolted via an adjustable adapter to a pile fixed to the river bottom
-the turbine will operate below 1 m/s but for economic efficiency it recommends velocities
greater than of 2 m/s and water depths of at least 6.5 m
Figure 73 Free flow turbine, Verdant Popwer
72
FREE FLOW TURBINE UEK Corporation Underwater Electric Kite
-the system employs two axial-flow turbines in a side-by-side configuration. Each turbine consists
of five blades that dive a single internal generator housed within the nacelle.
-the system incorporates an augmenter ring that is integral with rear edge of the shroud. The
augmenter ring extends outwardly with respect to the axial alignment of the turbine shafts and
deflects the flow of water about the shroud. This creates a low pressure zone of the rear of the
shroud that pulls water through the turbine blades at a velocity greater than that of the normal
surrounding flow of water.
-the unit is positively buoyant and is secured to the river bed by a single anchorage system using a
cable bridle. When flown as a kite the angle of attack is altered by a patented ballasting system
that shifts a weight forward and back in the keel. Keeping a controlled operational depth, the units
are not affected by the surface effect of the large waves or navigation. Lateral positioning controls
permit the units to stay in the core of current.
-the turbine is designed to operate in river, tidal or ocean currents
-various models exist from 2 m to 5 m and operate in extremely low velocities of 0.2 m/s or less
Figure 74 Free flow turbine, UEK Corporation Underwater Electric Kite
FREE FLOW TURBINE SWAN TURBINE
-the unit is a three-bladed axial flow turbine
-a gearless low speed generator offers a high efficiency over a range of speeds with minimal
maintenance demands through the use of novel structural electromagnetic topologies
-a simple, robust and serviceable yawing mechanism is used for maximum flow capture
73
Figure 75 Free flow turbine, Swan turbine
FREE FLOW TURBINE GORLOV HELICAL
-a cross-axis turbine consisting of one or more helical blades that run along an imaginary
cylindrical surface of rotation like a screw thread
-the helical airfoil blades provide a reaction thrust perpendicular to the leading edges of the blades
that can pull them faster than the fluid flow itself
-the GHT allows a large mass of slow water to flow through, capturing its kinetic energy and
utilizing a very simple rotor
-it can be assembled vertically, horizontally or in any other cross-axis combination using common
shaft and generator for an array of multiple turbines
-generating capacity is proportional to the number of modules
-in its vertical orientation the generator and gearing can easily be positioned above water
-the standard unit is 1 m diameter by 2.5 m in length
-it starts producing power at approximately 0.6 m/s.
74
Figure 76 Free flow turbine, Gorlov helical turbine
FREE FLOW TURBINE MILLAU VLH (VERY LOW HEAD TURBINE)
-installed capacity 410 kW
-commissioned on March 19, 2007
-a DN 4500 (runner diameter 4.5 m).
-the 410 kW max nominal power at grid was reached at the nominal speed of 37 rpm
-very smooth, vibration-free and silent operation (one must touch the machine to find out whether it
is operating)
Figure 77 Free flow turbine, Millau VLH
75
7. References
References (RO)
1. Cojocar, Mihai. Hidroconstructia 2005. 2005.
2. Pavel, Dorin. Masini Hidraulice. Bucuresti : Editura Energetica de stat, 1965.
3. Arbiter Systems. Model 1133A Power Sentinel GPS-Synchronized Power Quality/Revenue Standard -
Operation Manual.
4. www.elpros.si. UniFusion. ELPROS. [Online] ELPROS. http://www.elpros.si/eng/UF_base.htm.
5. National Instruments. NI cRIO-9022 Operating Instructions and Specifications.
6. Carlo Gavazzi. PQT H Smart Modular Network Power Quality Transducer - Instruction Manual.
7. Boyle, Godfrey. Renewable Energy - Power for a Sustainable Future. s.l. : Oxford University Press.
8. Silviu, Folea. LabVIEW - Practical Applications and Solutions. s.l. : InTech, Croatia, 2011.
9. Alexandru Fransua, Razvan Magureanu. Electrical machines and Drive Systems. s.l. : Technical Press,
Oxford, UK, 1984.
10. Aquatic Renewable Energy Technologies AQUARET, Leonardo da Vinci, Project No: IRL/06/B/F/PP-
153111, 2006, www.aquaret.com.
11. Websites: www.vlh-turbine.com, www.swanturbines.co.uk, www.uekus.com,
www.verdantpower.com, www.gcktechnology.com, en.wikipedia.org
References (MOLD)
Monographs:
1. I. Bostan, V. Dulgheru, I. Sobor, V. Bostan, A. Sochirean. Renewable energy conversion systems: / - Ch. :
Tehnica-Info, 2007. 592pp.
2. Bostan, V. Dulgheru, V. Bostan, R. Ciuperc. Anthology of inventions: renewable energy conversion
systems: / - Ch.: Bons Offices SRL, 2009. 458pp.
3. Jula A., Mogan Gh., Bostan I., Dulgheru V. et al. ECOMECA ECO- mechanical engineering. Braov,
Publ. House of Transilvania University, Braov, p.324.
4. BOSTAN, I.; GHEORGHE, A.; DULGHERU, V.; BOSTAN, V.; SOCHIREANU, A.; DICUSAR, I.
Conversion of Renewable Kinetic Energy of Water: synthesis, Theoretical Modeling, and Experimental
Evaluation. Energy Security: International and Local Issues, Theoretical Perspectives, and Critical Energy
Infrastructures (NATO Science for Peace and Security Series - C: Environmental Security). 2011. Published
by Springer, p. 125-177. ISBN 978-94-007-0718-4
Article:
5. Bostan I., Dulgheru V., Sobor I., Bostan V., Sochirean A. Valorisation of renewable energy // ENERG VI:
Energy, Environment, Economy, Resources, Globalization. Publ. AGIR, Bucureti. P. 152-205. ISBN 978-
973-720-263-5.
6. Bostan I., Bostan V., Dulgheru V. Microhidro stations / Seminar addressed to hydropower stakeholders in
Moldova - SEE HYDROPOWER - clear water, clean energy (SEE - South East Europe) 24.03.2010.
7. Dulgheru V. Utilisation of renewable energy sources - wind, solar and hydro in the Republic of Moldova.
Meridian Ingineresc, nr. 3, 2009. P. 63-69.
8. Bostan I., Dulgheru V. Renewable energy conversion systems one of basic element for sustenable
development of society. // Pro-Active Partnership in Creativity for the Next Generation. Proceedings: 32
st
ARA Congress. - Sibiu, 2009. - P. 78-82.
76
Patents:
9. PATENT 2981 (MD), CIB B 63 B 35/44; E 02 B 17/00. Hydraulic station / I. BOSTAN, V. Dulgheru, V.
Bostan. R. Ciuperc. Publ. BOPI 2006. - Nr. 2.
10. PATENT 2991 (MD), CIB F03 B 7/00. Hydroelectric station / I. BOSTAN, V. Dulgheru, V. Bostan, O.
Ciobanu, A. Sochireanu. Publ. BOPI 2006. - Nr. 2.
11. PATENT 2992 (MD), CIB F 03 B 7/00. Hydraulic station / I. BOSTAN, V. Dulgheru, A. Sochireanu, V.
Bostan, O. Ciobanu, R. Ciobanu. Publ. BOPI 2006. - Nr. 2.
12. PATENT 2993 (MD), CIB F 03 B 7/00; F 03 B 13/00. Hydraulic turbine / I. BOSTAN, V. Dulgheru, V.
Bostan, A. Sochireanu, N. Trifan. Publ. BOPI 2006. - Nr. 2.
13. PATENT 3104 (MD), CIB F 03 B 7/00: F 16 H 1/00. Hydraulic station / I. BOSTAN, V. Dulgheru, V.
Bostan, A. Sochireanu, O. Ciobanu; R. Ciobanu, I. Dicusar. Publ. BOPI2006. -Nr. 7.
14. PATENT 3845 (MD), CIB F 03 B 13/00; F 03 B 7/00; F 03 B 13/10; ; F 03 B 13/22; ; F 03 B 17/06.
Hydraulic station / I. BOSTAN, V. Dulgeru, V. Bostan, A. Sochireanu, O. Ciobanu, R. Ciobanu. Publ. BOPI
2009. - Nr. 2.
15. PATENT 3846 (MD), CIB F 03 B 13/00; F 03 B 7/00; F 03 B 13/18; ; F 03 B 13/22; ; F 03 B 17/06.
Hydraulic station with horizontal axle / I. BOSTAN, A. Gheorghe, V. Dulgheru, V. Bostan, A. Sochireanu, O.
Ciobanu, R. Ciobanu. Publ. BOPI 2009. - Nr. 2.
Presentation on the International Salon of Research and Innovations
16. Bostan Ion, Dulgheru Valeriu, Sobor Ion, Bostan Viorel, Sochireanu Anatol, Crudu Radu, Guu Marin,
Ciobanu Oleg, Ciobanu Radu, Trifan Nicolae. Industrial prototype of mini hidropower station for flow
water kinetic energy conversion. Salon des Inventions, Geneva,- Palexpo, 6 au 10 avril 2010 (Silver
medal).
17. Bostan Ion, Dulgheru Valeriu, Sobor Ion, Bostan Viorel, Sochireanu Anatol, Crudu Radu, Guu Marin,
Ciobanu Oleg, Ciobanu Radu,Trifan Nicolae. Industrial prototype of mini hidropower station for flow
water kinetic energy conversion. XIIIth Moskow International Salon of Research and Innovations
ARHIMED-2010. 30.03..02.04.2010 (Gold medal).
18. Bostan I. (MD), Dulgheru V. (MD), Bostan V. (MD), Ciobanu O. (MD), Ciobanu R. (MD), Sochireanu A.
(MD), Dicusar I. (MD), Trifan N. (MD). Industrial prototype of mini hidropower station for flow water
kinetic energy conversion into electrica land mechanical energy. (DIPLOM i Medalia de aur. Premiul
Special al Asociaiei Inventatorilor din Zagreb). EUROINVENT"-European Exhibition of Creativity and
Innovation-Iai, Romnia. 07..09.05.2010 (Gold medal).
19. Bostan Ion, Dulgheru Valeriu, Bostan Viorel, Sochireanu Anatol, Ciobanu Oleg, Ciobanu Radu, Dicusar
Ion, Trifan Nicolae. Industrial prototype of mini hidropower station for flow water kinetic energy
conversion into electrica land mechanical energy. International Salon of Research and Innovations,
INVENTICA 2010, XIVth edition, 9 - 11 June 2010 (Gold medal).
20. Bostan Ion, Dulgheru Valeriu, Bostan Viorel, Sochireanu Anatol, Ciobanu Oleg, Ciobanu Radu, Dicusar
Ion, Trifan Nicolae. Industrial prototype of mini-hydropower station for flow water kinetic energy conversion.
EURECA 2009, Bruxel (Gold medal).
21. DIPLOMA. Awarded to: I. Bostan, V. Dulgheru, V. Bostan, A. Sochirean, O. Ciobanu, R. Ciobanu, N.
Trifan for the Floatable micro-hydropower station. PRIZE ENVIRONMENT PROTECTION.
EUROINVENT2009 - Iai, 9/05/2009.
22. I. Bostan, V. Dulgheru, V. Bostan, A. Sochirean, O. Ciobanu, R. Ciobanu, Nicolae Trifan for the Floatable
micro-hydropower station with adjustable blades. EUROINVENT2009.- Iai, 9/05/2009. (Gold medal).
23. I. Bostan, I. Via, V. Dulgheru, V. Bostan, A. Sochireanu, O. Ciobanu, N. Trifan. Micro-hydropower station
for the rivers water kinetic energy conversion // - Cluj-Napoca, 2009 The International Salon of Research
and Innovations PROINVENT 2009. The DIPLOMA of EXCELLENCE and PROINVENT medal.
24. I. Bostan, A.Greorghe, V. Dulgheru, V. Bostan, A. Sochireanu, V. Cartofeanu, O. Ciobanu, R. Ciobanu, I.
Dicusar, N. Trifan. Flotable micro-hydropower station with self-oriented hydrodynamic blades // - Cluj-
Napoca, 2011 The International Salon of Research and Innovations PROINVENT 2011. The DIPLOMA of
EXCELLENCE and (Gold medal).
77
Authors Contact
Razvan Magureanu (POLI-B)
[email protected]
Telephone: +40 722228514
Fax: +40 214029342
Sergiu Ambrosi (POLI-B)
[email protected]
Telephone: +40 721761481
Fax: +40 214029342
Bogdan Popa (POLI-B)
e-mail: [email protected]
Telephone: +40 214029189
Fax: +40 214029189
Bostan Ion, Dulgheru Valeriu,
Bostan Viorel, Sochirean Anatol (MOLD)
e-mail:[email protected]
www.seehydropower.eu
Project Contact
Ing. Maximo Peviani
[email protected]
Telephone: +39 035 55771 (switchboard)