1. Einführung
In vielen ingenieurtechnischen CFD-Simulationen, insbesondere in der Außenaerodynamik und der Windtechnik, wird die Turbulenz in Wandnähe üblicherweise mithilfe von Standard-Wandfunktionen modelliert, anstatt die viskose Unterschicht direkt aufzulösen. Ein typisches Beispiel ist die in OpenFOAM verwendete (und in RWIND angewandte) Funktion nutkWallFunction, die die turbulente Viskosität an der Wand auf der Grundlage des logarithmischen Wandgesetzes und empirischer Annahmen zur Turbulenzmodellierung schätzt. Diese Modellierungsstrategie ermöglicht die Durchführung von Simulationen mit relativ groben Netzen, was zu einer erheblichen Reduzierung des Rechenaufwands und der Simulationszeit führt und dadurch für groß angelegte ingenieurtechnische Probleme wie Gebäude, Brücken und städtische Windsimulationen praktisch ist.
Die durch Standard-Wandfunktionen eingeführten Vereinfachungen bringen jedoch auch erhebliche Einschränkungen hinsichtlich der Genauigkeit der Strömungsvorhersage in Wandnähe mit sich. Diese Einschränkungen können Größen wie die Wandschubspannung, Grenzschichtentwicklung und lokale Druckverteilungen beeinflussen und müssen daher bei der Interpretation aerodynamischer Ergebnisse, der Beurteilung der Modellgenauigkeit oder der Durchführung von Validierungsstudien anhand experimenteller Daten sorgfältig berücksichtigt werden.
Dieser Artikel behandelt die wichtigsten Konsequenzen und technischen Auswirkungen der Verwendung von Standardansätzen mit Wandfunktionen in aerodynamischen CFD-Simulationen unter der Annahme sehr hoher y⁺-Werte, die außerhalb des empfohlenen oder akzeptierten Bereichs liegen. Außerdem werden Situationen aufgezeigt, in denen diese Modellierungsstrategie geeignet ist, sowie Fälle, in denen ihre Einschränkungen die Zuverlässigkeit der Ergebnisse beeinträchtigen können.
2. Hintergrund: Wandfunktionen in der aerodynamischen CFD
In der Nähe fester Wände weisen turbulente Grenzschichten mehrere Bereiche auf:
- Viskose Unterschicht (niedriges 𝑦+)
- Pufferschicht
- Logarithmischer Bereich (hohes 𝑦+)
Standard-Wandfunktionen gehen von folgenden Annahmen aus:
- Die erste Berechnungszelle liegt im logarithmischen Bereich (typischerweise 𝑦+>30).
- Das Geschwindigkeitsprofil in Wandnähe folgt dem logarithmischen Gesetz der Wand.
- Die turbulente Viskosität und die Schubspannung werden anhand empirischer Beziehungen geschätzt, anstatt kleine Strömungsstrukturen aufzulösen.
Diese Modellierungsstrategie findet breite Anwendung in:
- der Außenaerodynamik
- Simulationen im Bereich der Windtechnik
- groß angelegten Strömungssimulationen in der Umwelt
da sie stabile Lösungen mit relativ groben Netzen ermöglicht.
3. Vorteile des derzeitigen Wandfunktionsansatzes mit hohem y⁺
Die Verwendung von Standard-Wandfunktionen mit hohen 𝑦+-Werten bietet mehrere praktische Vorteile für ingenieurtechnische CFD-Simulationen. Da die viskose Unterschicht nicht aufgelöst werden muss, ermöglicht dieser Ansatz die Verwendung gröberer Netze, was zu geringeren Rechenaufwänden und schnelleren Simulationen führt, während das für viele groß angelegte aerodynamische Studien relevante allgemeine Strömungsverhalten weiterhin erfasst wird.
- Hohe Recheneffizienz
Ermöglicht die Verwendung relativ grober Netze in Wandnähe, wodurch die Anzahl der Rechenzellen und die Simulationszeit erheblich reduziert werden.
- Geringerer Rechenaufwand
Dadurch wird das extrem feine Netz vermieden, das zur Darstellung der viskosen Unterschicht (𝑦+≈1) erforderlich wäre und den CPU- und Speicherbedarf drastisch erhöhen würde.
- Geeignet für große Rechenbereiche
Effizient für Simulationen in großflächigen Umgebungen wie Stadtgebieten, Geländemodellen und komplexen Gebäudeclustern.
- Numerische Robustheit und Stabilität
Standard-Wandfunktionsansätze sind im Allgemeinen weniger empfindlich gegenüber Netzunregelmäßigkeiten und liefern in der Regel stabile Lösungen für komplexe Geometrien.
- Effizient für parametrische Studien
Ermöglicht einen schnellen Vergleich mehrerer Bemessungsalternativen, Gebäudeausrichtungen oder Umgebungsbedingungen.
- Gut geeignet für Bemessungsstudien im Frühstadium
Liefert schnelle aerodynamische Erkenntnisse in der Konzeptionsphase, wenn keine detaillierte Auflösung der Grenzschicht erforderlich ist.
- Effektiv zur Erfassung des globalen Strömungsverhaltens
Stellt großräumige Strömungsmuster wie Nachlaufgebiete, Strömungsbeschleunigungszonen und die allgemeine Windverteilung genau dar.
- Praktischer technischer Kompromiss
Schafft einen Ausgleich zwischen Recheneffizienz und akzeptabler technischer Genauigkeit, wenn der Schwerpunkt eher auf dem allgemeinen aerodynamischen Verhalten als auf detaillierten physikalischen Vorgängen in Wandnähe liegt.
4. Hauptnachteile in aerodynamischen Simulationen
4.1 Reduzierte Auflösung in Wandnähe
Standard-Wandfunktionen lösen die viskose Unterschicht nicht auf. Folglich:
- werden detaillierte Geschwindigkeitsgradienten in Wandnähe nicht erfasst
- werden Turbulenzstrukturen in Wandnähe nur approximiert statt berechnet
- kann die lokale Strömungsphysik zu stark vereinfacht werden
Diese Einschränkung wirkt sich direkt auf die aerodynamische Genauigkeit aus, wenn wandgebundene Effekte dominieren.
4.2 Ungenaue Wandschubspannung und Mantelreibung
Wandfunktionen schätzen die Wandschubspannung anhand empirischer logarithmischer Gesetzmäßigkeiten. Dies kann folgende Auswirkungen haben:
- Fehler beim Mantelreibungswiderstand
- Geringere Genauigkeit bei der Vorhersage des Gesamtwiderstands für stromlinienförmige Körper
- Ungenaue Schubspannungsverteilung entlang der Oberflächen
Für die Analyse der aerodynamischen Leistung können diese Abweichungen erheblich sein
4.3 Grenzen der Vorhersage von Strömungsablösungen
Wandfunktionen gehen von einem Gleichgewichtsverhalten der Grenzschicht aus, das in folgenden Fällen nicht mehr zutrifft:
- Starke negative Druckgradienten
- Schnelle Strömungsverlangsamung
- Komplexe geometriebedingte Ablösung
Zu den typischen Folgen zählen:
- Falsche Vorhersage des Ablösungsbeginns
- Fehler bei der Wiederanlegestelle
- Ungenaue Nachlaufstruktur
Da die Ablösung einen starken Einfluss auf aerodynamische Kräfte hat, stellt dies einen erheblichen Nachteil dar
4.4 Reduzierte Genauigkeit in Nachlauf- und Rezirkulationsbereichen
Da die Turbulenz in Wandnähe nicht vollständig aufgelöst wird:
- kann die Entwicklung der Scherschicht ungenau modelliert sein
- kann die Größe der Rezirkulationszonen falsch berechnet sein
- können Wirbelstrukturen übermäßig diffus dargestellt sein
Dies wirkt sich auf die Abschätzung der aerodynamischen Last und die Strömungsvisualisierung aus
4.5 Empfindlichkeit gegenüber y⁺-Platzierung
Die Gültigkeit der Wandfunktion hängt von der richtigen Platzierung der ersten Netzzelle in Wandnähe ab:
- Hohes 𝑦+ → Logarithmusgesetz gilt
- Niedriges 𝑦+ → Annahmen treffen nicht mehr zu
- Mittleres 𝑦+ (Pufferzone) → Modellierungsunsicherheit nimmt zu
Eine schlechte Netzgestaltung kann daher zu versteckten Fehlern führen, selbst wenn Simulationen numerisch stabil erscheinen.
4.6 Eingeschränkte Leistungsfähigkeit bei Übergangsströmungen
Standard-Wandfunktionen gehen von vollständig turbulenten Strömungsbedingungen aus. Das hat folgende Konsequenzen:
- Der Übergang von laminarer zu turbulenter Strömung kann nicht genau erfasst werden.
- Die Entwicklung der Grenzschicht kann falsch dargestellt werden.
Dies ist wichtig für aerodynamische Anwendungen, bei denen Folgendes eine Rolle spielt:
- Geringe Turbulenzintensität
- Glatte Oberflächen
- Strömungen an Flügelprofilen
4.7 Reduzierte Genauigkeit bei lokalen Druckspitzen
Die empirische Wandmodellierung kann Druckgradienten in Oberflächennähe glätten, was zu Folgendem führen kann:
- Unter- oder Überschätzung der Spitzendruckbeiwerte
- Verminderte Genauigkeit bei der lokalen Lastbewertung
- Mögliche Fehler bei statischen Bemessungslasten
4.8 Abhängigkeit von Annahmen bezüglich des empirischen Logarithmusgesetzes
Wandfunktionen basieren auf vereinfachten empirischen Modellen, die Folgendes annehmen:
- Vollentwickelte turbulente Grenzschichten
- Glatte Strömungsentwicklung
- Moderate Druckgradienten
Wenn die tatsächlichen aerodynamischen Bedingungen von diesen Annahmen abweichen, nimmt die Genauigkeit ab
5. Technischer Kompromiss: Genauigkeit vs. Effizienz
Trotz der Nachteile werden Standardwandfunktionen in der technischen CFD nach wie vor häufig verwendet, da sie mehrere wichtige praktische Vorteile bieten:
- Verbesserte numerische Stabilität
- Geringere Netzanforderungen
- Geringerer Rechenaufwand
- Bessere Zugänglichkeit für technische Anwendungen, wodurch Simulationen innerhalb realistischer Zeit- und Ressourcengrenzen durchgeführt werden können
In praktischen technischen Workflows sind extrem hochauflösende Simulationen, die die viskose Unterschicht vollständig auflösen, oft sehr rechenintensiv und für routinemäßige Bemessungsaufgaben nicht immer durchführbar. Daher stellen Wandfunktionsansätze nicht nur einen Kompromiss zwischen Genauigkeit und Aufwand dar, sondern auch eine praktische Modellierungsstrategie, die aerodynamische Simulationen für alltägliche technische Anwendungen zugänglich macht. Wie auch in dem Artikel hervorgehoben wird, sind hochdetaillierte und rechenintensive Simulationen in typischen technischen Projekten, bei denen Zeit und Rechenressourcen begrenzt sind, oft noch nicht praktikabel.
Aus diesen Gründen sind Standard-Wandfunktionsansätze oft akzeptabel für:
- Strömungsvisualisierung
- Bemessungsstudien in frühen Phasen
- Konzeptionelle aerodynamische Bewertung
- Topografische oder geländebedingte Windstudien
- Globale Lastabschätzung
Sie sind jedoch möglicherweise nicht für Anwendungen geeignet, bei denen eine detaillierte Strömungsphysik in Wandnähe entscheidend ist, wie zum Beispiel:
- Spitzendruckberechnungen
- Lokale Fassadendruckbemessung
- Instationäre Simulationen mit detaillierten Strömungsschwankungen
- Lokale Effekte in der Nähe von scharfen Kanten oder Ecken
- Präzise Analysen der Ablösung und des Wiederanlegens
Tabelle 1 klassifiziert aerodynamische CFD-Anwendungen nach ihrer Eignung bei Verwendung eines Standard-Wandfunktionsansatzes mit hohem 𝑦+-Wert, bei dem die viskose Unterschicht nicht aufgelöst wird. Anwendungen, die sich auf das globale Strömungsverhalten konzentrieren, werden empfohlen; einige Analysen, die eine moderate Genauigkeit erfordern, werden unter Vorbehalt empfohlen, während Untersuchungen, die von detaillierten physikalischen Phänomenen in Wandnähe abhängen, nicht empfohlen werden.
Tabelle 1: Anwendbarkeit von aerodynamischen RWIND-CFD-Simulationen unter der Annahme einer Standard-Wandfunktion für hohe y⁺-Werte
| Nr. | Anwendung | Hauptziel | Empfehlung | Technische Rechtfertigung |
|---|---|---|---|---|
| 1 | Globale Windlasten auf Gebäude | Gesamtkräfte und -momente | 🟢 Empfohlen | Von großräumiger Druckverteilung dominiert; detaillierte Wandschubauflösung weniger kritisch |
| 2 | Fassadenbemessung | Verkleidungsdruckzonierung | 🔴 Nicht empfohlen | Lokale Spitzendrücke in der Nähe von Kanten und Ecken werden möglicherweise zu niedrig eingeschätzt |
| 3 | Windkomfort für Fußgänger | Windgeschwindigkeit in 1,5–2 m Höhe | 🟢 Empfohlen | Schwerpunkt auf dem Geschwindigkeitsfeld; Schubspannung in Wandnähe hat nur begrenzten Einfluss |
| 4 | Untersuchungen zur Windströmung im städtischen Raum | Windverteilung in Stadtblöcken | 🟢 Empfohlen | Effizienz im großen Maßstab hat Vorrang vor detaillierter Auflösung der Wandschicht |
| 5 | Strömungsvisualisierung | Strömungslinien und Nachlaufstrukturen | 🟢 Empfohlen | Qualitative Trends bleiben trotz Vereinfachung in Wandnähe robust |
| 6 | Vergleich parametrischer Bemessungen | Entwicklungen der relativen Leistung | 🟢 Empfohlen | Einheitliche Modellannahmen ermöglichen eine zuverlässige vergleichende Bewertung |
| 7 | Dachauftrieb (mittlere Saugkraft) | Durchschnittliche Dachdrücke | 🟡 Bedingt empfohlen | Mittlere Saugkraft akzeptabel; lokale Eckwerte stark netzabhängig |
| 8 | Mittlerer Luftwiderstand und Auftrieb des Staukörpers | Globale aerodynamische Koeffizienten | 🟡 Bedingt empfohlen | Mittlere Kräfte werden angemessen erfasst; Ablösungsstelle ist weniger präzise |
| 9 | Detaillierte Ablösungsstudie | Genaue Ablösungs- und Wiederanlegestelle | 🔴 Nicht empfohlen | Die Annahme einer Gleichgewichts-Wandfunktion ist bei starken negativen Druckgradienten ungültig |
| 10 | Analyse des Mantelreibungswiderstands | Aufschlüsselung des Reibungswiderstands | 🔴 Nicht empfohlen | Ein hoher y+-Wert verhindert eine genaue Auflösung der Wandschubspannung und des Grenzschichtgradienten |
| 11 | Tragflächen- oder Fahrzeugaerodynamik | Auftriebs- und Widerstandsverhalten | 🔴 Nicht empfohlen | Erfordert einen y+-Wert von etwa 1 und eine aufgelöste viskose Unterschicht |
| 12 | Lokaler Spitzendruck | Cp-Spitzen an scharfen Kanten | 🔴 Nicht empfohlen | Ein hoher y+-Wert glättet Druckgradienten und unterschätzt lokale Extrema |
| 13 | Instationäre Wirbelablösung | Dominante Ablösungsfrequenz | 🔴 Nicht empfohlen | Fluktuationsamplitude weniger zuverlässig |
| 14 | Windstudie erhöhter Terrasse | Lokale Beschleunigungszonen | 🟡 Bedingt empfohlen | Allgemeine Trends werden erfasst; Randeffekte ungewiss |
| 15 | Validierung auf Forschungsniveau | Genauigkeit auf Publikationsniveau | 🔴 Nicht empfohlen | Eine hochgenaue Validierung erfordert eine wandaufgelöste oder verbesserte Modellierung in Wandnähe |
6. Fazit
Standard-Wandfunktionsansätze wie nutkWallFunction bieten eine praktische und effiziente Methode zur Modellierung turbulenter Grenzschichten in aerodynamischen Simulationen. Die Verwendung empirischer Näherungen nach dem Logarithmusgesetz bringt jedoch erhebliche Einschränkungen mit sich, darunter eine geringere Genauigkeit bei der Physik in Wandnähe, eine weniger zuverlässige Vorhersage der Ablösung sowie eine hohe Empfindlichkeit gegenüber der Netzgestaltung. Ingenieure sollten sich dieser Nachteile bewusst sein, wenn sie Simulationsergebnisse interpretieren, CFD-Validierungen durchführen oder Simulationen mit experimentellen Daten vergleichen.