Numerische Optimierung dreidimensional parametrisierter ...
Numerische Optimierung dreidimensional parametrisierter ... Numerische Optimierung dreidimensional parametrisierter ...
Impuls aus der Seitenwandsenke füllt die Ablöseblase auf dem Profil und wird ohne Rotationaufzubauen re-energetisiert.Durch die dreidimensionale Gestaltung von Profil und Seitenwand konnten die ursächlichenanfachenden Mechanismen der Sekundärströmungen wirkungsvoll reduziert werden. Des weiterenwird ein Aufrollen der Grenzschichten zu den bekannten Wirbelphänomenen vermieden.Die Sekundärströmungsphänomene bilden sich erst weit hinten aus. Es zeigt sich nicht dasübliche Verlustbild in der Profilabströmung. Die Sekundärphänomene bleiben trotz der divergentenSeitenwände in Seitenwandnähe.Die Güte der Optimierung läßt sich anhand der integralen Beurteilungsparameter überprüfen.Die wichtigsten integralen Größen sind in der folgenden Übersicht gegenübergestellt:T106D T106Dopt T106Doptnumerisch numerisch experimentellZuströmwinkel 127.7° 127.7° 127.7°Bezugs-Mach-Zahl 0.59 0.59 0.59Bezugs-Reynolds-Zahl 500000 500000 498735Abströmwinkel β 2 grobes Netz 25.15° 24.90° 24.76°Optimierung flußgemittelt alsRandbedingungAbströmwinkel β 2 feines NetzNachrechnung / massengemitteltExperiment24.84° 24.57° 24.76°Totaldruckverlustbeiwert feineNachrechnung / ExperimentmassengemitteltTotaldruckverlustbeiwert flußgemitteltfeine Nachrechnung0.0658 0.0591 0.05730.0881 0.0719Die integralen flußgemittelten Verluste konnten durch die Optimierung um 22 % relativ zurStartgeometrie reduziert werden. Das Auslegungsziel einer Reduzierung der Sekundärströmungen,bei gleichzeitiger integraler Verlustminimierung, konnte erreicht werden.9. Diskussion und AusblickDurch die Automatisierung der Auslegungskette wird die aerodynamische Berechnung mit derkonstruktiven Entwurfsebene von Profilen gekoppelt. Die Anwendung einer solchen Auslegungsmethodikauf eine vollständig parametrisierte Beschaufelungsgeometrie ist zum heutigen102
Zeitpunkt sicherlich noch aufwendig und neu. Deswegen sollen hier noch einige abschließendeAnmerkungen gegeben werden.Um mit einem Optimierungsalgorithmus bei der Auslegung von dreidimensionalen Beschaufelungen,nach Möglichkeit global, optimale Lösungen erreichen zu können, sind mehrere Voraussetzungennotwendig. Zuerst muß die Parametrisierung der Geometrie, im Rahmen ihrergesetzten Grenzen, es ermöglichen, die optimale Konfiguration überhaupt einzunehmen. Diesersehr einfach klingenden Voraussetzung versucht das oben beschriebene Verfahren durchseine sehr weitgehenden Freiheitsgrade, auch was mögliche Amplituden anbelangt, gerecht zuwerden. Die hier auf Erfahrungswerten und Literatur basierten geometrischen Freiheitsgrademüssen mittelfristig aber sicherlich noch einiges an Weiterentwicklung erfahren. EineMethode der diesbezüglichen Ideenwerkstatt könnte z. B. eine Weiterentwicklung der physikalischenice-formation Optimierungsmethode 57 von La Fleur darstellen. Eine mit Eisschichtenversehene Startgeometrie wird dabei im Wasserkanal bei gleichzeitiger Unterkühlung eingesetzt.Der sich einstellende Gleichgewichtszustand stellt ein gleichzeitiges aerodynamischesund thermisches Optimum dar. Bekannte Beispiele zeigen z. B. den gewellten Diffusor, oderzwei hintereinander liegende hufeisenförmige Vertiefungen vor einem umströmten Zylinder.Um eine solche Parametrisierung im Rahmen ingenieurmäßiger Verfahren einzusetzen, ist esaußerdem nötig, daß zur Darstellung dieser topologischen Möglichkeiten nur eine minimaleAnzahl voneinander unabhängiger problemspezifischer Variablen notwendig sind, die sichauch einfach mittels Grenzen im dreidimensionalen Raum limitieren lassen.Der zweite Punkt betrifft die physikalische Modellbildung der durch die gebildete Geometriefließenden Strömung. Die Optimierung der dreidimensionalen Strömungsphänomene erfordertes, daß die physikalische Modellbildung im eingesetzen numerischen Verfahren die dreidimensionalenturbulenten und reibungsbehafteten Strömungseffekte auch physikalisch korrekt auflösenkann. Eine besondere Schwierigkeit besteht darin, auch die vomOptimierungsalgorithmus vorgeschlagenen, an den Grenzen der Freiheitsgrade liegendenextremen Geometrien noch korrekt zu berechnen und eventuell auch als nicht geeignet zubewerten. Hierbei stoßen in einigen Fällen die unter anderen Vorraussetzungen erarbeitetennumerischen Modelle der Turbulenzmodellierung und der Transitionsbestimmung an ihreGrenzen. Es besteht dann die Gefahr, den Optimierungserfolg durch die Optimierung numerischerArtefakte zu unterbinden. Besonders in den letzten Jahren wurden allerdings erhebliche57. Bei der ice-formation Methode wird die Startgeometrie aus einer gekühlten Aluminiumbasis mit einer konstantenEisschicht gebildet. Wird diese Geometrie in einem Wasserkanal betrieben wird sich mit der Zeit eine,sowohl strömungsmechanisch, als auch thermisch, verlustoptimale Eistopologie als Gleichgewicht einstellen.Aus einer Reihe von Versuchen mit unterschiedlichen Unterkühlungsraten der Seitenwände, als auch bei verschiedenenReynolds-Zahlen, kann versucht werden die topologischen Gemeinsamkeiten für die höheren realenReynolds-Zahlen zu extrapolieren. Eine flexible Parametrisierung dieser Geometrie kann dann eine Basiszukünftiger Optimierungen darstellen.103
- Seite 63 und 64: erforderliche Detektierung der rich
- Seite 65 und 66: dimensionalen Auslegung aber nur in
- Seite 67 und 68: Die Optimierung stellt mathematisch
- Seite 69 und 70: den mehrdimensionalen Variablenraum
- Seite 71 und 72: In der Praxis finden sich zahlreich
- Seite 73 und 74: festzustellen, welche Variablen bzw
- Seite 75 und 76: sung nicht exakt gewährleistet. In
- Seite 77 und 78: Das systematische Ablaufschema des
- Seite 79 und 80: Zu- und Abströmgrößen des ebenen
- Seite 81 und 82: Zur Gewährleistung einer symmetris
- Seite 83 und 84: Die Parametrisierung der Schaufelsc
- Seite 85 und 86: dingungen wurde die maximale Machza
- Seite 87 und 88: Abb. 6.5:Gefertigtes Windkanalgitte
- Seite 89 und 90: .Abb. 7.1:Hochgeschwindigkeits-Gitt
- Seite 91 und 92: Schaufeln entfernt sind. Im Rahmen
- Seite 93 und 94: p x - p k, = ------------------c px
- Seite 95 und 96: Die genaue Lage der Profildruckmeß
- Seite 97 und 98: Optimierung durchgeführt wurde. Al
- Seite 99 und 100: Der Mittelschnitt des Profils T106D
- Seite 101 und 102: lichkeiten durch die Freigabe von V
- Seite 103 und 104: In der Abb. 8.5 und Abb. 8.6 sind d
- Seite 105 und 106: Eine Saugspitze an der Vorderkante
- Seite 107 und 108: Abb. 8.10: Ölanstrichbild T106Dopt
- Seite 109 und 110: 8.4 Totaldruckverluste in der Abstr
- Seite 111 und 112: Experimentell ist der s-förmige Ve
- Seite 113: Abb. 8.14: Umfangsgemittelter Total
- Seite 117 und 118: Basiskriterium für die Anwendung e
- Seite 119 und 120: Herausforderungen an eine gute Modu
- Seite 121 und 122: nach Kriterien abgeleitete umfangsg
- Seite 123 und 124: [12] Dennis, B. H.; Dulikravich, G.
- Seite 125 und 126: [34] Harvey, N. W.; Rose, M. G.; Ta
- Seite 127 und 128: [55] Pierret, S.: Designing Turboma
- Seite 129 und 130: [77] Wallis, A. M.; Denton, J. D.:
Impuls aus der Seitenwandsenke füllt die Ablöseblase auf dem Profil und wird ohne Rotationaufzubauen re-energetisiert.Durch die <strong>dreidimensional</strong>e Gestaltung von Profil und Seitenwand konnten die ursächlichenanfachenden Mechanismen der Sekundärströmungen wirkungsvoll reduziert werden. Des weiterenwird ein Aufrollen der Grenzschichten zu den bekannten Wirbelphänomenen vermieden.Die Sekundärströmungsphänomene bilden sich erst weit hinten aus. Es zeigt sich nicht dasübliche Verlustbild in der Profilabströmung. Die Sekundärphänomene bleiben trotz der divergentenSeitenwände in Seitenwandnähe.Die Güte der <strong>Optimierung</strong> läßt sich anhand der integralen Beurteilungsparameter überprüfen.Die wichtigsten integralen Größen sind in der folgenden Übersicht gegenübergestellt:T106D T106Dopt T106Doptnumerisch numerisch experimentellZuströmwinkel 127.7° 127.7° 127.7°Bezugs-Mach-Zahl 0.59 0.59 0.59Bezugs-Reynolds-Zahl 500000 500000 498735Abströmwinkel β 2 grobes Netz 25.15° 24.90° 24.76°<strong>Optimierung</strong> flußgemittelt alsRandbedingungAbströmwinkel β 2 feines NetzNachrechnung / massengemitteltExperiment24.84° 24.57° 24.76°Totaldruckverlustbeiwert feineNachrechnung / ExperimentmassengemitteltTotaldruckverlustbeiwert flußgemitteltfeine Nachrechnung0.0658 0.0591 0.05730.0881 0.0719Die integralen flußgemittelten Verluste konnten durch die <strong>Optimierung</strong> um 22 % relativ zurStartgeometrie reduziert werden. Das Auslegungsziel einer Reduzierung der Sekundärströmungen,bei gleichzeitiger integraler Verlustminimierung, konnte erreicht werden.9. Diskussion und AusblickDurch die Automatisierung der Auslegungskette wird die aerodynamische Berechnung mit derkonstruktiven Entwurfsebene von Profilen gekoppelt. Die Anwendung einer solchen Auslegungsmethodikauf eine vollständig parametrisierte Beschaufelungsgeometrie ist zum heutigen102