Method for reducing the dynamic forces and vortical structures acting on a wind turbine

Abstract

A method for reducing the dynamic forces and vortical structures acting on a wind turbine is proposed. The method comprises providing a suction pump at a given location of a wind turbine; generating a constant or periodic suction of the fluid at a blade tip by the suction pump; and decreasing a pressure at the blade tip until reducing the vortices up to a given threshold. The invention also relates to a method and computer program for maximizing the power generated by a wind turbine.

Description

[0001] METHOD FOR REDUCING THE DYNAMIC FORCES AND VORTICAL STRUCTURES ACTING ON A WIND TURBINE

[0002] TECHNICAL FIELD

[0003] The present invention generally relates to solutions for wind turbines. More specifically, the invention relates to a method for reducing the dynamic forces and vortical structures acting on a wind turbine. Two different strategies / solutions are proposed with that objective, which can be independent or complementary to one another.

[0004] BACKGROUND OF THE INVENTION

[0005] In recent years, the number of installed large-scale horizontal axis wind turbines (HAWT) is rapidly growing around the world. These modern large-scale wind turbines are affected by the unsteady characteristics of the wind turbine operating environment such as turbulence, atmosphere stability, gusts, misaligned yaw of the rotor, wind shear, tower shadow, and upwind wakes, all these phenomena will cause the local inflow angle on the rotor blades to change periodically or suddenly. As a result, several airfoils may operate under stall conditions, resulting in boundary layer separation, drag increase and loss of lift and ultimately the deterioration of the airfoils aerodynamic performance and a reduction of wind turbine power output.

[0006] Flow control can reduce such unfavorable behavior by delaying the boundary layer separation. Flow control techniques can generally be categorized as passive and active. Passive flow control techniques do not require input energy and usually utilize geometrical modifications like changing the shape of turbine airfoils via implementing Gurney flaps [2-5] or vortex generators [1 , 6-8] for example, but their performance is optimal only at design conditions. Wind turbine blades must face a variety of flow conditions during operation, thus it is doubtful whether such control strategies can have a positive effect under all operating conditions. In contrast, Active flow control (AFC) techniques (requiring power input) remove the drawback associated with passive control at off-design conditions [9-11], This form of control has been associated mainly with adding / subtracting momentum to / from the mean velocity in particular locations in order to interact with the boundary layer and delay / promote its separation, synthetic jet actuators (SJA) [11-15] or direct-barrier-discharge plasma actuators (DBD-PA) [16-18], are two of the main devices used in AFC implementations. Further information on flow control actuators can be found in review articles by Cattafesta and Steplak

[0019] , Johnson et al

[0040] and Barias and Vonkaila

[0021] ,

[0007] Among the different AFC techniques, synthetic jet actuators (SJA) have the advantage of generating pulsating flow which combines suction and blowing phases, and they can have associated large momentum coefficients. Their high pulsating velocity, fast time response, low power consumption and easiness of integration to the surface of the body, enables their practical use and maximizes the possible aerodynamic gains. Due to the fact that SJ’s control is a mature flow control technique of AFC, which is well established and has achieved good results in the field of aeronautics, it is logic to think that their application can be expanded as an effective means to control flow separation on HAWT blades. The application and influence of SJA-AFC on the aerodynamic performance of rotating blades of wind turbines cannot be traced until almost the last decade. Stalnov et al.

[0022] and Maldonado et al. [23, 24] were among the first researchers to illustrate experimentally the influence of SJA on the aerodynamic performance of HAWTs.

[0008] Stalnov et al.

[0022] experimentally investigated the performance of SJ based-controlled blades on a modified IAI pr8-SE profile at Reynolds numbers ranging between 2 x 105< Re < 8 x 105using a Hot-wire sensor. They demonstrated an increase of the overall efficiency of about 5-15%. In Maldonado et al. [23, 24], they initially applied SJs to an S809 airfoil in order to control the unsteady loading for a non-rotating and fixed-pitch wind turbine blade. They used 14 SJ’s in total, distributed in three locations, near the root, around the middle, and at the blade tip. A set of Reynolds numbers ranging between 7.1 x 104, and 2.38 x 105were considered, the momentum coefficient range employed was between Cp = 1.34 x 10'3and Cp = 5.97 x 1O’3.

[0009] In the research undertaken by Taylor et al.

[0025] , they studied experimentally a Wind Turbine (WT) blade using an S809 airfoil across the entire blade. They incorporated two arrays of 12 SJ’s placed along the blade and at the normalized chord locations of x / C = 0.1 and x / C = 0.2, measured from the leading edge. A lift to drag ratio hysteresis reduction of around 73% and a significant reduction in flow separation near the leading edge, was found in the numerical simulations performed by Tran et al.

[0026] on the NREL VI turbine with the S809 airfoil shape at Re = 1 x 106and Cp = 0.003 when the blade was undergoing pitching. To investigate the effect of the SJ on the delay of the boundary layer separation and its associated vortical structures, WT studies performed by Rice et al.

[0027] using stereo Particle Image Velocimetry (PIV) and employing the S817 airfoil type were conducted at Re = 3.75 x 105, the SJ was placed at x / C = 0.35, being the momentum coefficient Cp = 0.012. When AFC was applied, they observed a significant reduction in the hysteresis of the lift and pitching moment coefficient of 41% and 85%, respectively. A numerical study was conducted by Moshfeghi and Hur

[0028] on a wind turbine blade 5809 HAWT model for a range of angles of attack (AoA) between 0° and 25° at Re = 1 x 106under the excitation control of a SJ located at the chordwise location of x / C = 0.75. Three different SJ injection angles of 5°, 15° and 25° were considered.

[0010] Maldonado and Gupta

[0029] , experimentally investigated at Re = 1 x 106the performance of a real 3-bladed rotor based on the S809 airfoil model, they considered five turning speeds, 250, 500, 750, 1000 and 1250rpm, the pitching angles evaluated were 0°, 3°, and 6°. They placed 20 SJ actuators along each blade radius and they used laser Doppler Anemometry (LDA) to analyze the flow performance at the SJ’s exit. For an angular velocity of 500rpm and a pitch angle of 3°, the actuation of 20 SJs per blade decreased the power input to drive the rotor by 10.6%, the rotor efficiency increased by 28%.

[0011] Troshin and Seifert

[0030] studied experimentally at Reynolds numbers ranging from Re = 4 x 105to 5 x 1 o5, the performance of a thick airfoil AH93-W-300 model typically used at the root section of wind turbines. Three rows of synthetic jets consisting of 13 actuators in each row were considered, the momentum coefficients applied to each SJ were varied from zero (baseline case) to 0.024. A closed-loop control system was used to determine momentum coefficient applied at each instant, the feedback loop was based on the information provided by the pressure sensors and hot-film sensors placed on the airfoil upper surface. A significant increase in net-energy generated by the wind turbine, up to 60% when compared with the non-actuated performance, was obtained over a wide range of AoA and Reynolds numbers. Wu et al.

[0031] conducted a numerical study on an elliptic cross-section airfoil at Re = 1100, they placed a pair of SJs on the upper and lower surfaces and considered the combined effect of pitching and plunging motion. The effect of different parameters such as jet inclination, jet location and the phase angle between the jet and pitching motion on the energy production capability were analyzed. They observed an increase of up to 30% in energy harvesting along with a significant increase in lift force. The S809 profile employed in a small HAWT and considering suction to delay the boundary layer separation was studied by [32, 33], In the former research the suction groove was located at 0.15C versus the leading edge, four different wind speeds of 7, 10, 13 and 20m / s and 12 different momentum coefficients for each wind speed were considered. Improvements were particularly relevant at low radius sections, where large boundary layer separation exists. The power capability of the HAWT raised with the wind speed, obtaining an increase versus the base case of 225.56% for a wind speed of 20m / s. The effectiveness of using a single suction slot or two suctions slots combined with a blowing one was analyzed in

[0033] , Slots were located at 0.1 C and 0.5C for the dual case, just the former one was used for the single case. Two momentum coefficients were considered and they observed that when using the larger momentum coefficient and the dual suction configuration, the net power output had a maximum increase. Due to its high performance and efficiency, the SD7003 airfoil has been particularly investigated to be employed in small HAWTs and VAWTs (Vertical Axis Wind Turbines), since they operate at relatively low Reynolds numbers Re ~ 6 x 104, [34-37], One of the most comprehensive numerical studies to investigate the performance of SJs on this particular airfoil and Reynolds number was undertaken by Tousi et al. [11 , 12], They optimized at four AoA 4°, 6°, 8°, and 14° and using an in-house Genetic Algorithms optimizer, the five AFC parameters: jet position, jet with, momentum coefficient, frequency, and jet injection angle. They observed that aerodynamic efficiency could be increased by 251 % at AoA = 14° while a 39% increase was obtained at AoA = 8°. In fact, in the present invention the same optimizer and optimization process previously employed to obtain the five optimum AFC parameters for different airfoils cut along the DTU-10MW-RWT

[0038] is to be used.

[0012] Farhadi Asghar et al.

[0039] discloses a three-parameter numerical optimization of aerodynamic characteristics of an airfoil at high angle of attack. The separation region is reduced with two types of suction / blowing usage: single blowing jet and combined suction and blowing jets. Genetic algorithm is used to optimize the three parameters of the suction / blowing jet strength, jet angle and jet location. The objective of this investigation was to reduce the high separation region over the NACA0012 airfoil to enhance its aerodynamic characteristics. The time- averaged compressible Navier-Stokes equations with Spalart-Allmaras turbulence model are solved along with the optimization algorithm.

[0013] Ali Boudis et al.

[0040] provides a method in which synthetic jet parameters are optimized using response surface methodology (RSM) to improve aerodynamic performances of a S809 profile placed in a flow at a Reynolds number Re = 106 with an angle of attack of 15.2°. The synthetic jet actuator is positioned on the upper surface of profile at a distance 5% of the chord length from the leading edge. The considered optimization parameters are: jet amplitude, jet frequency and jet injection / suction angle. The RSMs are constructed based on a quadratic meta-model established from a BOX-BEHNKEN design. The genetic algorithm (GA) is used to explore the surface and determine the parameters giving best performances. Navier-Stokes equations are solved in unsteady calculation of a two-dimensional domain using ANSYS Fluent software. Turbulence is accounted using the k-w SST model.

[0014] Duvigneau et al.

[0041] deals with the optimization of the location of synthetic jet on the suction side of an airfoil to control stall. The optimal location is found by coupling a time-accurate flow solver with adaptive mesh refinement / coarsening techniques and an automatic optimization algorithm. The flow and jet are modeled by the unsteady Reynolds-averaged Navier-Stokes equations (LIRANS) with a near-wall low-Reynolds number turbulence closure. An unstructured grid refinement / coarsening method is used to automatically generate meshes adapted to the presence of the synthetic jet at a prescribed location. An optimization algorithm modifies the location of the synthetic jet to determine the best actuator location to increase the time-averaged lift for high angles of attack.

[0015] The primary distinction between existing solutions and the present invention lies in the fact that previous studies neither define nor consider the necessity of identifying boundary layer separation in a given airfoil section under varying operating conditions. Moreover, they do not recognize the need to optimize all Active Flow Control (AFC) parameters across different wind turbine operating conditions to maximize airfoil efficiency and / or lift, based on the vortex shedding frequency associated with boundary layer separation.

[0016] In contrast, the present invention establishes that, to maximize wind turbine power generation, it is essential to first determine the boundary layer separation point and its corresponding vortex shedding frequency for each airfoil (i.e. each blade section) along the blade span under all AoA and operating conditions. Once this information is obtained, the invention proposes optimizing all AFC parameters for each airfoil and then adapting (integrating) these optimized parameters along the entire blade span for all operating conditions of the wind turbine.

[0017] For instance, Farhadi Asghar et al.

[0039] numerically optimized three AFC parameters for a single airfoil at a fixed angle of attack (AoA). Ali Boudis et al.

[0040] positioned synthetic jet actuators (SJAs) at 5% of the chord using an S809 airfoil and optimized the jet injection / suction angle, frequency, and amplitude. Duvigneau et al.

[0041] focused on optimizing the location of the SJA for a given airfoil at stall AoA. Notably, all these studies considered only a single blade section, a limited number of AoA, and the optimization of one or a few AFC parameters. By contrast, the present invention proposes a comprehensive optimization of all AFC parameters across all blade sections, at all AoA, and under all operating conditions of the wind turbine. Once optimized, these AFC parameters must then be adapted (integrated) along the blade span for each operating condition, ensuring maximum efficiency and power generation.

[0018] Consequently, based on the existing prior art, it can be concluded that research on aerodynamic performance of HAWTs with SJs control has been generally done without a previous analysis of the boundary layer separation and its associated frequency.

[0019] New methods for maximizing the power generated by a wind turbine are therefore needed, specifically to enable the obtention of the boundary layer separation location and associated frequency for each airfoil.

[0020] Considering the reduction / minimization of the vortical structures generated on the blade tip via implementing constant / periodic suction (on the blade tip), not a single research work is found in the prior art.

[0021] References:

[0022] [1] C. Zhu, Y. Qiu, Y. Feng, T. Wang, and H. Li, “Combined effect of passive vortex generators and leading-edge roughness on dynamic stall of the wind turbine airfoil,” Energy Conversion and Management, vol. 251 , p. 115015, 2022.

[0023] [2] W. L. Zhao and X. L. Zheng, “Numerical simulation on the flow control by mounting indented gurney flaps for a wind turbine,” Advanced Materials Research, vol. 512, pp. 623-627, 2012.

[0024] [3] B. L. Storms and C. S. Jang, “Lift enhancement of an airfoil using a gurney flap and vortex generators,” Journal of Aircraft, vol. 31 , no. 3, pp. 542-547, 1994.

[0025] [4] D. Baldacchino, C. Ferreira, D. D. Tavernier, W. Timmer, and G. Van Bussel, “Experimental parameter study for passive vortex generators on a 30% thick airfoil,” Wind Energy, vol. 21 , no. 9, pp. 745-765, 2018.

[0026] [5] J. Wang, Y. Li, and K.-S. Choi, “Gurney flap — lift enhancement, mechanisms and applications,” Progress in Aerospace sciences, vol. 44, no. 1 , pp. 22-47, 2008.

[0027] [6] M. Kerho, S. Hutcherson, R. Blackwelder, and R. Liebeck, “Vortex generators used to control laminar separation bubbles,” Journal of aircraft, vol. 30, no. 3, pp. 315-319, 1993.

[0028] [7] A. Abbas, J. De Vicente, and E. Valero, “Aerodynamic technologies to improve aircraft performance,” Aerospace science and technology, vol. 28, no. 1 , pp. 100-132, 2013.

[0029] [8] J. C. Lin, “Review of research on low-profile vortex generators to control boundary-layer separation,” Progress in aerospace sciences, vol. 38, no. 4-5, pp. 389-420, 2002.

[0030] [9] G. Jones, M. Santer, and G. Papadakis, “Control of low reynolds number flow around an airfoil using periodic surface morphing: A numerical study,” Journal of Fluids and Structures, vol. 76, pp. 95-115, 2018.

[0031]

[0010] A. Rezaeiha, H. Montazeri, and B. Blocken, “Active flow control for power enhancement of vertical axis wind turbines: Leading-edge slot suction,” Energy, vol. 189, p. 1 16131 , 2019.

[0011] N. Tousi, M. Coma, J. Bergada, J. Pons-Prats, F. Mellibovsky, and G. Bugeda, “Active flow control optimization on sd7003 airfoil at pre and post-stall angles of attack using synthetic jets,” Applied Mathematical Modelling, vol. 98, pp. 435-464, 2021.

[0032]

[0012] N. Tousi, J. Bergad'a, and F. Mellibovsky, “Large eddy simulation of optimal synthetic jet actuation on a sd7003 airfoil in post-stall conditions,” Aerospace Science and Technology, p. 107679, 2022.

[0033]

[0013] M. Amitay, D. R. Smith, V. Kibens, D. E. Parekh, and A. Glezer, “Aerodynamic flow control over an unconventional airfoil using synthetic jet actuators,” AIAA journal, vol. 39, no. 3, pp. 361-370, 2001.

[0034]

[0014] N. Buchmann, C. Atkinson, and J. Soria, “Influence of znmf jet flow control on the spatiotemporal flow structure over a naca-0015 airfoil,” Experiments in fluids, vol. 54, pp. 1-14, 2013.

[0035]

[0015] W. Zhang and R. Samtaney, “A direct numerical simulation investigation of the synthetic jet frequency effects on separation control of low-re flow past an airfoil,” Physics of Fluids, vol. 27, no. 5, 2015.

[0036]

[0016] J. -J. Wang, K.-S. Choi, L.-H. Feng, T. N. Jukes, and R. D. Whalley, “Recent developments in dbd plasma flow control,” Progress in Aerospace Sciences, vol. 62, pp. 52-78, 2013.

[0037]

[0017] M. Sato, T. Nonomura, K. Okada, K. Asada, H. Aono, A. Yakeno, Y. Abe, and K. Fujii, “Mechanisms for laminar separated-flow control using dielectric-barrier-discharge plasma actuator at low reynolds number,” Physics of Fluids, vol. 27, no. 11 , 2015.

[0038]

[0018] M. G. De Giorgi, A. Ficarella, F. Marra, and E. Pescini, “Micro dbd plasma actuators for flow separation control on a low pressure turbine at high altitude flight operating conditions of aircraft engines,” Applied Thermal Engineering, vol. 114, pp. 511-522, 2017.

[0039]

[0019] L. N. Cattafesta III and M. Sheplak, “Actuators for active flow control,” Annual Review of Fluid Mechanics, vol. 43, pp. 247-272, 2011.

[0040]

[0020] S. J. Johnson, J. P. Baker, C. Van Dam, and D. Berg, “An overview of active load control techniques for wind turbines with an emphasis on microtabs,” Wind Energy: An International Journal for Progress and Applications in Wind Power Conversion Technology, vol. 13, no. 2-3, pp. 239-253, 2010.

[0041]

[0021] T. K. Barias and G. A. van Kuik, “Review of state of the art in smart rotor control research for wind turbines,” Progress in Aerospace Sciences, vol. 46, no. 1 , pp. 1-27, 2010.

[0042]

[0022] O. Stalnov, A. Kribus, and A. Seifert, “Evaluation of active flow control applied to wind turbine blade section,” Journal of Renewable and Sustainable Energy, vol. 2, no. 6, 2010.

[0043]

[0023] V. Maldonado, M. Boucher, R. Ostman, and M. Amitay, “Active vibration control of a wind turbine blade using synthetic jets.,” International Journal of Flow Control, vol. 1 , no. 4, 2009.

[0044]

[0024] V. Maldonado, J. Farnsworth, W. Gressick, and M. Amitay, “Active control of flow separation and structural vibrations of wind turbine blades,” Wind Energy: An International Journal for Progress and Applications in Wind Power Conversion Technology, vol. 13, no. 2-3, pp. 221-237, 2010.

[0045]

[0025] K. Taylor, C. M. Leong, and M. Amitay, “Load control on a dynamically pitching finite span wind turbine blade using synthetic jets,” Wind Energy, vol. 18, no. 10, pp. 1759-1775, 2015.

[0026] S. A. Tran, D. A. Corson, and O. Sahni, “Synthetic jet based active flow control of dynamic stall phenomenon on wind turbines under yaw misalignment,” in 32nd ASME wind energy symposium, p. 0871 , 2014.

[0046]

[0027] T. T. Rice, K. Taylor, and M. Amitay, “Wind tunnel quantification of dynamic stall on an s817 airfoil and its control using synthetic jet actuators,” Wind Energy, vol. 22, no. 1 , pp. 21-33, 2019.

[0047]

[0028] M. Moshfeghi and N. Hur, “Numerical study on the effects of a synthetic jet actuator on s809 airfoil aerodynamics at different flow regimes and jet flow angles,” Journal of Mechanical Science and Technology, vol. 31 , pp. 1233-1240, 2017.

[0048]

[0029] V. Maldonado and S. Gupta, “Increasing the power efficiency of rotors at transitional reynolds numbers with synthetic jet actuators,” Experimental Thermal and Fluid Science, vol. 105, pp. 356-366, 2019.

[0049]

[0030] V. Troshin and A. Seifert, “Performance recovery of a thick turbulent airfoil using a distributed closed-loop flow control system,” Experiments in fluids, vol. 54, pp. 1-19, 2013.

[0050]

[0031] J. Wu, M. Shen, and L. Jiang, “Role of synthetic jet control in energy harvesting capability of a semi-active flapping airfoil,” Energy, vol. 208, p. 118389, 2020.

[0051]

[0032] P. Wang, Q. Liu, C. Li, W. Miao, M. Yue, and Z. Xu, “Investigation of the aerodynamic characteristics of horizontal axis wind turbine using an active flow control method via boundary layer suction,” Renewable Energy, vol. 198, pp. 1032-1048, 2022.

[0052]

[0033] Y. Sun, Y. Qian, Y. Gao, T. Wang, and L. Wang, “Stall control on the wind turbine airfoil via the single and dual-channel of combining bowing and suction technique,” Energy, vol. 290, p. 130224, 2024.

[0053]

[0034] A. Tummala, R. K. Velamati, D. K. Sinha, V. Indraja, and V. H. Krishna, “A review on small scale wind turbines,” Renewable and Sustainable Energy Reviews, vol. 56, pp. 1351-1371 , 2016.

[0054]

[0035] F. Mlihle, J. Schottler, J. Bartl, R. Futrzynski, S. Evans, L. Bernini, P. Schito, M. Draper, A. Guggeri, E. Kleusberg, et al., “Blind test comparison on the wake behind a yawed wind turbine,” Wind Energy Science, vol. 3, no. 2, pp. 883-903, 2018.

[0055]

[0036] M. Selig, “Low reynolds number airfoil design lecture notes,” VKI Lecture Series, November, pp. 24-28, 2003.

[0056]

[0037] A. Ducoin, J.-C. Loiseau, and J.-C. Robinet, “Numerical investigation of the interaction between laminar to turbulent transition and the wake of an airfoil,” European Journal of Mechanics-B / Fluids, vol. 57, pp. 231-248, 2016.

[0057]

[0038] A. Nabhani, N. M. Tousi, M. Coma, G. Bugeda, JM. Bergada. “Large-scale horizontal axis wind turbine aerodynamic efficiency optimization using active flow control and synthetic jets”. Energy 319 (2025) 134940.

[0058]

[0039] Farhadi Asghar et al: “Aerodynamic Multi-Parameter Optimization of NACA0012 Airfoil Using Suction / Blowing Jet Technique”. 2016.

[0059]

[0040] Ali Boudis et al: “Aerodynamic optimization of active flow control over S809 airfoil using synthetic jet”. 2019.

[0041] Duvigneau R et al: “Optimal location of a synthetic jet on an Airfoil for Stall control”. 2007.

[0060] DESCRIPTION OF THE INVENTION

[0061] To solve the drawbacks of the prior art, the present invention proposes distinct strategies, which can be implemented independently or in combination, for reducing dynamic forces and vortical structures acting on a wind turbine.

[0062] To that end, a first aspect of the invention provides a method that can be applied to HAWTs and VAWTs for reducing the dynamic forces and vortical structures generated by a tip of a wind turbine blade. In such a case, the method particularly comprises providing a suction pump at a given location of a wind turbine; generating a constant or periodic suction of the fluid at the blade tip by the suction pump; and decreasing a pressure at the blade tip until reducing the vortices up to a given threshold.

[0063] The corresponding pressure required to minimize / suppress the blade tip vortical structures has to be optimized for each wind turbine, blade shape and operating conditions.

[0064] A second aspect of the invention provides a (computer implemented) method for maximizing the power generated by a wind turbine, particularly, a HAWT. The method comprises performing by one or more processors the following steps: dividing a length of each turbine blade into different airfoils; obtaining a boundary layer separation point, a vortex shedding frequency and an amplitude of the dynamic lift and drag forces for each airfoil by analyzing the flow around each airfoil considering a wind speed, a turning speed and an angle of attack; implementing Active Flow Control (AFC) technology on each airfoil; performing an AFC rotor performance optimization process for each airfoil, wind speed and angle of attack through a parametric optimization or using an optimizer implementing an optimization algorithm, using the obtained boundary layer separation point and vortex shedding frequency, the AFC rotor performance optimization process comprising obtaining different AFC actuation parameters including an AFC actuator location, an actuator width, a momentum coefficient, an inclination angle and a pulsating flow frequency; and integrating the different parameters along the turbine blade span.

[0065] Therefore, the optimization of the AFC actuation parameters is performed across all blade sections along the blade span and fora range of wind speeds (Reynolds numbers) and angles of attack. These parameters can be estimated based on the wind turbine’s operational envelope - i.e. the range of conditions within which the turbine can function safely and efficiently, for instance, the wind speed, temperature, power output, and mechanical limits. That is, unlike previous approaches that optimize a single airfoil and inflow angle, the proposed method anticipates the optimal actuation for every load case, ensuring enhanced performance across varying operational conditions.

[0066] According to the invention, the AFC technology can comprise Synthetic Jets, Fluidic Oscillators, Plasma Actuators, constant / periodic blowing and / or suctions devices.

[0067] In some embodiments, the AFC rotor performance optimization process comprises locating the actuator in the vicinity of the boundary layer separation point, and aligning the pulsating flow frequency with the vortex shedding frequency.

[0068] In some embodiments, the method also comprises establishing the actuator width to be between 0.05% - 2% of an airfoil chord length and establishing the inclination angle to range between 0 - 180 degrees, measured relative to a tangent of the airfoil section at the actuator location.

[0069] In case the parametric optimization is used, the method can comprise keeping two of the AFC actuation parameters constant and optimizing the other three AFC actuation parameters. For instance, the actuator location and actuator width can be maintained constant and the other three parameters can be optimized until obtaining an optimal momentum coefficient, optimal inclination angle and optimal pulsating flow frequency that maximize or minimize a desired integrated parameter of the airfoil.

[0070] In case an optimizer is used, the method can comprise defining maximum and lower values for each AFC actuation parameter, and randomly assigning values to each AFC actuation parameter until obtaining an optimal values combination of the AFC actuation parameters that maximizes or minimizes a desired integrated parameter of the airfoil section.

[0071] In some embodiments, the desired integrated parameter comprises an airfoil efficiency, an airfoil lift or an airfoil torque.

[0072] The optimization algorithm can be based in any optimization system / algorithm. In some embodiments, a genetic algorithm is used. In this case, the genetic algorithm provides a Pareto Front (i.e. the group of final results that maximizes or minimizes what is desired, for example, the wing efficiency, the lift, etc.

[0073] Therefore, the second solution provided by the present invention optimizes the power generated by any HAWT, at any wind speed. At very large wind speeds it can also decrease the time averaged torque generated by the blades, because it promotes the boundary layer separation, and therefore allows the wind turbine to work at higher wind speeds than the initially designed one. At very small wind speeds it allows the turbine to operate at smaller wind speeds than the minimum designed one, because it delays the boundary layer separation. All this is done via delaying the boundary layer separation along the entire blade (at small and medium wind speeds), or promoting the boundary layer separation (at large wind speeds). At small and medium wind speeds, and due to the delay of the boundary layer separation, the downstream vortical structures are drastically reduced; therefore allowing to considerably decrease the distance between two consecutive wind turbines.

[0074] Furthermore, the present invention can be directly employed for minimizing the downstream vortical structures of the wind turbines, allowing the wind turbines to be placed closer to each other, without having to consider the effect that the boundary layer separation has on the wind turbine's power output.

[0075] Moreover, the vortical structures generated by the tip of the blades at any wind speed, can be minimized or even suppressed when the AFC technology is applied in the form of constant or periodic suction.

[0076] Other embodiments of the invention that are disclosed herein also include software programs to perform the method embodiment steps and operations summarized above and disclosed in detail below. More particularly, a computer program product is one embodiment that has a computer-readable medium including computer program instructions encoded thereon that when executed on at least one processor in a computer system causes the processor to perform the operations indicated herein as embodiments of the invention.

[0077] BRIEF DESCRIPTION OF THE DRAWINGS

[0078] The previous and other advantages and features will be more fully understood from the following detailed description of embodiments, with reference to the attached figures, which must be considered in an illustrative and non-limiting manner, in which:

[0079] Fig. 1 is a flow chart illustrating a method for maximizing the power generated by a wind turbine, particularly a HAWT, according to an embodiment of the present invention.

[0080] Fig. 2 illustrates an example of blade airfoils (sections).

[0081] DETAILED DESCRIPTION OF THE INVENTION AND OF PREFERRED EMBODIMENTS Present invention provides different methodologies for wind turbines, which can be used in combination or independently.

[0082] Specifically, according to a first solution, the present invention provides a method for reducing or minimizing the vortices generated by the blade tip. This is applicable to all HAWT, as well as to all VAWT, regardless of the operating wind speed. Present invention’s approach to minimize or completely eliminate the blade tip vortical structures is based on generating a constant or periodic suction of the fluid at this particular location. A suction pump is placed somewhere in the wind turbine generating a constant or periodic suction of the fluid at the tip. Then, the pressure at the tip is reduced until reducing the vortices up to a given threshold.

[0083] Additionally, according to a second solution, the present invention also provides a method to maximize the power generated by any HAWT. The procedure is performed for each turbine blade section and is applicable to any operating wind speed.

[0084] Fig. 1 shows an embodiment of such method for maximizing the power. According to this embodiment, at step 101 , the entire blade length is divided in a sufficient set of different airfoils, as shown in Fig. 2. The higher the number of sections the higher the accuracy of the procedure.

[0085] At step 102, for each airfoil the following parameters are determined a) the boundary layer separation point, b) the vortex shedding frequency and c) the amplitude of the dynamic lift and drag forces by analyzing the flow around each airfoil taking into account the corresponding wind speed, turning speed and angle of attack (or inflow angle).

[0086] At step 103, AFC technology is applied to each airfoil. Possible AFC technologies include Synthetic Jets (SJ), Fluidic Oscillators (FO), Plasma Actuators, also known as Dielectric Barrier Discharge (DBD) actuators, constant or periodic blowing and / or sucking.

[0087] Then, at step 104, an AFC rotor performance optimization process is carried out, whether parametric or using an optimizer. The AFC rotor performance optimization process, for each particular airfoil, wind speed and angle of attack, uses the previously obtained boundary layer separation point and vortex shedding frequency and, regardless of the AFC technology employed, it focuses on five AFC actuation parameters (actuator location, actuator width, momentum coefficient, inclination angle and pulsating flow frequency). The output of this step is the optimum AFC actuation parameters associated with each blade section and wind speed. Finally, at step 105, the different parameters are integrated along the turbine blade span. This allows a smooth evolution of each AFC actuation parameter along the blade span in order to make sure that the boundary layer remains attached (or separately, if required) in all airfoils (blade sections).

[0088] Consequently, the method designs and implements a feedback control loop, or closed loop control system, that can read a fluid parameter (for example the pressure), on the trailing edge of each airfoil, and can modify one or several parameters associated to the implemented AFC technology to obtain the desired effect on the boundary layer along the entire blade span.

[0089] In an embodiment, the AFC rotor performance optimization process considers locating the AFC actuator around the boundary layer separation point, ideally, locating it precisely at this point. The pulsating flow frequency is aligned with the vortex shedding frequency.

[0090] For the rest of the AFC actuation parameters, the actuator width is established to be between 0.05% - 2% of the airfoil chord length based on the inventor’s and other researchers experience. On the other hand, the inclination angle is established to range between 0 - 180 degrees, measured relative to a tangent of the airfoil section at the actuator location

[0091] The momentum coefficients required greatly differ from different applications, usual values fall between 0.001 and 0.2, the tendency being of using greater momentum coefficients as the Reynolds number increases. At this point, is relevant to remember that the Reynolds number involved in the present large wind turbine applications are of the order of 12000000, and there are no studies of AFC applications at these very high Reynolds numbers, except for

[0038] , Therefore to start the AFC rotor performance optimization process a possible momentum coefficient is estimated and the entire range of momentum coefficients studied to obtain the optimum one for each application.

[0092] Parametric optimization (small and medium wind speeds)

[0093] In an embodiment, the location of the actuator (corresponding to the previously determined boundary layer separation point) and the actuator width are kept constant throughout the process. This approach is not restrictive; other embodiments may keep different AFC actuation parameters constant while optimizing the remaining ones.

[0094] The process begins by establishing a constant jet inclination angle of approximately 10 degrees and a constant pulsation flow frequency matching the previously determined vortex shedding frequency. With these four parameters (actuator location, actuator width, jet inclination angle, and pulsation flow frequency) fixed, a range of momentum coefficients is evaluated to identify the one that maximizes the airfoil efficiency.

[0095] Once the optimal momentum coefficient is determined, it is kept constant along with the actuator location, actuator width, and pulsation flow frequency. Next, the inclination angle is varied until the angle that maximizes airfoil efficiency is found. At this stage, with both the optimum momentum coefficient and inclination angle established, the method evaluates different pulsating flow frequencies to find the one that maximizes airfoil efficiency. The parametric optimization process can stop here, having optimized the momentum coefficient, inclination angle, and pulsating flow frequency.

[0096] It is important to note that the method could also initiate the parametric optimization process by evaluating the inclination angle first, followed by the momentum coefficient, and so on. Furthermore, the actuator location and width, which were initially kept constant, could also be optimized. In other words, the parametric optimization process can begin with any parameter and can be extended to optimize all five AFC actuation parameters instead of fixing two and optimizing the other three.

[0097] Using an optimizer (small and medium wind speeds)

[0098] The advantage of using an optimizer resides in the fact that the entire optimization process is automatic. This means that the five AFC actuation parameters can be automatically modified after each set of simulations (each generation). To begin, the limits of each parameter are defined, specifying the maximum and minimum values to be explored for each of the five AFC actuation parameters. For example, the actuator location can move upstream and downstream, 10% or 20% of the airfoil chord length from the estimated initial position, which is the boundary layer separation point. The range of the pulsating frequency to be evaluated falls between 0.1 and ten times the vortex shedding frequency. The boundaries of the inclination angle can be placed between one and 25 degrees versus the airfoil tangent position, the actuator width estimated range can be between 0.05% and 0.3% of the airfoil chord length and finally the momentum coefficients evaluated can be between 0.001 and 0.5 for example.

[0099] Once the limits of each AFC actuation parameter in the optimization process have been defined, the optimizer keeps changing the different parameters until obtaining the optimum values for each of the five AFC actuation parameters, this process is long, it has to be applied for each airfoil and for each wind turbine operating condition, and requires a large number of computational fluid mechanics simulations.

[0100] In some embodiments, the optimization algorithm executed by the optimizer is a genetic algorithm because it’s good at finding the best solutions when multiple goals are had. Moreover, it’s also strong enough to handle situations where some configurations might not work due to technical issues. Nevertheless, the use of the genetic algorithm is not limitative as any optimization algorithm could be equally used.

[0101] Once the optimization process is finished (whether parametric or using and optimizer), the resulting lift and drag dynamic peak to peak amplitudes can be compared with the ones obtained in step 102 (baseline-case). The peak to peak amplitudes of the forces obtained after the optimization process, should be much smaller than the ones obtained in the baseline-case.

[0102] Case optimization for the maximum wind speed

[0103] The process is similar to the one previously described for small and medium wind speeds; however, in this case, the aim is not to maximize airfoil efficiency, lift, or torque. Instead, the objective is to reduce these factors to prevent damage to the rotor blades. That is, for high wind speeds, the boundary layer is separated selectively whereas for medium and low speeds it is attached. Both parametric optimization and optimization using an optimizer can be used for this purpose.

[0104] The present invention has been described in particular detail with respect to specific possible embodiments. Those of skill in the art will appreciate that the invention may be practiced in other embodiments. For example, the nomenclature used for components, capitalization of component designations and terms, the attributes, data structures, or any other programming or structural aspect is not significant, mandatory, or limiting, and the mechanisms that implement the invention or its features can have various different names, formats, and / or protocols. Further, the system and / or functionality of the invention may be implemented via various combinations of software and hardware, as described, or entirely in software elements. Also, particular divisions of functionality between the various components described herein are merely exemplary, and not mandatory or significant. Consequently, functions performed by a single component may, in other embodiments, be performed by multiple components, and functions performed by multiple components may, in other embodiments, be performed by a single component. Certain aspects of the present invention include process steps or operations and instructions described herein in an algorithmic and / or algorithmic-like form. It should be noted that the process steps and / or operations and instructions of the present invention can be embodied in software, firmware, and / or hardware, and when embodied in software, can be downloaded to reside on and be operated from different platforms used by real-time network operating systems.

[0105] The scope of the present invention is defined in the following set of claims.

Claims

CLAIMS1 . A method for reducing the dynamic forces and vortical structures acting on a wind turbine, the method comprising: providing a suction pump at a given location of a wind turbine; generating a constant or periodic suction of the fluid at a blade tip by the suction pump; and decreasing a pressure at the blade tip until reducing the vortices up to a given threshold.

2. The method of claim 1 , further comprising maximizing a power generated by a wind turbine by means of a processor performing the following steps: dividing a length of each turbine blade into different airfoils; obtaining a boundary layer separation point, a vortex shedding frequency and an amplitude of the dynamic lift and drag forces for each airfoil by analyzing the flow around each airfoil considering a wind speed, a turning speed and an angle of attack; implementing Active Flow Control, AFC, technology on each airfoil; performing an AFC rotor performance optimization process for each airfoil, wind speed and angle of attack through a parametric optimization or using an optimizer implementing an optimization algorithm, using the obtained boundary layer separation point and vortex shedding frequency, the AFC rotor performance optimization process comprising obtaining different AFC actuation parameters including an actuator location, an actuator width, a momentum coefficient, an inclination angle and a pulsating flow frequency; and integrating the different parameters along the wind turbine blade span.

3. The method of claim 2, wherein the AFC rotor performance optimization process comprises: locating the actuator in a vicinity of the boundary layer separation point; and aligning the pulsating flow frequency with the vortex shedding frequency.

4. The method of claim 3, further comprising establishing the actuator width to be between 0.05% - 2% of an airfoil chord length and establishing the inclination angle to range between 0 - 180 degrees, measured relative to a tangent of the airfoil section at the actuator location.

5. The method of claim 4, wherein the parametric optimization is used, the method comprising keeping two of the AFC actuation parameters constant and optimizing the other three AFC actuation parameters.

6. The method of claim 5, wherein the actuator location and actuator width are maintained constant and the other three parameters are optimized until obtaining an optimal momentum coefficient, optimal inclination angle and optimal pulsating flow frequency that maximize or minimize a desired integrated parameter of airfoil section.

7. The method of claim 4, wherein an optimizer is used, the method comprising: defining maximum and lower values for each AFC actuation parameter; and randomly assigning values to each AFC actuation parameter until obtaining an optimal values combination of the AFC actuation parameters that maximizes or minimizes a desired integrated parameter of the airfoil considered.

8. The method of claim 7, wherein the desired integrated parameter comprises an airfoil efficiency or an airfoil lift.

9. The method of any one of the previous claims, wherein the AFC technology comprises Synthetic Jets, Fluidic Oscillators, Plasma Actuators, a constant or periodic blowing, or a suction device.

10. The method of claim 1 , wherein the wind turbine comprises a Horizontal Axis Wind Turbine, HAWT, or a Vertical Axis Wind Turbine, VAWT.

11. A non-transitory computer readable medium comprising executable instructions that, when executed by one or more processors of a computer system cause the computer system to perform the method of any one of the preceding claims 2 to 9.

Images