AP05_3.vp Notation Symbol Unit Description a Interference factor c m Blade chord length Cd Airfoil drag coefficient Cl Airfoil lift coefficient Cmc/4 Airfoil quarter chord pitching moment coefficient CD Blade drag coefficient CL Blade lift coefficient Cp Turbine performance coefficient Cq Turbine torque coefficient D N Drag dF N Elementary force acting on the elemen- tary actuator disk Ip kg m 2 Blade moment of inertia IT kg m 2 Turbine moment of inertia L N Lift M Nm Instantaneous turbine torque MC Nm Instantaneous load torque Mc/4 Nm Quarter chord pitching moment Mm Nm Average turbine torque N N Blade radial force Nb Number of blades P W Instantaneous turbine mechanical power Pm W Average turbine mechanical power R m Turbine radius Re Blade Reynolds number S m2 Turbine frontal area T Nm Blade tangential force V m /s Local velocity VR m /s Tip speed V � m /s Asymptotic velocity xc/4 %blade chord Blade aerodynamic centre position xhinge %blade chord Floating hinge position � rad Blade angle of attack �tan rad Angle between the local velocity and the local tangent at the blade �zv rad Blade pitch angle ��� rad /s2 Blade pitch angle acceleration � Tip speed ratio � �R/V� � kg /m3 Fluid density � Solidity � Nb c /R � rad Blade azimuth angle ��� rad /s 2 Turbine acceleration � rad /s Turbine angular velocity Subscript d Conditions at the downwind actuator disk u Conditions at the upwind actuator disk h Hinge 1 Introduction Marine current energy is a type of renewable energy resource that has been less exploited than wind energy. Only recent years, have some countries devoted funds to research aimed at developing tidal current power stations. Tidal cur- rent turbines, as in the wind community, can be divided into vertical-axis and horizontal axis types. Although horizontal axis turbines have been more widely used than vertical axis types for wind energy exploitation, vertical axis turbines could present significant advantages for tidal current ex- ploitation, because they are simple to build and reliable in working conditions. Therefore, at beginning of the studies © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 77 Czech Technical University in Prague Acta Polytechnica Vol. 45 No. 3/2005 Dynamic Behaviour of the Patented Kobold Tidal Current Turbine: Numerical and Experimental Aspects D. P. Coiro, A. De Marco, F. Nicolosi, S. Melone, F. Montella This paper provides a summary of the work done at DPA on numerical and experimental investigations of a novel patented vertical axis and variable pitching blades hydro turbine designed to harness energy from marine tidal currents. Ponte di Archimede S.p.A. Company, located in Messina, Italy, owns the patented KOBOLD turbine that is moored in the Messina Strait, between the mainland and Sicily. The turbine has a rotor with a diameter of 6 meters, three vertical blades of 5 meters span with a 0.4 m chord ad hoc designed curved airfoil, producing high lift with no cavitation. The rated power is 160 kW with 3.5 m/s current speed, which means 25% global system efficiency. The VAWT and VAWT_DYN computer codes, based on Double Multiple Steamtube, have been developed to predict the steady and dynamic performances of a cycloturbine with fixed or self-acting variable pitch straight-blades. A theoretical analysis and a numerical prediction of the turbine performances as well as experimental test results on both a model and the real scale turbine will be presented and discussed. Keywords: vertical-axis-hydro-turbine, variable pitch, Double-Multiple-Streamtube, tidal currents, tidal energy. vertical axis wind turbines were taken as models for hy- dro-turbines. The blades of Darrieus-type vertical axis wind turbines are fixed, and they perform well when the blade so- lidity is low and the working speed is high. For this reason, the first hydro-turbines were impossible to start. A variable-pitch blade system can be a solution to this problem. Some proto- types with different variants of this system have therefore been developed around the world: the Kobold turbine in the Strait of Messina, Italy; the cycloidal turbine in Guanshan, China; the moment-control turbine at Edinburgh University, UK; and the mass-stabilised system turbine by Kirke and Lazauskas in Inman Valley, South Australia. The Kobold tur- bine has been under development since 1997: the rotor has a self-acting variable pitch and the Kobold blades have an ad hoc designed airfoil, called HLIFT, to be cavitation free and to have high lift performance. The methods for calculat- ing the hydrodynamic performances of vertical axis turbines also come from wind turbines: in the 1970s Templin devel- oped the Single-Disk Single-Tube model, and then Strickland put forward the Single-Disk Multi-Tube model. In the 1980s Paraschivoiu introduced the Double-Disk Multi-Tube model. The VAWT and VAWT_DYN computer codes, based on this theory, have been developed to predict the steady and dy- namic performances of a cycloturbine with fixed or self-acting variable pitch straight-blades. The numerical results have been compared with two sets of experimental data: one set is obtained from wind tunnel tests on a scaled model, and the other set is based on field data from the Kobold prototype. 2 Double multiple streamtube In order to analyze the flow field around a vertical axis turbine, a DMS model was used. The DMS model is an evolu- tion of the previous “momentum models”: the single streamtube model, the multiple streamtube model and the double streamtube model [1]. The DMS model [2] assumes that the flow through the rotor can be modelled by examining the flow through several streamtubes, and the flow disturbance, produced by the rotor is determined by equating the aerody- namic forces on the turbine rotor to the time rate of change in momentum through the rotor as depicted in Fig. 1. In the DMS model, the flow velocities vary in both the upwind and downwind regions of the streamtube, as well as varying from streamtube to streamtube. So DMS is able to analyse the inter- ference between the downwind blade and the upwind blade’s wake in order to evaluate more accurately the local value of the velocity and the instantaneous blade load. As shown in Fig. 1., the rotor is modelled as a series of elementary streamtubes, and each streamtube is modelled with two actua- tor disks in series. Across the actuator disk the pressure drops and this drop is equivalent to the streamwise force dF on the actuator disk divided by the actuator disk area dA. The elementary force dFu and dFd, respectively on the up- wind and downwind disk, given by the momentum principle, are d du u uF V A V V� ��� ( )2 (1) d dd d dF V A V V� �� ( )2 3 (2) where Vd, which is the velocity on the downwind actuator disk, is influenced by the velocity Vu on the upwind actuator disk. The elementary forces dF on the actuator disks may be calculated using Blade Element Theory. If the upwind and downwind interference factors are de- fined as a V V V a V V Vu u d d� � � �� � � � (3) the mathematical problem can be reduced to the calculation of au and ad. Because of the non-linearity of the equations, the problem must be resolved iteratively. If the rotor blades have a fixed pitch angle or an assigned pitch variation (i.e. sinusoi- dal like in Pinson, cycloidal, etc.), the mathematical model is reduced, for each elementary streamtube, to an equation for the momentum balance for the upwind actuator disk and an equation for the momentum balance for the downwind actua- tor disk. � ��a a sen V V C senu u Ru lu u( ) ( ) tan1 1 8 2 2 � � � � � � � � �� � � � � � � ��C a a a sen V V du u d d u Rd cos ( ) ( )( ) tan2 1 2 1 8 � � � � � � � � � � � � � � � � � � � � 2 C sen C ld d dd d ( ) cos( ) tan tan � � � � (4) If the rotor blades have a self-acting variable pitch angle [3], [4], [5], another equation is also necessary for each actua- tor disk: the hinge moment equilibrium. In this case, in fact, the blade is partially free to pitch under the action of the aero- dynamic and inertia forces so as to reduce the angle of attack 78 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 45 No. 3/2005 Czech Technical University in Prague Downwind disk d�d � i j u Upwind disk Elementary streamtube Circular path of the blade V V Downwind disk dAddAu V V Upwind disk V dFu dFd V V V V V V V d� Fig. 1: Double multiple streamtube model and hence the tendency of the blade to stall. The allowed an- gular swinging of the blade is limited by the presence of two blocks. In this way the mathematical model is represented by two systems of equations, each constituted of two equations: momentum balance and hinge moment equilibrium. For one blade and for the upwind actuator disk a a sen V V Cu u u Ru lu u zvu u( ) ( , ,Retan1 1 8 2 � � � � � �� � � � �� � ) ( ) ( , ,Re ) cos( ) tan tan tan � � � � � sen C C � � � � � � u u du u zvu u u u mc u zvu u nu u zvu u c cer 4 4 ( , ,Re ) ( , ,Re ) ( ) tan tan� � � �� � � C x x � � � � � � cos ( , ,Re ) ( ) tan � � � � zvu tu u zvu u c cer zvu C x x sen4 0 (5) For the downwind actuator disk ( )( ) ( ,tan 1 2 1 8 2 � � � � � � � a a a sen V V C d d u d Rd ld d z � � � � �� vd d d d dd d zvd d d ,Re ) ( ) ( , ,Re ) cos( tan tan� � � � � sen C� � � � � � �tan tan tan ) ( , ,Re ) ( , ,Re ) ( d mc d zvd d nd d zvd d c C C x 4 4 � � � �� � � � � � � � x C x x sen cer zvd td d zvd d c cer ) cos ( , ,Re ) ( ) tan � � � �4 zvd � 0 (6) The instantaneous torque and power produced by the blade are given by the moment equilibrium around the tur- bine axis � � M N x x T R x x sen P M � � � � � � � ( ) cos ( ) c hinge zv c hinge zv 4 4 � � � (7) To obtain the mean torque and mechanical power pro- duced by Nb blades in a revolution it is necessary to average the instantaneous values. �M N x x T R x x sen m b c hinge zv c hinge � � � � � �2 4 0 2 4 � � � � ( ) cos ( )� ��zv d ,� (8) �P N x x T R x x sen m b c hinge zv c hinge � � � � � � � 2 4 0 2 4 � � � ( ) cos ( )� ��� �zv d . (9) To simulate dynamic performances, we have to resolve only the equation of the moment equilibrium around turbine axis (10) for fixed blade or Nb� 1 equations for Nb floating blades around their hinge axis (11). I M MT i i N �� � � � � C b 1 (10) I M I M I P zv1 h1 P zv h2 P zv (�� ��) (�� ��) (�� ��) � � � � � � � � 2 3 � � � M I M h3 P zvn hn � (�� ��)� (11) 3 Airfoils and aerodynamic characteristics The rotor performances were tested using five airfoils: three symmetrical airfoils NACA 0012, 0015, 0018 and two cambered airfoils NACA 4415 and HLIFT18. The last named was designed at DPA; it is a high lift [6], [7] and no cavitating airfoil. It has in fact been designed to work in the water on the KOBOLD turbine. For the NACA airfoils, data was taken from the literature [2], [8] while for the HLIFT18 airfoil, TBVOR [9], [10], [11] codes were used to generate the values of the aerodynamic coefficients. The airfoil 2D data is corrected, in VAWT and VAWT_DYN codes [12], to take into account the three-dimensional effects due to the blade finite aspect ratio. To take into account the three-dimensional effects, Prandtl’s lifting line theory, extended to treat high lift flow, has been used, evaluating in this way the 3D lift curve beginning from the 2D data. This theory is valid only in the linear zone of the lift curve but with care it can also be extended to non linear conditions. The total blade drag coefficient is the sum of the airfoil drag coefficient (Cd), due to skin friction, and the induced drag coefficient. To take into account the inter- ference between the blade and the support arms, a further drag coefficient increment, �CD, has been introduced. More- over, 2D post-stall modelling, based on the Viterna-Corrigan correlation method, has been introduced to extend the 2D aerodynamic coefficients to this angle range. © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 79 Czech Technical University in Prague Acta Polytechnica Vol. 45 No. 3/2005 V arm � �R V V n n t c L D Mc/4 �zv Fig. 2: Hinge moment equilibrium j i Mh1 Mh2 Mh3 M � MC TI �zv : Fig. 3: Torques on a Kobold turbine 4 Experimental models The reported experimental data is divided into two parts: � Experimental data measured in the DPA wind-tunnel at the University of Naples on a small straight-bladed cyclo- turbine, which was designed, developed and assembled at DPA [13]; � Experimental data measured in water on the “Kobold” prototype (real scale) [12], [14]. Both the DPA straight-bladed cycloturbine (Model A) and the “Kobold” prototype (Model B) will be described, as follows. Both turbines have variable pitch blades with a self- -acting system, made up of two balancing masses for each blade. In this way, the centre of gravity of the blade can be moved into its optimal position in order to optimize the global performance of the rotor, and, using two stops, the pitch blade range can be limited, as shown in Fig. 5. Model A in the DPA wind-tunnel is shown in fig. 6. Using different stop positions, it was possible to test different pitch angle ranges, and while using different numbers of blades it was possible to take into account different solidity [Nc/R] val- ues. Model A has the following geometric parameters: number of blades tested 2, 3, 4, 6 blade chord 0.15 m blade airfoil NACA 0018 blade span 0.8 m Aspect Ratio 5.33 radius 1.05 m number of radial arm 4, 6, 8, 12 arm chord 0.05 m solidity 1 0.286 solidity 2 0.428 solidity 3 0.571 solidity 4 0.857 80 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 45 No. 3/2005 Czech Technical University in Prague 0 0.2 0.4 0.6 0.8 1 x/c -0.2 -0.1 0 0.1 0.2 y/c NACA 0012 NACA 0015 NACA 0018 0 0.2 0.4 0.6 0.8 1 x/c -0.2 -0.1 0 0.1 0.2 y/c NACA 4415 HLIFT18 -30 -20 -10 0 10 20 30 A lfa [de g ] -1.5 -1 -0.5 0 0.5 1 1.5 2 C L NA CA 00 12 NA CA 00 15 NA CA 00 18 NA CA 44 15 HLI FT 18 -30 -20 -10 0 10 20 30 A lfa [d e g ] 0 0.1 0.2 0.3 0.4 0.5 0.6 C D N A C A 0 0 12 N A C A 0 0 15 N A C A 0 0 18 N A C A 4 4 15 H L IFT 1 8 Fig. 4: Airfoils tested and aerodynamic characteristics (Re � 10e6) � i Stop V� Hinge Stop Pitch Range j � Fig. 5: Balancing mass and blade stops Fig. 6: DPA straight-bladed cycloturbine (Model A) The “Kobold” prototype (Model B) lies out in the Strait of Messina, close to the Sicilian shore, facing a village called Ganzirri, close to the lake of the same name, as shown in Fig. 7. In this site the peak current speed is 2 m/s (4 knots), the sea depth is 20 meters and the plant has been moored 150 meters offshore. The current changes direction every 6 hrs and 12 minutes, and the amplitude period is equal to 14 days. A high lift airfoil, called H-LIFT18, is used for the blade sections and has been specially designed at DPA to be cavitation free and to optimise the turbine performance. Two arms sustain each blade and the arms have been stream- lined using another ad hoc designed symmetrical airfoil. The turbine has a very high starting torque, being able in this way to start spontaneously, also with electrical load connected, without the need for any starting devices. The ENERMAR plant is composed of the turbine rotor hanging under a floating buoy that contains the remaining mechanical and electrical parts to deliver energy to the grid, as shown in Fig. 7. The rotor has a diameter of 6 meters with 6 radial arms holding three blades with a five-meter span and with a chord of 0.4 m employing the H-LIFT18 airfoil leading to an aspect ratio of 12.5 and to solidity � � 0.4. The acquisition data sys- tem is made up of a torque-meter, a tidal current speed-meter and an RPM counter all connected to a PLC that converts analog signals to digital data transferring them to a PC. Data handling software has been developed to monitor the ac- quired data in real time. Fig. 8 shows some components of the acquisition data system. The PLC also acts as an electrical load controller to keep the turbine working always at its maxi- mum efficiency independently from the current speed. 5 Experimental tests Fig. 9 and Fig. 10 compare the VAWT code numerical re- sults with the experimental data measured on Model A, for © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 81 Czech Technical University in Prague Acta Polytechnica Vol. 45 No. 3/2005 Fig. 7: The Strait of Messina, the position of the plant and a picture of the plant Fig. 8: The torque-meter, the PLC and the PC 0 100 200 300 400 Rotor Angular Velocity (rpm) 0 50 100 150 200 250 M e c h a n ic a l R o to r P o w e r (W ) pitch range = 0° - 10° blades = 2 Solidity = 0.256 experimental data V AWT r esults 0 100 200 300 400 Rotor Angular Velocity (rpm) 0 50 100 150 200 250 M e c h a n ic a l R o to r P o w e r (W ) pitch range = 0° - 10° blades = 3 Solidity = 0.428 experim ental data VAWT results Fig. 9: Experimental data and VAWT results (Model A) 82 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 45 No. 3/2005 Czech Technical University in Prague 0 100 200 300 400 Rotor Angular Velocity (rpm) 0 50 100 150 200 250 M e c h a n ic a l R o to r P o w e r (W ) pitch range = 0° - 10 ° blad es = 4 Solidity = 0.571 experimental data V AWT results 0 100 200 300 Rotor Angular Velocity (rpm) 0 50 100 150 200 250 M e c h a n ic a l R o to r P o w e r (W ) pitch range = 0° - 10° blades = 6 Solidity = 0.857 experim ental data VAWT results Fig. 10: Experimental data and VAWT results (Model A) 1.6 1.8 2 2.2 2.4 2.6 2.8 TSR 0.12 0.16 0.2 0.24 0.28 0.32 0.36 N e t ro to r p o w e r c o e ff ic ie n t experimental data VAWT results Fig. 11: Experimental data and VAWT results (Model B) 0 10 20 30 40 50 2 4 6 8 10 12 numerical experimental 0 10 20 30 40 50 -40 0 40 80 120 160 T o rq u e (N m ) numerical experimental Fig. 12: Experimental data and VAWT_DYN results (Model B) different blade numbers and different pitch range. The pitch angle is measured between the local tangent to the blade circular path and the blade trailing edge, starting from the tangent in clockwise direction. The power measured is the net rotor power, and the tests are carried out with 9 m/s air speed in the DPA wind tunnel for Model A. A comparison of the numerical results with the experimental data shows good agreement, especially for 3 and 4 blades; in the case of 6 blades, the solidity is very high and there is a strong wake effect on each blade. Fig. 11 compares the experimental data measured during field tests in water for the Kobold turbine prototype – model B – with the VAWT code numerical predic- tion. The net rotor power coefficient is measured, and the tests are carried out with a presumed current speed of 1.4 m/s, but there is uncertainty around 25 % on the real “undis- turbed” current speed value, which is strongly influenced by the location of the current speed meter: this is at present be- ing investigated. Fig. 12 and Fig. 13 show the experimental data and VAWT_DYN code numerical results referred to the starting condition for the Kobold turbine prototype (Model B). The rotor angular velocity variation in time predicted by the code seems to be very accurate, while the rotor torque and power amplitudes are in good agreement only in the first part of the time range: this is probably due to the uncertain value of the numerical predicted losses. The frequency of the torque and power is, however, very well predicted. 6 Conclusions The Double Multiple Steamtube model seems to be good in predicting the vertical axis turbine performances, with fixed or floating blades, expecially for solidity values less than 0.5. This model has been implemented in VAWT and VAWT_DYN codes, which are then capable of predicting both static and dynamic performances of the turbine with very low computing time. In these codes blade 3D effects have been included as well as arm losses, while it is difficult to predict the losses due to other effects. In the field tests, the accuracy of the tidal speed current measurement is a problem, because it is difficult to set the current speed-meter in the “real” un- disturbed flow, so the measurements have an uncertainty level around 20 %. On the Kobold turbine, the HLIFT18 non-symmetrical and non-cavitating airfoil gives better per- formance than a symmetrical airfoil, and the self-acting vari- able pitch system with balancing mass has proven to be simple to build and more reliable than other more complex systems. References [1] Strickland, J. H.: “A Review of Aerodynamic Analysis Methods for Vertical-Axis Wind Turbine.” In: Fifth ASME Wind Energy Symposium, SED – Vol. 2 edited by A. H. P. Swift. [2] Paraschivoiu, I.: Wind Turbine Design with Emphasis on Darreius Concept. Ecole Polytechnique de Montreal, Poly- technic International Press, 2002. [3] Kentfield, J. A. C.: “Cycloturbines with Freely Hinged Blades or Freely Hinged Leading Edge Slats.” In: Alter- native Energy Sources V. Part C: Indirect Solar/Geo- thermal. (Editor: T. N. Veziroglu). Amsterdam: Elsevier Science Publishers B. V., 1983, p. 71–86. [4] Lazauskas, L.: “Three Pitch Control Systems for Vertical Axis Wind Turbines Compared.” Wind Engineering, Vol. 16 (1992), No. 5, p. 269–281. [5] Kirke, B. K., Lazauskas, L.: “Experimental Verification of a Mathematical Model for Predicting the Perfor- mance of a Self-Acting Variable Pitch Vertical Axis Wind Turbine.” Wind Engineering, Vol. 17 (1993), No. 2, p. 58–66. [6] Healy, J. V.: “The Influence of Blade Camber on the Output of Vertical Axis Wind Turbine.” Wind Engineer- ing, Vol. 2 (1978), No. 3, p. 146–155. [7] Healy, J. V.: “The Influence of Blade Thickness on the Output of Vertical Axis Wind Turbine.” Wind Engineer- ing, Vol. 2 (1978), No. 1, p. 1–9. [8] Reuss, R. L. et al.: Effects of Surface Roughness and Vor- tex Generators on the NACA 4415 Airfoil. Report: NREL/TP-442-6472. Golden, Colorado: National Re- newable Energy Laboratory, December 1995. [9] Coiro, D. P., de Nicola, C.: Prediction of Aerodynamic Per- formance of Airfoils in Low Reynolds Number Flows. In: Low © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 83 Czech Technical University in Prague Acta Polytechnica Vol. 45 No. 3/2005 0 10 20 30 40 50 ti m e ( s ) 0 4 8 12 M e c h a n ic a l p o w e r (k W ) n u m e r i c a l e x p e r im e n t a l Fig. 13: Experimental data and VAWT_DYN results (Model B) Reynolds Number Aerodynamics Conference, Notre Dame, Indiana, U.S.A., June 1989. [10] Coiro, D. P., de Nicola, C.: “Low Reynolds Number Flows: The Role of The Transition.” In: X Congresso Nazionale AIDAA, Pisa, Ottobre 1989. [11] Coiro, D. P.: “Convergence Acceleration Procedure for a Viscous/Inviscid Coupling Approach for Airfoil Perfor- mances Prediction.” In: XI AIDAA National Congress, Forli’, Italy, October 1991. [12] Montella, F., Melone S.: “Analisi Sperimentale e Nume- rica del Comportamento Statico e Dinamico di una Cicloturbina ad Asse Verticale.” Aerospace Engineering Bachelor Thesis. Napoli, Italy: Dipartimento di Proget- tazione Aeronautica, December 2003. [13] Coiro, D. P., Nicolosi, F.: “Numerical and Experimental Tests for the Kobold Turbine.” In: SINERGY Sympo- sium, Hangzhou, Republic of China, November 1998. [14] Coiro, D. P. et al.: “Exploitation of Marine Tidal Currents: Design, Installation and Experimental Results for the Patented Kobold Vertical Axis Hydro Turbine.” In: Poster-Session OWEMES 2003. Napoli, Italy, 10–12 April 2003. Prof. D. P. Coiro phone: +39 081 7683322 fax: +39 081 624609 e-mail: coiro@unina.it Dr. A. De Marco e-mail: agodemar@unina.it Dr. F. Nicolosi e-mail: fabrnico@unina.it Dr. S. Melone e-mail: stmelone@unina.it Dr. F. Montella e-mail: francmont@genie.it Dipartimento di Progetazione Aeronautica (DPA) University of Naples “Federico II” Via Claudio 21, 80125 Naples – Italy 84 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 45 No. 3/2005 Czech Technical University in Prague