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 #2 del Duomo di Milano Aula Politecnico Magna di Rettorato Milano Mercoledì 27 maggio 2015 Department of Mechanical Engineering x = 0:.01:8*pi; s = x.*sin(x); c = x.*cos(x); plot3(s,c,x)

2 Presentation overview Project outline Details of participants numerical codes Phase IA(1) results Open discussion and updates

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 UPSALA Sweden Academia Mendoza Others?????

4 Phase IA active participants Institution Reference Politecnico di Milano (IT) Ecole centrale de Nantes (FR) Ilmas Bayati Luca Bernini Jean Christophe Gilloteaux Antoine Ducoin ec-nantes.fr EOLFI (FR) Benoit Paillard Cranfield University (UK) NENUPHAR (FR) Pierre-Luc Delafin Simon Horb UPSALA UNIV. (SE) Victor Mendoza

5 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

6 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 )

7 Participants codes summary (Phase IA1) Institution Name Type Added Info Blade Modelling 2D/3D/tip loss Pitch Control *please confirm Not used in this phase Fluid Dynamics Equations Politecnico di Milano (IT) OpenFoam Actuator Line Actuator Line with effective velocity method EVM for the definition of the AoA. No Dynamic stall LookUp Tables 2D CFD-LES Navier Stokes Ecole centrale de Nantes (FR) Innwave/ Cactus starccm+ Panel Code - LookUp Table 3D? CFD- URANS - Body Fitted Mesh 2D? Euler (Vorticity) CFD-URANS Navier Stokes EOLFI (FR) OpenFoam CFD- URANS / DES - Body Fitted Mesh 2D? Cactus Vortex - LookUp Table 3D? Navier Stokes Euler Cranfield University (UK) TM4E - BEM CFD RANS- DES Gormont-Berg dynamic stall model, 3D effects, rigid Platform Motion LookUp Table 2D? - Body Fitted Mesh 2D? Conservation of momentum Navier Stokes NENUPHAR (FR) PHARWEN3 D Panel Code Coupled with structural beam code (NeSTor), rigid platform motion. Also coupled with inviscid with dynamic stall (Beddoes-Leishman) Viscous Polars (?) 3D Euler (Vorticity)

8 Participants codes summary (Phase IA1) *please confirm Not used in this phase Institution Name Type Added Info Blade Modelling 2D/3D/tip loss Pitch Control Fluid Dynamics Equations - Vortex Lifting Line Dynamic stall (model) LookUp Table 3D? Euler (Vorticity) UPSALA - Actuator Line Dynamic Stall (model) LookUp Table 3D? CFD Navier Stokes (LES/DES/URAN S) UCL - 3D Vortex Particle Mesh (VPM) Immersed lifting line. Dynamic Stall (Leishmann- Beddoes) LookUp Table 3D? Lagrangian Approach

9 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 μ ρ 1.45e-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 - [-]

10 Phase IA results: Cp

11 Phase IA results: Cp Airfoil geometry Cl, Cd lookup table

12 Phase IA results: Cp Viscous Inviscid

13 Phase IA results: Ct

14 Phase IA results: Ct Airfoil geometry Cl, Cd lookup table

15 Phase IA results: Ct Viscous Inviscid

16 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 Further simulations: Fn, Ft, AoA as function of azimuthal angle θ only for λ =1.5 and 3 to better evaluate the differences of the phase (1) Cl and Cd given by PoliMi CASES (1) Fixed Pitch φ = 0 (2) Fixed Pitch φ = 5/+5 Cp(λ), Ct(λ) 15 May June 2017 Further outputs Fn, Ft, AoA as function of azimuthal angle θ

17 Open Discussion New simulation at 0 pitch? (IA1) Agreeing on the details of the Phase IA(2) Next Net Meetings Updates Website: (+FTP) ilmasandrea.bayati@polimi.it vodca@polimi.it Wind Tunnel days booked: November: VODCA participants are more than welcome if they can join PoliMi during the experiments