Modelowanie i kwalifikacja przepływów wewnętrznych

Teza: wiarygodność modelu miniaturowego silnika turboodrzutowego nie wynika z dokładnego odwzorowania pojedynczego komponentu ani jednego punktu projektowego. Powstaje dopiero wtedy, gdy model jednocześnie zamyka bilanse masy, energii i momentu wału oraz wykazuje — dla jednoznacznej konfiguracji, przy jawnym statusie danych i rachunku niepewności — zdolność przewidywania zachowania całej drogi gazowej w stanach ustalonych, przejściowych, zdegradowanych i po zabudowie.

Sprężarka nie wybiera swojego punktu pracy. Punkt ten jest skutkiem zdolności przepływowej wlotu, strat w dyfuzorze i komorze, temperatury przed turbiną, charakterystyki turbiny, pola dyszy, warunków otoczenia oraz działania sterowania. Zmiana jednego z tych elementów przemieszcza stan wszystkich pozostałych. Mapa komponentu jest więc warunkiem wejściowym do rozwiązania układu sprzężonych równań, a nie samodzielnym dowodem osiągów lub stateczności silnika.

Tekst dotyczy fizyki, modelowania, metrologii i kwalifikacji dowodu. Nie podaje geometrii wykonawczej współczesnego napędu wojskowego, jego map komponentowych, harmonogramów paliwa, progów zabezpieczeń, nastaw rozruchowych ani procedur wywoływania niestateczności. Szczegółowe skutki samej miniaturyzacji oraz ogólny proces przechodzenia od CFD do próby stanowiskowej omówiono w odrębnych monografiach; tutaj przedmiotem jest węższy problem: jak z rozproszonych charakterystyk komponentów utworzyć sprawdzalny argument o zachowaniu kompletnej drogi gazowej.

Ostrzeżenie EOD i stanowiskowe: nieznanego napędu, jego wirnika, przewodów ani osprzętu nie wolno zasilać, obracać, przedmuchiwać, napełniać paliwem, rozłączać ani uruchamiać „na próbę”. Układ może zawierać paliwo, ciśnienie, gorące pozostałości, czynny zapłonnik oraz energię zmagazynowaną w wirniku. Należy odizolować miejsce i powiadomić właściwe służby. Identyfikacja, pozbawienie energii, badanie i demontaż należą do uprawnionego personelu, a w przypadku nieznanego uzbrojenia — do EOD.

Granica układu, stacje i status danych

Analizę trzeba rozpocząć nie od wyboru solvera, lecz od narysowania granicy układu. Dla jednowałowego silnika turboodrzutowego ogólny łańcuch ma postać:

otoczenie
   │
   ▼
wlot → sprężarka → dyfuzor/kanał nawrotny → komora spalania
                                                   │
                                                   ▼
                              dysza ← kanał wylotowy ← turbina
                                │
                                ▼
                              struga

Do głównego strumienia dochodzą paliwo, upusty, powietrze chłodzące, przecieki uszczelnień i ewentualne przepływy pomocnicze. Część z nich opuszcza granicę, część omija wybrany komponent, a część wraca do głównej drogi w innej stacji. Zapis „przepływ silnika wynosi $0{,}25 \mathrm{kg/s}$” jest niewystarczający, dopóki nie wiadomo, czy chodzi o powietrze na wlocie, strumień przez wirnik sprężarki, gaz przed turbiną czy masę wypływającą przez dyszę.

W tej monografii przyjęto umowne stacje:

Stacja	Lokalizacja
0	niezaburzone otoczenie przed zespołem napędowym
2	płaszczyzna wejściowa sprężarki
3	płaszczyzna wyjściowa sprężarki
4	wylot komory i wlot turbiny
5	wylot turbiny
8	charakterystyczny przekrój dyszy, na przykład gardziel
9	płaszczyzna wylotowa

Numeracja nie jest normą uniwersalną. W innym źródle stacja 3 może leżeć przed kanałem nawrotnym, za dyfuzorem albo przed komorą. Przed porównaniem dwóch zbiorów danych należy odwzorować fizyczne przekroje, a nie tylko zgodność indeksów.

Indeks $t$ oznacza wielkość całkowitą, czyli stagnacyjną. Brak indeksu oznacza wielkość statyczną. Każda karta stacji powinna zawierać co najmniej:

fizyczne położenie i pole przekroju;
układ osi oraz dodatni kierunek przepływu;
ciśnienie i temperaturę statyczną oraz całkowitą, jeśli są potrzebne;
ciśnienie i temperaturę całkowitą;
strumień masy i skład gazu;
profile przestrzenne;
metodę uśredniania;
konfigurację sond i ich blokadę przepływu;
pasmo pomiarowe i wspólną podstawę czasu;
niepewność oraz flagi jakości danych;
identyfikator geometrii, egzemplarza i kampanii.

Równie ważny jest status każdej liczby. W dalszych studiach stosuje się pięć etykiet:

Etykieta	Znaczenie
pomiar	wynik łańcucha pomiarowego o opisanej spójności i niepewności
dane publikowane	wartość z katalogu, raportu albo literatury, bez twierdzenia o niezależnym potwierdzeniu
wynik modelu	wyjście zależne od równań, map, geometrii, solvera i warunków brzegowych
wartość wyznaczona pośrednio	rezultat redukcji pomiarów lub połączenia pomiaru z modelem
dane syntetyczne autora	liczby dobrane do demonstracji metody, nieopisujące istniejącego napędu

Zgodność modelu z pomiarem nie zamienia wyniku modelu w pomiar. Dane użyte do narzucenia warunku brzegowego albo kalibracji nie są niezależną walidacją tego samego wyniku.

Wielkości statyczne, całkowite i uśrednianie przekroju

W przepływie jednowymiarowym łatwo przypisać przekrojowi pojedyncze ciśnienie, temperaturę i prędkość. Rzeczywisty kanał zawiera warstwy przyścienne, zawirowanie, ślady łopatek, separację, gradienty promieniowe i obwodowe oraz lokalne gorące pasma. „Średnia temperatura przed turbiną” nie jest wobec tego wielkością samowyjaśniającą się.

Dla gazu doskonałego, przepływu adiabatycznego i lokalnego stanu o liczbie Macha $M$ relacje między wielkościami statycznymi i całkowitymi mają postać

$$\frac{T_t}{T}=1+\frac{\gamma-1}{2}M^2,$$

$$\frac{p_t}{p} = \left(1+\frac{\gamma-1}{2}M^2\right)^{\frac{\gamma}{\gamma-1}}.$$

Relacji tych nie wolno bez kontroli stosować przez obszar spalania, falę uderzeniową, silną wymianę ciepła albo przepływ o zmiennym składzie. Ciśnienie całkowite jest szczególnie użyteczne jako miara nieodwracalnych strat: w adiabatycznym przepływie bez pracy użytecznej jego spadek sygnalizuje wytwarzanie entropii.

Dla wielkości $\phi$ średnia powierzchniowa wynosi

$$\overline{\phi}_A = \frac{1}{A}\int_A \phi\,\mathrm dA,$$

natomiast średnia ważona strumieniem masy

$$\overline{\phi}_m = \frac{\int_A \rho V_n\phi\,\mathrm dA} {\int_A \rho V_n\,\mathrm dA}.$$

Druga definicja zwykle lepiej odpowiada bilansom transportowanej masy lub energii, ale także nie jest automatycznie poprawna dla każdego zadania. Do bilansu energii właściwe może być uśrednienie zachowujące strumień entalpii:

$$\overline{h_t} = \frac{\int_A \rho V_n h_t\,\mathrm dA} {\int_A \rho V_n\,\mathrm dA}.$$

Do rachunku siły potrzebne jest natomiast zachowanie strumienia pędu. Jedna „uniwersalna średnia” nie zachowuje jednocześnie masy, energii i pędu, jeśli profil jest silnie niejednorodny.

Równie rozmieszczone promieniowo punkty sondy grzebieniowej nie reprezentują równych pól pierścieni. Prosta średnia arytmetyczna może nadmiernie ważyć okolice piasty albo wierzchołka. Sonda ciśnienia całkowitego ma ograniczony kąt akceptacji, termopara ma współczynnik odzysku, bezwładność cieplną, przewodzenie po osłonie i wymianę radiacyjną. Nazwa kanału w systemie akwizycji nie dowodzi jeszcze, że zmierzono właściwy mezurand.

Podobieństwo i wielkości skorygowane

Mapy maszyn przepływowych wykorzystują wielkości skorygowane, aby ograniczyć pierwszorzędny wpływ stanu na wlocie. Przy jawnie przyjętych warunkach odniesienia $T_\mathrm{ref}$ i $p_\mathrm{ref}$ można zdefiniować

$$\theta=\frac{T_t}{T_\mathrm{ref}}, \qquad \delta=\frac{p_t}{p_\mathrm{ref}},$$

$$\dot m_c = \dot m\frac{\sqrt{\theta}}{\delta}, \qquad N_c = \frac{N}{\sqrt{\theta}}.$$

Niektóre źródła stosują inną normalizację lub jednostki. Dlatego zapis „corrected flow” bez definicji $T_\mathrm{ref}$, $p_\mathrm{ref}$, stacji wlotowej i użytych własności gazu jest niekompletny. Dokumentacja NPSS przedstawia skorygowany przepływ, skorygowaną prędkość, spręż oraz sprawność jako podstawowe współrzędne map, a jednocześnie ostrzega, że mapa konkretnej maszyny ma ograniczoną przydatność po przeniesieniu do innego układu.¹

Korekcja nie usuwa wpływu:

liczby Reynoldsa;
chropowatości i względnej grubości warstwy przyściennej;
luzu wierzchołkowego;
wymiany ciepła;
wilgotności i składu gazu;
odkształceń cieplnych;
niejednorodności wlotowej;
historii przejścia dynamicznego.

Skorygowana prędkość służy lokalizacji na mapie aerodynamicznej. Rzeczywista prędkość obrotowa nadal określa naprężenia, energię kinetyczną wirnika i częstotliwości wymuszeń. Sterownik może zatem używać współrzędnej skorygowanej do interpretacji mapy, a równocześnie chronić limit rzeczywistych obrotów.

Szczegółowy konflikt podobieństwa Macha, Reynoldsa, luzu względnego, wymiany ciepła i czasu procesu należy analizować osobno. Dla dopasowania drogi gazowej najważniejszy jest wniosek: mapy dwóch geometrycznie podobnych maszyn nie muszą być podobne w stopniu wystarczającym do decyzji kwalifikacyjnej.

Równania zamykające drogę gazową

Bilans masy

Dla objętości kontrolnej w stanie ustalonym suma strumieni wchodzących musi być równa sumie strumieni wychodzących:

$$\sum \dot m_\mathrm{in} - \sum \dot m_\mathrm{out} =0.$$

Dla uproszczonego rdzenia od stacji 2 do 9:

$$\dot m_9 = \dot m_2+\dot m_f-\dot m_b-\dot m_l,$$

gdzie $\dot m_f$ oznacza paliwo, $\dot m_b$ upust opuszczający granicę, a $\dot m_l$ pozostałe strumienie tracone poza nią. Powietrze chłodzące pobrane przed komorą i zwrócone przed turbiną powinno wystąpić w księdze przepływów w obu miejscach, lecz nie może zostać podwójnie odjęte od bilansu całego rdzenia.

Dla objętości o istotnej pojemności gazowej potrzebny jest człon akumulacji:

$$\frac{\mathrm d m_\mathrm{CV}}{\mathrm dt} = \sum \dot m_\mathrm{in} - \sum \dot m_\mathrm{out}.$$

Pominięcie tego członu w modelu przejściowym może sztucznie wymusić natychmiastową zgodność przepływów i ukryć mechanizm pompażu.

Praca sprężarki

Rzeczywista jednostkowa praca sprężarki wynosi

$$w_c=h_{t3}-h_{t2}.$$

Dla kalorycznie doskonałego gazu:

$$w_c=c_{p,a}(T_{t3}-T_{t2}).$$

Sprawność izentropowa całkowita–całkowita jest zdefiniowana jako

$$\eta_c = \frac{h_{t3s}-h_{t2}} {h_{t3}-h_{t2}},$$

a przy stałych $c_p$ i $\gamma$:

$$T_{t3s} = T_{t2}\pi_c^{(\gamma-1)/\gamma}, \qquad \pi_c=\frac{p_{t3}}{p_{t2}}.$$

Stąd

$$T_{t3} = T_{t2} + \frac{T_{t3s}-T_{t2}}{\eta_c}.$$

Moc pobierana przez sprężarkę:

$$P_c=\dot m_c^\mathrm{phys} w_c,$$

gdzie $\dot m_c^\mathrm{phys}$ jest rzeczywistym strumieniem przez komponent, a nie współrzędną skorygowaną mapy.

Komora: bilans energii i strata ciśnienia

Ogólny bilans energii komory w stanie ustalonym można zapisać jako

$$\dot m_3 h_{t3} + \eta_b\dot m_f H_u = \dot m_4 h_{t4} + \dot Q_\mathrm{loss},$$

gdzie $H_u$ jest przyjętą wartością opałową, $\eta_b$ sprawnością spalania, a $\dot Q_\mathrm{loss}$ stratą ciepła przez ściany w granicy modelu. Trzeba jawnie rozstrzygnąć, czy użyto dolnej czy górnej wartości opałowej oraz czy entalpie zawierają energię tworzenia składników.

Komora jest zarazem elementem przepływowym. Jej charakterystyka strat

$$\sigma_b = \frac{p_{t4}}{p_{t3}}$$

nie powinna być automatycznie stałą z punktu nominalnego. Zależy od strumienia, gęstości, rozdziału powietrza, stanu cieplnego i geometrii. Średnia $T_{t4}$ nie opisuje lokalnego gorącego pasma, które może ograniczać turbinę.

Turbina i bilans wału

Rzeczywista jednostkowa praca oddana przez gaz:

$$w_t=h_{t4}-h_{t5}.$$

Sprawność izentropowa turbiny:

$$\eta_t = \frac{h_{t4}-h_{t5}} {h_{t4}-h_{t5s}}.$$

W stanie ustalonym jednowałowego zespołu:

$$P_t = P_c+P_\mathrm{aux}+P_\mathrm{loss},$$

czyli

$$\dot m_t(h_{t4}-h_{t5}) = \dot m_c(h_{t3}-h_{t2}) + P_\mathrm{aux} + P_\mathrm{loss}.$$

NASA w dydaktycznym omówieniu dopasowania sprężarki i turbiny również rozpoczyna od warunku równości pracy na wspólnym wale.² Jest to warunek konieczny, ale nie wystarczający: nadal trzeba uzgodnić zdolności przepływowe turbiny i dyszy.

W stanie przejściowym właściwe równanie ma postać momentową:

$$J\frac{\mathrm d\omega}{\mathrm dt} = \tau_t-\tau_c-\tau_\mathrm{aux}-\tau_\mathrm{loss},$$

lub energetyczną:

$$\frac{\mathrm d}{\mathrm dt} \left(\frac12J\omega^2\right) = P_t-P_c-P_\mathrm{aux}-P_\mathrm{loss}.$$

Chwilowo stała prędkość nie zawsze dowodzi pełnego bilansu: krótki odcinek może być skutkiem filtracji sygnału, małej rozdzielczości lub kompensujących się błędów modelu.

Dysza: bilans pędu i ciąg

Dysza zamienia część entalpii całkowitej na energię kinetyczną. Dla całej objętości kontrolnej zespołu, przy jednym wylocie i z pominięciem pędu doprowadzanego paliwa:

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

$F_\mathrm{inst}$ grupuje tylko oddziaływania instalacyjne nieujęte innymi powierzchniami kontrolnymi. Przy innej granicy zmieni się jawny podział na ciąg strugi, człony ciśnieniowe i opór instalacyjny, lecz wynik fizyczny powinien pozostać zgodny. Oficjalne omówienie NASA wyprowadza tę samą strukturę równania z bilansu pędu i członów ciśnieniowych.³

Pole gardzieli oraz ciśnienie otoczenia zamykają drogę gazową. Jeśli przepływ w krytycznym przekroju jest zdławiony, dalsza zmiana ciśnienia za gardzielą nie zwiększa prostym sposobem strumienia. W przepływie niezdławionym ciśnienie wsteczne silniej przenosi się w górę układu. Dlatego dysza nie jest biernym zakończeniem rachunku: jej zdolność przepływowa współwyznacza ciśnienie za turbiną, stopień rozprężania turbiny i punkt sprężarki.

Mapa sprężarki jako część układu

Mapa sprężarki wiąże zwykle skorygowany strumień $\dot m_c$, spręż $\pi_c$, skorygowaną prędkość $N_c$ i sprawność $\eta_c$. Linie stałej prędkości leżą pomiędzy obszarem małego przepływu a dławieniem. Kontury sprawności, nazywane wyspami, są interpolacją danych, nie odrębnym prawem fizycznym.

Schematycznie:

spręż πc
  ▲
  │         granica niestateczności
  │        /
  │   N3  /───────────────┐
  │      /   • punkt      │
  │ N2  /──── pracy ──────┤  obszar dławienia
  │    /                  │
  │ N1/───────────────────┘
  └────────────────────────────► skorygowany strumień
       mały przepływ       duży przepływ

Kierunki przesunięć trzeba interpretować systemowo:

wzrost oporu za sprężarką przy zbliżonej prędkości zwykle przesuwa punkt ku mniejszemu przepływowi i większemu sprężowi;
zwiększenie zdolności przepływowej dalszej drogi może przesunąć punkt ku większemu przepływowi;
zmiana pola dyszy wpływa przez turbinę, ciśnienie przed nią i równowagę wału;
zaburzenie wlotowe może zmniejszyć lokalny zapas stateczności mimo podobnej średniej wartości odzysku ciśnienia;
zmiana paliwa oddziałuje najpierw na temperaturę i pracę turbiny, następnie na przyspieszenie wału, przepływ sprężarki i położenie punktu na mapie.

Wirujące oderwanie, pompaż i dławienie

Wirujące oderwanie (rotating stall) jest lokalnym, nieosiowosymetrycznym stanem sprężarki. Jedna lub więcej komórek oderwanego przepływu przemieszcza się obwodowo. Średni strumień przez maszynę może pozostać dodatni, choć lokalne obciążenia i pulsacje rosną.

Pompaż (surge) jest niestatecznością całego układu sprężarka–przewody–objętości–opór. Może obejmować duże oscylacje ciśnienia i strumienia, a niekiedy chwilowe odwrócenie przepływu. Jego dynamika zależy od pojemności objętości, przewodów, charakterystyki dławienia i bezwładności wirnika, nie tylko od mapy sprężarki.

Dławienie (choke) oznacza osiągnięcie lokalnego warunku krytycznego ograniczającego dalszy wzrost strumienia. W złożonym komponencie lokalne dławienie może współistnieć z separacją i innymi stratami. Prawa granica mapy bywa również ustalana kryterium spadku sprawności, dlatego każda „linia dławienia” wymaga definicji.

Granica niestateczności z mapy stanowiska komponentowego nie jest automatycznie granicą kompletnego silnika. Inne przewody, objętości, zaburzenie wlotowe, dynamika sterowania i sposób zbliżania się do granicy zmieniają argument. Badania NASA nad zaburzeniami wlotowymi pokazują, że nierównomierność obwodowa wpływa zarówno na osiągi, jak i zakres statecznego przepływu, a reakcja musi być analizowana jako zjawisko nieustalone w układzie wirnika.⁴

Zapas stateczności

Nie istnieje jedna uniwersalna definicja zapasu do pompażu. W używanym raporcie trzeba podać wzór. Dla demonstracji można zdefiniować ciśnieniowy zapas na tej samej linii skorygowanej prędkości:

$$M_\pi = \left( \frac{\pi_{c,\mathrm{lim}}}{\pi_{c,\mathrm{op}}} -1 \right)100\%.$$

Inne definicje uwzględniają jednocześnie zmianę sprężu i strumienia albo rzutują punkt na granicę w inny sposób. Wartości obliczonych różnymi definicjami nie wolno zestawiać bez przeliczenia.

Zapas nominalny trzeba rozdzielić na co najmniej trzy składniki:

$$M_\mathrm{available} = M_\mathrm{nominal} - M_\mathrm{transient} - M_\mathrm{uncertainty} - M_\mathrm{degradation}.$$

Nie jest to równanie probabilistyczne, jeżeli składniki „transient” i „degradation” są rezerwami projektowymi. Nie należy ich bez uzasadnienia łączyć pierwiastkiem z sumy kwadratów z niepewnością pomiarową.

Od wlotu do dyszy: pełny łańcuch przyczynowy

Wlot i zaburzenie

Wlot wyznacza $p_{t2}$, $T_{t2}$, profil prędkości oraz przestrzenny rozkład ciśnienia całkowitego. Utrata odzysku ciśnienia obniża gęstość i potencjalny spręż bezwzględny, lecz sama średnia wartość odzysku nie opisuje stateczności. Dwa wloty o tym samym średnim $p_{t2}$ mogą wytworzyć inny sektor deficytu, wir albo widmo pulsacji.

Sprężarka obracająca się za nierównomiernym wlotem doświadcza okresowego obciążenia. Wolny czujnik może pokazywać spokojną średnią, podczas gdy łopatki widzą szybkozmienne wymuszenie. Kwalifikacja zabudowanego układu wymaga więc zgodnego pasma pomiarowego, synchronizacji oraz deskryptorów odnoszących się do rzeczywistej wrażliwości komponentu.

Sprężarka, dyfuzor i kanał nawrotny

W małych silnikach często stosuje się jednostopniową sprężarkę promieniową, ponieważ w zwartym układzie może uzyskać znaczący spręż. Miniaturyzacja zwiększa jednak względne znaczenie luzu, chropowatości i warstw przyściennych. Monografia Mariana Gierasa pokazuje, że charakterystyka wirnika nie może być oddzielona od dyfuzora i kanału nawrotnego: zmiana tych elementów zmienia straty, kąt napływu i granice całego zespołu.⁵

Dyfuzor odzyskuje część ciśnienia statycznego z energii kinetycznej, lecz nadmierna dyfuzja sprzyja separacji. Kanał nawrotny wytwarza przepływy wtórne oraz nierównomierność na wlocie komory. Mapa „sprężarki” musi zatem określać, czy obejmuje sam wirnik, wirnik z dyfuzorem, czy również dalszy kanał. Zamiana dyfuzora przy zachowaniu mapy starej konfiguracji jest zmianą modelu, nie drobną korektą geometrii.

Komora w łańcuchu przyczynowym

Komora ustala trzy sprzężone wielkości:

spadek ciśnienia całkowitego;
przyrost entalpii;
przestrzenny i czasowy rozkład temperatury przed turbiną.

Wzrost spadku ciśnienia obniża dostępne ciśnienie przed turbiną i dyszą. Zmiana rozdziału powietrza może poprawić stabilizację płomienia, ale zmienić stratę oraz profil temperatury. Lokalny gorący obszar ogranicza dopuszczalną średnią temperaturę, nawet jeśli model zerowymiarowy widzi wyłącznie bezpieczną wartość uśrednioną.

W badaniach cywilnego GTM-120 bezpośredni pomiar ciśnienia całkowitego w małym kanale uznano za nadmiernie zaburzający przepływ i część wyników wyznaczano z pomocą CFD.⁵ Jest to poprawny powód do zastosowania estymatora hybrydowego, ale nie powód do nazwania wyjścia CFD „pomiarem”.

Turbina

Turbina musi jednocześnie przepuścić gaz i dostarczyć pracę sprężarce. Jeżeli jej gardziel jest zdławiona, skorygowany strumień staje się względnie mało wrażliwy na dalsze obniżanie ciśnienia. Niewielka zmiana efektywnego pola gardzieli może jednak przesunąć zdolność przepływową, stopień rozprężania i linię pracy sprężarki.

Łańcuch dla zmniejszenia zdolności przepływowej turbiny wygląda następująco:

mniejsza efektywna gardziel turbiny
        ↓
mniejszy strumień możliwy przy danym pₜ4 i Tₜ4
        ↓
wzrost ciśnienia przed turbiną lub zmiana stopnia rozprężania
        ↓
zmiana ciśnienia za sprężarką i oporu dalszej drogi
        ↓
przesunięcie punktu sprężarki ku mniejszemu przepływowi
        ↓
zmiana zapasu stateczności, pracy sprężarki i temperatur
        ↓
nowa równowaga wału albo przyspieszenie/zwolnienie wirnika

Nie każdy układ zareaguje identyczną wielkością lub nawet identycznym chwilowym kierunkiem, ponieważ sterowanie i dysza mogą przejąć część zaburzenia. Łańcuch wskazuje zmienne, które trzeba rozwiązać, a nie gotowy wynik konkretnego silnika.

Dysza

Zmniejszenie efektywnego pola dyszy zwiększa opór dalszej drogi. Może podnieść ciśnienie za turbiną, zmniejszyć jej stopień rozprężania i pracę, a następnie wpłynąć na prędkość wału oraz punkt sprężarki. W układzie z dyszą regulowaną pole staje się wejściem sterującym, lecz korzyść w jednym wymiarze może pogorszyć ciąg, temperaturę albo obciążenie aktuatora.

Próba w komorze wysokościowej lub kanale odprowadzającym spaliny wymaga korekty wpływu stanowiska. Eżektor, przewód wylotowy, recyrkulacja oraz ssanie instalacji mogą zmienić ciśnienie wsteczne i wskazanie wagi ciągu. Wynik „ciągu silnika” bez opisanej objętości kontrolnej, tarowania i korekty stanowiska jest niepełny.

Iteracyjne wyznaczanie punktu pracy

Dopasowanie drogi gazowej jest rozwiązaniem układu nieliniowych równań. Nie wystarczy przejść przez komponenty jednokrotnie od wlotu do dyszy, ponieważ warunek z ostatniego komponentu wpływa na pierwszy.

Wektor niewiadomych może zawierać na przykład:

$$\mathbf x = \begin{bmatrix} \beta_c & \pi_t & T_{t4} & \dot m & \omega \end{bmatrix}^{T},$$

gdzie $\beta_c$ jest współrzędną pomocniczą na mapie sprężarki. Wektor reszt może mieć postać

$$\mathbf r(\mathbf x) = \begin{bmatrix} r_{\dot m,t}\\ r_{\dot m,n}\\ r_P\\ r_p\\ r_N \end{bmatrix},$$

obejmując niezgodność przepływu przez turbinę, niezgodność dyszy, bilans mocy, warunek ciśnieniowy oraz zadany stan wału. Przykładowa znormalizowana reszta przepływu:

$$r_{\dot m} = \frac{\dot m_\mathrm{upstream}-\dot m_\mathrm{component}} {\dot m_\mathrm{ref}}.$$

Reszta mocy:

$$r_P = \frac{P_t-P_c-P_\mathrm{aux}-P_\mathrm{loss}} {P_\mathrm{ref}}.$$

Solver zmienia tyle niezależnych niewiadomych, ile istnieje niezależnych warunków zgodności. Dokumentacja historycznego programu NNEP89 nazywa takie rozbieżności odpowiednio błędami przepływu i pracy oraz podkreśla konieczność wyboru zmiennych, które mogą je usunąć.⁶

Typowy tok obliczenia jest następujący:

ustalenie otoczenia, konfiguracji i wejść sterujących;
przeliczenie stanu wlotowego do stacji 2;
odczyt sprężu i sprawności z mapy sprężarki;
obliczenie stanu w stacji 3;
zastosowanie charakterystyki straty i bilansu komory;
odczyt zdolności przepływowej i sprawności turbiny;
obliczenie pracy turbiny oraz stanu w stacji 5;
obliczenie zdolności przepływowej dyszy;
utworzenie reszt masy, mocy i ciśnień;
aktualizacja niewiadomych;
powtarzanie do uzyskania zbieżności wszystkich reszt;
sprawdzenie ograniczeń fizycznych, zakresu map i niepewności.

Zbieżność numeryczna nie jest równoznaczna z poprawnością fizyczną. Solver może znaleźć matematyczny punkt w ekstrapolowanym obszarze mapy, z ujemną stratą, nierealną sprawnością albo niezamkniętym strumieniem pomocniczym. Po zakończeniu iteracji potrzebna jest osobna kontrola praw zachowania i domeny ważności.

Studium I — domknięcie nominalnego punktu pracy

Cel i status danych

Studium pokazuje pełną drogę od danych stacji do decyzji. Wszystkie liczby są danymi syntetycznymi autora i nie opisują istniejącego napędu. Przykład nie jest mapą, procedurą uruchomienia ani podstawą doboru geometrii.

Zamierzone użycie modelu: sprawdzenie, czy uproszczony model jednowymiarowy logicznie zamyka nominalny punkt ustalony oraz czy jego księga przepływów i rachunek mocy są spójne.

Warunki odniesienia:

$$T_\mathrm{ref}=288{,}15 \mathrm K, \qquad p_\mathrm{ref}=101{,}325 \mathrm{kPa}.$$

Dane wejściowe:

Wielkość	Wartość	Status
$p_{t2}$	$99{,}2 \mathrm{kPa}$	dane syntetyczne
$T_{t2}$	$291{,}0 \mathrm K$	dane syntetyczne
$\dot m_2$	$0{,}2450 \mathrm{kg/s}$	dane syntetyczne
$N$	$92\,000 \mathrm{obr/min}$	dane syntetyczne
$\dot m_f$	$0{,}0045 \mathrm{kg/s}$	dane syntetyczne
upust opuszczający granicę	$0{,}0030 \mathrm{kg/s}$	dane syntetyczne
przeciek opuszczający granicę	$0{,}0010 \mathrm{kg/s}$	dane syntetyczne
powietrze chłodzące pobierane i zwracane przed turbiną	$0{,}0060 \mathrm{kg/s}$	dane syntetyczne
$T_{t4}$	$1010 \mathrm K$	wynik iterowanego modelu
$P_\mathrm{aux}+P_\mathrm{loss}$	$1{,}20 \mathrm{kW}$	dane syntetyczne

Współrzędne mapy

Skorygowany strumień na wlocie sprężarki:

$$\dot m_c = 0{,}2450 \frac{\sqrt{291{,}0/288{,}15}} {99{,}2/101{,}325} = 0{,}2515 \mathrm{kg/s}.$$

Skorygowana prędkość:

$$N_c = \frac{92\,000} {\sqrt{291{,}0/288{,}15}} = 91\,548 \mathrm{obr/min}.$$

Interpolacja syntetycznej mapy w tej współrzędnej zwraca

$$\pi_c=2{,}25, \qquad \eta_c=0{,}72.$$

Są to wyniki modelu mapy, a nie pomiar.

Stan za sprężarką

Dla $\gamma_a=1{,}4$ i $c_{p,a}=1005 \mathrm{J/(kg\,K)}$:

$$T_{t3s} = 291{,}0\cdot2{,}25^{0{,}285714} = 366{,}87 \mathrm K,$$

$$T_{t3} = 291{,}0+ \frac{366{,}87-291{,}0}{0{,}72} = 396{,}38 \mathrm K.$$

Jednostkowa praca sprężarki:

$$w_c = 1005(396{,}38-291{,}0) = 105{,}91 \mathrm{kJ/kg}.$$

Moc sprężarki:

$$P_c = 0{,}2450\cdot105{,}91 = 25{,}95 \mathrm{kW}.$$

Ciśnienie całkowite za sprężarką:

$$p_{t3}=2{,}25\cdot99{,}2=223{,}2 \mathrm{kPa}.$$

Księga przepływów

Powietrze chłodzące zostaje pobrane z głównego strumienia, ale wraca przed turbiną. Nie zmienia więc masy przecinającej końcową granicę rdzenia, choć zmienia lokalne strumienie przez część komory.

Strumień przed turbiną:

$$\dot m_4 = 0{,}2450+0{,}0045-0{,}0030-0{,}0010 = 0{,}2455 \mathrm{kg/s}.$$

Księga lokalna powinna zarazem pokazać, że przez strefę reaktywną komory przepływa mniej powietrza niż przez turbinę, ponieważ $0{,}0060 \mathrm{kg/s}$ wraca dopiero przed nią.

Pominięcie tego zapisu tworzy charakterystyczny fałszywy wynik: solver może dopasować pracę turbiny przez zmianę sprawności albo „przecieku kalibracyjnego”, mimo że masa w stacjach jest niespójna.

Komora i turbina

Przy syntetycznym spadku ciśnienia całkowitego komory równym 7%:

$$p_{t4}=0{,}93p_{t3}=207{,}58 \mathrm{kPa}.$$

Wymagana moc turbiny:

$$P_t = 25{,}95+1{,}20 = 27{,}15 \mathrm{kW}.$$

Wymagana jednostkowa praca turbiny:

$$w_t = \frac{27{,}15}{0{,}2455} = 110{,}58 \mathrm{kJ/kg}.$$

Dla syntetycznych własności gorącego gazu $c_{p,g}=1150 \mathrm{J/(kg\,K)}$, $\gamma_g=1{,}33$ i $\eta_t=0{,}76$:

$$T_{t5} = 1010- \frac{110\,580}{1150} = 913{,}84 \mathrm K.$$

Z definicji sprawności turbiny wynika przybliżony stosunek ciśnień

$$\frac{p_{t5}}{p_{t4}} = \left[ 1- \frac{w_t} {\eta_t c_{p,g}T_{t4}} \right]^{\frac{\gamma_g}{\gamma_g-1}} = 0{,}583.$$

Stąd

$$p_{t5}=121{,}04 \mathrm{kPa}.$$

Domknięcie dyszy i ciągu

Dla $p_0=101{,}325 \mathrm{kPa}$ syntetyczna dysza pracuje w tym punkcie bez dławienia. Przy $\eta_n=0{,}94$ temperatura po izentropowym rozprężeniu do otoczenia wynosi

$$T_{9s} = T_{t5} \left(\frac{p_0}{p_{t5}}\right)^{(\gamma_g-1)/\gamma_g} = 874{,}41 \mathrm K.$$

Prędkość wylotowa:

$$V_9 = \sqrt{ 2\eta_n c_{p,g}(T_{t5}-T_{9s}) } = 291{,}97 \mathrm{m/s}.$$

Dla próby statycznej, idealnego dopasowania ciśnienia w płaszczyźnie 9 i bez osobnego członu instalacyjnego:

$$F_N = \dot m_9V_9 = 0{,}2455\cdot291{,}97 = 71{,}68 \mathrm N.$$

Wartość $71{,}68 \mathrm N$ jest wynikiem modelu syntetycznego. Nie jest przewidywaniem dla GTM-120 ani innego rzeczywistego silnika.

Historia iteracji

Właściwe kryterium zbieżności obejmuje wszystkie niezależne reszty:

Iteracja	reszta przepływu turbiny	reszta dyszy	reszta mocy
0	$+4{,}8\%$	$-3{,}2\%$	$-7{,}6\%$
1	$+1{,}2\%$	$-0{,}8\%$	$+2{,}0\%$
2	$-0{,}18\%$	$+0{,}09\%$	$+0{,}32\%$
3	$+0{,}03\%$	$+0{,}02\%$	$-0{,}06\%$

Zbieżność jedynie reszty mocy w iteracji 2 nie byłaby wystarczająca. Dopiero iteracja 3 spełnia syntetyczne kryterium $|r_i|<0{,}1\%$ dla wszystkich trzech warunków.

Budżet niepewności

Model pomiaru ciągu ma postać

$$F_N=f(\dot m,p_{t5},T_{t5},p_0,\eta_n,\text{mapy},\text{stanowisko}).$$

Zgodnie z GUM niepewność złożoną dla linearyzowanego modelu oblicza się z uwzględnieniem współczynników wrażliwości i kowariancji:⁷

$$u_c^2(F) = \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).$$

Syntetyczny budżet standardowych wkładów do ciągu:

Źródło	Standardowy wkład do $u(F)$
strumień masy i jego redukcja	$0{,}65 \mathrm N$
ciśnienia oraz model rozprężania	$1{,}05 \mathrm N$
temperatury i własności gazu	$0{,}55 \mathrm N$
interpolacja i postać map komponentów	$1{,}30 \mathrm N$
tarowanie oraz oddziaływanie stanowiska	$0{,}80 \mathrm N$
niepewność konfiguracji i przecieków	$0{,}75 \mathrm N$

Przy syntetycznym założeniu braku korelacji:

$$u_c(F)=2{,}17 \mathrm N, \qquad U(F)=2u_c(F)=4{,}34 \mathrm N.$$

Wynik należy więc zapisać jako

$$F_N=(71{,}7\pm4{,}4) \mathrm N, \qquad k=2.$$

Niepewność mocy sprężarki, po uwzględnieniu strumienia, temperatur i mapy sprawności, wynosi w tym przykładzie $U(P_c)=1{,}3 \mathrm{kW}$ dla $k=2$.

Historyczna analiza NASA dla badań turbojetu wykazała, że w wynikach takich jak strumień wlotowy, ciąg netto i jednostkowe zużycie paliwa wkład błędu systematycznego może być dwu- lub trzykrotnie większy od wkładu precyzji.⁸ Duża liczba próbek nie usuwa więc błędnej kalibracji, modelu stanowiska ani wspólnego obciążenia kanałów.

Zapas i decyzja

Syntetyczna mapa daje w tym punkcie $\pi_{c,\mathrm{lim}}=2{,}47$. Dla przyjętej definicji:

$$M_\pi = \left(\frac{2{,}47}{2{,}25}-1\right)100\% = 9{,}8\%.$$

Rozszerzona niepewność położenia granicy i punktu wynosi syntetycznie $2{,}3$ punktu procentowego. Rezerwa degradacyjna wynosi $2{,}0$ p.p., a rezerwa przejściowa $3{,}0$ p.p. Konserwatywny pozostały zapas:

$$M_\mathrm{res} = 9{,}8-2{,}3-2{,}0-3{,}0 = 2{,}5 \mathrm{p.p.}$$

Decyzja: model zamyka nominalny punkt i nadaje się do demonstracji procedury. Nie kwalifikuje napędu, ponieważ mapy, rezerwy i budżet są syntetyczne, a przypadek nie został zwalidowany w stanach przejściowych ani po zabudowie.

Regresja punktu nominalnego

Ponowne wykonanie niezmienionej konfiguracji SYN-NOM-01 powinno odtworzyć:

$\dot m_c=0{,}25148 \mathrm{kg/s}$ z tolerancją $5\cdot10^{-5} \mathrm{kg/s}$;
$P_c=25{,}947 \mathrm{kW}$ z tolerancją $0{,}010 \mathrm{kW}$;
$T_{t5}=913{,}84 \mathrm K$ z tolerancją $0{,}05 \mathrm K$;
$F_N=71{,}68 \mathrm N$ z tolerancją $0{,}05 \mathrm N$;
każdą znormalizowaną resztę poniżej $0{,}1\%$.

To test implementacji i konfiguracji, nie walidacja fizyki.

Punkt projektowy, praca poza projektem i dynamika wału

Punkt projektowy służy do doboru rozmiarów, spręży, temperatur i pól charakterystycznych. Nie musi być najczęściej używanym stanem ani punktem największego ciągu. Optymalizacja wyłącznie jednego punktu może pogorszyć rozruch, stan małej mocy, działanie na wysokości, zdolność przyspieszania i tolerancję zaburzeń wlotowych.

Model pracy poza punktem projektowym wykorzystuje mapy do ponownego rozwiązania dopasowania przy zmianie:

ciśnienia i temperatury otoczenia;
prędkości lotu;
wejść sterujących;
pola dyszy lub innej geometrii zmiennej;
upustów i poborów mocy;
strat, luzów i chropowatości;
stanu cieplnego;
konfiguracji instalacyjnej.

Region bez danych nie powinien być ukrywany przez gładką interpolację. Na mapie i w wynikach należy rozróżnić:

interpolację między danymi;
kontrolowaną ekstrapolację;
obszar opisany modelem zastępczym;
obszar niedopuszczony do użycia.

W stanie przejściowym punkt sprężarki nie porusza się po ustalonej linii pracy. Zwiększenie dopływu energii podnosi temperaturę przed turbiną i jej moment, ale wał, przepływ przez objętości, ciśnienia i czujniki reagują z różnymi stałymi czasu. Trajektoria może chwilowo zbliżyć się do granicy niestateczności bardziej niż dowolny punkt ustalony.

Narzędzia dynamiczne NASA służą między innymi do oceny kompromisu między czasem odpowiedzi a zapasem sprężarki; ich dokumentacja wprost odróżnia analizę ustaloną od zamkniętej pętli przejściowej.⁹ Nie oznacza to, że konkretne wyniki dla dużego turbofanu można przenieść ilościowo do miniaturowego turbojetu. Przenoszalny jest sposób argumentacji.

Studium II — przyspieszenie przy zaburzeniu wlotowym

Pytanie i hipotezy

Studium wykorzystuje wyłącznie dane syntetyczne autora. Nie przedstawia harmonogramu paliwa ani progów rzeczywistego sterowania.

Pytanie: czy zapas $9{,}8\%$ wyznaczony w nominalnym stanie ustalonym wystarcza do decyzji o dopuszczeniu przejściowego przyspieszenia przy zaburzonym wlocie?

Hipotezy:

$H_S$: linia stanów ustalonych wystarcza, ponieważ trajektoria przejściowa pozostaje blisko niej;
$H_T$: bezwładność wału i pojemności gazowe tworzą odrębną trajektorię, zużywając część zapasu;
$H_D$: zaburzenie wlotowe oraz degradacja nakładają się na efekt przejściowy i mogą zmienić decyzję;
$H_0$: niepewność granicy jest tak duża, że wynik pozostaje nierozstrzygający.

Model dynamiczny

Minimalny model zawiera:

$$J\dot\omega = \tau_t-\tau_c-\tau_\mathrm{loss},$$

$$\frac{\mathrm d m_v}{\mathrm dt} = \dot m_\mathrm{in}-\dot m_\mathrm{out},$$

$$\frac{\mathrm d U_v}{\mathrm dt} = \sum\dot m h_{t,\mathrm{in}} - \sum\dot m h_{t,\mathrm{out}} + \dot Q - \dot W,$$

a także modele zwłoki spalania, sensorów, aktuatorów i logiki ograniczającej. Mapa ustalona nadal dostarcza charakterystyki komponentów, ale nie zastępuje równań akumulacji i wału.

Syntetyczne zaburzenie jest opisane przez czasową zmianę średniego odzysku ciśnienia oraz niezależny deskryptor nierównomierności obwodowej. Użycie tylko średniego odzysku odpowiadałoby innej hipotezie fizycznej.

Wyniki

Stan	Zapas według jawnej definicji	Status
początkowy punkt ustalony	$9{,}8\%$	wynik modelu syntetycznego
minimum przejścia bez zaburzenia	$6{,}7\%$	wynik modelu syntetycznego
minimum z zaburzeniem wlotowym	$4{,}1\%$	wynik modelu syntetycznego
minimum po dodaniu scenariusza degradacji	$2{,}8\%$	wynik modelu syntetycznego

Łańcuch przyczynowy jest następujący:

wzrost doprowadzanej energii
        ↓
wzrost Tₜ4 i momentu turbiny
        ↓
przepływ sprężarki reaguje wolniej niż obciążenie dalszej drogi
        ↓
punkt odchodzi od ustalonej linii pracy
        ↓
sektor niższego pₜ2 zwiększa lokalne obciążenie części sprężarki
        ↓
efektywna granica stateczności zbliża się do trajektorii
        ↓
degradacja mapy zmniejsza pozostały zapas

Wynik wspiera $H_T$ i $H_D$: utrata zapasu nie jest możliwa do odtworzenia samą serią stanów ustalonych.

Niepewność i decyzja

Rozszerzona syntetyczna niepewność minimalnego zapasu wynosi $1{,}9$ p.p. Obejmuje położenie granicy, błąd mapy, czasową synchronizację, deskryptor zaburzenia i niepewność modelu dynamicznego. Konserwatywna dolna wartość dla przypadku zdegradowanego:

$$M_{\mathrm{lower}} = 2{,}8-1{,}9 = 0{,}9 \mathrm{p.p.}$$

Przy syntetycznym kryterium decyzji równym $2{,}0$ p.p. przypadek nie przechodzi bramki.

Decyzja: nie dopuszczać ocenianego przejścia na podstawie modelu ustalonego. Potrzebne są zwalidowany model dynamiczny, reprezentatywne zaburzenie, kontrola wspólnej podstawy czasu i niezależny dowód położenia granicy. Wynik nie wskazuje, jak zmienić rzeczywisty harmonogram sterowania.

Regresja modelu przejściowego

Dla konfiguracji SYN-TR-02 zamrożone są dane wejściowe, krok czasu, interpolacja map i przebieg zaburzenia. Test regresyjny powinien odtworzyć:

minimum bez zaburzenia $6{,}7\%\pm0{,}05$ p.p.;
minimum z zaburzeniem $4{,}1\%\pm0{,}05$ p.p.;
minimum zdegradowane $2{,}8\%\pm0{,}05$ p.p.;
czas minimum w granicy jednego kroku obliczeniowego;
brak użycia komórek mapy oznaczonych jako niedopuszczone.

Zmiana modelu zaburzenia, dynamiki wału albo granicy stateczności jest zmianą tezy fizycznej i wymaga ponownej pełnej walidacji, a nie tylko regresji liczbowej.

Sterowanie jako część dopasowania drogi gazowej

Sterowanie nie działa obok termodynamiki. Wybiera trajektorię w przestrzeni map, wykorzystując mierzone lub estymowane wielkości. Typowe wejścia funkcjonalne obejmują żądany stan, prędkość wału, ciśnienia, temperatury i stan ograniczeń. Wyjściem jest polecenie dla układu wykonawczego dozującego energię, a w odpowiednich architekturach także dla dyszy, upustu albo geometrii zmiennej.

Każdy tor zawiera:

błąd i dryft sensora;
filtrację oraz opóźnienie;
ograniczenie szybkości aktuatora;
nieliniowość i histerezę;
niepewność estymatora;
arbitraż między odpowiedzią a ochroną limitów.

Jeżeli model zakłada natychmiastowy pomiar i wykonanie polecenia, może przewidywać bezpieczną trajektorię niemożliwą do zrealizowania. Jeśli limit jest sprawdzany dopiero po obliczeniu kroku, model może przeoczyć krótkie przekroczenie. Kwalifikacja drogi gazowej obejmuje więc również interfejs z regulatorem, choć nie musi ujawniać konkretnej logiki sterowania.

Hierarchia modeli i podział odpowiedzialności

Różne poziomy modelu odpowiadają na odmienne pytania:

Poziom	Założenia i wejścia	Główne wyniki	Typowy koszt	Typowy błąd	Dopuszczalne użycie
obieg zerowymiarowy	stany średnie, sprawności lub mapy	bilanse, ciąg, zużycie, temperatury	niski	ukrycie pól lokalnych	przesiew koncepcji i dopasowanie
model średnioliniowy/1D	geometria średnia, trójkąty prędkości, korelacje strat	wstępne mapy, rozkład pracy	niski–średni	zależność od korelacji	trendy i budowa hipotez mapy
model systemowy mapowy	mapy komponentów, objętości, wały, sterowanie	obwiednia ustalona i przejściowa	średni	błąd map i ekstrapolacji	analiza współpracy komponentów
RANS pojedynczego kanału	periodyczność, model turbulencji, uśrednienie wirnik–stator	lokalne pola i trendy strat	średni–wysoki	wykluczenie modów niesymetrycznych	mechanizmy lokalne z dala od niestateczności
CFD pełnoobwodowe nieustalone	pełna geometria, krok czasu, ruch względny	mody obwodowe i wymuszenia	bardzo wysoki	model turbulencji, koszt i warunki brzegowe	wybrane krytyczne mechanizmy
stanowisko komponentowe	rzeczywisty komponent, sztuczne granice	mapa i lokalna diagnostyka	wysoki	wpływ stanowiska i brak integracji	identyfikacja charakterystyki komponentu
silnik zintegrowany	pełna droga gazowa	dopasowanie i osiągi zespołu	bardzo wysoki	ograniczona obserwowalność	walidacja systemowa
próba po zabudowie	reprezentatywna instalacja	wpływ wlotu, wylotu i cieplny	najwyższy	ograniczony zakres i ryzyko	końcowy dowód integracji

Model wyższego rzędu nie jest automatycznie bardziej wiarygodny. Pełnoobwodowe CFD ze źle określonym wlotem może być mniej użyteczne niż prosty model z dobrze zmierzonymi warunkami. Stanowisko komponentowe może dokładnie wyznaczyć mapę w swojej konfiguracji, lecz nie dowieść zachowania po zabudowie.

NASA zademonstrowała sprzężenie zerowymiarowego modelu cyklu z trójwymiarowymi modelami komponentów, w którym warunki brzegowe przepływały z modelu systemowego do CFD, a częściowe mapy wracały do modelu cyklu.¹⁰ Najważniejszą lekcją nie jest konkretna implementacja, lecz konieczność iteracji: wynik CFD zmienia rozwiązanie dopasowania, a ono tworzy nowe warunki brzegowe dla CFD.

Kiedy modele mogą się korygować

Model średnioliniowy może generować kształt mapy przed próbą. CFD może wskazać mechanizm odpowiedzialny za zmianę nachylenia lub straty. Próba komponentowa może skorygować zdolność przepływową i sprawność. Model systemowy sprawdza, jak korekta przemieszcza linię pracy całego silnika.

Nie wolno jednak liczyć zgodności dwóch modeli jako dwóch niezależnych dowodów, jeśli oba używają:

tej samej geometrii nominalnej;
tej samej korelacji strat;
tego samego warunku brzegowego;
tej samej mapy bazowej;
tego samego pomiaru do kalibracji.

Zbieżność modelu obiegowego i CFD może wynikać ze wspólnego założenia, a nie z poprawności fizyki. Monografia Piotra Tarnawskiego trafnie pokazuje różnicę między koncepcyjnym potencjałem wynikającym z CFD a potwierdzeniem eksperymentalnym: zgodność dwóch obliczeń nie zastępuje próby rzeczywistego mechanizmu.¹¹

Budowa, skalowanie i ekstrapolacja map

Mapa ze stanowiska

Mapa stanowiskowa wymaga znanych warunków wlotowych, prędkości, przeciwciśnienia, geometrii oraz instrumentacji. Powinna przechowywać surowe punkty i niepewność, a nie tylko wygładzone kontury. Definicja granicy niestateczności musi wskazywać metodę jej rozpoznania i bezpieczny zakres dowodu.

Mapa jest ważna dla konkretnego:

wirnika i statora;
dyfuzora i kanałów;
luzu;
chropowatości;
stanu cieplnego;
gazu;
stanowiska;
sposobu redukcji danych.

Mapa średnioliniowa

Model średnioliniowy wykorzystuje trójkąty prędkości, korelacje poślizgu, strat, blokady i przecieków. Jest szybki i dobry do analizy wrażliwości, ale jego korelacje reprezentują określone rodziny maszyn. Kalibracja w jednym punkcie nie ustala poprawnie nachylenia linii prędkości, granicy dławienia ani granicy niestateczności.

Mapa z CFD

Serię punktów można uzyskać przez zmianę przeciwciśnienia lub strumienia. Każdy punkt wymaga kontroli:

bilansu masy, energii i momentu;
zbieżności iteracyjnej;
niezależności od siatki i kroku czasu;
rozmiaru domeny;
modelu przejścia i turbulencji;
szczelin, chropowatości i wymiany ciepła;
sposobu uśredniania;
kryterium uznania rozwiązania za fizycznie ustalone.

Stacjonarne RANS pojedynczego kanału okresowego nie może samodzielnie wyznaczyć systemowego pompażu. Taka domena z definicji wyklucza część modów obwodowych i dynamikę objętości całego układu.

Skalowanie

Proste skalowanie mapy może przesunąć przepływ, spręż i sprawność tak, aby trafić w punkt projektowy:

$$\dot m_c' = S_W\dot m_c,$$

$$\pi_c'-1=S_\pi(\pi_c-1),$$

$$\eta_c'=S_\eta\eta_c.$$

Taka operacja nie gwarantuje poprawnego:

nachylenia linii prędkości;
odstępu między liniami;
przebiegu dławienia;
granicy niestateczności;
zmiany liczby Reynoldsa;
zachowania poza punktem kalibracji.

Skalowanie jest hipotezą. Powinno otrzymać zakres ważności, analizę wrażliwości i osobny składnik niepewności postaci modelu.

Ekstrapolacja

Model rozruchu, zatrzymania lub awarii może potrzebować charakterystyk poza obszarem pomiarowym. Mechaniczne przedłużenie splajnu jest niebezpieczne: może wytworzyć sprawność większą od jedności, ujemną stratę, niefizyczny zwrot nachylenia albo sztuczną stateczność.

Każda komórka ekstrapolowana powinna mieć:

flagę pochodzenia;
model zastępczy;
ograniczenia termodynamiczne;
rosnącą niepewność;
warunek niedopuszczenia do decyzji, jeśli mechanizm zmienia reżim.

Przycięcie wartości do ostatniego punktu mapy nie rozwiązuje problemu. Ukrywa brak danych i może sztucznie zatrzymać trajektorię na granicy.

Metrologia i rachunek niepewności

Od mezurandu do wyniku

Łańcuch dowodowy powinien mieć cztery warstwy:

mezurand
   ↓
czujnik i instalacja
   ↓
akwizycja, czas i kalibracja
   ↓
model redukcji oraz wynik pochodny

Dla każdego wyniku trzeba oddzielić:

precyzję — rozrzut powtórzeń;
obciążenie systematyczne — kalibrację, dryft, model stanowiska;
niepewność instalacyjną — profil, kąt, blokadę i przewody;
niepewność modelu redukcji — własności gazu, korekcje, interpolację i CFD;
niepewność konfiguracji — niepełną wiedzę o geometrii, luzie i przeciekach.

Mateusz Turkowski przedstawia ogólny model pomiaru $Y=f(X_1,\ldots,X_N)$, podział na oceny typu A i B oraz propagację przez współczynniki wrażliwości.¹² Dla granicy $\pm a$ o rozkładzie prostokątnym:

$$u=\frac{a}{\sqrt3},$$

a dla rozkładu trójkątnego:

$$u=\frac{a}{\sqrt6}.$$

Niepewność rozszerzona ma postać

$$U=ku_c.$$

Współczynnik $k=2$ bywa używany dla pokrycia zbliżonego do 95%, lecz nie jest automatyczną gwarancją takiego prawdopodobieństwa przy małej liczbie stopni swobody, silnej nieliniowości lub rozkładzie nienormalnym.

Korelacje

Jeżeli dwa ciśnienia mierzy ten sam skaner, dzielą odniesienie albo korzystają z tej samej kalibracji, nie są niezależne. Dla ilorazu

$$\pi=\frac{p_{t3}}{p_{t2}}$$

względna wariancja w przybliżeniu wynosi

$$\left(\frac{u(\pi)}{\pi}\right)^2 \approx \left(\frac{u(p_{t3})}{p_{t3}}\right)^2 + \left(\frac{u(p_{t2})}{p_{t2}}\right)^2 - 2\rho \frac{u(p_{t3})}{p_{t3}} \frac{u(p_{t2})}{p_{t2}},$$

gdzie $\rho$ jest współczynnikiem korelacji. Pominięcie korelacji może zawyżyć albo zaniżyć wynik.

Dla silnie nieliniowego modelu dyszy, granicy dławienia albo interpolacji blisko krawędzi mapy linearyzacja powinna zostać sprawdzona propagacją numeryczną lub metodą Monte Carlo. BIPM publikuje zarówno podstawowy GUM, jak i suplement poświęcony propagacji rozkładów.⁷

Instrumentacja drogi gazowej

Sonda ciśnienia całkowitego. Jej odpowiedź zależy od kąta napływu. Trzon wytwarza ślad i blokadę. W małym kanale kilka sond może zmienić ten sam profil, który mają mierzyć.

Termopara. Odczyt zależy od odzysku temperatury, promieniowania, przewodzenia i stałej czasu. Osobne wygładzanie temperatury i strumienia może zniekształcić chwilowy strumień entalpii.

Pomiar prędkości wału. Licznik impulsów potrzebuje poprawnej liczby impulsów na obrót, detekcji brakujących impulsów i znacznika czasu. W przejściu opóźnienie filtru zmienia obliczone $\dot\omega$.

Pomiar przepływu. Zwężka, dysza, venturi lub wzorcowany wlot wymagają pomiaru ciśnienia, temperatury i współczynnika przepływu. Pulsacja oraz zaburzenie profilu mogą zmienić wynik. Przepływ odczytany z mapy nie jest niezależnym pomiarem do walidacji tej samej mapy.

Pomiar ciągu. Waga widzi silnik, mocowanie, przewody, siły poprzeczne i dryft cieplny. Tarowanie zimne może nie odpowiadać gorącemu. Korekta instalacyjna powinna mieć własny model i niepewność.

Synchronizacja. Moc $P=\tau\omega$, bilans entalpii i trajektoria na mapie wymagają wspólnej osi czasu. Kanały o różnych opóźnieniach mogą utworzyć pozorne źródło albo stratę energii.

Tolerancje, degradacja i populacja

Nominalna mapa opisuje nominalną reprezentację, nie flotę. Rozrzut egzemplarzy wynika z tolerancji kątów, gardzieli, luzów, chropowatości, współosiowości, uszczelnień, montażu i czujników. Te zmienne bywają skorelowane przez proces produkcyjny.

Złożenie wszystkich tolerancji w jednym geometrycznym „najgorszym przypadku” może utworzyć egzemplarz niemożliwy do wyprodukowania. Z kolei niezależne losowanie wszystkich wymiarów może zaniżyć ryzyko, jeśli jeden proces przesuwa kilka cech równocześnie. Potrzebne są dane produkcyjne, macierz korelacji i model hierarchiczny rozdzielający:

zmienność egzemplarzy;
zmienność pomiaru;
zmienność między kampaniami;
dryft w cyklu życia.

Luz

Luz wierzchołkowy zależy od montażu na zimno, rozszerzalności cieplnej, wzrostu odśrodkowego, ruchu łożysk i owalizacji obudowy. Mniejszy luz może poprawić sprawność, ale zwiększyć ryzyko przytarcia. Pomiar szczelinomierzem w stanie zimnym nie opisuje gorącego stanu przejściowego.

Chropowatość, erozja i osady

Chropowatość zmienia tarcie, przejście warstwy przyściennej i separację. Erozja zaokrągla krawędzie oraz zwiększa luz. Osady mogą zmniejszyć efektywne pole gardzieli i zmienić kąty. Degradacji nie należy sprowadzać do jednego mnożnika sprawności. Model stanu powinien, na ile dane pozwalają, oddzielać:

zdolność przepływową;
spręż;
sprawność;
przeciek;
stratę komory;
pole dyszy.

Parametry te mogą być słabo identyfikowalne. Ten sam spadek ciągu może wynikać z kilku kombinacji zmian. Bez dodatkowej instrumentacji odwrotne „wyznaczenie zużycia” nie musi mieć jednoznacznego rozwiązania.

Uszkodzenie obcym ciałem

Lokalne uszkodzenie może zmienić przepływ, wyważenie i odpowiedź drganiową. Krótka próba funkcjonalna po zdarzeniu nie dowodzi odzyskania zapasu stateczności ani trwałości. Ocena, naprawa i przywrócenie do użytkowania muszą wynikać z zatwierdzonej dokumentacji; nie wolno ich improwizować na podstawie ogólnego modelu.

Studium III — rozbieżność stanowisko–CFD–silnik zintegrowany

Dane i konfiguracje

Studium jest zanonimizowaną demonstracją diagnostyczną. Wszystkie liczby są syntetyczne, ale kolejność kontroli odpowiada prawidłowej praktyce.

Porównywane konfiguracje:

MAP-R3 — mapa stanowiskowa komponentu;
CFD-C17 — model stacjonarny komponentu;
TEST-E04 — próba zintegrowanego silnika.

W pozornie jednakowym punkcie otrzymano:

Źródło	Spręż $\pi_c$	Status
stanowisko komponentowe	$2{,}25\pm0{,}04$	pomiar syntetyczny, $k=2$
CFD	$2{,}31$	wynik modelu; $U_\mathrm{num}=0{,}02$, niepełne $U_\mathrm{form}$
silnik zintegrowany	$2{,}12\pm0{,}05$	wartość syntetyczna wyznaczona pośrednio

Pierwsza pokusa polega na skalibrowaniu mapy do $2{,}12$. Byłoby to przedwczesne, ponieważ trzy liczby nie odnoszą się jeszcze do tego samego mezurandu.

Hipotezy

$H_1$: stacje mają różne fizyczne położenie;
$H_2$: profile uśredniono odmiennie;
$H_3$: sondy na stanowisku zmieniły pole przepływu;
$H_4$: użyto innej wersji dyfuzora lub mapy;
$H_5$: stanowisko zintegrowane ma inny odzysk wlotu albo przeciwciśnienie;
$H_6$: CFD pomija przeciek, chropowatość lub wymianę ciepła;
$H_7$: rzeczywisty komponent odbiega od modelu mimo poprawnej redukcji.

Krok 1 — definicja stacji

MAP-R3 definiowała wylot bezpośrednio za dyfuzorem, natomiast TEST-E04 raportowała ciśnienie po kanale nawrotnym. Strata kanału została więc przypisana sprężarce tylko w jednym zbiorze.

Po przeniesieniu danych testu do wspólnej granicy spręż rośnie z $2{,}12$ do $2{,}15$. Jest to zmiana definicji mezurandu, nie korekta wyniku na lepszy.

Krok 2 — uśrednianie

Stanowisko komponentowe stosowało średnią ważoną strumieniem masy. Pierwsza redukcja próby silnika używała średniej arytmetycznej trzech punktów w przekroju z silnym profilem promieniowym.

Ponowne przeliczenie zachowujące strumień podnosi wartość z $2{,}15$ do $2{,}17$. Surowe profile pozostają niezbędne; same średnie nie pozwoliłyby wykonać korekty.

Krok 3 — blokada i współrzędna mapy

Grzebień stanowiska komponentowego zajmował niepomijalną część małego kanału. Jego korekta była zastosowana do ciśnienia, ale nie do efektywnego pola użytego przy redukcji strumienia. Po poprawce skorygowany strumień punktu zmienia się z $0{,}2515$ do $0{,}2480 \mathrm{kg/s}$.

Na pochyłej linii prędkości odpowiada to zmianie wartości odniesienia mapy z $2{,}25$ do $2{,}22$. Nie należy korygować samego sprężu bez przesunięcia współrzędnej.

Krok 4 — wersja geometrii

MAP-R3 dotyczyła wcześniejszej wersji dyfuzora, natomiast TEST-E04 i CFD-C17 — wersji późniejszej. Po użyciu zgodnej mapy oczekiwanie dla punktu testu wynosi $2{,}20$, nie $2{,}22$.

Zestawienie po kontroli konfiguracji:

Źródło	Wartość na wspólnej stacji i współrzędnej
mapa zgodnej konfiguracji	$2{,}20\pm0{,}04$
silnik zintegrowany po ponownej redukcji	$2{,}17\pm0{,}05$
CFD	$2{,}31$

Różnica mapa–silnik wynosi około $1{,}4\%$ i mieści się w złożonej niepewności tych dwóch wyników. Nie dowodzi idealnej zgodności, ale usuwa podstawę do arbitralnego skalowania mapy.

CFD pozostaje wyżej o około 5–6%. Przegląd modelu ujawnia gładkie ściany, brak jednego przecieku i nominalny luz. Kierunek błędu jest zgodny z hipotezą $H_6$, lecz sama zgodność kierunku nie wyznacza wartości korekty. Potrzebne są obliczenia wrażliwości oraz niezależne pomiary lub metrologia geometrii.

Decyzja

Nie zmieniać mapy bazowej wyłącznie po to, aby wymusić zgodność z pierwszą redukcją testu.
Przyjąć ponownie zredukowany punkt $2{,}17\pm0{,}05$ jako wynik zintegrowany dla oznaczonej konfiguracji.
Ograniczyć CFD-C17 do interpretacji mechanizmu i trendu; nie używać jego bezwzględnego sprężu do kwalifikacji.
Utworzyć nową wersję modelu dopiero po wprowadzeniu fizycznie uzasadnionych zmian.
Zachować pierwotny zbiór oraz pełny ślad redukcji.

Regresja i kontrola zmiany

Pakiet regresyjny powinien obejmować:

mapowanie nazw stacji na fizyczne płaszczyzny;
reprodukcję średniej masowej $2{,}17$;
reprodukcję skorygowanego strumienia $0{,}2480 \mathrm{kg/s}$;
zgodność identyfikatora dyfuzora;
bilans masy stanowiska;
test syntetycznego profilu, dla którego średnia arytmetyczna i masowa celowo różnią się;
kontrolę, że dane kalibracyjne nie trafiają do zbioru niezależnej walidacji.

Zmiana modelu turbulencji, szczeliny, chropowatości, wersji mapy lub definicji stacji tworzy nową konfigurację. Jeżeli zmienia zasadniczą tezę o mechanizmie rozbieżności, wymaga nowej walidacji, a nie tylko porównania różnic plików.

Od wymagania do decyzji kwalifikacyjnej

Kwalifikacja modelu nie jest pochwałą jego realizmu. Jest argumentem, że oznaczony wynik nadaje się do określonej decyzji. Aktywny NASA-STD-7009B wiąże praktyki modelowania z kryteriami akceptacji ustalanymi przez program oraz z ograniczaniem ryzyka decyzji zależnych od modelu.¹³ Standard jest źródłem ogólnej struktury wiarygodności, nie normą kwalifikacji konkretnego silnika wojskowego.

Proces powinien przebiegać następująco:

wymaganie
   ↓
wielkość oceniana i kryterium
   ↓
domena oraz konfiguracja
   ↓
model i plan eksperymentu
   ↓
weryfikacja numeryczna i metrologiczna
   ↓
walidacja na danych niezależnych
   ↓
propagacja niepewności i marginesów
   ↓
decyzja, ograniczenia i regresja

Zamierzone użycie

Przykładowe zamierzone użycia różnią się wymaganym dowodem:

Użycie	Główna wielkość	Minimalny dowód
porównanie dwóch koncepcji	różnica osiągów	poprawny trend i niepewność różnicy
przewidywanie ciągu	ciąg netto	model stanowiska, tarowanie, walidacja integralna
ocena temperatury turbiny	lokalne pole temperatury	reprezentatywna instrumentacja i model transferu ciepła
projektowanie sterowania	odpowiedź przejściowa i limity	dynamika wału, objętości, sensorów i aktuatorów
planowanie próby	oczekiwane obciążenia i zakres	konserwatywne obwiednie oraz jawne ekstrapolacje
kwalifikacja drogi gazowej	marginesy w całej obwiedni	dowody komponentowe, zintegrowane, środowiskowe i instalacyjne

Model zwalidowany dla ustalonego ciągu na poziomie morza nie jest przez to zwalidowany dla trajektorii pompażu, rozruchu, lokalnej temperatury łopatki ani lotu na wysokości.

Macierz dowodów

Końcowy argument powinien być zapisany w macierzy:

Twierdzenie	Wielkość oceniana	Źródło modelowe	Źródło pomiarowe	Niepewność	Domena	Decyzja
masa jest rozliczona	reszta bilansu	księga przepływów	niezależne pomiary wlotu, paliwa i odpływów	pomiar + redukcja	stany ustalone i przejściowe	zaliczone/niezaliczone
wał jest zbilansowany	reszta mocy lub momentu	model komponentów i strat	prędkość, temperatury, ewentualny moment	model + synchronizacja	pełny zakres prędkości	decyzja
zachowano stateczność	minimalny jawnie zdefiniowany zapas	mapa i model dynamiczny	szybkie ciśnienia, reprezentatywne zaburzenie	granica + trajektoria	nominalnie i po degradacji	decyzja
ciąg spełnia wymaganie	ciąg netto	model objętości kontrolnej	waga z tarowaniem	bias + precision + stanowisko	otoczenie i zabudowa	decyzja
temperatura jest dopuszczalna	średnia i lokalne maksimum	komora–turbina–termika	grzebień, pirometria lub inne metody	przestrzeń + czas + odzysk	pełna obwiednia	decyzja
wynik dotyczy populacji	rozkład wielkości	analiza tolerancji	seria egzemplarzy i inspekcje	hierarchiczna	produkcja i cykl życia	decyzja

„Zaliczone” musi oznaczać spełnienie wcześniej ustalonego kryterium. Samo podobieństwo wykresów nie jest kryterium walidacji.

Weryfikacja

Weryfikacja odpowiada na pytanie, czy równania i ich implementację rozwiązano poprawnie. Obejmuje:

testy jednostkowe właściwości i komponentów;
przypadki analityczne;
bilanse zachowania;
zbieżność iteracyjną;
badanie siatki i kroku czasu;
testy interpolacji map;
testy granic oraz jednostek;
regresję wyników;
niezależny przegląd kodu i konfiguracji.

Niski poziom reszt solvera nie zastępuje kontroli siatki, domeny ani bilansu wielkości ocenianej.

Walidacja

Walidacja pyta, czy model reprezentuje rzeczywistość wystarczająco dobrze dla zamierzonego użycia. Metrykę należy określić przed obejrzeniem wyniku. Może to być:

$$E_y=y_\mathrm{model}-y_\mathrm{test},$$

znormalizowany błąd:

$$E_N = \frac{y_\mathrm{model}-y_\mathrm{test}} {\sqrt{U_\mathrm{model}^2+U_\mathrm{test}^2}},$$

albo wielowymiarowa metryka profilu i czasu. Próg akceptacji zależy od konsekwencji decyzji.

Zbiór użyty do kalibracji nie może pełnić roli jedynego zbioru walidacyjnego. Poprawna praktyka rozdziela co najmniej:

dane do identyfikacji;
dane do walidacji;
dane do końcowego testu predykcyjnego.

Domena walidacji

Domena powinna być opisana współrzędnymi fizycznymi, a nie zdaniem „model działa w zakresie prób”. Obejmuje:

zakres $N_c$ i $\dot m_c$;
spręż i temperaturę;
warunki otoczenia;
stan ustalony albo przejściowy;
konfigurację wlotu i dyszy;
geometrię oraz stan powierzchni;
poziom degradacji;
typ zaburzenia;
wersję sterowania;
właściwości gazu.

Interpolacja wewnątrz tej domeny i ekstrapolacja poza nią powinny mieć różne statusy.

Kwalifikacja komponentu, silnika i zabudowy

Komponent

Dowód komponentowy obejmuje mapę, granice, obciążenia termiczne i mechaniczne, tolerancje, trwałość oraz stan po próbie. Nie każdy element tej kwalifikacji może zostać wyprowadzony z badań przepływowych. Udana mapa sprężarki nie zastępuje dowodu integralności wirnika, a udana charakterystyka turbiny nie zastępuje oceny życia cieplnego.

Silnik zintegrowany

Na poziomie kompletnej drogi gazowej trzeba wykazać:

zgodność masy i energii;
bilans wału;
osiągi w obwiedni;
działanie poza punktem projektowym;
stany przejściowe;
współpracę ze sterowaniem;
marginesy do wszystkich granic;
wpływ tolerancji i degradacji;
zachowanie po ekspozycjach środowiskowych;
identyfikowalność konfiguracji.

Próba kompletnego silnika nie daje automatycznie dokładnej mapy każdego komponentu, ponieważ wiele parametrów jest nieobserwowalnych albo współzależnych. Z kolei próby komponentów nie odtwarzają pełnych warunków integracji. Potrzebna jest piramida dowodów.

Zabudowa

Wlot, korpus platformy, sąsiednie powierzchnie, kanał wylotowy oraz stan cieplny po zabudowie zmieniają warunki graniczne. Model silnika z jednorodnym wlotem nie dowodzi kompatybilności z zabudowanym wlotem. Próba stanowiskowa z wolnym wylotem nie dowodzi zachowania w zamkniętej przestrzeni lub przy innym przeciwciśnieniu.

Badanie po zabudowie powinno rozdzielić:

odzysk i zaburzenie wlotowe;
opór instalacyjny;
przeciwciśnienie wydechu;
oddziaływanie cieplne;
drgania;
wpływ elementów pomiarowych;
różnice między stanowiskiem naziemnym a warunkami lotu.

Trwałość i środowisko

Cykl trwałościowy powinien reprezentować liczbę uruchomień, przejścia, postoje, wyłączenia i cykle cieplne wynikające z profilu życia. Długi bieg w jednym punkcie nie zastępuje wielu krótkich cykli. Wynik obejmuje zarówno dryft osiągów, jak i inspekcję powykonawczą.

Po magazynowaniu, wibracji, wilgoci i temperaturze napęd musi nadal spełniać wymagania. Krótka kontrola funkcjonalna nie dowodzi zachowania całej mapy. Każde znalezisko z inspekcji powinno mieć lokalizację, metodę, rozmiar, prawdopodobieństwo wykrycia i decyzję techniczną.

Konfiguracja i cyfrowy ślad dowodu

Wynik modelu powinien wskazywać jednoznacznie:

wersję modelu systemowego;
wersje map;
źródła i sumy kontrolne danych;
geometrię CAD i zasady uproszczeń;
siatkę oraz ustawienia solvera;
własności gazu i materiałów;
model sterowania;
skrypt redukcji;
kalibracje instrumentów;
egzemplarz testowy;
konfigurację stanowiska;
kampanię i przebieg;
zbiór danych użyty do strojenia;
zbiór walidacyjny;
kryteria decyzji.

Przykładowy rekord wyniku:

analysis_id: GASPATH-QOI-017
system_model: SYS-4.2.1
compressor_map: CMP-R3
turbine_map: TUR-R5
geometry: ENG-E04
controller: CTRL-2.7
test_configuration: CELL-B17
reduction_script: RED-1.9
calibration_set: CAL-2026-04
validation_set: VAL-2026-07
qoi: minimum_transient_margin
domain: DOM-TR-03
decision_rule: DR-12

Nazwa final_map.csv nie jest kontrolą konfiguracji. Po aktualizacji modelu należy zachować poprzednią wersję, uzasadnienie zmiany, listę wyników dotkniętych zmianą oraz regresję.

Pełna ponowna walidacja jest potrzebna, gdy zmienia się:

zasadnicza fizyka modelu;
definicja stacji lub wielkości ocenianej;
architektura komponentu;
kluczowa mapa;
model zaburzenia;
kryterium granicy;
model sterowania;
klasa bezpieczeństwa;
domena zamierzonego użycia.

Drobna korekta kodu może podlegać regresji różnicowej, jeśli dane, teza i konfiguracja fizyczna pozostają niezmienione.

Końcowa karta oceny

Końcowy przegląd nie powinien powtarzać rozdziałów w formie długiej ankiety. Osiem bramek wiąże pytanie z artefaktem, który można sprawdzić i zachować:

Bramka	Pytanie rozstrzygające	Wymagany artefakt	Warunek zatrzymania decyzji
1. granica i stacje	czy model, CFD i próba opisują te same objętości oraz przekroje?	mapa granic, stacji, dopływów, upustów, przecieków i powrotów	niezgodne definicje fizyczne
2. definicje i redukcja	czy wielkości statyczne, całkowite, uśrednione i skorygowane mają jedną definicję?	słownik wielkości, warunki odniesienia i wersjonowany skrypt redukcji	wynik zmienia się wraz z nieudokumentowaną metodą uśredniania
3. mapy i konfiguracja	skąd pochodzi każda mapa i gdzie wolno ją interpolować lub ekstrapolować?	surowe punkty, identyfikator geometrii, flagi domeny i niepewność postaci	trajektoria używa komórek niedopuszczonych albo innego wariantu
4. zamknięcie fizyczne	czy jednocześnie zamknięto masę, energię, moment wału, ciśnienia i odpływ?	wektor niezależnych reszt oraz kontrola praw zachowania po zbieżności	mała reszta solvera przy niezbilansowanej wielkości fizycznej
5. dynamika i sterowanie	czy model rozdziela linię ustaloną od trajektorii wału, objętości, sensorów i aktuatorów?	przebiegi czasowe, jawny model zaburzenia i zapasy nominalny, przejściowy oraz degradacyjny	decyzja przejściowa oparta wyłącznie na mapie ustalonej
6. pomiar i niepewność	czy wynik pochodny ma pełny model pomiaru, wspólny czas, blokadę sond i kowariancje?	budżet niepewności od surowego kanału do wielkości ocenianej	brak identyfikowalności dominującego wkładu albo korelacji
7. weryfikacja i walidacja	czy testy implementacji są oddzielone od niezależnego sprawdzenia rzeczywistości w określonej domenie?	przypadki analityczne, badanie siatki i kroku, rozdział danych kalibracyjnych i walidacyjnych, metryka ustalona z góry	kalibracja i walidacja używają jedynego wspólnego zbioru
8. decyzja i zmiana	czy zamierzone użycie, kryterium, populacja, ograniczenia i wyzwalacze ponownej walidacji są jawne?	macierz dowodów, rekord konfiguracji, reguła decyzji i pakiet regresyjny	decydent otrzymuje liczbę bez domeny, niepewności lub pochodzenia

Niespełnienie bramki nie dowodzi automatycznie, że wszystkie obliczenia są błędne. Oznacza, że wynik nie ma jeszcze kompletnego prawa do wskazanej decyzji. Dzięki temu brak nie ginie pod estetyką mapy ani ogólnym stwierdzeniem o „dobrej zgodności”, lecz ma właściciela, artefakt zamykający i skutek dla zakresu użycia.

Bezpieczeństwo prób i stan bezpieczny

Stanowisko sprężarki lub turbiny magazynuje znaczną energię obrotową. Gorący silnik dodaje paliwo, płomień, ciśnienie, toksyczne produkty spalania, hałas i możliwość uwolnienia części wirnika. Próby prowadzi wyłącznie uprawniona organizacja w zatwierdzonym obiekcie z zabezpieczeniem odłamkowym, zdalnym sterowaniem, blokadami, strefą wyłączenia, ochroną przeciwpożarową i planem reagowania właściwymi dla konkretnego urządzenia.

Polecenie wyłączenia nie oznacza potwierdzonego stanu bezpiecznego. Dostęp może nastąpić dopiero według zatwierdzonej procedury, po potwierdzeniu między innymi zatrzymania wirnika, izolacji paliwa i pozostałych źródeł energii oraz bezpiecznych ciśnień i temperatur. Nietypowy dźwięk, drgania, zapach, wyciek albo utrata telemetrii wymagają zachowania stanu wyłączenia i reakcji personelu stanowiska, nie samodzielnej diagnostyki przy urządzeniu.

Model nie może być jedyną podstawą uznania obiektu za bezpieczny. Nie obserwuje niezamodelowanych pęknięć, obcego ciała, błędu montażu, utraty sygnału ani niezgodnej konfiguracji.

Synteza

Dopasowanie drogi gazowej zamienia charakterystyki komponentów w rozwiązanie kompletnego silnika. Jego rdzeniem nie jest pojedyncza mapa, lecz jednoczesne spełnienie ciągłości masy, bilansu energii, bilansu momentu wału, zgodności ciśnień oraz zdolności przepływowej turbiny i dyszy.

Pierwsze studium wykazało, że nominalny punkt wymaga jawnej księgi przepływów, iteracji wszystkich reszt i propagacji niepewności do ciągu oraz zapasu. Drugie pokazało, dlaczego trajektoria przejściowa przy zaburzeniu wlotowym nie może być zastąpiona linią stanów ustalonych. Trzecie wykazało, że rozbieżności stanowisko–CFD–silnik należy najpierw badać przez definicje stacji, uśrednianie, blokadę, wersję mapy i wpływ stanowiska, a dopiero później przez kalibrację fizyki modelu.

Dojrzała kwalifikacja zachowuje profile zamiast samych średnich, trajektorie zamiast samych punktów, rozkład populacji zamiast jednego nominalnego egzemplarza oraz historię konfiguracji zamiast pliku nazwanego „final”. Każdy wynik jest przekazywany razem z pochodzeniem danych, niepewnością, domeną walidacji, kryterium decyzji i ograniczeniami użycia.

Najważniejszym produktem analizy nie jest maksymalna wyspa sprawności ani najmniejszy błąd dopasowania. Jest nim kontrolowana obwiednia pracy, w której wiadomo, jakie prawa zostały domknięte, co rzeczywiście zmierzono, czego model nie obejmuje, jaki pozostaje margines i jaka zmiana wymusi ponowne otwarcie dowodu.

Źródła

Przypisy

Thomas M. Lavelle, An Introduction to Thermodynamic Performance Analysis of Aircraft Gas Turbine Engine Cycles Using the Numerical Propulsion System Simulation Code, NASA/TM-2007-214690, 2007, dodatek C „Turbomachinery Maps”, NASA Technical Reports Server, dostęp 2026-09-10. mapy, wielkości skorygowane i skalowanie. Źródło ogólne, niedotyczące konkretnego współczesnego napędu wojskowego.
↩
NASA Glenn Research Center, Nancy Hall, Compressor–Turbine Matching, aktualizacja 13 maja 2021 r., NASA Glenn, dostęp 2026-09-10. bilans pracy sprężarki i turbiny. Nie jest normą kwalifikacyjną.
↩
NASA Glenn Research Center, Thrust Equation, NASA Glenn, dostęp 2026-09-10. źródło ogólne.
↩
H. C. Melick Jr., Analysis of Inlet Flow Distortion and Turbulence Effects on Compressor Stability, NASA-CR-114577, 1973, NASA Technical Reports Server, dostęp 2026-09-10. wpływ obwodowego zaburzenia ciśnienia całkowitego na stateczność sprężarki. Źródło historyczne, nie mapa współczesnego napędu.
↩
Marian Gieras, Miniaturowe silniki turboodrzutowe, wyd. 2, 2025, zwłaszcza s. 72–160 i 169–185. sprężarka promieniowa, dyfuzor, komora, turbina, CFD i stanowiska małych silników. Dane GTM-120 dotyczą cywilnego obiektu historycznego i nie są użyte jako parametry studiów syntetycznych.
↩1 ↩2
NASA, Navy/NASA Engine Program NNEP89 User’s Manual, 1991, NASA Technical Reports Server — dokument 19910020827, dostęp 2026-09-10. błędy zgodności przepływu i pracy oraz dobór zmiennych iterowanych. Źródło historyczne.
↩
Joint Committee for Guides in Metrology, Evaluation of Measurement Data — Guide to the Expression of Uncertainty in Measurement, JCGM 100:2008, DOI 10.59161/JCGM100-2008E, BIPM, dostęp 2026-09-10. propagacja niepewności, współczynniki wrażliwości i kowariancje. BIPM publikuje także JCGM 101:2008 dotyczący propagacji rozkładów metodą Monte Carlo.
↩1 ↩2
Mahmood Abdelwahab, Thomas J. Biesiadny, Dean Silver, Measurement Uncertainty for the Uniform Engine Testing Program Conducted at NASA Lewis Research Center, NASA-TM-88943, 1987, NASA Technical Reports Server, dostęp 2026-09-10. niepewność strumienia wlotowego, ciągu netto i jednostkowego zużycia paliwa. Źródło historyczne dotyczące programu porównań stanowisk.
↩
Jeffrey T. Csank, Dynamic Systems Analysis for Turbine Based Aero Propulsion Systems, NASA/TM-2016-219133, 2016, NASA Technical Reports Server, dostęp 2026-09-10. analiza osiągów przejściowych, sterowania i zapasu sprężarki. Wyniki konkretnego modelu turbofanu nie są przenoszone ilościowo do miniaturowego turbojetu.
↩
Mark G. Turner, John A. Reed, Robert Ryder, Joseph P. Veres, Multi-Fidelity Simulation of a Turbofan Engine with Results Zoomed into Mini-Maps for a Zero-D Cycle Simulation, NASA/TM-2004-213076, 2004, NASA Technical Reports Server, dostęp 2026-09-10. sprzężenie modelu cyklu NPSS z modelami trójwymiarowymi i mapami częściowymi. Dotyczy demonstracji turbofanu GE90-94B, nie miniaturowego napędu wojskowego.
↩
Piotr Tarnawski, Symulacyjne badania koncepcyjne nietypowego silnika turbowałowego pracującego w obiegu termodynamicznym Humphreya, 2024, zwłaszcza rozdz. 2–3 i 8. bilanse obiegu, CFD oraz jawne ograniczenia wnioskowania bez próby stanowiskowej. Źródło dotyczy odmiennej, koncepcyjnej architektury turbowałowej.
↩
Mateusz Turkowski, Metrologia przepływów, 2018, s. 253–266. model pomiaru, oceny typu A i B, współczynniki wrażliwości, niepewność złożona oraz CFD w metrologii przepływów. Dane katalogowe przyrządów z książki nie są traktowane jako aktualne wymagania.
↩
NASA Office of the Chief Engineer, Standard for Models and Simulations, NASA-STD-7009B, 5 marca 2024 r., NASA Technical Standards System, dostęp 2026-09-10. cykl życia modelu, kryteria akceptacji i komunikowanie wiarygodności wyników. Standard nie jest specyfikacją kwalifikacyjną konkretnego wojskowego systemu napędowego.
↩