Numerische Strömungsmechanik

Visualisierung einer CFD-Simulation der Boeing X-43 bei Mach 7

Die numerische Strömungsmechanik (englisch: computational fluid dynamics, CFD) ist eine etablierte Methode der Strömungsmechanik. Sie hat das Ziel, strömungsmechanische Probleme approximativ mit numerischen Methoden zu lösen. Die benutzten Modellgleichungen sind meist Navier-Stokes-Gleichungen, Euler-Gleichungen oder Potentialgleichungen. Die Motivation hierzu ist, dass wichtige Probleme wie zum Beispiel die Berechnung des Widerstandsbeiwerts sehr schnell zu nichtlinearen Problemen führen, die nur in Spezialfällen exakt lösbar sind. Die numerische Strömungsmechanik bietet dann eine kostengünstige Alternative zu Versuchen im Windkanal oder Wasserkanal.

Die international gebräuchliche Abkürzung CFD wird etwa seit einer Konferenz der American Institute of Aeronautics and Astronautics (AIAA) 1973 benutzt. Dort wurde auch die Verwendung von CFD als Werkzeug zum Design von Flugzeugen etabliert.

Modelle

Das umfassendste Modell sind die Navier-Stokes-Gleichungen. Es handelt sich hierbei um ein System von nichtlinearen partiellen Differentialgleichungen 2. Ordnung, die die meisten Fluide komplett beschreiben. Insbesondere sind auch Turbulenz und die hydrodynamische Grenzschicht enthalten, was allerdings zu höchsten Ansprüchen an Rechnerleistung, Speicher und die numerischen Verfahren führt.

Ein einfacheres Modell sind die Euler-Gleichungen, die aufgrund der vernachlässigten Reibung die Grenzschicht nicht abbilden und auch keine Turbulenz enthalten, womit beispielsweise Strömungsabriss nicht über dieses Modell simuliert werden kann. Dafür sind wesentlich gröbere Gitter geeignet, um die Gleichungen sinnvoll zu lösen. Für diejenigen Teile der Strömung, in denen die Grenzschicht keine wesentliche Rolle spielt, sind die Euler-Gleichungen sehr gut geeignet.

Die Potentialgleichungen schließlich sind vor allem nützlich, wenn schnell grobe Vorhersagen gemacht werden sollen. Bei ihnen wird die Entropie als konstant vorausgesetzt, was bedeutet, dass keine starken Schockwellen auftreten können, da an diesen die Entropie sogar unstetig ist. Weitere Vereinfachung über konstante Dichte führt dann zur Laplace-Gleichung.

Bei Mehrphasenströmungen spielen Wechselwirkungskräfte zwischen den Phasen eine Rolle, wobei geeignete Vereinfachungen durchgeführt werden können.

CFD-Verfahren bilden auch die Grundlage für die numerische Aeroakustik, die sich mit der Berechnung von Strömungsgeräuschen befasst.

Verfahren

Die verbreitetsten Lösungsmethoden der numerischen Strömungsmechanik sind

Die FEM ist für viele Probleme geeignet, insbesondere für elliptische und parabolische im inkompressiblen Bereich, weniger für hyperbolische. Sie zeichnet sich durch Robustheit und solide mathematische Untermauerung aus. FVM ist für Erhaltungsgleichungen geeignet, insbesondere für kompressible Strömungen. FDM ist sehr einfach und deswegen vor allem von theoretischem Interesse.

Weitere sind

die Spektralmethode
die Lattice-Boltzmann-Methode (LBM)
die Smoothed Particle Hydrodynamics (SPH)
die Randelementmethode (boundary element method, BEM)
die Finite Pointset Methode (FPM)

Bei allen Methoden handelt es sich um numerische Näherungsverfahren, die zur Validierung mit quantitativen Experimenten verglichen werden müssen. Mit Ausnahme der partikelbasierten Methoden ist der Ausgangspunkt der oben genannten Methoden die Diskretisierung des Problems mit einem Rechengitter.

Zeitabhängige Gleichungen

Bei zeitabhängigen Gleichungen führt die Reihenfolge von Orts- und Zeitdiskretisierung auf zwei verschiedene Lösungsansätze:

Vertikale Linienmethode: Hier wird zunächst im Ort diskretisiert, sodass man ein System von gewöhnlichen Differentialgleichungen in der Zeit erhält.
Horizontale Linienmethode (oder Rothe-Methode): Die Zeitdiskretisierung erfolgt zuerst, und die Gleichungen reduzieren sich auf die Lösung eines Randwertproblems in jedem Zeitschritt.

Die erste Methode wird vor allem bei hyperbolischen Gleichungen und kompressiblen Strömungen, letztere bei inkompressiblen Strömungen eingesetzt. Außerdem ist die Rothe-Methode flexibler im Hinblick auf eine Implementierung einer adaptiven Gitterverfeinerung im Ort während der Zeitevolution der Strömungsgleichungen.

Turbulente Strömungen

Bei turbulenten Strömungen gibt es für die numerische Strömungssimulation noch viele offene Fragen: Entweder man verwendet sehr feine Rechengitter wie bei der direkten numerischen Simulation oder man verwendet mehr oder weniger empirische Turbulenzmodelle, bei denen neben numerischen Fehlern zusätzliche Modellierungsfehler auftreten. Einfache Probleme können auf Highend-PCs in Minuten gelöst werden, während komplexe 3D-Probleme selbst auf Großrechnern teilweise kaum zu lösen sind.

Basierend auf Artikeln in:

Wikipedia.de