Wind Turbine Aeromechanics Research at UMD

Transcription

1 Wind Turbine Aeromechanics Research at UMD Anholt Wind Farm PI: James D. Baeder Alfred Gessow Rotorcraft Center of Excellence Department of Aerospace Engineering Wind Energy Thrust Leader, UMERC University of Maryland, College Park NREL/NTWC June 30, 2015

2 Offshore Wind Energy is Inherently Multi-Disciplinary Assembled a team across five departments at three universities to support eight faculty members and eight students in research! Leverages existing expertise at UMERC, ESSIC, CALCE and AGRC! Significant cost share from UMCP campus! Research focused on two major thrust areas! Builds a foundation for long-term OSW research in Maryland Research Thrusts Co-PIs GRA/URA Refined Wind Resource Characterization For Next-Gen Prognostics & Health Management 1.1 Offshore Wind Power Resource Evaluation and Wind Forecasting 1.2 Development of a Maintenance Option Model to Enable Optimized OSW Farm Sustainment 1.3 Frostburg Testbed and Data Analysis 1.4 Prognostics and Health Management of Offshore Wind Turbine Innovative Wind Turbine Aeromechanics to Aid Energy Capture 2.1 Wake Alleviation Devices for OSW Turbines 2.2 CFD Modeling of OSW Turbines and Wake Interactions 2.3 3D Rotor Design for Maximum Mechanical Performance and Safety 2.4 Advanced Composite Couplings for Passive Fatigue Loads Reduction in OSW Turbines List of Research Tasks and Co-PIs and Student Support Zeng Sandborn Eltayeb Azarian Jones Baeder Goloubev Chopra GRA GRA URA GRA GRA GRA URA GRA

3 Thrust 2: Innovative Wind Turbine Aeromechanics to Aid Energy Capture Thrust 2 Objectives: This thrust addresses key offshore wind turbine aeromechanic issues through a balanced experimental / computational / analytical program to: (1) evaluate the effect of blade-mounted devices on the structure, strength and dissipation of the turbine wake at laboratory scales; and (2) validate CFD modeling of laminar/turbulent transition and other turbulence modeling issues for investigating 2-D airfoil characteristics as well as 3-D slotted tips and tubercle shaped leading edges; and (3) develop state-of-the-art comprehensive aeromechanics analysis for investigating the effects of tailored composite couplings on vibratory loads. This thrust is a collaboration between AGRC (UMCP Department of Aerospace Engineering) and Bowie State University.

4 Task 2.1 Wake Alleviation Devices for OSW Turbines (UMCP/AGRC: Anya Jones, Vera Klimchenko) Objec+ve(( Evaluate(effec+veness(of(blade(mounted( devices( on( the( dissipa+on( of( wind( turbine(wake( U"" " Methods( Sub=scale(wind(tunnel(tes+ng( PIV(used(to(iden+fy(and(analyze(wake(structures(behind(wind( turbine( Evaluate(effect(of(blade(mounted(devices(on(energy(produc+on( Baseline( Winglet( Serrated(

5 Small Scale Wind Turbine Design Linear(Taper( Root"Chord"=(1.173(in( Diameter(=(8.4(in 2( Low(Reynolds(airfoil(SG6040" Tip"Chord"=((1/4)(Root(Chord( Tower"Height"=(10(in( Op+mum(Rotor(Theory(used(to(calculate(twist(necessary( to(achieve(design(tsr( Twist"at"the"root"is(30.2(degrees( Twist"at"the";p"is(=3.7(degrees( Design(TSR""of"6(

6 Wind Tunnel and PIV Setup FOV(1( FOV(4( FOV(2( FOV(3( U"" " Four(Fields(of(View((FOV)(for(PIV( S;tched"together"

7 Time=Averaged(Total(Velocity(for(3(Cases( Baseline(Case(Average(Velocity( U"" " U"" " Velocity(Cut(at(x/D(=(=1.8(( U"" "

8 Snapshot(of(Total(Velocity((and(Vor+city(for(3(Cases( "" " "" " "" "

9 Current Work 1. Characterize(the(strength(of(the(+p(vor+ces( 2. Validate(the(theory(that(blade(mounted(devices(can(weaken(+p(vor+ces( 3. Analyze(wake(recovery(by(comparing(the(+me=averaged(data(for(the(three(cases(

10 Task 2.2 CFD Modeling of OSW Turbines and Wake Interactions (UMCP/AGRC: James Baeder, Taylor Rinehart) Introduction - Motivation Sandia 100-meter all-glass HAWT blade design (13.2 MW) for offshore wind turbine applications! High-resolution CFD CSD simulations expensive for design and analysis of multiple new configurations! Accurate aeroelastic analysis tools needed for quick estimates of aerodynamic/structural loads, and rotor thrust/torque Aeroelastic analysis tools depend on reduced-order aerodynamic modeling (2-D airfoil characteristics, stall models)! Lift / Drag/ Pitching moment queried from look-up tables! Experiments costly to generate airfoil tables for new airfoil designs over wide range of angles of attack & Reynolds numbers! CFD simulations provide an affordable alternative for generating look-up tables if accurate! Use of GPGPU can provide 40-50x speedup over single core on CPU

11 CFD Modeling of OSW Turbines and Wake Interactions Introduction Turbulence Model Limitations S809 Airfoil, Re = 10 6 Near Copenhagen (Wikipedia) Accuracy of design estimates depends on accuracy of airfoil characteristics Issues with conventional RANS turbulence models in CFD simulations Over-predict drag for partially laminar boundary layers at low AoA Over-predict maximum lift and stall onset angle (incipient separation) Need accurate laminar/turbulent transition; adverse pressure gradient

12 CFD Modeling of OSW Turbines and Wake Interactions 2-D Airfoil Results on S809 Airfoil at Re = 1,000,000 Added Laminar/Turbulent Transition & Adverse Pressure Gradient Corrections SA-Transition-APG model combines best of both Predicts the lower drag at lower AoA and earlier stall onset at higher AoA Looking at incorporating surface roughness effects; cross-flow; DDES

13 CFD Modeling of OSW Turbines and Wake Interactions Sandia 100m w/ and w/o tubercles 3-D URANS simulations 10 m/s Excellent comparison with previous analysis (without tubercles) No detrimental effect on performance to include tubercles

14 Subdivision and Hamiltonian Paths Triangular unstructured grid Divide triangle into three quadrilaterals Loops constructed by connecting the midpoint of edges Loops formed through all triangles connected by a triangular node Each face part of only one distinct loop Each cell centroid is intersected by loops of different colors Makes the algorithms very efficient!!!

15 CFD Modeling of OSW Turbines and Wake Interactions Plans Continue to validate CFD for Sandia 100-meter all-glass HAWT blade design (13.2 MW) for offshore wind turbine applications! Continue improving physical modeling in CFD; CFD/CSD coupling Continue to investigate tubercle shaped leading-edge, but for desensitizing turbines to upstream disturbances! Complimentary to experiments Investigate slotted tip for reducing tip vortex swirl using CFD tools! Momentum and turbulence increases vortex diffusion CFD/CSD to aid in studying tailored composite coupling Developing Hamiltonian/Strand solver HAMSTR for DOD

16 Task 2.4 Advanced Composite Couplings for Passive Fatigue Load Reduction in OSW (UMCP/AGRC: Chopra; Ananthan) Composite tailoring technology: Intentional distribution of fiber orientation and layup Meet specific structural requirements Achieve desired elastic couplings Lag Bending Nose Down Twist Successful applications in aerospace: Forward sweep wing of Grumman X-29 Composite tailored couplings on the AgustaWestland AW101 rotor blade for vibratory loads reduction Potential benefits of composite tailored couplings on wind turbines: Blade loads reduction Rotor/tower aeromechanical stability improvement Tower vibration reduction

17 Comprehensive Analysis AGRC Modeling Features " Multiple configurations (SMR/TR, Coax, Tilt-Rotor, Wind Turbine) " Flexible blades with flap, lag, torsion and axial DOF ( Euler-Bernoulli ) " Multibody-type kinematics and rotating reference frames " Coupled with free wake, panel* method, CFD " Modular inclusion of phenomena into dynamics and trim " All geometric nonlinearities large deflections " Modal reduction or full nodal DOF analysis " Composite coupling, gravity loads and material nonlinearities " FET, harmonic balance, time integration for trim " Free-flight, wind-tunnel mode for single/multi-rotor configurations

18 Power Regulation with RPM, Pitch 50% RPM 100% RPM 75% RPM 2.5 o pitch 0 o pitch Power limit to reduce loads 5 o pitch RPM reduction used to extract more power at low wind speeds Blade pitch used to limit power at high speeds

19 Results Tip Deflection w/r v/r 50% stiffness Elastic blade With atm. Boundary layer Rigid blade Rigid blade 50% stiffness Tip deflection scales with inverse of bending stiffness Steady out-of-plane deflection = 5 m Good agreement with AcuSolve predictions (Corson et. al, AIAA 2012)

20 Results Flap Bending Moment Less Ft-lb 50% stiffness Elastic blade With atm. boundary layer Rigid blade More Blade elasticity alleviates root flap bending moment Minor effect from atmospheric boundary layer

21 Results Torsion Moment Ft-lb 50% stiffness Elastic blade With atm. boundary layer Rigid blade Elastic motions introduce torsion loads into the blade through flap bending, chordwise force and dynamic twist

22 Summary " Applied multi-body UMARC version to turbines " Implemented " Axial elasticity in blade dynamics " Coupling with Maryland Free Vortex Wake for turbines " Wind velocity gradient in wake " Gravity loading on blades " Performance dominated by external geometry (aerodynamics) " Flap bending loads, blade stiffness and atmospheric boundary layer (oscillatory forcing) drive oscillations in hub moments " Implemented basic composite coupling for turbine blades " Modeling of axial-flap-lag-torsion structural interactions! Rotor- tower couplings (in progress)! Targets for Year-2 Characterizing cross-section properties # Fiber orientation $ blade properties (and cross-couplings) # Validation of composite beam models Blade load reduction with bend-twist, extension twist coupling

23 Summary Assembled a team across five departments at three universities to support eight faculty members and eight students in research! Leverages existing expertise at UMERC, ESSIC, CALCE and AGRC! Significant cost share from UMCP campus! Research focused on two major thrust areas! Builds a foundation for long-term OSW research in Maryland Developing facilities and analytical/computational tools! Many now in place generating preliminary results! Frostburg and Bowie became more involved this past summer with URA Already presenting results at conferences! Posters at Offshore EWEA 2013/15; Three papers at NAWEA 2015! Papers at AIAA Scitech 32 nd ASME Wind Energy Symposium January 2014/15(2-D CFD;Sustainment;Flatback); Several other AIAA; PHM

24 Synergistic Activities at UMCP in Wind Energy Wind Energy Thrust as part of UMERC Three faculty members visited Europe! Attended: Offshore EWEA and met with vendors; O&M Conference! Made contacts with researchers at DTU and OSW industry! Visited Anholt Offshore Wind Farm Wind Energy Theory Class (ENAE788I) taught last fall and in 2014! 15 grads and 2 undergrads (3 and 2 this fall); Guest lectures by specialists! Final exam/project looked at NREL 5MW Offshore Design (2009 Jonkman) at Buoy St The Maryland Wind TERPines formed last fall! Student organization looking at promoting wind energy on campus! Developing facilities / analysis for participating in DOE Wind Energy Challenge Offshore Wind Energy Seminar Series started last fall! Continued last spring; brought in faculty from Delaware and Maine this fall Eight UMCP faculty investigating Offshore Wind Energy Issues! Hope to expand; make long-term; and work with Maryland companies (BizMDOSW)

25 University of Maryland College Park Wind Turbine Aeromechanics Potential Collaboration GPGPU for Aeromechanic AnalysisResults and Discussion: Lift vs.! Laminar/Turbulent transition with APG correction and DDES! Coupling isolated blade with grid motion to FV-W, CSD and multiple GPGPU cards Hamiltonian/Strand Solver Halfway cut 3/4 flatback! Unstructured triangles turned into quads that form Hamiltonian chains! Strands in the third direction! Can then examine complicated nacelle / tower interactions Halfway cut 1/2 flatback Large Blades!!!! (a) Surface unstructured mesh (b) Hamiltonian loops on surface Aeroelastic instabilities at extreme scales? Individual blade control (using LIDAR for feedback) Tubercles, slotted tips to desensitize turbines 90C cut 3/4 flatback Flatback airfoils with wavy trailing edges The Maryland Wind TERPines! Examining cyclo-turbine VAWT! Examining shroud for DAWT (c) Volume mesh (d) Curved strand grids Fig. 23. Longitudinal slice of the mesh system around the Robin fuselage highlighting the curved strands in concave fuselage surface. Quad-Rotor Bi-Plane Tailsitter Stagnation point Stagnation point Flow deceleration Flow Zero velocity on fuselage surface