Numerical simulations of transition and long-term response of a wind turbine airfoil

Diese Arbeit präsentiert numerische Simulationen unter Verwendung von Nek5000 und EllipSys zur Analyse des Übergangs und der langfristigen Strömungsreaktion eines FFA-W3-Tragflügels bei niedrigen Reynolds-Zahlen, wobei die Fähigkeit des EllipSys-Codes zur Erfassung von Kelvin-Helmholtz-Instabilitäten validiert und eine langsame, periodische Modulation des Normalkraftkoeffizienten identifiziert wird, die mit niederfrequenten Oszillationen und potenziellem Aufbrechen von Blasen zusammenhängt.

Das Wesentliche

Stellen Sie sich ein Windkraftanlagen-Rotorblatt wie einen riesigen, rotierenden Flügel vor. Genau wie ein Flugzeugflügel benötigt es einen glatten Luftstrom, um effizient zu arbeiten. Doch wenn der Wind in bestimmten Winkeln auf den Flügel trifft, kann die Luft „hängenbleiben“ und von der Oberfläche abreißen, was ein chaotisches, wirbelndes Durcheinander erzeugt. Diese Arbeit ist wie ein hochmoderner Windkanal-Versuch, aber anstatt eines physischen Modells haben die Forscher ein virtuelles Modell in einem Supercomputer gebaut, um genau zu beobachten, wie sich diese Luft verhält.

Hier ist die Geschichte ihrer Entdeckung, unterteilt in einfache Konzepte:

1. Der virtuelle Windkanal

Die Forscher wollten einen spezifischen Ausschnitt eines massiven Windkraftanlagen-Rotorblatts (einer 10-MW-Turbine) untersuchen. Sie verwendeten zwei verschiedene Computerprogramme, NEK5000 und ELLIPSYS, um die Luftströmung über dieses Blatt zu simulieren.

Stellen Sie sich NEK5000 als eine hochklassige, ultrapräzise Kamera vor, die jedes winzige Detail erfasst, aber sehr langsam und teuer in der Ausführung ist. ELLIPSYS ist wie eine etwas schnellere, effizientere Kamera. Das Team musste zuerst beweisen, dass die „schnellere“ Kamera (ELLIPSYS) dieselben Dinge sehen kann wie die „hochklassige“. Sie fanden heraus, dass ELLIPSYS zwar ein paar winzige, schwache Kräuselungen in der glatten Luft übersah (weil es die Dinge ein wenig zu stark glättete), aber exzellent darin war, die großen, chaotischen Wirbel zu erfassen, die tatsächlich für die Leistung des Blattes entscheidend sind.

2. Wie breit muss der Tunnel sein?

Bevor sie die langen Simulationen starteten, mussten sie herausfinden, wie breit ihr virtueller „Windkanal“ sein müsste. Wenn der Tunnel zu schmal ist, könnte er die Luft zusammenpressen und falsche Ergebnisse erzeugen. Wenn er zu breit ist, verschwendet er Rechenleistung.

Sie testeten einen „schmalen“ Tunnel (10 % der Breite des Flügels) gegen einen „breiten“ Tunnel (20 % der Breite).

• Die Analogie: Stellen Sie sich vor, Sie beobachten einen Fluss. Wenn Sie nur einen schmalen Streifen des Flusses betrachten, übersehen Sie dann die großen Wellen?
• Das Ergebnis: Sie fanden heraus, dass der schmale Tunnel tatsächlich ausreichte. Die großen Wellen und Wirbel bildeten sich im engen Raum perfekt aus. Dies bedeutete, dass sie viel Computerzeit sparen konnten, indem sie den kleineren, schmaleren Simulationsraum nutzten.

3. Die „Blase“ und die „Klappe“

Der interessanteste Teil der Studie ereignete sich auf der Oberseite des Flügels (der Saugseite).

• Die Ablöseblase: Während die Luft über den Flügel strömt, löst sie sich für einen Moment ab und erzeugt eine kleine, rezirkulierende Lufttasche, die sogenannte „Laminare Ablöseblase“ (Laminar Separation Bubble, LSB). Denken Sie an einen winzigen, vorübergehenden Wirbel auf der Oberfläche des Flügels.
• Die Instabilität: Innerhalb dieser Blase ruht die Luft nicht einfach still; sie vibriert und rollt sich zu Wellen auf (wie Kräuselungen auf einem Teich). Die Forscher beobachteten, wie diese Wellen wuchsen. Sie fanden heraus, dass der Haupt„Roller“ in dieser Blase eine Art von Instabilität ist, die als Kelvin-Helmholtz-Mode bezeichnet wird.
• Die Entdeckung: Sie bestätigten, dass das „schnellere“ Computerprogramm (ELLIPSYS) diese Wellenbewegungen und das Verhalten der Blase genau vorhersagen konnte, was mit den Ergebnissen des ultrapräzisen Programms übereinstimmte.

4. Der langsame, rhythmische Puls (Die große Überraschung)

Nachdem sie ihre Werkzeuge validiert hatten, ließen sie die Simulation über eine sehr lange Zeit laufen (entsprechend 50 Mal dem Durchströmen des Flügels durch die Luft). Hier entdeckten sie etwas Besonderes.

Während die Luft mit schnellen, chaotischen Bewegungen wirbelte, bemerkten sie einen sehr langsamen, rhythmischen Puls in der Kraft, die auf den Flügel drückt.

• Die Analogie: Stellen Sie sich einen Trommelschlag vor. Das schnelle Wirbeln der Luft ist wie ein schnelles, hochfrequentes Trommelfeuer. Der langsame Puls, den sie fanden, ist wie ein tiefer, langsamer Herzschlag, der alle 48 Sekunden (in der Simulationszeit) auftritt.
• Der Effekt: Dieser langsame Herzschlag verursachte ein Auf und Ab der Kraft auf dem Flügel um etwa 10,5 %.
• Die Verbindung zu echten Windturbinen: Als sie dies auf eine echte, rotierende Windkraftanlage übertrugen, erkannten sie, dass dieser langsame Puls einmal alle 7,7 vollen Umdrehungen des Blattes auftritt.

5. Warum passiert das?

Die Forscher glauben, dass dieser langsame Puls durch einen Zyklus verursacht wird, in dem die Luft am Flügel „stallt“ (feststeckt/abreißt) und wieder „unstallt“ (loslässt).

• Der Zyklus: Die Luft bleibt hängen und erzeugt eine große Blase. Dann wird etwas ausgelöst, das die Blase zum Platzen bringt, und die Luft strömt wieder glatt an. Dann baut sich der Druck erneut auf, die Blase bildet sich, und der Zyklus wiederholt sich.
• Der Auslöser: Sie vermuten, dass dies geschieht, weil die Luft auf dem Flügel so stark rückwärts wirbelt, dass sie einen Zustand der „absoluten Instabilität“ erzeugt – eine fachsprachliche Art zu sagen, dass die Luft so turbulent ist, dass sie zwangsläufig von selbst oszilliert.

6. Das Fazit

Diese Studie ist eine Erfolgsgeschichte für die Computermodellierung. Sie haben bewiesen, dass ein schnelleres, effizienteres Computerprogramm (ELLIPSYS) darauf vertraut werden kann, komplexe Windturbinen-Physik zu untersuchen, sofern man es zuerst gegen den „Goldstandard“ prüft.

Sie entdeckten, dass selbst auf einem dicken Windkraftanlagen-Rotorblatt eine sehr langsame, rhythmische „Atmung“ der Luftströmung stattfindet, die etwa alle 8 Umdrehungen des Blattes auftritt. Diese Atmung verursacht ein deutliches Ansteigen und Abfallen des Auftriebs (der Kraft, die die Turbine antreibt). Das Verständnis dieses langsamen Rhythmus ist entscheidend, denn obwohl es die Turbine nicht sofort beschädigen wird, könnten diese langsamen, großen Kraftschwankungen die Materialien über viele Jahre des Betriebs hinweg ermüden.

Kurz gesagt: Sie bauten einen virtuellen Windkanal, bewiesen, dass ein schneller Computer die Aufgabe erfüllen kann, und entdeckten, dass Windkraftanlagen-Rotorblätter einen langsamen, rhythmischen „Herzschlag“ haben, der durch das Entstehen und Platzen von Luftblasen auf ihrer Oberfläche verursacht wird.

Technische Zusammenfassung

Technisches Resümee: Numerische Simulationen des Übergangs und des Langzeitverhaltens eines Windkraftanlagen-Profils

Problemstellung
Strömungsablösung an Windkraftanlagen-Blättern, hervorgerufen durch erhöhte Anströmwinkel und dicke Profilgeometrien, beeinflusst die Leistung und die strukturelle Belastung erheblich. Vorangegangene Studien haben mehrere Eigenfrequenzen in abgelösten Strömungen identifiziert, einschließlich Wirbelablösung (Wake Shedding), Kelvin-Helmholtz (KH)-Moden und niederfrequenter Oszillationen (LFOs). LFOs sind durch starke Schwankungen der Auftriebskraft gekennzeichnet und werden häufig auf das periodische Stalling und Unstalling der Strömung zurückgeführt, was potenziell mit der absoluten Instabilität innerhalb laminarer Ablöseblasen (LSBs) zusammenhängt. Während LFOs in verschiedenen Profilstudien beobachtet wurden, erfordern ihr Auftreten bei niedrigen Anströmwinkeln (AoA) an dicken Windkraftanlagen-Profilen sowie die spezifischen Mechanismen, die sie unter diesen Bedingungen antreiben, weitere Untersuchungen. Diese Studie konzentriert sich auf den Übergang und das Langzeitverhalten eines Abschnitts des DTU 10-MW Referenz-Windkraftanlagen-Profils (FFA-W3 Serie) bei einer Sehnen-Reynolds-Zahl ($Re_c$) von $1 \times 10^5$.

Methodik
Die Forschung verwendet wandaufgelöste Large-Eddy-Simulationen (LES) ohne Subgrid-Scale-Modelle, um die inkompressiblen Navier-Stokes-Gleichungen zu lösen. Es werden zwei Solver eingesetzt:

1. NEK5000: Ein Spectral-Element-Methode (SEM) Code, der für minimale numerische Dissipation und hohe Genauigkeit bekannt ist.
2. ELLIPSYS: Ein Finite-Volumen-Code, der bei DTU/Risø entwickelt wurde, traditionell für RANS verwendet, hier jedoch als Implicit LES angewendet.

Die Studie beinhaltet eine systematische Validierung von ELLIPSYS gegenüber NEK5000. Die wesentlichen methodischen Schritte umfassen:

• Domänensensitivität: Bewertung der spanweiten Domänenbreite ($L_z$) durch Vergleich von Simulationen mit $L_z = 0,1c$ und $L_z = 0,2c$, um die Auflösung relevanter Strömungsstrukturen sicherzustellen.
• Solver-Validierung: Vergleich der mittleren Strömungsfelder, der Perturbationsentwicklung und der spektralen Charakteristika zwischen den beiden Solvern.
• Stabilitätsanalyse: Verwendung von Parabolisierten Stabilitätsgleichungen (PSE) und Spectral Proper Orthogonal Decomposition (SPOD), um das Wachstum von Instabilitäten (Tollmien-Schlichting- und KH-Moden) sowie die räumliche Struktur kohärenter Moden zu analysieren.
• Langzeitintegration: Nutzung des recheneffizienteren ELLIPSYS-Solvers, um die Strömungsevolution über 50 Durchströmzeiten ($T = 50c/U_\infty$) zu simulieren, um niederfrequente Phänomene zu erfassen.

Kernergebnisse

• Domänenbreite: Eine spanweite von 10 % der Sehne ($L_z = 0,1c$) erwies sich als ausreichend, um die Entwicklung der Hauptstörungen zu erfassen und die zeitgemittelten Strömungsstatistiken zu reproduzieren, was mit den Ergebnissen einer breiteren Domäne ($L_z = 0,2c$) übereinstimmt.
• Solver-Validierung: ELLIPSYS zeigte eine enge Übereinstimmung mit NEK5000 hinsichtlich der mittleren Strömungsfelder und der Entwicklung der am stärksten verstärkten Perturbationen. ELLIPSYS unterschätzte jedoch die Amplitude von Tollmien-Schlichting (TS)-Wellen in der anliegenden Grenzschicht aufgrund höherer numerischer Dissipation. Im Gegenzug sagte ELLIPSYS die Entwicklung der KH-Mode innerhalb der laminaren Ablöseblase (LSB) genau voraus, was gut mit den PSE-Vorhersagen übereinstimmte.
• Übergangsmechanismus: Der Übergang auf der Saugseite wird durch die KH-Instabilität getrieben, die sich auf der abgelösten Scherschicht entwickelt. Die am stärksten verstärkte KH-Mode hat eine Frequenz von $f \approx 15$, und ihre stromabwärtige Entwicklung wird von beiden Solvern sowie durch PSE gut erfasst.
• Niederfrequente Oszillationen (LFOs): Langzeitsimulationen ergaben eine langsame Modulation des Normalkraftkoeffizienten ($C_N$) mit einer Amplitude von 10,5 % und einer Periode von 48 Durchströmzeiten. Dies entspricht einer Frequenz von $f = 0,021$ und einer Strouhal-Zahl von $St = 0,0012$.
  • Diese Frequenz ist etwa 720 Mal niedriger als die primäre KH-Instabilitätsfrequenz.
  • Im Kontext der DTU 10-MW Windkraftanlage entspricht diese Periode etwa 7,7 Blattrotationen.
  • Die beobachtete $St$ ist niedriger als die in der Literatur für andere Profile berichteten Werte und tritt bei einem relativ kleinen effektiven Anströmwinkel von $3,1^\circ - 3,3^\circ$ auf.

Bedeutung und Behauptungen
Die Arbeit etabliert, dass ELLIPSYS ein lebensfähiges Werkzeug für wandaufgelöste LES von Übergangsströmungen ist und eine recheneffiziente Alternative zu Spectral-Element-Codes für Langzeitintegrationen bietet. Die Studie bestätigt, dass das periodische Stalling und Unstalling der Strömung, das potenziell durch absolute Instabilität aufgrund hoher Rückströmung (bis zu -17 % auf der Saugseite) ausgelöst wird, für die beobachteten LFOs verantwortlich ist.

Ein bedeutender Befund ist die Identifizierung von LFOs an einem dicken Profil (24 % Dicke) bei einem niedrigen Anströmwinkel, eine Bedingung, die sich von den für dünne Profile typischen Stalls unterscheidet, die mit solchen Phänomenen in der Literatur assoziiert werden. Die Autoren merken an, dass die Rückströmung auf beiden Seiten des Profils ausreichend ist, um eine absolute Instabilität zu ermöglichen, welche möglicherweise das periodische Aufbrechen der LSB in eine voll abgelöste Strömung antreibt. Die Studie trägt zum Verständnis darüber bei, wie die dicke Profilgeometrie den Beginn und die Charakteristika niederfrequenter Oszillationen beeinflusst, und hebt hervor, dass diese Phänomene bei niedrigeren Anströmwinkeln und Strouhal-Zahlen auftreten können als bisher dokumentiert für Standardprofile.