Lehrstuhl Verbrennungskraftmaschinen und Flugantriebe ...
Lehrstuhl Verbrennungskraftmaschinen und Flugantriebe ... Lehrstuhl Verbrennungskraftmaschinen und Flugantriebe ...
Das k-ε Modell wird auch in speziellen Modifikationen für den Niedrig-Reynoldszahlbereich verwendet, wenn die Gl. 2.62 mit der Dämpfungsfunktion f µ ergänzt wird. 2 = C ⋅ f ⋅ ρ ⋅ k ~ ε Gl. 2.66 µ t µ µ Des weiteren werden in Abhängigkeit vom modifizierten Turbulenzmodell zusätzliche Funktionen in die reynoldsgemittelten Navier-Stokes-Gleichungen eingefügt. Die Beschreibung der Funktionen entnehmen wir NUMECA [30]. Eine besondere Form der Turbulenzmodelle sind Zweischichtenmodelle, die auf modifizierten Versionen des Eingleichungsmodells bzw. der k-ε Modelle beruhen. Bei einer Art der Zweischichtenmodelle wird in der viskosen Unterschicht auf ein einfacheres Eingleichungsmodell umgeschaltet, das keine ε-Gleichung löst, sondern die Längenmaßverteilung empirisch vorgibt, Rodi [33]. Das Eingleichungsmodell wird nur in direkter Wandnähe eingesetzt, während die wandfernen Strömungsgebiete mit dem Standardk-ε Modell berechnet werden. Das Umschalten vom Eingleichungsmodell in Wandnähe auf das Standard k-ε Modell im äußeren Bereich erfolgt an einer Stelle, bei der die Dämpfungsfunktion f µ den Wert 0,95 erreicht, d.h. wo die viskosen Effekte vernachlässigt werden können. Eine weitere Art der Zweischichtmodelle ist die Verwendung eines modifizierten Niedrig- Reynoldszahl-k-ε Modells für die viskose Unterschicht und des Standard k-ε Modells bei einem normierten Wandabstands y + . In FINE/TURBO V3.0 ist die Verwendung der letztbeschriebenen Art des Zweischichtenmodells unter Verwendung des k-ε Modells der Version Yang und Shi im Bereich der viskosen Unterschicht und das Standard k-ε Modell im äußeren Bereich ab einem normierten Wandabstand von y + (f µ =0,95)=260,5 möglich. 2.9.3 Turbulenzmodelle von FINE/TURBO V3.0 Folgende Turbulenzmodelle sind im Euranus/Turbo-Code von FINE/TURBO V3.0 implementiert. • Baldwin-Lomax-Modell • k-ε Modell Standard • k-ε Modell Chien • k-ε Modell Launder & Sharma • k-ε Modell Yang & Shi • Zweigleichungsmodell k-ε Modell Yang & Shi Das Beeinflussen der spezifischen Turbulenzparameter ist ebenfalls möglich. 30
2.10 Anfangs- und Randbedingungen Zur Lösung der Erhaltungsgleichungen ist die Vorgabe von bestimmten Anfangs- und Randbedingungen notwendig. Dies umfaßt sämtliche konservative Größen. 2.10.1 Anfangsbedingungen Der zeitgenauen oder instationären Berechnung geht eine Linearisierung der Erhaltungsgleichungen zum Zeitpunkt t = 0 voraus, um den exakten Lösungsvektor der Erhaltungsgleichungen im betrachtenen Rechenraum zu erhalten. Eine triviale Lösung wäre die Null-Strömung-Annahme, die jedoch aufgrund ihrer schwierigen numerischen Beherrschbarkeit nicht angewendet wird, Null-Geschwindigkeitskomponenten, konstanter Druck und Temperaturverteilung. Die Null-Strömung-Annahme erfordert einen hohen numerischen Aufwand, denn je näher die Anfangsbedingung letztendlich der endgültigen Lösung ist, ums so schneller erfolgt die numerische Approximation. Stationäre Strömungen werden allein durch die Randbedingungen bestimmt. Passen wir die Anfangsbedingungen dem erwarteten Strömungsverlauf an, so erhalten wir eine stabilere und schnellere Konvergenz der Rechnungen. Somit ist die Angabe einer approximierten Anfangsbedingung immer sinnvoll. Des weiteren bieten sich Rechnungen im Rahmen des full multi grid modus auf einem groben Netz an, um deren Ergebnisse dann als Startbedingung in ein feineres Netz interpolieren zu können. Diese Funktion ist in FINE/TURBO V3.0 jedoch nur für der Lösung der Euler-Gleichungen sowie für die turbulenten Navier-Stokes-Gleichungen nur bei Verwendung des Baldwin- Lomax-Modells vorgesehen. 2.10.2 Randbedingungen An den Rändern des Rechenraums müssen sämtliche konservative Größen bekannt sein. Physikalisch lassen sich die Randbedingungen durch Formulierungen von Dirichlet und Neumann beschreiben. Die Dirichlet-Randbedingung legt die Strömungsgröße am Rand fest. Die Neumann-Randbedingug definiert den Fluß dieser Größe, d.h. die zeitliche Abbleitung der Strömungsgröße senkrecht zur Berandung. Wir unterscheiden Ein- und Ausströmrand, Festkörperwand sowie Festlegungen für innere Randbedingungen, wie periodische Randbedingungen, Symmetrierandbedingungen, Verbindungsrandbedingungen sowie externe Randbedingungen für die freie Umströmung. 31
- Seite 1 und 2: Brandenburgische Technische Univers
- Seite 3 und 4: II. INHALTSVERSZEICHNIS I. EIDESSTA
- Seite 5 und 6: III. VERWENDETE FORMELZEICHEN LATEI
- Seite 7 und 8: µ dynamische Viskosität -1 Kg·m
- Seite 9 und 10: Mach- oder Reynoldszahl. Die Meßte
- Seite 11 und 12: konvertierten CAD-Dateien verträgt
- Seite 13 und 14: ∫∫ ρ v ⋅ dA = ∫∫∫ ∇
- Seite 15 und 16: Einen Zusammenhang zwischen der Vol
- Seite 17 und 18: sind eindeutig und voneinander unab
- Seite 19 und 20: Wir unterscheiden zwei wichtige Dis
- Seite 21 und 22: einen bestimmten kritischen Betrag
- Seite 23 und 24: Für stationäre Strömungen sowie
- Seite 25 und 26: Die Turbulenz einer Strömung kann
- Seite 27 und 28: F KLEB ( y) 6 −1 ⎡ ⎛ CKLEB
- Seite 29: Der Dissipationsterm ist ~ ∂ v r
- Seite 33 und 34: 2.10.2.3 Innere Randbedingungen Inn
- Seite 35 und 36: Turbulenzmodellen bessere Ergebniss
- Seite 37 und 38: 2.13 Residuum Mit dem Residuum best
- Seite 39 und 40: Referenzgrößen Referenzlänge l 1
- Seite 41 und 42: Abb. 3.3 - Stromlinien in der Meßk
- Seite 43 und 44: Abb. 3.7 - Isolinien der Machzahl M
- Seite 45 und 46: Neben den einfachen zweidimensional
- Seite 47 und 48: 2 v c p ⋅ T0 = c p ⋅ T + Gl. 3.
- Seite 49 und 50: unterschiedliche Netze für die Ber
- Seite 51 und 52: In der Genzschicht werden stark ver
- Seite 53 und 54: Abb. 3.16 Geschwindigkeitsfeld in d
- Seite 55 und 56: In den Abb. 3.20 bis Abb. 3.22 sind
- Seite 57 und 58: Abb. 3.22 - Isolinien des statische
- Seite 59 und 60: In den Abb. 3.26 bis Abb. 3.28 sind
- Seite 61 und 62: Abb. 3.29 - Isolinien der statische
- Seite 63 und 64: Abb. 3.33 - Isolinien der Dichte ρ
- Seite 65 und 66: Abb. 3.37 - Isolinien der kinetisch
- Seite 67 und 68: Unser Ziel ist es, in Abhängigkeit
- Seite 69 und 70: Den Anstellwinkel α stellen wir ü
- Seite 71 und 72: Abb. 3.43 - Verteilung des statisch
- Seite 73 und 74: Abb. 3.47 - Geschwindigkeitsfeld um
- Seite 75 und 76: Abb. 3.51 - Geschwindigkeitsfeld um
- Seite 77 und 78: Abb. 3.55 - Geschwindigkeitsfeld um
- Seite 79 und 80: Abb. 3.59 - Geschwindigkeitsfeld um
Das k-ε Modell wird auch in speziellen Modifikationen für den Niedrig-Reynoldszahlbereich<br />
verwendet, wenn die Gl. 2.62 mit der Dämpfungsfunktion f µ ergänzt wird.<br />
2<br />
= C ⋅ f ⋅ ρ ⋅ k ~ ε<br />
Gl. 2.66<br />
µ t µ µ<br />
Des weiteren werden in Abhängigkeit vom modifizierten Turbulenzmodell zusätzliche<br />
Funktionen in die reynoldsgemittelten Navier-Stokes-Gleichungen eingefügt. Die<br />
Beschreibung der Funktionen entnehmen wir NUMECA [30].<br />
Eine besondere Form der Turbulenzmodelle sind Zweischichtenmodelle, die auf modifizierten<br />
Versionen des Eingleichungsmodells bzw. der k-ε Modelle beruhen.<br />
Bei einer Art der Zweischichtenmodelle wird in der viskosen Unterschicht auf ein einfacheres<br />
Eingleichungsmodell umgeschaltet, das keine ε-Gleichung löst, sondern die<br />
Längenmaßverteilung empirisch vorgibt, Rodi [33]. Das Eingleichungsmodell wird nur in<br />
direkter Wandnähe eingesetzt, während die wandfernen Strömungsgebiete mit dem Standardk-ε<br />
Modell berechnet werden. Das Umschalten vom Eingleichungsmodell in Wandnähe auf<br />
das Standard k-ε Modell im äußeren Bereich erfolgt an einer Stelle, bei der die<br />
Dämpfungsfunktion f µ den Wert 0,95 erreicht, d.h. wo die viskosen Effekte vernachlässigt<br />
werden können.<br />
Eine weitere Art der Zweischichtmodelle ist die Verwendung eines modifizierten Niedrig-<br />
Reynoldszahl-k-ε Modells für die viskose Unterschicht <strong>und</strong> des Standard k-ε Modells bei<br />
einem normierten Wandabstands y + . In FINE/TURBO V3.0 ist die Verwendung der<br />
letztbeschriebenen Art des Zweischichtenmodells unter Verwendung des k-ε Modells der<br />
Version Yang <strong>und</strong> Shi im Bereich der viskosen Unterschicht <strong>und</strong> das Standard k-ε Modell im<br />
äußeren Bereich ab einem normierten Wandabstand von y + (f µ =0,95)=260,5 möglich.<br />
2.9.3 Turbulenzmodelle von FINE/TURBO V3.0<br />
Folgende Turbulenzmodelle sind im Euranus/Turbo-Code von FINE/TURBO V3.0<br />
implementiert.<br />
• Baldwin-Lomax-Modell<br />
• k-ε Modell Standard<br />
• k-ε Modell Chien<br />
• k-ε Modell La<strong>und</strong>er & Sharma<br />
• k-ε Modell Yang & Shi<br />
• Zweigleichungsmodell k-ε Modell Yang & Shi<br />
Das Beeinflussen der spezifischen Turbulenzparameter ist ebenfalls möglich.<br />
30