Numerische Untersuchung einer Düsenströmung mit Schiebewinkel

Numerisches Verfahren 18
Parameter Typ Bedeutung
my Integer Zahl der Gitterpunkte in y Richtung
nx Integer Zahl der Gitterpunkte in x Richtung
nwel Integer Punkte pro Wellenlänge in x-Richtung → ∆x = 2π/(α · nwel)
alpha Real Wellenzahl α in x-Richtung
dy Real Schrittweite ∆y in y-Richtung
Re Real Reynolszahl
Pr Real Prandtlzahl
p unend Real Anströmdruck p∞ [bar]
rb ausstrom Integer Randbedingung im Ausstromrand
Machzahl Real Machzahl
T unend Real Anströmtempertaur T∞ [K]
x0 Grenzschicht Real Grenzschicht Stelle x0 am Anfang der Paltte
Gebiete in y Integer Anzahl der Gebiete in y-Richtung
rb einstrom Integer Einstrom Randbedingung
rb freistrom Integer Freistrom Randbedingung
rb wand Integer Randbedingung an der Wand
T wand Real Temperatur an der Wand [K]
Gebiete in x Integer Anzahl der Gebiete in x-Richtung
Gittermodus Integer Art der Gittertransformation
3.1.3 Programm Beschreibung
Tab. 3.1: Bedeutung der Wichtigsten Parameter
Das modifizierte Programm Shearwake besteht aus viele Subroutine, die zur Berechnung des Nach-
laufes inkl. Scherschicht nach einer ebenen Platte dienen.
Um einen grundlegenden Überblick über die Arbeitsweise des Programms zu geben, sind in Folgen-
den die wichtigsten Schritte für die Berchnung des Nachlaufes erklärt.
Zuerst findet das Einlesen der in dem Inputfile Shearwake.i eingegebenen Parameter statt. Weiterhin
werden die Speicher für die verschiedenen Gebiete allokiert.
Danach werden Daten, die für die Berechnung in dem gesamten Programm nötig sind wie z.b.
die Schrittweite in x-Richtung, die Geschwindigkeitsverhältniss zwischen obere und untere Außen-
strömung aber auch die nötigen Daten für die Berechnung des Profils der Grenzschicht berechnet
Als nächster Schritt werden Die Grenzschichten oberhalb und unterhalb der Paltte mithilfe des Mo-
duls PROFKOM bestimmt. Hier werden die verschiedenen Grenzschichtdaten abgespeichert um sie
in den anderen Subroutinen benutzen zu können.
Mit dem Modul PROFKOM [5] wird das Grundströmungsprofil für die ebene kompressible zweidi-
mensionale Plattengrenzschicht für eine Machzahl zwischen 0 und 10 berechnet.
Danach wird das Gitter für alle Gebiete erzeugt. Der Nullpunkt ist dabei die Hinterkante der Platte.
Das Gitter hier können äquidistant oder gestreckt in y-Richtung sein. Dieses hängt von dem angegebe-
nen Gittermodus ab. Wenn der Gittermodus gleich 0 ist , ist der Abstand zwischen alle Gitterpunkte