Firma convenzione. del Duomo di Milano

Transcription

1 VODCA Firma convenzione Politecnico VAWT di Open Milano Data e for Veneranda Code Assessment Fabbrica Phase I Net Meeting #1 del Duomo di Milano Ilmas Bayati and Marco Belloli Aula Magna Rettorato Mercoledì 27 maggio 2015 Politecnico di Milano Department of Mechanical Engineering x = 0:.01:8*pi; s = x.*sin(x); c = x.*cos(x); plot3(s,c,x)

2 Presentation overview Participants PoliMi IPC machine Long Term: Phase I, II and III Short Term: Phase I proposed outline Open discussion

3 VODCA Participants Institution Country Type Contact Person Politecnico di Milano Italy Academia Bayati DELFT University Cranfield University The Netherlands United Kindom Academia Ferreira Academia Collu EOLFI France Industry Paillard Nenuphar France Industry Kluczewska-Bordier DTU Denmark Academia Paulsen Sandia National Laboratories Ecole Centrale de Nantes United States Government Griffith France Industry/Academia Gilloteaux VUB Bruxelles Belgium Academia Runacres UCL Louvain Belgium Academia Chatelain Others?????

4 PoliMi IPC-VAWT Reference machine for VOCA project: VAWT Open data for Code Assessment Airfoil NACA0021 Chord C 0.25 m Height H 2.7 m Radius R 1.0 m Blades N 3 Solidity σ 0.375

5 PoliMi IPC-VAWT Blades NACA 0021, C = 0.25 m Pitch rotation at 22% of the chord

6 PoliMi s Proposal Phase Validation Activity Management issues I II Preliminary Numerical Vs Numerical PoliMi reference VAWT Investigating the differences among different codes Fixed Pitch φ = 0 Fixed Pitch φ 0 Variable Pitch: IPC simple control laws (e.g 1P) consistent with PoliMi VAWT Numerical Vs Experimental Wind tunnel tests: cases choosen among the ones investigated in the Phase I PoliMi will provide the Person Months for Project coordination Data processing Writing papers For Phase I, II and III PoliMi will provide 5 wind tunnel and person days for the tests III Numerical Vs Experimental (Control) Wind tunnel tests: implementation of IPC advanced control laws, based on the results of Phase I and Phase II Required N wind tunnel and person days for the tests. N to be defined Control laws must be public

7 PoliMi s Proposal Phase Validation Activity Time schedule I Preliminary Numerical Vs Numerical PoliMi reference VAWT Investigating the differences among different codes Fixed Pitch φ = 0 1 Year II Fixed Pitch φ 0 Variable Pitch: IPC simple control laws (e.g 1P) consistent with PoliMi VAWT Numerical Vs Experimental Wind tunnel tests: cases choosen among the ones investigated in the Phase I TBD III Numerical Vs Experimental (Control) Wind tunnel tests: implementation of IPC advanced control laws, based on the results of Phase I and Phase II TBD (IEA task?)

8 Participants list for the Phase I Institution Country Type Contact Person Phase I confirmation Politecnico di Milano Italy Academia Bayati ilmasandrea.bayati@polimi.it TU Delft Cranfield University The Netherlands United Kindom Academia Ferreira C.J.SimaoFerreira@tudelft.nl - Academia Collu maurizio.collu@cranfield.ac.uk TBC EOLFI France Industry Paillard benoit.paillard@eolfi.com Nenuphar France Industry Kluczewska- Bordier joanna.kluczewska-bordier@nenupharwind.com DTU Denmark Academia Paulsen uwpa@dtu.dk - Sandia National Laboratories Ecole Centrale de Nantes United States Government Griffith dgriffi@sandia.gov TBC France Industry/Acade mia Gilloteaux jean-christophe.gilloteaux@ec-nantes.fr VUB Belgium Academia Runacres mark.runacres@vub.ac.be TBC UCL Belgium Academia Chatelain philippe.chatelain@uclouvain.be

9 Participants codes summary Institution Name Type Politecnico di Milano (IT) Ecole centrale de Nantes (FR) Active Pitch OpenFoam CFD-LES Innwave/ Cactus Vortex Panel? Reference Ilmas Bayati Luca Bernini Jean Christophe Gilloteaux ec-nantes.fr starccm+ CFD-URANS? Antoine Ducoin EOLFI (FR) OpenFoam CFD-URANS / DES Cactus Vortex Panel? Benoit Paillard benoit.paillard@eolfi.com Cranfield University (UK) TM4E BEM?? CFD RANS-DES? Pierre-Luc Delafin p.p.delafin@cranfield.ac.uk NENUPHAR (FR) PHARWEN3D Vortex Simon Horb simon.horb@nenupharwind.com

10 Participants codes summary Institution Name Type Added Info Input Output Inputs Outputs Politecnico di Milano (IT) Ecole centrale de Nantes (FR) OpenFoam Innwave/ Cactus CFD-LES Vortex Panel Actuator Line with effective velocity method EVM for the definition of the AoA. No Dynamic stall 1,2,3 1,2-1,2,3 1,(2) starccm+ CFD-URANS - 1,2,4,5 1,2 1. Geometry, 2. Operational Conditions and Inflow field 3. Aero coefficients 4. Turbulence intensity and length scale 5. Experimental test section geometry 1. Aerodynamic Force and Pressure on blades 2. Wake velocity and pressure field EOLFI (FR) OpenFoam Cactus CFD-URANS / DES Vortex Panel - 1,2,4 1,2-1,2,3 1,(2) Cranfield University (UK) TM4E? BEM CFD RANS- DES Gormont-Berg dynamic stall model, 3D effects, rigid Platform Motion 1,2,3 1,2-1,2 1,2 NENUPHAR (FR) PHARWEN3 D Vortex Coupled with structural beam code (NeSTor), rigid platform motion. Also coupled with inviscid with dynamic stall (Beddoes-Leishman) 1,2,3,4 1,2

11 Phase I: parameters, symbols and conventions Symbol Parameter Definition Value Unit R Rotor radius - 1 [m] N Number of blades - 3 [-] U c Chord length [m] ω Angular Velocity - - [rad/s] ρ Air density [kg/m3] ν Kinematic Viscosity μ ρ 1e-5 [m^2/s] U Undisturbed wind speed - - [m/s] λ Tip Speed Ratio ω R U - [-] σ Solidity N c 2R [-] Cp Power coefficient Power 0.5ρ 2RU^3 - [-] The plane OXY is supposed at half wind turbine height H Ct Thrust coefficient Thrust 0.5ρ 2RU^2 - [-]

12 Phase I: general scheme Phase A: steady warm up Phase B: deeper analysis at max Cp 1D 3D 10D Phase C: with control Phase D : final simulations and 1 comparison with exp. Warm up for Phase II φ t = φ 0 sin(θ t )

13 Phase I simulations: steady warm up Phase Simulations setup Deliverables DUE DATE A No tower, No braces V = 5 m/s Different Tip Speed Ratios λ : from 1 to 5, with step 0.5 Cl and Cd given by PoliMi Cp(λ), Ct(λ) CASES (1) Fixed Pitch φ = 0 15 Feb 2017 (2) Fixed Pitch φ = 5/+5 31 March 2017 Recommendations Cp and Ct are considered mean values at each λ (difference between mean value between a complete rotation and the previous one is less than at least 2-3% ). Provide also max and min values of Cp and Ct within a complete rotation ϑ=[1:360] for each λ

14 Phase I simulations: deeper analysis Phase Simulations setup Deliverables DUE DATE B No tower, No braces, only blades V = 5 m/s Cl and Cd given by PoliMi Cp(ϑ) Ct(ϑ) Ft(ϑ) Detailed analysis for a single Tip Speed Ratio Cp max (e.g. 3) Probes in the wake: distance 1D, 3D and 10 D Fn(ϑ) α ϑ Vr ϑ Wake velocity u,v,w 31 May 2017 Recommendations ϑ range:1 to 360, increment 1 Wake probes sampling: x = -3.75m to 3.75 m, increment 0.25m y = -1D, -3D, -10D z = 0 m α ϑ and Vrel ϑ are requested from those participants who can easily extract these info from their simulations. The forces on the blade are intended at mid-span

15 Phase I simulations: with control Phase Simulations setup Deliverables DUE DATE C V = 5 m/s Cp(ϑ) Cl and Cd from measurements (PoliMi and VUB) Ct(ϑ) Ft(ϑ) Detailed analysis for a single Tip Speed Ratio Cp max (e.g.3) Fn(ϑ) α ϑ 1 Sept 2017 Vr ϑ Probes in the wake: distance 1D, 3D and 10 D Wake velocity u,v,w Control laws: φ t = φ 0 sin θ t with φ 0 =2.5 / 5 / 7.5

16 Phase I simulations: final simulations and 1 comparison with exp. Phase Simulations setup Deliverables DUE DATE D Tower included (200mm) Cl and Cd from measurements (PoliMi and VUB) Static curves (Based on experimental tests) Cp(λ) Ct(λ) Cp(ϑ) Ct(ϑ) Ft(ϑ) 15 November 2017 Detailed analysis for a single TSR Cp max (Based on experimental tests) Probes in the wake distance xd (Based on experimental tests) Fn(ϑ) α ϑ Vr ϑ Wake velocity u,v,w 1P Control laws (Based on experimental tests)

17 Open Discussion Proposed Phase I outline Agreeing on the details of the Phase I/A Next Net Meetings References Website: To be included in the VODCA project, send an to ilmasandrea.bayati@polimi.it reporting the details of the numerical code adopted