Wie beeinflusst die asymptotische Struktur der Grenzschicht die Gitterauslegung? CFD für Pilotanlagen optimieren

Die asymptotische Struktur der Grenzschicht bestimmt ein stark nichtgleichförmiges Gitter: Verwenden Sie sehr kleine Zellen normal zu festen Wänden, wo Geschwindigkeits- und Skalargradienten steil sind, und zunehmend größere Zellen in der Außenströmung, wo die Gradienten schwächer sind. Die erforderliche wandnahe Auflösung hängt vom Löser und der Wandbehandlung ab: DNS löst die viskosen und turbulenten Skalen direkt auf, LES erfordert eine ausreichende Auflösung wandnaher Strukturen oder ein Wandmodell, und RANS zielt üblicherweise auf eine bestimmte Wandkoordinate wie (y^+) ab. In einer Pilotanlage liefern gemessene Durchflussraten, Druckverluste, Temperaturen und Wandbedingungen die erforderlichen Nachweise zur Validierung dieser Gitterstrategie.

Bei hohen Reynolds-Zahlen ist die Strömung nicht überall gleichermaßen schwer aufzulösen. Die Grenzschicht erzeugt eine dünne viskose Region, die Wandschubspannung, Druckverlust, Ablösung sowie Wärme- oder Stoffübertragung bestimmt; daher sollte die Gitterauslegung die Auflösung entsprechend den lokalen Grenzmaßstäben verteilen und sie anhand von Messungen an der Pilotanlage überprüfen.

Warum die Grenzschicht die Gitterauslegung bestimmt

Die Außenströmung benötigt keine gleichmäßig feine Auflösung

Bei hohen Reynolds-Zahlen kann die Strömung oft konzeptionell in eine annähernd reibungsfreie Außenregion und eine dünne, an die Wand angrenzende viskose Grenzschicht unterteilt werden.

Die Geschwindigkeitsgradienten normal zur Wand sind in weiten Teilen der Außenregion klein. Ein überall gleichmäßig feines Gitter verschwendet dort Rechenressourcen, ohne notwendigerweise die interessierenden Größen zu verbessern.

Die Wandregion enthält die wichtigen Gradienten

Die Haftbedingung erfordert, dass die Geschwindigkeit an der Wand null erreicht. Anschließend steigt sie innerhalb der Grenzschicht rasch auf die Geschwindigkeit der Außenströmung an.

Dadurch entstehen große Werte von (\partial u/\partial y), die direkt die Wandschubspannung bestimmen:

[ \tau_w = \mu \left.\frac{\partial u}{\partial y}\right|_w ]

Ein in wandnormaler Richtung zu grobes Gitter glättet diesen Gradienten numerisch. Die daraus entstehenden Fehler können sich als falsche Reibungsbeiwerte, Druckverluste, Widerstandskräfte, Verweilzeitverteilungen oder Wärmeübergangskoeffizienten zeigen.

Die Grenzschichtdicke nimmt mit der Reynolds-Zahl ab

Eine charakteristische Skalierung lautet

[ \frac{\delta}{L} \sim \frac{1}{\sqrt{Re}} ]

für eine laminare Grenzschicht unter geeigneten Bedingungen. Eine zunehmende Reynolds-Zahl macht die Grenzschicht daher im Allgemeinen relativ zur Anlagenlänge dünner.

Die praktische Konsequenz ist unmittelbar: Simulationen mit höheren Reynolds-Zahlen erfordern im Allgemeinen eine stärkere wandnormale Verfeinerung, sofern keine geeignete Wandfunktions- oder Wandmodellstrategie verwendet wird.

Übertragung der Struktur auf ein Gitter

Zuerst normal zur Wand verfeinern

Die erste Priorität ist der wandnormale Abstand. Über die Grenzschicht sollten mehrere Zellen verteilt sein, wobei die Platzierung nahe der Wand, wo der Geschwindigkeitsgradient am größten ist, besonders sorgfältig erfolgen muss.

Ein typisches wandnahes Gitter verwendet Inflations- oder Prismenschichten, deren Dicke mit zunehmendem Abstand von der Wand allmählich wächst. Abrupte Übergänge zwischen winzigen Wandzellen und großen Außenzellen können die Genauigkeit verringern und Konvergenzprobleme verursachen.

Das Gitter von der Wand weg strecken

Das Gitter sollte sich von der Wand zur Außenströmung hin gleichmäßig vergrößern. Dies entspricht dem abnehmenden Auflösungsbedarf, während sich das Geschwindigkeitsprofil dem Wert der Außenströmung annähert.

Das Wachstumsverhältnis muss ausreichend moderat bleiben, um die Diskretisierungsqualität zu erhalten. Übermäßige Streckung kann Interpolationsfehler, schlechte Gradientenschätzungen und eine erhöhte Empfindlichkeit des Lösers verursachen.

Die Entwicklung in Strömungsrichtung auflösen

Grenzschichten wachsen stromabwärts, sodass sich die erforderliche wandnormale Ausdehnung und Zellverteilung entlang der Anlage ändern.

Das Gitter in Strömungsrichtung muss außerdem Folgendes auflösen:

• Die Entwicklung im Einlaufbereich.
• Änderungen des Druckgradienten.
• Ablösung und Wiederanlagerung.
• Vermischung, Bögen, Verengungen und Erweiterungen.
• Lokale Änderungen der Wärme- oder Stoffübertragung.

Ein Gitter, das ausschließlich anhand der Grenzschicht am Einlass ausgelegt wurde, kann in der Nähe eines Bogens, Ventils, Rührers, Leitblechs oder einer plötzlichen Querschnittsänderung unzureichend werden.

Bei relevanten dreidimensionalen Effekten eine Spannweitenauflösung hinzufügen

Bei geraden, idealisierten Strömungen kann die Auflösung in Spannweitenrichtung weniger anspruchsvoll sein als die wandnormale Auflösung. In Pilotanlagen sind dreidimensionale Effekte jedoch häufig.

Bögen, Verteiler, Verteilerplatten, Rohrbündel, rotierende Ausrüstung und asymmetrische Formstücke können Sekundärströmungen und Quergradienten erzeugen. Diese erfordern eine ausreichende Auflösung in den Querrichtungen, statt sich ausschließlich auf eine zweidimensionale Grenzschichtannahme zu verlassen.

Das Gitter an die CFD-Methode anpassen

RANS-Simulationen

RANS-Modelle verwenden entweder wandaufgelöste Gitter oder Wandfunktionen.

Ein wandaufgelöstes RANS-Gitter platziert die erste Zelle ausreichend nahe an der Wand, um das wandnahe Geschwindigkeitsprofil aufzulösen. Ein Wandfunktionsgitter platziert die erste Zelle in einem geeigneten Bereich der logarithmischen Schicht und modelliert die wandnähere Region.

Das richtige Ziel ist nicht einfach „so fein wie möglich“. Entscheidend ist der (y^+)-Bereich, der vom gewählten Turbulenzmodell und der Wandbehandlung verlangt wird. Ein mit dieser Behandlung nicht konsistentes Gitter kann unzuverlässige Vorhersagen der Wandschubspannung und des Wärmeübergangs liefern.

LES-Simulationen

LES löst große turbulente Strukturen auf und modelliert kleinere. In der Nähe fester Wände werden die kleinsten turbulenten Skalen besonders aufwendig in der Auflösung.

Daher kann LES entweder Folgendes erfordern:

• Eine sehr feine wandnahe Auflösung für wandaufgelöste LES.
• Ein validiertes Wandmodell für Simulationen von Anlagen im Maßstab von Pilotanlagen bei hohen Reynolds-Zahlen.

Das Gitter muss außerdem in Strömungs- und Spannweitenrichtung ausreichend fein sein, um die aufgelösten turbulenten Strukturen darzustellen. Eine Verfeinerung ausschließlich in wandnormaler Richtung reicht nicht aus.

DNS-Simulationen

DNS löst alle relevanten räumlichen und zeitlichen Skalen ohne Turbulenzmodell auf. Die Anforderungen an das Gitter sind wesentlich höher, insbesondere in Wandnähe und bei hohen Reynolds-Zahlen.

DNS eignet sich im Allgemeinen für Grundlagenuntersuchungen oder kleine repräsentative Geometrien, nicht für Anlagen im Vollmaßstab einer Pilotanlage unter industriellen Betriebsbedingungen. Dennoch kann DNS Referenzdaten zur Validierung reduzierter Modelle oder von Wandbehandlungen liefern.

Messungen an der Pilotanlage zur Definition und Validierung des Gitters verwenden

Mit messbaren Betriebsbedingungen beginnen

Die Gitterauslegung sollte auf dem tatsächlichen Betriebsbereich und nicht nur auf der nominalen Geometrie basieren. Zu den wichtigen Eingangsgrößen gehören:

• Durchflussrate und Geschwindigkeit.
• Fluiddichte, Viskosität und thermische Eigenschaften.
• Turbulenz- oder Störungsniveau am Einlass.
• Wandtemperatur oder Wärmestrom.
• Oberflächenrauheit.
• Betriebsdruck und -temperatur.
• Transiente Anfahr- oder Zyklusbedingungen.

Diese Größen bestimmen lokale Reynolds-Zahlen, Grenzschichtdicken und die Regionen, in denen die Verfeinerung besonders wichtig ist.

Wandauflösungsabhängige Größen validieren

Daten aus der Pilotanlage sollten mit Simulationsergebnissen verglichen werden, die empfindlich auf die Grenzschicht reagieren, darunter:

• Druckverlust entlang von Rohrleitungen und Anlagenkomponenten.
• Reibungsbeiwert.
• Wandschubspannung.
• Geschwindigkeitsprofile.
• Positionen von Ablösung und Wiederanlagerung.
• Wärmeübergangskoeffizienten.
• Konzentrations- oder Stoffübertragungsprofile.
• Mischungs- und Verweilzeitverhalten.

Eine Übereinstimmung allein bei der globalen Durchflussrate belegt nicht, dass das wandnahe Gitter ausreichend ist.

Gitterkonvergenzuntersuchungen durchführen

Eine Gitterstudie sollte mindestens drei systematisch miteinander verknüpfte Gitter vergleichen, beispielsweise ein grobes, mittleres und feines Gitter. Der Vergleich sollte sich auf technische Ergebnisgrößen und nicht nur auf die Verringerung der Residuen konzentrieren.

Eine Simulation kann numerisch konvergieren und dennoch gitterabhängig bleiben. Wenn sich Vorhersagen von Druckverlust, Wandschubspannung oder Wärmeübergang bei einer Verfeinerung deutlich ändern, zeigt das Gitter noch keine Lösungsunabhängigkeit.

Numerische und experimentelle Unsicherheit kombinieren

Messungen an Pilotanlagen enthalten Unsicherheiten durch Sensoren, Kalibrierung, Strömungsstörungen und betriebliche Schwankungen. Numerische Ergebnisse enthalten Unsicherheiten durch Diskretisierung, Modellierung und Eingabedaten.

Die entscheidende Frage ist nicht, ob Simulation und Experiment Punkt für Punkt identisch sind. Entscheidend ist, ob die Abweichung mit der kombinierten Unsicherheit vereinbar ist und ob sich unter verschiedenen Betriebsbedingungen dieselben Trends zeigen.

Ränder, Ablösung und transiente Vorgänge berücksichtigen

Anliegende und abgelöste Strömungsbereiche unterscheiden

Grenzschichtannahmen sind für überwiegend anliegende Strömungen mit allmählichen geometrischen Änderungen am zuverlässigsten. Starke ungünstige Druckgradienten können eine Ablösung verursachen, bei der die Grenzschicht ihre Haftung an der Wand verliert und die Strömung stark instationär wird.

Ablösebereiche erfordern eine lokale Verfeinerung in mehr als einer Richtung. Nachlauf, Scherschicht, Rezirkulationsbereich und Wiederanlagerungspunkt können jeweils den Druckverlust und die Vermischung bestimmen.

Um geometrische Störungen herum verfeinern

Pilotanlagen enthalten Merkmale, die das einfache Bild einer „dünnen Schicht entlang einer ebenen Wand“ ungültig machen. Verfeinern Sie um:

• Einlässe und Auslässe.
• Bögen und T-Stücke.
• Ventile und Formstücke.
• Plötzliche Erweiterungen und Verengungen.
• Leitbleche und Einbauten.
• Rohrmündungen.
• Rührer und rotierende Teile.
• Poröse Verteiler.

Diese Bereiche erzeugen starke Gradienten, Sekundärströmungen oder Rezirkulationen, die das Ergebnis im Maßstab der Anlage dominieren können.

Relevante Zeitskalen auflösen

Bei Anfahren, Abschalten, Pulsation oder periodischem Betrieb entwickelt sich die Grenzschicht zeitlich. Der Zeitschritt muss den relevanten transienten Vorgang auflösen und mit dem räumlichen Gitter kompatibel sein.

Dimensionslose Größen wie die Strouhal-Zahl und die Reynolds-Zahl helfen zu bestimmen, ob die Prozessänderung langsam im Verhältnis zur Strömungsentwicklung ist oder schnell genug abläuft, um eine ausgeprägte instationäre Reaktion der Grenzschicht zu erzeugen.

Die Zielkonflikte verstehen

Feinere Gitter erhöhen die Kosten stark

Die wandnahe Verfeinerung erhöht die Anzahl der Zellen, den Speicherbedarf und die Rechenzeit. Dreidimensionale Inflationsschichten können diese Kosten in langen Rohrleitungen, Wärmetauschern und komplexen Prozessanlagen besonders stark erhöhen.

Das Ziel ist daher eine gezielte Auflösung und nicht die maximale Auflösung überall.

Wandfunktionen senken die Kosten, fügen jedoch Annahmen hinzu

Wandfunktionen können Simulationen bei hohen Reynolds-Zahlen praktikabel machen, indem sie die nicht aufgelöste wandnahe Region modellieren. Sie sind weniger zuverlässig, wenn die Strömung starke Ablösung, rasche Druckgradientenänderungen, Rauheitseffekte, starke Krümmung oder eine bedeutende Wärme- und Stoffübertragung aufweist.

Ihre Verwendung muss mit der Lage des Gitters und dem physikalischen Regime konsistent sein.

Grenzschichtvereinfachungen sind nicht universell gültig

Die reduzierten Grenzschichtgleichungen können bei geeigneten Strömungen mit hohen Reynolds-Zahlen, die überwiegend stromabwärts gerichtet sind, als parabolisch behandelt werden. Dies ermöglicht effiziente Marschierverfahren, da Informationen hauptsächlich stromabwärts übertragen werden.

Eine vollständige Behandlung mit den Navier-Stokes-Gleichungen ist jedoch erforderlich, wenn der Einfluss stromaufwärts, Rezirkulation, starke Dreidimensionalität, Ablösung oder eine elliptische Druckkopplung von Bedeutung sind. Eine Pilotanlage sollte nicht allein deshalb vereinfacht werden, weil ihre globale Reynolds-Zahl hoch ist.

Die formale Ordnung garantiert keine technische Genauigkeit

Eine Diskretisierung höherer Ordnung kann den Abschneidefehler reduzieren, aber sie kann keine unzureichende Wandplatzierung, schlechte Geometriedarstellung, ungeeignete Turbulenzmodellierung oder unzureichende Auflösung der Ablösung kompensieren.

Ebenso beweist das Hinzufügen von Zellen ohne Durchführung einer Konvergenzuntersuchung nicht, dass das Ergebnis vertrauenswürdig ist.

So wenden Sie dies auf Ihr Projekt an

Ein praxisnaher Arbeitsablauf besteht darin, das Gitter anhand des Betriebsbereichs und der wandempfindlichen Ergebnisgrößen auszulegen, die für den jeweiligen Prozessschritt relevant sind.

• Wenn Ihr Hauptaugenmerk auf Druckverlust oder Wandschubspannung liegt: Priorisieren Sie die wandnormale Verfeinerung, wählen Sie eine konsistente (y^+)-Strategie und validieren Sie anhand gemessener Druckverluste und Reibungsbeiwerte.
• Wenn Ihr Hauptaugenmerk auf Wärme- oder Stoffübertragung liegt: Lösen Sie sowohl die hydrodynamischen als auch die thermischen oder konzentrationsbezogenen Grenzschichten auf, die unterschiedliche Dicken aufweisen können, und vergleichen Sie die vorhergesagten Übertragungskoeffizienten mit Daten aus der Pilotanlage.
• Wenn Ihr Hauptaugenmerk auf Mischung oder Verweilzeit liegt: Verfeinern Sie Nachläufe, Rezirkulationszonen, Einlässe, Einbauten und dreidimensionale Querströmungsbereiche, statt die Auflösung nur an geraden Wänden zu konzentrieren.
• Wenn Ihr Hauptaugenmerk auf Ablösung oder instationärem Betrieb liegt: Verwenden Sie eine lokale multidirektionale Verfeinerung, ausreichend kleine Zeitschritte und Messungen, die transientes Druck-, Geschwindigkeits- oder Wiederanlagerungsverhalten erfassen.
• Wenn Ihr Hauptaugenmerk auf Recheneffizienz liegt: Verwenden Sie ein geschichtetes, lokal verfeinertes Gitter und wählen Sie RANS, LES oder DNS entsprechend der erforderlichen Physik, statt überall die aufwendigste Methode anzuwenden.
• Wenn Ihr Hauptaugenmerk auf Modellvalidierung liegt: Führen Sie Gitterkonvergenzuntersuchungen durch und vergleichen Sie mehrere Betriebsbedingungen, nicht nur einen einzigen Datenpunkt der Pilotanlage.

Ein durch die asymptotische Grenzschichtstruktur bestimmtes Gitter konzentriert den Rechenaufwand dort, wo die Physik – und die technischen Entscheidungen – am empfindlichsten sind.

Übersichtstabelle:

Planen Sie eine zuverlässige CFD-Studie für Ihre Prozessanlage? Kontaktieren Sie LABPARK, um Ihre Anforderungen zu besprechen. LABPARK bietet Ausbildungs- und berufsbildende Pilotanlagen für verfahrenstechnische Grundoperationen in den Bereichen Chemieingenieurwesen, Bioprozess- und Biotechnologie sowie Umwelt- und Wasseraufbereitung. Unsere Pilotanlagen unterstützen Universitäten, Forschungsinstitute und Unternehmen bei der Bewertung von Druckverlust, Wärme- und Stoffübertragung, Mischung, Verweilzeit, Trennung und transientem Betriebsverhalten unter messbaren Bedingungen. Wir helfen Ihnen bei der Auswahl geeigneter Anlagen, der Definition von Betriebsbedingungen, der Unterstützung wandgrenzschichtbasierter Gitterstrategien und der Verknüpfung von Simulationsergebnissen mit einer praxisnahen experimentellen Validierung.