0
Research Papers

Enhancement of Free Vortex Filament Method for Aerodynamic Loads on Rotor Blades

[+] Author and Article Information
Hamidreza Abedi

Division of Fluid Dynamics,
Department of Applied Mechanics,
Chalmers University of Technology,
Göteborg SE-412 96, Sweden
e-mail: abedih@chalmers.se

Lars Davidson

Division of Fluid Dynamics,
Department of Applied Mechanics,
Chalmers University of Technology,
Göteborg SE-412 96, Sweden
e-mail: lada@chalmers.se

Spyros Voutsinas

Fluid Section,
School of Mechanical Engineering,
National Technical University of Athens,
Athens 15780, Greece
e-mail: spyros@fluid.mech.ntua.gr

Contributed by the Solar Energy Division of ASME for publication in the JOURNAL OF SOLAR ENERGY ENGINEERING: INCLUDING WIND ENERGY AND BUILDING ENERGY CONSERVATION. Manuscript received February 12, 2016; final manuscript received January 26, 2017; published online March 16, 2017. Assoc. Editor: Douglas Cairns.

J. Sol. Energy Eng 139(3), 031007 (Mar 16, 2017) (12 pages) Paper No: SOL-16-1072; doi: 10.1115/1.4035887 History: Received February 12, 2016; Revised January 26, 2017

The aerodynamics of a wind turbine is governed by the flow around the rotor, where the prediction of air loads on rotor blades in different operational conditions and its relation to rotor structural dynamics is one of the most important challenges in wind turbine rotor blade design. Because of the unsteady flow field around wind turbine blades, prediction of aerodynamic loads with high level of accuracy is difficult and increases the uncertainty of load calculations. An in-house vortex lattice free wake (VLFW) code, based on the inviscid, incompressible, and irrotational flow (potential flow), was developed to study the aerodynamic loads. Since it is based on the potential flow, it cannot be used to predict viscous phenomena such as drag and boundary layer separation. Therefore, it must be coupled to tabulated airfoil data to take the viscosity effects into account. Additionally, a dynamic approach must be introduced to modify the aerodynamic coefficients for unsteady operating conditions. This approach, which is called dynamic stall, adjusts the lift, the drag, and the moment coefficients for each blade element on the basis of the two-dimensional (2D) static airfoil data together with the correction for separated flow. Two different turbines, NREL and MEXICO, are used in the simulations. Predicted normal and tangential forces using the VLFW method are compared with the blade element momentum (BEM) method, the GENUVP code, and the MEXICO wind tunnel measurements. The results show that coupling to the 2D static airfoil data improves the load and power predictions while employing the dynamic stall model to take the time-varying operating conditions into consideration is crucial.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.

References

Hansen, M. O. , 2008, Aerodynamics of Wind Turbines, 2nd ed., EarthScan, London.
van Garrel, A. , 2003, “ Development of a Wind Turbine Aerodynamics Simulation Module,” Energy Research Centre of the Netherlands (ECN), Petten, The Netherlands, Report No. ECN-C–03-079.
Vermeer, L. , Sørensen, J. , and Crespo, A. , 2003, “ Wind Turbine Wake Aerodynamics,” Prog. Aerosp. Sci., 39(6–7), pp. 467–510. [CrossRef]
Opoku, D. G. , Triantos, D. G. , Nitzsche, F. , and Voutsinas, S. G. , 2002, “ Rotorcraft Aerodynamic and Aeroacoustic Modeling Using Vortex Particle Methods,” 23rd International Congress of the Aeronautical Sciences (ICAS), Toronto, ON, Canada, Sept. 8–13, Paper No. ICAS 2002-2.1.3.
Voutsinas, S. G. , Beleiss, M. A. , and Rados, K. G. , 1995, “ Investigation of the Yawed Operation of Wind Turbines by Means of a Vortex Particle Method,” AGARD Conference Proceedings, Vol. 552, pp. 11.1–11.
Zhao, J. , and He, C. , 2010, “ A Viscous Vortex Particle Model for Rotor Wake and Interference Analysis,” J. Am. Helicopter Soc., 55(1), p. 12007. [CrossRef]
Egolf, T. A. , 1988, “ Helicopter Free Wake Prediction of Complex Wake Structures Under Blade-Vortex Interaction Operating Conditions,” 44th Annual Forum of the American Helicopter Society, Washington, DC, June 16–18, pp. 819–832.
Rosen, A. , and Graber, A. , 1988, “ Free Wake Model of Hovering Rotors Having Straight or Curved Blades,” J. Am. Helicopter Soc., 33(3), pp. 11–21. [CrossRef]
Bagai, A. , 1995, “ Contribution to the Mathematical Modeling of Rotor Flow-Fields Using a Pseudo-Implicit Free Wake Analysis,” Ph.D. thesis, University of Maryland, College Park, MD.
Gupta, S. , 2006, “ Development of a Time-Accurate Viscous Lagrangian Vortex Wake Model for Wind Turbine Applications,” Ph.D. thesis, University of Maryland, College Park, MD.
Pesmajoglou, S. , and Graham, J. , 2000, “ Prediction of Aerodynamic Forces on Horizontal Axis Wind Turbines in Free Yaw and Turbulence,” J. Wind Eng. Ind. Aerodyn., 86(1), pp. 1–14. [CrossRef]
Voutsinas, S. , 2006, “ Vortex Methods in Aeronautics: How To Make Things Work,” Int. J. Comput. Fluid Dyn., 20(1), pp. 3–18. [CrossRef]
Chattot, J. , 2007, “ Helicoidal Vortex Model For Wind Turbine Aeroelastic Simulation,” Comput. Struct., 85(11–14), pp. 1072–1079. [CrossRef]
Chattot, J. , 2003, “ Optimization of Wind Turbines Using Helicoidal Vortex Model,” ASME J. Sol. Energy Eng., 125(4), pp. 418–424. [CrossRef]
Holierhoek, J. , de Vaal, J. , van Zuijlen, A. , and Bijl, H. , 2013, “ Comparing Different Dynamic Stall Models,” J. Wind Energy, 16(1), pp. 139–158. [CrossRef]
Leishman, J. , 2002, “ Challenges in Modeling the Unsteady Aerodynamics of Wind Turbines,” AIAA Paper No. 2002-0037.
Bierbooms, W. , 1992, “ A Comparison Between Unsteady Aerodynamic Models,” J. Wind Eng. Ind. Aerodyn., 39(1–3), pp. 23–33. [CrossRef]
Schepers, J. , and Boorsma, K. , 2006, “ Description of Experimental Setup MEXICO Measurements,” Energy Research Centre of the Netherlands (ECN), Petten, The Netherlands, Report No. ECN-X–09-0XX.
Anderson, J. , 2001, Fundamentals of Aerodynamics, 3rd ed., McGraw-Hill, New York.
van Garrel, A. , 2001, “ Requirements for a Wind Turbine Aerodynamics Simulation Module,” 1st ed., Energy Research Centre of the Netherlands (ECN), Petten, The Netherlands, Report No. ECN-C–01-099.
Leishman, J. , and Bagai, M. J. , 2002, “ Free Vortex Filament Methods for the Analysis of Helicopter Rotor Wakes,” J. Aircr., 39(5), pp. 759–775. [CrossRef]
Bhagwat, M. , and Leishman, J. , 2002, “ Generalized Viscous Vortex Model for Application to Free-Vortex Wake and Aeroacoustic Calculations,” 58th Annual Forum and Technology Display of the American Helicopter Society International, Montreal, Canada, June 11–13, pp. 2042–2057.
Vasitas, G. , Kozel, V. , and Mih, W. , 1991, “ A Simpler Model for Concentrated Vortices,” Exp. Fluids, 11(1), pp. 73–76. [CrossRef]
Bagai, A. , and Leishman, J. , 1993, “ Flow Visualization of Compressible Vortex Structures Using Density Gradient Techniques,” Exp. Fluids, 15(6), pp. 431–442. [CrossRef]
Katz, J. , and Plotkin, A. , 2001, Low-Speed Aerodynamics, 2nd ed., Cambridge University Press, New York.
Jonkman, J. , Butterfield, S. , Musial, W. , and Scott, G. , 2009, “ Definition of a 5-MW Reference Wind Turbine for Offshore System Development,” National Renewable Energy Laboratory, Golden, CO, Report No. NREL/TP-500-38060.
Abedi, H. , 2013, “ Development of Vortex Filament Method for Aerodynamic Loads on Rotor Blades,” Licentiate thesis, Chalmers University of Technology, Gothenburg, Sweden.
Chasapogiannis, P. , and Voutsinas, S. , 2014, “ Aerodynamic Simulations of the Flow Around a Horizontal Axis Wind Turbine Using the GAST Software,” National Technical University of Athens, Athens, Greece, Ref. No. 68/1110 (in Greek).
Abedi, H. , Davidson, L. , and Voutsinas, S. , 2015, “ Numerical Studies of the Upstream Flow Field Around a Horizontal Axis Wind Turbine,” AIAA Paper No. 2015-0495.
Reddy, T. S. R. , and Kaza, K. R. V. , 1989, “ Analysis of an Unswept Propfan Blade With a Semi Empirical Dynamic Stall Model,” NASA Lewis Research Center, Cleveland, OH, Report No. NASA-TM-4083.

Figures

Grahic Jump Location
Fig. 1

Schematic for the Biot–Savart law

Grahic Jump Location
Fig. 2

Lifting surface and vortex panels construction

Grahic Jump Location
Fig. 3

Schematic of vortex lattice free wake

Grahic Jump Location
Fig. 4

Numbering procedure

Grahic Jump Location
Fig. 5

Schematic of generation and moving of wake panels at each time step

Grahic Jump Location
Fig. 6

Schematic of wake evolution at the first time step

Grahic Jump Location
Fig. 7

Schematic of wake evolution at the second time step

Grahic Jump Location
Fig. 8

Potential load decomposition

Grahic Jump Location
Fig. 9

Viscous load decomposition

Grahic Jump Location
Fig. 10

Radial distribution of blade elements for NREL-5 MW turbine

Grahic Jump Location
Fig. 11

Distribution of angle of attack along the blade for the NREL turbine: (a) case 2, (b) case 6, and (c) case 10; — potential and – – viscous

Grahic Jump Location
Fig. 12

Distribution of tangential force along the blade for the NREL turbine: (a) case 2, (b)case 6, and (c) case 10; — potential and – – viscous

Grahic Jump Location
Fig. 13

Distribution of normal force along the blade for the NREL turbine: (a) case 2, (b) case 6, and (c) case 10; — potential and – – viscous

Grahic Jump Location
Fig. 14

(a) power and (b) thrust curves for the NREL turbine; — VLFW potential, –– VLFW viscous, – – BEM, – ∘ – GENUVP potential, and – × – GENUVP viscous

Grahic Jump Location
Fig. 15

Radial distribution of blade elements for MEXICO turbine

Grahic Jump Location
Fig. 16

Distribution of normal force along the MEXICO turbine's blade, nonyawed flow: (a) case 1, (b) case 2, and (c) case 3, — potential, – – viscous, and ∘ experiment

Grahic Jump Location
Fig. 17

Distribution of tangential force along the MEXICO turbine's blade, nonyawed flow: (a) case 1, (b) case 2 and (c) case 3, — potential, – – viscous, and ∘ experiment

Grahic Jump Location
Fig. 18

Azimuthal variation of (a) tangential and (b) normal forces at 0.25R of radial position, MEXICO turbine, yawed flow, case 1, — potential, – – viscous, – - – dynamic stall, and ∘ experiment

Grahic Jump Location
Fig. 19

Azimuthal variation of (a) tangential and (b) normal forces at 0.35R of radial position, MEXICO turbine, yawed flow, case 1, — potential, – – viscous, – - – dynamic stall, and ∘ experiment

Grahic Jump Location
Fig. 20

Azimuthal variation of (a) tangential and (b) normal forces at 0.60R of radial position, MEXICO turbine, yawed flow, case 1, — potential, – – viscous, – - – dynamic stall, and ∘ experiment

Grahic Jump Location
Fig. 21

Azimuthal variation of (a) tangential and (b) normal forces at 0.82R of radial position, MEXICO turbine, yawed flow, case 1, —: potential, – – viscous, – - – dynamic stall, and ∘ experiment

Grahic Jump Location
Fig. 22

Azimuthal variation of (a) tangential and (b) normal forces at 0.92R of radial position, MEXICO turbine, yawed flow, case 1, — potential, – –: viscous, – - – dynamic stall, and ∘ experiment

Grahic Jump Location
Fig. 23

Hysteresis loop around the stall angle

Grahic Jump Location
Fig. 24

Definition of the lift coefficient parameters in the ONERA model

Grahic Jump Location
Fig. 25

Definition of the drag coefficient parameters in the ONERA model

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In