Tutaj prezentujemy eksperymentalną metodę rozdzielania współzależnej siły Coriolisa i wyporu obrotowego na rozkłady wymiany ciepła w pełnym polu obracającego się kanału.
Artykuł metodologiczny
Tutaj prezentujemy eksperymentalną metodę rozdzielania współzależnej siły Coriolisa i wyporu obrotowego na rozkłady wymiany ciepła w pełnym polu obracającego się kanału.
Zaproponowano eksperymentalną metodę badania charakterystyki wymiany ciepła w osiowo obracającym się kanale. Decydujące parametry przepływu, które charakteryzują zjawiska transportu w obracającym się kanale, są identyfikowane za pomocą analizy parametrycznej równań pędu i energii odnoszących się do obracającego się układu odniesienia. W oparciu o te bezwymiarowe równania przepływu, sformułowano strategię eksperymentalną, która łączy projekt modułu testowego, program eksperymentalny i analizę danych, z próbą ujawnienia izolowanego wpływu siły Coriolis i wyporu na wydajność wymiany ciepła. Wpływ siły Coriolisa i wyporu wirującego zilustrowano za pomocą selektywnych wyników mierzonych z obracających się kanałów o różnych geometriach. Podczas gdy siła Coriolis i obrotowa pływalność mają kilka wspólnych cech dla różnych obracających się kanałów, unikalne sygnatury wymiany ciepła występują w związku z kierunkiem przepływu, kształtem kanału i rozmieszczeniem urządzeń poprawiających wymianę ciepła. Niezależnie od konfiguracji przepływowych wirujących kanałów, przedstawiona metoda eksperymentalna umożliwia opracowanie fizycznie spójnych korelacji wymiany ciepła, które pozwalają na ocenę izolowanych i współzależnych efektów siły Coriolis i wyporu obrotowego na właściwości wymiany ciepła obracających się kanałów.
Podczas gdy prawa termodynamiczne dyktują poprawę mocy właściwej i sprawności cieplnej silnika turbiny gazowej poprzez podniesienie temperatury wejścia do turbiny, kilka gorących elementów silnika, takich jak łopatki turbiny, jest podatnych na uszkodzenia termiczne. Wewnętrzne chłodzenie łopaty wirnika turbiny gazowej pozwala na temperaturę wejścia do turbiny przekraczającą granice temperaturowe oporu pełzania materiału łopatki. Jednak konfiguracje wewnętrznych kanałów chłodzących muszą być zgodne z profilem łopatek. W szczególności płyn chłodzący obraca się w łopatce wirnika. Przy tak trudnych warunkach termicznych dla pracującej łopaty wirnika turbiny gazowej, skuteczny schemat chłodzenia łopat ma kluczowe znaczenie dla zapewnienia integralności konstrukcji. W związku z tym lokalne właściwości wymiany ciepła dla obracającego się kanału są ważne dla efektywnego wykorzystania ograniczonego dostępnego przepływu chłodziwa. Pozyskiwanie użytecznych danych dotyczących wymiany ciepła, które mają zastosowanie do projektowania wewnętrznych kanałów chłodziwa w realistycznych warunkach silnika, ma zasadnicze znaczenie przy opracowywaniu eksperymentalnej metody pomiaru właściwości wymiany ciepła symulowanego kanału chłodzenia wewnątrz łopatki wirnika turbiny gazowej.
Obrót z prędkością powyżej 10 000 obr./min znacząco zmienia wydajność chłodzenia obracającego się kanału wewnątrz łopaty wirnika turbiny gazowej. Określenie warunków pracy silnika dla takiego obracającego się kanału jest dopuszczalne przy użyciu prawa podobieństwa. Dzięki rotacji bezwymiarowe grupy, które kontrolują zjawiska transportowe wewnątrz promieniowo obracającego się kanału, można ujawnić poprzez wyprowadzenie równań przepływu względem obracającego się układu odniesienia. Morris1 wyprowadził równanie zachowania pędu przepływu względem obracającego się układu odniesienia jako:
(1)
W równaniu (1) lokalna prędkość płynu, v̄, z wektorem położenia, r̄, względem układu odniesienia obracającego się z prędkością kątową, ω, jest pod wpływem przyspieszenia Coriolisa pod względem 2(ω×v̄), odsprzęgniętej dośrodkowej siły wyporu, β(T-Tref)(ω×ω×r̄), napędzany gradient ciśnienia piezometrycznego,
, oraz lepkość dynamiczna płynu, ν. Odnosząca się gęstość płynu, ρref, odnosi się do wstępnie zdefiniowanej temperatury odniesienia płynu Tref, która jest typowa dla lokalnej temperatury objętości płynu do eksperymentów. Jeżeli nieodwracalna konwersja energii mechanicznej na energię cieplną jest pomijalna, równanie zachowania energii sprowadza się do:
(2)
Pierwszy składnik równania (2) otrzymuje się przez traktowanie entalpii właściwej jako bezpośrednio związanej z lokalną temperaturą płynu, T, poprzez stałe ciepło właściwe, Cp. Ponieważ zaburzenie gęstości płynu spowodowane zmianą temperatury płynu w ogrzewanym kanale wirującym ma znaczny wpływ na ruch płynów, gdy łączy się z przyspieszeniem dośrodkowym w równaniu (1), pola prędkości i temperatury płynu w kanale obracającym się osiowo są sprzężone. Ponadto zarówno przyspieszenia Coriolisa, jak i dośrodkowe zmieniają się jednocześnie, gdy prędkość obrotowa jest regulowana. Tak więc wpływ siły Coriolisa i wyporu wirującego na pola prędkości i temperatury płynu jest naturalnie sprzężony.
Równania (1) i (2) w formach bezwymiarowych ujawniają parametry przepływu, które rządzą konwekcją ciepła w obracającym się kanale. Przy zasadniczo równomiernym strumieniu ciepła nałożonym na obracający się kanał, lokalna temperatura objętościowa płynu, Tb, wzrasta liniowo w kierunku strumienia, s, od referencyjnego poziomu wlotowego, Tref. Lokalną temperaturę objętości płynu określa się jako Tref + τs, gdzie τ jest gradientem temperatury objętości płynu w kierunku przepływu. Podstawienia następujących parametrów bezwymiarowych:
(3)
(4)
(5)
(6)
(7)
do równań (1) i (2), gdzie średnia V, N i d oznaczają odpowiednio średni przepływ przez prędkość, prędkość obrotową i średnicę hydrauliczną kanału, bezwymiarowe równania pędu przepływu i energii są wyprowadzane odpowiednio jako równania (8) i (9).
(8)
(9)
Najwyraźniej η w równaniu (9) jest funkcją Re, Ro i Bu = Ro2βτdR, które są odpowiednio określane jako liczby Reynoldsa, rotacji i wyporu. Liczba Rossby'ego, która określa ilościowo stosunek sił bezwładności i sił Coriolisa, jest równoważna odwrotnej liczbie rotacji w równaniu (8).
Gdy Tb jest obliczane jako Tref + τs w obracającym się kanale podlegającym równomiernemu strumieniowi ciepła, wartość τ może być alternatywnie obliczona jako Qf/(mCpL), w której Qf, m i L to odpowiednio konwekcyjna moc grzewcza, masowe natężenie przepływu chłodziwa i długość kanału. Zatem bezwymiarowa lokalna temperatura objętości płynu, ηb, jest równa s / d, a bezwymiarowa temperatura na ścianie kanału, ηw, daje [(Tw-T b)/Qf][mCp][L/d]+s / d. Przy konwekcyjnym współczynniku wymiany ciepła zdefiniowanym jako Qf/(Tw-T b), bezwymiarowa różnica temperatur między ścianą a cieczą, ηw-η b, jest przekształcana na lokalną liczbę Nusselta za pomocą równania (10), w którym ζ jest bezwymiarową funkcją kształtu powierzchni grzewczej i pola przekroju kanału.
(10)
Z zestawem predefiniowanych geometrii oraz hydrodynamicznymi i termicznymi warunkami brzegowymi, bezwymiarowe grupy kontrolujące lokalną liczbę Nusselta obracającego się kanału są identyfikowane jako:
(11)
(12)
(13)
Z eksperymentalnymi testami, dostosowanie prędkości obrotowej, N, dla zmiennego Ro w celu wygenerowania danych o wymianie ciepła przy różnych natężeniach sił Coriolisa, nieuchronnie zmienia przyspieszenie dośrodkowe, a tym samym względną siłę wirującej wyporności. Co więcej, zestaw danych dotyczących wymiany ciepła zebranych z obracającego się kanału zawsze podlega skończonemu stopniowi efektu wyporu obrotowego. Ujawnienie indywidualnego wpływu siły Coriolis i wyporu na wydajność wymiany ciepła przez obracający się kanał wymaga rozłączenia wpływu Ro i Bu na właściwości Nu poprzez procedurę przetwarzania danych końcowych, która jest zawarta w niniejszej metodzie eksperymentalnej.
Warunki przepływu silnika i laboratorium dla obracającego się kanału wewnątrz łopaty wirnika turbiny gazowej mogą być określone przez zakresy Re, Ro i Bu. Typowe warunki silnika dla przepływu płynu chłodzącego przez łopatę wirnika turbiny gazowej, a także budowa i uruchomienie obrotowego obiektu testowego, który pozwolił na przeprowadzanie eksperymentów w pobliżu rzeczywistych warunków silnika, zostały opisane przezMorrisa 2. Opierając się na realistycznych warunkach pracy silnika podsumowanych przez Morrisa2, Rysunek 1 przedstawia realistyczne warunki pracy w zakresie zakresów Re, Ro i Bu dla obracającego się kanału chłodzącego w łopatce wirnika turbiny gazowej. W Rysunek 1, wskazanie najgorszego stanu silnika jest określane jako stan pracy silnika przy najwyższej prędkości obrotowej wirnika i najwyższym współczynniku gęstości. W Rysunek 1, dolna granica i najgorsze warunki pracy silnika pojawiają się odpowiednio przy najniższych i najwyższych prędkościach obrotowych silnika. Niezwykle trudno jest zmierzyć pełnopolowy rozkład Nu w kanale wirującym pracującym z rzeczywistą prędkością obrotową silnika w zakresie od 5000 do 20 000 obr./min. Jednak w oparciu o prawo podobieństwa, testy na skalę laboratoryjną przeprowadzono przy zmniejszonych prędkościach obrotowych, ale z kilkoma próbami zapewnienia pełnego pokrycia zakresów Re, Ro i Bu dla rzeczywistych silników. Jako innowacyjną metodę eksperymentalną, program NASA HOST3,4,5,6 przyjął testy wysokociśnieniowe w celu zwiększenia gęstości płynu w predefiniowanym Re w celu rozszerzenia zakresu Ro poprzez zmniejszenie średniej prędkości płynu. W związku z tym specyficzne zależności między Re, Ro i Bu dla gazu doskonałego o stałej gazowej Rc i lepkości μ są powiązane jako:
(14)
(15)
Aby dostosować warunki laboratoryjne do nominalnej zgodności z warunkami silnika widzianymi w Rysunek 1, prędkość obrotowa, N, ciśnienie płynu chłodzącego, P, średnica hydrauliczna kanału, d, promień obrotu, R, oraz różnica temperatur między ścianą a cieczą, Tw-T b, muszą być kontrolowane w celu dopasowania realistycznych Re, Ro i Bu Zakresy. Oczywiste jest, że jednym z najskuteczniejszych podejść do rozszerzenia zakresu Ro jest zwiększenie średnicy hydraulicznej kanału, ponieważ Ro jest proporcjonalne do d2. Ponieważ laboratoryjny test wymiany ciepła przy realistycznym N jest niezwykle trudny, ciśnienie chłodziwa, P, jest technicznie łatwiejsze do podniesienia w celu rozszerzenia zakresu Ro; nawet jeśli Ro jest tylko proporcjonalne do P. Opierając się na tym zapleczu teoretycznym, filozofia projektowania obecnej metody eksperymentalnej polega na zwiększeniu Ro poprzez zwiększenie ciśnienia w obracającym się kanale testowym przy użyciu maksymalnej średnicy hydraulicznej kanału, jaka może zmieścić się w obracającej się platformie. Po zwiększeniu zakresu Ro, zakres Bu jest odpowiednio rozszerzany, ponieważ Bu jest proporcjonalne do Ro2. W Rysunek 1, uwzględniono również warunki badań laboratoryjnych przyjęte do generowania danych o wymianie ciepła w obracających się kanałach3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22, 23,24,25,26,27,28,29. Jak wskazano w Rysunek 1, pokrycie realistycznych warunków pracy silnika przez dostępne dane dotyczące wymiany ciepła jest nadal ograniczone, szczególnie dla wymaganego zakresu Bu. Symbole otwartych i kolorowych brył przedstawione w Rysunek 1 to odpowiednio eksperymenty wymiany ciepła w spiczastym i pełnym polu. Jak zebrano w Rysunek 1, większość danych dotyczących wymiany ciepła z aplikacjami chłodzącymi do łopat wirnika turbiny gazowej1,2,3,4,5,6,7,8,9, 10,11,12,13,14,15,16,17,18,20,21,22,23,24,25,26 są pomiarami punktowymi wykonanymi metodą termopary. Wpływ przewodzenia ściennego na pomiar strumienia ciepła przewodzącego ścianę oraz temperatury na granicy faz ciecz-ściana pogarszają jakość danych dotyczących wymiany ciepła przekształconych z pomiarów termopary. Ponadto pomiary przenikania ciepła1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,20,21,22,23,24, 25,26 przy użyciu metody termopary nie może wykryć dwuwymiarowych zmian wymiany ciepła na obracającej się powierzchni. Przy obecnej metodzie eksperymentalnej29,30,31,32, wykrywanie pełnopolowych rozkładów liczby Nusselta nad obracającą się ścianą kanału jest dopuszczalne. Minimalizacja efektu przewodzenia ścian przy użyciu folii ze stali nierdzewnej o grubości 0,1 mm o liczbie Biot >>1 do generowania mocy grzewczej niniejszą metodą eksperymentalną pozwala na jednowymiarowe przewodzenie ciepła od folii grzewczej do przepływu chłodziwa. W szczególności niedopuszczalne jest pozyskiwanie danych dotyczących wymiany ciepła w pełnym polu, obejmujących zarówno efekty Ro, jak i Bu, przy użyciu techniki ciekłokrystalicznej w warunkach nieustalonych i metody termopary. Przy obecnej metodzie termografii ciekłokrystalicznej w stanie ustalonym19, wykrywalny zakres temperatur 35-55 °C uniemożliwia generowanie danych o wymianie ciepła z realistycznymi współczynnikami gęstości.
Używanie parametrów przepływu rządzących konwekcją ciepła w obracającym się kanale w celu wykazania, że pełne pokrycie realistycznych warunków silnika widoczne w Rysunek 1 nie zostało jeszcze osiągnięte, więc potrzeba uzyskania danych o wymianie ciepła w pełnym polu w realistycznych warunkach silnika była stale naglona. Niniejsza metoda eksperymentalna umożliwia generowanie wymiany ciepła w pełnym polu z wykrytymi efektami zarówno siły Coriolisa, jak i wyporu obrotowego. Protokoły mają na celu pomoc badaczom w opracowaniu strategii eksperymentalnej istotnej dla realistycznego pomiaru wymiany ciepła w całym polu obracającego się kanału. Wraz z metodą analizy parametrycznej, która jest unikalna dla niniejszej metody eksperymentalnej, dozwolone jest generowanie korelacji wymiany ciepła w celu oceny izolowanych i współzależnych efektów Ro i Bu na Nu.
Artykuł ilustruje eksperymentalną metodę mającą na celu generowanie dwuwymiarowych danych o wymianie ciepła w wirującym kanale o warunkach przepływu podobnych do realistycznych warunków silnika turbiny gazowej, ale działającego przy znacznie niższych prędkościach obrotowych w laboratoriach. We wstępie przedstawiono metodę opracowaną w celu wyboru prędkości obrotowej, średnicy hydraulicznej kanału testowego oraz zakresu różnic temperatur między ścianką a cieczą w celu uzyskania danych dotyczących wymiany ciepła w realistycznych warunkach silnika. Przedstawiono testy kalibracyjne systemu termografii w podczerwieni, testy kalibracji strat ciepła oraz działanie obrotowego stanowiska do badania wymiany ciepła. W artykule opisano czynniki powodujące znaczne niepewności dla pomiarów wymiany ciepła oraz procedury rozdzielania wpływu siły Coriolis i wyporu na właściwości wymiany ciepła kanału wirującego wraz z wynikami selektywnymi w celu zademonstrowania niniejszej metody doświadczalnej.
Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.
UWAGA: Szczegóły dotyczące rotacji obiektów testowych, pozyskiwania danych, przetwarzania danych i modułu testowego wymiany ciepła emulującego wewnętrzny kanał chłodzenia łopaty wirnika turbiny gazowej znajdują się w naszych poprzednich pracach29,30,31,32.
1. Przygotowanie testów wymiany ciepła
2. Wyznaczanie współczynnika emisyjności cieplnej dla systemu termografii w podczerwieni
3. Dynamiczna równowaga obracającego się zestawu
4. Ocena współczynników strat ciepła
5. Podstawowe testy wymiany ciepła
6. Testy wymiany ciepła przy rotacji
7. Analiza parametryczna
Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.
Realistyczne warunki pracy wewnętrznego przepływu chłodziwa wewnątrz obracającej się łopatki turbiny gazowej pod względem Re, Ro i Bu są porównywane z emulowanymi warunkami laboratoryjnymi w Rysunek 1. Punkty danych mieszczą się w realistycznych warunkach silnika przy użyciu obecnej metody eksperymentalnej podsumowanej w protokołach11,14,17,20...
Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.
Podczas gdy temperatury ścianek czołowych obracającego się kanału są wykrywane przez system termografii w podczerwieni, temperatury płynu są mierzone przez termopary. Ponieważ alternatywne pole magnetyczne silnika prądu przemiennego, który napędza obracającą się platformę, indukuje potencjał elektryczny, który zakłóca pomiary termopary, silnik prądu stałego musi być przystosowany do napędzania obracającego się stanowiska testowego.
Rozkład temperatury płynu w płaszczyźnie wyjściowej ogrzewaneg...
Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.
Autorzy nie mają nic do ujawnienia.
Obecna praca badawcza była finansowana przez Ministerstwo Nauki i Technologii Tajwanu w ramach grantu NSC 94-2611-E-022-001, NSC 95-2221-E-022-018, NSC 96-2221-E-022-015MY3 i NSC 97-2221-E-022-013-MY3.
Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.
| Nazwa | Firma | Numer katalogowy | Komentarze |
|---|---|---|---|
| Obrotowe stanowisko testowe | Wykonane | zakresie Projekt przez tę grupę badawczą | |
| Moduł testowy wymiany ciepła | Wykonany | we własnym zakresie Projekt przez tę grupę badawczą | |
| Przepływomierz masowy | Eldride Product, Inc. | 3100301-01-01 359-1007 | |
| System termografii w podczerwieni | NEC P384A-8 | 3100401-04 3127A-4 | |
| Pierścień ślizgowy oprzyrządowania | Michigan Scientific SR36M | 3100506-62 3553-372 |
Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.
Poproś o pozwolenie na ponowne wykorzystanie tekstu lub ilustracji tego artykułu JoVE
Poproś o pozwolenie