kth.sePublications
Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
Numerical study of the hydrodynamic stability of a wind-turbine airfoil with a laminar separation bubble under free-stream turbulence
KTH, School of Engineering Sciences (SCI), Engineering Mechanics, Fluid Mechanics and Engineering Acoustics, Fluid Physics. KTH, Centres, SeRC - Swedish e-Science Research Centre.ORCID iD: 0009-0007-8056-6109
Mechanical Engineering Department, Kiel University of Applied Sciences 2 , D-24149 Kiel, Germany.ORCID iD: 0000-0002-7568-8757
Aerospace Engineering, Monash University 3 , 3800 Clayton, Australia.ORCID iD: 0000-0001-7831-8121
Mechanical Engineering Department, Kiel University of Applied Sciences 2 , D-24149 Kiel, Germany.ORCID iD: 0000-0002-9357-3232
Show others and affiliations
2023 (English)In: Physics of fluids, ISSN 1070-6631, E-ISSN 1089-7666, Vol. 35, no 8Article in journal (Refereed) Published
Abstract [en]

The interaction of several instabilities and the influence of free-stream turbulence on laminar-turbulent transition on a 20% thick wind-turbine blade section with a laminar separation bubble (LSB) are investigated with wall-resolved large-eddy simulations (LES). Turbulence intensities (TI) of 0%, 2.2%, 4.5%, 8.6%, and 15.6% at chord Reynolds number 100,000 are considered. Linear receptivity occurs for the most energetic disturbances; high-frequency perturbations are excited via non-linear mechanisms for  TI≥8.6%⁠. Unstable Tollmien–Schlichting (TS) waves appear in the inflectional flow region for  TI≤4.5%⁠, shifting to inviscid Kelvin–Helmholtz (KH) modes upon separation and forming spanwise rolls. Sub-harmonic secondary instability occurs for  TI=0%⁠, with rolls intertwining before transition. Streaks spanwise modulate the rolls and increase their growth rates with TI for  TI≤4.5%⁠, reducing separation and shifting transition upstream. The  TI=4.5% case presents the highest perturbations, leading to the smallest LSB and most upstream transition. Earlier inception of TS/KH modes occurs on low-speed streaks, inducing premature transition. However, for  TI=8.6%⁠, the effect of the streaks is to stabilize the attached mean flow and front part of the LSB. This occurs due to the near-wall momentum deficit alleviation, leading to the transition delay and larger LSB than  TI=4.5%⁠. This also suppresses separation and completely stabilizes TS/KH modes for  TI=15.6%⁠. Linear stability theory predicts well the modal evolution for  TI≤8.6%⁠. Optimal perturbation analysis accurately computes the streak development upstream of the inflectional flow region but indicates higher amplification than LES downstream due to the capture of low-frequency, oblique modal instabilities from the LSB. Only low-amplitude [ O(1%)] streaks displayed exponential growth in the LES since non-linearity precludes the appearance of these modes.

Place, publisher, year, edition, pages
AIP Publishing , 2023. Vol. 35, no 8
National Category
Fluid Mechanics
Identifiers
URN: urn:nbn:se:kth:diva-337104DOI: 10.1063/5.0159783ISI: 001041183200006Scopus ID: 2-s2.0-85166950584OAI: oai:DiVA.org:kth-337104DiVA, id: diva2:1800147
Funder
StandUp for WindStandUp
Note

QC 20230926

Available from: 2023-09-25 Created: 2023-09-25 Last updated: 2025-02-09Bibliographically approved
In thesis
1. Stability and transition on wind turbine blades
Open this publication in new window or tab >>Stability and transition on wind turbine blades
2023 (English)Doctoral thesis, comprehensive summary (Other academic)
Abstract [en]

Transition on wind turbine blades is a highly complex phenomenon due to the myriad effects influencing the process. This thesis studies some of them, namely free-stream turbulence (FST), rotation, and three-dimensionality. The investigations employ large eddy simulations (LES) with and without (implicit or wall-resolved LES) a subgrid-scale model. The role of FST in the modal and non-modal stability of the flow on the suction side of a wind turbine section at a Reynolds number 𝑅𝑒𝑐 = 100,000 is studied. This involved several simulations at varying turbulence intensity (𝑇𝐼) and primary and secondary linear stability analyses. The separated shear layers strongly govern the flow stability with the characteristic Kelvin-Helmholtz (KH) modes. Low FST levels increase the growth rates of the secondary instability of Tollmien-Schlichting (TS) and KH modes, leading to an upstream shift of transition and shrinking of the LSB. High enough 𝑇𝐼 stabilizes the flow to these modes, leading to an unexpected increase in the LSB. However, further rises in the turbulence level suppress separation. The spanwise-averaged part of the mean-flow distortion causes the stabilizing effect. The increase in the turbulence intensity also leads to a monotonic drop in the energy of coherent structures, shed from the separated shear layer, passing near the trailing edge. In the case of 𝑅𝑒𝑐 = 1,000,000, streak growth is much more intense, and even low levels of 𝑇𝐼 are enough to suppress the LSB present in the absence of FST. For 𝑇𝐼 ≤ 2.4%, transition is caused by inner modes, which in the limit of zero FST tend to TS waves. This range of 𝑇𝐼 presents linear receptivity, good agreement of the 𝑁 factor from Mack’s correlation with simulation data, and an exponential dependency of the transition location with 𝑇𝐼. For 2.4% < 𝑇𝐼 ≤ 7.0%, bypass transition occurs, characterized by the predominance of the outer varicose mode. In this regime, the transition location displays a variation ∝ 𝑇𝐼−2. A low-frequency cut-off for the free-stream turbulence is proposed to allow the computation of an effective turbulence intensity for wind turbine blades. Regarding the role of rotation, a model is developed to compute the quasi-three-dimensional base flow for stability analyses over a blade. The flow in the inboard region is highly three-dimensional and significantly affected by rotation. Highly oblique modes are the most unstable in this region, leading to a transition up to 19% earlier than the widely used two-dimensional semi-empirical 𝑒𝑁 transition model of Drela and Giles, used in the RANS simulations. A transition-prediction framework based on the boundary layer and parabolized stability equations accounting for these effects was developed. It indicates that rotation shifts transition upstream if the Reynolds number is allowed to increase with the reference velocity. Subsequent LES indicated that rotation stabilizes the flow for a fixed Reynolds number in the attached flow region and front part of the LSB for low rotation rates, delaying transition and reattachment. Even though rotation delays these phenomena, rotation may act as an adverse pressure gradient after separation occurs, leading to an increase in the growth rates of the KH modes and reverse flow. Furthermore, crossflow transition may be triggered for higher rotation rates and towards the inboard blade region, leading to an upstream shift of the transition point. Crossflow transition leads to a rise in the pressure difference between the two sides of the airfoil, generating a higher lift. In the outboard blade region, a self-excited type of instability may occur in an LSB forming near the leading edge, promoting an early transition that may cause a sudden shift of the separation line to the leading edge after a certain critical radius, as observed in experiments. Finally, a low-frequency oscillation in the normal force coefficient, with an amplitude of 10.5% around the mean, was identified in a wind turbine airfoil. The period of these oscillations was long, corresponding to several turns of a wind turbine at rated rotation speed. The occurrence of such a phenomenon in real wind turbines should be assessed and considered in the structural design of the rotor.

Abstract [sv]

Laminär-turbulent omslag på vindturbinblad är ett mycket komplext fenomen på grund av de många faktorer som påverkar processen. I denna avhandling studerar vi några av dem, nämligen hur friströmsturbulens (FST), rotation och tredimensionalitet påverkar strömningen. Analysen är baserad på det så kallade large eddy simulation (LES) med och utan en modell för subgrid-skala (implicit eller väggupplöst LES). FST:s roll i den modala och icke-modala stabiliteten av flödet på sugsidan av en vindturbinsektion vid ett Reynolds nummer Rec=100.000 studeras. Detta involverade flera simuleringar med varierande turbulensintensitet (TI) och primära och sekundära linjära stabilitetsanalyser. De separerade skjuvskikten styr flödesstabiliteten via de karakteristiska Kelvin-Helmholtz (KH) virvlarna. Låga FST-nivåer ökar tillväxthastigheten för den sekundära instabiliteten i Tollmien-Schlichting (TS) och KH-fallet, vilket leder till en uppströmsförskjutning av omslaget till turbulens och en minskning av LSB. Tillräckligt hög TI stabiliserar flödet i dessa lägen, vilket leder till en oväntad ökning av LSB. Ytterligare höjningar av turbulensnivån undertrycker emellertid avlösningen. Den spännvidds medelvärdesbildade delen orsakar den stabiliserande effekten. Ökningen av turbulensintensiteten leder också till en monoton minskning av energin hos koherenta strukturer, bakom det separerade skjuvskiktet. I fallet med Rec=1.000.000 är stråktillväxten mycket mer intensiv, och även låga nivåer av TI är tillräckligt för att undertrycka LSB som finns i frånvaro av FST. För TI≤ 2,4% orsakas övergången av inre moder som liknar TS-vågor i gränsen mot noll FST. Detta intervall av TI visar god överensstämmelse mellan N-faktorn från Macks korrelation med simuleringsdata och ett exponentiellt beroende av omslaget med TI. För 2,4%<TI≤7,0% sker sk  bypass omslag, kännetecknad av dominansen symmetriska stråk. En lågfrekvent cut-off för friströmsturbulensen möjliggör beräkning av en effektiv turbulensintensitet för vindturbinblad. När det gäller rotationens roll utvecklas en modell för att beräkna det kvasi-tredimensionella basflödet för stabilitetsanalyser över ett blad. Flödet är mycket tredimensionellt och påverkas avsevärt av rotation. Mycket sneda vågor är de mest instabila i denna region, vilket leder till ett omslag upp till 19% tidigare än den allmänt använda tvådimensionella semi-empiriska eN omslagsmodellen av Drela och Giles, som används i RANS-simuleringarna. Ett ramverk för omskagsprediktion baserat på gränsskiktet och paraboliserade stabilitetsekvationer som står för dessa effekter utvecklades. Det indikerar att rotation skiftar övergång uppströms om Reynolds-talet tillåts öka med referenshastigheten. Efterföljande LES indikerade att rotation stabiliserar flödet för ett fast Reynolds-tal i flödesområdet och främre delen av LSB för låga rotationshastigheter, vilket fördröjer övergång och återlaminarisering. Även om rotation fördröjer dessa fenomen, kan rotation verka som en negativ tryckgradient efter att separation inträffar, vilket leder till en ökning av tillväxthastigheterna för KH-moden och återcircirkulerande flöde. Vidare kan omslag till turbulens utlösas för högre rotationshastigheter och mot det inre bladområdet, vilket leder till en uppströmsförskjutning av omslagspunkten. Crossflow-övergång leder till en ökning av tryckskillnaden mellan de två sidorna av vingytan, vilket genererar en högre lyftkraft. I andra områden kan en självexciterad typ av instabilitet uppstå i en LSB som bildas nära framkanten, vilket främjar en tidigt omslag som kan orsaka en plötslig förskjutning av separationslinjen till framkanten efter en viss kritisk radie, som också observerats i experiment. Slutligen identifierades en lågfrekvent oscillation i normalkraftskoefficienten, med en amplitud på 10,5% runt medelvärdet, i ett vindkraftverks vingprofil. Perioden för dessa svängningar var lång, motsvarande flera varv av ett vindturbin med nominell rotationshastighet. Förekomsten av ett sådant fenomen i verkliga vindkraftverk bör bedömas och beaktas i den strukturella designen av rotorn.

Place, publisher, year, edition, pages
Stockholm: KTH Royal Institute of Technology, 2023. p. 61
Series
TRITA-SCI-FOU ; 2023:47
Keywords
Boundary layer stability, free-stream turbulence, laminar-turbulent transition, low-frequency oscillations, rotation effects, wind turbine blades., Gränsskiktsstabilitet, friströmsturbulens, laminär-turbulent omslag, lågfrekventa svängningar, rotationseffekter, vindkraftverksblad.
National Category
Fluid Mechanics
Research subject
Engineering Mechanics
Identifiers
urn:nbn:se:kth:diva-337136 (URN)978-91-8040-701-4 (ISBN)
Public defence
2023-10-23, Kollegiesalen, Brinellvägen 8, Stockholm, 10:00 (English)
Opponent
Supervisors
Funder
StandUpStandUp for Wind
Note

QC 231002

Available from: 2023-10-02 Created: 2023-09-26 Last updated: 2025-02-09Bibliographically approved

Open Access in DiVA

No full text in DiVA

Other links

Publisher's full textScopus

Authority records

Coelho Leite Fava, ThalesHenningson, Dan S.Hanifi, Ardeshir

Search in DiVA

By author/editor
Coelho Leite Fava, ThalesLobo, B. A.Nogueira, P. A. S.Schaffarczyk, A. P.Breuer, M.Henningson, Dan S.Hanifi, Ardeshir
By organisation
Fluid PhysicsSeRC - Swedish e-Science Research CentreFluid Mechanics and Engineering AcousticsLinné Flow Center, FLOW
In the same journal
Physics of fluids
Fluid Mechanics

Search outside of DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetric score

doi
urn-nbn
Total: 200 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf