Power Performance Assessment of Vertical-Axis Tidal Turbines Using an Experimental Test Rig

energies Article Power Performance Assessment of Vertical-Axis Tidal Turbines Using an Experimental Test Rig Aitor Fernández-Jiménez 1,* , Eduardo Álvarez-Álvarez 1, Mario López 2, Mateo Fouz 3, Iván López 3, Ahmed Gharib-Yosry 4, Rubén Claus 2and Rodrigo Carballo 3   Citation: Fernández-Jiménez, A.; Álvarez-Álvarez, E.; López, M.; Fouz, M.; López, I.; Gharib-Yosry, A.; Claus, R.; Carballo, R. Power Performance Assessment of Vertical-Axis Tidal Turbines Using an Experimental Test Rig. Energies 2021, 14, 6686. https://doi.org/10.3390/ en14206686 Academic Editor: Elhoussin Elbouchikhi Received: 7 September 2021 Accepted: 12 October 2021 Published: 15 October 2021 Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/). 1Hydraulic R&D Group, University of Oviedo, Gonzalo Gutiérrez Quirós St, 33600 Mieres, Spain; [email protected] 2DyMAST R&D Group, Department of Construction and Manufacturing Engineering, University of Oviedo, Gonzalo Gutiérrez Quirós St, 33600 Mieres, Spain; [email protected] (M.L.); [email protected] (R.C.) 3 Área de Ingeniería Hidráulica, University of Santiago de Compostela, Benigno Ledo St, 2, 27002 Lugo, Spain; [email protected] (M.F.); [email protected] (I.L.); [email protected] (R.C.) 4 Mechanical Power Department, Port Said University, Port Said 42526, Egypt; [email protected] *Correspondence: [email protected] Abstract: This article presents the characteristic curves of a vertical-axis hydrokinetic tidal turbine of the Darrieus subtype aimed at meeting the electricity demand of port facilities located at harbors and estuaries with low water-speed conditions. The turbine was tested in the water-current flume of the University of Santiago de Compostela for several flow conditions with different water heights and water speeds. Blockage conditions were tested by examining the results from two groups of tests: with and without an accelerator device that restricts the flow around the rotor. The tip speed ratio and the power coefficient were used to characterize the performance of the turbine for each test. Finally, the results for open-field conditions were obtained by applying empirical expressions, which allowed us to assess the performance of the device in estuaries and harbors with known water-flow regimes. Keywords: marine renewable energy; tidal energy; laboratory tests; tidal turbine 1. Introduction Port facilities utilize a series of key services which promote the economic and social development of regions. There are a number of key activities which demand significant consumption of energy, existing important activities that demand big amounts of energy. Therefore, improvements in energy consumption are essential in order to ensure the most favorable use of electricity. It should be borne in mind that, due to their nature, there are diverse renewable hypo-carbon energies that could be used at these facilities, such as tidal currents or wave energy [1]. The total amount of energy consumed at ports represents 75% of all the energy used by the transport sector worldwide, which equates to 221 million tons of oil equivalent (MTOE) [ 2 ]. This represents 2.7% of global CO 2 emissions. Unless drastic measures are taken to reduce these emissions, it is expected that in the coming years, as a consequence of the growth of international trade, they will continue to increase, reaching 7% of all CO 2 emissions by 2040 [2]. The International Maritime Organization (IMO) has proposed different measures to reduce today’s emissions by 50% up to 2050. Some such proposals include the use of renewable fuels, such as hydrogen or ammonia, or electric generation from low-carbon resources [ 3 ]. It should be noted that even some private concerns, such as the Maersk shipping company, have set a zero-emissions target for the year 2050 [4]. In the case of ports, the incorporation of electricity generation technologies from low-carbon resources will be a reality in 2030 [ 5 ]. Some examples are the port of Houston Energies 2021,14, 6686. https://doi.org/10.3390/en14206686 https://www.mdpi.com/journal/energies Energies 2021,14, 6686 2 of 12 (United States), where, in 2020, a solar panel farm was set up to supply electricity for harbor activity [ 6 ] or the port of Gothenburg (Sweden), where 80% of the electricity generated comes from a mixed plant of solar panels and windmills [7]. In this context, the development and implementation of other electric generation technologies from renewable resources, such as vertical-axis hydrokinetic turbines to harness water currents, are an excellent alternative to provide green energy to port facilities and reach the zero-emissions target set by the IMO for 2050. These devices are low-pressure turbines which generate electricity from the kinetic term (velocity) of a flowing stream with low potential energy, where the difference between upstream and downstream height is less than 0.2 m [ 8 ]. The generated energy can be used directly in the port, optimizing the environmental benefits of renewable energy, constituting a nearby point of consumption, and thereby reducing transport, operation, and maintenance costs. The use of hydrokinetic turbines for their application in harbor areas is based on the electrical generation from the energy of tidal currents. These streams are a renewable, clean, and largely foreseeable resource which could be used to contribute to the decarbonization of port facilities. It has been estimated that tidal currents have a global energy potential of 800 TWh/year [ 9 ], meaning that their development and implementation would have a considerable niche in the world market. To date, this technology has been studied from various perspectives. These include: energy storage [ 10 ], their effects in the electrical market [ 11 ], and their economic optimization of operation [ 12 ]. Additionally, on many occasions, auxiliary mechanical elements, such as flow accelerators, need to be installed to increase velocity and ensure satisfactory operation [13]. Hitherto, the two main types of hydrokinetic turbines used to supply energy to harbor areas have been axial turbines (flow and axis parallel) and crossflow turbines (flow and axis perpendicular). Some of the real projects carried out using axial turbines include the RITE system based in New York (USA) that can generate 70 MWh/year [ 14 ] and the Kvalsund project installed in Finnmark (Norway) with a total power of 500 kW [ 15 ]. As regards crossflow rotor-based projects, the most noteworthy are the Blue Energy system settled in Canada with a total power of 2 kW [ 16 ] and the Oryon Watermill project located in the Netherlands with 150 kW of power [ 17 ], both with diameters greater than 2 m. In addition, the OceanQuest project should be highlighted, in which the testing of a 1 MW H-Darrieus prototype was carried out on the coast of Brittany (France) [ 18 ]. All of them are still in an embryonic stage, but recent results have been successful, which augers well for the future. From a scientific point of view, the study of crossflow turbines has been based mainly on the development of numerical models. This is due to existing difficulties in the experimental characterization of these devices, such as problems with starting rotation at low water velocities and obtaining measurements of their power stage since, at these low speeds, torque is small [ 19 ]. In this sense, some studies stand out, such as the study carried out by Behroui et al. [ 20 ] where a Savonius rotor was characterized using Computational Fluid Dynamic (CFD) techniques, obtaining the mechanical performance. Another study of interest is the work of Kumar and Saini [ 21 ], which presents a numerical methodology for optimizing the blade size of a Savonius rotor taking into account the hydrodynamic parameters of the water stream. As regards Darrieus typology, the studies presented by Kirke [ 22 ] and Benchikh Le Hocine et al. [ 23 ] stand out. Both publications considered the influence of the blockage effect on the power stage of a crossflow turbine. Thus, in the initial study, a deflector was numerically tested to observe its influence on the energy production of a rotor, while in the second publication, two numerical models were used together (CFD and Blade Element Method, BEM) to study the influence of the position of the water-free sheet on the turbine’s performance. This article presents an analysis and study of a Darrieus vertical-axis crossflow turbine in a current water flume, under blocking effects induced by a flow-accelerator device and with low water-velocity conditions (<1 m/s). The study was carried out in a water-current Energies 2021,14, 6686 3 of 12 flume installed in the laboratory of the Hydraulic Engineering Area at the Polytechnic School of Engineering of Lugo (University of Santiago de Compostela, USC). 2. Materials and Methods 2.1. Experimental Set-Up A Darrieus rotor (Hydraulic R&D Group, Mieres, Spain) design was used, geometrical characteristics of which are shown in Figure 1. Figure 1. Darrieus rotor design tested, with the geometrical characteristics. Using Equation (1), the solidity of this design was 0.34. This solidity value was selected because, since the rotor would work under low water-velocity conditions, having a high solidity made it easier to start the rotation of the turbine [ 24 ]. In fact, some other studies such as Singh et al. [ 25 ] and Mohammed et al. [ 26 ] show that Darrieus rotors improve their starting capacity when solidity is high since the static torque is increased, although there may be a downturn in performance. σ=N·c 2·π·R(1) where σ is the solidity (dimensionless), N is the number of blades (dimensionless), c is the blade chord (m), and Ris the rotor’s radius (m). The rotor was printed using 3D additive technology, attending to lift designs criteria to maximize performance. The material used was nylon filament (PA) due to its resistance against stresses, lightness, and non-degradation under water presence. The surface roughness of the elements that made up the rotor was reduced due to the high precision and excellent final finish of the printing equipment used (BCN3D SigmaxR19 ) (BCN3D, Barcelona, Spain). Additionally, after printing, the elements were subjected to sanding and sealing primer with plastic material. In consequence, the surface was completely smooth and there were no effects of turbulence in the flow due to the roughness of the turbine. The turbine was placed in a stainless-steel metal case, in the lower part of which a radial bearing was located to ensure rotation, and the rotors were coupled to a 1.20 m-long stainless-steel shaft with a diameter measuring 10 mm. This structure was 0.62 m wide, 0.70 m long, and 1 m high, and could be found inside the water-current flume (Figure 2). In addition, two lateral windows of 0.24 m 2 were built so that the tests could be filmed. Energies 2021,14, 6686 4 of 12 The turbine was located at 0.175 m above the floor of the channel, using the lower rotor cover as a reference. Figure 2. Structure used to couple the different rotor designs. The structure was designed in such a way that it was possible to add and easily remove the rotor and the accelerator device by using two sets of sprats. The accelerator device was based on a symmetrical NACA profile (0015), being placed on each side of the turbine. This element allowed the incident speed to be increased by 1.6 times as the cross-section was reduced (from 0.65 m to 0.39 m). Furthermore, the use of this standardized profile for the accelerator allowed a gradual reduction in the passage section, avoiding detachments and directing the flow of water towards the rotor, minimizing the appearance of turbulent effects [ 27 ]. Figure 3shows the accelerator device (Hydraulic R&D Group, Mieres, Spain) geometry. Figure 3. Geometry of the accelerator device. Energies 2021,14, 6686 5 of 12 A 20 m-long water-current flume, which was 0.65 m wide and 1 m high with transparent crystal walls, was used to carry out the power characterization of the selected rotor. The movement of the water inside the channel was forced by a hydraulic pump (Pedrollo Spa, San Bonifacio VR, Italy) with 22 kW of power (IDEAL-RNL-150-250/264 series). There was no free-discharge point on the circuit (Figure 4). The flow inside the channel was controlled by means of several valves that were calibrated with an ultrasonic flow meter. In this manner, different water speeds and heights could be reproduced under controlled conditions. Figure 4. Water-current flume with the hydraulic control devices. To monitor the hydraulic parameters during the tests such as flow rate, water velocity, and height variation, different sensors were used, including water-pressure devices, electric and ultrasonic water-height sensors, and an Acoustic Doppler Velocimeter (ADV) (Sontek, San Diego, CA, USA). The ADV was located just at the inlet zone of the turbine so that the speed upstream the rotor was known (Figure 5). The power characterization simulation process of the rotor was based on variation of the resistive torque of the electric brake. In this way, for each value induced on the electric brake, a rotational speed (n) and mechanical torque (T) in the turbine were obtained, so that the rotor’s power stage could be fully and satisfactorily simulated. Figure 6shows an outline of the procedure adopted. The characterization began with an initial measurement, called “no load”, where the electric brake did not offer resistance so that the turbine would rotate at its maximum speed but would not generate mechanical torque (Step 1). Then, by using a power-supply unit, the brake voltage was progressively varied, increasing the braking torque and obtaining, for each voltage, an nand T value of the turbine (Step 2). The test ended when the turbine stopped its rotation (Step 3). The measurements were recorded with a high-precision torque meter (MAGTROL TS 103 series, Bufalo, NY, USA) and an electrical brake (MAGTROL HB140M series, Bufalo, NY, USA). This system measured instant torque (N · m) and rotational speed (rpm), storing and processing all data using the manufacturer’s software (MAGTROL TORQUE Version 10, Bufalo, NY, USA). Energies 2021,14, 6686 6 of 12 Figure 5. Hydraulic sensors’ locations during tests. Figure 6. Scheme of the experimental characterization procedure. 2.2. Performance Indicators The power characterization of the turbine was carried out from the raw data obtained by the torque meter. Subsequently, these data were processed, and several performance indicators were obtained. The relation of the speed at the tip of the blade with the upstream flow speed—i.e., the Tip Speed Ratio (TSR), was obtained by applying: TSR =ω·R U(2) Energies 2021,14, 6686 7 of 12 with ω the rotational speed of the rotor (rad/s), and Uthe average water speed in the uniform zone upstream of the turbine at the channel (m/s). In addition, the power coefficient of the turbine (dimensionless) was obtained by applying: Cp =Pm Pt(3) where Pm is the average mechanical power in a complete turn of the turbine (W), obtained using: Pm=ω·T(4) with T the torque force at the shaft (N · m), and Pt the average hydraulic power of the current during a given test (W), obtained via: Pt=1 2·ρ·At·U3(5) with ρ the fluid density (kg/m 3 ), and At the area swept by the blades (m 2 ), which was calculated by multiplying the height of the turbine ( h ) by the transverse width occupied by the rotor (2R). 2.3. Open-Field Conditions In the presence of any blockage, the fluid is forced to pass through the rotor and the flow is accelerated [ 28 ]. This may be applied to experiments in current flumes such as those carried out in this study, but not to estuaries where the flow is unrestricted by walls near the rotor—i.e., open-field conditions. On this basis, the performance of the turbine for open-field conditions is required to assess the power production of the turbine for a possible location of interest in the vicinity of a port. To characterize the performance of the turbine in open-field conditions, the modified values of Cp and TSR were obtained by applying: CpF=CpU UF3 (6) TSRF=TSRU UF(7) respectively, where UF is the open-field flow speed corresponding to each value of water speed tested in the flume under blocked conditions (m/s). The relation between velocities, applied velocity relationship, is based on the expression proposed by Gauvin et al. [ 28 ]. This Equation (8) was selected since it performs well with blockage coefficients greater than 20%, as indicated in the case studied [28]: UF U2 =1 1−(m·B)(8) where B(dimensionless) is the blockage coefficient provided by: B=2·R·h b·y(9) with hthe height of the rotor (m), bthe flume width (m), ythe water height at the flume (m), and man empirical coefficient provided by [28]: m=8.14·B2−7.31·B+3.23 (10) A complete power characterization of the rotor was carried out to study the influence of blockage phenomena during the power stage of the turbine under low water-speed Energies 2021,14, 6686 8 of 12 conditions (<1 m/s). For each velocity, three water heights (0.5, 0.55, and 0.6 m) were induced with and without an accelerator device, carrying out a total of six tests. Under these conditions, values of B of 0.21, 0.23, and 0.25 were obtained for the tests without an accelerator and of 0.38, 0.41, and 0.45 when this element was placed. In both cases, the variation of blockage was made by varying the water level, although each situation was different because of changes in the hydrodynamic behavior due to the presence of the accelerator. Figure 7shows an outline of the tests carried out. Figure 7. Outline showing how the tests were carried out. An initial water velocity of U1 = 0.174 m/s was selected because it was the cut-in speed of the rotor, while U3 = 0.246 m/s was the maximum water speed of the current flume. U2 = 0.210 m/s was selected as an intermediate speed. It should be noted that the position of the rotor, water-height sensors, and the ADV speed meter were invariably located in the same place during all tests. More specifically, all sensors were located at a height of 0.30 m from the bottom of the channel, and the ADV located at 0.70 m upstream of the rotor axis. Figure 8shows a diagram with the main distances. Figure 8. Diagram with the position of the different elements used in tests. 3. Results and Discussion The proposed Darrieus water-current turbine was studied following the power characterization procedure described in previous sections. On the one hand, Figure 9shows the variation of Pm versus the rotational speed of the turbine ( n ) for tests without an accelerator represented as being (a–c) whilst the different water speeds used were of 0.174, 0.210, and 0.246 m/s, respectively. On the other hand, Figure 10 shows the same tests but Energies 2021,14, 6686 9 of 12 using the accelerator device represented as (a–c) with water speeds of 0.174, 0.210, and 0.246 m/s, respectively. Figure 9. Variation of power vs. nfor the different tests without an accelerator device. Figure 10. Variation of power vs. nfor the different tests with an accelerator device. The characteristic curves show how power values were affected by variations in the current speed and in water height. More specifically, when the current velocity increased or the water height decreased, the power of the rotor went up. This produces an increase of blockage conditions, which has direct consequences on energy extraction. In fact, the effects of blockage were noticeable when the accelerator device was installed, producing a significant increase in the rotational speed and power, even under the low water velocities of the tests. It could also be observed how the recorded values reached a point of maximum power (maximum power point, MPP) and then the turbine stopped its rotation. The reason was that the left part of the Cp vs. TSR curve corresponds to an unstable operating zone of the turbine caused by the loss of lift force on blades (the turbine was the Darrieus NACA 0015 type that works mainly under lift forces). In the case that the rotor had a design based on drag forces (for example a Savonius type) it would be possible to obtain those points of the instability zone. [29]. Using the dimensionless coefficients of TSR and Cp it was possible to obtain an initial approximation of the operation of the designed turbine. Figure 11 shows these power parameters only with the results of the tests where the accelerator device was used since these tests conditions allow for a better visualization of the behavior of the turbine. For the three water velocities, maximum TSR values higher than 4.5 were attained, which indicates a significant acceleration of the speed at the tip of the blade with respect to the current velocity. Furthermore, in all tests the maximum Cp values were greater than 1 due to the great blockage existing during the tests (B > 0.38). Therefore, it is evident that the results obtained were strongly influenced by blockage phenomenon but were not consistent with the behavior of the turbine under open-field conditions. To