Pressure-gain combustion i obieg Humphreya

Spalanie z przyrostem ciśnienia (pressure-gain combustion, PGC) obiecuje coś, czego klasyczna komora turbiny gazowej zwykle nie daje: doprowadzenie ciepła bez utraty całkowitego ciśnienia, a w idealizacji nawet z jego wzrostem. Sama fala wysokiego ciśnienia nie jest jednak jeszcze zyskiem napędu. Trzeba wykazać, że większa zdolność gazu do wykonania pracy przetrwała napełnianie, izolację wlotu, opróżnianie, mieszanie, współpracę ze sprężarką, turbiną lub dyszą, chłodzenie oraz koszt masy i trwałości.

Teza artykułu brzmi: przewaga PGC nie jest właściwością samego procesu spalania, lecz wynikiem bilansu egzergii całego niestacjonarnego układu. Wywód zaczyna się więc od idealnych obiegów Braytona–Joule’a i Humphreya, następnie wprowadza bilans okresowej objętości kontrolnej, a na końcu sprawdza, jak idealny zysk przechodzi — albo nie przechodzi — przez rzeczywiste interfejsy.

Artykuł służy analizie termodynamicznej i ocenie dowodów. Nie zawiera geometrii komór, składu mieszanin, warunków inicjacji detonacji, nastaw wtrysku, sekwencji zaworów ani procedury uruchomienia urządzenia.

Bezpieczeństwo: stanowiska PGC łączą gorące gazy, gwałtowne zmiany ciśnienia, hałas impulsowy, toksyczne produkty spalania i często szybko wirujące części. Próby wymagają uprawnionej organizacji, zdalnej obsługi, osłon, strefy wyłączenia i zatwierdzonej analizy zagrożeń. Zanik dźwięku nie oznacza stanu bezpiecznego: mogą pozostać ciśnienie, temperatura, reagenty i energia elektryczna.

Pięć poziomów, których nie wolno mieszać

W dyskusji o PGC słowo „zysk” bywa odnoszone do różnych granic. Na poziomie fali można obserwować lokalny wzrost ciśnienia statycznego. Na poziomie porcji gazu pyta się, jak zmieniły się entalpia, entropia i zdolność do wykonania pracy. Na poziomie komory dochodzą straty napełniania, przecieków i opróżniania. Na poziomie urządzenia granica obejmuje wlot, komorę wyrównawczą (plenum), mieszalnik, turbinę albo dyszę. Dopiero poziom napędu i misji uwzględnia masę, chłodzenie, sterowanie, opór, niezawodność i zużycie paliwa.

Twierdzenie może być prawdziwe na jednym poziomie i fałszywe na kolejnym. Dodatnie ciśnienie za czołem fali nie dowodzi dodatniego średniego ciśnienia całkowitego na wylocie. Dodatni bilans komory nie dowodzi większej mocy wału. Większa moc stanowiskowa nie dowodzi większego zasięgu, jeśli urządzenie wymaga ciężkiej komory wyrównawczej i intensywnego chłodzenia.

Dlatego każde porównanie powinno podawać: powierzchnię kontrolną, stan odniesienia, sposób uśredniania, wspólne ograniczenia oraz to, czy wynik pochodzi z obiegu idealnego, modelu, demonstratora czy kompletnego zespołu.

Brayton–Joule jako punkt odniesienia

Idealny zamknięty obieg Braytona–Joule’a składa się z izentropowego sprężania, izobarycznego doprowadzenia ciepła, izentropowego rozprężania i izobarycznego odprowadzenia ciepła. Nie jest dokładnym modelem otwartego silnika, ale porządkuje wpływ sprężu i temperatury.

Przy stałych właściwościach gazu zdefiniujmy

$$a=\frac{\gamma-1}{\gamma}, \qquad \tau=r_p^a=\frac{T_2}{T_1},$$

gdzie $r_p=p_2/p_1$ jest sprężem, a $\gamma=c_p/c_v$. Dla idealnego sprężania i rozprężania sprawność cieplna wynosi

$$\eta_B=1-\frac{1}{\tau} =1-r_p^{-(\gamma-1)/\gamma}.$$

W realnej turbinie gazowej komora nazywana „stałociśnieniową” ma spadek ciśnienia całkowitego. Ciepło jest doprowadzane podczas przepływu, występują tarcie, mieszanie, chłodzenie ścian i nierównomierność profilu. Sprężarka musi dostarczyć pracę, turbina pokrywa jej zapotrzebowanie, a dopiero nadwyżka może zasilać wał lub dyszę. Monografia Mariana Gierasa pokazuje ten pełny łańcuch dla małego silnika: obieg → modele 1D i CFD → hamownia → walidacja → zmiana konstrukcji.¹

Brayton jest zatem nie tylko konkurencyjnym obiegiem. Jest także dyscypliną bilansową: każda obietnica PGC powinna zostać porównana z konwencjonalnym układem o tych samych warunkach wlotowych, doprowadzonym cieple, limitach temperatury i definicji produktu użytecznego.

Wyprowadzenie idealnego obiegu Humphreya

Obieg Humphreya zastępuje izobaryczne doprowadzenie ciepła przemianą izochoryczną. Po izentropowym sprężeniu $1\rightarrow2$ zamknięta porcja gazu otrzymuje ciepło przy stałej objętości $2\rightarrow3$. Ciśnienie rośnie wraz z temperaturą. Następnie gaz rozpręża się izentropowo $3\rightarrow4$ do ciśnienia początkowego, a ciepło jest odprowadzane izobarycznie $4\rightarrow1$.

Niech

$$\theta=\frac{T_3}{T_2}=\frac{p_3}{p_2}.$$

Dla idealnego gazu temperatura po rozprężeniu do $p_1$ spełnia

$$\frac{T_4}{T_1}=\theta^{1/\gamma}.$$

Ciepło doprowadzone i odprowadzone wynoszą odpowiednio

$$q_{in,H}=c_vT_1\tau(\theta-1),$$

$$q_{out,H}=c_pT_1\left(\theta^{1/\gamma}-1\right).$$

Stąd idealna sprawność Humphreya:

$$\eta_H =1-\frac{q_{out,H}}{q_{in,H}} =1-\frac{\gamma\left(\theta^{1/\gamma}-1\right)} {\tau(\theta-1)}.$$

Różnica względem Braytona wynika z drogi procesu, a nie z „darmowej energii”. Przy stałej objętości doprowadzenie ciepła podnosi ciśnienie przed rozprężaniem, ale rzeczywiste urządzenie musi najpierw zamknąć lub gazodynamicznie odizolować porcję, a potem zamienić jej okresowe opróżnianie na użyteczny strumień.

Na wykresie $p$–$v$ różnicę widać w odcinku $2\rightarrow3$: izobara Braytona przesuwa stan w prawo, natomiast izochora Humphreya podnosi ciśnienie pionowo. Na wykresie $T$–$s$ oba idealne odcinki sprężania i rozprężania są izentropowe, lecz przyrost entropii podczas doprowadzania tej samej relacji temperatur jest mniejszy dla przemiany izochorycznej, ponieważ występuje $c_v$, a nie $c_p$. Diagram nie rozstrzyga jeszcze wyniku, bo przy wspólnym $T_3$ obie drogi pobierają różne ilości ciepła.

Uczciwy protokół porównania

Wynik zależy od tego, co uznano za wspólne. Można narzucić ten sam spręż, tę samą temperaturę maksymalną, tę samą ilość ciepła, ten sam stosunek paliwa do powietrza albo tę samą moc. Nie są to równoważne protokoły.

Rozważmy neutralny przykład gazu o $\gamma=1{,}4$, sprężu $r_p=4$ i stosunku temperatur $\theta=2$. Wtedy $\tau\approx1{,}486$. Idealne równania dają

$$\eta_B\approx0{,}327, \qquad \eta_H\approx0{,}396.$$

Różnica wynosi około $6{,}9$ punktu procentowego. Nie jest prognozą silnika. Przy tym samym $r_p$ i $\theta$ oba obiegi przyjmują inną ilość ciepła: w Braytonie występuje $c_p$, a w Humphreyu $c_v$. Znormalizowana praca netto przy tych właśnie ograniczeniach wynosi

$$\frac{w_B}{c_vT_1} =\gamma(\theta-1)(\tau-1) \approx0{,}680,$$

$$\frac{w_H}{c_vT_1} =\tau(\theta-1) -\gamma\left(\theta^{1/\gamma}-1\right) \approx0{,}589.$$

Humphrey ma w tym przykładzie większą sprawność, lecz mniejszą pracę na przyjętą jednostkę $c_vT_1$, ponieważ przy stałym $\theta$ przyjął mniej ciepła. To prosty dowód, że ranking „lepszy–gorszy” bez protokołu porównania jest bez znaczenia.

W realnym studium należy dodatkowo zrównać sprawności sprężarki i turbiny, limit temperatury materiału, straty ciśnienia, przepływy chłodzące, moc urządzeń pomocniczych oraz stan otoczenia. Wyniki idealne pozostają górną wskazówką, nie obietnicą osiągów.

Ciśnienie całkowite, egzergia i znaczenie przyrostu

Ciśnienie statyczne opisuje lokalne obciążenie termodynamiczne, a ciśnienie całkowite — stan po hipotetycznym izentropowym wyhamowaniu jednorodnego strumienia. W PGC przepływ może jednocześnie zmieniać się w czasie i przestrzeni, mieć składowe poprzeczne oraz niejednorodny skład. Nie istnieje wtedy jeden oczywisty „średni” stan całkowity.

Bardziej fundamentalną miarą jest właściwa egzergia przepływu względem otoczenia:

$$\psi=(h-h_0)-T_0(s-s_0) +\frac{V^2}{2}+gz.$$

Egzergia mierzy maksymalną pracę użyteczną możliwą przy doprowadzeniu strumienia do równowagi z otoczeniem. PGC dodaje energię chemiczną, lecz zarazem niszczy część egzergii przez nieodwracalność reakcji, fal, mieszania, tarcia i wymiany ciepła.

NASA opisuje pressure gain jako praktyczną miarę dostępności płynu do wytworzenia pracy lub ciągu, podkreślając trudność pomiaru w gorącym i niejednorodnym wylocie.² Poprawna publikacja musi zatem powiedzieć, czy raportuje lokalny pik, średnie ciśnienie statyczne, średnie ciśnienie całkowite, stan po wymieszaniu, strumień egzergii czy rzeczywistą pracę turbiny.

Mieszanie niszczy część zdolności do pracy

Neutralny rachunek można wykonać bez modelu konkretnego urządzenia. Dla dwóch nieruchomych strug tego samego idealnego gazu, przy wspólnym ciśnieniu otoczenia $p_0$, pomińmy energię kinetyczną i przyjmijmy stałe $c_p$. Temperaturowa część bezwymiarowej egzergii wynosi wtedy

$$\psi_T^*=\frac{\psi_T}{c_pT_0} =\frac{T}{T_0}-1-\ln\!\left(\frac{T}{T_0}\right).$$

Niech struga A niesie 25% masy i ma $T_A/T_0=1{,}6$, a struga B — 75% masy i $T_B/T_0=1{,}1$. Przed wymieszaniem średnia egzergia właściwa na jednostkę całkowitej masy wynosi

$$\bar\psi_{przed}^* =0{,}25(0{,}6-\ln1{,}6) +0{,}75(0{,}1-\ln1{,}1) \approx0{,}0360.$$

Bilans entalpii daje po wymieszaniu $T_m/T_0=0{,}25(1{,}6)+0{,}75(1{,}1)=1{,}225$. Dla jednorodnego stanu końcowego

$$\psi_{po}^*=0{,}225-\ln1{,}225\approx0{,}0221.$$

Różnica około $0{,}0140\,c_pT_0$ została zniszczona przez nieodwracalne mieszanie, mimo że energia została zachowana. Stan o poprawnej średniej temperaturze nie zachowuje więc średniej zdolności niejednorodnego strumienia do wykonania pracy.

Ta sama nieliniowość pojawia się w uśrednianiu przebiegu. Rozważmy cztery równe przedziały czasu z umownymi parami $(\dot m,T/T_0)$ równymi $(1;1{,}0)$, $(2;1{,}3)$, $(0{,}5;1{,}6)$ i $(1{,}5;1{,}1)$. Otrzymujemy:

Operator	Wynik	Znaczenie
średnia arytmetyczna temperatury	1,250	każdy moment ma tę samą wagę, niezależnie od transportowanej masy
średnia masowo-strumieniowa temperatury	1,210	zachowuje średni strumień entalpii przy stałym (c_p)
średnia masowo-strumieniowa ψ*	0,0295	uśrednia zdolność każdej porcji do wykonania pracy
ψ* obliczona ze średniej temperatury 1,210	0,0194	traci informację o niejednorodności przed zastosowaniem funkcji nieliniowej

Ostatnie dwie wartości różnią się o około 34%. Nie oznacza to, że jeden operator jest zawsze właściwy: średnia temperatury służy bilansowi entalpii, a średnia egzergii — bilansowi zdolności do pracy. Błąd polega na użyciu operatora bez wskazania zachowywanej wielkości. W rzeczywistym przekroju dochodzą ciśnienie, prędkość, skład i korelacje między nimi.

Bilans okresowej objętości kontrolnej

PGC osiąga nie stan ustalony, lecz powtarzalny cykl graniczny. Chwilowa masa $m$ i energia $E$ w objętości kontrolnej spełniają

$$\frac{dm}{dt}=\sum\dot m_{in}-\sum\dot m_{out},$$

$$\frac{dE}{dt} =\sum\dot m_{in}\left(h+\frac{V^2}{2}\right) -\sum\dot m_{out}\left(h+\frac{V^2}{2}\right) +\dot Q-\dot W.$$

Po scałkowaniu przez okres $\mathcal T$ w ustalonym cyklu granicznym powinno być

$$\Delta m_{\mathcal T}=0, \qquad \Delta E_{\mathcal T}=0.$$

Jeżeli te domknięcia nie zachodzą, średnia może obejmować rozruch, dryft cieplny albo błąd modelu. Samo powtarzanie się częstotliwości fali nie dowodzi okresowości masy, energii i stanu ścian.

Dla wielkości $\phi$ przenoszonej przez przekrój sensowna średnia strumieniowa ma postać

$$\langle\phi\rangle_m= \frac{ \int_0^{\mathcal T}\!\int_A \rho u_n\phi\,dA\,dt }{ \int_0^{\mathcal T}\!\int_A \rho u_n\,dA\,dt }.$$

Średnia arytmetyczna z czujnika nie musi jej odtwarzać. Co więcej, obliczenie turbiny ze średnich wejść nie musi dać średniej pracy turbiny, ponieważ charakterystyka maszyny jest nieliniowa. Raport NASA/TM-2012-217443 wykazał na modelu, że sposób uśredniania wylotu PGC istotnie wpływa na przewidywaną pracę, a miary związane ze strumieniem pędu lepiej odtwarzały wynik okresowy niż część prostszych średnich.³

Rodziny PGC rozwiązują różne problemy

PGC nie jest synonimem detonacji. Wspólnym celem jest ograniczenie rozprężania podczas doprowadzania ciepła, lecz architektury różnią się mechanizmem izolacji i charakterem wylotu.

Rodzina	Mechanizm organizujący proces	Charakter wylotu	Dominujący problem integracji
komora pulsacyjna rezonansowa	sprzężenie reakcji z akustyką i elementem zwrotnym	okresowy w czasie	zawór, hałas, mieszanie i trwałość
wirnik falowy	fale sprężania i rozprężania w obracających się kanałach	okresowy względem portów	synchronizacja portów, przecieki i straty końcowe
PDE	kolejne napełnianie, reakcja falowa i opróżnianie komór	silnie impulsowy	izolacja wlotu, przepłukanie i wygładzenie
RDE	fala obiega komorę, a dopływ i odpływ są quasi-ciągłe	niejednorodny w czasie i obwodzie	wlot o dużej diodowości, tryby falowe i turbina
cykliczna komora mechaniczna	fizyczne odcinanie porcji i etapowe rozprężanie	okresowy, zależny od rozrządu	uszczelnienie, tarcie, bezwładność i chłodzenie

RDE zasilana powietrzem i sprężarką nie ma tej samej granicy systemu co rakietowa RDRE z przechowywanym utleniaczem. Wspólna fizyka falowa nie upoważnia do przenoszenia sprawności, masy ani dojrzałości między zastosowaniami. Przegląd NASA obejmuje wirniki falowe, komory pulsacyjne, PDE i RDE oraz ich interfejsy z wlotami, dyszami, eżektorami i turbinami.⁴

Komora pulsacyjna: eżektor jest częścią wyniku

Eksperyment NASA Glenn z 2005 r. połączył zaworową komorę pulsacyjną z eżektorem i osłoną, tworząc zespół badany jako komora dla turbiny gazowej.¹⁰ Na granicy całego zespołu raportowano około 3,5% wzrostu ciśnienia całkowitego przy stosunku temperatur odpowiadającym rozważanemu zastosowaniu. Za eżektorem odchylenie standardowe pulsacji ciśnienia całkowitego wynosiło około 5% średniej i było mniejsze niż w samej strudze pulsacyjnej.

Wartość tego przykładu nie polega na przeniesieniu dwóch procentów do innej architektury. Pokazuje on równoczesne działanie dwóch mechanizmów: komora wytwarza niestacjonarny potencjał ciśnieniowy, a eżektor miesza strugę z przepływem obejściowym, zmniejszając niejednorodność widzianą przez odbiornik. To wygładzenie nie jest darmowe, ponieważ mieszanie niszczy część egzergii. Granicę należy więc postawić za eżektorem, jeżeli to on czyni przepływ użytecznym dla turbiny.

Przypadek pulsacyjny równoważy studium RDE z dalszej części artykułu. W obu eksperymentach lokalny proces i wynik zintegrowanego urządzenia są odrębnymi twierdzeniami, ale skutek integracji może mieć przeciwny znak: dobrze dobrany eżektor zachowuje niewielki dodatni bilans, podczas gdy ograniczenie wlotu i długość badanego układu RDE mogą zużyć lokalny przyrost w całości.

Pięć faz tworzy jeden sprzężony cykl

W urządzeniu porcjowym można wyróżnić napełnianie, izolację, doprowadzenie ciepła, opróżnianie i przepłukanie. Udziały czasu spełniają

$$f_A=f_{nap}+f_{izol}+f_{ciep}.$$

Łącznie z opróżnianiem i przepłukaniem udziały te wypełniają cały okres:

$$f_A+f_{wyl}+f_{płuk}=1.$$

Nie można poprawić jednej fazy w izolacji. Dłuższe napełnianie może zwiększyć masę porcji, ale zmniejszyć częstotliwość cyklu. Mocniejsza izolacja wlotu może ograniczyć przepływ zwrotny, lecz zwiększyć stratę w kierunku roboczym. Szybsze opróżnianie zwiększa chwilowy strumień i obciążenie odbiornika. Przepłukanie zmniejsza ryzyko niepożądanego zapłonu pozostałości, ale zużywa medium i obniża udział czasu roboczego.

Bilans masy jednego okresu można zapisać jako

$$m_{nap}=m_{wyl}+m_{strat}+\Delta m.$$

W cyklu granicznym $\Delta m=0$. Jeżeli model zakłada $m_{strat}=0$, idealną izolację i natychmiastowe przełączenie faz, to oszacował potencjał procesu, a nie działanie mechanizmu. Te założenia powinny występować w tabeli bilansu obok wyniku, nie w odległym przypisie.

W RDE fazy nie są rozdzielone globalnie w czasie. Współistnieją w różnych miejscach obwodu, a liczba i prędkość fal mogą się zmieniać. Nadal obowiązuje jednak ten sam wymóg: po pełnym okresie obserwacji masa, energia, modalność i obciążenie ścian muszą być opisane wspólnym bilansem.

Deflagracja, detonacja i przełączanie trybów

Deflagracja rozchodzi się głównie dzięki transportowi ciepła i składników, detonacja łączy reakcję z falą uderzeniową. Nazwanie sygnału ciśnienia „detonacją” wymaga dowodu zgodnego z rozdzielczością czasową i przestrzenną pomiaru. Pojedynczy pik nie rozstrzyga mechanizmu.

Układ może pracować z jedną albo wieloma falami, falą przeciwnie biegnącą, nieregularnym spalaniem lub okresowym wygasaniem. Dwa przebiegi o podobnej średniej mogą mieć inną amplitudę, fazę i widmo, a więc inaczej obciążać wlot, turbinę i ścianę.

Stan cyklu warto opisywać wektorem

$$x=[\bar p,\bar T,\dot m,N_w,A_p,\varphi,\Theta_w],$$

gdzie $N_w$ oznacza liczbę rozpoznanych fal, $A_p$ amplitudę pulsacji, $\varphi$ relację fazową, a $\Theta_w$ stan cieplny ściany. Częstotliwość sama nie wystarcza do identyfikacji trybu.

Wlot i sprężarka: izolacja nie może być korkiem

Sprężarka pracuje na mapie zależnej od skorygowanego przepływu, prędkości i sprężu. Pulsujące przeciwciśnienie może okresowo przesuwać punkt pracy ku niestateczności, nawet gdy średnia leży w dozwolonym obszarze. Komora wyrównawcza tłumi część pulsacji, lecz dodaje objętość i masę oraz może utworzyć własny rezonans.

Idealny wlot PGC ma małą stratę przy przepływie do komory i dużą odporność na przepływ zwrotny. Tę asymetrię można nazwać diodowością przepływową. Raport NASA o wpływie cofki na RDE pokazuje, że przepływ zwrotny może jednocześnie pogarszać osiągi, wzbudzać oscylacje wysokości fali i doprowadzić do jej zaniku; zdolność wlotu do blokowania cofki pozostaje więc częścią wyniku PGC, a nie pomocniczym detalem.⁵

Łańcuch przyczynowy ma postać

$$\text{fala ciśnienia} \rightarrow\text{cofka lub modulacja przepływu} \rightarrow\text{ruch punktu sprężarki} \rightarrow\text{zmiana masy porcji} \rightarrow\text{zmiana następnej fali}.$$

To sprzężenie może zamknąć pętlę dodatnią. Model samej komory z narzuconym stałym dopływem nie wykryje ryzyka, bo usunął z granicy systemu element, który je tworzy.

Turbina i dysza odbierają strumień, nie średnią tabelę

Na wlocie turbiny liczą się nie tylko średnie $p_t$ i $T_t$, ale ich rozkład promieniowy, obwodowy i czasowy, kąt napływu, skład oraz turbulencja. Gorący sektor może cyklicznie obciążać te same łopatki. Zmienny kąt może zwiększać straty profilu, a impuls ciśnienia — wymuszenia łopatkowe.

Chwilową moc stopnia można schematycznie zapisać jako

$$P_t(t)=\dot m(t)\,\Delta h_0(t)\,\eta_t[x(t)],$$

gdzie sprawność zależy od bieżącego stanu. Zwykle

$$\left\langle P_t(t)\right\rangle \ne P_t\!\left(\langle\dot m\rangle, \langle\Delta h_0\rangle, \langle x\rangle\right).$$

Neutralny przebieg czterech równych faz pokazuje skalę błędu. Przyjmijmy znormalizowane $\dot m$, $\Delta h_0$ i współrzędną mapy $q$, a sprawność opiszmy wyłącznie dydaktyczną funkcją

$$\eta_t(q)=0{,}90-0{,}08(q-1)^2.$$

Faza	$\dot m$	$\Delta h_0$	$q$	$\eta_t$	$P_t=\dot m\,\Delta h_0\,\eta_t$
1	1,0	1,0	0,7	0,8928	0,8928
2	1,4	0,8	1,3	0,8928	0,9999
3	0,8	1,3	0,8	0,8968	0,9327
4	1,2	0,9	1,2	0,8968	0,9685

Średnia mocy z przebiegu wynosi $\langle P_t\rangle=0{,}9485$. Tymczasem średnie wejścia są równe $\langle\dot m\rangle=1{,}1$, $\langle\Delta h_0\rangle=1{,}0$ i $\langle q\rangle=1{,}0$, więc podstawienie ich do mapy daje

$$P_t(\langle\dot m\rangle,\langle\Delta h_0\rangle, \langle q\rangle)=1{,}1\cdot1{,}0\cdot0{,}90=0{,}9900.$$

Wynik ze średnich zawyża rzeczywistą średnią o około 4,4%. Powodem jest zarówno nieliniowość mapy sprawności, jak i korelacja strumienia masy ze spadkiem entalpii. Zgodność samych średnich $\dot m$, $\Delta h_0$ i $q$ nie jest zatem wystarczającym warunkiem równoważnego obciążenia turbiny; potrzebne są ich zsynchronizowane przebiegi albo model zachowujący odpowiednie momenty i korelacje.

Dysza także odpowiada na chwilowe ciśnienie, temperaturę i skład. W komorze opróżnianej optymalne dopasowanie zmienia się podczas cyklu. Stała geometria wybiera kompromis; eżektor albo mieszalnik może wygładzić strumień, lecz mieszanie niejednorodnych strug niszczy egzergię.

Pełny wynik należy więc liczyć w co najmniej dwóch postaciach: jako bilans komory oraz jako pracę wału lub ciąg na zewnętrznej granicy urządzenia. Różnica ujawnia koszt interfejsu.

Temperatura, konstrukcja i czas życia

Silniejsza nieustaloność może podnieść teoretyczny zysk, ale równocześnie zwiększa szczytowy strumień ciepła i amplitudę naprężeń. Ściana filtruje szybkie wymuszenie zależnie od dyfuzyjnego czasu cieplnego; średnia temperatura nie pokazuje gradientu ani lokalnego cyklu.

Ścieżka uszkodzenia powinna być jawna:

$$\text{fala i gorący sektor} \rightarrow q''(x,t) \rightarrow T_w(x,t) \rightarrow \sigma(x,t) \rightarrow\text{zmęczenie, pełzanie lub pękanie}.$$

Skumulowane uszkodzenie w neutralnym modelu blokowym można zapisać

$$D=\sum_i\frac{n_i}{N_i}.$$

Równanie Minera nie dowodzi trwałości PGC; przypomina jedynie, że rozkład amplitud i liczba cykli są ważniejsze niż pojedyncza wartość maksymalna. Rzeczywista ocena potrzebuje widma, temperatury, średniego naprężenia, karbów, utleniania, rozrzutu materiału i właściwego modelu pękania.

Chłodzenie zabiera medium sprężone, zmienia lokalny skład i może zakłócać falę. Powłoka termiczna zmniejsza temperaturę podłoża, ale ma własną przyczepność i życie. W urządzeniach mechanicznie cyklicznych dochodzą uszczelnienia, tarcie, zużycie i dynamika wirnika. To nie „problemy wykonawcze po termodynamice”, lecz składniki tego samego bilansu netto.

Sterowanie i diagnostyka muszą rozpoznawać tryb

Sterownik nie obserwuje bezpośrednio całego pola. Dostaje opóźnione i filtrowane sygnały ciśnienia, światła, temperatury, przepływu oraz drgań. Musi odróżnić rozruch, prawidłowy cykl graniczny, przełączenie modów, zanik reakcji, niestabilność wlotu i przeciążenie cieplne.

Minimalna logika diagnostyczna łączy trzy rodzaje informacji:

1. stan procesu — energia pasma, prędkość fazowa, korelacja czujników i powtarzalność okresu;
2. stan urządzenia — temperatura ściany, drgania, przepływ zwrotny i przeciek;
3. stan odbiornika — zapas sprężarki, odpowiedź turbiny, moc albo ciąg.

Sygnał o prawidłowej częstotliwości może pochodzić od niepożądanego trybu. Z kolei chwilowe zniknięcie jednej składowej nie zawsze oznacza wygaśnięcie reakcji. Progi muszą wynikać z walidacji na znanych stanach i uwzględniać niepewność czujnika.

Publiczny opis bezpiecznego wyłączenia powinien kończyć się na celach: odcięcie energii doprowadzanej, zapobieżenie niekontrolowanej reakcji resztkowej, kontrolowane rozładowanie ciśnienia, chłodzenie oraz potwierdzenie stanu. Konkretna kolejność i czasy należą do zatwierdzonej procedury danego stanowiska.

Hierarchia modeli jest łańcuchem pytań

Model powinien być najprostszy, który zachowuje zjawisko rozstrzygające dane pytanie. Zwiększanie liczby komórek nie naprawia złej granicy systemu ani braku walidacji.

Poziom	Pytanie, na które odpowiada	Główne stany	Czego sam nie dowodzi
obieg idealny	czy droga stałej objętości ma potencjał	stany termodynamiczne	strat, masy i wykonalności
model 0D	jak zamyka się masa i energia w okresie	$m(t)$, $E(t)$, $p(t)$	fal i lokalnych obciążeń
model 1D/hybrydowy	jak przebiegają napełnianie i opróżnianie	fale osiowe, strumienie, stan porcji	trójwymiarowej modalności
CFD 2D/3D	jak geometria tworzy fale i niejednorodność	pola przepływu i reakcji	osiągów bez warunków i walidacji
model zredukowany	jak układ zachowuje się w pętli sterowania	wybrane mody i stany cieplne	zjawisk pominiętych przy redukcji
model całego napędu	czy zysk przechodzi na wał, ciąg i misję	mapy komponentów, masa, ograniczenia	trwałości bez danych materiałowych

NASA/TM-2010-216821 jest dobrym przykładem modelu hybrydowego: globalna relacja wyznacza stan po idealizowanym doprowadzeniu ciepła, a prosty model jednowymiarowy śledzi ciągłe opróżnianie i ponowne napełnianie; wylot jest uśredniany do miar przydatnych w analizie korzyści.⁶ Jego siłą jest dopasowanie do pytania, a nie udawanie pełnej chemii i geometrii.

Weryfikacja pyta, czy równania rozwiązano poprawnie: zbieżność czasowa i siatkowa, zachowanie masy i energii, testy kodu oraz czułość algorytmu numerycznego. Walidacja pyta, czy model reprezentuje rzeczywistość w zadanym zastosowaniu: porównuje niezależne pomiary, podaje ich niepewność i zachowuje dane kontrolne niewykorzystane do strojenia.

Zgodność dwóch modeli nie jest automatycznie walidacją, gdy dzielą założenie idealnej szczelności, tę samą korelację strat albo te same warunki brzegowe. Może jedynie potwierdzić spójność ich wspólnego rdzenia.

Metrologia: od czujnika do bilansu

Pulsacyjny przepływ narzuca wymagania dotyczące pasma, odpowiedzi dynamicznej, temperatury, kierunkowości i synchronizacji. Przewód impulsowy może tłumić i przesuwać fazę; czujnik może mieć rezonans; sonda zaburza mały kanał. Kalibracja statyczna nie wystarcza do rekonstrukcji szybkiego przebiegu.

Łańcuch pomiarowy można zapisać jako

$$y(t)=H\{x(t)\}+b(t)+n(t),$$

gdzie $H$ jest odpowiedzią całego toru, $b$ dryftem, a $n$ szumem. Odwracanie $H$ może wzmacniać szum, więc korekcja dynamiczna potrzebuje własnej niepewności.

Minimalny bilans PGC wymaga zsynchronizowanych strumieni masy, energii i pędu na wlocie oraz wylocie, doprowadzonej energii chemicznej, strat ciepła, produktu użytecznego i akumulacji. Ciśnienie mierzone w jednym punkcie nie zamyka powierzchni kontrolnej.

Raport powinien pokazać:

• definicję okresu i kryterium osiągnięcia cyklu granicznego;
• lokalizacje i orientacje czujników oraz ich odpowiedź dynamiczną;
• sposób rekonstrukcji wielkości całkowitych i uśredniania strumieniowego;
• propagację niepewności, brakujące obszary przekroju i korekty modelowe;
• bilans masy oraz energii przed publikacją wskaźnika pressure gain;
• osobno średnią i amplitudę, widmo, fazę oraz zmienność między cyklami.

Optyczna obserwacja fali może wesprzeć identyfikację trybu, termometria ściany — bilans cieplny, a pomiar siły lub momentu — produkt końcowy. Żaden z tych kanałów osobno nie zastępuje pozostałych.

Studium 1: od idealnej przewagi do wyniku netto

Wróćmy do neutralnego przykładu, w którym idealna różnica sprawności Humphreya i Braytona wyniosła $0{,}069$. Zbudujmy nie prognozę urządzenia, lecz dydaktyczny bilans utraty przewagi:

$$\Delta\eta_{net} =\Delta\eta_{ideal} -\delta_{nap} -\delta_{miesz} -\delta_{int} -\delta_{chł} -\delta_{pom}.$$

Przypiszmy umownie $0{,}018$ stratom napełniania i izolacji, $0{,}014$ mieszaniu, $0{,}011$ niedopasowaniu odbiornika, $0{,}009$ chłodzeniu i $0{,}006$ urządzeniom pomocniczym. Pozostaje

$$\Delta\eta_{net}=0{,}011.$$

Około 84% idealnej przewagi zniknęło, choć każda kara z osobna wyglądała umiarkowanie. Liczby są wyłącznie neutralnym przykładem rachunkowym. Ich wartość dydaktyczna polega na pokazaniu, jakie pomiary musiałyby zastąpić symbole.

Same wartości średnie nadal nie wystarczają. Załóżmy następujące standardowe niepewności składników, oszacowane w umownym modelu:

Składnik	Wartość średnia	Niepewność standardowa	Wspólna przyczyna w przykładzie
idealna przewaga	+0,069	0,0020	protokół własności gazu
napełnianie i izolacja	−0,018	0,0025	amplituda pulsacji i cofka
mieszanie	−0,014	0,0020	niejednorodność przekroju
interfejs odbiornika	−0,011	0,0025	amplituda pulsacji i cofka
chłodzenie	−0,009	0,0015	obciążenie urządzeń pomocniczych
urządzenia pomocnicze	−0,006	0,0010	obciążenie chłodzenia

Gdyby błędnie przyjąć pełną niezależność, standardowa niepewność wyniku byłaby równa

$$u_{net,ind}=\sqrt{0{,}0020^2+0{,}0025^2+0{,}0020^2 +0{,}0025^2+0{,}0015^2+0{,}0010^2} \approx0{,}00487.$$

W tym fikcyjnym bilansie straty napełniania i interfejsu rosną wraz z tą samą pulsacją, dlatego przyjmijmy $\rho_{nap,int}=0{,}8$. Koszt chłodzenia i urządzeń pomocniczych ma wspólny napęd pomp, więc niech $\rho_{chł,pom}=0{,}6$. Przy pozostałych korelacjach równych zero wariancja zwiększa się o

$$2(0{,}8)(0{,}0025)(0{,}0025) +2(0{,}6)(0{,}0015)(0{,}0010) =11{,}8\cdot10^{-6},$$

co daje $u_{net,kor}\approx0{,}00596$. Dla $k=2$ wynik ma postać $0{,}011\pm0{,}0119$, więc przedział obejmuje brak korzyści. Model niezależny dawałby $0{,}011\pm0{,}0097$ i pozornie dodatnią dolną granicę. Korelacji nie wolno jednak dobierać po to, by uzyskać pożądany werdykt: muszą wynikać ze wspólnych zmiennych wejściowych, wspólnej aparatury albo danych z prób.

Jeżeli masa dodatkowej komory wyrównawczej zwiększa zużycie energii misji o równoważność większą niż $0{,}011$, przewaga napędu znika mimo dodatniego wyniku termodynamicznego urządzenia. Jeżeli natomiast poprawa wlotu zmniejsza $\delta_{nap}$ bez wzrostu innych kar, zespół wie, gdzie istnieje realna dźwignia projektowa.

Studium wymusza trzy osobne wykresy: idealny potencjał obiegu, bilans urządzenia i bilans misji. Jedna kolumna „sprawność PGC” nie potrafi ich zastąpić.

Studium 2: czego uczy koncepcja Tarnawskiego

Monografia Piotra Tarnawskiego opisuje mechanicznie cykliczny silnik turbowałowy pracujący według przybliżenia Humphreya.⁷ Najcenniejszy nie jest pojedynczy deklarowany wynik, lecz ciąg odrzucanych wariantów: kolejne mechanizmy rozrządu zmieniają przeciwne momenty, straty przepływu, synchronizację, uszczelnienie i możliwość doprowadzenia mediów.

W jednym etapie urealnienie paliwa, chłodzenia, promieniowania i rekuperacji obniżyło wynik modelowy z około $0{,}399$ do $0{,}281$. Spadek o $0{,}118$ jest praktyczną demonstracją bilansu z poprzedniego studium: idealizacje nie tworzą małej korekty na końcu, lecz mogą zmienić ranking architektur.

Finalna koncepcja z wirującymi komorami uzyskała w modelu CFD około $36{,}6\%$, a w obliczeniu analitycznym około $37{,}5\%$. Różnica $0{,}9$ punktu procentowego świadczy o spójności rachunków w ich wspólnym zakresie, nie o walidacji eksperymentalnej. Obie metody dzielą część warunków i założeń, w tym idealizowaną szczelność i przyjęte sprawności mechaniczne.

Najbardziej ryzykowny węzeł — segmentowe uszczelnienie pracujące cieplnie i mechanicznie — autor proponuje zbadać osobno przed pełnym prototypem. To właściwa logika redukcji ryzyka: najpierw eksperyment rozstrzygający o wykonalności interfejsu, potem zintegrowana hamownia.

Studium prowadzi do czterech reguł:

1. zachowywać tabelę założeń przy każdym wariancie;
2. nie porównywać modeli z produkcyjnymi silnikami jak wyników tej samej dojrzałości;
3. rozdzielać zgodność modeli od walidacji pomiarem;
4. budować program prób wokół węzła, który może unieważnić całą korzyść.

Studium 3: RDE, turbina i lokalny zysk

NASA przeanalizowała stanowisko łączące RDE, przepływ obejściowy z eżektorem i turbinę.⁸ Model CFD komory połączono z algebraicznym modelem mieszania, aby uzupełnić skąpe pomiary stanowiskowe. Celem próby była przede wszystkim interakcja RDE–turbina, nie optymalizacja osiągów komory.

Analiza wykazała trzy różne prawdy jednocześnie:

1. samo zdarzenie detonacyjne zwiększało ciśnienie całkowite względem stanu przed reakcją;
2. przepływ na wylocie komory był w analizowanym układzie silnie uwarunkowany geometrią i pozostawał nietypowy dla odbiornika;
3. ograniczenie wlotu oraz długość urządzenia powodowały całkowitą stratę ciśnienia na granicy całego urządzenia.

Nie ma tu sprzeczności. Pierwszy wniosek dotyczy procesu, trzeci — większej powierzchni kontrolnej. Gdyby w publikacji pokazano tylko lokalny przyrost, czytelnik otrzymałby prawdziwy pomiar i fałszywą odpowiedź systemową.

Do wykazania korzyści wałowej potrzebne byłyby jeszcze: pełny bilans obejścia, zsynchronizowany stan przed turbiną, rzeczywista praca stopnia, porównanie z konwencjonalnym źródłem o tych samych wejściach oraz kara masy i chłodzenia. Studium jest więc wzorcem uczciwej negatywnej informacji: demonstrator może potwierdzić fizykę fali, a jednocześnie nie potwierdzić zysku urządzenia.

Dojrzałość mierzy się zamykaniem ryzyk

Krótka praca gorąca potwierdza zapłon i istnienie trybu, ale nie ustalony cieplnie stan, liczbę ponownych uruchomień ani trwałość materiału. Program NASA dotyczący ciągłego obiegu detonacyjnego wskazywał ekstremalne obciążenia cieplne i potrzebę aktywnego chłodzenia, a zarazem udokumentował próby dochodzące do stanu cieplnie ustalonego oraz wielokrotne uruchomienia.⁹ To inny poziom dowodu niż pojedynczy obraz fali, choć nadal nie jest to pełny dowód resursu.

Dojrzałość PGC można śledzić przez kolejne pytania:

1. Czy proces daje powtarzalny cykl i domyka masę oraz energię?
2. Czy poprawnie zdefiniowany zysk jest dodatni na granicy komory?
3. Czy pozostaje dodatni po wlocie, mieszaniu, turbinie lub dyszy?
4. Czy urządzenie utrzymuje wynik w zakresie warunków i przejść?
5. Czy ściana, uszczelnienia, łożyska i sensory osiągają wymagane życie?
6. Czy masa, chłodzenie, sterowanie i obsługa pozostawiają korzyść misji?

Negatywna odpowiedź nie zawsze kończy program. Wskazuje poziom, na którym powinien powstać następny eksperyment. Najgorszym krokiem jest przejście do większego demonstratora bez zamknięcia bilansu, bo skala mnoży równocześnie energię, koszt i ryzyko.

Jak czytać komunikat o przełomie PGC

Rzetelna ocena nie zaczyna się od procentu w nagłówku. Najpierw trzeba ustalić rzeczownik w mianowniku: sprawność czego, ciśnienie gdzie i praca na jakiej granicy?

Wynik jest wiarygodniejszy, gdy publikacja:

• pokazuje wspólny protokół odniesienia do Braytona;
• oddziela stan idealny, prognozę modelu i pomiar;
• raportuje cykl graniczny, rozrzut między cyklami i stan cieplny;
• zamyka masę, energię i niepewność;
• podaje zarówno średnie, jak i niejednorodność oraz widmo;
• obejmuje wlot i odbiornik w powierzchni kontrolnej;
• ujawnia koszt chłodzenia, przepłukania i urządzeń pomocniczych;
• opisuje liczbę prób, ponownych uruchomień i stan urządzenia po badaniu.

Czerwoną flagą jest lokalny pik ciśnienia nazwany pressure gain, średnia bez definicji, porównanie przy różnych temperaturach maksymalnych, CFD bez bilansu albo procent przeniesiony z obiegu idealnego na misję. Podobnie należy traktować słowo „ciągły” użyte wobec urządzenia, którego pole jest w istocie okresowe: ciągły dopływ nie oznacza stanu ustalonego.

Najkrótsza poprawna konkluzja brzmi: PGC ma realną podstawę termodynamiczną, lecz korzyść jest hipotezą systemową aż do zamknięcia bilansu na właściwej granicy i w reprezentatywnym czasie życia. Obieg Humphreya pokazuje, gdzie może powstać przewaga. Dopiero pomiar egzergii, pracy lub ciągu pokazuje, czy urządzenie jej nie zużyło.

Źródła

pominięto geometrie, mieszaniny, warunki inicjacji, częstotliwości, nastawy wtrysku, sekwencje zaworów, progi sterowania i procedury prób współczesnych urządzeń PGC. Wszystkie przykłady liczbowe poza jawnie oznaczonymi wynikami źródłowymi są neutralnymi demonstracjami rachunku, a nie danymi projektowymi.