Two-way coupling of gravity waves and wind farm wakes: a reduced-order boundary-layer model

Dieser Beitrag stellt ein rechnerisch effizientes Reduced-Order-Modell vor, das die wechselseitige Kopplung zwischen Schwerewellen und Windparknachläufen erfolgreich erfasst, indem die nicht-hydrostatischen Boussinesq-Gleichungen linearisiert und die Dynamik der Grenzschicht und der freien Atmosphäre über eine Deckungsinversion gekoppelt werden, wobei eine Validierung gegenüber Large-Eddy-Simulationen die Fähigkeit des Modells bestätigt, wesentliche Strömungsmerkmale wie die Aufstauung stromaufwärts und eine beschleunigte Nachlaufwiederherstellung wiederzugeben.

Das Wesentliche

Stellen Sie sich einen riesigen Windpark nicht nur als Ansammlung sich drehender Turbinen vor, sondern als eine riesige, unsichtbare Hand, die in den Himmel greift und versucht, eine Handvoll Wind zu ergreifen, um Strom zu erzeugen. Diese Arbeit handelt davon zu verstehen, was passiert, wenn diese „Hand" so groß wird, dass sie den Wind nicht nur abfängt, sondern tatsächlich gegen die Atmosphäre selbst drückt und einen komplexen Tanz zwischen dem Wind und der Luft darüber erzeugt.

Hier ist die Geschichte dieses Tanzes, aufgeschlüsselt in einfache Konzepte:

Das Problem: Ein Windpark zu groß für seine Brisen

In der Vergangenheit waren Windparks klein genug, um wie ein paar Kieselsteine in einem Fluss zu wirken. Das Wasser (der Wind) floss leicht um sie herum, und der Fluss bemerkte sie kaum. Doch heute sind Windparks riesig – manchmal so hoch wie die gesamte Luftschicht, in der wir leben (die atmosphärische Grenzschicht).

Wenn ein so großer Windpark versucht, Energie aus dem Wind zu stehlen, verlangsamt er die Luft. Da Luft nicht einfach verschwinden kann, zwingt diese Verlangsamung die Luft, sich auf- und abzubewegen, um Platz zu schaffen. Denken Sie daran wie an eine überfüllte U-Bahn: Wenn alle plötzlich aufhören, sich vorwärts zu bewegen, müssen sie sich nach oben oder unten verschieben, um nicht gegeneinander zu stoßen.

Der „Trampolin"-Effekt (Schwerewellen)

Die Atmosphäre ist nicht nur leerer Raum; sie hat Schichten. Direkt über dem Windpark befindet sich eine deutliche „Decke", die Abschattungsinversion (capping inversion) genannt wird. Man kann sich diese Decke wie ein Trampolin oder eine schwere Decke vorstellen, die über dem Windpark gespannt ist.

Wenn der Windpark die Luft verlangsamt, drückt er die Luft nach oben und stößt dabei eine Wölbung in diese Trampolin-Decke.

1. Der Stoß: Der Windpark drückt die Luft nach oben.
2. Das Abprallen: Das „Trampolin" (die stabile Luft darüber) möchte zurückfedern. Dieses Zurückfedern erzeugt Wellen, die als Schwerewellen bekannt sind.
3. Das Feedback: Diese Wellen bleiben nicht einfach liegen; sie drücken zurück auf den Windpark. Es ist, als würde das Trampolin gegen Ihre Füße drücken. Dies erzeugt Druckänderungen, die den Wind entweder daran hindern können, die Turbinen zu erreichen (was sie weniger effizient macht), oder die Windgeschwindigkeit hinter dem Park beschleunigen können (was die Nachlaufströmung schneller erholen lässt).

Der alte Weg vs. der neue Weg

Der alte Weg (Der schwere Hammer):
Wissenschaftler verwendeten früher hochkomplexe Computersimulationen namens „Large Eddy Simulations" (LES), um dies zu untersuchen. Stellen Sie sich vor, Sie versuchen, jedes einzelne Luftmolekül und jede winzige Welle im Trampolin zu simulieren. Es ist unglaublich genau, erfordert aber so viel Rechenleistung, dass es so ist, als würde man versuchen, jedes Sandkorn an einem Strand zu zählen, nur um zu sehen, wie sich die Gezeiten bewegen. Es ist zu langsam für die Planung neuer Windparks oder deren Optimierung in Echtzeit.

Der neue Weg (Die intelligente Skizze):
Die Autoren dieser Arbeit haben ein „reduziertes Modell" (reduced-order model) erstellt. Denken Sie daran wie an eine intelligente Skizze anstelle eines fotorealistischen Gemäldes.

• Sie haben die Mathematik vereinfacht, indem sie sich nur auf die wichtigsten Teile konzentrierten: die vertikale Bewegung der Luft und die Wellen auf dem „Trampolin".
• Sie behandelten den Windpark als eine kontinuierliche Kraft, anstatt jede einzelne Turbinenschaufel zu simulieren.
• Sie verwendeten einen cleveren mathematischen Trick (eine Mischung aus spektralen und Finite-Differenzen-Methoden), um die Gleichungen schnell zu lösen.

Was sie herausfanden

Sie testeten ihre „intelligente Skizze" gegen den „schweren Hammer" (die superkomplexen Simulationen) und reale Daten. Hier ist, was sie entdeckten:

1. Die Blockade: Wenn sich der Windpark in einer stabilen Atmosphäre befindet (wie an einem ruhigen Tag mit einer klaren „Decke"), erzeugen die Schwerewellen einen „Gegenwind", bevor der Park überhaupt beginnt. Es ist, als würde man versuchen, gegen einen starken Gegenwind zu laufen, der sich bevor man das Hindernis erreicht bildet. Dies verlangsamt den Wind erheblich, bevor er die Turbinen erreicht.
2. Die Erholung: Hinter dem Windpark federt das „Trampolin" zurück und erzeugt einen „Rückenwind", der die Luft vorwärts drückt. Dies hilft, dass sich die Windgeschwindigkeit viel schneller erholt als an einem ruhigen, neutralen Tag.
3. Genauigkeit: Ihr vereinfachtes Modell stimmte fast perfekt mit den Ergebnissen der superkomplexen Simulationen überein, lief aber tausende Male schneller.

Das Fazit

Diese Arbeit bietet Ingenieuren ein schnelles, zuverlässiges Werkzeug, um vorherzusagen, wie riesige Windparks mit dem Himmel interagieren werden. Anstatt Tage darauf zu warten, dass ein Supercomputer ihnen sagt, wie ein Park funktionieren wird, können sie nun dieses Modell verwenden, um in Sekunden zu sehen, wie der „Trampolin"-Effekt der Atmosphäre dem Windpark helfen oder ihn behindern wird. Es überbrückt die Lücke zwischen einfachen Schätzungen und unmöglich ausführbaren Super-Simulationen und hilft uns, bessere Windparks zu entwerfen, die mit der Atmosphäre arbeiten und nicht nur gegen sie.

Technische Zusammenfassung

Technische Zusammenfassung: Zweiwege-Kopplung von Schwerewellen und Windpark-Wirbeln: ein reduziertes Grenzschichtmodell

Problemstellung
Die rasche Expansion von Windparks im Versorgungsmaßstab hat zu Anlagen geführt, die in ihrer Größe mit der Dicke der atmosphärischen Grenzschicht (ABL) vergleichbar sind. Dieser Skalenshift erzeugt eine komplexe Zweiwege-Wechselwirkung zwischen dem Windpark und der Atmosphäre, insbesondere unter thermischer Schichtung. In einer konventionell neutralen Grenzschicht (CNBL), die durch eine Deckungsinversion und eine geschichtete freie Atmosphäre gekennzeichnet ist, induziert die Energieentnahme durch Turbinen vertikale Bewegungen, um die Kontinuität zu erfüllen. Diese Bewegungen verdrängen die Deckungsinversion und erzeugen großskalige Schwerewellen (sowohl Grenzflächen- als auch interne Wellen). Diese Wellen induzieren wiederum Druckgradienten, die auf die Grenzschichtströmung zurückwirken, was zu einer signifikanten Aufstauung stromaufwärts und einer Veränderung der Nachlauf-Wiederherstellungsraten stromabwärts führt.

Während Large-Eddy-Simulationen (LES) diese Phänomene erfassen können, sind sie aufgrund der massiven räumlichen Domänen, die für die Ausbreitung von Schwerewellen erforderlich sind, und der Empfindlichkeit dieser Wellen gegenüber Randbedingungen rechnerisch sehr aufwendig. Bisherige reduzierte Modelle haben versucht, dies zu adressieren, leiden jedoch unter Einschränkungen: Einige gehen von vertikal einheitlichen Strömungsfeldern aus, während andere linearisierte Wellenmodelle mit technischen Nachlaufmodellen koppeln, was zu hochparametrisierten Systemen führt, die auf empirischen Entscheidungen für turbulente Diffusion, Nachlauf-Wiederherstellungsraten und Überlagerungsmethoden beruhen. Es besteht weiterhin ein Bedarf an einem einheitlichen Rahmenwerk, das die kontinuierliche vertikale Struktur der ABL auflöst und gleichzeitig die Recheneffizienz beibehält und die Parameterabhängigkeit reduziert.

Methodik
Die Autoren entwickeln ein reduziertes, linearisiertes numerisches Modell, das die kontinuierliche vertikale Struktur der ABL auflöst und die Zweiwege-Kopplung zwischen Nachlaufströmungen und Schwerewellen erfasst. Der Ansatz umfasst folgende Schritte:

1. Governing Equations (Grundgleichungen): Das Modell linearisiert die stationären, zweidimensionalen, nicht-hydrostatischen Boussinesq-Gleichungen um eine Hintergrundströmung mit horizontaler Geschwindigkeit $U(z)$. Die Gleichungen werden in zwei Bereiche unterteilt: die ABL (wo Turbulenz und Park-Antriebskräfte existieren) und die freie Atmosphäre (wo Schichtung existiert, Turbulenz und Antriebskräfte jedoch vernachlässigt werden).
2. Grenzschichtformulierung: Innerhalb der ABL leitet das Modell eine Grundgleichung für die vertikale Geschwindigkeitsstörung ($w$) her, indem es die horizontale Geschwindigkeit und den Druck eliminiert. Die turbulente Scherspannung wird unter Verwendung der Boussinesq-Hypothese mit einem Mischungsansatz modelliert.
3. Freie-Atmosphäre-Formulierung: In der freien Atmosphäre reduzieren sich die Gleichungen auf eine Helmholtz-Gleichung für die vertikale Geschwindigkeit, die mittels Fourier-Transformation gelöst wird, um Lösungen für interne und Grenzflächenwellen zu erhalten.
4. Kopplungsmechanismus: Die beiden Bereiche werden an der Deckungsinversion ($z=H$) durch kinematische und dynamische Randbedingungen gekoppelt.
   • Kinematisch: Die Inversionsschicht bewegt sich mit der Flüssigkeit, was die vertikale Geschwindigkeit mit der Inversionsverdrängung ($\eta$) verknüpft.
   • Dynamisch: Eine Druckkontinuitätsbedingung wird durch den Vergleich der Drücke vor und nach der Verdrängung der Inversion hergeleitet. Dies ergibt eine wesentliche Randbedingung (Gleichung 2.19), die den vertikalen Geschwindigkeitsgradienten an der Spitze der ABL mit der vertikalen Geschwindigkeit selbst in Beziehung setzt. Diese Bedingung fasst die Effekte sowohl von Grenzflächenwellen (über reduzierte Schwerkraft $g_r$) als auch von internen Wellen (über die Brunt-Väisälä-Frequenz $N$) zusammen.
5. Numerische Diskretisierung: Das System wird mit einer gemischten Spektral-Finite-Differenzen-Methode gelöst. Der vertikale Bereich wird mittels Finite-Differenzen diskretisiert, während die horizontale Richtung über Fourier-Transformationen behandelt wird. Dieser Ansatz ist rechnerisch effizient und passt sich natürlicherweise periodischen Domänen an.
6. Windpark-Antriebskraft: Der Park wird als halbunendliche Anordnung von Aktuatorscheiben modelliert. Der Antriebsbegriff wird iterativ berechnet, um den lokalen Nachlaufdefizit und das Verhältnis lokaler zu lateral gemittelter Nachlaufeinflüsse ($\beta_n$) zu berücksichtigen.

Hauptergebnisse
Das Modell wurde gegen drei verschiedene LES-Datensätze validiert, die unterschiedliche atmosphärische Regime (variierende Schichtung und Froude-Zahlen) repräsentieren:

• Qualitatives Verhalten: Vergleiche zwischen einer truly neutralen Grenzschicht (TNBL) und einer CNBL zeigen, dass Schwerewellen die Strömungsdynamik erheblich verändern. In der CNBL reproduziert das Modell:
  • Aufstauung stromaufwärts: Ein adverser Druckgradient entwickelt sich stromaufwärts des Parks und verursacht starke globale Blockage-Effekte, die in der TNBL nicht vorhanden sind.
  • Nachlauf-Wiederherstellung: Ein günstiger Druckgradient bildet sich innerhalb und stromabwärts des Parks, was die Nachlauf-Wiederherstellung im Vergleich zum TNBL-Fall beschleunigt.
  • Vertikale Verdrängung: Das Modell erfasst die Verdrängung der Deckungsinversion und zeigt, dass sich bei subkritischen Strömungen (starke Schichtung) die maximale Verdrängung zum Beginn des Parks verschiebt, während sie bei superkritischen Strömungen nahe dem Ende auftritt.
• Quantitative Übereinstimmung:
  • Das Modell zeigt eine gute Übereinstimmung mit LES-Daten für stromabwärts gerichtete Geschwindigkeitsstörungen in Nabenhöhe und vertikal gemittelte Geschwindigkeiten.
  • Es reproduziert erfolgreich die Größe der Aufstauung stromaufwärts und die Beschleunigung der Nachlauf-Wiederherstellung in subkritischen Fällen.
  • Im stark geschichteten subkritischen Fall (Q00) überschätzt das Modell die Nachlauf-Wiederherstellung leicht (was zu einer kleinen positiven Geschwindigkeitsstörung führt), was die Autoren auf die 2D-Vereinfachung zurückführen, die den durch Schwerewellen induzierten günstigen Druckgradienten möglicherweise überschätzt.
• Recheneffizienz: Das Modell ist rechnerisch effizient; jede Iteration dauert auf einem Standard-PC ohne Parallelisierung etwa 3 Sekunden, was es für Parametersweeps und Optimierungsaufgaben geeignet macht, bei denen LES untragbar wäre.

Bedeutung und Behauptungen
Die Arbeit behauptet, ein einheitliches Rahmenwerk zu liefern, das die Lücke zwischen rechnerisch aufwendigen LES und hochparametrisierten technischen Modellen überbrückt. Ihre Hauptbeiträge sind:

1. Physikalische Genauigkeit: Es löst die kontinuierliche vertikale Struktur der ABL auf, ohne auf die „Schichtdiskretisierung" (z. B. Aufteilung der ABL in separate Park- und obere Schichten) zurückzugreifen, die in früheren reduzierten Ansätzen verwendet wurde.
2. Reduzierte Parametrisierung: Durch die dynamische Herleitung der Kopplung aus ersten Prinzipien (Physik der Deckungsinversion) anstatt der Kopplung linearer Wellenmodelle mit empirischen Nachlaufmodellen reduziert das Rahmenwerk die Anzahl willkürlicher technischer Parameter (wie Nachlauf-Wiederherstellungsraten oder Überlagerungsmethoden).
3. Zweiwege-Kopplung: Das Modell erfasst die Rückkopplung von Schwerewellen auf die Grenzschichtströmung auf natürliche Weise durch die hergeleitete Randbedingung und reproduziert sowohl die Aufstauung stromaufwärts als auch die Beschleunigung des Nachlaufs stromabwärts genau.

Die Autoren positionieren dieses Rahmenwerk als Grundlage für zukünftige Parametrisierungen von Windpark-Effekten in Wetter- und Klimamodellen sowie als Werkzeug für die effiziente Bewertung der Windpark-Leistung und von Wechselwirkungen zwischen Windparks. Sie weisen darauf hin, dass zukünftige Arbeiten das Rahmenwerk auf drei Dimensionen erweitern und zusätzliche Physik einbeziehen werden, einschließlich nichtlinearer Effekte, Coriolis-Kräfte und verbesserter Turbulenzdarstellungen.