Komory spalania, stateczność procesu i obciążenia termoakustyczne

Problem inżynierski, teza i zakres

Komora jest reaktorem przepływowym, rezonatorem, granicą ciśnieniową i gorącą strukturą. Zaburzenie spalania zmienia ciśnienie, prędkość, temperaturę i skład; pola te wracają przez akustykę, przepływ, zasilanie i konstrukcję do następnego wydzielania energii. Przedmiotem analizy jest ta zamknięta pętla.

Monografia prowadzi dwie niezależne osie oceny. Pierwsza pyta, czy energia małego zaburzenia rośnie, czy zanika. Druga pyta, czy nawet pole ograniczone, przejściowe albo wymuszone tworzy dopuszczalną historię temperatury, naprężenia i uszkodzenia. Ujemne tempo wzrostu odpowiada tylko na pierwsze pytanie: nie dowodzi dopuszczalnego obciążenia cieplno-mechanicznego. Obie osie spotykają się dopiero w argumentacji kwalifikacyjnej.^1,17

Teza tej monografii brzmi: znak lokalnego bilansu jest tylko jedną własnością określonego stanu i granicy układu spalanie–akustyka–przepływ–struktura. Dowód przydatności wymaga osobnego wykazania odpowiedzi na perturbacje, obciążeń cieplno-mechanicznych, niepewności oraz pokrycia konfiguracji. Pierwsza oś dowodu dotyczy wzrostu energii zaburzenia, druga — dopuszczalności obciążenia; żadna nie zastępuje drugiej.^1,6,8

Tekst prowadzi od skal, energii i odpowiedzi źródła przez nienormalność oraz nieliniowość do pomiarów, obciążeń i kwalifikacji. Formalizm zaawansowany można odłożyć, zachowując rozróżnienie modu odseparowanego i podprzestrzeni modów nakładających się. Nierozwiązane pola wirowe lub entropowe wymagają rozszerzenia stanu, gdy ich energia przekracza niepewność modelu.¹

Spalanie nierównomierne (ang. rough combustion) oznacza podwyższoną, zwykle szerokopasmową lub przerywaną aktywność, nie mechanizm. Odpowiedź przejściową od samowzbudzenia odróżnia historia energii i zachowanie po ustaniu wymuszenia. Granica obejmuje wszystkie dynamicznie istotne magazyny i interfejsy, także stanowisko, jeśli jego impedancja zmienia dynamikę.^9,10

Opis nie jest instrukcją budowy, modyfikacji ani badania silnika. do ręcznego uzupełnienia wyłącznie w kontrolowanej dokumentacji konkretnego programu: wykonawcze geometrie komory, wtryskiwacza, powierzchni materiału pędnego, przegród i absorberów; składy materiałów pędnych; parametry oraz nastawy robocze; własne częstotliwości i marginesy konkretnego systemu; progi kwalifikacyjne, zabezpieczenia i procedury celowego wzbudzania. Nieznanego lub uszkodzonego silnika ani jego fragmentu nie wolno przemieszczać ani poddawać próbom; miejsce przekazuje się służbom EOD i HAZMAT.

Dwa modele bazowe: komora ciekłopędna i stałopędna

Napęd ciekły. Zaburzenie ciśnienia przez drogę zasilania zmienia strumień masy, pęd strug i przygotowanie mieszaniny. Bezwładność, podatność, przemiana fazowa, mieszanie, transport i chemia nadają odpowiedzi amplitudę oraz fazę: zasilanie → wtrysk → opóźnione wydzielanie ciepła → pole akustyczne i hydrodynamiczne → zasilanie oraz płomień.¹ Wirowość i entropia mogą po czasie konwekcji zmienić spalanie lub wytworzyć dźwięk, a chłodzenie i dysza zmieniają straty, odbicie i podatność konstrukcji.

Napęd stały. Cofająca się powierzchnia jest źródłem masy, a przewodzenie, przemiany powierzchniowe i płomień tworzą pamięć cieplno-chemiczną. Ciśnienie oraz prędkość równoległa do powierzchni wyznaczają dwie drogi sprzężenia: transport ciepła i odpowiedź powierzchni → generacja gazu → przepływ oraz akustyka → transport przy powierzchni.⁷ Cząstki rozpraszają energię, a cofanie powierzchni i erozja dyszy zmieniają objętość, pole i odpływ.⁸

W napędzie ciekłym główna pamięć prowadzi przez hydraulikę, wtrysk i przygotowanie reagentów; w stałym — przez odpowiedź powierzchni i ewolucję objętości. Wspólna księgowość energetyczna wymaga mocy ze znakiem, czyli dodatniej dla dopływu do wybranej granicy i ujemnej dla odpływu; dalej skrót „moc” zachowuje tę konwencję. Przekazanie energii do konstrukcji jest transferem wewnętrznym, gdy konstrukcja należy do analizowanego układu; stratą jest dopiero $D_s$ albo odpływ poza granicę.

Stan średni, zaburzenia i hierarchia skal czasu

Analiza dynamiczna wymaga rozdzielenia stanu, wokół którego pracuje układ, od zmian niosących energię oscylacji. Dla dowolnej mierzalnej wielkości pola $f(\mathbf{x},t)$ przyjmujemy rozkład

$$f(\mathbf{x},t)=\overline{f}(\mathbf{x},t)+f'(\mathbf{x},t), \qquad \overline{f'}=0.$$

Symbol $\mathbf{x}$ oznacza położenie, $t$ — czas, kreska średnią w oknie długim wobec oscylacji, lecz krótkim wobec zmiany punktu pracy, a prim odchylenie. Wielkość, średnia i zaburzenie mają ten sam wymiar. Warunek $\overline{f'}=0$ nie zeruje energii, zależnej m.in. od $p'^2$ i $|\mathbf u'|^2$; poprawne średnie nie wykluczają dużego obciążenia dynamicznego.

Okno średniej musi być długie wobec oscylacji i krótkie wobec zmiany punktu pracy; inaczej miesza tło z zaburzeniem. Linearyzacja małych $f'/\overline f$ daje postać i tempo, lecz nie przewiduje cyklu granicznego, histerezy ani przełączania modów przy dużej amplitudzie.

Znaczenie opóźnień ujawnia bezwymiarowa relacja

$$\Pi_i=\omega\tau_i=2\pi f\tau_i.$$

Tutaj $\tau_i[\mathrm s]$, $f[\mathrm{Hz}]$, a $\omega=2\pi f[\mathrm{rad\,s^{-1}}]$; $\Pi_i$ jest bezwymiarowe. Dla $\Pi_i\ll1$ proces nadąża, przy wartości rzędu jedności silnie zmienia fazę, a dla $\Pi_i\gg1$ zwykle filtruje szybkie zmiany. Jest to porządkowanie pasm, nie próg kwalifikacyjny.

Przygotowanie i reakcja obejmują w cieczy rozpad, przemianę fazową, mieszanie i chemię, a w napędzie stałym — przewodzenie oraz przemiany powierzchni; model równowagowy nie wyznacza ich fazy. Transport rozdziela konwekcję $L/U$ od propagacji akustycznej $L/c$, a rzeczywiste mody zależą od impedancji.⁶ Zegary konstrukcji, dyfuzji ciepła, zasilania, chłodzenia, sterowania i stanowiska pozostają odrębne; zbieżność częstości jest tylko hipotezą do sprawdzenia polem i fazą.

Dla każdego mechanizmu źródło, zwłoka, odpływ, postać i pomiar wyznaczają tempo modu albo odpowiedź podprzestrzeni.

Bilans energii zaburzenia i definicja marginesu

Bilans mocy nadaje mechanizmom znak, jednostkę i granicę, dzięki czemu modele wyjaśniają zmianę energii, a nie tylko odtwarzają widmo.^1,6,8

Od równań zaburzeń do energii akustycznej

Najprostszy przypadek ustala sens energii i jednostek. Dla małych zaburzeń jednorodnego, nieruchomego gazu doskonałego równania pędu i ciśnienia z objętościowym źródłem ciepła mają postać

$$\overline{\rho}\,\frac{\partial \mathbf{u}'}{\partial t}+\nabla p'=\mathbf 0,$$

$$\frac{\partial p'}{\partial t} +\overline{\rho}\,\overline c^{\,2}\nabla\!\cdot\!\mathbf{u}' =(\overline\gamma-1)\dot q'.$$

Prim oznacza zaburzenie; $p'[\mathrm{Pa}]$, $\mathbf u'[\mathrm{m\,s^{-1}}]$, $\overline\rho[\mathrm{kg\,m^{-3}}]$, $\overline c[\mathrm{m\,s^{-1}}]$ i $\dot q'[\mathrm{W\,m^{-3}}]$ oznaczają odpowiednio ciśnienie, prędkość, średnią gęstość, prędkość dźwięku i szybkość wydzielania ciepła.

Po pomnożeniu równań przez $\mathbf u'$ i $p'/(\overline\rho\overline c^2)$, dodaniu oraz użyciu tożsamości iloczynu i $\overline\rho\overline c^2=\overline\gamma\overline p$ otrzymujemy

$$\frac{\partial e_a}{\partial t}+\nabla\!\cdot(p'\mathbf u') =\frac{\overline\gamma-1}{\overline\gamma\,\overline p}\,p'\dot q',$$

gdzie

$$e_a= \frac{p'^2}{2\overline\rho\,\overline c^{\,2}} +\frac{\overline\rho\,|\mathbf u'|^2}{2}$$

jest gęstością energii akustycznej: pierwszy składnik odpowiada ściśliwości, drugi ruchowi. Po scałkowaniu po $V$

$$E_a(t)=\int_V e_a\,\mathrm dV.$$

Ponieważ $\overline\rho\,\overline c^{\,2}$ ma wymiar paskala,

$$\left[\frac{p'^2}{\overline\rho\,\overline c^{\,2}}\right] =\mathrm{Pa}=\mathrm{J\,m^{-3}}, \qquad [\overline\rho|\mathbf u'|^2] =\mathrm{kg\,m^{-1}s^{-2}}=\mathrm{J\,m^{-3}}.$$

Całka ma jednostkę dżula. $p'\mathbf u'$ ma wymiar $\mathrm{W\,m^{-2}}$, a $p'\dot q'/\overline p$ — $\mathrm{W\,m^{-3}}$; po całkowaniu są to waty, zgodnie z $\mathrm dE_a/\mathrm dt$.

Miara dotyczy małych zaburzeń stanu quasi-stałego. Gorący przepływ niejednorodny dodaje wymianę ze stanem średnim, entropię i wirowość, a duża amplituda człony nieliniowe; $E_a$ pozostaje odniesieniem dla bilansu i wymiarów.

Intuicję daje fala stojąca w prostej objętości: antywęzły ciśnienia i prędkości leżą w innych miejscach, a źródło ciepła może pokrywać się z jednym z nich tylko częściowo. Ta sama amplituda źródła zasila albo osłabia pole zależnie od fazy i przestrzennego nakładania. Dlatego przed nadaniem znaków mocom trzeba określić postać pola, granicę magazynowanej energii i samą miarę energii; formalizm poniżej zapisuje właśnie tę księgowość.

Bilans mocy, granica i konwencja znaków

Niech $\mathbf n$ będzie normalną skierowaną na zewnątrz gazu, $S_o$ zbiorem otwartych interfejsów, a $S_s$ podatną granicą konstrukcji. Dla $\mathbf z_a=(p',u'_1,u'_2,u'_3)^{\mathsf T}$ wyprowadzona energia ma postać $e_a=\tfrac12\mathbf z_a^H\mathbf H_E\mathbf z_a$, gdzie

$$\mathbf H_E=\operatorname{diag}\!\left( \frac{1}{\overline\rho\,\overline c^{2}}, \overline\rho,\overline\rho,\overline\rho \right).$$

Ten pomost operatorowy jest szablonem dla klasy symetryzowalnej linearyzowanych układów, nie wyprowadzeniem uniwersalnym. Obowiązuje w oknie quasi-stałym, w którym stan bazowy, $\mathbf H_E$ i współczynniki operatora są zamrożone. Wymaga dodatniego hermitowskiego symetryzatora, regularności potrzebnej do całkowania przez części, zgodnych warunków granicznych oraz hermitowskich macierzy strumienia po symetryzacji. Wtedy można zapisać

$$\mathbf H_E\frac{\partial\mathbf x}{\partial t} +\sum_j\mathbf A_j\frac{\partial\mathbf x}{\partial x_j} +\mathbf B\mathbf x=\mathbf f_q, \qquad E_b=\frac12\int_{\mathcal D_b}\mathbf x^H\mathbf H_E\mathbf x\,\mathrm d\mathcal D,$$

gdzie $\mathcal D_g=V$, a $\mathcal D_c$ obejmuje także konstrukcję. Macierze $\mathbf A_j$ zawierają strumienie, w tym adwekcję przez $\overline{\mathbf U}$, a $\mathbf B$ człony bez pochodnych, w tym gradienty stanu średniego.

Aby wskazać podstawę tego wzorca bez odwołania do modelu szczególnego, pomnóżmy przyjęte równanie z lewej strony przez $\mathbf x^H$ i weźmy część rzeczywistą. Dla zamrożonej hermitowskiej $\mathbf H_E$ oraz hermitowskich $\mathbf A_j$ zachodzą tożsamości

$$\operatorname{Re}(\mathbf x^H\mathbf H_E\partial_t\mathbf x) =\frac12\partial_t(\mathbf x^H\mathbf H_E\mathbf x),$$

$$\operatorname{Re}(\mathbf x^H\mathbf A_j\partial_j\mathbf x) =\frac12\partial_j(\mathbf x^H\mathbf A_j\mathbf x) -\frac12\mathbf x^H(\partial_j\mathbf A_j)\mathbf x,$$

oraz $\operatorname{Re}(\mathbf x^H\mathbf B\mathbf x)=\tfrac12\mathbf x^H(\mathbf B+\mathbf B^H)\mathbf x$. Wynika stąd lokalna tożsamość

$$\begin{aligned} \partial_t\!\left(\frac12\mathbf x^H\mathbf H_E\mathbf x\right) &+\frac12\sum_j\partial_j(\mathbf x^H\mathbf A_j\mathbf x)\\ &=\operatorname{Re}(\mathbf x^H\mathbf f_q) +\frac12\mathbf x^H\!\left[ \sum_j\partial_j\mathbf A_j-\mathbf B-\mathbf B^H \right]\!\mathbf x . \end{aligned}$$

Po scałkowaniu twierdzenie Gaussa zamienia dywergencję na $\tfrac12\int_{\partial V}\mathbf x^H\mathbf A_n\mathbf x\,\mathrm dS$. Przeniesienie jej na prawą stronę daje człon graniczny ze znakiem minus, natomiast pochodne $\partial_j\mathbf A_j$ i hermitowska część $\mathbf B$ dają człon objętościowy. Wydzielenie z nich części zależnych od przepływu średniego prowadzi do

$$\begin{aligned} P_{\mathrm{adv}}^{\mathrm{ex}} &=-\frac12\int_{\partial V} \mathbf x^H\mathbf A_n^{(\overline U)}\mathbf x\,\mathrm dS,\\ P_{\mathrm{grad}}^{\mathrm{ex}} &=\frac12\operatorname{Re}\int_V \mathbf x^H\!\left[ \sum_j\frac{\partial\mathbf A_j^{(\overline U)}}{\partial x_j} -\mathbf B^{(\overline U)} -\mathbf B^{(\overline U)H} \right]\mathbf x\,\mathrm dV . \end{aligned}$$

Superskrypt „ex” oznacza człon wyprowadzony z równań, $\mathbf A_n^{(\overline U)}=\sum_j n_j\mathbf A_j^{(\overline U)}$, a $H$ transpozycję sprzężoną. Konkretne $P_{\mathrm{adv}}^{\mathrm{ex}}$ i $P_{\mathrm{grad}}^{\mathrm{ex}}$ trzeba wyprowadzić z wybranego zestawu zmiennych, symetryzatora i granic. Model AVPE może ilustrować znaczenie przepływu średniego w układzie rakietowym, lecz nie jest dowodem powyższej ogólnej tożsamości.⁶

Zakres wyprowadzenia kończy się na zamrożonym symetryzatorze i nieruchomej dziedzinie. Jeżeli $\mathbf H_E=\mathbf H_E(\boldsymbol\theta(t))$, to

$$\frac{\mathrm d}{\mathrm dt}\left(\frac12\mathbf x^H\mathbf H_E\mathbf x\right) =\operatorname{Re}(\mathbf x^H\mathbf H_E\dot{\mathbf x}) +\frac12\mathbf x^H\dot{\mathbf H}_E\mathbf x,$$

więc ostatni składnik musi wejść jawnie do bilansu. Dla ruchomej dziedziny $\mathcal D(t)$ twierdzenie transportowe Reynoldsa daje ponadto

$$\frac{\mathrm d}{\mathrm dt}\int_{\mathcal D(t)}e\,\mathrm dV =\int_{\mathcal D(t)}\partial_t e\,\mathrm dV +\int_{\partial\mathcal D(t)}e\,\mathbf v_b\!\cdot\!\mathbf n\,\mathrm dS,$$

gdzie $\mathbf v_b$ jest prędkością granicy. Człon z $\dot{\mathbf H}_E$, strumień transportowy ruchomej granicy oraz możliwa zależność pozostałych współczynników od czasu wymagają osobnego bilansu czasowo zmiennego; przedstawione wyprowadzenie nie dowodzi jego pełnej postaci i nie pozwala użyć stałego marginesu z jednego okna quasi-stałego.

Jeżeli stan nie zawiera entropii, wirowości lub innych pól hydrodynamicznych wpływających na wynik, ich wymianę zapisuje się jako model zamknięcia $P_h^{\mathrm{cl}}$ z własną niepewnością. Model musi określać wejścia, wyjście, znak i fazę, pasmo oraz parametry identyfikowane na danych innych niż końcowy margines. Musi też przewidzieć zmianę w co najmniej jednym niewykorzystanym stanie; strukturalna reszta skorelowana z jego wejściem falsyfikuje zamknięcie. Nie wolno dopasować $P_h^{\mathrm{cl}}$ wyłącznie tak, aby wyzerować resztę. Kanoniczny bilans gazowy brzmi

$$\boxed{ \frac{\mathrm dE_a}{\mathrm dt} =P_q+P_{\mathrm{adv}}^{\mathrm{ex}} +P_{\mathrm{grad}}^{\mathrm{ex}}+P_h^{\mathrm{cl}} -P_o-P_{\mathrm{abs}}-D_w-D_v-P_{a\to s} +R_{\mathrm{num},g} }$$

Znaki odnoszą się do tej samej granicy: dodatnie $P_q,P_{\mathrm{adv}}^{\mathrm{ex}},P_{\mathrm{grad}}^{\mathrm{ex}},P_h^{\mathrm{cl}}$ zwiększają energię gazu; dodatnie $P_o,P_{\mathrm{abs}},D_w,D_v,P_{a\to s}$ ją zmniejszają. Dyssypację objętościową rozkłada się tylko raz,

$$D_v=D_{v,\mathrm{rest}}+D_p,$$

gdzie $D_p\geq0$ jest częścią cząstkową, a $D_{v,\mathrm{rest}}\geq0$ pozostałą. Reszta numeryczna mierzy zachowanie po dyskretyzacji i nie zastępuje mechanizmu.

W najprostszym modelu wkład spalania wynosi

$$P_q(t)=\int_V \frac{\overline\gamma-1}{\overline\gamma\,\overline p} \,p'(\mathbf x,t)\dot q'(\mathbf x,t)\,\mathrm dV.$$

Jest to treść kryterium Rayleigha: dodatnia średnia zasila pole, ujemna je osłabia. Przeciwne fazy mogą się przestrzennie znosić, więc para sygnałów punktowych nie zastępuje całki. Dla okresu $T$ lub odpowiedniego okna

$$\langle P_q\rangle_T= \frac{1}{T}\int_{t_0}^{t_0+T}P_q(t)\,\mathrm dt.$$

Dla harmonicznych amplitud zespolonych średnia moc na otwartej granicy przyjmuje postać

$$\langle P_o\rangle =\frac{1}{2}\operatorname{Re}\!\int_{S_o} \widehat p\,\widehat{\mathbf u}^{\,*}\!\cdot\!\mathbf n\,\mathrm dS,$$

gdzie gwiazdka oznacza sprzężenie zespolone. Dysza może wynosić energię i odbijać falę, a zasilanie ją odbierać lub dostarczać; znak wymaga amplitudy i fazy obu pól.⁶

Jeśli $\mathbf v_s'$ jest prędkością podatnej ścianki, to

$$P_{a\to s}(t)=\int_{S_s}p'\mathbf v_s'\!\cdot\!\mathbf n\,\mathrm dS.$$

Dodatni wynik oznacza pracę pola nad konstrukcją; sprężysta ścianka może ją później oddać. Po dodaniu kinetycznej i sprężystej energii $E_s$ transfer gaz–ścianka kasuje się i powstaje kanoniczny bilans sprzężony:

$$\boxed{ \frac{\mathrm d(E_a+E_s)}{\mathrm dt} =P_q+P_{\mathrm{adv}}^{\mathrm{ex}} +P_{\mathrm{grad}}^{\mathrm{ex}}+P_h^{\mathrm{cl}} +P_{s,\mathrm{ext}} -P_o-P_{\mathrm{abs}}-P_{s,\mathrm{out}} -D_w-D_v-D_s+R_{\mathrm{num},c} }.$$

$P_{s,\mathrm{ext}}$ jest zewnętrznym dopływem do konstrukcji, a $P_{s,\mathrm{out}}$ odpływem przez mocowania lub inne interfejsy. Po włączeniu pola hydrodynamicznego lub absorbera do stanu jego dotychczasowy odpływ staje się transferem wewnętrznym, więc bilans wyprowadza się ponownie.

Konwencja modu i podprzestrzeni. Indeksu $k$ oraz wielkości $E_{b,k},P_{*,k},m_{b,k}$ używa się tylko dla modu energetycznie odseparowanego: ortogonalnego względem $\mathbf W_E$ albo objętego kontrolowanym przybliżeniem, w którym wszystkie wyrazy krzyżowe są mniejsze od niepewności decyzji. W przeciwnym razie obiektem rachunku jest podprzestrzeń $K$. Dla $\mathbf x=\boldsymbol\Phi_K\mathbf a$

$$E_{b,K}=\frac12\mathbf a^H \boldsymbol\Phi_K^H\mathbf W_E\boldsymbol\Phi_K\mathbf a, \qquad N_{b,K}(\mathbf a)=\left\langle\frac{\mathrm dE_{b,K}}{\mathrm dt}\right\rangle-R_{\mathrm{num},b},$$

z pełnymi wyrazami krzyżowymi. Pojedynczy skalar opisuje całą $K$ tylko wtedy, gdy forma mocy ze znakiem jest jednostajnie określona względem energii; konserwatywny odpowiednik to

$$m_{b,K}^{\min}=\inf_{\mathbf a\ne0} \frac{-N_{b,K}(\mathbf a)}{2E_{b,K}(\mathbf a)}.$$

Jeżeli infimum nie jest dodatnie albo wynik zależy od kierunku $\mathbf a$, raportuje się ograniczoną odpowiedź lub pełną formę kwadratową, nie „margines modu”. Ta konwencja obowiązuje w całej monografii; późniejsze dopiski „mod odseparowany” są tylko odsyłaczem.

Poniższy słownik zbiera symbole używane dalej. Znaki mocy wynikają z dwóch bilansów powyżej.

Symbol	Znaczenie	Granica / warunek użycia	Jednostka
$b\in\{g,c\}$, $g$, $c$	wybór bilansu; gaz; gaz z konstrukcją	jedna zadeklarowana granica w całym liczniku i energii	—
$k$, $K$	mod odseparowany; podprzestrzeń modów nakładających się	$k$ tylko przy ortogonalności energetycznej lub pomijalnych wyrazach krzyżowych	—
$E_a,E_s,E_b$	energia akustyczna, konstrukcyjna i energia dla granicy $b$	$E_g=E_a$, $E_c=E_a+E_s$; wszystkie jawne magazyny	$\mathrm J$
$N_g,N_c,N_{b,K}$	fizyczny licznik mocy gazu, układu sprzężonego i podprzestrzeni	bez reszt numerycznych i bez dopasowywania brakującej fizyki	$\mathrm W$
$P_q,P_o,P_{\mathrm{abs}}$	moc spalania; odpływ przez otwarcie; odpływ do absorbera	$P_q$ ze znakiem; dwa odpływy nie nakładają się	$\mathrm W$
$P_h^{\mathrm{cl}},P_{a\to s}$	zamknięcie pól hydrodynamicznych; praca gazu nad konstrukcją	pierwsze tylko dla pól poza stanem; drugie znika w bilansie $c$	$\mathrm W$
$P_{\mathrm{adv}}^{\mathrm{ex}},P_{\mathrm{grad}}^{\mathrm{ex}}$	adwekcja i wymiana z gradientami stanu bazowego	wyłącznie po wyprowadzeniu z wybranego modelu symetryzowalnego	$\mathrm W$
$D_w,D_v,D_s$	dyssypacja przyścienna, objętościowa i konstrukcyjna	średnio nieujemne; $D_v=D_{v,\mathrm{rest}}+D_p$	$\mathrm W$
$R_{\mathrm{num}},R_{\mathrm{exp}}$	reszta zachowania po dyskretyzacji; reszta pomiarowego bilansu	test jakości, nigdy człon fizyczny poprawiający licznik	$\mathrm W$
$\mathbf H_E,\mathbf W_E$	ciągły symetryzator i dyskretna macierz tej samej energii	dodatnie na dopuszczalnym stanie; $\mathbf W_E$ zawiera kwadraturę	zależna od skalowania
$\mathbf M_d,\mathbf M_s$	macierz uogólnionego zagadnienia; macierz masy konstrukcji	$\mathbf M_d\not\equiv\mathbf W_E$ bez dowodu; $\mathbf M_s$ dla stopni swobody struktury	zależna od dyskretyzacji
$\sigma,m$	tempo amplitudy i przeciwny mu margines zaniku	skalar dla odseparowanego $k$; dla $K$ forma lub dolne ograniczenie	$\mathrm{s^{-1}}$
$\Omega$	zbiór stanów i konfiguracji objętych twierdzeniem	tylko obszar z jawną mapą danych, modeli i niepewności	—

Bilans pomiarowy definiuje się jawnie jako

$$R_{\mathrm{exp}}=\frac{\mathrm dE_{\mathrm{obs}}}{\mathrm dt}-\sum_i P_{i,\mathrm{est}}, \qquad [R_{\mathrm{exp}}]=\mathrm W.$$

$E_{\mathrm{obs}}$ jest energią odtworzoną z sygnałów pomiarowych, a $P_{i,\mathrm{est}}$ estymatami wszystkich mocy dla tej granicy. Raport podaje kowariancję składników i $R_{\mathrm{exp}}$; domknięcie wymaga zgodności reszty z zerem w niepewności oraz braku struktury względem stanu i wejść modeli. $R_{\mathrm{num}}$ sprawdza zachowanie bilansu po dyskretyzacji, natomiast $R_{\mathrm{exp}}$ — domknięcie obserwowanego układu.

Przed analizą definiuje się skalę $S_P$, na przykład większą z sumy modułów dominujących mocy i $2E/T_{\mathrm{ref}}$, oraz dopuszczony udział $\varepsilon_{\mathrm{cl}}$. Domknięcie wymaga jednocześnie

$$|R_\ell|\leq U_{R_\ell}, \qquad \frac{U_{R_\ell}}{S_P}\leq\varepsilon_{\mathrm{cl}}, \qquad \ell\in\{\mathrm{num},\mathrm{exp}\},$$

Pierwszy warunek sprawdza zgodność reszty z zerem, drugi zdolność wykrycia braku wpływającego na wynik. Program ustala $S_P$ i $\varepsilon_{\mathrm{cl}}$; reszta zdolna odwrócić znak licznika blokuje wniosek.^13,14

Od przyrostu energii do liniowego tempa wzrostu

W krótkim odcinku quasi-stałym amplitudę energetycznie odseparowanego modu zapisujemy dla granicy $b\in\{g,c\}$:

$$A_{b,k}(t)=A_{b,k,0}\exp(\sigma_{b,k}t).$$

Dla granicy gazowej $E_{g,k}\equiv E_{a,k}$, a dla sprzężonej $E_{c,k}\equiv E_{a,k}+E_{s,k}$. Zatem

$$E_{b,k}(t)=E_{b,k,0}\exp(2\sigma_{b,k}t), \qquad \sigma_{b,k}=\frac{1}{2E_{b,k}} \frac{\mathrm dE_{b,k}}{\mathrm dt}.$$

Współczynnik ma jednostkę $\mathrm{s^{-1}}$. Jego znak można estymować z nachylenia $\ln|A_{b,k}|$ albo z właściwego licznika mocy podzielonego przez $2E_{b,k}$. Użycie $E_{a,k}$ z licznikiem sprzężonym albo $E_{a,k}+E_{s,k}$ z członem $-P_{a\to s,k}$ miesza granice i daje pozorny wynik.

Gdy warunek odseparowania nie zachodzi, obowiązuje wcześniejszy bilans podprzestrzeni $K$, z pełną formą $\boldsymbol\Phi_K^H\mathbf W_E\boldsymbol\Phi_K$; tę samą $\mathbf W_E$ stosuje się później w $G_W(t)$.

Chwilowa pochodna energii może zawierać odwracalną wymianę; $\sigma_{b,k}$ opisuje obwiednię po uśrednieniu okresu. W przejściu raportuje się historie energii, stanu i mocy zamiast jednej stałej.

Margines lokalny i margines wykazywany

Dla odseparowanego modu i stanu $\boldsymbol\theta$ niech $N_{g,k}$ i $N_{c,k}$ oznaczają uśrednione prawe strony dwóch kanonicznych bilansów po usunięciu $R_{\mathrm{num}}$. W pierwszym pozostaje $-P_{a\to s,k}$; w drugim transfer ten kasuje się. Margines dla $b\in\{g,c\}$ ma jedną postać

$$M_{P,b,k}=-N_{b,k}, \qquad m_{b,k}=\frac{M_{P,b,k}}{2\langle E_{b,k}\rangle_k} =-\sigma_{b,k}.$$

Dodatni znak oznacza zanik, a reszty jedynie poszerzają niepewność. Dla nakładających się modów stosuje się $m_{b,K}^{\min}$ lub pełną odpowiedź podprzestrzeni.

Niech $\mathcal J(\boldsymbol\theta)$ zawiera odseparowane mody $k$ i konieczne podprzestrzenie $K$, a $\mu_{c,j,L}$ oznacza odpowiednio dolną granicę $m_{c,k}$ albo $m_{c,K}^{\min}$. Dla zbioru stanów i konfiguracji $\Omega$

$$m_{\mathrm{wyk}}= \inf_{\boldsymbol\theta\in\Omega} \inf_{j\in\mathcal J(\boldsymbol\theta)} \mu_{c,j,L}(\boldsymbol\theta).$$

Zapis nie ustanawia progu. $\boldsymbol\theta$ obejmuje stan cieplny, rozrzut, interfejsy, zużycie i otoczenie; wynik dotyczy tylko zbadanej części $\Omega$ i stanów dopuszczających linearyzację. Przejścia, wzrost niemodalny i skończona perturbacja wymagają osobnych kryteriów.

Skończony zbiór prób wyznacza dolną granicę tylko przez jawną mapę pokrycia. W poniższym przykładzie współrzędne $\boldsymbol\theta$ są bezwymiarowe i znormalizowane do dopuszczonych zakresów, a $\rho=0{,}177$ jest promieniem w normie euklidesowej tej przestrzeni. Dla punktowej dolnej granicy $m_{c,k,L}=1{,}9 \mathrm{s^{-1}}$ i zwalidowanego ograniczenia Lipschitza $L_\theta=1{,}2 \mathrm{s^{-1}}$ obwiednia wynosi

$$\underline m_C=m_{c,k,L}-L_\theta\rho =1{,}9-1{,}2\times0{,}177 =1{,}69 \mathrm{s^{-1}}.$$

Inna normalizacja zmienia $L_\theta$, dlatego skalowanie i stałą waliduje się łącznie. Wynik obejmuje tylko tę komórkę; emulator musi dodatkowo odjąć zwalidowany błąd.^13,14,15

Nadrzędna macierz dowodowa

Macierz porządkuje późniejsze wyniki; nie zastępuje lokalnego wniosku fizycznego.

Kategoria dowodu	Minimalna podstawa	Warunek decyzji	Wynik przy braku
stan, granica i konfiguracja	stan średni, okno, normalne, interfejsy, wersje	wszystkie człony dotyczą tego samego układu	luka zakresu
energia, postać i norma	pola, właściwa energia, pełna forma dla $K$	postać rozróżniona; wyrazy krzyżowe ujęte lub ograniczone	nieznana tożsamość pola
źródło i przepływ	sygnał pomiarowy związany z $\dot q'$, dokładne człony przepływu i zwalidowane $P_h^{\mathrm{cl}}$	znak, faza i całki zachowane w paśmie	znak licznika nierozstrzygnięty
granice, absorber i konstrukcja	fazowe pola lub impedancje dla $P_o,P_{\mathrm{abs}},P_{a\to s},P_{s,\mathrm{ext}},P_{s,\mathrm{out}}$	każdy transfer sklasyfikowany raz; jawne warunki zerowe	bilans niedomknięty
dyssypacja	$D_w,D_{v,\mathrm{rest}},D_p,D_s$ z modelem stanu	$D_v=D_{v,\mathrm{rest}}+D_p$, bez podwójnego liczenia	zapas nierozstrzygnięty
tempo, przejście i obciążenie	obwiednia lub energia podprzestrzeni, $G_W$, historia cieplno-mechaniczna	zgodność niezależnej estymaty i osobne kryteria dla przejścia oraz obciążenia	brak decyzji o danym mechanizmie
domknięcie i pokrycie	reszty, niepewność i mapa $\Omega$	reszty akceptowalne; wniosek mieści się w danych	niewykazany zakres

pominięto lokalnie: liczbowe marginesy, progi decyzji oraz parametry celowej perturbacji konkretnego programu, potrzebne do zastosowania bilansu w jego kontrolowanej dokumentacji.

Pole akustyczne, rodziny modów i impedancje graniczne

Bilans wymaga przestrzennych pól ciśnienia i prędkości. Ich dopuszczalne postacie wybierają geometria, ośrodek i warunki brzegowe, a źródła oraz straty nadają im tempo. Mod jest więc rozwiązaniem przestrzenno-czasowym określonego modelu, nie nazwą piku. Poniżej model zimnego, jednorodnego gazu służy jako odniesienie; impedancje, przepływ, stan gorący i model toru stopniowo przywracają fizykę rzeczywistego układu.

Od równań zaburzeń do równania falowego

W ogólnym modelu odniesienia przyjmijmy małe zaburzenia w jednorodnym, nieruchomym, bezstratnym gazie o średniej gęstości $\overline\rho$ i prędkości dźwięku $\overline c$. Po pominięciu niestacjonarnego źródła ciepła liniowe równania użyte w bilansie energii mają postać

$$\overline\rho\frac{\partial \mathbf u'}{\partial t}+\nabla p'=\mathbf 0, \qquad \frac{\partial p'}{\partial t} +\overline\rho\,\overline c^{\,2}\nabla\!\cdot\!\mathbf u'=0.$$

Obliczenie rozbieżności pierwszego równania i pochodnej czasowej drugiego, a następnie eliminacja $\nabla\!\cdot\!\mathbf u'$, prowadzą do równania falowego

$$\boxed{ \frac{\partial^2 p'}{\partial t^2} -\overline c^{\,2}\nabla^2p'=0 }.$$

Oba wyrazy mają jednostkę $\mathrm{Pa\,s^{-2}}$. Częstości dopuszczalne powstają dopiero po określeniu dziedziny i wszystkich warunków brzegowych.

Dla składowej harmonicznej stosujemy zapis

$$p'(\mathbf x,t)=\operatorname{Re}\!\left\{ \widehat p(\mathbf x)\exp(\mathrm i\omega t) \right\},$$

gdzie $\widehat p$ jest zespoloną amplitudą ciśnienia, $\mathrm i^2=-1$, a $\omega$ jest częstością kołową w $\mathrm{rad\,s^{-1}}$. Radian jest bezwymiarowy. Po podstawieniu otrzymujemy równanie Helmholtza

$$\boxed{ \nabla^2\widehat p+k^2\widehat p=0, \qquad k=\frac{\omega}{\overline c} }.$$

Liczba falowa $k[\mathrm{m^{-1}}]$ daje rozwiązanie niezerowe tylko dla warunków zgodnych z granicami. W modelu zachowawczym wartości własne są rzeczywiste; odpływ, straty lub aktywne spalanie zwykle czynią je zespolonymi. Część oscylacyjną i tempo interpretuje się razem z przyjętą konwencją czasową.

Model odniesienia wyjaśnia rodziny modów. Tam, gdzie przepływ i jego gradienty zmieniają propagację oraz odpływ przez dyszę, potrzebny jest operator linearyzowany wokół przepływu średniego.^1,6

Warunek brzegowy jako relacja ciśnienia i prędkości

W modelu idealnie sztywnej, nieprzepuszczalnej ściany normalna składowa prędkości znika:

$$\widehat u_n=\widehat{\mathbf u}\!\cdot\!\mathbf n=0.$$

Z harmonicznego równania pędu

$$\mathrm i\omega\overline\rho\,\widehat{\mathbf u} =-\nabla\widehat p$$

wynika wtedy warunek Neumanna $\partial\widehat p/\partial n=0$. Jest on przybliżeniem granicy o bardzo dużej impedancji, a nie uniwersalną własnością metalu, wtryskiwacza czy dyszy. Bardziej ogólny lokalny warunek zapisujemy przez zespoloną impedancję akustyczną

$$Z_b(\omega)=\frac{\widehat p}{\widehat u_n}, \qquad \frac{\partial\widehat p}{\partial n} =-\mathrm i\omega\overline\rho\, \frac{\widehat p}{Z_b}.$$

Przy normalnej skierowanej na zewnątrz impedancja ma jednostkę

$$[Z_b]=\frac{\mathrm{Pa}}{\mathrm{m\,s^{-1}}} =\mathrm{Pa\,s\,m^{-1}} =\mathrm{kg\,m^{-2}s^{-1}}.$$

Admitancja $Y_b=1/Z_b\,[\mathrm{m\,s^{-1}Pa^{-1}}]$ opisuje ruch wywołany ciśnieniem. Mała wartość odpowiada ścianie prawie sztywnej; duża może oznaczać podatną powierzchnię lub interfejs do innej objętości, niekoniecznie otwarcie geometryczne.

Składowa reaktywna opisuje odwracalne magazynowanie energii i zmienia fazę oraz częstość; rezystancyjna odpowiada za średni transfer lub dyssypację. Moc wychodząca przez $S_b$ wynosi

$$\left\langle P_b\right\rangle =\frac12\operatorname{Re} \int_{S_b}\widehat p\,\widehat u_n^{,*}\,\mathrm dS,$$

Znak mocy zależy od fazy i orientacji normalnej, nie od samego $|Z_b|$. Granica sprzężona z zasilaniem lub spalaniem może dostarczać energię; jej klasyfikację daje bilans całej pętli.

Jeżeli ciśnienie w jednym miejscu wywołuje przepływ lub ruch w innym, granicę reprezentuje operator przestrzenny albo macierz admitancji. Redukcja wielokanałowego interfejsu do jednej liczby wymaga zachowania postaci modu i bilansu mocy.

Jak geometria tworzy rodziny modów

Nazwy „osiowy”, „styczny” i „promieniowy” wskazują kierunek zmian pola. Ich pochodzenie pokazuje idealizowany cylinder, choć rzeczywista postać obejmuje całą objętość.

Niech $x$ oznacza współrzędną wzdłuż osi, $r$ odległość od osi, a $\theta$ kąt azymutalny. W jednorodnym cylindrycznym modelu szukamy rozwiązań rozdzielonych w postaci

$$\widehat p(r,\theta,x)=R(r)\,\Theta(\theta)\,X(x).$$

Po rozdzieleniu zmiennych całkowita liczba falowa spełnia

$$k^2=k_x^2+k_\perp^2.$$

Relacja pokazuje współistnienie zmian osiowych i poprzecznych oraz wpływ warunku końcowego na każdą postać.

Rodzina osiowa. Przy sztywnej ścianie bocznej i najniższej postaci poprzecznej zmienność skupia się w $X(x)$:

$$\frac{\mathrm d^2X}{\mathrm dx^2}+k_x^2X=0.$$

Sztywne końce dają zerowy gradient, a kolejne rozwiązania dodają węzły. Rzeczywiste impedancje przesuwają je wraz z częstością i stratami; reguła połówek fali jest tylko przybliżeniem.

Rodzina styczna, czyli azymutalna. Okresowość $\Theta(\theta+2\pi)=\Theta(\theta)$ dopuszcza całkowite $m$ i pary $\cos(m\theta),\sin(m\theta)$. Dla $m>0$ tworzą sektory przeciwnych faz; w osiowej symetrii zdegenerowana para składa się w falę stojącą lub biegnącą.

Niejednorodne granice wybierają orientację lub rozszczepiają parę. Dwa bliskie tony rozróżnia dopiero rozłożona faza.

Rodzina promieniowa. $R(r)$ spełnia równanie Bessela. Regularność usuwa gałąź osobliwą, a sztywna ściana zeruje pochodną radialną; kolejne liczby falowe tworzą współosiowe powierzchnie węzłowe.

„Promieniowy” oznacza zmianę w $r$: dla $m=0$ między pierścieniami, dla $m>0$ razem ze zmianą azymutalną. W dziedzinie pierścieniowej obie gałęzie mogą spełniać warunki na dwóch granicach.

Ogólna postać łączy trzy kierunki, które mieszają zwężenie, dysza, niejednorodny gaz i impedancje. Nazwa wskazuje składową dominującą, lecz identyfikacja opiera się na całym polu.¹

Co naprawdę obciąża granice komory

Wtryskiwacz lub powierzchnia spalania łączy granicę mechaniczną z aktywnym źródłem; w napędzie stałym trzeba oddzielić bierną obudowę od cofającej się powierzchni.^1,8 Ściana przewodzi i rozprasza ciepło, a jej impedancja może sprzęgać gaz z konstrukcją. Dysza wynosi, odbija i przekształca energię zależnie od przepływu.⁶ Objętość boczna magazynuje albo odbiera energię przez admitancję zależną od częstości. Każda z tych granic może przesunąć postać i węzły pola.

Trzy różne znaczenia częstości własnej

Wynik	Co obejmuje	Właściwe użycie i granica
zimny model bez przepływu	znaną geometrię, nieruchomy ośrodek i idealizowane impedancje	topologia, symetria i degeneracja; bez konwekcji, gradientów, rzeczywistego odpływu i $\dot q'$
model z przepływem	pola $\overline{\mathbf U},\overline\rho,\overline p,\overline c$, transport i gradienty; odpływ może dać zespolone wartości własne	postać oraz tempo właściwe temu polu; nie istnieje uniwersalna korekta modelu zimnego 6
stan gorący	zmienny skład i temperaturę, rzeczywiste impedancje, straty i zespoloną odpowiedź spalania	wartość własna sprzężenia spalanie–przepływ–akustyka, ewentualnie konstrukcja

Jedna temperatura średnia koryguje tylko $c^2=\gamma p/\rho$, pomijając odbicie, konwekcję, podatność i źródło. Wyniki między stanami łączy ciągłość postaci, nie sama częstość.⁸

Poziom zaawansowany — operator niesamosprzężony i wzrost niemodalny

Operator z przepływem, stratami, otwartymi granicami i aktywnym źródłem jest zwykle niesamosprzężony, a w normie energii często nienormalny. Wartości własne opisują zachowanie asymptotyczne, lecz nie ograniczają koniecznie krótkotrwałego wzrostu: nieortogonalne składowe mogą się rozfazować mimo ujemnych temp modalnych. Jest to dynamika liniowa, nie dowód nieliniowego atraktora.²⁰

Po dyskretyzacji zapiszmy uogólnione zagadnienie

$$\mathbf M_d\dot{\mathbf x}=\mathbf L\mathbf x, \qquad \mathbf A=\mathbf M_d^{-1}\mathbf L.$$

Prawe mody spełniają $\mathbf L\boldsymbol\phi_j=\lambda_j\mathbf M_d\boldsymbol\phi_j$, a lewe $\boldsymbol\psi_i^H\mathbf L=\lambda_i\boldsymbol\psi_i^H\mathbf M_d$. Przy prostych wartościach własnych można wybrać biortogonalną normalizację

$$\boldsymbol\psi_i^H\mathbf M_d\boldsymbol\phi_j=\delta_{ij}.$$

$\mathbf M_d$ nie jest automatycznie wagą energii $\mathbf W_E$. Lewy mod określa czułość współczynnika na wymuszenie, prawy przestrzenną odpowiedź, a duża norma pary sygnalizuje słabe uwarunkowanie projekcji.²⁰ Biortogonalność względem $\mathbf M_d$ nie uchyla wcześniejszej reguły podprzestrzeni $K$.

Maksymalne liniowe wzmocnienie energii w czasie $t$ mierzy propagator:

$$\boxed{ G_W(t)= \max_{\mathbf x_0\ne0} \frac{\|\exp(\mathbf A t)\mathbf x_0\|_{W_E}^2} {\|\mathbf x_0\|_{W_E}^2} = \left\| \mathbf W_E^{1/2}\exp(\mathbf A t)\mathbf W_E^{-1/2} \right\|_2^2 }, \qquad \|\mathbf x\|_{W_E}^2=\mathbf x^H\mathbf W_E\mathbf x=2E_b.$$

Warunek $\operatorname{Re}\lambda_j<0$ zapewnia zanik asymptotyczny, ale nie wymusza $G_W(t)\leq1$. Dowód lokalny obejmuje więc tempa modalne i wzrost przejściowy w tej samej normie; osobne badanie nieliniowe rozstrzyga wyjście poza obszar linearyzacji.²⁰

Miniaturowy model bez wymiaru fizycznego pokazuje ten efekt rachunkowo. Niech $\mathbf W_E=\mathbf I$ i

$$\dot{\mathbf x}=\begin{bmatrix}-1&4\\0&-2\end{bmatrix}\mathbf x, \qquad \exp(\mathbf At)= \begin{bmatrix} e^{-t}&4(e^{-t}-e^{-2t})\\0&e^{-2t} \end{bmatrix}.$$

Wartości własne $-1$ i $-2$ gwarantują zanik asymptotyczny. Ich prawe wektory $(1,0)^T$ oraz $(-4,1)^T$ nie są ortogonalne. Dla $\mathbf x_0=(0,1)^T$ te składowe początkowo znoszą się w pierwszej współrzędnej, lecz zanikają z różnymi tempami. W chwili $t=\ln2$ otrzymujemy $\mathbf x=(1,0{,}25)^T$, a zatem w normie euklidesowej $\|\mathbf x\|_I^2=1{,}0625>\|\mathbf x_0\|_I^2=1$; stąd $G_I(\ln2)\geq1{,}0625$. Jest to przejściowy wzrost jednego małego warunku początkowego, nie cykl graniczny ani twierdzenie o odpowiedzi na dowolną perturbację skończoną. Liczby nie reprezentują parametrów komory.²⁰

Jak rozpoznać postać modu w danych rozłożonych

Pojedynczy kanał może ujawnić ton i obwiednię, lecz przy węźle przeoczy mod, a rezonans przewodu może utworzyć pozorne maksimum. Postać wymaga rozłożonych, zsynchronizowanych pomiarów z zespolonymi transmitancjami torów; pełny model pomiaru podaje rozdział o instrumentacji.⁴

Hipoteza postaci	Minimalna informacja rozróżniająca	Typowa luka
osiowa	próbki w więcej niż jednym położeniu osiowym, wzór amplitudy i przejść fazy	jeden przekrój może widzieć ten sam lokalny znak kilku postaci
azymutalna	wiele kierunków obwodu; przyrost fazy odróżnia falę biegnącą od stojącej	symetria rozmieszczenia może aliasować liczby azymutalne
radialna	więcej niż jeden promień albo zwalidowana rekonstrukcja pola	jedna ściana nie rozróżnia postaci podobnych na brzegu, innych wewnątrz
sprzężona	wspólny czas dla ciśnienia, drgań/odkształceń i sygnału spalania	zgodna częstość bez fazy nie określa postaci ani kierunku mocy

Zespolony wektor danych porównuje się z próbkowanymi modami, zachowując względną amplitudę, fazę, wspólną skalę i niepewność toru. Dla par zdegenerowanych dopasowuje się podprzestrzeń i śledzi jej współczynniki, aby odróżnić obrót fali od przełączania. Rekord identyfikacji podaje stan, operator, impedancje, wartość własną, postać, transmitancje, miarę zgodności i obszary nieobserwowalne.

pominięto lokalnie: wykonawczą geometrię granic, mapy impedancji, własne postacie i częstości oraz dokładne rozmieszczenie torów konkretnego napędu.

Dynamiczna odpowiedź spalania w napędzie ciekłym

Aktywne sprzężenie powstaje, gdy pole zmienia późniejsze wydzielanie ciepła. Zasilanie, wtrysk, przygotowanie faz, transport i reakcja nadają odpowiedzi amplitudę oraz fazę; o mocy rozstrzyga ich przestrzenne nakładanie z modem.^1,3

Łańcuch od zaburzenia do wydzielania ciepła

Wejścia $p'_c,p'_f,\dot m'$ i $\mathbf u'$ nie są zamienne: różnica ciśnień pobudza hydraulikę, strumień jest odpowiedzią drogi, a prędkość może zmieniać warstwę ścinającą bez dużej zmiany całkowitego dopływu.

Etap	Pamięć i skutek dynamiczny
zasilanie	bezwładność, ściśliwość i podatność magazynują energię; opory ją rozpraszają, a odgałęzienia dodają mody; impedancja wiąże fazę $\dot m'$ z $p'_c$ 5
wtrysk	zmienny pęd i spadek ciśnienia kształtują strugi oraz warstwy ścinające; interfejs jest granicą impedancyjną i przetwornikiem, więc jego moc oddziela się od późniejszego $P_q$
przygotowanie i mieszanie	rozpad, nagrzewanie i odparowanie tworzą rozkład czasów; po zaniku klasycznej granicy faz zastępują je mieszanie gęstych strug i własności rzeczywiste; różne drogi konwekcji mogą wnosić przeciwne fazy
reakcja i płomień	skład, temperatura, powierzchnia oraz położenie reakcji tworzą przestrzenne $\dot q'$; o mocy decyduje jego korelacja z $p'$, także gdy sama chemia jest szybka

Pełny łańcuch można zatem zapisać symbolicznie jako

$$\begin{aligned} (p'_c,p'_f,\dot m',\mathbf u') &\longrightarrow \text{zasilanie i różnica ciśnień}\\ &\longrightarrow \text{niestacjonarny wtrysk}\\ &\longrightarrow \text{przygotowanie faz i mieszanie}\\ &\longrightarrow \text{transport oraz reakcja}\\ &\longrightarrow \dot q'(\mathbf x,t) \longrightarrow p',\mathbf u'. \end{aligned}$$

Ostatnia strzałka zamyka pętlę; równoległe oddziaływanie prędkości, zasilania i entropii pozostaje w modelu, jeśli wpływa na fazę lub moc.

Funkcja przenoszenia jako model lokalnej odpowiedzi

Aby oddzielić stan średni od odpowiedzi dynamicznej, wybierzmy w małym obszarze lub w kontrolnej części strefy spalania bezwymiarowe wejście

$$\chi'(t)=\frac{s'(t)}{\overline s},$$

gdzie $s$ może oznaczać jednoznacznie zdefiniowany strumień masy, prędkość, ciśnienie albo inną wielkość przyczynowo poprzedzającą lokalne spalanie. Wielkości $s'$ i $\overline s$ muszą mieć te same jednostki, zatem $\chi'$ jest bezwymiarowe. Wyjściem jest względne zaburzenie objętościowej szybkości wydzielania ciepła

$$y_q'(\mathbf x,t)=\frac{\dot q'(\mathbf x,t)}{\overline{\dot q}(\mathbf x)}.$$

W dziedzinie częstotliwości lokalną funkcję przenoszenia definiujemy jako model liniowy

$$\boxed{ \hat y_q(\mathbf x,\omega) =F_{\chi q}(\mathbf x,\omega) \hat\chi(\mathbf x,\omega) }, \qquad F_{\chi q}(\mathbf x,\omega) =\frac{G(\mathbf x)}{1+\mathrm i\omega\tau_f(\mathbf x)} \exp[-\mathrm i\omega\tau_d(\mathbf x)].$$

Bezwymiarowe $G$ jest małosygnałowym wzmocnieniem, a $\tau_f,\tau_d[\mathrm s]$ czasem filtracji i zwłoką. Iloczyny $\omega\tau$ i $F_{\chi q}$ są bezwymiarowe, więc $\hat{\dot q}'=\overline{\dot q}F_{\chi q}\hat\chi$ ma jednostkę $\mathrm{W\,m^{-3}}$.

Przy konwencji harmonicznej $\exp(\mathrm i\omega t)$ czyste opóźnienie ma czynnik $\exp(-\mathrm i\omega\tau_d)$. Dla dodatniego rzeczywistego $G$ faza modelu wynosi

$$\arg F_{\chi q} =-\arctan(\omega\tau_f)-\omega\tau_d.$$

Ujemna faza oznacza opóźnienie. Dla wejścia innego niż ciśnienie dochodzi faza drogi od $\hat p$, a zwłokę odtwarza się z ciągłości fazy modulo $2\pi$.

Kilka porównywalnych etapów wymaga sumy lub iloczynu członów, a kilka wejść — modelu wielowejściowego:

$$\begin{aligned} \hat y_q={}&F_{pq}\frac{\hat p}{\overline p} +F_{\dot m q}\frac{\hat{\dot m}}{\overline{\dot m}}\\ &+\mathbf F_{uq}\!\cdot \frac{\hat{\mathbf u}}{U_{\mathrm{ref}}} +\hat r, \end{aligned}$$

Jawne wejścia i $\hat r$ są bezwymiarowe. Ich korelacja ogranicza identyfikowalność, a zależność współczynników od amplitudy zawęża model do zbadanego cyklu i stanu.

Gdy $\overline{\dot q}(\mathbf x)$ jest bliskie zeru, stosuje się wymiarowe $\hat{\dot q}'$ albo dodatnią skalę całkową, aby uniknąć sztucznie nieograniczonego wzmocnienia.

Od odpowiedzi lokalnej do mocy modu odseparowanego

Niech zespolona amplituda ciśnienia rozpatrywanego modu $k$ będzie zapisana jako

$$\hat p_k(\mathbf x,\omega_k) =A_k\,\phi_{p,k}(\mathbf x),$$

gdzie $A_k$ ma jednostkę paskala, a postać $\phi_{p,k}$ jest bezwymiarowa przy jawnie wybranej normalizacji. Lokalne bezwymiarowe wejście odpowiedzi spalania zapisujemy w postaci

$$\hat\chi_k(\mathbf x,\omega_k) =B_k(\mathbf x,\omega_k)A_k,$$

więc $B_k[\mathrm{Pa^{-1}}]$ opisuje drogę od modalnego ciśnienia przez granicę, zasilanie i wtrysk do wejścia funkcji $F_{\chi q}$. Rozdzielenie tych funkcji lokalizuje niepewność fazy.

W obszarze, w którym dozwolona jest normalizacja lokalna,

$$\begin{aligned} \hat{\dot q}'_k(\mathbf x)={}& \overline{\dot q}(\mathbf x) F_{\chi q}(\mathbf x,\omega_k)\\ &\times B_k(\mathbf x,\omega_k)A_k. \end{aligned}$$

Po wstawieniu do harmonicznej postaci członu Rayleigha z części „Bilans energii zaburzenia i definicja marginesu” otrzymujemy

$$\begin{aligned} \langle P_{q,k}\rangle &=\frac12\operatorname{Re}\!\int_V \frac{\overline\gamma-1} {\overline\gamma\,\overline p} \hat p_k\,\hat{\dot q}'^{\,*}_k\,\mathrm dV\\ &=\frac{|A_k|^2}{2} \operatorname{Re}\!\int_V \frac{\overline\gamma-1} {\overline\gamma\,\overline p} \phi_{p,k}\, \overline{\dot q}\, F_{\chi q}^{*}\, B_k^{*}\,\mathrm dV. \end{aligned}$$

Jednostki całki redukują się do wata. Moc wymaga $F_{\chi q}$, wejścia $B_kA_k$ i nakładania z $\phi_{p,k}$; obszary przeciwnych faz mogą się znosić.

Jak dane ograniczają niepewność modelu

Model zasilania ogranicza drogę ciśnienie–przepływ, obserwacje wtrysku — transport, a zsynchronizowana próba ze spalaniem — zespoloną odpowiedź. Optyka bez kalibracji pozostaje wskaźnikiem.^3,4

Rozważmy syntetyczne mini-studium bez parametrów silnika. Po korekcji dwa tory o różnej dynamice i trzeci kanał odtwarzają tę samą postać, a dwa przekroje ciśnienia zachowują jej ciągłość; ogranicza to hipotezy artefaktu i przełączenia. Cross-PSD wskaźnika ciepła z ciśnieniem pokazuje zgodność fazową z możliwym zasilaniem rozpoznanej postaci, a niezależna obwiednia $E_b$ określa znak pełnego tempa. Sama para sygnałów nie wyznacza jednak $\langle P_{q,k}\rangle$: potrzebne są przestrzenna rekonstrukcja $\dot q'$, kwadratura całki modalnej, określone pasmo i normalizacja oraz propagacja kalibracji. Pomiar nie rozdziela też $B_k$ od $F_{\chi q}$, bo współdzielą skorelowane wejście; wymaga to dodatkowego stanu, niezależnej zmiany fazy wejść lub sygnału wtrysku. Wniosek ogranicza się zatem do zgodności danych z hipotezą zasilania postaci przez źródło; zarówno wartość mocy, jak i przypisanie pierwotnej przyczyny dopływowi albo wtryskowi pozostają nierozstrzygnięte.^1,3,4

pominięto lokalnie: wykonawcze charakterystyki wtrysku i zasilania, liczbowe funkcje odpowiedzi oraz dane celowej perturbacji potrzebne do strojenia i kwalifikacji konkretnego napędu.

Dynamiczna odpowiedź powierzchni i spalania w napędzie stałym

W napędzie stałym pole oddziałuje na rozłożoną, cofającą się granicę, która przewodzi ciepło, ulega przemianom i zasila komorę gazem. Ciśnienie zmienia procesy powierzchniowe, a prędkość styczna warstwę przyścienną; statyczna zależność średniej szybkości spalania od ciśnienia nie opisuje obu dróg.^7,8

Pamięć cieplna, odpowiedź zespolona i moc źródła

Lokalne $p'$ i $u_t'$ zmieniają transport ciepła w warstwie przyściennej. Zaburzenie temperatury wnika na głębokość zależną od częstości, po czym zmienia przemiany powierzchniowe; pamięć tej drogi rozsuwa maksima ciśnienia, prędkości i strumienia produktów.⁷

W małym otoczeniu stanu średniego lokalną odpowiedź powierzchni zapisujemy wymiarowo:

$$\widehat{\dot m}_s'(\boldsymbol\xi,\omega) =R_p(\boldsymbol\xi,\omega)\widehat p(\boldsymbol\xi,\omega) +R_u(\boldsymbol\xi,\omega)\widehat u_t(\boldsymbol\xi,\omega),$$

gdzie $\boldsymbol\xi$ leży na powierzchni, a $\dot m_s'[\mathrm{kg\,m^{-2}s^{-1}}]$. Jednostki $R_p$ i $R_u$ wynikają z wejść. Obejmują one czułość cieplno-chemiczną i transportową, są zespolone, zależne od stanu, pasma i amplitudy, nie są zaś stałymi materiałowymi.⁷

Transport produktów i rozkład czasów reakcji tworzą dopiero objętościowe pole ciepła. Nazwijmy ten krok operatorem przyczynowym $\mathcal T_{sq}$:

$$\widehat{\dot q}'(\mathbf x,\omega) =\mathcal T_{sq}(\omega;\boldsymbol\theta)\widehat{\dot m}_s' =\int_{S_s}K_{sq}(\mathbf x,\boldsymbol\xi,\omega;\boldsymbol\theta) \widehat{\dot m}_s'(\boldsymbol\xi,\omega)\,\mathrm dS_{\xi}.$$

Jądro $K_{sq}$ zapisuje transport, mieszanie i reakcję; pojedyncze opóźnienie wymaga walidacji modułu i fazy.⁷

Cały równorzędny łańcuch ma więc postać

$$(p',u_t') \longrightarrow q_{\mathrm{conv},s}''{}' \longrightarrow (T_s',\mathcal C_s') \longrightarrow \dot m_s' \overset{\tau_{\mathrm{konw}},\tau_{\mathrm{chem}}}{\longrightarrow} \dot q'(\mathbf x,t) \longrightarrow \langle P_{q,k}\rangle .$$

$\mathcal C_s'$ oznacza zbiorcze zaburzenie stanu chemicznego powierzchni, nie skład ani recepturę. Dla modu odseparowanego ostatni etap wyznacza modalna całka mocy:

$$\left\langle P_{q,k}\right\rangle =\frac12\operatorname{Re} \int_V \frac{\overline\gamma-1} {\overline\gamma\,\overline p} \widehat p_k\, \widehat{\dot q}_k'^{\,*}\,\mathrm dV .$$

Próbka ogranicza odpowiedź tylko dla zachowanych temperatury, pasma, amplitudy i stanu powierzchni; nie dostarcza pola potrzebnego do $P_{q,k}$. Węzeł ciśnienia może współistnieć z dużym $u_t'$, a współliniowe wejścia nie rozdzielają $R_p$ i $R_u$.^7,8

Cząstki, dysza i ewolucja wolnej objętości

Cząstki rozpraszają energię przez poślizg i wymianę ciepła; dla modu odseparowanego $D_{p,k}$ zależy od relaksacji, położenia i częstości. Gdy uczestniczą w reakcji, ich wkład wchodzi też do opóźnionego $\dot q'$. Dysza tworzy zależny od przepływu odpływ i odbicie. Cofanie powierzchni zmienia $V_g(t)$, pole i granice, więc postać śledzi się po ciągłości. Dla odseparowanego modu i powolnej ewolucji wystarcza ciąg stanów quasi-ustalonych:^7,8

$$\begin{aligned} \mathcal L[\boldsymbol\theta(t)]\,\phi_k(t) &=\lambda_k(t)\,\phi_k(t),\\ \sigma_{c,k}(t) &=\frac{ P_{q,k}+P_{\mathrm{adv},k}^{\mathrm{ex}} +P_{\mathrm{grad},k}^{\mathrm{ex}}+P_{h,k}^{\mathrm{cl}} +P_{s,\mathrm{ext},k}-P_{o,k}-P_{\mathrm{abs},k} -P_{s,\mathrm{out},k}-D_{w,k}-D_{v,\mathrm{rest},k} -D_{p,k}-D_{s,k}} {2(E_{a,k}+E_{s,k})}. \end{aligned}$$

$P_{s,\mathrm{ext},k}$ wolno wyzerować tylko bez zewnętrznej mocy konstrukcyjnej. Gdy geometria zmienia się w oknie identyfikacji albo zbliżają się wartości własne, przechodzi się do operatora czasowo zmiennego i bilansu $K$. Dane historyczne potwierdzają możliwość czasowego zbliżenia mechanizmu przepływowego i modu podłużnego, nie jego parametry dla innego silnika.¹⁸

Mechanizmy tłumienia i warunki brzegowe

Tłumienie jest sumą mocy dla jednej postaci i granicy, nie stałą komory. Element stabilizujący może zmieniać źródło lub pole; NASA SP-8113 dlatego rozdziela przegrody ograniczające organizację sprzężenia od absorberów zwiększających odpływ.²

Księgowość mocy na granicach

Równanie mocy i znaki zdefiniowano w nadrzędnym bilansie. Dysza jest odpływem względem komory, lecz interfejsem wewnętrznym po włączeniu strugi do stanu. Straty gorące zależą od stanu i pasma; próba zimna może potwierdzić topologię albo impedancję, nie ich wartość podczas pracy.^6,8

Przegrody, absorbery i odporność marginesu

Mechanizmy różnią się miejscem w księgowości i selektywnością. Poniższa synteza nie podaje geometrii ani strojenia urządzeń.^2,6,7,8

Mechanizm	Człon bilansu	Charakter pasmowy lub modalny	Możliwy skutek uboczny	Dowód po zmianie konfiguracji
przegroda	zmiana $P_q$, postaci, $D_w$ i strat przepływowych	selektywna względem postaci; może rozszczepić lub przeorientować mod	lokalne ogrzewanie i obciążenie przegrody, przeniesienie energii do innej postaci	ponowne wartości własne, pola fazowe, bilans mocy i lokalna odpowiedź cieplno-konstrukcyjna 2
absorber lub wnęka celowa	$P_{\mathrm{abs}}$, a po rozszerzeniu granicy energia magazynu i jego dyssypacja	zwykle ograniczone pasmo; zależy od nakładania z modem	przesunięcie częstości, nowe maksimum pola lub obciążenie elementu	zmierzona albo zwalidowana admitancja z fazą, bilans i próba odpowiedzi dla objętych stanów 2
dysza i inne otwarte granice	$P_o$, $P_{\mathrm{adv}}^{\mathrm{ex}}$, $P_{\mathrm{grad}}^{\mathrm{ex}}$	zależne od kierunku fali, przepływu i postaci przy interfejsie	odbicie albo konwersja składowych; obciążenie przesuwa się do strugi lub mocowania	fazowe $p',u_n'$, model przepływu średniego i domknięcie strumienia 6
straty ścienne	$D_w=D_\nu+D_\kappa$	rosną lub maleją ze stanem warstwy przyściennej, częstością i amplitudą	zmienione pole cieplne, lokalne gradienty i trwałość ściany	walidacja gorącej impedancji/warstwy oraz transferu do temperatury i naprężenia 1,2
straty objętościowe i cząstkowe	$D_v=D_{v,\mathrm{rest}}+D_p$	zależne od relaksacji, rozkładu fazy i położenia względem modu	zmiana źródła lub transportu; ryzyko podwójnego policzenia wymiany	rozdzielone składniki, zgodność fazowa i bilans dla aktualnego stanu 7,8
dyssypacja konstrukcyjna	$D_s$ oraz $P_{s,\mathrm{out}}$	modalna; największa przy sprzężeniu postaci gazu i konstrukcji	większe drganie, zmęczenie albo przeniesienie obciążenia do mocowań	sprzężone postacie, kierunek $P_{a\to s}$, dyssypacja i historia naprężeń 2,17

Dla odseparowanego modu na granicy gazowej $m_{\mathrm{abs},g,k}=P_{\mathrm{abs},k}/(2E_{a,k})$, lecz po dołączeniu magazynu absorbera do stanu $P_{\mathrm{abs}}$ staje się transferem wewnętrznym. Większy odpływ może zarazem zmienić częstość, postać, $P_q$, pole cieplne lub obciążenie konstrukcji. Po modyfikacji powtarza się dowody dotkniętych mechanizmów, a na końcu aktualizuje ich kategorie w macierzy.^2,10

Nieliniowość, rozruch, wyłączenie i przełączanie modów

Liniowy wynik lokalny nie rozstrzyga odpowiedzi na skończony bodziec ani odległego cyklu granicznego. Praca spontaniczna, odpowiedź po perturbacji, amplituda ograniczona i przejścia są odrębnymi twierdzeniami.^1,4,11

Równanie amplitudy i pięć różnych twierdzeń

Minimalny model jakościowy jednego modu można zapisać jako

$$\begin{aligned} \frac{\mathrm dA}{\mathrm dt} &=g(A)A+F(t),\\ g(A)&=\sigma_{c,0}+\beta|A|^2-\gamma|A|^4,\\ E_c&=C_E|A|^2, \qquad C_E>0 . \end{aligned}$$

$\sigma_{c,0}$ opisuje tempo liniowe, $\beta,\gamma$ zależność od amplitudy, a $F(t)$ wymuszenie. Ujemne $\sigma_{c,0}$ może współistnieć ze stabilnym rozwiązaniem niezerowym, gdy nieliniowość tworzy drugi obszar przyciągania.

Dla przypadku autonomicznego $F=0$, rzeczywistych współczynników i $r=|A|$ faza pozostaje stała, a amplituda spełnia

$$\dot r=r(\sigma_{c,0}+\beta r^2-\gamma r^4), \qquad \dot y=2y(\sigma_{c,0}+\beta y-\gamma y^2), \quad y=r^2.$$

Niezerowe stany równowagi wymagają dodatniego pierwiastka równania $\sigma_{c,0}+\beta y-\gamma y^2=0$:

$$y_\pm=\frac{\beta\pm\sqrt{\beta^2+4\gamma\sigma_{c,0}}}{2\gamma},$$

o ile wyróżnik jest nieujemny, $\gamma\ne0$ i $y_\pm>0$. Dla $\gamma>0$ zero jest stabilne przy $\sigma_{c,0}<0$, a stabilność $r_*>0$ określa znak $2y_*(\beta-2\gamma y_*)$.

Przypadek przy $\gamma>0$	Gałęzie autonomiczne	Wniosek
nadkrytyczny, $\beta<0$	po przejściu przez $\sigma_{c,0}=0$ pojawia się mała stabilna gałąź, $y\simeq-\sigma_{c,0}/\beta$ dla $\sigma_{c,0}>0$	zero traci stabilność bez lokalnego obszaru bistabilnego
podkrytyczny, $\beta>0$	dla $-\beta^2/(4\gamma)<\sigma_{c,0}<0$ istnieją dodatnie $y_-$ niestabilne i $y_+$ stabilne, obok stabilnego zera	skończony bodziec może przekroczyć granicę basenu; punkt zwrotny jest przy $\sigma_{c,0}=-\beta^2/(4\gamma)$
wymuszony, $F(t)\ne0$	amplituda zależy także od widma i fazy wymuszenia	niezerowa odpowiedź podczas działania $F$ nie dowodzi samowzbudnego cyklu

Po odjęciu znanego wymuszenia rodzina trajektorii pozwala dopasować $\dot y/(2y)=\sigma_{c,0}+\beta y-\gamma y^2$: odcinek małosygnałowy ogranicza $\sigma_{c,0}$, a krzywizna wielu przebiegów — $\beta,\gamma$. Dryf stanu, postaci lub toru może udawać zależność od amplitudy, dlatego potrzebna jest walidacja na innych warunkach początkowych.¹

Przełączanie modów wymaga modelu co najmniej dwóch amplitud:

$$\begin{aligned} \dot A_1&=\sigma_{c,1}A_1 -\mu_{11}|A_1|^2A_1 -\mu_{12}|A_2|^2A_1,\\ \dot A_2&=\sigma_{c,2}A_2 -\mu_{22}|A_2|^2A_2 -\mu_{21}|A_1|^2A_2 . \end{aligned}$$

Wzajemne człony oznaczają zmianę dostępnego źródła; zanik jednego tonu przy pojawieniu się drugiego wymaga śledzenia obu postaci, energii i obciążenia sumarycznego.

Przejścia i historia stanu

Podczas rozruchu zmieniają się przepływy, temperatura, impedancje i płomień. Wzrost z wymuszeniem $F(t)$ może przejść w gałąź samopodtrzymującą, a jedno $\sigma_{c,k}$ z całego przebiegu miesza tło z zaburzeniem. J-2X oceniał przejścia osobno.^4,11

Wyłączenie nie jest odwróconym rozruchem: strumienie zanikają w różnym tempie, struktura oddaje ciepło, a opróżniane objętości zmieniają impedancje. Krótka odpowiedź wymuszona nadal może tworzyć niedopuszczalne obciążenie.

Histerezę rozpoznaje się po zależności od historii przy tym samym stanie obserwowanym. Najpierw ogranicza się dryf temperatury, zmianę konfiguracji i toru; zapis historii musi czynić twierdzenie falsyfikowalnym.

pominięto lokalnie: progi przejść, energia i sposób wprowadzania perturbacji oraz nastawy ochronne konkretnego programu. Bez tych danych artykuł nie jest procedurą próby.

Od pola termoakustycznego do obciążenia i uszkodzenia

To druga z zapowiedzianych osi: pole może być dynamicznie ograniczone, a mimo to wytworzyć niedopuszczalną historię lokalną. Łańcuch biegnie od $p'$, $\mathbf u'$ i $\dot q'$, przez trakcję oraz strumień ciepła, do temperatury, odkształcenia, naprężenia i miary uszkodzenia.^1,17

Po stronie gorącej rozdzielamy stan średni i oscylację:

$$\begin{aligned} q_{g\to w}''&=\overline q_{g\to w}''+q_w''{}', &q_w''{}'&=q_{\mathrm{conv}}''{}'+q_{\mathrm{chem}}''{}'+q_{\mathrm{rad}}''{}',\\ T&=\overline T+T', &\boldsymbol\sigma&=\overline{\boldsymbol\sigma}+\boldsymbol\sigma'. \end{aligned}$$

Konwekcję moduluje prędkość, turbulencja i warstwa przyścienna; człon chemiczny obejmuje zmiany składu, reakcji i rekombinacji; promieniowanie zależy od temperatury i składu ośrodka. Ich udział jest właściwy stanowi. Pomiary lokalne potwierdzają niejednorodność ogrzewania, nie dostarczają jednak mapy innej komory.^19,21

W jednorodnym przewodniku harmoniczna składowa wnika na skalę

$$\boxed{\delta_T=\sqrt{\frac{2\alpha}{\omega}}}, \qquad \alpha=\frac{k_m}{\rho_m c_m}.$$

Amplituda maleje w przybliżeniu jak $\exp(-z/\delta_T)$, a faza opóźnia się o $z/\delta_T$. Wzrost częstości zmniejsza $\delta_T$: szybkie wymuszenie lokalizuje zmianę temperatury i jej gradient przy powierzchni, zamiast równomiernie ogrzewać ściankę. Nie wynika z tego automatycznie lokalizacja $\boldsymbol\sigma'$. Położenie i znak maksimum naprężenia zależą od więzów, geometrii, temperaturowo zależnych własności oraz pełnego rozwiązania konstrukcyjnego; koncentracja naprężenia przy powierzchni jest tylko hipotezą wymagającą obliczenia i walidacji. Dla grubości porównywalnej z $\delta_T$ odpowiedź cieplna zależy już od obu granic i nie jest półprzestrzenią.^17,21

Drogi mają różne pasma bez uniwersalnej częstości granicznej. Konwekcję filtruje warstwa przyścienna, reakcję czasy chemiczne i transport, promieniowanie stan źródła, a przewodzenie dyfuzja w ściance. Chłodziwo dodaje bezwładność cieplną i dynamikę przepływu; kolejność tych skal ustala amplitudę i fazę lokalnego cyklu.^17,21

Po stronie chłodziwa przyjmujemy dwie kierunkowe normalne fizycznego transportu: $\mathbf n_g$ biegnie z gazu do ścianki, a $\mathbf n_c$ ze ścianki do chłodziwa. Dodatni strumień biegnie zgodnie z odpowiednią normalną, zatem

$$-k_m\nabla T\!\cdot\!\mathbf n_g=q_{g\to w}'', \qquad -k_m\nabla T\!\cdot\!\mathbf n_c =h_c(t)\,[T_w-T_c(t)].$$

Oscylacje $h_c'$ i $T_c'$ zmieniają amplitudę oraz fazę $T'$, więc stała temperatura chłodziwa wymaga uzasadnienia.²¹ Stan średni steruje własnościami, pełzaniem i relaksacją, a szybka składowa wyznacza cykl i może propagować wadę. Ich faza określa lokalną pętlę naprężenie–odkształcenie oraz ratcheting.^17,22

Operator przejścia od modu do odpowiedzi konstrukcji

Dla zespołu modów rekonstrukcja pola może mieć postać

$$\begin{aligned} p'(\mathbf x,t) &=\operatorname{Re} \sum_k A_k(t)\phi_{p,k}(\mathbf x) \exp(\mathrm i\omega_k t),\\ \mathbf u'(\mathbf x,t) &=\operatorname{Re} \sum_k A_k(t)\boldsymbol\phi_{u,k}(\mathbf x) \exp(\mathrm i\omega_k t),\\ \dot q'(\mathbf x,t) &=\operatorname{Re} \sum_k A_k(t)\phi_{q,k}(\mathbf x) \exp(\mathrm i\omega_k t). \end{aligned}$$

Postacie muszą mieć spójną normalizację i fazę. Globalna energia $E_{a,k}+E_{s,k}$ dotyczy tu składnika odseparowanego; przy nakładaniu obciążenie wyznacza pełna superpozycja $K$. Ani energia, ani poziom widma nie wskazują miejsca największego obciążenia. Antywęzeł ciśnienia jest kandydatem na dużą trakcję $\mathbf t_p=-p'\mathbf n$, lecz maksimum temperatury może leżeć gdzie indziej. Zmienny przepływ i spalanie można powiązać ze strumieniem za pomocą zwalidowanego operatora

$$\widehat q_w''(\mathbf x,\omega) =H_{qp}(\mathbf x,\omega)\widehat p +\mathbf H_{qu}(\mathbf x,\omega) \!\cdot\widehat{\mathbf u} +H_{qq}(\mathbf x,\omega) \widehat{\dot q}' .$$

Strumień $q_w''$ ma jednostkę $\mathrm{W\,m^{-2}}$, więc każda transmitancja ma jednostkę wynikającą z wejścia i zachowuje fazę. Taki operator redukuje wielowejściowy mechanizm tylko w paśmie oraz stanie objętym walidacją.^2,19,21

Identyfikacja łączy zsynchronizowane wejścia $p',\mathbf u',\dot q'$ z odpowiedziami $q_w'',T',\boldsymbol\varepsilon',\boldsymbol\sigma'$ w wielu miejscach. FRF lub cross-PSD zachowują moduł, fazę, koherencję i dynamikę torów. Skorelowane wejścia rozdziela model wielowejściowy wsparty dodatkowym stanem, niezależnym pobudzeniem albo sygnałem pomiarowym zmieniającym współliniowość; kowariancję identyfikacji propaguje się do temperatury, naprężenia i zapasu.^13,14,19

Operator sprawdza się na stanach niewykorzystanych do identyfikacji, kontrolując przyczynowość, pasmo, fazę, chłodzenie i mapowanie pole–siatka zachowujące strumień, siłę lub pracę. Prognoza obciążenia wymaga przewidzenia amplitudy, fazy i lokalizacji maksimum z niepewnością mniejszą niż dostępny zapas; samo odtworzenie punktów uczących jest jedynie wskaźnikiem.^13,14,21

Odpowiedź cieplna wynika z równania przewodzenia z warunkami na obu stronach ścianki; odpowiedź konstrukcyjna — z równań ruchu z obciążeniem ciśnieniowym i odkształceniem cieplnym. W zapisie symbolicznym:

$$\rho_m c_m\frac{\partial T}{\partial t} =\nabla\!\cdot(k_m\nabla T)+s_T,$$

$$\mathbf M_s\ddot{\mathbf d} +\mathbf C\dot{\mathbf d} +\mathbf K\mathbf d=\mathbf f(t),$$

$$\mathbf f(t)=\mathbf f_p(t)+\mathbf f_T[T(t)] .$$

Dywergencja $\nabla\!\cdot(k_m\nabla T)$ opisuje lokalną nierównowagę strumienia przewodzonego, natomiast $s_T[\mathrm{W\,m^{-3}}]$ jest rzeczywistym objętościowym źródłem ciepła, na przykład wydzielaniem wewnętrznym, i nie wolno utożsamiać tych wielkości. W rozważanym przypadku bez wewnętrznej generacji przyjmujemy $s_T=0$. Dla składowej średniej przewodzenie ustala powolne pole $\overline T$, a dla harmonicznej daje zanik i przesunięcie fazy opisane przez $\delta_T$. $\mathbf f$ jest sumą obciążenia ciśnieniowego i cieplnego.

Sprzężenie termoelastyczne obejmuje rozszerzenie cieplne, zależność sztywności i tłumienia od temperatury oraz pracę ciśnienia nad poruszającą się ścianką. Łańcuch jednokierunkowy pole → temperatura → naprężenie traci ważność, gdy odkształcenie cieplne przesuwa postać lub impedancję, ruch ścianki wnosi istotny $P_{a\to s,k}$, zmienia się chłodzenie albo zbliżają się postacie gazu i konstrukcji. Wtedy rozwiązanie dwukierunkowe musi zachować energię i fazę.^2,17

Reprezentatywnym mechanizmem jest zmęczenie termomechaniczne. Sama amplituda naprężenia nie wyznacza trwałości: znaczenie mają naprężenie średnie i temperatura średnia, kolejność cykli, lokalna plastyczność i relaksacja oraz rozmiar i położenie wady początkowej. Historia $\boldsymbol\sigma(t),\boldsymbol\varepsilon(t),T(t)$ tworzy pętlę naprężenie–odkształcenie; jej miarą może być zakres odkształcenia albo energia rozproszona w cyklu, zależnie od skalibrowanego mechanizmu. Dwie historie o tej samej amplitudzie naprężenia dadzą inną pętlę, gdy w jednej maksimum obciążenia mechanicznego przypada na gorącą, zmiękczoną część cyklu, a w drugiej na chłodniejszą; podwyższona temperatura średnia może dodatkowo zwiększyć relaksację i odkształcenie nieodwracalne.^17,22

Model życia jest ważny tylko dla mechanizmu i zakresu temperatur objętych kalibracją. Pełzanie wymaga czasu ekspozycji i zmiennych opisujących odkształcenie zależne od czasu, ratcheting — ewolucji odkształcenia średniego, a propagacja pęknięcia — stanu wady i miary napędu szczeliny; nie wolno zastępować ich jedną krzywą amplitudową. Proces, NDE i starzenie określają stan początkowy, zaś pomiar musi obejmować przewidywane maksimum albo walidować przeniesienie do niego.^2,10,17,22

Jawny przykład liczbowy z kontrolą jednostek

Rozważmy bezwymiarowo wyskalowany, syntetyczny układ niepowiązany z uzbrojeniem. Skale dobrano wyłącznie tak, aby wkłady ciśnieniowy i cieplny były porównywalne; nie reprezentują materiału ani stanu roboczego. Dla $|A|=60 \mathrm{kPa}$, $|\phi_p|=0{,}40$ i $|H_{qp}|=0{,}30 \mathrm{W\,m^{-2}Pa^{-1}}$ otrzymujemy lokalnie $|\widehat p|=24 \mathrm{kPa}$ oraz $|\widehat q_w''|=7{,}2 \mathrm{kW\,m^{-2}}$.

Dla półprzestrzeni o $\alpha=1{,}0\times10^{-5} \mathrm{m^2s^{-1}}$, $k_m=20 \mathrm{W\,m^{-1}K^{-1}}$ i $f=500 \mathrm{Hz}$ otrzymujemy $\delta_T=\sqrt{2\alpha/(2\pi f)}=7{,}98\times10^{-5} \mathrm m$.

Przy harmonicznym strumieniu rozwiązanie ma postać

$$\widehat T(z)= \frac{\widehat q_w''\delta_T}{k_m(1+\mathrm i)} \exp\!\left[-(1+\mathrm i)\frac{z}{\delta_T}\right].$$

Temperatura powierzchni opóźnia się o $45^\circ$, a $|\widehat T_s|=|\widehat q_w''|\delta_T/(\sqrt2k_m)=0{,}0203 \mathrm K$. Mała amplituda nie oznacza małego gradientu: zmiana skupia się na skali $\delta_T$.

W ustalonym punkcie $\mathbf x_*$ i dla stałego wektora jednostkowego $\mathbf n_*$ definiujemy podpisaną, liniową składową tensora

$$\sigma_{nn}(t)=\mathbf n_*^{\mathsf T} \boldsymbol\sigma(\mathbf x_*,t)\mathbf n_*, \qquad \widehat\sigma_{nn} =\mathbf n_*^{\mathsf T}\widehat{\boldsymbol\sigma}(\mathbf x_*)\mathbf n_*.$$

Dodatni znak oznacza rozciąganie w kierunku $\mathbf n_*$, a ujemny ściskanie. Niech $H_{\sigma_{nn}p}=40 \mathrm{Pa/Pa}$ odwzorowuje ciśnienie na tę samą składową w tym samym punkcie, a lokalny operator cieplny wynosi

$$\widehat\sigma_{nn,T}=H_{\sigma_{nn}T}\widehat T_s, \qquad |H_{\sigma_{nn}T}|=34{,}5 \mathrm{MPa\,K^{-1}},$$

co nie jest stałą materiałową. Zespolona składowa ciśnieniowa $\widehat\sigma_{nn,p}$ ma amplitudę $0{,}96 \mathrm{MPa}$, a zespolona składowa cieplna $\widehat\sigma_{nn,T}$ — $0{,}700 \mathrm{MPa}$. Ponieważ są to podpisane składowe liniowe tego samego tensora, punktu i kierunku, wolno najpierw utworzyć $\widehat\sigma_{nn}=\widehat\sigma_{nn,p}+\widehat\sigma_{nn,T}$. Dwa warianty fazy pokazują graniczne przypadki tej sumy:

$$\begin{aligned} \Delta\varphi_{\sigma_{nn}}=0^\circ: &&|\widehat\sigma_{nn}| &=0{,}96+0{,}70=1{,}66 \mathrm{MPa},\\ \Delta\varphi_{\sigma_{nn}}=180^\circ: &&|\widehat\sigma_{nn}| &=|0{,}96-0{,}70|=0{,}26 \mathrm{MPa}. \end{aligned}$$

Bez fazy cieplno-mechanicznej ciśnienie nie wyznacza amplitudy cyklu. W drugim, nieobserwowanym miejscu $|\widehat p|=15 \mathrm{kPa}$, lecz $|H_{qp}|=0{,}75 \mathrm{W\,m^{-2}Pa^{-1}}$, więc strumień wynosi $11{,}25 \mathrm{kW\,m^{-2}}$: mniejsze ciśnienie nie wyklucza większego lokalnego obciążenia cieplnego. W pełnym przypadku najpierw sumuje się tensorowe historie $\boldsymbol\sigma_p(t)+\boldsymbol\sigma_T(t)$, a dopiero z historii wypadkowej oblicza nieliniową lub zredukowaną miarę zmęczeniową. Nie wolno dodawać naprężeń równoważnych, amplitud zredukowanych ani wyników miar obliczonych osobno dla obu obciążeń. Przejście do trwałości wymagałoby pełnej pętli $\boldsymbol\sigma$–$\boldsymbol\varepsilon$ i modelu skalibrowanego dla temperatury, nie arbitralnego progu amplitudy.^17,22

Fazę i położenie maksimum zawężają zsynchronizowane pomiary temperatury lub strumienia ciepła oraz odkształcenia w przewidywanym maksimum, wraz ze stanem chłodziwa, wykonane w co najmniej dwóch stanach rozdzielających wejścia. Liczby pozostają wyłącznie demonstracją fazy i nieobserwowanej lokalizacji.^13,17,21,22

Hierarchia modeli, walidacja i niepewność

Kolejne poziomy łączą impedancje, postacie, źródło, wymianę ciepła i konstrukcję. Redukcja jest ważna, gdy zachowuje mechanizm decyzji.^3,8

Weryfikacja, walidacja i zakres użycia

Weryfikacja sprawdza implementację, jednostki, znaki, bilanse i rafinację. Walidacja porównuje fazę, postać, tempo, lokalny strumień i odpowiedź konstrukcji z niezależnymi danymi.¹³ Zakres użycia wymienia stany, konfiguracje, pasma, amplitudy i pytania; model zimny, próbka, komponent i zespół nie dzielą automatycznie jednego zakresu.¹⁵

Budżet niepewności i identyfikowalność

Niepewność $m_{c,k}$ dla modu odseparowanego — albo funkcjonału odpowiedzi $K$ — obejmuje pomiar, parametry, model, numerykę i reprezentatywność konfiguracji. Dla skalarnego $m=m(\mathbf z)$ pierwsze przybliżenie ma postać

$$u_m^2 \approx \mathbf J\,\mathbf U_z\,\mathbf J^{\mathsf T} +u_{\mathrm{model}}^2 +u_{\mathrm{num}}^2, \qquad \mathbf J=\frac{\partial m}{\partial\mathbf z}.$$

Macierz $\mathbf U_z$ zachowuje kowariancje, a brak wiedzy o modelu nie staje się losowym rozrzutem wejść.¹⁴ Równoważne pary opóźnienia i impedancji rozdziela dodatkowy sygnał pomiarowy, stan albo badanie komponentu. Przykład $m_{c,k}=3 \mathrm{s^{-1}}$ przy niepewności rozszerzonej $5 \mathrm{s^{-1}}$ nie daje dodatniej dolnej granicy.

Instrumentacja, obserwowalność i identyfikacja przyczynowa

Przewód, wnęka, membrana, okablowanie, filtr, przetwornik i zegar tworzą $H_i(\omega)$. Rejestr wiąże się z polem dopiero po ustaleniu amplitudy, fazy, liniowości i stanu toru.⁴

Od kanału do obserwowalnego modu

Dla kilku kanałów model pomiaru można zapisać

$$\widehat{\mathbf y} =\mathbf H(\omega)\, \mathbf C\boldsymbol\Phi(\omega)\, \widehat{\mathbf A} +\widehat{\mathbf n}.$$

Macierze $\boldsymbol\Phi,\mathbf C,\mathbf H$ opisują postacie, próbkowanie i tory. Rząd oraz uwarunkowanie $\mathbf H\mathbf C\boldsymbol\Phi$ wyznaczają obserwowalność; kanały we wspólnym węźle lub ze wspólną nieznaną dynamiką nie są niezależne.

Przewód pomiarowy może utworzyć własne maksimum lub opóźnienie, zatem korekcja samej amplitudy nie wystarcza do obliczania mocy. Błąd synchronizacji $\Delta t$ daje

$$\Delta\varphi=\omega\Delta t=2\pi f\Delta t.$$

Dla dydaktycznych $f=5 \mathrm{kHz}$ i $\Delta t=10 \mu\mathrm s$ błąd wynosi $0{,}314 \mathrm{rad}$, około $18^\circ$, więc zegar i latencja należą do budżetu fazy.^4,14

Wskaźnik ciepła ogranicza $P_q$, ciśnienie z prędkością odpływ, drganie transfer konstrukcyjny, a historia cieplna skutek. Emisja optyczna wymaga kalibracji; wspólna częstość nie określa kierunku mocy.

Widmo, obwiednia, koherencja i falsyfikacja

Widmo mocy lokalizuje aktywność w częstotliwości, spektrogram pokazuje jej ewolucję, a widmo wzajemne — fazę. Koherencja

$$\gamma_{xy}^2(f) =\frac{|S_{xy}(f)|^2} {S_{xx}(f)S_{yy}(f)}$$

mierzy liniową współzależność w wybranych oknach, lecz nie kierunek przyczynowy. Wspólne wymuszenie stanowiska może dać wysoką wartość bez sprzężenia spalania, a zależność nieliniowa — niską wartość mimo związku fizycznego.

Małosygnałowe tempo wyznacza się tylko z odcinka o stacjonarnym tle, rozdzielonym modzie i amplitudzie poza nasyceniem:

$$\begin{aligned} \ln |A_k(t)|={}&\ln|A_{k,0}|\\ &+\sigma_{c,k}t+\varepsilon(t). \end{aligned}$$

Niepewność obejmuje okno, szum, nakładanie modów, dryf i tor; nachylenie sprawdza się niezależną sumą mocy oraz postacią.

Samowzbudzenie, wymuszenie, artefakt, drganie stanowiska, dudnienie i dryf otrzymują przed analizą falsyfikowalne przewidywania. Samowzbudzenie wymaga zgodnego znaku mocy i obwiedni; artefakt nie powinien odtwarzać skorygowanej postaci w niezależnych kanałach.

Studium rozwijane I — historyczny przypadek F-1

Etap / konfiguracja	Obserwacja	Hipotezy nadal żywe	Próba lub dowód rozróżniający	Decyzja źródłowa	Lekcja dla macierzy
wczesny wzór wtryskiwacza; wiele prób przy ograniczonym poziomie pracy	[źródło] większość prób nie ujawniła zdarzenia, lecz odnotowano siedem przypadków samorzutnych, około 10% prób 23	[rekonstrukcja] zależność od stanu; sprzężenie zasilania, wtrysku lub pola	przejście do pełnego stanu zamiast ekstrapolacji	wczesny materiał nie dowodził odporności	wynik ograniczony do zbadanych stanów
pełny poziom pracy	[źródło] wróciły zdarzenia samorzutne z automatycznym odcięciem; 28 czerwca 1962 r. utracono silnik 12	odpowiedź źródła, droga zasilania i organizacja pola	reprodukcja w stanie, w którym zdarzenie wróciło; publiczny zapis nie daje rozłożonych pól	problem uznano za krytyczny 12	„tempo, przejście i obciążenie” otwarte; [interpretacja energetyczna] znak pełnego licznika dodatni podczas zdarzenia, bez rozdziału członów
„Project Go”; równoległe zmiany wtrysku i przegród	[źródło] badano 15 układów przegród i 14 konfiguracji wtryskiwacza 23	zmiana $P_q$, postaci, strat ściennych albo ich kombinacji	porównania konfiguracji, lecz zmiany równoległe pozostawiają degenerację przyczyn	iteracyjne zawężanie do konfiguracji końcowej 23	otwarte razem „źródło i przepływ” oraz „granice, absorber i konstrukcja”
trzeci rok; brak zdarzeń samorzutnych	[źródło] spokojna praca nie wyczerpała oceny 23	układ odporny albo metastabilny i podatny na skończony bodziec	dziewięć ocen z kontrolowaną perturbacją przeszło w niestateczność 23	utrzymano test powrotu, nie poprzestano na braku zdarzenia	osobne komórki „brak samorzutnego zdarzenia” i „powrót po perturbacji”
konfiguracja końcowa	[źródło] pracowała spokojnie i wygaszała badaną odpowiedź 23	publicznie nierozdzielone wkłady wtrysku i przegród	różnica między przebiegiem spontanicznym a odpowiedzią po perturbacji	przyjęcie konfiguracji wyłonionej w iteracjach 12,23	wynik obejmuje zbadane stany i perturbacje; brak surowych pól blokuje tylko liczbowe $P_q$, straty i $m_{\mathrm{wyk}}$

[źródło] Chronologia F-1 pokazuje, że spokojna praca i odporność po perturbacji są odrębnymi twierdzeniami: program zachował próbę powrotu także po zaniku zdarzeń samorzutnych.^12,23 [rekonstrukcja] Równoległe zmiany wtrysku i przegród nie pozwalają przypisać wyniku jednemu mechanizmowi. [interpretacja energetyczna] Źródła potwierdzają zmianę zachowania końcowej konfiguracji, lecz bez surowych pól nie rozdzielają $P_q$, strat ani postaci; dowód pozostaje ograniczony do zbadanej konfiguracji i stanów.^1,2

Studium rozwijane II — publiczny dowód dla J-2X

Etap / konfiguracja	Obserwacja	Hipotezy nadal żywe	Próba lub dowód rozróżniający	Decyzja źródłowa	Lekcja dla macierzy
plan i metoda E10001 od czerwca 2011 r.	[źródło] raport 2013 opisał komorę z wnękami i przegrodami, perturbacje, kryteria zdarzeń samorzutnych, spalania nierównomiernego i przejść oraz dynamikę kanałów 4	[rekonstrukcja] niestateczność po perturbacji, odpowiedź szerokopasmowa, dyskretna, przejściowa albo artefakt toru	rozdzielenie klas obserwacji oraz charakterystyka portów i doboru czujników 4	ustanowiono metodę wielokryterialną, nie jedno kryterium widmowe	„stan, granica i konfiguracja”, „tempo, przejście i obciążenie”, model toru
E10001: 21 prób, czerwiec 2011–wrzesień 2012	[źródło] pierwsza część pokrycia programowego 11	hipotezy właściwe obserwowanym stanom; kilka warunków niepobudzonych	powrót po wykonanych perturbacjach i osobna ocena przebiegów samorzutnych 11	brak zaobserwowanej dynamicznej niestateczności w objętym zbiorze 11	wyniki perturbacyjne dodatnie tylko w pobudzonych komórkach; pozostałe są luką pokrycia
E10002: 13 prób, luty–wrzesień 2013; E10003: 12 prób, listopad 2013–kwiecień 2014	[źródło] pokrycie rozszerzono do trzech silników i łącznie 46 prób 11	rozrzut między egzemplarzami oraz zależność od warunku nadal podlegają zakresowi prób	zgodność klas wyników między trzema silnikami 11	wszystkie trzy spełniły wymaganie szerokopasmowe 11	„domknięcie i pokrycie”: większa reprezentatywność, nie dowód stanów niebadanych
obserwacje dyskretne	[źródło] w pewnych warunkach wystąpiły małe składowe dyskretne 11	własny mod, wymuszenie lub resztkowa odpowiedź toru; źródło publiczne nie rozdziela ich energetycznie	ocena amplitudy i ryzyka oddzielona od wyniku szerokopasmowego 11	minimalne ryzyko bezpiecznej pracy w objętych warunkach 11	brak pól ogranicza „energię, postać i normę”; decyzję o obciążeniu zachowuje się jako źródłową
rozruch i wyłączenie	[źródło] oscylacje przejściowe analizowano osobno i nie wykazano zasadniczych problemów 11	wymuszenie tłem albo chwilowa odpowiedź własna	historia czasowa zamiast uśrednienia z odcinkiem ustalonym	osobny pozytywny wniosek dla badanych przejść 11	osobna komórka „przejście”; nie zastępuje tempa quasi-stałego

[źródło] J-2X rozdzielał wymaganie szerokopasmowe, małe składowe dyskretne, historie przejściowe oraz odpowiedź po perturbacji, uwzględniając dynamikę portów i czujników.^4,11 [rekonstrukcja] Są to cztery różne twierdzenia: brak jednej klasy sygnału nie rozstrzyga pozostałych. [interpretacja energetyczna] Publiczny materiał pokazuje spójność wyniku między trzema egzemplarzami, lecz nie daje pól do ponownego obliczenia $N_{c,K}$, $G_W$ ani pełnej macierzy transmitancji; jawnie niepokryte pozostają warunki bez wykonanej perturbacji i stany poza zbiorem prób.^4,11

Studium rozwijane III — syntetyczny napęd stały podczas pracy

To studium modelowe jednego odseparowanego modu sprawdza rachunek i logikę falsyfikacji, nie opisuje istniejącego silnika ani materiału. Fizyczną inspiracją jest obserwacja, że częstości modów podłużnych i mechanizmów przepływowych mogą zbliżać się podczas pracy.¹⁸

Model, przewidywania i syntetyczny zapis

Obserwacja początkowa. Ton dryfuje z około $417$ do $333 \mathrm{Hz}$, a obwiednia rośnie. Konkurują: mod wolnej objętości, rezonans przewodu, wymuszenie przepływowe i zmiana odpowiedzi powierzchni.

W przybliżeniu jednowymiarowym przyjmujemy wyłącznie efektywną długość akustyczną $L_a$, nie geometrię urządzenia. Dla $c=1000 \mathrm{m\,s^{-1}}$:

$$\begin{aligned} f_1(t)&=\frac{c}{2L_a(t)},\\ f_1(0)&=\frac{1000}{2\times1{,}20} =417 \mathrm{Hz},\\ f_1(t_2)&=\frac{1000}{2\times1{,}50} =333 \mathrm{Hz}. \end{aligned}$$

Jednostką jest herc. Mod powinien zachować podłużną postać fazową podczas dryfu; rezonans przewodu zależałby od typu toru.

Cross-PSD wskaźnika ciepła z $p'$, po korekcji torów, jest tu wyłącznie jakościowym wskaźnikiem fazy i koherencji; nie stanowi operatora kalibracji do watów. Do syntetycznego bilansu jawnie przyjmujemy jako założone wejścia $P_{q,1}=180 [155;205] \mathrm W$ oraz $P_{q,2}=320 [285;355] \mathrm W$. Wartości te nie są wyznaczone przez parę sygnałów. Ich rzeczywiste wyznaczenie wymagałoby przestrzennej rekonstrukcji $\dot q'(\mathbf x,t)$, poprawnej kwadratury całki modalnej z $p'(\mathbf x,t)$, określonego pasma i normalizacji oraz propagacji kalibracji i jej kowariancji. Współliniowość $p'$ i $u_t'$ pozostawia ponadto nierozdzielone $R_p$ i $R_u$; cross-PSD nie identyfikuje ich wynikowej mocy.

Energia układu sprzężonego wynosi w obu oknach $E_{c,k}=20 \mathrm J$, przedział $19$–$21 \mathrm J$. Pozostałe człony mocy dają łącznie odpływ i dyssypację $260 \mathrm W$, przedział $240$–$280 \mathrm W$. Centralnie $P_o=90$, $P_{s,\mathrm{out}}=20$, $P_{\mathrm{abs}}=30$, $D_w=40$, $D_{v,\mathrm{rest}}=30$, $D_p=20$ i $D_s=30 \mathrm W$; $D_v=50 \mathrm W$, więc $D_p$ nie jest liczone drugi raz. Pozostałe estymaty są zgodne z zerem w łącznej niepewności. Wartości centralne dają:

$$\sigma_{c,k,1} =\frac{180-260}{2\times20} =-2{,}0 \mathrm{s^{-1}}.$$

$$\sigma_{c,k,2} =\frac{320-260}{2\times20} =+1{,}5 \mathrm{s^{-1}}.$$

Ponieważ $\mathrm W/\mathrm J=\mathrm{s^{-1}}$, jednostki są poprawne. Skrajne kombinacje założonych przedziałów dają odpowiednio $\sigma_{c,k,1}\in[-3{,}29;-0{,}83] \mathrm{s^{-1}}$ i $\sigma_{c,k,2}\in[0{,}12;3{,}03] \mathrm{s^{-1}}$: warunkowy znak zmienia się między oknami także po uwzględnieniu tej syntetycznej niepewności. Jest to rachunek dla założonych wejść mocy, nie estymacja mocy z cross-PSD; niezależny test obwiedni sprawdza jedynie zgodność przewidywanego tempa.

Test hipotez, ograniczenia i decyzja

Syntetyczny pomiar. Cztery kanały ciśnienia i dwa tory z różnymi przewodami, po korekcji $H_i$, odtwarzają postać oraz ciągły dryf. Niezależna od cross-PSD regresja $\ln|A_k|$ daje $-1{,}9\pm0{,}3 \mathrm{s^{-1}}$ i $+1{,}6\pm0{,}4 \mathrm{s^{-1}}$, zgodnie z warunkowym znakiem bilansu wykorzystującego założone wejścia $P_q$.

Hipoteza	Sygnatura w fazie, postaci, obwiedni i kanałach niezależnych	Pomiar rozstrzygający i wynik syntetyczny	Rozstrzygnięcie	Kategoria dowodu
zmiana odpowiedzi powierzchni	zmiana fazy $p'$–wskaźnik ciepła przy zachowanej postaci; trend obwiedni zgodny ze zmianą założonego bilansu; odpowiedź obecna poza jednym torem	skorygowana cross-PSD ciepło–ciśnienie zmienia relację fazową, a niezależne $\ln|A_k|$ zmienia znak; brak przestrzennego $\dot q'$ nie pozwala wyznaczyć $P_q$	zgodna z danymi, nierozstrzygnięta; $R_p$ i $R_u$ nie są osobno zidentyfikowane	źródło i przepływ
wymuszenie lub sprzężenie przepływowe	faza związana z markerem przepływu, maksimum obwiedni przy zbliżeniu częstości; pole może zachować postać akustyczną	marker przepływu potwierdza przecięcie, ale nie wyjaśnia sam znaku mocy ani wzrostu po odjęciu wymuszenia	nieodrzucona jako wyzwalacz lub droga sprzężenia; odrzucona jako samodzielne wyjaśnienie danych	źródło i przepływ; tempo i przejście
dryf albo przejście modu	dryf zachowuje ciągłość fazy i postaci; przejście daje zmianę podprzestrzeni lub skok korelacji	cztery kanały zachowują podłużny wzór podczas $417\to333 \mathrm{Hz}$	dryf wspierany, przełączenie modu odrzucone w tych oknach	energia, postać i norma; pokrycie
artefakt toru	faza i maksimum zależą od typu przewodu, bez wspólnej postaci i niezależnej obwiedni	po korekcji dwa różne tory i pozostałe kanały dają ten sam dryf oraz tempo	odrzucona w syntetycznym przykładzie	stan, granica i konfiguracja; domknięcie

Ciągłość postaci odróżnia dryf od przełączenia, różne tory ograniczają artefakt, a obwiednia sprawdza tempo przewidziane przez bilans z jawnie założonym $P_q$. Cross-PSD ogranicza możliwą relację fazową, lecz bez rekonstrukcji przestrzennej nie lokalizuje przyczyny po stronie źródła i nie wyznacza watów. Marker przepływu nie rozstrzyga jego roli, a $R_p$ i $R_u$ pozostają nierozdzielone.^7,8 Warunkowy wniosek dla dwóch okien brzmi: przy przyjętych wejściach dodatnie $\sigma_{c,k,2}$ oznacza lokalne narastanie; hipoteza zmiany odpowiedzi powierzchni pozostaje zgodna z danymi, lecz nierozstrzygnięta. Zmiana objętości, powierzchni, cząstek, dyszy lub materiału wymaga ponownej oceny właściwych członów.¹⁸

Oś cieplno-konstrukcyjna tego samego studium

Drugi werdykt wymaga przeprowadzenia tego samego zrekonstruowanego pola do lokalnej miary uszkodzenia. Poniższy ciąg pozostaje celowo syntetyczny: liczby nie opisują materiału, geometrii ani progu rzeczywistego napędu. W drugim oknie rekonstrukcja postaci daje w ustalonym punkcie $\mathbf x_*$ amplitudę $|\widehat p(\mathbf x_*)|=24 \mathrm{kPa}$. Niezależnie od cross-PSD przyjmujemy jako wejście lokalnego bilansu demonstracyjny operator $|H_{qp}(\mathbf x_*)|=0{,}30 \mathrm{W\,m^{-2}Pa^{-1}}$, wraz z jego fazą. Stąd założony lokalny strumień wynosi

$$|\widehat q_w''(\mathbf x_*)| =|H_{qp}\widehat p|=7{,}2 \mathrm{kW\,m^{-2}}.$$

Nie jest to wartość uzyskana z pary sygnałów ani parametr rzeczywistej ścianki. W dowodzie rzeczywistym $H_{qp}$ musiałby pochodzić z niezależnie zwalidowanego mapowania pola na strumień, z zachowaniem położenia, fazy i niepewności.

Dla wyłącznie demonstracyjnego modelu półprzestrzeni przyjmijmy $f=333 \mathrm{Hz}$, $\alpha=1{,}0\times10^{-5} \mathrm{m^2s^{-1}}$ i $k_m=20 \mathrm{W\,m^{-1}K^{-1}}$. Wówczas

$$\delta_T=\sqrt{\frac{2\alpha}{2\pi f}} =9{,}78\times10^{-5} \mathrm m, \qquad |\widehat T_s| =\frac{|\widehat q_w''|\delta_T}{\sqrt2k_m} =0{,}0249 \mathrm K.$$

Wynik lokalizuje zmianę temperatury i jej gradient, nie maksimum naprężenia. To ostatnie musi wynikać z rozwiązania konstrukcyjnego z rzeczywistymi więzami, geometrią i własnościami zależnymi od temperatury.

Odpowiedź konstrukcyjną śledzimy w tym samym $\mathbf x_*$ i w stałym kierunku jednostkowym $\mathbf n_*$, przez podpisaną składową

$$\sigma_{nn}(t)=\mathbf n_*^{\mathsf T} \boldsymbol\sigma(\mathbf x_*,t)\mathbf n_*.$$

Demonstracyjne liniowe operatory, niebędące własnościami materiału ani progiem dopuszczenia, dają amplitudy $a_p=|\widehat\sigma_{nn,p}|=0{,}96 \mathrm{MPa}$ i $a_T=|\widehat\sigma_{nn,T}|=0{,}70 \mathrm{MPa}$. Dla zachowanej różnicy faz $\Delta\varphi=60^\circ$ zespolone, podpisane składowe można dodać przed obliczeniem amplitudy:

$$a_{nn}=|\widehat\sigma_{nn,p}+\widehat\sigma_{nn,T}| =\sqrt{a_p^2+a_T^2+2a_pa_T\cos\Delta\varphi} =1{,}44 \mathrm{MPa}.$$

Nie jest to naprężenie równoważne. W analizie wieloosiowej należałoby najpierw zsumować pełne historie tensorowe ciśnieniową i cieplną, a dopiero potem zastosować wybraną miarę uszkodzenia.

Niepewność przeprowadzamy przez oba etapy, zamiast dopisywać ją dopiero do werdyktu. Przy niezależnych standardowych niepewnościach względnych $u_q/q=10\%$, $u_\alpha/\alpha=10\%$, $u_k/k=5\%$ i $u_f/f=3/333$ liniowa propagacja dla temperatury daje

$$\left(\frac{u_T}{T_s}\right)^2 =\left(\frac{u_q}{q}\right)^2 +\frac14\left(\frac{u_\alpha}{\alpha}\right)^2 +\left(\frac{u_k}{k_m}\right)^2 +\frac14\left(\frac{u_f}{f}\right)^2, \qquad u_T=0{,}0031 \mathrm K.$$

Dla $u_{a_p}=0{,}10 \mathrm{MPa}$, $u_{a_T}=0{,}12 \mathrm{MPa}$ i $u_\varphi=10^\circ=0{,}175 \mathrm{rad}$, bez znanych kowariancji, różniczkowanie powyższego wzoru daje

$$\begin{aligned} u_{a_{nn}}^2={}& \left(\frac{a_p+a_T\cos\Delta\varphi}{a_{nn}}u_{a_p}\right)^2 +\left(\frac{a_T+a_p\cos\Delta\varphi}{a_{nn}}u_{a_T}\right)^2\\ &+\left(\frac{a_pa_T\sin\Delta\varphi}{a_{nn}}u_\varphi\right)^2, \qquad u_{a_{nn}}=0{,}15 \mathrm{MPa}. \end{aligned}$$

Brak kowariancji jest tu elementem założenia syntetycznego; w dowodzie rzeczywistym wspólna rekonstrukcja pola zwykle je tworzy i macierz kowariancji musi zostać zachowana. Otrzymany zapis $a_{nn}=1{,}44\pm0{,}15 \mathrm{MPa}$ dotyczy jednej składowej i nie stanowi zapasu trwałości.

Jawnie wybraną miarą uszkodzenia byłaby energia niesprężyście rozproszona w ustabilizowanym cyklu,

$$W_d=\oint\boldsymbol\sigma:\mathrm d\boldsymbol\varepsilon^{\mathrm{inel}}, \qquad [W_d]=\mathrm{J\,m^{-3}\,cykl^{-1}}.$$

Nie można jej jednak wyznaczyć z dostępnych danych syntetycznych: brakuje pełnej historii tensora odkształcenia niesprężystego, temperaturowo zależnego i zwalidowanego prawa konstytutywnego, danych o stanie początkowym i liczbie cykli, kalibracji relacji $W_d$–uszkodzenie oraz zatwierdzonego kryterium dopuszczenia. Niepewność końcowej miary nie jest więc skończonym przedziałem, który można porównać z progiem; łańcuch dowodowy urywa się przed $W_d$.

Werdykt osi cieplno-konstrukcyjnej: niewykazane. Zrekonstruowane pole i założony lokalny strumień prowadzą przez rozwiązanie cieplne do jawnie zdefiniowanej, podpisanej składowej naprężenia z propagacją niepewności, ale nie dostarczają miary uszkodzenia ani podstawy jej akceptacji. Brakującym dowodem jest zwalidowana pełna historia tensorowa $(\boldsymbol\sigma,\boldsymbol\varepsilon^{\mathrm{inel}},T)$ w rzeczywistym maksimum, wraz z modelem uszkodzenia, jego niepewnością i progiem programu. Warunkowy wynik stateczności z poprzedniej części nie zmienia tego werdyktu.

Piramida dowodu, kwalifikacja i zarządzanie zmianą konfiguracji

Od wymagania do argumentu zgodności

Dowód rośnie poziomami; wynik wyższy korzysta z niższego, ale go nie zastępuje.^9,10,13,15

Poziom	Pytanie	Typ dowodu	Zależność od poziomu niższego	Czego nie wolno wnioskować
kalibracja toru lub parametru	czy wskazanie ma poprawną skalę, fazę i niepewność?	wzorzec, identyfikowalność, model toru, kowariancja	ustanawia znaczenie danych wejściowych	że skalibrowany kanał obserwuje całe pole lub przyczynę
eksperyment elementarny	czy wydzielony mechanizm działa i model oddaje jego trend?	próba materiałowa lub fizyczna w kontrolowanym paśmie i stanie	używa skalibrowanych torów i parametrów	że mechanizm dominuje w podzespole albo zespole
podzespół	czy mechanizm zachowuje znak, fazę i skalę po dodaniu interfejsów?	bilans, FRF, pola i walidacja na niewykorzystanym stanie	łączy modele elementarne i ich niepewność	że otoczenie zespołu nie zmieni granic ani sprzężeń
zespół w reprezentatywnym środowisku	czy sprzężony układ spełnia wymaganie w zbadanych stanach?	próby zespołu, bilans domknięcia, perturbacje i historia obciążenia	opiera interpretację na zwalidowanych podzespołach	że pokryto całą obwiednię lub każdy egzemplarz
kwalifikacja obwiedni	czy reprezentatywna konfiguracja wykazała wymaganie w zadanym zakresie użycia?	mapa pokrycia, dolne granice zapasu, przejścia i kontrola konfiguracji	scala dowody zespołu z ekstrapolacją ograniczoną niepewnością	że każda sztuka wykonana później jest zgodna
akceptacja konkretnego egzemplarza	czy ta sztuka odpowiada konfiguracji i kryteriom odbioru?	zapisy produkcji, NDE, badania odbiorcze i ślad konfiguracji	korzysta z uprzednio zakwalifikowanej obwiedni	że rozszerzono kwalifikację lub zwalidowano nowy mechanizm

Weryfikacja sprawdza, czy rozwiązanie spełnia specyfikację i czy obliczenia wykonano poprawnie; walidacja — czy model wystarczająco reprezentuje rzeczywistość w zadanym zakresie użycia. Kwalifikacja wykazuje zgodność reprezentatywnej konfiguracji w określonej obwiedni, a akceptacja potwierdza zgodność konkretnego egzemplarza z kryteriami odbioru i konfiguracją zakwalifikowaną.^13,15

Konfiguracja, cykl życia i regresja

Ślad ECSS prowadzi od wersji sprzętu i oprogramowania przez model oraz dane do decyzji.¹⁶ Materiał, proces, NDE, magazynowanie, montaż lub czujnik wskazują człony bilansu i transferu wymagające regresji.^9,10,16

Regresja zaczyna się od analizy wpływu i powtarza dowód zmienionego mechanizmu. Zmiana tezy, architektury, dominującego sprzężenia, klasy obciążenia lub zakresu użycia wymaga pełnej bramki. Werdykt podaje konfigurację, obwiednię, przejścia i niewzbudzone obszary.^15,16

Synteza — od obserwacji do werdyktu o marginesie

Drabina wnioskowania i werdykt

Ciąg dowodu prowadzi od obserwacji przez pola, moc i tempo do historii cieplno-mechanicznej. Wspólna księgowość nie usuwa różnicy przyczynowej między klasami napędu:

Oś porównania	Napęd ciekły	Napęd stały
magazyny i opóźnienia	zasilanie, wtrysk, przygotowanie faz, mieszanie i reakcja	pamięć powierzchni, warstwa przyścienna, cząstki, dysza i zmienna objętość
model źródła	wielowejściowe $B_kF_{\chi q}$ prowadzące do przestrzennego $\dot q'$	zespolone $R_p,R_u$, operator transportu i całka $P_q$
sygnały rozróżniające	zsynchronizowany dopływ, wtrysk/ciepło, postać i obwiednia	$p',u_t'$, wskaźnik powierzchni, ewolucja postaci i niezależne tory
typowa degeneracja	zasilanie i spalanie dają podobną fazę; mieszają się mody komory, dopływu i kanału	współliniowe $p',u_t'$ nie rozdzielają $R_p/R_u$; dryf przypomina zmianę źródła

Werdykt wyrasta z jednego łańcucha. Po korekcji toru i opisie stanu dane rozłożone wskazują postać $k$ albo podprzestrzeń $K$; faza wyznacza moc ze znakiem, a niezależna obwiednia sprawdza jej wynik. Odtworzona dynamika zasila historię cieplno-mechaniczną, a niepewność i pokrycie ograniczają zakres twierdzenia. Przerwanie łańcucha oznacza wynik węższy albo niewykazany.^13,14,15

Ostatnia kontrola spójności i bezpieczeństwa

Przyjęta granica publikacji pomija instrukcje wykonawcze, mapy słabych miejsc i procedury wymuszania współczesnego uzbrojenia. Nieznanego, uszkodzonego lub odnalezionego elementu uzbrojenia ani napędu nie wolno dotykać, przenosić ani próbować samodzielnie „zabezpieczać”. Należy odsunąć ludzi, zadzwonić pod 112 i wykonywać polecenia służb. Nie wolno go rozłączać, wiercić, ciąć, ogrzewać, chłodzić, płukać ani pobierać próbek.

Margines staje się twierdzeniem dopiero wtedy, gdy bilans, obserwacje, obciążenie i zakres dotyczą tego samego układu.

Źródła