Benard, Pierre

Bricteux, Laurent

Houtin-Mongrolle, Felix

Tene Hedje, Patrick

Vigny, Ulysse

Zeoli, Stéphanie

ORBi UMONS

[1] Mertens S 2002 Refocus 3 22-24.[2] Ge et al. 2021 Ren. En. 163 1063-77[3] Moureau et al. 2011 Comptes Rendus Mécanique[4] Bénard et al. 2018 Computers & Fluids[5] Nicoud et al. 2011. Physics of Fluids[6] Li et al. 2015 J. Wind Eng. Ind. Aerodyn. 146 51-58.[7] Vasaturo et al. 2018 J. Wind Eng. Ind. Aerodyn. 173 241–61.[8] Hrastinski et al. 2015 Croatian Meteorological Journal 50 105-120..5-REFERENCES & ACKNOWLEDGEMENTS2- METHODOLOGY & SETUPLARGE EDDY SIMULATION OF HAWT AND VAWT PERFORMANCES IN THE VICINITY OF A BUILDINGP Tene Hedje1, S Zeoli1, U Vigny1,2, F Houtin-Mongrolle2, P Benard2, L Bricteux11 Fluids-Machines Unit, Faculty of engineering, UMONS2 CNRS UMR 6614, INSA and University of Rouen, CORIApatrick.tenehedje@umons.ac.beThis study has received funding from the POPE research project Pôle Mecatech, infrastructure funded by the Walloon Region under the grantagreement No.7961. This research benefited from computational resources were made available on the Tier-1 supercomputer of theFédération Wallonie-Bruxelles, infrastructure funded by the Walloon Region under the grant agreement No.1117545. Computational resourceswere also provided by the Consortium des Equipements de Calcul Intensif, funded by the Fonds de la Recherche Scientifique de Belgiqueunder Grant No.2.5020.11 and by the Walloon Region.1- CONTEXT 3- RESULTS4- CONCLUSION & PERSPECTIVESObjectivesProblematicsTurbulent inflow generationNumerical model• 𝑢 𝑧 = !!"log #$#"#"• 𝑢%!& = 5𝑚/𝑠 ;	κ = 0.37Channel flow Precursor Database Method [6,7] Wall Modeled Large Eddy Simulations (WMLES) • Low Mach number N-S equations solver• 4th order central finite volume method, 4th order time integration• Turbulence model: 𝝈-model [5]The wind flow in urban environment has an unpredictable behaviour [2]:• complex turbulent phenomena at the bottom of the boundary layer• constantly changing wind directionWhich turbine (VAWT vs HAWT) performs best in this context?→ high fidelity numerical simulation to better understand how wind turbineswork in complex environments.Cases(color)WindTurbineBuilding WinddirectionIncident angle[°]V0 (-----)VAWTNo - -V1 (-----) Yes SW 45V2 (-----) Yes SSW 67.5H1 (      )HAWTNo - -H1 (      ) Yes SW 45H2 (      ) Yes SSW 67.5V0 H0 V1 H1 V2 H2CP 0.4 0.42 0.42 0.46 0.42 0.45CT 0.72 0.71 0.75 0.75 0.75 0.74The current energetic and environmental challenges induce a significantgrowth of the use of renewable energies. Wind energy production isundergoing intense research and development work. Progress in thedevelopment of wind energy in urban areas [1] has led to competitionbetween Vertical Axis Wind Turbines (VAWTs) and Horizontal Axis WindTurbines (HAWTs).[3,4]• HAWT : NREL5MW scaled; 𝜆/01 = ⁄𝜔 𝑟 𝑢234 = 7.55• VAWT : 3-bladed H-rotor; generic NACA (𝑐 = 0.45𝑚) 𝜆/01 = 3.65• ALM : 32 points per blade; ⁄𝜖 ∆𝑥 = 2Wind TurbineBuilding & simulation casesFlow statisticsPower and ThrustTime-averaged streamwise velocity and Turbulent Kinetic Energy (TKE)Flow topology• Horizontal profiles (𝑧 = 20𝑚 ) at several downstream locations of the rotor• Local overspeed up to	 ⁄𝑥 𝐷 = 3 (~4% for SW and ~2.4% for SSW cases): low-rise building• Higher TKE Levels for SW and SSW cases• In the near wake 𝑇𝐾𝐸5671 ≤ 𝑇𝐾𝐸8671• Vertical slices colored by time-averaged streamwise velocity and TKE• Slight flow deviation in SW and SSW cases (greater in SW cases)• Higher production for HAWT overall.• 𝐶9 (𝑉1 and 𝑉2) = 𝐶9 (H0) • Higher production for SW and SSW cases→ positive effect of the buildingPerspectives :Methodology :• WMLES to access aerodynamic performances• Modelling the flow in a realistic way• Highlight the effect of the buildingThe building locally disturbs the flow. It induces overspeed, which improves the power production of the wind turbines, but also leads to high turbulence levels, so the blades are subjected to more loads. HAWTs however have better performances                   Industrial building1zhub =20mDHawt = 9.388 m DVawt = 11.2 mHVawt =9mH = 8 my/D Hawt[!]u/U ref[!]y/D Vawt[!]u/U ref[!]u/U ref[!]x /DVawt [ ! ] x /DHawt [ ! ]y/D Vawt[!]y/D Vawt[!]y/D Hawt[!]y/D Hawt[!]x /DVawt [ ! ]x /DVawt [ ! ]x /DHawt [ ! ]x /DHawt [ ! ]V0V1V2H0H1H2TKE/U2 ref[!]TKE/U2 ref[!]TKE/U2 ref[!]V0V1V2H0H1H2y/D Hawt[!]y/D Vawt[!]x /DVawt [ ! ] x /DHawt [ ! ]y/D Vawt[!]y/D Vawt[!]y/D Hawt[!]y/D Hawt[!]x /DVawt [ ! ]x /DVawt [ ! ]x /DHawt [ ! ]x /DHawt [ ! ]Faculté PolytechniqueThermique - FondementsBAB3 - Bachelier Ingénieur civilV. Feldheim, P. LybaertService de Thermique et combustionS. Zeoli (stephanie.zeoli@umons.ac.be)Service de Fluides-MachinesFaculté P lytechniq eTh rmique - FondementsBAB3 - Bachelier Ingénieur civilV. Feldheim, P. LybaertService de Thermique et c mbustionS. Zeoli (stephanie.zeoli@umons.ac.be)Service de Fluides-MachinesCONFIDENTIEL	 	 Page	1	                       RAPPORT SEMESTRIEL     POPE - Convention n° 7961 Pôle MecaTech –Appel N° 23  Semestre n° 6  -  Période couverte : du 01/06/2021 au 30/11/2021        WaPT UCLouvain 3E UMons Fairwind   Coordinateur du projet Ivan DE VISSCHER (WaPT)  Ivan.devisscher@wapt.be  +32 10 392 245   Date de démarrage du projet 01/12/2018 Durée du projet 36 mois a supprimé: ¶CONFIDENTIEL	 	 Page	1	                       RAPPORT SEMESTRIEL     POPE - Convention n° 7961 Pôle MecaTech –Appel N° 23  Semestre n° 6  -  Période couverte : du 01/06/2021 au 30/11/2021        WaPT UCLouvain 3E UMons Fairwind   Coordinateur du projet Ivan DE VISSCHER (WaPT)  Ivan.devisscher@wapt.be  +32 10 392 245   Date de démarrage du projet 01/12/2018 Durée du projet 36 mois a supprimé: ¶→ Investigates of other wind directions→ Influences of the upstream buildings→ Noise pollution and fatigue effects• Low-rise building (H	=	8	m) • Typical industrial area (TI		~	17	%	)• Wind rose obtained from IRM meteorological model ALADIN [8] • Two relevant wind directions: SW and SSW0.4 0.6 0.8 1.0!1.5!1.0!0.50.00.51.01.5y/DVawt[!]x/DV awt = 00.6 0.8 1.0x/DV awt = 10.6 0.8 1.0x/DV awt = 30.6 0.8 1.0 1.2x/DV awt = 5u/Uref [!]0.4 0.6 0.8 1.0!1.5!1.0!0.50.00.51.01.5y/DVawt[!]x/DV awt = 00.6 0.8 1.0x/DV awt = 10.6 0.8 1.0x/DV awt = 30.6 0.8 1.0 1.2x/DV awt = 5u/Uref [!]0.00 0.02 0.04!1.5!1.0!0.50.00.51.01.5y/DVawt[!]x/DV awt = 00.02 0.04x/DV awt = 10.02 0.04x/DV awt = 30.02 0.04 0.06x/DV awt = 5TKE/U2ref [!]0.00 0.02 0.04!1.5!1.0!0.50.00.51.01.5y/DHawt[!]x/DHawt = 00.02 0.04x/DHawt = 10.02 0.04x/DHawt = 30.02 0.04 0.06x/DHawt = 5TKE/U2ref [!]Comparison of the aerodynamic performance of a VAWT and a HAWT located close to a full-scale low-rise industrial building, for various wind directions• 𝑧%!& = 20𝑚 ; 𝑧: = 0.06𝑚• TI	[%]	~	15	- 200.0 0.5 1.0 1.5u/Uref0.00.20.40.60.81.0z+z 0z max+z 0Log ! lawPrecursor5 10 15 20TI [%]0.00.20.40.60.81.0z+z 0z max+z 0• Design criteria :  𝑃3;3< 8=>? = 𝑃3;3< 5=>? for 𝑢@= 𝑢2340 5 10 15 20U! [m/s]051015202530Pelec[kW]V AWTHAWTuref

Large Eddy Simulation of HAWT and VAWT performances in the vicinity of a building

https://orbi.umons.ac.be/bitstream/20.500.12907/45051/1/Poster_Tene_Hedje_et_al_Torque_2022.pdf

Large Eddy Simulation of HAWT and VAWT performances in the vicinity of a building

Abstract

Similar works

Full text

Available Versions

ORBi UMONS