THE DEVELOPMENT OF A SYSTEM TO OBTAIN ENERGY FROM TIDAL CURRENTS

V. Ramirez (1) , D. Uribe (2), E. Chandomi(3).

(1) Renewable Energy Unit, Centro de Investigación Científica de Yucatán, Calle 43 No. 130, Col. Chuburná de Hidalgo, Mérida Yucatán, Mexico, Post Code: 97200, Phone: +5219997382748 e-mail: [email protected]

Abstract – In this research is described the design and development of a system capable of harnessing ocean current for the generation of energy. For this purpose, a vertical axis helical turbine will be employed. One of its advantages is that the design is unidirectional, it can be positioned vertically or horizontally, it has auto start; this means that compared among other similar models it does not need any external force to be activated and is 35% more efficient. The objective of the research is designing a turbine capable to maximize its torque, being aided with the design software SolidWorks, later with a physical prototype get the experimental data and proceed with the creation of the final device for the generation of electric energy.

Keywords: Turbine, Helical, Solidity, Gorlov, Blade.

1. Introduction – Global warming caused by greenhouse gases has made us reconsider new alternatives of energy production and consumption. Currently the increase of CO2 on the atmosphere is mainly caused by the combustion of fossil fuels which are employed for transportation and generation of electric energy [1]. According to IEA, in 2016 in Mexico the main energy generation sources were gas, oil and carbon, with a production goes from 33926 GWh for the carbon to 192259 GWh for the gas. As an alternative for the usage of fossil fuels and to help reducing the greenhouse gases it is important the exploitation of alternative energy sources which are those that generate renewable energies. One of this energy sources is the one that we can find on the ocean currents, which research and development of new devices capable of harnessing its potential approximated to 800 THW/year has increased the last years [2]. Some of the advantages worth mentioning of this energy source is the capacity to predict its availability, its lack of visual or noise pollution and its requirement for less space compared to , this is due to the higher density of the water [3]. One way to harness this kind of energy is trough hydrokinetic turbines which extract mechanical energy from the movement of water, these work under an essentially identical principle applied for wind power [4]. Gorlov helical cross-flow turbine is an example of a hydrokinetic turbine, which is a vertical type turbine. Some of its advantages are its auto start function, its unidirectional design and its superior 35% efficiency compared to other similar models [5]. For this research we focused on the design of a Gorlov-type turbine capable of generating 1.7 Kw of energy.

GEET-19, Paris, 24-26 July 2019 Pag. 83 Nomenclature

P: Power 퐶퐷: Drag coefficient η: Efficiency 퐶퐿: coefficient 휌: Density D: Drag force (N) 2 퐴푇: Cross-sectional area (푚 ) L: Lift force (N) V: Current Velocity (푚 ∕ 푠) F: Force (N) L: Length (m) T: Torque (Nm) D: Diameter (m) R: Radio (m) AR: Aspect Ratio ω: Angular Velocity (rad/s) σ: Solidity TSR: Tip Speed Ratio n: Number of blades 휑0: Torsion Angle (°) c: Chord (m) δ: Lift Angle (°) 퐶푇: Torque coefficient α: Attack Angle (°)

2. Experimental

2.1. Calculation for the Turbine design - For the dimensioning of the vertical axis turbine, we started from the power equation 1, which depicts the amount of energy that can be extracted from a flow [6].

1 3 푃 = 휂휌퐴 푉 (1) 2 푇

Where P represents the desired power to be extracted from the flow, η the efficiency (which is 35 % for this type of turbine [7], ρ the flow density, 퐴푇 is the cross-sectional area of flow that gets across the turbine and V is the flow speed that goes through the turbine. Considering that an amount power of 2 Kw is desired, a cross-sectional area of 1.39 m2 will be required, given a 2 m/s speed, with a density of 1024 Kg/m3. According to the required area, the length and diameter dimensions would be 1.4 m and 1 m respectively which will give us an aspect relation of 1.4:1, as it showed in Table I. Having the previous measurements, the next step is designing the blades. For this, a profile of 0021 from NACA series is chosen which is a symmetric profile, these are recommended for low power generators [8]. These show an increased lift coefficient at low attack angles; this is advantageous because the turbine makes use of the lift forces to generate its angular movement; in simulations shows a good start performance [9]. Solidity is a non-dimensional parameter, which is related with the blade width and is defined as the ratio between the total area of the blades and its swept area [4], since solidity is directly proportional to the force of the turbine (as shown in equation 2) and at the same time to torque; a higher solidity will imply a higher torque [10], however, values close to 1 should be considered, these will cause the turbine to close as if it was cylinder and its efficiency will drop. A solidity of 0.8 is chosen.

1 2 F = C ρσAV (2) 2 d

Solidity will be used for the calculation of the chord, which could be defined as the straight line drawn from the frontal to the rear part of the winged profile, this is calculated solving equation 3 of solidity σ [11][12], which is related to the number of blades n, the chord c, the turbine diameter D and π is a constant.

nc σ = (3) Dπ

GEET-19, Paris, 24-26 July 2019 Pag. 84 For the torsion calculation of the blades, the equation 4 [13]:

2

2 퐿 퐿 푇 = √1 + ( 푐푛 ) ⋅ sin (4) 휑 ( ) 0 2휎 푐푛 퐿 ( ) 2휎 퐿 휑0( 푐푛 ) 휑표( ) [ ( ( 2휎 ))]

This equation will be employed to maximize the torque, and considering the torsion angle φo as variable, we proceed to derive the function with the purpose of making it close to zero, then using the Newton- Raphson method, after three iterations we get the torsion value for which the torque is maximum. Given this, we continue to calculate the blade elevation angle, based on equation 5; which considers the length, radius and torsion angle [13].

퐿 tan⁡(훿) = (5) 휑0푅 Making the dimensioning of the turbine, we produce a prototype made out of nylon with a 3D printer, which is a 10% scaled down version of the turbine previously dimensioned. As the previous one, this turbine is designed with the SolidWorks software, then we proceed to experimentation.

Table I. Constant parameters ______Turbina Escala 1:1 Turbina Escala 1:10

Diameter (m) 1 0.1 Length (m) 1.4 0.14 Cross-sectional area (m2) 1.4 0.014 Radius (m) 0.5 0.05 Number of blades 3 3 Solidity 0.8 0.8 Blade profile NACA 0021 NACA 0021 Blade chord (m) 0.8377 0.08377 Attack Angle (°) 15 15 Elevation angle (°) 66 66 Torsion angle (°) 71 71 Power (W) 2007 2.4 Angular velocity (rpm) 38 190 ______

2.2. Gear train and generator – For the calculation of the gear train, it is considered the value of the angular velocity estimated at table I for the 1:1 scale turbine and the required angular velocity for the generator (180 rpm), which is capable of producing up to 2 kw. With this data we make use of equation 6 [14], which considers angular velocity, and the number of gear teeth. A calculation for 2 gears and 2 pinions of 62 and 22 teeth respectively is made, getting from this relation the required revolutions for the chosen generator; the gears are manufactured with onyx aided with a 3D printer.

휔퐺 푅푝 푁푝 = = (6) 휔푃 푅퐺 푁퐺

GEET-19, Paris, 24-26 July 2019 Pag. 85 2.3. – Experimental setup – For the experimentation described below we based on the Rachmat et al work, in which an engine was adapted in a turbine and a flow was produced to measure the torque [15]. An aluminum structure is designed, as Image 1 shows, in which the turbine is fixed vertically, on the upper part of the structure a torque Futek® sensor is placed, this has a maximum sensing capacity of 20 Nm and above this component there is an engine. Both are attached to the axis of the turbine; the engine varies the revolutions per minute of the turbine when flux goes through and test its performance at every assay. For the visualization of the torque a software developed by Madgeteck was used, it includes a data logger to read the data on the torquemeter. In Image 2 and 3 we can see an image showing the test in the varying slope water channel, and the prototype itself.

Image 1 Image 2 Image 3

Table II. Experimental set up in the channel. ______Parámetros Valores______

Current Velocity (m/s) 1 Temperature (°C) 30 Water density (Kg/m3) 996 Channel height (m) 0.3 Turbine material Nylon ______

For the graphs obtained after the assay the calculation of tip speed ratio TSR, drag coefficient 퐶푑, torque coefficient 퐶푇 and power coefficient 퐶푝. was required. For the calculation of efficiency 퐶푝⁡equation 7 is used, which considers the power generated by the turbine divided by the power that can be extracted from the water flow [16]. Tip speed ratio TSR is a non-dimensional parameter, this one relates the speed at the tip of a blade with the speed of the water flow [17], and it is calculated with equation 8. The drag coefficient 퐶푑 is calculated trough the relation of the velocity and efficiency, divided by the solidity, the radius and the angular velocity, as its shown at equation 9 [10]. To obtain torque coefficient 퐶푇 it was required to use equation 10 which relates the power coefficient with tip speed ratio [18].

푃푡 푐푝 = (7) 푃∞

휔⋅푅 푇푆푅 = 훥 = (8) 푉

푉휂 푐 = (9) 푑 휎푅휔

퐶푝 퐶 = (10) 푇 푇푆푅

GEET-19, Paris, 24-26 July 2019 Pag. 86 3. Results and Discussion – In order to obtain the efficiency vs tip speed ratio chart 퐶푝⁡vs ∆, as its shown at Image 4, a water flow with a speed of 1 m/s was calculated in the channel with a variable slope, it was employed a rotating turbine powered with a motor that varied at different angular velocity; the parameters at which the engine was shifted were 40, 60, 80, 100, 120 and 140 rpm respectively. In Image 4 we can see this chart, where in the x axis we have ∆ and in the y axis 퐶푝, then we obtain for this chart a maximum efficiency of 21.18 % when the turbine rotates at 100 rpm. This results are compared with the maximum efficiency at the same speed of the water current for a bigger turbine with solidity of 0.5, resulting in a smaller efficiency compared with the one presented here [9], we can conclude then that a higher solidity harnesses better the water flow, which can be considered as a better performance. For assays the turbine was rotated without load for the purpose of measuring its angular velocity, which was 160 rpm, a similar valued compared to the previously calculated parameter at Table I (190 rpm). In Image 5 we got the graph of power vs tip speed ratio P vs ∆, with a maximum value of 1.48 w, when ∆ is 0.52 In figure 6 we can see the drag coefficient vs tip speed ratio graph 퐶푑 vs ∆; this indicates that the resistance to movement decreases when ∆ increases, which also creates a bigger lift. The Torque Coefficient vs tip speed ratio graph 퐶푇 vs ∆ that we have at Image 7 shows that the torque coefficient diminishes when ∆ increases, we can find this relation in the works of Parag, in which an experiment with 3 turbines, larger than the one we employed results in a similar coefficient [17].

Turbine Turbine

25% 20% 21% 20% 2.00 18% 18% 1.39 1.48 1.42 20% 1.26 1.26 13% 1.50 15% 0.92 p

C 1.00

10% P (w) 5% 0.50 0% 0.00 0 0.2 0.4 0.6 0.8 0 0.2 0.4 0.6 0.8 ∆ ∆

Image 4 Image 5

Turbine Turbine

1 0.8

0.8 0.6 0.6 T d C C 0.4 0.4 0.2 0.2 0 0 0 0.2 0.4 0.6 0.8 0 0.2 0.4 0.6 0.8 ∆ ∆

Image 6 Image 7

The power calculated for the prototype (Table I), results larger than the maximum value obtained in the assays. However, this is because the calculation was made with the maximums efficiency for this kind of turbines, if we recalculate for an efficiency of 21% this results in 1.4 w, which matches with the value previously obtained at the tests. When making the tests, we noticed that there were some losses as a result of the collision between the flow and the structure. This let us conclude that a fully developed turbine will have an equal or even superior efficiency compared with the one obtained at the prototype.

GEET-19, Paris, 24-26 July 2019 Pag. 87 4. Conclusions – Given de data here presented, we estimate that the final turbine will make better use of the water currents, compared to other similar models, this is due to its geometry here described. Moreover, it is believed that the device will produce a higher torque, thanks to the solidity used in the blades and the torsional angle employed, which was fundamental for a maximum torque. The production of the final turbine is made with a 3D printer. Analyzing the real environment of use it is recommended to build the pieces in aluminum for properties such as corrosion resistance [19], and low density; the material is suitable for maritime environments and would contribute to a better performance. Considering the efficiency losses caused by many factors such as material or generator, it is estimated that the final version of the turbine will generate 1.7 kw.

4. References

[1] J. Mart, M. R. Hern, and V. Carrillo, “Simulación computacional de fluidos en micro turbina eólica de eje vertical tipo helicoidal,” 2015. [2] U. N. A. Rodriguez Rivas, Carlos. Plataforma and D. E. E. Marítima, “Análisis de viabilidad económica de una plataforma de energía marítima,” 2010. [3] U. P. de C. Graciá Ribes María, “Estudio de las diferentes formas de conseguir energìa con el mar y su aplicabilidad en el litoral español,” 2014. [4] P. Bachant and M. Wosnik, “Performance measurements of cylindrical- and spherical-helical cross-flow marine hydrokinetic turbines, with estimates of exergy efficiency,” Renew. Energy, vol. 74, pp. 318–325, 2015. [5] A. N. Gorban’, A. M. Gorlov, and V. M. Silantyev, “Limits of the Turbine Efficiency for Free Fluid Flow,” J. Energy Resour. Technol., vol. 123, no. 4, p. 311, 2001. [6] B. A. Gorlov, “Helical Turbine and Fish Safety,” pp. 1–14, 2010. [7] J. Anderson, B. Hughes, C. Johnson, and N. Stelzenmuller, “Capstone Project Report: Design and Manufacture of a Cross-Flow Helical Tidal Turbine,” p. 143, 2011. [8] A. J. Gonz, L. Jos, G. Coronado, and Y. E. Gonz, “Selección del perfil alar simétrico óptimo para un aerogenerador de eje vertical utilizando la dinámica de flujos computacional Selection of the optimal symmetrical , for a vertical axis using computational fluid dynamics,” no. 22, pp. 83–91, 2017. [9] R. Keough, V. Mullaley, H. Sinclair, and G. Walsh, “Design , Fabrication and Testing of a Water Current Energy Device,” pp. 1–75, 2014. [10] E. Abril, “Diseño de una turbina para una pico central hidroeléctrica para las condiciones del rio Vaupés en Mitú,” pp. 1–66, 2016. [11] P. E. Conference, P. Engineers, M. Shiono, K. Suzuki, and S. Kiho, “Output Characteristics of Darrieus with Helical Blades for Tidal Current Generations,” vol. 3, pp. 859–864, 2002. [12] S. Pongduang, C. Kayankannavee, and Y. Tiaple, Experimental Investigation of Helical Tidal Turbine Characteristics with Different Twists, vol. 79, no. November 2015. Elsevier B.V., 2015. [13] D. A. Gorlov, “Development of the Helical Reaction Hydraulic Turbine,” 1998. [14] R. L. Mott, Diseño de Elementos de Maquinas. 2006. [15] R. Firdaus, T. Kiwata, T. Kono, and K. Nagao, “Numerical and experimental studies of a small vertical-axis wind turbine with variable-pitch straight blades,” vol. 10, no. 1, pp. 1–15, 2015. [16] M. A. Al-dabbagh and M. I. Yuce, “Simulation and Comparison of Helical and Straight- Bladed Hydrokinetic Turbines,” no. March, 2018. [17] P. K. Talukdar, V. Kulkarni, and U. K. Saha, “Field-testing of model helical-bladed hydrokinetic turbines for small- scale power generation,” Renew. Energy, vol. 127, pp. 158–167, 2018. [18] S. Kaprawi, D. Santoso, and R. Sipahutar, “Performance of Combined Water Turbine Darrieus- Savonius with Two Stage Savonius Buckets and Single Deflector,” vol. 5, no. 1, pp. 1–5, 2015. [19] E. Mario, “Diseño de un rotor eólico tipo Darrieus helicoidal,” vol. 1, no. 2, pp. 34–41, 2017.

GEET-19, Paris, 24-26 July 2019 Pag. 88