Miniaturowy silnik turboodrzutowy nie jest pomniejszonym egzemplarzem większej maszyny w tym samym sensie, w jakim pomniejsza się rysunek. Jeżeli wszystkie długości zmniejszyć w jednej proporcji, podobieństwo geometryczne zostanie zachowane, ale podobieństwo przepływu, wymiany ciepła, spalania, dynamiki wirnika i pomiaru — nie. Utrzymanie zbliżonej liczby Macha wymaga innej drogi skalowania niż utrzymanie liczby Reynoldsa; tolerancja, sonda i warstwa przyścienna nie muszą maleć wraz z kanałem; czas gazodynamiczny skraca się, podczas gdy czas chemiczny i czas odpowiedzi czujnika nie podlegają temu samemu prawu.

To jest zasadniczy problem tej monografii: miniaturyzacja silnika turboodrzutowego jest problemem niezgodnych podobieństw. Wiarygodny wniosek nie brzmi więc „małe silniki mają większe straty”, lecz: która niezgodność podobieństwa dominuje w oznaczonej konfiguracji, jak zmierzono jej skutek, jaka jest niepewność wyniku i dokąd wolno ten wynik przenieść? Odpowiedź jest rozwijana na trzech łańcuchach dowodowych: bezwymiarowym studium syntetycznym, historycznym cywilnym GTM-120 oraz przypadku, w którym wynik jest hybrydą pomiaru i modelu.

Tekst dotyczy fizyki, metrologii i oceny dowodu. Nie jest instrukcją projektowania, budowy, uruchamiania, strojenia ani prób napędu współczesnego uzbrojenia. Nie podaje geometrii wykonawczej komory, tolerancji produkcyjnych, materiałów konkretnego wirnika, map paliwowych, nastaw rozruchu, progów ochronnych ani kompletnej procedury stanowiskowej.

Pytanie badawcze, definicja miniaturyzacji i stała granica systemu

Roboczo przez miniaturowy silnik turboodrzutowy rozumiany jest tu jednoprzepływowy silnik turbinowy o tak małej skali liniowej i strumieniu masy, że co najmniej dwa efekty nieskalujące się geometrycznie — na przykład liczba Reynoldsa, luz względny, bezwładność cieplna, czas procesu albo blokada sondą — wpływają na wynik w stopniu porównywalnym ze zmianą badanego parametru. Nie jest to definicja normatywna ani próg masy lub ciągu. Celowo nie przypisuje etykiety „miniaturowy” na podstawie jednej liczby katalogowej, ponieważ ta sama maszyna może być mała względem lotniczego turbojetu, a duża względem układu mikroturbinowego.

Miniaturyzacja oznacza w tej analizie przejście między konfiguracjami różniącymi się charakterystycznym wymiarem $L$, przy jawnym wskazaniu tego, co próbuje się zachować: geometrię bezwymiarową, stan całkowity na wlocie, prędkość końcówki, skorygowany przepływ, obciążenie stopnia, liczbę Reynoldsa lub jeszcze inną wielkość. Sama relacja $L_B/L_A$ nie definiuje eksperymentu. Dwa obiekty o tej samej skali mogą różnić się architekturą, technologią luzów, chropowatością, sposobem chłodzenia i instalacją tak głęboko, że porównanie „duży–mały” nie izoluje żadnego mechanizmu.

W całym tekście obowiązuje następująca, niezmienna księga granic:

Obiekt Granica przyjęta w tekście Co należy do obiektu Czego nie wolno dopisywać przez analogię
rdzeń silnika przekroje od wejścia do sprężarki do wylotu z dyszy oraz wspólny zespół wirujący sprężarka, dyfuzor/kanał nawrotny, komora, turbina, wał i dysza objęte bilansem instalacyjnego spadku ciśnienia wlotu, masy zbiornika, sterownika i płatowca
kompletny zespół napędowy z instalacją osłona obejmująca rdzeń, rozruch, zasilanie paliwem, sterowanie, przewody, czujniki i zabudowany wlot/wylot masy, pobory i straty wszystkich wymienionych podzespołów oporu i zniekształcenia powodowanego przez cały płatowiec, jeśli nie zostały zmierzone
integracja z płatowcem objętość kontrolna wokół zespołu napędowego w rzeczywistym kanale i polu opływu odzysk ciśnienia wlotu, opór pędu strumienia wlotowego (ang. ram drag), upusty, oddziaływanie dyszy i opór instalacyjny osiągów „gołego” rdzenia jako osiągów lotnych
GTM-120 historyczny cywilny egzemplarz/rodzina opisana przez Gierasa oraz dwie wskazane publikacje tylko wersja, stanowisko, stacje i warunki jawnie udokumentowane cech SO-3, silnika NASA, innego GTM ani współczesnego napędu wojskowego
SO-3 historyczny silnik lotniczy omawiany przez Gierasa w odrębnym przykładzie diagnostycznym tylko dane przypisane w źródle do SO-3 charakterystyki tarcia, sterowania lub dynamiki GTM-120
NASA TM X-3463 historyczny eksperymentalny turbojet o sprężarce osiowej opisany w raporcie z 1976 r. jego własna architektura, stanowisko i zakres wskaźnika Reynoldsa cech GTM-120 albo współczesnego napędu wojskowego
konfiguracje dydaktyczne A–D jawnie syntetyczne, bezwymiarowe obiekty zdefiniowane w Studium I wyłącznie relacje wynikające z podanych równań i współczynników autora wymiarów, osiągów lub tolerancji realnego silnika
współczesny napęd wojskowy dokładnie oznaczona wersja sprzętu wraz z zatwierdzoną dokumentacją tej wersji tylko twierdzenia udokumentowane dla tej wersji i konfiguracji jakiegokolwiek ilościowego transferu z obiektu cywilnego, historycznego lub syntetycznego

Słowo „silnik” bez dopowiedzenia oznacza rdzeń tylko tam, gdzie równanie ma granicę od stacji 2 do 9. „Napęd” nie jest automatycznie synonimem rdzenia. Ta dyscyplina języka zapobiega częstemu błędowi: dodaniu katalogowego ciągu rdzenia do równania lotu i pominięciu wlotu, oporu pędu strumienia wlotowego, poborów oraz oporu zabudowy.

Silnik turboodrzutowy pobiera utleniacz z atmosfery, natomiast silnik rakietowy niesie go ze sobą. Zmienia to granicę bilansu masy i zależność napędu od własności oraz działania wlotu: w turbojecie strumień powietrza przecina granicę układu, a w silniku rakietowym zapas utleniacza należy do pojazdu. To rozróżnienie jest potrzebne do bilansów, lecz nie wymaga odtwarzania katalogu odmiennych architektur napędów.

Każdej liczbie przypisana jest jedna z pięciu etykiet pochodzenia:

  • pomiar — odczyt wielkości przez opisany łańcuch pomiarowy;
  • dane publikowane — wartość zadeklarowana w katalogu, tabeli producenta albo publikacji bez twierdzenia, że autor artykułu niezależnie ją zmierzył;
  • wynik modelu — wyjście obliczenia zależne od równań, siatki, warunków brzegowych i kalibracji;
  • wartość wyznaczona pośrednio — rezultat redukcji innych pomiarów albo połączenia pomiaru z modelem;
  • dane syntetyczne autora — liczby dobrane wyłącznie do demonstracji metody, nieopisujące żadnego istniejącego silnika.

„Zgodność modelu z pomiarem” nie zmienia automatycznie wyniku modelu w pomiar. „Dane publikowane” również nie znaczą „niezależnie potwierdzone”. Te dwa rozróżnienia są konieczne do interpretacji GTM-120.

Stacje termodynamiczne i dwie objętości kontrolne

Do dalszych rachunków wystarcza siedem umownych stacji. Indeks $t$ oznacza wielkość całkowitą, a brak indeksu — statyczną. Stacja 0 leży w niezaburzonym otoczeniu przed zespołem; 2 na wejściu do sprężarki; 3 za sprężarką; 4 za komorą, przed turbiną; 5 za turbiną; 8 w charakterystycznym przekroju dyszy; 9 w płaszczyźnie wylotowej. Numeracja jest konwencją tej monografii. Przy porównaniu z publikacją najpierw trzeba odwzorować jej fizyczne przekroje, a dopiero potem przenieść symbole. Jednakowe cyfry w dwóch źródłach nie gwarantują jednakowego miejsca.

Pierwsza objętość kontrolna obejmuje rdzeń 2–9. Jej księga masy w stanie ustalonym ma postać

$$\dot m_9=\dot m_2+\dot m_f-\dot m_{b}-\dot m_{l},$$

gdzie $\dot m_2$ jest strumieniem powietrza na wejściu sprężarki [kg/s], $\dot m_f$ strumieniem paliwa [kg/s], $\dot m_b$ strumieniem wyprowadzonym poza granicę przez upust [kg/s], a $\dot m_l$ pozostałym strumieniem opuszczającym granicę inną drogą [kg/s]. Jeżeli źródło nie opisuje upustów i przecieków, nie wolno bez śladu przyjąć ich za zero. Można jedynie zadeklarować uproszczenie i zbadać jego wpływ.

Bilans energii wirującego zespołu, z jednoznaczną granicą mocy gazowej, zapisuje się jako

$$\frac{\mathrm d}{\mathrm dt}\left(\frac{1}{2}J\omega^2\right) =\dot m_g(h_{t4}-h_{t5})-\dot m_a(h_{t3}-h_{t2})-P_{\mathrm{loss}}-P_{\mathrm{aux}}.$$

$J$ jest biegunowym momentem bezwładności zespołu [kg·m²], $\omega$ prędkością kątową [rad/s], $h_t$ entalpią całkowitą [J/kg], $P_{\mathrm{loss}}$ mocą wszystkich strat mechanicznych przekraczających zdefiniowaną granicę [W], a $P_{\mathrm{aux}}$ mocą odbiorników [W]. Strumień $\dot m_g$ przecina fizyczny przekrój turbiny od stacji 4 do 5, a $\dot m_a$ — przekrój sprężarki od stacji 2 do 3; z powodu dopływu paliwa, upustów i przecieków nie muszą być sobie równe. Różnice entalpii są rzeczywistymi zmianami energii gazu; nie mnoży się ich ponownie przez sprawność izentropową. Sprawności wykorzystuje się przy przejściu między stanem rzeczywistym i izentropowym, nie jako drugi potrącony składnik. W stanie ustalonym lewa strona jest równa zeru. W przejściu dynamicznym nie jest, dlatego równość mocy turbiny i sprężarki nie może zastępować równania przyspieszenia.

Druga objętość kontrolna obejmuje kompletny zespół napędowy i przecina niezaburzony strumień w stacji 0 oraz wylot w stacji 9. Dla ustalonego, jednowymiarowo uśrednionego przepływu, z pominięciem pędu paliwa i osobnych wylotów pomocniczych, siła osiowa wywierana na zespół ma wartość

$$F_N=\dot m_9V_9-\dot m_0V_0+(p_9-p_0)A_9-F_{\mathrm{inst}},$$

gdzie $V$ oznacza średnią osiową prędkość [m/s], $p$ ciśnienie statyczne [Pa], $A_9$ pole wylotu [m²], a $F_{\mathrm{inst}}$ jest dodatnią wartością wypadkowej siły instalacyjnej skierowanej przeciwnie do dodatniej osi ciągu i grupuje tylko te oddziaływania, których nie objęto wybranymi powierzchniami kontrolnymi. Jeżeli objętość kontrolna zostanie poprowadzona inaczej, zmienią się jawne człony ciśnieniowe i oporowe, lecz wynik fizyczny nie powinien. To równanie służy później do wykazania, dlaczego spadek ciśnienia całkowitego wewnątrz rdzenia nie przekłada się na ciąg jednym uniwersalnym współczynnikiem.10

W próbie statycznej $V_0\approx0$, ale w locie człon $-\dot m_0V_0$, czyli opór pędu strumienia wlotowego, nie znika. Jednocześnie zabudowany wlot może zmienić $p_{t2}$, profil prędkości i margines stateczności. Dlatego wynik stanowiskowy i wynik lotny należą do dwóch różnych objętości kontrolnych, nawet gdy rdzeń jest ten sam.

Dwa często mylone wskaźniki mają różne wymiary. $F/m$ [N/kg] jest ciągiem przypadającym na jednostkę jawnie wskazanej masy $m$, na przykład suchego rdzenia albo kompletnego zespołu napędowego. $F/W$ jest bezwymiarowym stosunkiem ciągu do ciężaru tej samej masy, ponieważ $W=mg_0$, zatem

$$\frac{F}{W}=\frac{F/m}{g_0}.$$

Bez określenia, czy mianownik obejmuje suchy rdzeń, instalację, czy cały zespół, obu wskaźników nie da się uczciwie porównywać.

Droga skalowania: co zostaje stałe, a co musi się zmienić

Niech konfiguracja B powstaje z A przez zmianę wszystkich idealnie skalowanych długości według $L_B=\lambda L_A$, gdzie $0<\lambda<1$. Dla podobnych brył $A_B/A_A=\lambda^2$, $V_B/V_A=\lambda^3$. Jest to dopiero geometria. Aby zdefiniować drogę aerodynamiczną, przyjmijmy ten sam gaz, tę samą temperaturę całkowitą i ciśnienie całkowite na wlocie oraz zbliżoną prędkość końcówki $U=\omega r$. Wówczas prędkość dźwięku $a=\sqrt{\gamma RT}$ pozostaje zbliżona, więc podobny lokalny udział $U/a$ pomaga utrzymać podobieństwo Macha. Gieras opisuje te jakościowe skutki miniaturyzacji — zmianę znaczenia Reynoldsa, warstwy przyściennej, wymiany ciepła, czasu procesu i prędkości obrotowej — na s. 43–54 druku / 44–55 PDF.11

Cena tego wyboru jest natychmiastowa:

$$\frac{\omega_B}{\omega_A}=\frac{1}{\lambda},\qquad \frac{n_B}{n_A}=\frac{1}{\lambda},$$

gdzie $n$ jest prędkością obrotową [obr/min]. Dwukrotne zmniejszenie promienia przy tym samym $U$ wymaga dwukrotnie większej liczby obrotów na minutę. Sama wartość $n$ nie jest jednak podobieństwem aerodynamicznym. Jest konsekwencją wybranej prędkości końcówki i wymiaru.

Lokalna liczba Reynoldsa

$$Re=\frac{\rho VL}{\mu}$$

przy tym samym stanie wlotowym, gazie i podobnej prędkości maleje jak $\lambda$. Aby jednocześnie utrzymać $Re$, trzeba by zwiększyć $\rho$, zwiększyć $V$, zmienić $\mu$ albo użyć instalacji ciśnieniowej. Zwiększenie $V$ niszczy podobieństwo Macha; zwiększenie gęstości zmienia warunki mechaniczne, wymianę ciepła i często punkt pracy; zmiana gazu narusza podobieństwo właściwości i spalania. Nie istnieje jeden ruch pokrętła przywracający wszystkie warunki.4,11

Strumień masy dla podobnego pola przepływu skaluje się w tej drodze jak

$$\dot m\sim \rho V A\propto \lambda^2.$$

Charakterystyczny czas przepływu $\tau_f=L/V$ maleje jak $\lambda$. Energia kinetyczna wirnika przy idealnym podobieństwie gęstości materiału i stałym $U$ skaluje się w przybliżeniu jak masa, czyli $\lambda^3$, ale częstotliwości własne, tłumienie, sztywność łożysk i niewyważenie nie wynikają z tego jednego wykładnika. W realnej miniaturyzacji wał, łożysko, przewód, powłoka i minimalna grubość nie zachowują idealnego podobieństwa.

Luz nad końcówką należy zapisać jako wielkość względną

$$\varepsilon=\frac{c}{h},$$

gdzie $c$ jest lokalnym prześwitem, a $h$ wysokością kanału lub łopatki. Gdy $h\propto\lambda$, lecz absolutny luz jest ograniczony technologią, rozszerzalnością, biciem i dynamiką, $\varepsilon$ rośnie szybciej niż wynikałoby z podobieństwa geometrycznego. Przeciek nie jest wówczas „tym samym przepływem w mniejszej szczelinie”. Zajmuje większą część użytecznego przekroju i silniej oddziałuje z warstwą przyścienną.5,6

Analogicznie chropowatość względna $k_s/L$, średnica sondy do wysokości kanału $d_p/h$, grubość ściany do promienia oraz bezwładność cieplna ściany do czasu próby mogą rosnąć. Rzeczywista droga skalowania musi więc być tabelą założeń, nie pojedynczą strzałką:

Wielkość Na przyjętej drodze $U,T_{t2},p_{t2}\approx \mathrm{const}$ Status podobieństwa
lokalny Mach oparty na $U/a$ w przybliżeniu stały celowo zachowany, ale nie gwarantuje identycznych trójkątów prędkości
liczba obrotów $n$ $\lambda^{-1}$ musi wzrosnąć
Reynolds $Re$ $\lambda$ naruszony
strumień masy $\dot m$ $\lambda^2$ zmniejszony
czas przepływu $\tau_f$ $\lambda$ skrócony
luz względny $c/h$ około $\lambda^{-1}$, jeśli $c\approx \mathrm{const}$ naruszony
blokada powierzchniowa sondą około $\lambda^{-2}$, jeśli pole sondy nie maleje gwałtownie naruszona
czas odpowiedzi czujnika / czas procesu około $\lambda^{-1}$, jeśli czujnik się nie zmienia naruszony

Ostatnia kolumna jest właściwym początkiem analizy. Mówi, jakie niezgodności trzeba osobno zmierzyć lub ograniczyć, zanim wynik z A zostanie użyty dla B.

Konflikt Mach–Reynolds–końcówka–luz–ciepło–czas

Liczba Macha $M=V/a$ kontroluje znaczenie ściśliwości. Liczba Reynoldsa porównuje transport pędu przez bezwładność i lepkość. Prędkość końcówki $U$ wiąże pracę stopnia, względne prędkości i naprężenia z ruchem wirnika. Luz względny $\varepsilon$ otwiera drogę przecieku omijającego użyteczną wymianę pracy. Żadna z tych wielkości nie jest zamiennikiem pozostałych.

Próba odzyskania $Re$ przez wzrost prędkości o $1/\lambda$ podnosi Mach mniej więcej o $1/\lambda$, jeśli temperatura pozostaje stała. Próba odzyskania go przez wzrost gęstości o $1/\lambda$ wymaga innego ciśnienia stanowiska i zmienia obciążenia oraz przepływy przeciekowe. Próba odzyskania luzu względnego przez zmniejszenie absolutnego luzu do $\lambda c$ może zderzyć się z nieskalującą się rozszerzalnością, dokładnością wykonania i przemieszczeniami wirnika. Nie jest zatem poprawne twierdzenie, że jeden eksperyment „zachował skalę”, bez listy zachowanych grup.

Konflikt dotyczy również ciepła. Sama proporcja $A/V\propto1/L$ informuje jedynie, że powierzchnia przypadająca na jednostkę objętości rośnie. Nie wyznacza straty cieplnej. Dla uproszczonego, ustalonego, jednowymiarowego przejścia przez dwie warstwy konwekcyjne i płaską ścianę o stałych własnościach strumień ciepła na jednostkę powierzchni wynosi

$$q''=\frac{T_g-T_\infty}{\displaystyle \frac{1}{h_g}+\frac{t_w}{k_w}+\frac{1}{h_o}}, \qquad \dot Q=q''A,$$

gdzie $h_g$ i $h_o$ są współczynnikami przejmowania po stronie gazu i otoczenia [W/(m²·K)], $t_w$ grubością ściany [m], a $k_w$ jej przewodnością [W/(m·K)]. Współczynnik $h_g$ zależy przez liczbę Nusselta $Nu=h_gL/k_g$ od $Re$, liczby Prandtla $Pr=\mu c_p/k_g$, geometrii, stanu warstwy przyściennej i nierównomierności. Grubość ściany i opór zewnętrzny nie muszą się skalować. Dlatego wniosek z samego $A/V$ jest niezamknięty.12

Wymiarowo poprawną miarą znaczenia utraty ciepła dla gazu jest na przykład

$$\chi_Q=\frac{\dot Q}{\dot m c_p\Delta T_{\mathrm{procesu}}}.$$

Przy podobnej geometrii stosunek pól powierzchni wymiany i przepływu jest stały, a $h_g/(\rho Vc_p)$ można powiązać z $Nu/(RePr)$. Dopiero ta kombinacja pokazuje, czy ściana odbiera istotną część entalpii strumienia. Dla stanu nieustalonego dochodzą pojemność cieplna ściany $C_w=\rho_wc_wV_w$, opór cieplny $R_{th}$, czas $\tau_w=R_{th}C_w$, liczba Biota $Bi=h_gt_w/k_w$ oraz liczba Fouriera $Fo=\alpha_wt/t_w^2$, gdzie $\alpha_w$ jest dyfuzyjnością cieplną ściany [m²/s], a $t$ czasem [s]. Próba trwająca ułamek $\tau_w$ i próba po nasyceniu cieplnym nie mierzą tego samego stanu, nawet przy tej samej prędkości obrotowej.12

Konflikt czasów jest równie fundamentalny. Lokalny czas przebywania najlepiej definiować przez inwentarz masy w oznaczonej objętości $V_c$:

$$\tau_{res}=\frac{\displaystyle\int_{V_c}\rho\,\mathrm dV}{\dot m_c}.$$

Jest to definicja uśredniona; rozkład czasów przebywania może być szeroki i zawierać recyrkulację. Czas chemiczny $\tau_{chem}$ jest czasem charakterystycznym wybranego etapu reakcji dla określonej temperatury, ciśnienia i składu. W rzeczywistej komorze wybór właściwej długości strefy recyrkulacji, prędkości odniesienia i czasu chemicznego sam stanowi część modelu, a nie oczywistą operację. Liczba Damköhlera

$$Da=\frac{\tau_{res}}{\tau_{chem}}$$

porównuje czas dostępny z czasem reakcji. Na przyjętej drodze geometrycznej $\tau_{res}\sim L/V\propto\lambda$, lecz $\tau_{chem}$ nie musi maleć; co więcej, większa względna utrata ciepła może obniżyć temperaturę i przez kinetykę reakcji wydłużyć czas chemiczny nieliniowo. Miniaturyzacja sprzęga więc czas i bilans cieplny, zamiast tworzyć dwie niezależne poprawki.13

Czas obrotu $\tau_{rot}=2\pi/\omega$ maleje jak $\lambda$ przy stałym $U$. Jeżeli czas odpowiedzi sondy $\tau_s$ pozostaje stały, stosunek $\tau_s/\tau_{rot}$ rośnie jak $1/\lambda$. Czujnik może wygładzić szybką fluktuację, a jego trzon może jednocześnie zmienić pole, które miał obserwować. Metrologia jest zatem kolejnym uczestnikiem konfliktu podobieństw, nie neutralnym oknem.

Od lokalnej straty do sprężarki, turbiny, dyszy i ciągu

Lokalną stratę ciśnienia całkowitego między przekrojami $a$ i $b$ można opisać przez współczynnik

$$\zeta_{a-b}=\frac{p_{ta}-p_{tb}}{p_{ta}}, \qquad \sigma_{a-b}=\frac{p_{tb}}{p_{ta}}=1-\zeta_{a-b}.$$

Definicja wymaga jednak takiego samego sposobu uśredniania ciśnienia całkowitego w obu przekrojach. Średnia powierzchniowa, masowa i strumieniowa mogą się różnić w przepływie niejednorodnym. Pomiar jednopunktowy nie staje się średnią masową tylko dlatego, że został wykonany w „reprezentatywnym” miejscu.

Bilans komory, przy jawnie pominiętej energii kinetycznej paliwa, entalpii paliwa na wlocie względem przyjętego stanu odniesienia oraz przy wszystkich upustach uwzględnionych w strumieniach, można zapisać jako

$$\dot m_g h_{t4}=\dot m_a h_{t3}+\eta_b\dot m_fH_u-\dot Q_w,$$

gdzie $\eta_b$ oznacza wyłącznie chemiczny stopień uwolnienia energii paliwa w objętości kontrolnej, $H_u$ wartość opałową [J/kg], a $\dot Q_w\geq0$ ciepło opuszczające objętość kontrolną przez ściany [W]. Strata ścienna nie jest zawarta w $\eta_b$, dlatego odejmuje się ją osobno i tylko raz. Ta równość pokazuje, że taka sama temperatura $T_{t4}$ może wynikać z innej kombinacji strumienia paliwa, niecałkowitego spalania i strat do ściany. Temperatura nie identyfikuje samodzielnie mechanizmu.11,12

Strata $\zeta_{3-4}$ obniża $p_{t4}$. Turbina nadal musi pokryć rzeczywistą pracę sprężarki i odbiorników. W punkcie dopasowania zmieniają się więc jednocześnie stosunek rozprężania turbiny, $T_{t5}$, prędkość skorygowana i warunki dyszy. Przyjmijmy przepływ ustalony, adiabatyczną dyszę bez pracy wału, jednowymiarowe uśrednienie wylotu i rozprężenie idealnego stanu odniesienia od stanu całkowitego w stacji 8 do zadanego ciśnienia statycznego $p_9$. Jeżeli $h_{9s}$ jest entalpią statyczną po takim izentropowym rozprężeniu, a $\eta_n$ oznacza stosunek rzeczywistej energii kinetycznej wylotu do idealnego izentropowego spadku entalpii od $h_{t8}$ do $h_{9s}$, to prędkość wylotowa ma postać

$$V_9=\sqrt{2\eta_n(h_{t8}-h_{9s})}.$$

Energia kinetyczna odpowiadająca prędkości w stacji 8 jest już zawarta w entalpii całkowitej $h_{t8}$; nie dodaje się więc osobnego składnika $V_8^2$. Mniejsze dostępne ciśnienie całkowite może zmniejszyć możliwy spadek entalpii, człon pędu i człon ciśnieniowy w równaniu ciągu. Jeżeli układ sterowania próbuje utrzymać prędkość lub ciąg przez zwiększenie dopływu paliwa, wzrasta $\dot m_f$, temperatura albo oba parametry — aż do ograniczenia właściwego dla danej konfiguracji. Nie ma uniwersalnego przelicznika „jeden punkt procentowy straty ciśnienia = x procent ciągu”. Potrzebne jest ponowne dopasowanie całego układu.14

Zużycie paliwa właściwe dla ciągu

$$TSFC=\frac{\dot m_f}{F_N}$$

rośnie zarówno wtedy, gdy dla podobnego $\dot m_f$ maleje ciąg, jak i wtedy, gdy sterowanie zwiększa $\dot m_f$ szybciej niż odzyskuje $F_N$. Wniosek o TSFC wymaga wspólnego czasu, tej samej granicy paliwa i ciągu oraz ich kowariancji. Oddzielne wykresy z różnych przebiegów nie tworzą automatycznie jednego punktu.

Pełny łańcuch ma zatem postać: lokalny mechanizm → rozkład straty i nierównomierności → uśredniona strata całkowita → nowe dopasowanie sprężarki, turbiny i dyszy → reakcja sterowania → ciąg i zużycie paliwa w określonej objętości kontrolnej. Pominięcie któregokolwiek ogniwa zamienia mechanizm w sugestię, nie w dowód osiągu.

Co rzeczywiście wiadomo o niskim Reynoldsie, luzie i niestateczności

Historyczny eksperyment opisany w NASA TM X-3463 jest wartościowy dlatego, że zmieniano warunki wlotowe tego samego silnika i analizowano skorygowaną charakterystykę względem zdefiniowanego wskaźnika liczby Reynoldsa. Dla wyższej części badanego zakresu przepływ i sprawność sprężarki dawały się uogólniać, natomiast przy niższych wartościach wskaźnika załamanie podobieństwa stawało się wyraźne; stosunek ciśnień zachowywał się inaczej niż sprawność i przepływ.4 To pomiar historycznego silnika osiowego, nie prawo progowe dla GTM-120. Dowodzi natomiast ważnej zasady: dwie wielkości wyjściowe tej samej sprężarki mogą utracić podobieństwo w różnym stopniu.

Xiang, Schluter i Duan rozdzielili numerycznie dwa czynniki w sprężarkach miniaturowych: najpierw zmieniali luz przy utrzymaniu liczby Reynoldsa, a następnie liczbę Reynoldsa przy ustalonym luzie absolutnym. Wynik badania numerycznego wskazał, że spadek sprawności związany z prześwitem nasila się przy miniaturyzacji.5 Nie jest to pomiar ani uniwersalna korelacja. Jego wartość polega na projekcie porównania, który nie miesza od razu $Re$ z $c/h$.

Cheng, Lu, Zhao, Huang i Zhu zbadali numerycznie wpływ zmian luzu w modelu 1,5-stopniowej sprężarki osiowej przy wysokiej i niskiej liczbie Reynoldsa, stosując trójwymiarowy model RANS.6 Tylko ten zakres twierdzenia został tu zweryfikowany z abstraktu. Badanie uzasadnia rozpatrywanie luzu i Reynoldsa jako dwóch jawnych zmiennych, lecz bez pełnego tekstu nie stanowi dowodu postaci ich sprzężenia ani poprawności addytywnego modelu Studium I.

Należy też rozdzielić wirujące oderwanie (rotating stall) od pompażu (surge). W ujęciu Greitzera wirujące oderwanie jest przestrzennie nieosiowym stanem sprężarki: jedna lub więcej komórek obniżonego przepływu przemieszcza się obwodowo, a średni przepływ przez układ może pozostać dodatni. Pompaż jest niestatecznością całego układu sprężającego, z dużymi oscylacjami strumienia i ciśnienia, niekiedy z odwróceniem przepływu. O tym, który tryb rozwinie się po przekroczeniu granicy, współdecydują objętość odbiorcza, długość przewodów, charakterystyka dławienia i bezwładność całego układu, nie tylko lokalna mapa wirnika.7,8

Te zjawiska mogą następować po sobie i współistnieć, lecz nie są synonimami. Pojedyncza pulsacja ciśnienia, dźwięk albo ślad na mapie bez rozdzielczości przestrzennej i czasowej nie uprawnia do dowolnego nazwania trybu. W tej monografii nie podaje się progów ani procedur wywoływania niestateczności. Dla uprawnionego zespołu każda nieoczekiwana, narastająca oscylacja jest przesłanką do wykonania zatwierdzonej reakcji ochronnej, a nie okazją do samodzielnego „sprawdzenia, czy to surge”.

Studium I — konflikt Macha, Reynoldsa i luzu względnego

Pytanie, wymaganie i konkurencyjne hipotezy — Studium I

Pytanie brzmi: czy zachowanie podobnej prędkości końcówki i tego samego stanu wlotowego wystarcza, aby po zmniejszeniu skali uznać podobną stratę aerodynamiczną? Wymaganie dydaktyczne jest umyślne i nieoperacyjne: syntetyczna konfiguracja może zostać uznana za zgodną z konfiguracją odniesienia tylko wtedy, gdy górna granica przedziału rozszerzonego dla przyrostu umownego współczynnika straty nie przekracza $0{,}020$, czyli dwóch punktów procentowych. Jest to kryterium autora służące pokazaniu reguły decyzji, a nie wymaganie dla prawdziwego silnika.

Rozpatrujemy trzy hipotezy. $H_M$: podobny Mach końcówki dominuje, zatem po zachowaniu $U/a$ zmiana straty jest zaniedbywalna. $H_{R\varepsilon}$: spadek $Re$ i wzrost $c/h$ dają mierzalny przyrost straty mimo podobnego Macha. $H_I$: interakcja Reynoldsa i luzu jest na tyle duża, że prosty model addytywny nie wystarcza do rozstrzygnięcia. Hipotezy nie są twierdzeniami o realnym silniku; określają, co ma odróżnić eksperyment lub model.

Model, źródła i tabela danych

Konfiguracje A, B i C są jawnie syntetyczne i bezwymiarowe. Nie mają przypisanych średnic, liczby łopatek, geometrii kanału, materiału ani paliwa. Skala wynosi odpowiednio $\lambda=1$, $0{,}5$ i $0{,}25$. Wszystkie mają ten sam syntetyczny stan wlotowy oraz $U/a=\mathrm{const}$. Zatem $n/n_A=1/\lambda$ i $Re/Re_A=\lambda$.

Nie zakładamy idealnego skalowania luzu. Syntetyczny model jego ograniczenia ma postać

$$\frac{c}{c_A}=0{,}55+0{,}45\lambda, \qquad \frac{\varepsilon}{\varepsilon_A}=\frac{0{,}55+0{,}45\lambda}{\lambda}.$$

Współczynnik $0{,}55$ reprezentuje umowną część absolutnego luzu, która nie maleje z geometrią, a $0{,}45$ część skalowalną. Nie pochodzi z produkcji żadnego urządzenia. Przesiewowy model przyrostu straty jest zdefiniowany jako

$$\Delta Y_s= a_R\left[\left(\frac{Re}{Re_A}\right)^{-m}-1\right] +a_\varepsilon\left(\frac{\varepsilon}{\varepsilon_A}-1\right),$$

z syntetycznymi współczynnikami autora $a_R=0{,}020$, $m=0{,}20$ i $a_\varepsilon=0{,}012$. $Y_s$ jest umownym, bezwymiarowym współczynnikiem przesiewowym. Równanie nie jest prawem konstrukcyjnym, mapą stopnia ani dopasowaniem do GTM-120. Kierunki zależności są zgodne z rozdzieleniem efektów w badaniu Xianga i współautorów.5 Addytywność jest wyłącznie syntetycznym założeniem autora; praca Chenga i współautorów nie jest używana do jego walidacji ani falsyfikacji.6

Wielkość A B C Status pochodzenia
skala $\lambda$ 1,000 0,500 0,250 dane syntetyczne autora
$U/a$ 1,000 1,000 1,000 dane syntetyczne autora, wartość znormalizowana
$n/n_A$ 1,000 2,000 4,000 wartość wyznaczona pośrednio z $1/\lambda$
$Re/Re_A$ 1,000 0,500 0,250 wartość wyznaczona pośrednio z przyjętej drogi skalowania
$c/c_A$ 1,000 0,775 0,6625 wynik modelu syntetycznego luzu
$\varepsilon/\varepsilon_A$ 1,000 1,550 2,650 wartość wyznaczona pośrednio
$\Delta Y_s$ 0 0,00957 0,02619 wynik modelu syntetycznego straty

W tabeli nie ma pomiaru. Nie wolno jej cytować jako charakterystyki klasy silników. Jej funkcją jest umożliwienie sprawdzenia rachunku i wykazanie, że podobieństwo jednego parametru nie zamyka pozostałych.

Rachunek — Studium I

Dla konfiguracji B prędkość obrotowa rośnie dwukrotnie, ale Reynolds maleje o połowę. Stosunek luzu względnego wynosi

$$\frac{\varepsilon_B}{\varepsilon_A} =\frac{0{,}55+0{,}45\cdot0{,}5}{0{,}5} =1{,}55.$$

Składnik reynoldsowski przyjmuje wartość

$$0{,}020\left(0{,}5^{-0{,}20}-1\right) =0{,}00297,$$

a składnik luzu

$$0{,}012(1{,}55-1)=0{,}00660.$$

Suma to $\Delta Y_{s,B}=0{,}00957$, czyli umownie $0{,}957$ punktu procentowego. Około 69% przyrostu w tym konkretnym modelu pochodzi ze składnika luzowego. Nie oznacza to, że w rzeczywistym stopniu 69% straty powoduje luz; jest to wyłącznie wynik zadanych współczynników.

Dla C liczba obrotów rośnie czterokrotnie, $Re/Re_A=0{,}25$, a

$$\frac{\varepsilon_C}{\varepsilon_A} =\frac{0{,}55+0{,}45\cdot0{,}25}{0{,}25} =2{,}65.$$

Odpowiednie składniki wynoszą

$$\Delta Y_{R,C}=0{,}020(0{,}25^{-0{,}20}-1)=0{,}00639,$$

$$\Delta Y_{\varepsilon,C}=0{,}012(2{,}65-1)=0{,}01980,$$

stąd $\Delta Y_{s,C}=0{,}02619$, czyli $2{,}619$ punktu procentowego. Zachowany $U/a$ nie zapobiegł ani spadkowi Reynoldsa, ani ponad dwukrotnemu wzrostowi luzu względnego.

Przewidywane i rzeczywiste obserwacje

Gdyby $H_M$ wystarczała, A, B i C powinny mieć $\Delta Y_s\approx0$ w granicach niepewności. $H_{R\varepsilon}$ przewiduje monotoniczny wzrost przy malejącym $\lambda$, z rosnącym udziałem luzu. $H_I$ przewiduje, że błąd postaci modelu będzie rósł wraz z oddaleniem od A, a addytywna suma może nie zachować nawet prawidłowej krzywizny trendu.

„Rzeczywiste obserwacje” w tym studium nie istnieją, ponieważ konfiguracje są syntetyczne. Obserwacją rachunkową jest tylko to, że zadeklarowany model zwraca 0,00957 dla B i 0,02619 dla C. Literatura pierwotna wspiera jakościowo kierunek $H_{R\varepsilon}$, ale nie dostarcza walidacji tych współczynników dla konfiguracji A–C. Fakt ten zamyka możliwość podszywania się demonstracji pod eksperyment.

Budżet niepewności i decyzja

Stosujemy liniową propagację GUM. Dla niezależnych wejść

$$u_c^2(\Delta Y_s)=\sum_i \left(\frac{\partial \Delta Y_s}{\partial x_i}\right)^2u^2(x_i)+u_{form}^2.$$

Przyjęte syntetyczne standardowe niepewności typu B wynoszą: $u(a_R)=0{,}004$, $u(m)=0{,}05$, $u(a_\varepsilon)=0{,}003$; dla $\varepsilon/\varepsilon_A$: 0,10 w B i 0,20 w C; rezerwa postaci modelu $u_{form}=0{,}0025$. Skala i stosunek Reynoldsa są w tym ćwiczeniu definicjami, więc nie dodano im niepewności. W realnym badaniu musiałyby ją mieć.2

Źródło wkład standardowy do B wkład standardowy do C Charakter
współczynnik $a_R$ 0,00060 0,00128 syntetyczny typ B
wykładnik $m$ 0,00080 0,00183 syntetyczny typ B
współczynnik $a_\varepsilon$ 0,00165 0,00495 syntetyczny typ B
względny luz 0,00120 0,00240 syntetyczny typ B
postać modelu 0,00250 0,00250 syntetyczny typ B
$u_c$, suma kwadratowa 0,00338 0,00644 wynik propagacji
$U=2u_c$ 0,00675 0,01288 przedział rozszerzony, $k=2$

Dla B przedział wynosi $0{,}00957\pm0{,}00675$, czyli od 0,00282 do 0,01633. Jego górna granica jest niższa niż wymaganie 0,020, więc B spełnia syntetyczne kryterium z marginesem 0,00367. Dla C przedział $0{,}02619\pm0{,}01288$ rozciąga się od 0,01331 do 0,03907 i przecina granicę. Nie wykazano zgodności C. Decyzja jest konserwatywna: nie wolno przenieść ilościowego wyniku z A do C na podstawie zachowanego Macha końcówki. Nie jest to twierdzenie, że rzeczywisty obiekt C na pewno nie spełniłby wymagania; jest to stwierdzenie, że przedstawiony dowód tego nie wykazuje.

Użycie $k=2$ nie gwarantuje automatycznie dokładnie 95% pokrycia. GUM wymaga podania założeń o rozkładach i stopniach swobody; tu $k=2$ jest jawnie dydaktycznym mnożnikiem, a przedział służy regule decyzji, nie certyfikacji.2

Kwestie nierozstrzygnięte i test regresyjny

Nierozstrzygnięte pozostają: interakcja $Re\times\varepsilon$, wpływ chropowatości względnej, przesunięcie punktu pracy, zmiana profilu warstwy przyściennej oraz to, czy podobieństwo $U/a$ rzeczywiście utrzyma lokalne trójkąty prędkości. Ich brak nie jest „małą poprawką”; określa domenę modelu jako przesiewową.

Test regresyjny dla implementacji rachunku używa dodatkowej, syntetycznej konfiguracji D o $\lambda=0{,}75$, niewykorzystanej w tabeli. Z równań wynika $\varepsilon_D/\varepsilon_A=1{,}18333$ i $\Delta Y_{s,D}=0{,}00338$. Każda ponowna implementacja powinna odtworzyć wynik z dokładnością bezwzględną 0,00005. Jest to test kodu i zapisu równań, nie walidacja fizyki. Walidacja fizyczna wymagałaby niezależnych danych dla oznaczonej geometrii i warunków.

Karta dowodowa historycznego cywilnego GTM-120

GTM-120 jest użyteczny jako dobrze opisany, cywilny obiekt badawczy, lecz właśnie dlatego wymaga ochrony przed nadinterpretacją. Nazwa obejmuje konkretną architekturę z jednostopniową sprężarką promieniową, układem dyfuzora i nawrotu, komorą oraz jednostopniową turbiną. Nie jest nazwą ogólnej klasy miniaturowych turbojetów. Poniższa karta rozdziela deklarację katalogową, pomiar, CFD i redukcję pośrednią.

Twierdzenie lub wartość Etykieta pochodzenia Dokładny lokalizator Co rzeczywiście ustanawia Ograniczenie
masa 1,5 kg; maksymalny ciąg 120 N; maksymalna prędkość 120 000 obr/min; deklarowany spręż 2; przepływ powietrza 0,4 kg/s; paliwo 400 ml/min dane publikowane, katalog producenta przytoczony przez autorów Gieras, wyd. 2, s. 55–56 druk / 56–57 PDF, tab. 5.1; Gieras–Stańkowski, s. 69, tab. 1 deklarowaną charakterystykę historycznego GTM-120 nie jest niezależnym pomiarem; nie podano wspólnego punktu i pełnej niepewności wszystkich maksimów
zewnętrzna średnica 110 mm i długość 265 mm dane publikowane Chmielewski–Gieras, s. 48, tab. 1 gabaryt badanego typu nie ujawnia definicji kompletnej instalacji ani geometrii wykonawczej; nie wolno z tego rekonstruować kanałów
sześć prób; warunki otoczenia 17–19°C, 1005–1016 hPa, wilgotność 41–56%; punkty ustalone 40–120 tys. obr/min; 2,5–5 min na przebieg pomiar — protokół serii Chmielewski–Gieras, s. 49 zakres i powtarzanie kampanii stanowiskowej krótki czas i jeden opisany typ stanowiska; stan cieplny zależy od procedury chłodzenia
ciąg, statyczne ciśnienie za sprężarką, strumień powietrza, paliwa, temperatury i emisje przedstawione na wykresach pomiar Chmielewski–Gieras, s. 49–51, rys. 3–12 trendy w oznaczonej kampanii wartości odczytane z wykresu miałyby dodatkową niepewność digitalizacji; artykuł jej nie dodaje
około 90 N w punkcie użytym do porównania z obiegiem pomiar siłomierzem przywołany w pierwotnej analizie Gieras–Stańkowski, s. 74 jeden integralny punkt ciągu stanowiskowego nie jest bezpośrednim pomiarem przepływu przez komorę ani spadku ciśnienia całkowitego
około 20,6% straty ciśnienia całkowitego dla wariantu z założonym przepływem 0,4 kg/s wartość wyznaczona pośrednio z dwóch wyników CFD Modelu I Gieras–Stańkowski, s. 72, model I czułość wyniku konkretnego modelu na katalogowy warunek przepływu nie jest pomiarem; 0,4 kg/s jest daną katalogową
przepływ około 0,25 kg/s wybrany na podstawie trendu dla silników o zbliżonej masie; obieg z założoną stratą 6% dał 88 N wobec około 90 N z pomiaru; około 10% straty uzyskano następnie w CFD wartość scenariuszowa, wynik modelu obiegowego i wynik CFD Gieras–Stańkowski, s. 73–74, model III, rys. 11–14; s. 78, wnioski scenariusz zgodny integralnie z jednym punktem ciągu i wynik zimnego CFD dla opisanej domeny przepływ nie był w tej pracy bezpośrednio zmierzony; 10% zależy od geometrii, przekrojów, RANS i warunku brzegowego
sprawność spalania około 86% przy 80 000 obr/min wartość wyznaczona pośrednio z pomiarów i założeń cieplnych Gieras, s. 93–94 druk / 94–95 PDF wynik bilansu autora dla jednego stanu nie jest bezpośrednim odczytem; zależy m.in. od ciepła właściwego, wartości opałowej, profilu temperatury i strat do ścian
ciśnienie całkowite za sprężarką wyznaczone z modelu, ponieważ sonda nieproporcjonalnie zaburzała mały kanał wartość wyznaczona pośrednio, hybryda pomiar–CFD Gieras, s. 174–175 druk / 175–176 PDF dokumentuje ograniczenie metrologiczne zgodność modelu do około 0,5% na innych wielkościach nie jest niezależną walidacją niezmierzonej stacji

W źródłach pojawiają się różne średnice opisane jako średnica wlotu. Nie należy ich „uśredniać”: odnoszą się do odmiennie nazwanych przekrojów albo różnych sposobów opisu. Jest to przykład problemu semantycznego stacji, nie błąd do naprawienia arbitralną liczbą. W artykule nie wykorzystuje się tych średnic do geometrii wykonawczej.

Kampania Chmielewskiego i Gierasa miała wartościową strukturę: sześć przebiegów, punkty podczas zwiększania i zmniejszania prędkości, tysiące próbek kanałów szybkich i mniej próbek emisji. Autorzy połączyli niepewność pomiaru i odchylenie standardowe sumą kwadratową.3 Nie oznacza to jednak, że wyznaczyli pełny budżet każdego wyniku. Brakujące składniki — profil przestrzenny, wpływ sondy, dryft stanu cieplnego, korelacje i niepewność modelu redukcji — nadal trzeba nazwać.

Katalogowy maksymalny ciąg i maksymalny przepływ nie muszą należeć do tego samego ustalonego punktu, jeżeli dokument tego nie mówi. Wartość „maksymalna” jest cechą tabeli, nie równaniem stanu. Karta dowodowa pozwala wykorzystać dane bez nadawania im wyższego statusu niż mają.

Studium II — czy strata w dyfuzorze, nawrocie i komorze GTM-120 jest prawem skali

Pytanie, wymaganie i hipotezy — Studium II

Pytanie brzmi: czy publikowane około 10% utraty ciśnienia całkowitego można uznać za ilościowy skutek miniaturyzacji, czy jedynie za wynik konkretnej geometrii GTM-120 i modelu? Wymaganie dla transferu jest ostre: ogólna teza o skali wymagałaby co najmniej dwóch skal geometrycznych porównanych przy kontrolowanych grupach podobieństwa, zgodnych definicji stacji oraz niezależnego pomiaru albo zwalidowanego modelu integralnej straty. Jedna geometria i jeden scenariusz RANS nie spełniają tego wymagania.

Konkurują cztery hipotezy. $H_G$: dominuje konkretna geometria dyfuzora, nawrotu i komory. $H_S$: dominuje efekt małej skali, przede wszystkim niższy Reynolds i większy udział warstw przyściennych. $H_B$: pozornie duża strata jest w znacznej mierze skutkiem warunku brzegowego, wyboru stacji i modelu turbulencji. $H_{GSB}$: mechanizmy są sprzężone, a dostępne dane nie identyfikują osobno ich udziałów.

Model i źródła

Pierwotną podstawą jest praca Gierasa i Stańkowskiego z 2012 r. Autorzy modelowali pełną trójwymiarową domenę dyfuzora i komory GTM-120 jako ustalony, ściśliwy przepływ gazu doskonałego. Zastosowali model RANS realizable $k\!-\!\varepsilon$ w Fluent 12.1 oraz siatkę około 1,835 mln elementów tetraedrycznych o charakterystycznym rozmiarze 0,5–0,7 mm.1 To wystarczająco precyzyjny opis domeny dowodu, lecz nie instrukcja odtworzenia komory; nie podaje się tu jej wykonawczej geometrii.

Model I używał katalogowego przepływu 0,4 kg/s. W Modelu III autorzy wybrali około 0,25 kg/s na podstawie zestawienia przepływu z masą innych małych silników; dopiero potem obieg z założoną stratą 6% dał 88 N, wartość bliską około 90 N z pomiaru siłomierzem. Późniejsza praca stanowiskowa Chmielewskiego i Gierasa mierzyła przepływ wlotowy w serii sześciu prób, ale nie jest niezależnym pomiarem rozkładu ciśnienia całkowitego przez wszystkie przekroje domeny z pracy CFD.1,3

NASA-STD-7009B wymaga, by zamierzone użycie modelu, kryteria akceptacji, domena walidacji, źródła danych empirycznych i niezwalidowane aspekty były zapisane, a wynik przekazywany wraz z niepewnością, zastrzeżeniami i ryzykiem użycia.9 Zastosowane tu zamierzone użycie jest węższe niż w oryginalnej publikacji: model może pokazać mechanizm i czułość GTM-120 na przyjęty przepływ, lecz nie może wyznaczyć prawa miniaturyzacji.

Dane i ich status

Pozycja Wartość Status Znaczenie dowodowe
przepływ wariantu I 0,400 kg/s dane publikowane, katalogowe; wejście modelu warunek scenariusza, nie pomiar w domenie CFD
$p_{t,in}$ wariantu I 0,233493 MPa średnia ważona, wynik CFD/Fluent Modelu I stan wejściowy uzyskany w domenie CFD
$p_{t,out}$ wariantu I 0,185405 MPa średnia ważona, wynik CFD/Fluent Modelu I stan wyjściowy uzyskany w domenie CFD
strata wariantu I 20,6% wartość wyznaczona pośrednio z dwóch wyników CFD Modelu I niepewność katalogowego warunku przepływu jest dominująca
przepływ wariantu III około 0,250 kg/s wartość scenariuszowa wybrana z trendu przepływu względem masy małych silników nie jest bezpośrednim pomiarem przepływu w pracy z 2012 r.; obieg dał 88 N wobec około 90 N z pomiaru
ciąg porównawczy około 90 N pomiar integralny ogranicza kombinację wielu parametrów, nie samą stratę komory
założona strata w obiegu wariantu III 6,0% założenie modelu obiegowego służy do wyznaczenia warunków brzegowych; nie jest jego niezależnym wynikiem
strata w Modelu III około 10% wynik zimnego CFD dotyczy straty wewnątrz modelowanej komory; przekrój $y=0{,}01$ m tuż za dyfuzorem jest jej przekrojem wejściowym; nie jest pomiarem ani prawem skali
typowa strata większych komór 4–8% dane publikowane w literaturze porównawczej przywołanej przez autorów punkt odniesienia konfiguracje i definicje nie są dopasowaną serią skalową

Rozróżnienie 6% i około 10% nie jest kosmetyczne. Pierwsza wartość pochodzi z założeń modelu obiegowego między jego stacjami, druga z interpretacji rozkładu CFD w określonej części domeny. Traktowanie ich jako dwóch pomiarów tego samego mezurandu byłoby błędem.

Rachunek — Studium II

Dla wariantu I współczynnik straty wynikający z opublikowanych stanów modelu wynosi

$$\zeta_I=\frac{0{,}233493-0{,}185405}{0{,}233493} =\frac{0{,}048088}{0{,}233493}=0{,}20595.$$

Po zaokrągleniu jest to 20,6%. Ten rachunek jest odtwarzalny, ale obie liczby wejściowe nadal są średnimi ważonymi z CFD/Fluent Modelu I. Nie nabywają statusu pomiaru przez wykonanie dzielenia.

Dla wariantu III opublikowane wartości modelu obiegowego, obliczone po przyjęciu straty 6%, wynoszą $p_{t,in}=0{,}202655$ MPa i $p_{t,out}=0{,}190496$ MPa. Zatem

$$\zeta_{III,cykl}=\frac{0{,}202655-0{,}190496}{0{,}202655} =0{,}059999\approx6{,}0\%.$$

Rachunek odtwarza założenie obiegu, a nie niezależny rezultat modelu. Autorzy uzyskali następnie około 10% utraty w zimnym CFD Modelu III wewnątrz modelowanej komory; przekrój $y=0{,}01$ m tuż za dyfuzorem był przekrojem wejściowym komory, a nie końcem ocenianego odcinka. Różnica czterech punktów procentowych pokazuje, że przed porównaniem trzeba uzgodnić stacje, uśrednienie i zakres domeny. Nie wolno wybrać wartości bliższej oczekiwaniu i nazwać jej „zmierzoną stratą GTM-120”.1

Zmiana wejściowego przepływu z 0,4 do około 0,25 kg/s, czyli o 37,5%, przesunęła raportowaną skalę straty z około 20% w kierunku około 10%. To analiza czułości warunku brzegowego, nie eksperyment skali: geometria i wymiar pozostały te same. Już ten fakt wystarcza, aby odrzucić interpretację „10% jest stałą miniaturowej komory”.

Przewidywane i zaobserwowane skutki

$H_G$ przewiduje skupienie strat w obszarach silnego skrętu, dyfuzji, recyrkulacji i niejednorodnego zasilania tej konkretnej geometrii. $H_S$ przewiduje, że po zmianie skali przy zachowaniu bezwymiarowej geometrii strata zmieni się systematycznie z $Re$ i chropowatością względną. $H_B$ przewiduje dużą czułość na przyjęty przepływ, stacje, turbulencję i cieplny model. $H_{GSB}$ przewiduje, że dane z jednej skali nie pozwolą na identyfikację udziałów.

Rzeczywiste obserwacje z publikacji są następujące. Model wykazał silnie trójwymiarowy przepływ i nierównomierności; zmiana warunku masowego znacząco zmieniła stratę. Dodanie uproszczonego źródła ciepła zmieniło profile termiczne, lecz nie odwróciło głównej topologii przepływu w analizie autorów. Integralny pomiar ciągu był zgodny z obiegiem dopiero po zmianie założonego przepływu. Nie przedstawiono niezależnego, przestrzennie rozdzielonego pomiaru $p_t$ przed i za całą domeną, który bez istotnej blokady potwierdzałby 10%.1

Obserwacje są zgodne z $H_B$ i $H_{GSB}$, nie wykluczają $H_G$, a nie wystarczają do wydzielenia $H_S$. Porównanie z 4–8% w większych komorach jest przesłanką, nie kontrolowanym eksperymentem skali, ponieważ zmieniają się również geometria, obciążenie, stacje i metoda.

Budżet niepewności

Publikacja nie daje danych potrzebnych do obliczenia liczbowego przedziału GUM dla 10%. Pełny budżet należy więc przedstawić jako budżet identyfikowalnych, lecz niekwantyfikowalnych składników, zamiast wymyślić cyfry.

Składnik Typ/status Informacja dostępna Skutek dla decyzji
przepływ wejściowy dominująca niepewność warunku brzegowego dwa scenariusze 0,4 i ok. 0,25 kg/s zmienia wynik straty o skalę porównywalną z samym wynikiem
stacje i uśrednianie niepewność definicji mezurandu 6% w modelu cyklu i ok. 10% w interpretacji CFD wartości nie są bezpośrednio zamienne
model turbulencji i ustalenie przepływu niepewność postaci modelu jeden model RANS; przepływ rzeczywisty może być nieustalony brak liczbowego oszacowania
siatka i zbieżność niepewność numeryczna porównano warianty siatki, rozkłady były zbliżone ogranicza jeden składnik, nie waliduje fizyki
geometria i stan powierzchni niepewność reprezentacji konkretna geometria historyczna; brak serii tolerancji nieznany wpływ na transfer między egzemplarzami
model cieplny i spalanie niepewność modelu fizycznego przepływ zimny i uproszczone źródło ciepła nie ustanawia straty w pełnym reaktywnym stanie
dane empiryczne niepewność walidacji integralny ciąg, ciśnienia statyczne w innych miejscach problem nieidentyfikowalności: wiele kombinacji daje ten sam ciąg
efekt skali brak projektu porównawczego tylko jedna skala nie można wyznaczyć współczynnika skali

Zakres 10–20,6% jest rozpiętością dwóch scenariuszy modelowych, nie przedziałem ufności. Nie znamy rozkładu prawdopodobieństwa ani korelacji składników, więc obliczenie $U=2u_c$ byłoby pozorne. Poprawny wynik metrologiczny brzmi: wartość modelowa zależy silnie od warunku brzegowego, a ilościowa niepewność całkowita nie została wykazana.

Wniosek, sprawy otwarte, decyzja i test regresyjny

Publikowane około 10% jest właściwością konkretnego wyniku modelu konkretnej geometrii GTM-120 przy określonych założeniach, nie udowodnioną stałą skali. Dane nie rozstrzygają ilościowo, jaka część pochodzi z miniaturyzacji, jaka z układu dyfuzora i nawrotu, a jaka z warunków brzegowych i postaci modelu. Sprawa ilościowego rozdziału $H_G$, $H_S$ i $H_B$ pozostaje otwarta.

Decyzja jest jednak zamknięta: wolno użyć wyniku do sformułowania hipotezy dla historycznego GTM-120; nie wolno użyć go jako współczynnika projektowego ani korekty osiągów innego silnika. Brak podstaw do wyniku ilościowego poza domeną. Można podać przedział scenariuszy modelowych, ale nie przedział pomiarowy i nie prawdopodobieństwo.

Test regresyjny przyszłej analizy powinien zamrozić definicje stacji i sposób uśredniania, a następnie przewidzieć wielkości w punkcie lub przebiegu niewykorzystanym do dostrajania. Minimalne kryterium logiczne brzmi: model dopasowany do ciągu w jednym stanie ma odtworzyć jednocześnie niezależny przepływ i różnicę $p_t$ w drugim stanie wraz z ich niepewnościami; zgodność wyłącznie statycznych ciśnień, które były warunkami wejściowymi, nie wystarcza. Dopóki takiego wyniku nie ma, regresja ma status „niewykonana”, a nie „zaliczona przez podobieństwo wykresu”.

Metrologia jest częścią badanego układu

W dużym kanale trzon sondy może być małym zaburzeniem. Po miniaturyzacji jego pole czołowe nie maleje automatycznie jak pole kanału. Jeśli charakterystyczny wymiar sondy $d_p$ jest stały, a wysokość kanału $h$ maleje jak $\lambda$, liniowy stosunek $d_p/h$ rośnie jak $1/\lambda$, a prosta miara blokady powierzchniowej $A_p/A_c$ może rosnąć jak $1/\lambda^2$. Sonda zmienia lokalną prędkość, warstwę przyścienną i ciśnienie, a w kanale z zakrzywieniem może generować separację. Odczyt nie jest wtedy wartością „sprzed włożenia sondy”.

Mezurand trzeba zdefiniować przed doborem instrumentu. „Ciśnienie za sprężarką” może znaczyć lokalne ciśnienie statyczne na ścianie, lokalne ciśnienie całkowite sondy, średnią powierzchniową, średnią masową albo wartość z modelu obiegowego. Każda odpowiada na inne pytanie. Dla bilansu mocy i straty całkowitej potrzebna jest wielkość uśredniona w sposób zgodny ze strumieniem energii; dla obciążenia ściany — lokalne ciśnienie statyczne. Zastąpienie jednego drugim powinno być modelem redukcji z własną niepewnością.

Łańcuch pomiarowy ma co najmniej cztery warstwy:

  1. czujnik — rozdzielczość, szum, histereza, czas odpowiedzi, czułość poprzeczna;
  2. kalibracja i akwizycja — wzorzec, dryft zera i skali, liniowość, synchronizacja, filtracja;
  3. stanowisko i stacja — blokada sondą, przewody impulsowe, przewodzenie cieplne, pole niejednorodne, drgania, siły przewodów i oddziaływanie ścian;
  4. model redukcji — korekcje, własności gazu, uśrednianie, mapy, CFD oraz interpolacja.

Niepewność z powtórzeń jest oceną typu A, lecz nie obejmuje automatycznie stałego błędu kalibracji ani wpływu sondy. Te składniki ocenia się typem B na podstawie świadectw, specyfikacji, wiedzy o metodzie i porównań. GUM podkreśla, że „A” i „B” opisują sposób oceny, nie naturę błędu, a kowariancje trzeba zachować, gdy te same kanały wejściowe występują w kilku wynikach.2

Dla wyniku $y=f(x_1,\ldots,x_N)$ liniowa propagacja z korelacjami ma postać

$$u_c^2(y)= \sum_i\left(\frac{\partial f}{\partial x_i}\right)^2u^2(x_i) +2\sum_{i<j}\frac{\partial f}{\partial x_i} \frac{\partial f}{\partial x_j}u(x_i,x_j).$$

Równanie wykonuje konkretną pracę w Studium III: pokazuje, dlaczego iloraz dwóch ciśnień dziedziczy nie tylko szum czujników, lecz także niepewność modelu stacji i ich wspólne dane wejściowe. Jeżeli nieliniowość lub rozkłady są istotne, linearyzacja wymaga sprawdzenia numerycznego albo propagacji Monte Carlo; samo wypisanie pierwiastka z sumy kwadratów nie wystarcza.

Wynik hybrydowy powinien mieć etykietę wprost w tabeli i na wykresie. Poprawny zapis to na przykład: „$p_{t3}$: wartość wyznaczona pośrednio z CFD ograniczonego pomiarami w stacjach 2 i 4”. Niepoprawny: „zmierzono $p_{t3}$”, jeśli żadna sonda go nie obserwowała. Zasada jest szczególnie ważna przy walidacji: dane użyte do warunków brzegowych lub kalibracji nie są niezależnym testem tej samej wielkości.

Studium III — sonda i wynik hybrydowy dla GTM-120

Pytanie, wymaganie i konkurencyjne hipotezy — Studium III

Pytanie brzmi: czy numerycznie wyznaczone ciśnienie całkowite w stacji, w której sonda nadmiernie ingeruje w przepływ, może służyć do ilościowej decyzji o sprężu i stracie? Wymaganie ma trzy poziomy. Wynik ilościowy jest dopuszczalny, gdy rozszerzona niepewność nie zmienia decyzji względem ustalonego progu. Przedział jest dopuszczalny, gdy wartość środkowa nie rozstrzyga, ale zakres jest obliczalny i użyteczny. Brak podstaw do rozstrzygnięcia zachodzi, gdy dominująca niepewność postaci modelu lub korelacji nie jest oszacowana.

$H_P$: wprowadzenie sondy daje mniejszy błąd niż zastąpienie pomiaru modelem. $H_M$: sonda zaburza kanał nieproporcjonalnie, więc mniej obciążający jest model ograniczony innymi pomiarami. $H_H$: hybryda jest najlepszym dostępnym estymatorem, lecz jej niepewność jest na tyle duża, że pozwala tylko na przedział. $H_0$: brak niezależnego odniesienia nie pozwala ocenić błędu modelu, więc nawet przedział byłby spekulacją.

Model, źródło przypadku i dane

Gieras opisuje stanowisko GTM-120 na s. 169–177 drukowanych (170–178 PDF). Wielkości statyczne i całkowite mierzono w kilku stacjach, ale pomiar ciśnienia całkowitego bezpośrednio za sprężarką uznano za nadmiernie ingerujący; wartość uzyskano numerycznie. Autor podaje zgodność obliczeń z mierzonymi ciśnieniami statycznymi przed i za sprężarką oraz ciśnieniem całkowitym na wylocie na poziomie do około 0,5%.11 Jest to rzeczywisty, udokumentowany przypadek wyboru hybrydy. Nie ma jednak opublikowanego pełnego budżetu niepewności postaci modelu dla nieobserwowanej stacji.

Aby pokazać rachunek bez przypisywania źródłu nieistniejących liczb, używamy poniższego syntetycznego, znormalizowanego i nieoperacyjnego zestawu autora. Ciśnienia podzielono przez umowną wartość odniesienia; nie są to kPa, MPa ani dane GTM-120.

Wielkość Wartość znormalizowana Status Rola
$p_{t2}^*$ 1,000 dane syntetyczne autora, reprezentujące pomiar mianownik sprężu
$p_{t3}^*$ 2,000 dane syntetyczne autora, reprezentujące wynik modelu ograniczonego pomiarami licznik sprężu i mianownik zachowania ciśnienia
$p_{t4}^*$ 1,860 dane syntetyczne autora, reprezentujące pomiar w innej stacji licznik zachowania ciśnienia
zgodność do ok. 0,5% na innych wielkościach dane publikowane o przypadku rzeczywistym Gieras, s. 174–175 druk / 175–176 PDF przesłanka lokalna, nie bezpośredni błąd $p_{t3}$

Rachunek — Studium III

Syntetyczny spręż całkowity sprężarki jest równy

$$\pi_c^*=\frac{p_{t3}^*}{p_{t2}^*} =\frac{2{,}000}{1{,}000}=2{,}000.$$

Syntetyczny współczynnik zachowania ciśnienia między 3 i 4 wynosi

$$\sigma_{3-4}^*=\frac{p_{t4}^*}{p_{t3}^*} =\frac{1{,}860}{2{,}000}=0{,}930,$$

a odpowiadająca mu strata

$$\zeta_{3-4}^*=1-\sigma_{3-4}^*=0{,}070.$$

Te trzy wartości są wynikami redukcji syntetycznych danych. Liczba 0,070 nie pochodzi z książki ani z pomiaru komory. Została dobrana tak, aby unaocznić, że umiarkowana niepewność ciśnienia modelowego może stać się bardzo dużą względną niepewnością małej różnicy $1-\sigma$.

Przewidywania i obserwacje

$H_P$ przewidywałaby, że po wprowadzeniu sondy jej poprawki i blokada będą mniejsze od niepewności modelu. Dokumentowany przypadek mówi przeciwnie na poziomie decyzji eksperymentalnej: sonda była na tyle duża względem kanału, że odstąpiono od bezpośredniego pomiaru $p_{t3}$. Nie znamy jednak kontrfaktycznego, niezaburzonego $p_{t3}$, zatem nie można liczbowo stwierdzić, o ile sonda byłaby gorsza.

$H_M$ przewiduje zgodność modelu w stacjach obserwowanych i możliwość interpolacji do stacji nieobserwowanej. Rzeczywistą obserwacją publikowaną jest zgodność do około 0,5% dla wybranych statycznych ciśnień i całkowitego ciśnienia na wylocie.11 To wspiera spójność, lecz nie wystarcza do przyjęcia 0,5% jako niepewności $p_{t3}$: inny mezurand, inna stacja i możliwe wspólne dane kalibracyjne ograniczają niezależność walidacji.

$H_H$ przewiduje, że po dodaniu niepewności modelu wynik 7% zamieni się w szeroki przedział. Poniższy syntetyczny budżet potwierdza tę własność matematyczną. $H_0$ pozostaje właściwą oceną rzeczywistej publikacji w zakresie, w którym nie podano pełnej niepewności modelowej: nie wolno podmienić braku danych budżetem autora.

Pełny syntetyczny budżet niepewności

Wartości w tabeli są względnymi standardowymi niepewnościami typu B, dobranymi przez autora. Osobno pokazano czujnik, kalibrację, stanowisko i model.

Wejście Czujnik / dane wejściowe Kalibracja Stanowisko / przestrzeń Numeryka Postać modelu Reprezentacja konfiguracji Złożona $u_r$
$p_{t2}^*$, pomiar 0,2% 0,3% 0,4% — — — 0,539%
$p_{t3}^*$, model hybrydowy 0,3% 0,4% 0,7% 1,0% 3,0% 1,0% 3,426%
$p_{t4}^*$, pomiar 0,3% 0,4% 1,2% — — — 1,300%

W tym wyłącznie syntetycznym modelu przyjęto zerowe korelacje między wszystkimi wyszczególnionymi wejściami. Składniki z tabeli złożono więc pierwiastkiem z sumy kwadratów. Do każdego ilorazu dodano ponadto osobny, niezależny syntetyczny składnik niezgodności modelu o względnej niepewności standardowej 1,0%; nie jest on kowariancją ani zastępstwem macierzy korelacji. W rzeczywistym badaniu przy nieznanych korelacjach nie wolno obliczać przedziału GUM metodą RSS: decyzją pozostaje brak podstaw, dopóki nie powstanie uzasadniona macierz korelacji i nie zostaną policzone człony mieszane. Przy zadanych tu zerowych korelacjach względna standardowa niepewność sprężu wynosi

$$u_r(\pi_c^*)= \sqrt{0{,}539^2+3{,}426^2+1{,}000^2}\%=3{,}6097\%.$$

Stąd $u(\pi_c^*)=2{,}000\cdot0{,}036097=0{,}07219$, a dla $k=2$:

$$\pi_c^*=2{,}000,\qquad U=0{,}14439, \qquad I_k=[1{,}85561;\,2{,}14439].$$

Dla zachowania ciśnienia

$$u_r(\sigma_{3-4}^*)= \sqrt{1{,}300^2+3{,}426^2+1{,}000^2}\%=3{,}799\%,$$

więc $u(\sigma^*)=0{,}930\cdot0{,}03799=0{,}03533$ i $U=0{,}07066$. Ponieważ $\zeta=1-\sigma$, bezwzględna niepewność jest taka sama:

$$\zeta_{3-4}^*=0{,}07000,\qquad U(\sigma^*)=U(\zeta^*)=0{,}07066.$$

Formalny przedział liniowy dla straty wynosi $[-0{,}00066;\,0{,}14066]$; po nałożeniu fizycznego warunku dla pasywnej straty otrzymujemy około $[0;\,0{,}1407]$. To obcięcie zmienia rozkład, dlatego nie wolno przypisać mu automatycznie poziomu 95%. Najważniejszy skutek jest decyzyjny: wartość środkowa 7% nie rozróżnia wiarygodnie 6% od 10%, chociaż ciśnienia wejściowe wyglądają precyzyjnie.

Wniosek, kwestie otwarte, decyzja i test regresyjny

Wynik hybrydowy może być lepszy od destrukcyjnie zaburzonego pomiaru, ale musi pozostać wynikiem hybrydowym. Jego wiarygodność jest właściwością zamierzonego użycia: model wystarczający do odtworzenia trendu statycznego ciśnienia może nie być wystarczający do małej różnicy dwóch ciśnień całkowitych.

Dla syntetycznego przykładu decyzja ma dwa poziomy. Spręż można raportować liczbowo jako 2,000 z rozszerzonym przedziałem 1,8556–2,1444, jeśli ten zakres odpowiada potrzebie dydaktycznej. Dla straty wolno raportować jedynie szeroki przedział około 0–14,07%; brak podstaw do rozstrzygnięcia, czy wynosi ona 6%, 7% czy 10%. Dla rzeczywistego GTM-120 budżet syntetyczny nie obowiązuje: z publikacji wynika status hybrydowy, lecz bez pełnego budżetu i macierzy korelacji nie wynika liczbowy przedział.

Otwarte pozostają: błąd postaci modelu w nieobserwowanej stacji, kowariancja wspólnych warunków brzegowych, przestrzenne uśrednienie $p_t$, zmienność między przebiegami i wpływ stanu cieplnego. Decyzja dla realnego przypadku brzmi: nie wykorzystywać samego $p_{t3}$ do wąskiego kryterium ilościowego bez dodatkowego odniesienia.

Test regresyjny musi użyć pomiaru, punktu albo stacji niewłączonej do kalibracji. Przed obliczeniem zapisuje się przewidywany przedział reszt. Test jest zaliczony tylko wtedy, gdy niezależny wynik mieści się w tym przedziale i nie ujawnia systematycznego trendu z prędkością lub temperaturą. Ponowne porównanie z tym samym ciśnieniem statycznym, które wyznaczało warunek brzegowy, jest testem odtwarzania wejścia, nie walidacją $p_{t3}$.

Bilans cieplny, czas spalania i znaczenie sprawności wyznaczonej pośrednio

Wartość około 86% podana dla GTM-120 przez Gierasa jest dobrym przykładem wyniku, który wygląda jak pojedyncza cecha komponentu, lecz w istocie jest bilansem. Dla stanu opisanego przy 80 000 obr/min wykorzystano publikowane wielkości pomiarowe: $T_2=345$ K, $T_3=897$ K, zużycie paliwa 237 ml/min przeliczone na około 3,25 g/s oraz przepływ gazu około 203 g/s, wraz z przyjętymi ciepłami właściwymi i wartością opałową 44 MJ/kg.11 Status wyniku to wartość wyznaczona pośrednio przez autora książki, nie bezpośredni pomiar sprawności płomienia.

W uproszczonej redukcji książkowej wielkość wynikowa ma postać

$$\eta_{b,\mathrm{red}}= \frac{\dot m_g c_{pg}T_{t4}-\dot m_a c_{pa}T_{t3}} {\dot m_fH_u}.$$

Indeks „red” odróżnia tę wielkość zredukowaną od zdefiniowanego wcześniej chemicznego stopnia uwolnienia energii $\eta_b$. Ponieważ redukcja nie wydziela osobno $\dot Q_w$, jej wartość może łączyć niepełne chemiczne uwolnienie energii ze stratą cieplną i pozostałymi uproszczeniami bilansu; nie należy jej interpretować jako czystej sprawności chemicznej. Symbole stacji w książce nie są automatycznie tożsame z numeracją tej monografii; równanie pokazuje sens redukcji, a nie przenosi bez sprawdzenia nazw przekrojów. Jeżeli $\dot m_a=\dot m_g-\dot m_f$, to dla przytoczonych przybliżeń różnica entalpii daje wynik rzędu 0,85–0,86. Rachunek zamyka arytmetykę, lecz nie zamyka metrologii: przestrzenny profil temperatury, zmienność $c_p(T)$, gęstość paliwa, $H_u$, utrata ciepła i niepełna jednoczesność kanałów tworzą budżet.11,12

W późniejszej pracy stanowiskowej autorzy sami zastrzegli jakościowy charakter sprawności obiegu ze względu na silnie niejednorodną temperaturę za komorą.3 Jedna termopara w profilu może ważyć zimną lub gorącą smugę inaczej niż średnia entalpijna. Błąd temperatury przechodzi przez różnicę dwóch dużych strumieni entalpii, dlatego względna niepewność $\eta_{b,\mathrm{red}}$ może być znacznie większa niż względna niepewność samego czujnika.

Miniaturyzacja wpływa na ten bilans dwiema drogami naraz. Wzrost $\chi_Q$ odbiera energię gazowi, a spadek $Da$ może zwiększyć udział niespalonych produktów. Te mechanizmy mogą dać podobny spadek temperatury wyjściowej. Bez emisji, bilansu cieplnego ściany i informacji o rozkładzie temperatury nie są identyfikowalne z jednego $T_{t4}$. Z kolei dodanie emisji nie rozwiązuje problemu automatycznie: miejsce poboru, czas transportu próbki i odpowiedź analizatora muszą należeć do tej samej osi czasu.

Skalowanie czasu wymaga również oddzielenia trzech reżimów. W stanie zimnego przepływu bada się straty aerodynamiczne bez reakcji. W przejściu rozgrzewania zmieniają się temperatura ścian, luzy i właściwości gazu. W quasi-ustalonym stanie gorącym pojemność cieplna może już nie dominować, ale rozkład reakcji nadal jest nieustalony lokalnie. Dane z tych reżimów nie powinny być łączone w jeden punkt tylko dlatego, że wskaźnik prędkości obrotowej jest podobny.

Konsekwencją szkoleniową nie jest nakaz „dodać poprawkę cieplną”, lecz pytanie identyfikacyjne: czy obserwowany spadek wyniku można przypisać aerodynamice, reakcji czy wymianie ciepła przy dostępnych kanałach? Jeśli dwie hipotezy przewidują ten sam odczyt, potrzebny jest inny mezurand lub należy pozostawić przyczynę nierozstrzygniętą.

Trzy historyczne obiekty, trzy odrębne domeny

GTM-120

GTM-120 pozostaje cywilnym, historycznym obiektem głównych studiów. Jego wartość wynika z połączenia opisu architektury, CFD oraz kampanii stanowiskowej. Granica transferu obejmuje jednak tylko udokumentowaną konfigurację i wskazane punkty. Nawet dwa egzemplarze o tej samej nazwie mogą różnić się stanem powierzchni, zużyciem, instalacją pomiarową i wersją podzespołów. Bez identyfikacji konfiguracji nie wiadomo, czy powtórzenie porównuje model z tym samym systemem.

Wnioskiem z GTM-120 nie jest „komora odwrócona traci 10%”, lecz: w opisanym modelu dyfuzora/nawrotu/komory strata była silnie zależna od przyjętego przepływu, a problem sondy wymusił wynik hybrydowy. Są to dwa dowody tezy o niezgodnych podobieństwach — aerodynamicznym i metrologicznym — z ograniczoną domeną.

SO-3

SO-3 pojawia się u Gierasa jako odrębny, historyczny przypadek analizy dynamicznej układu paliwowego. Portret fazowy zależności $\mathrm{d}n/\mathrm{d}t$ od $n$, omawiany na s. 31–33 drukowanych (32–34 PDF), dotyczy SO-3.11 Nie jest mapą tarcia małego wirnika, wynikiem GTM-120 ani uniwersalnym algorytmem diagnostycznym.

Tego przykładu wolno użyć do ogólnej lekcji: dynamika przyspieszania zawiera informację o całym bilansie momentów i sterowaniu, a nie tylko o jednej stracie. Nie wolno przenosić kształtu portretu, progów, nastaw ani interpretacji usterki do innego silnika bez dokumentacji jego dokładnej wersji. Równanie energii wirnika wyjaśnia dlaczego: taki sam $\mathrm{d}n/\mathrm{d}t$ może wynikać z innej mocy turbiny, sprężarki, odbiorników lub strat mechanicznych.

Silnik z NASA TM X-3463

Silnik z NASA TM X-3463 był inną maszyną: historycznym, tanim eksperymentalnym turbojetem o czterostopniowej sprężarce osiowej, pierścieniowej komorze i jednostopniowej turbinie. Raport podaje około 29 cm średnicy, 96,5 cm długości, 59 kg masy i badanie do około 37 000 obr/min.4 To dane publikowane i pomiary tego programu, nie dane GTM-120.

Raport jest przydatny, ponieważ eksperymentalnie analizuje wpływ wskaźnika Reynoldsa na kilka parametrów i wykazuje niejednakową zdolność do uogólnienia przepływu, sprawności i sprężu. Nie wolno jednak przenieść jego liczbowej granicy wskaźnika Reynoldsa do sprężarki promieniowej ani do innej instalacji. Architektura osiowa, liczba stopni, geometria, stanowisko i definicja wskaźnika są inne. Traktowanie raportu jako „większego odpowiednika GTM” byłoby dokładnie tym niekontrolowanym transferem, przed którym ostrzega teza.

Granica transferu: od rdzenia do lotu i do współczesnego napędu wojskowego

Transfer jest dozwolony tylko wtedy, gdy zachowane lub skorygowane zostały parametry istotne dla konkretnego wniosku. Nie trzeba zachować każdej liczby podobieństwa, jeśli wrażliwość wyniku na daną niezgodność jest wykazana jako mała. Trzeba jednak tę małość udowodnić w domenie, a nie założyć.

Przejście z rdzenia stanowiskowego do kompletnego zespołu napędowego dodaje spadki instalacji paliwowej i elektrycznej, masę i pobory układów pomocniczych, rzeczywisty wlot oraz warunki dyszy. Przejście do płatowca dodaje odzysk ciśnienia, zniekształcenie obwodowe i czasowe na wlocie, opór pędu strumienia wlotowego, upusty, oddziaływanie strugi z kadłubem, obciążenia i środowisko. Równanie ciągu zmienia się nie dlatego, że fizyka pędu jest inna, lecz dlatego, że granica kontrolna przecina inne strumienie i siły.

Transfer Minimalny dowód Co można przenieść Czego nie można
punkt GTM-120 → inny punkt tego samego egzemplarza mapa lub model przejścia zwalidowany poza punktem strojenia, wspólne stacje i stan cieplny trend z przedziałem w domenie walidacji maksimum katalogowego jako stałą
GTM-120 → inny mały turbojet seria skalowa lub niezależne dane oraz zgodność $M$, $Re$, $c/h$, cieplna i czasowa istotna dla wyniku mechanizm jakościowy jako hipoteza 10% straty, 86% sprawności albo reakcję dynamiczną
rdzeń statyczny → zespół lotny model wlotu, dyszy, poborów i objętości kontrolnej, zwalidowany w warunkach instalacji bilans z jawnymi korektami ciąg statyczny jako ciąg w locie
obiekt historyczny/cywilny → współczesny napęd wojskowy dokumentacja dokładnej wersji wojskowej, autoryzowane dane konfiguracji i zatwierdzony proces kwalifikacji tylko twierdzenia potwierdzone tą dokumentacją geometria, tolerancje, mapy paliwowe, progi ochronne i osiągi przez analogię

Podobna masa, ciąg, średnica lub prędkość obrotowa nie tworzą domeny walidacji. Dwa silniki o podobnym ciągu mogą mieć inny strumień masy i prędkość strugi; podobny gabaryt może ukrywać inną sprężarkę i komorę; podobne $n$ przy innym promieniu oznacza inną prędkość końcówki. Nazwa marketingowa rodziny także nie zastępuje numeru wersji i konfiguracji.

Dla współczesnego napędu wojskowego obowiązuje zasada zerowego transferu ilościowego bez dokumentacji dokładnej wersji. Niniejsza monografia może pomóc sformułować pytania o Reynolds, luz, temperaturę, czas i pomiar. Nie może dostarczyć brakującej mapy ani nastawy. Jeżeli dokument wersji nie definiuje granicy masy lub ciągu, odpowiedź brzmi „nie wiadomo”, a nie „jak w GTM-120”.

Wynik modelu lotnego musi mieć osobne zamierzone użycie. Model zwalidowany dla ustalonego przepływu rdzenia nie jest przez to zwalidowany dla zniekształconego wlotu, przejścia dynamicznego ani oporu instalacyjnego. NASA-STD-7009B nazywa tę granicę domeną walidacji i wymaga jawnego ostrzeżenia, gdy założenia lub limity modelu są naruszone.9

CFD, stanowisko, łożyska i sterowanie — tylko kontekst konieczny

CFD rozdziela lokalne pola, których nie daje pojedyncza sonda, ale wynik pozostaje zależny od siatki, turbulencji, przejścia laminarnego, ścian, warunków brzegowych i uśrednienia. Zbieżność siatkowa ogranicza błąd dyskretyzacji; nie dowodzi poprawności modelu fizycznego. Zgodność jednego integralnego ciągu nie identyfikuje jednocześnie przepływu, straty komory i sprawności komponentów. Szerszy warsztat budowania i kwalifikacji modeli opisano w artykule Modelowanie i kwalifikacja przepływów wewnętrznych.

Hamownia również nie jest „rzeczywistością bez modelu”. Tarowanie stanowiska, siły przewodów, temperatura otoczenia, odzysk wlotu, stacje, częstotliwość próbkowania i stan cieplny definiują wynik. Droga od obliczenia do bezpiecznego programu dowodowego została rozwinięta osobno: Od modelu i CFD do próby stanowiskowej małego napędu. Tutaj potrzebny jest tylko wniosek: próba waliduje oznaczoną konfigurację i mezurand w określonym zakresie, nie wszystkie możliwe użycia modelu.

Łożyska, smarowanie i uszczelnienia wpływają na $P_{\mathrm{loss}}$, ruch wirnika, luz oraz bilans cieplny. Uszczelnienie labiryntowe jest uszczelnieniem bezkontaktowym w prawidłowym stanie pracy; ogranicza przeciek przez sekwencję dławień i komór, a nie przez „tarcie zębów”. Kontakt może wystąpić jako niepożądane ocieranie, lecz nie jest zasadą działania. Szczegółowy kontekst trybologiczny znajduje się w Trybologii, łożyskach i smarowaniu mechanizmów.

Sterowanie zamyka pętlę między paliwem, prędkością, temperaturą i punktem sprężarki. To, że po lokalnej stracie układ wraca do podobnego $n$, nie znaczy, że wrócił do tego samego stanu: mógł zwiększyć paliwo i zmienić temperaturę oraz TSFC. Bez zapisanej logiki sterowania i ograniczeń nie wolno przypisać różnicy wyłącznie aerodynamice. Niniejszy tekst nie publikuje nastaw ani progów ochronnych.

Wnioski szkoleniowe prowadzące do decyzji

Zadaniem żołnierza, analityka lub inżyniera nie jest rozpoznanie wszystkich części po wyglądzie, lecz ustalenie, czy dostępny dowód pozwala podjąć konkretną decyzję. Proces zaczyna się od zdania decyzyjnego, na przykład: „czy opublikowany ciąg może być użyty do oceny tej konfiguracji lotnej?” albo „czy różnica straty między przebiegami przekracza niepewność?”. Dopiero potem wybiera się mezurand, granicę systemu i źródła.

Pierwszy filtr dotyczy tożsamości obiektu. Jeśli nie potwierdzono wersji, konfiguracji i granicy instalacji, dane analogiczne mogą służyć tylko do tworzenia hipotez. Drugi filtr dotyczy pochodzenia: katalog, pomiar, model, redukcja pośrednia i dane syntetyczne nie są wymienne. Trzeci dotyczy niezgodności podobieństwa: trzeba wskazać tę, która może zmienić decyzję, i oszacować jej skutek. Czwarty dotyczy niepewności i marginesu.

Praktyczna reguła decyzji ma trzy wyniki:

  • wynik ilościowy — mezurand, stacje i warunki są zgodne, a przedział niepewności pozostaje po jednej stronie kryterium;
  • wynik przedziałowy — środek jest niepewny, lecz uczciwie wyznaczony zakres nadal ogranicza możliwe decyzje;
  • brak podstaw — brak dominującego składnika budżetu, wersji, stacji, niezależnej walidacji albo zgodności skali uniemożliwia nawet przedział.

„Brak podstaw” jest wynikiem analizy, nie porażką. W Studium II chroni przed zamianą scenariusza CFD w prawo skali. W Studium III chroni przed nazwaniem wartości hybrydowej pomiarem. W zastosowaniu wojskowym chroni przed użyciem danych cywilnego obiektu jako cechy nieznanej wersji.

Jeżeli badany parametr jest różnicą dwóch dużych wielkości, jak strata $1-p_{t4}/p_{t3}$, budżet powinien być oceniony bezpośrednio dla tej różnicy. Dobra dokładność pojedynczych kanałów nie gwarantuje dobrej względnej dokładności małej straty. Jeżeli decyzja dotyczy przejścia dynamicznego, średnia ustalona nie zastępuje trajektorii, a wirujące oderwanie nie zastępuje pompażu.

W warunkach eksploatacyjnych jakikolwiek nieoczekiwany wzrost drgań, nietypowy dźwięk, dym, zapach paliwa, wyciek, widoczne uszkodzenie lub niekontrolowana zmiana wskazań jest wyłącznie przesłanką do przerwania działań przez uprawniony personel zgodnie z zatwierdzoną procedurą. Lista ta nie jest instrukcją samodzielnej inspekcji, zbliżania się ani diagnozowania pracującego urządzenia.

Bezpieczeństwo prób i postępowanie EOD

OSTRZEŻENIE. Wirnik może ulec wysokoenergetycznej fragmentacji. Paliwo, gorące gazy, zapłon i hałas są niebezpieczne. Samodzielna budowa i próby bez zatwierdzonego stanowiska oraz uprawnionego zespołu są niedopuszczalne.

Energia obrotowa nie znika natychmiast po odcięciu polecenia. Odłamki mogą przebić lekką osłonę, gorące części pozostają niebezpieczne po zatrzymaniu, paliwo może tworzyć palną atmosferę, a hałas uszkadzać słuch. Ocena osłony i strefy powstrzymującej fragmenty (ang. containment), strefy wyłączenia, zdalnej obsługi, ochrony pożarowej, wentylacji, izolacji energii i reakcji awaryjnej należy do zatwierdzonej organizacji dla konkretnego stanowiska. Ten artykuł nie stanowi procedury próby.

Uprawniony personel nie powinien traktować widocznego zatrzymania jako potwierdzenia stanu bezpiecznego. Dostęp może nastąpić wyłącznie według procedury obiektu, po potwierdzeniu wymaganych stanów energii, paliwa, ciśnienia i temperatury. Nietypowych symptomów nie bada się przez ponowne zasilenie, zwiększenie obciążenia ani zbliżenie do płaszczyzny wirnika.

Znalezionego lub nieznanego napędu albo fragmentu uzbrojenia nie dotykać, nie przemieszczać, nie zasilać, nie otwierać ani nie demontować. Należy oddalić się tą samą drogą, ostrzec innych i ograniczyć dostęp, oznaczyć lub zapamiętać miejsce z bezpiecznej odległości, zawiadomić 112, Policję albo właściwe służby i wykonywać ich polecenia.

Nieznany fragment może zawierać materiał energetyczny, źródło zasilania, ciśnienie resztkowe, gorące lub toksyczne pozostałości i urządzenia inicjujące. Brak płomienia, dźwięku albo widocznego przewodu niczego nie rozstrzyga. Postępowanie EOD nie jest rozwinięciem diagnostyki silnika; ma pierwszeństwo przed ciekawością techniczną i zabezpieczaniem materiału do analizy.

Synteza: jak zamknąć wniosek o miniaturyzacji

Miniaturyzacja nie tworzy jednej kary sprawności. Tworzy zbiór niezgodności, które mogą dominować w różnych miejscach i różnych mezurandach. Zachowanie $U/a$ podtrzymuje część podobieństwa ściśliwego, ale podnosi liczbę obrotów i nie zachowuje Reynoldsa. Nieskalujący się luz zwiększa $c/h$; nieskalująca się sonda zwiększa blokadę; skrócenie czasu przepływu zmienia $Da$ i stosunek czasu czujnika do procesu; bilans cieplny zależy od pełnego oporu i czasu ściany, nie od samego $A/V$.

Skutek lokalny trafia do osiągu dopiero przez dopasowanie całego rdzenia. Strata ciśnienia zmienia warunki turbiny i dyszy, reakcję sterowania, ciąg oraz TSFC. Bez map, warunków i wspólnej objętości kontrolnej nie da się nadać jej uniwersalnego przelicznika. W próbie lotnej dochodzą opór pędu strumienia wlotowego, odzysk i zniekształcenie wlotu oraz opór instalacyjny.

Trzy studia ustanawiają trzy różne rodzaje decyzji. Studium I pokazuje na danych syntetycznych, że podobny Mach nie wystarcza, a margines musi uwzględniać niepewność. Studium II ogranicza około 10% straty do modelu konkretnej geometrii GTM-120 i pozostawia przyczynę skali otwartą. Studium III pokazuje, że modelowe ciśnienie w nieobserwowanej stacji może być użyteczne, ale ma status hybrydowy; dla małej różnicy jego niepewność może odebrać zdolność rozstrzygnięcia.

Ostateczna reguła brzmi: nazwij obiekt i granicę, zdefiniuj drogę skalowania, wskaż dominującą niezgodność, zmierz jej skutek, przeprowadź niepewność przez wielkość decyzyjną i zapisz granicę transferu. Jeżeli któregoś ogniwa brakuje, należy obniżyć wynik z liczby do przedziału albo z przedziału do braku podstaw. Nie wolno uzupełniać brakującego dowodu podobieństwem nazwy, gabarytu ani przeznaczenia.

Bibliografia

Artykuł opiera się na oficjalnych raportach państwowego laboratorium oraz technicznej literaturze pierwotnej, normatywnej i monograficznej. Agregatory i Wikipedia nie stanowią dowodu dla przedstawionych twierdzeń.