Akumulator jako zależna od stanu impedancja boczna
Akumulator cieczowo-gazowy nie ma jednej „pojemności”, która odpowiada na wszystkie pytania projektowe. Jest zależną od stanu impedancją boczną: jego odpowiedź zmieniają ciśnienie i temperatura gazu, położenie separatora, amplituda oraz częstotliwość wymuszenia, opory przyłącza, tarcie, ograniczniki skoku i impedancja pozostałej sieci. Ta sama podatność może magazynować energię, przyjmować powolną zmianę objętości, przejmować okresowy strumień albo łagodzić pojedynczy udar, lecz są to cztery różne zadania. Korzystają z niej w innych skalach czasu, amplitudy i sprzężenia z siecią.
Dlatego każda analiza w tym artykule zachowuje jeden łańcuch dowodowy:
wymaganie → stan i granica układu → model → przewidywana sygnatura → pomiar → niepewność → kryterium → decyzja → regresja.
Przerwanie łańcucha daje pozorną pewność. Sama zgodność objętości nie dowodzi tłumienia. Sama mała pulsacja w jednym punkcie nie dowodzi, że energia nie przeniosła się w inne miejsce. Sam dodatni odczyt ciśnienia gazu nie dowodzi dostępności ciśnienia hydraulicznego. Sama próba szczelności nie dowodzi pełnego skoku separatora.
Bezpieczeństwo. Nieznany akumulator lub tłumik należy traktować jak urządzenie naładowane i potencjalnie skażone. Nie wolno go nacinać, ogrzewać, przebijać, odkręcać ani odpowietrzać bez identyfikacji, zatwierdzonej procedury, odizolowania obu stron i niezależnego potwierdzenia stanu bezpiecznego. Element znaleziony przy amunicji, wraku lub niezidentyfikowanym sprzęcie pozostawia się specjalistom EOD; podejrzenie skażenia chemicznego, biologicznego lub promieniotwórczego wymaga procedur CBRN. Artykuł wyjaśnia fizykę i sposób dowodzenia zgodności, nie podaje parametrów ani procedur ingerencji we współczesne uzbrojenie.
Cztery funkcje, cztery miary powodzenia
Pierwszym wymaganiem nie jest „zastosować akumulator”, lecz nazwać usługę, którą ma wykonać układ. Dokument JPL o akumulatorach pompowanych pętli cieplnych rozdziela przyjmowanie rozszerzalności cieczy, ochronę pompy przed kawitacją i zapobieganie wrzeniu; już ten cywilny przykład pokazuje, że wspólny element nie scala kryteriów funkcjonalnych.jpl-loop-functions
| Funkcja | Dominująca skala | Wielkość wejściowa | Właściwa metryka | Typowa porażka ukryta przez inną metrykę |
|---|---|---|---|---|
| magazynowanie energii | od sekund do godzin | krzywa absolutnego ciśnienia i objętości | użyteczna praca między dwoma stanami, przy zadanym ciśnieniu odbiornika | duża objętość geometryczna, lecz mała objętość oddana przed osiągnięciem ciśnienia minimalnego |
| kompensacja objętości | minuty, godziny, cykle cieplne | bilans rozszerzalności całej pętli | pozostały zapas skoku i przedział ciśnienia w całej obwiedni temperatur | poprawny stan w temperaturze napełnienia, lecz dojście do ogranicznika na gorąco albo na zimno |
| tłumienie pulsacji | wiele okresów wymuszenia | widmo źródła i impedancje gałęzi | transmitancja zespolona ciśnienia/przepływu w paśmie oraz średnia moc rozpraszana | mała amplituda w jednym punkcie, lecz rezonans albo wzrost amplitudy w innym |
| ograniczanie udaru | pojedyncze przejście fali i odbicia | historia zmiany prędkości, czas zamknięcia, granice przewodu | szczyt, czas narastania, impuls i brak kawitacji/rozdzielenia kolumny | wystarczająca objętość statyczna, lecz zbyt duża inertancja przyłącza, aby objętość była dostępna na czas |
Magazyn energii rozlicza się w dżulach i dostępnej objętości. Kompensator — z bilansu objętości i marginesu skoku. Tłumik — z zespolonej odpowiedzi częstotliwościowej. Ogranicznik udaru — z przebiegu przejściowego w rozłożonym przewodzie. Można spełnić jedną z tych funkcji i jednocześnie zawieść trzy pozostałe.
Granica układu i jednoznaczny bilans objętości
Przyjmijmy sztywną w danej chwili komorę roboczą o objętości $V_c$, rozdzieloną na objętość gazu $V_g$ i objętość cieczy $V_l$. Objętości przewodów poza komorą nie należą do $V_c$. Dla separatora bez istotnej objętości własnej:
$$V_c(T,p)=V_g+V_l.$$
Jeżeli separator ma objętość materiału $V_s$ zmienną w zakresie dokładności modelu, bilans brzmi $V_c=V_g+V_l+V_s$. Pominięcie mieszka, trzpienia lub tłoka jest dopuszczalne dopiero po oszacowaniu błędu, a nie na mocy nazwy urządzenia.
Ustalamy znaki:
- $Q_A>0$ oznacza strumień cieczy z węzła sieci do akumulatora;
- $dV_l/dt=Q_A$ i, dla stałego $V_c$, $dV_g/dt=-Q_A$;
- $V_d(t)=V_l(t_1)-V_l(t)$ oznacza objętość oddaną odbiornikowi od chwili $t_1$, zatem podczas rozładowania $dV_d>0$ i $dV_d=dV_g$;
- przemieszczenie separatora $x$ rośnie umownie w kierunku zwiększania $V_l$; jeśli rzeczywista geometria ma znak przeciwny, zmienia się definicję $x$, nie równania w połowie obliczeń.
Bilans całej zamkniętej pętli musi obejmować masę cieczy $m_l$, jej gęstość $\rho_l(T,p)$, objętość wszystkich przestrzeni poza kompensatorem $V_e(T,p)$, przeciek lub ubytek $V_{loss}$ oraz objętość cieczy w kompensatorze. W poniższej konwencji $m_l$ jest masą zapisaną dla stanu napełnienia przed osobno ewidencjonowanym ubytkiem, a $V_{loss}$ jest skumulowaną objętością utraconą, sprowadzoną do tej samej temperatury i gęstości co $m_l/\rho_l$:
$$V_l(T,p)=\frac{m_l}{\rho_l(T,p)}-V_e(T,p)-V_{loss}.$$
Jeżeli zamiast tego $m_l(t)$ oznacza bieżącą, rzeczywiście pozostałą masę po ubytku, człon $V_{loss}$ należy usunąć: $V_l=m_l(t)/\rho_l-V_e$. Tych dwóch konwencji nie wolno łączyć, bo policzyłoby to ubytek dwa razy. Równanie odróżnia rozszerzanie cieczy od rozszerzania obudowy. Współczynnik rozszerzalności cieczy użyty bez objętości i rozszerzalności całej pętli nie tworzy jeszcze bilansu. W układzie z wieloma temperaturami należy sumować odcinki: $\sum_i m_i/\rho_i(T_i,p_i)$ i $\sum_j V_{e,j}(T_j,p_j)$, zamiast podstawiać jedną „temperaturę układu”.
Trzy rodzaje ciśnienia
Równanie stanu gazu, kryterium wrzenia i kryterium kawitacji wymagają ciśnienia absolutnego $p_{abs}$, liczonego od próżni. Ciśnienie manometryczne $p_{gauge}=p_{abs}-p_{amb}$ zależy od aktualnego ciśnienia otoczenia $p_{amb}$. Różnica użyteczna $\Delta p_u=p_{up}-p_{down}$ jest ciśnieniem dostępnym do pokonania strat i wykonania pracy między dwiema granicami. Te trzy wartości bywają liczbami podobnego rzędu, ale nie są wymienne.
Manometr po stronie gazowej nie zastępuje pomiaru hydraulicznego. Separator przenosi ciśnienie przez równowagę sił, a więc $p_l-p_g$ może obejmować sprężystość, tarcie, ciężar, przyspieszenie i kontakt z ogranicznikiem. Załącznik Revision B with Change 1 do NASA-STD-8719.24 w kontekście bezpieczeństwa ładunków NASA wprost zabrania traktowania wskazania gazu jako równoważnego wskazaniu ciśnienia hydraulicznego. Nie jest to ogólna norma wszystkich akumulatorów, lecz trafna kontrola błędu obserwowalności.nasa-annex-gauge
Gazowa sprężyna: stan, praca i podatność
Dla gazu idealnego
$$p_gV_g=m_gR_sT_g=n_{mol}R_uT_g,$$
gdzie $m_g$ jest masą gazu, $n_{mol}$ liczbą moli, $T_g$ temperaturą bezwzględną, a $R_s=R_u/M$ indywidualną stałą gazową danego gazu. $M$ oznacza masę molową, a uniwersalna stała gazowa ma wartość $R_u=8{,}314\,462\,618 \mathrm{J/(mol\,K)}$. Nie wolno wpisywać uniwersalnej i indywidualnej stałej pod tym samym symbolem bez zdefiniowania jednostek. Własności termofizyczne, w tym gęstość, $c_p$, $c_v$, prędkość dźwięku i krzywe nasycenia, powinny pochodzić z danych dla właściwego medium i zakresu stanu, nie z pamięciowej liczby katalogowej.nist-fluid-data
Model idealny pomija objętość własną cząsteczek i oddziaływania między nimi. Jego błąd rośnie przy dużej gęstości, w pobliżu obszaru nasycenia i tam, gdzie współczynnik ściśliwości $Z_{comp}=pV/(mR_sT)$ odchodzi od jedności. Wtedy stosuje się zweryfikowane równanie stanu rzeczywistego gazu. Korekta nie usuwa potrzeby pomiaru temperatury: gaz ściskany szybko może być cieplejszy od korpusu, choć oba czujniki później pokażą tę samą temperaturę.
Model politropowy i jego granice czasowe
Na odcinku przemiany, dla którego można przyjąć stały wykładnik $n$,
$$p_gV_g^n=K.$$
$n=1$ opisuje granicę izotermiczną gazu idealnego, gdy wymiana ciepła nadąża za zmianą objętości. $n=\gamma=c_p/c_v$ opisuje odwracalną granicę adiabatyczną. W rzeczywistym urządzeniu $n$ jest parametrem zastępczym: zależy od czasu procesu względem czasu wymiany ciepła, amplitudy, kształtu komory, przewodnictwa separatora, częstotliwości i historii. Nie musi mieć jednej wartości dla powolnego napełniania, krótkiego udaru i małej oscylacji wokół tego samego średniego stanu.
Praca wykonana przez gaz podczas rozszerzania od $V_{g1}$ do $V_{g2}$ wynosi, dla $n\ne1$:
$$W_g=\int_{V_{g1}}^{V_{g2}}p_g\,dV_g =\frac{p_{g2}V_{g2}-p_{g1}V_{g1}}{1-n} =\frac{p_{g1}V_{g1}-p_{g2}V_{g2}}{n-1}.$$
Dla $n=1$ granica jest skończona i ma osobną postać:
$$W_g=p_{g1}V_{g1}\ln\!\left(\frac{V_{g2}}{V_{g1}}\right).$$
To praca gazu, nie automatycznie energia użyteczna odbiornika. Jeżeli $p_{port}$ jest ciśnieniem cieczy na porcie akumulatora, $p_{ret}$ ciśnieniem odniesienia/strony powrotnej odbiornika, a $\Delta p_{line}$ stratą między portem i odbiornikiem, użyteczna praca hydrauliczna ponad ciśnienie powrotne wynosi:
$$W_u=\int_0^{V_{d,2}} [p_{port}(V_d)-p_{ret}(V_d)-\Delta p_{line}(V_d)]\,dV_d.$$
$p_{port}$ i $p_{ret}$ muszą być mierzone w zgodnej konwencji. Jeśli $p_{ret}$ jest atmosferą, można całkować ciśnienie manometryczne na odbiorniku; dla zamkniętej drogi trzeba użyć rzeczywistej różnicy. Gdy $p_{port}$ jest estymowane z $p_g$, należy dodatkowo uwzględnić pracę tarcia i histerezy separatora; przy bezpośrednim pomiarze $p_{port}$ wpływ ten jest już zawarty w ciśnieniu cieczy. Jednostka $\mathrm{Pa\,m^3}=\mathrm J$ nie naprawia błędnej granicy układu ani znaku.
Pochodna politropy jest ujemna po stronie gazu:
$$\frac{dp_g}{dV_g}=-\frac{np_g}{V_g}.$$
Ponieważ wzrost ciśnienia hydraulicznego zwiększa $V_l$ i zmniejsza $V_g$, małosygnałowa podatność hydrauliczna idealnego separatora jest dodatnia:
$$C_h\equiv\left.\frac{\partial V_l}{\partial p_l}\right|_0 =-\left.\frac{\partial V_g}{\partial p_g}\right|_0 =\frac{V_{g0}}{n p_{g0}}, \qquad [C_h]=\mathrm{m^3/Pa}.$$
Indeks 0 oznacza punkt pracy, nie „stan pusty”. Podatność maleje ze wzrostem ciśnienia i rośnie z objętością gazu. Akumulator jest zatem nieliniowy nawet przy idealnym separatorze. Po dojściu do ogranicznika skoku $C_h$ gwałtownie dąży do podatności samych ścian i cieczy, o wiele mniejszej. Przed tym zdarzeniem nieliniowość objawia się harmonicznymi i zależnością odpowiedzi od amplitudy.
Separator dodaje własną sztywność. Jeżeli przyrost siły sprężystej wynosi $k_s\,dx$, a $dV_l=A_{eff}\,dx$, jego lokalna podatność ciśnieniowa ma rząd $A_{eff}^2/k_s$. Łączy się ona z podatnością gazu zgodnie z rzeczywistym układem sił; nie wolno po prostu dodać dwóch „pojemności” bez wyprowadzenia. Zmienna powierzchnia czynna $A_{eff}(x)$ i histereza sprawiają, że krzywa ładowania nie musi pokrywać się z rozładowaniem.
Liniaryzacja i impedancja widziana z węzła
Wokół ustalonego punktu $(p_0,Q_0,V_0,T_0)$ zapisujemy $p=p_0+\delta p$ i $Q=Q_0+\delta Q$. Liniaryzacja jest lokalna: jej parametry należy identyfikować dla stanu, temperatury, amplitudy i kierunku ruchu reprezentatywnych dla wymagania. Historyczne analizy NASA opisywały akumulator trzema skupionymi parametrami — inertancją, podatnością i oporem — oraz jego admitancją przepływową; jest to użyteczna struktura modelu, nie uniwersalny komplet liczb.nasa-passive-rlc
Przyjmujemy konwencję zespoloną $e^{j\omega t}$ i amplitudy szczytowe. Krótki port lub przewód o długości $\ell$, polu $A$ i gęstości cieczy $\rho$ ma w pierwszym przybliżeniu inertancję hydrauliczną
$$L_h=\frac{\rho\ell}{A}, \qquad [L_h]=\mathrm{Pa\,s^2/m^3}.$$
Lokalny opór różniczkowy w punkcie pracy definiujemy jako
$$R_h=\left.\frac{\partial(\Delta p_{str})}{\partial Q}\right|_{Q_0}, \qquad [R_h]=\mathrm{Pa\,s/m^3}.$$
Dla laminarnej gałęzi może być prawie stały; dla kryzy z $\Delta p\sim Q|Q|$ zależy od $Q_0$, a przy oscylacji wokół zera pojawia się nieliniowość i harmoniczne. $R_h$ jest parametrem dyssypacji, nie dowolnym „współczynnikiem tłumienia”.
Element widziany z węzła ma kolejno szeregowe: opór i inertancję przyłącza oraz podatność komory. Przy $Q_A>0$ do akumulatora:
węzeł p ── R_h ── L_h ── C_h ── odniesienie małosygnałowe
Q_A > 0 do komoryRównania gałęzi są
$$\delta p_{node}-\delta p_c=R_h\delta Q_A+L_h\frac{d\delta Q_A}{dt}, \qquad \delta Q_A=C_h\frac{d\delta p_c}{dt}.$$
Stąd impedancja i admitancja:
$$Z_A(\omega)\equiv\frac{\widehat{\delta p}_{node}}{\widehat{\delta Q}_A} =R_h+j\omega L_h+\frac{1}{j\omega C_h}, \quad [Z_A]=\mathrm{Pa\,s/m^3}, \qquad Y_A=\frac{1}{Z_A}, \quad [Y_A]=\mathrm{m^3/(Pa\,s)}.$$
Częstotliwość własna nierozłożonej gałęzi i współczynnik tłumienia mają postać
$$\omega_0=\frac{1}{\sqrt{L_hC_h}}, \qquad f_0=\frac{\omega_0}{2\pi}, \qquad \zeta=\frac{R_h}{2}\sqrt{\frac{C_h}{L_h}}.$$
Przy $\omega_0$ reaktancje inertancji i podatności znoszą się, $Z_A=R_h$, a przepływ gałęzi może być największy. To rezonans szeregowej gałęzi, który może równocześnie tworzyć minimum — antyrezonans — transmitancji ciśnienia w węźle zasilanym źródłem strumienia. Nazwa zależy więc od obserwowanej wielkości i punktów wejścia/wyjścia.
Akumulator jest bocznikiem, a nie samodzielnym filtrem. Jeśli pozostała droga z węzła do odniesienia ma impedancję $Z_N$, źródło strumienia $\hat Q_s$ wytwarza
$$\hat p=\hat Q_s\frac{Z_NZ_A}{Z_N+Z_A}.$$
Bez gałęzi $\hat p_0=\hat Q_sZ_N$, zatem transmitancja znormalizowana wynosi
$$H(\omega)=\frac{\hat p}{\hat p_0}=\frac{Z_A}{Z_N+Z_A}.$$
Zmiana długości przewodu, warunku brzegowego albo regulatora zmienia $Z_N$, a więc również wynik tej samej części. Badania NIST inertancji rur mierzyły równocześnie amplitudę i fazę chwilowego ciśnienia oraz przepływu na obu końcach, aby wyznaczyć impedancję zespoloną i moc akustyczną; model linii transmisyjnej odtwarzał rezonans podatności zbiornika z inertancją rury.nist-inertance
Faza, amplituda i moc
Sama informacja „przepływ jest przesunięty o 90°” nie mówi, ile energii przepływa. Dla amplitud szczytowych moc zespolona wpływająca do gałęzi wynosi
$$S_A=\frac12\hat p\hat Q_A^*, \qquad \overline P_A=\frac12\operatorname{Re}(\hat p\hat Q_A^*), \qquad P_{reakt}=\frac12\operatorname{Im}(\hat p\hat Q_A^*).$$
Dla fazorów skutecznych czynnik $1/2$ znika. Idealna podatność i inertancja wymieniają energię z siecią, lecz ich średnia moc czynna jest zerowa. Dodatnia $\overline P_A$ w $R_h$ oznacza dyssypację. Mała faza przy bardzo małych amplitudach może odpowiadać znikomej mocy, a faza bliska $90°$ przy dużych amplitudach może oznaczać duży strumień energii reaktywnej i znaczny skok separatora.
Wymaganie dynamiczne powinno więc obejmować przynajmniej $|H|$, fazę $\arg H$, zakres przepływu/skoku gałęzi i moc czynną w paśmie. Te wielkości ujawniają odpowiednio: skuteczność lokalną, ryzyko sprzężenia z innym elementem, obciążenie separatora i nagrzewanie dławika.
Węzeł ciśnienia nie jest węzłem przepływu
Boczna podatność pobiera strumień, gdy pobudza ją różnica ciśnienia. W węźle fali stojącej, w którym $\delta p\approx0$, gałąź ma małe wymuszenie i może być nieskuteczna. Węzeł przepływu/prędkości bywa natomiast strzałką ciśnienia; właśnie tam bocznik może przejąć znaczny strumień. Nie ma zatem zasady „umieść tłumik w węźle przepływu” ani zasady przeciwnej bez określenia fali, warunków brzegowych i transmitancji. Najpierw wyznacza się zespolone pole ciśnienia i przepływu sieci, potem wybiera lokalizację i ponownie liczy sieć z dodaną gałęzią.
Gdy przewód staje się układem rozłożonym
Model skupiony zakłada, że informacja o ciśnieniu dociera przez element praktycznie natychmiast w interesującej skali. Dla odcinka długości $\ell$ i prędkości fali $a$ czas propagacji wynosi $\tau_p=\ell/a$, długość fali $\lambda=a/f$, a liczba falowa $k=\omega/a$. Gdy $\omega\tau_p$ przestaje być małe — praktyczna reguła wstępna to $\ell/\lambda\gtrsim0{,}1$, nie granica fizyczna — faza i odbicia wzdłuż przewodu nie mieszczą się w jednym $L_h$ i $R_h$.
W bezstratnym przybliżeniu charakterystyczna impedancja dla strumienia objętościowego wynosi $Z_{char}=\rho a/A$. Rzeczywista linia ma opory zależne od częstotliwości, podatność ścian, rozgałęzienia, objętości martwe i odbicia od zbiorników, zaworów i pomp. Wtedy stosuje się równania linii transmisyjnej, metodę charakterystyk albo zweryfikowany model objętości skończonych sieci. NASA wskazywała metodę charakterystyk jako klasyczne narzędzie udaru, ale zarazem jej ograniczenia przy złożonych obwodach; wybór modelu jest częścią walidacji, nie kwestią przyzwyczajenia.nasa-transient-models
Udar i czas zamknięcia
Dla nagłej zmiany średniej prędkości cieczy $\Delta u$ w jednorodnym, wypełnionym przewodzie liniowa relacja Joukowsky’ego daje
$$\Delta p=\rho a\Delta u.$$
Pełny skok tej wartości jest właściwym przybliżeniem szybkiego zamknięcia, gdy czas działania $t_c$ jest krótszy od czasu powrotu fali do granicy i z powrotem, w prostym układzie zbiornik–zawór $t_c\lt2L/a$. Przy wolnym zamknięciu $t_c>2L/a$ fala wielokrotnie wraca podczas ruchu zaworu, a maksimum zależy od prawa zamykania, strat i warunków brzegowych. Monografia USBR wyprowadza równoważną wysokość udaru $h_{max}=aV_0/g$ dla zamknięcia natychmiastowego i rozróżnia czas zamykania względem czasu przebiegu fali.usbr-joukowsky
Relacja ta nie obejmuje kawitacji, rozdzielenia kolumny, gazu uwięzionego, nieliniowej podatności ścian, dwufazowości, ruchu konstrukcji ani złożonych odbić. Po spadku ciśnienia do granicy parowania model jednofazowy może przewidzieć nieistniejące ciśnienia rozciągające. Wtedy trzeba modelować powstanie i zapadanie obszarów parowych, a kryterium obejmuje zarówno dodatni szczyt, jak i minimalne ciśnienie absolutne.
Akumulator ograniczy udar tylko wtedy, gdy jego admitancja w skali frontu fali jest dostatecznie duża. Duża objętość gazu za długim, wąskim portem może być statycznie imponująca, lecz dynamicznie odcięta przez $j\omega L_h$. Przewidywaną sygnaturą jest wówczas opóźniony przepływ do komory: pierwszy szczyt pozostaje wysoki, choć późniejsze oscylacje zanikają.
Dwie różne granice przepływu: ciecz i gaz
Kawitacja i wrzenie cieczy
Kawitacja jest powstawaniem i dynamiką jam parowych w cieczy wskutek lokalnego obniżenia ciśnienia; wrzenie opisuje przemianę fazową powiązaną z temperaturą i warunkami zarodkowania. W obu ocenach porównuje się lokalne ciśnienie absolutne z ciśnieniem pary nasyconej $p_v(T)$ właściwego płynu. Odczyt manometryczny przy zbiorniku nie jest takim porównaniem.
Dla punktu krytycznego $x$ rekonstruuje się przykładowo
$$p_{abs,x}=p_{abs,s}+\rho g(z_s-z_x)-\Delta p_{s\rightarrow x} -L_{eq}\frac{dQ}{dt},$$
z zachowaniem znaków konkretnej drogi. $\Delta p_{s\rightarrow x}$ obejmuje straty przewodów, filtrów, zaworów i wejścia do pompy, a składnik bezwładnościowy ma znaczenie przejściowe. Temperatura $T_x$ musi pochodzić z punktu reprezentatywnego dla gorącego obszaru, nie z najwygodniejszego czujnika.
Dla pompy dostępna nadwyżka antykawitacyjna może być zapisana jako
$$NPSH_A=\frac{p_{t,in,abs}-p_v(T_{in})}{\rho g},$$
gdzie $p_{t,in,abs}$ jest absolutnym ciśnieniem całkowitym na zdefiniowanej płaszczyźnie wlotowej. $NPSH_A$ wynika z instalacji. $NPSH_R$ jest wynikiem umownej próby danej pompy przy określonym strumieniu, prędkości i medium; nie jest termodynamiczną stałą ani gwarancją braku wszystkich postaci kawitacji. Warunek decyzyjny ma postać $NPSH_A-U_{NPSH}>NPSH_R+M$, gdzie $M$ jest uzasadnionym marginesem projektu. Eksperymentalna praca NASA o przewidywaniu kawitacyjnych charakterystyk pomp wiąże NPSH z mierzoną degradacją osiągów i warunkami próby, co uzasadnia oddzielenie wielkości dostępnej od wymaganej.nasa-pump-cavitation
Samo chwilowe $p_{abs}\lt p_v$ nie zawsze wyznacza dokładnie chwilę obserwowalnej kawitacji: znaczenie mają jądra, czas przebywania, napięcie powierzchniowe i metastabilność. Analiza Roberta B. Jacobsa wskazuje wprost, że przewidywanie mierzalnych objawów kawitacji wymaga dodatkowych informacji o metastabilności, zarodkowaniu i dynamice fazy parowej.nist-cavitation
Krytyczny przepływ gazu
Gaz w małym porcie podlega innej granicy. Dla idealnego, stacjonarnego, izentropowego przepływu przez dyszę o najmniejszym polu $A_t$ maksimum strumienia masy występuje przy liczbie Macha równej jeden. Dla dostatecznie niskiego ciśnienia za przewężeniem strumień nie rośnie już proporcjonalnie do dalszego obniżania tego ciśnienia. Dla gazu doskonałego krytyczny stosunek ciśnień wynosi
$$\frac{p_d}{p_0}\le \left(\frac{2}{\gamma+1}\right)^{\gamma/(\gamma-1)},$$
a idealny strumień graniczny, po uwzględnieniu współczynnika wypływu $C_d$, można zapisać
$$\dot m_{max}=C_dA_t p_0\sqrt{\frac{\gamma}{R_sT_0}} \left(\frac{2}{\gamma+1}\right)^{(\gamma+1)/[2(\gamma-1)]}.$$
$p_0,T_0$ są wielkościami całkowitymi przed dyszą. Równanie wymaga kalibracji $C_d$ i nie opisuje krótkiego, lepkiego, ogrzewanego ani dwufazowego portu bez poprawek. NASA Glenn wyprowadza maksimum z równania ciągłości i relacji izentropowych oraz lokalizuje dławienie przy $M=1$ w gardzieli.nasa-gas-choking Krytyczny przepływ gazu nie jest synonimem kawitacji cieczy; różnią się fazą, kryterium i sygnaturą.
Separator jako część modelu, nie etykieta konstrukcji
Idealna, bezmasowa i bezoporowa granica faz jest punktem startowym. Rzeczywisty separator ma masę, sztywność, powierzchnię czynną zależną od położenia, tarcie lub histerezę materiałową i skończony zakres ruchu. W osi $x$ jego bilans sił można zapisać schematycznie:
$$[p_l-p_g]A_{eff}(x)=m_{eff}\ddot x+F_s(x,T)+F_f(\dot x,x,\mathcal H,T) +m_{eff}a_x+F_{contact}(x)+F_{side}.$$
$\mathcal H$ oznacza historię ruchu. Równanie nie jest katalogiem dodatków, lecz mapą sygnatur. Masa daje przesunięcie fazowe i rezonans mechaniczny. Sztywność zmienia krzywą $P$–$V$. Tarcie daje martwą strefę, skok po ruszeniu i pole pętli. Kontakt z ogranicznikiem daje nagły spadek podatności. Obciążenie poprzeczne zmienia siłę normalną, a więc tarcie. Jeśli model nie przewiduje kanału pomiarowego zdolnego odróżnić te zjawiska, urządzenie nie jest obserwowalne w wymaganym sensie.
Pęcherz
Pęcherz jest zamkniętą elastyczną przegrodą, zwykle zawierającą gaz. Ma małą masę ruchomą i może reagować szybko, ale jego objętość nie jest jednoznaczną funkcją jednego przemieszczenia osiowego. Fałdy przemieszczają się i stykają, materiał pracuje dwuosiowo, a lokalne promienie oraz odkształcenia zależą od historii napełnienia. Przy skrajnym skoku pęcherz może zakryć port, zostać do niego wciągnięty, ulec ściśnięciu albo ścieraniu o korpus.
Przewidywane sygnatury są konkretne: częściowe zasłonięcie portu zwiększa $R_h$ i opóźnia przepływ bez koniecznej utraty wstępnego ciśnienia gazu; zmiana układu fałd powoduje niepowtarzalność krzywej $P$–$V$; permeacja lub mikronieszczelność daje powolny trend masy/ciśnienia gazu; pęknięcie miesza fazy i może spowodować gaz w cieczy. Badania mechaniczne elastomeru w jednej temperaturze nie wystarczają. Trzeba znać konkretną mieszankę, medium, czas ekspozycji, temperaturę, zakres odkształceń i wpływ starzenia. Publiczny program NASA dotyczący konkretnego separatora wykazał narastanie gazu w pętli wskutek dyfuzji przez membranę polimerową i zweryfikował zmianę na spawany mieszek metalowy próbą szczelności funkcjonalnej; to przykład metody, nie uniwersalny ranking materiałów.nasa-diaphragm-permeation
Membrana
Membrana jest przegrodą zamocowaną na obwodzie. Usuwa uszczelnienie ślizgowe i może mieć dobrze kontrolowaną granicę faz, lecz zwykle oferuje mniejszy skok niż tłok albo pęcherz. Powierzchnia czynna zmienia się z ugięciem; naprężenia koncentrują się przy zamocowaniu, przetłoczeniach i strefach kontaktu. Membrana profilowana może przechodzić między konfiguracjami z lokalnym przeskokiem, a kontakt z podporą zmienia jej sztywność.
Sygnaturą narastającego uszkodzenia może być wzrost pola histerezy, przesunięcie ciśnienia ruszenia, zmiana harmonicznych albo przeciek między stronami jeszcze przed utratą integralności korpusu. Skrajny skok trzeba traktować jako osobny przypadek obciążenia. Polimerowe pętle histerezy opisują energię rozproszoną w cyklu, a samonagrzewanie i częstotliwość zmieniają odpowiedź zmęczeniową; lokalny podręcznik laboratoryjny wspiera te metody badań, ale nie jest źródłem danych kompatybilności ani kwalifikacji gotowego separatora.local-polymers
Tłok
Tłok daje duży, geometrycznie czytelny skok i umożliwia bezpośredni pomiar położenia. W zamian wprowadza uszczelnienie ślizgowe, masę i prowadzenie. Równowaga statyczna obejmuje tarcie spoczynkowe; małe zmiany ciśnienia mogą nie poruszyć tłoka, aż siła przekroczy próg. Po ruszeniu tarcie może zmaleć, co tworzy skok i histerezę. Niewspółosiowość, boczne obciążenie od przewodu, cząstka w szczelinie lub zużycie prowadzenia zwiększają siłę normalną i tarcie.
Przeciek przez uszczelnienie nie musi natychmiast obniżyć ciśnienia powłoki. W krzywej długookresowej objawia się dryfem położenia przy stałym bilansie cieczy; w próbie dynamicznej — zmianą fazy i oporu; w oględzinach medium — migracją faz. Orientacja i przyspieszenie dodają $m_{eff}a_x$ do równania, więc próba pozioma nie dowodzi ruchu w innej osi.
Mieszek
Mieszek metalowy oddziela fazy bez uszczelnienia ślizgowego. Jego zaletą jest hermetyczna bariera, a kosztami są sztywność osiowa, naprężenia cykliczne fałd, spoiny, możliwość lokalnego kontaktu, przechylenia i utraty stateczności. Pole czynne, długość i sztywność nie są stałe w całym skoku. Mieszek może być szczelny i utrzymywać właściwe ciśnienie gazu, a zarazem nie wykonać wymaganego przemieszczenia wskutek zakleszczenia.
Model musi łączyć ciśnienie z mechaniką prowadzenia. Zwiększenie różnicy ciśnień nie zawsze „przełamuje” tarcie: jeśli przechylenie zwiększa siłę normalną szybciej niż siłę przywracającą, układ staje się samowzmacniający. Ten mechanizm zostanie prześledzony w studium MSL. Dla trwałości nie wystarcza nominalne naprężenie błonowe; istotne są lokalne odkształcenia fałd, historia formowania, spoiny, karby, środowisko i pełne widmo cykli.
Wolny interfejs gaz–ciecz
Wolny interfejs nie jest „pęcherzem bez gumy”. Nie ma materialnej przegrody. Gaz styka się bezpośrednio z cieczą, więc do modelu wchodzą rozpuszczanie, dyfuzja, wynoszenie pęcherzy, falowanie, orientacja, przyspieszenia i geometria poboru. W warunkach grawitacyjnych pionowy zbiornik może utrzymywać rozdział dzięki wyporowi; przy zmianie orientacji albo przyspieszenia położenie faz może się odwrócić. W małej skali znaczenia nabierają zwilżanie i siły kapilarne.
Spadek objętości widocznego gazu może wynikać zarówno ze wzrostu ciśnienia, jak i z rozpuszczania; sam pomiar $p$ nie rozdziela tych hipotez. Obecność gazu w przewodzie za poborem, zmiana gęstości pozornej, niestabilny przepływ i wolna odbudowa stanu po postoju są sygnaturami fazowymi, których nie ma w idealnym modelu z pęcherzem. Architektury nie wolno mieszać z pęcherzem ani przenosić między nimi tych samych trybów awarii. Badania NASA zachowania płynów w obniżonej grawitacji wiążą kształt wolnej powierzchni, reorientację i pobór bez gazu z napięciem powierzchniowym, zwilżaniem, geometrią i przyspieszeniem.nasa-free-interface
Mapa mechanizmu do obserwacji
| Mechanizm | Przewidywana sygnatura | Pomiar rozstrzygający | Możliwa pomyłka |
|---|---|---|---|
| utrata masy gazu | przesunięta cała krzywa $P$–$V$, większy skok dla tego samego ciśnienia | $p_g,T_g,V_g$ po ustaleniu cieplnym | chwilowe chłodzenie gazu |
| przywieranie/tarcie | martwa strefa, różne ciśnienia ruszenia w obu kierunkach, skok po ruszeniu | zsynchronizowane $p_l,p_g,x,Q$ | opóźnienie czujnika lub dławik |
| wzrost inertancji portu | rosnąca różnica faz i słaba odpowiedź przy wysokiej częstotliwości, poprawna próba statyczna | ciśnienie po obu stronach portu i przepływ gałęzi | spadek podatności gazu |
| ogranicznik skoku | gwałtowny wzrost sztywności i zanik przyjmowanej objętości | $x$ lub niezależny bilans objętości | zawór odcinający albo zablokowany port |
| przeciek separatora | dryf masy faz, gaz w cieczy lub ciecz po stronie gazu | analiza medium, masa, długotrwały bilans | rozpuszczanie gazu przy wolnym interfejsie |
| zmęczenie/przerwanie | zmiana histerezy, następnie przeciek; niekiedy brak wcześniejszej sygnatury makro | NDE, próba szczelności, trend $P$–$V$ | zmiana temperatury lub własności medium |
| przechylenie/zakleszczenie | powtarzalny przedwczesny koniec skoku, rosnąca siła normalna | położenie, objętość, obrazowanie, próba obciążeń bocznych | niewłaściwe napełnienie |
Tabela nie zastępuje drzewa usterek. Jej funkcją jest wymuszenie co najmniej jednej konkurencyjnej hipotezy. Jeśli anomalię „wyjaśnia” jedyny dostępny czujnik, należy podejrzewać brak obserwowalności.
Od modelu do próby i decyzji
Próba ma obalić błędny model, nie tylko wytworzyć ładny wykres. W metrologii przepływów pulsujących wskazanie zależy od dynamiki przyrządu, profilu strugi, nieliniowości przetwarzania i miejsca instalacji. Lokalna monografia Turkowskiego omawia identyfikację pulsacji na s. 177–181, pomiary i skutki dynamiczne na s. 182–203 oraz praktyczny rachunek niepewności na s. 253–262; te części są podstawą projektu toru pomiarowego, nie źródłem dopuszczalnych parametrów konkretnej instalacji.local-flow-metrology
Minimalny zapis stanu i sygnatury
Dla próby quasi-statycznej zapis obejmuje co najmniej absolutne ciśnienie po obu stronach separatora, temperaturę gazu, cieczy i korpusu, objętość lub masę wtłoczonej cieczy, położenie separatora, czas ustalania i orientację. Dla próby dynamicznej dochodzą zsynchronizowane przebiegi ciśnienia i przepływu, częstotliwość/prędkość źródła, faza wyzwalania i stan średni. Jeden czujnik ciśnienia nie rozdzieli podatności komory od inertancji rurki impulsowej.
Przed próbą zapisuje się przewidywanie liczbowe z tolerancją: np. zakres $p(T)$, częstotliwość minimum $|H|$, znak fazy poniżej i powyżej rezonansu, maksymalny skok oraz energię rozpraszaną. Po próbie porównuje się reszty $r=y_{meas}-y_{model}$, a nie tylko wybrane punkty. Systematyczny kształt reszty jest informacją o brakującym mechanizmie: przesunięcie osi ciśnienia sugeruje kalibrację lub masę gazu, pętla — tarcie/termikę, błąd rosnący z częstotliwością — dynamikę sensora albo model rozłożony.
Pasmo, próbkowanie i synchronizacja
Pasmo czujnika musi obejmować nie tylko pasmo kryterium, lecz również harmoniczne potrzebne do rozpoznania nieliniowości i antyaliasingu. Podana przez producenta „częstotliwość własna” nie jest automatycznie płaskim pasmem pomiarowym. Czujnik wraz z przyłączem, objętością martwą i rurką tworzy własną gałąź $RLC$. Kalibracja statyczna nie dowodzi poprawnej amplitudy i fazy dynamicznej.
Częstotliwość próbkowania większa od dwukrotności najwyższej częstotliwości jest warunkiem Nyquista, a nie kompletnym projektem pomiaru. Potrzebne są filtr antyaliasingowy, wspólny zegar, znany błąd synchronizacji i wystarczająco długi rekord. Błąd czasu $\Delta t$ daje błąd fazy $\Delta\phi=360° f\Delta t$; ten sam zegar, który jest nieistotny przy 5 Hz, może zafałszować moc zespoloną przy wyższych częstotliwościach.
Niepewność jako część kryterium
Wielkość mierzona $y=f(x_1,\ldots,x_m)$ wymaga budżetu obejmującego wzorcowanie, rozdzielczość, dryf, temperaturę, montaż, synchronizację, powtarzalność oraz niepewność modelu użytego do przeliczenia. Dla liniaryzacji:
$$u_c^2(y)=\sum_i\left(\frac{\partial f}{\partial x_i}\right)^2u^2(x_i) +2\sum_{i\lt j}\frac{\partial f}{\partial x_i}\frac{\partial f}{\partial x_j}u(x_i,x_j).$$
Nie wolno dodawać korelacji dwukrotnie ani traktować wielkości wyprowadzonej z tych samych czujników jako niezależnej. Dla silnej nieliniowości, ogranicznika skoku lub rozkładu asymetrycznego właściwsza bywa propagacja Monte Carlo. NIST TN 1297 rozdziela składniki typu A i B, łączy standardowe niepewności z uwzględnieniem kowariancji i definiuje niepewność rozszerzoną $U=ku_c$.nist-uncertainty
Kryterium powinno działać na przedziale, nie na wartości zaokrąglonej. Dla górnej granicy $y_{max}$ bezpieczna decyzja wymaga $y+U\le y_{max}$; dla dolnej $y-U\ge y_{min}$. Wynik przecinający granicę jest nierozstrzygający lub niezgodny zgodnie z ustaloną regułą decyzyjną — nie „prawie spełniony”. Niepewność modelu i niepewność pomiaru zapisuje się osobno, aby kalibracja czujnika nie maskowała błędu fizyki.
Regresja
Próba regresyjna powtarza tę samą wielkość, stan i regułę decyzji po zdarzeniu mogącym zmienić urządzenie: cyklach ciśnienia, temperaturze, drganiach, naprawie, czyszczeniu, zmianie dostawcy lub integracji z przewodem. „Urządzenie nadal działa” jest słabym wynikiem. Silny wynik brzmi: krzywa $P$–$V$, częstotliwość $f_0$, $|H|$, faza, przeciek i próg ruszenia pozostają w zdefiniowanych korytarzach wraz z niepewnością.
Metoda badawcza powinna zachować łańcuch od parametrów przez model i stanowisko do walidacji, a modyfikacja otwiera ponownie te ogniwa, których dotyczy. Taką strukturę — choć dla innej klasy maszyn — konsekwentnie stosuje Marian Gieras, łącząc geometrię, obliczenia, stanowisko i porównanie wyników; jest to użyteczny wzorzec redakcyjny i inżynierski, nie źródło parametrów akumulatora.local-gieras-method
Integralność i cykl życia wynikają z przewidywanego uszkodzenia
Powłoka nie jest osobnym „tematem bezpieczeństwa” doczepionym po analizie funkcji. Gaz magazynuje energię, ciśnienie obciąża korpus, port przenosi reakcje przewodów, separator cyklicznie obciąża spoiny i prowadzenie, a zanieczyszczenie zmienia tarcie lub blokuje mały przekrój. Każdy proces zapewnienia ma zatem zamykać konkretny człon modelu.
Granica ciśnieniowa i nadzór odporności na pękanie
Analiza naprężeń cienkiej lub grubej powłoki daje pole nominalne, lecz rzeczywista granica zawiera króćce, spoiny, gwinty, zmiany grubości i obciążenia montażowe. Modele powłok, cylindrów grubościennych, dynamiki i zmęczenia w rozdziałach 11–13 i 15 podręcznika Misiaka pomagają wyznaczyć idealizację i krytyczne przekroje; autor sam nie zastępuje aktualnego kodu ciśnieniowego, danych materiałowych ani mechaniki pękania.local-strength
Nadzór odporności na pękanie (fracture control) rozpoczyna się od skutku hipotetycznej wady: jeśli pęknięcie może spowodować katastrofalną utratę granicy, określa się krytyczne części, widmo obciążeń, początkowy rozmiar wady wynikający ze zdolności badań nieniszczących (NDE), wzrost pęknięcia i pozostałe życie. NASA-STD-5019A z Change 4 ustanawia taki program dla sprzętu lotnego z udziałem człowieka; nie jest uniwersalnym kodem każdej instalacji, ale poprawnie wiąże projekt, analizę, inspekcję, próbę, dokumentację i śledzenie części krytycznych.
Wynika stąd wymaganie produkcyjne dla metalowej części krytycznej: jeżeli ocena tolerancji uszkodzeń analizuje wadę o określonym rozmiarze i kształcie, inspekcja NDE sprzętu lotnego ma przesiewać wady co najmniej tego rozmiaru i kształtu. NASA-STD-5019A odsyła przy tym do NASA-STD-5009 i objaśnia, że zdolność wykrywania przyjmowana w analizie opiera się na prawdopodobieństwie detekcji 90% z ufnością 95% lub lepszym.nasa-fracture „Sprawdzono penetrantem” bez wykazanej zdolności detekcji, dostępu i kryterium akceptacji nie zamyka modelu. Spoiny mieszka wymagają innej demonstracji niż cylindryczny korpus; złożona geometria lub różne rodzaje wad mogą wymagać dodatkowych metod.
Czystość, materiały i wytwarzanie
Cząstka w szczelinie prowadzącej może zwiększyć $F_{side}$, zadra może przeciąć pęcherz, pozostałość procesu może skleić powierzchnie, a film chemiczny może zmienić tarcie. Czystość ma więc poziom wynikający z najmniejszego krytycznego prześwitu, wrażliwości materiału i dopuszczalnej zmiany oporu, a nie arbitralne hasło „czyste”. NASA-STD-6016C z Change 1 wymaga planu kontroli zanieczyszczeń i programu zapobiegania obcym przedmiotom i zanieczyszczeniom (FOD) w projektowaniu, wytwarzaniu, montażu, naprawie, próbach i operacjach sprzętu kosmicznego; zastosowanie do innego projektu wymaga dostosowania zakresu.nasa-cleanliness
Materiał separatora kwalifikuje się jako konkretną konfigurację: receptura, partia, proces, grubość, orientacja, złącze, medium, temperatura, dawka środowiska i historia cykli. Próba rozciągania próbki nie dowodzi odporności fałdu pęcherza ani krawędzi membrany. Wytwarzanie utrwala naprężenia resztkowe, chropowatość, odchyłkę współosiowości i stan powierzchni — parametry w równaniu sił, a nie tylko pozycje na rysunku.
Zakres norm i obowiązek dostosowania
ECSS-E-ST-32-02C Rev.2, datowany 15 października 2025 r., ustanawia minimalne wymagania projektowania i weryfikacji sprzętu ciśnieniowego. Sam dokument wskazuje jednak obszary nieobjęte w pełni: urządzenia wypierające z pęcherzami lub membranami, wymagania funkcjonalne takie jak szybkie wypieranie i czystość oraz elementy obciążone znacząco inaczej niż samym ciśnieniem, w tym mieszki. Te wyłączenia oznaczają potrzebę dodatkowych wymagań i analiz, nie zwolnienie z dowodu.ecss-scope
NASA-STD-8719.17D należy cytować jako zatwierdzony 26 kwietnia 2022 r. — taka data widnieje w normatywnym PDF; karta katalogowa może prezentować inną datę. Standard dotyczy naziemnych naczyń i systemów ciśnieniowych oraz naziemnego przetwarzania lotnych systemów, nie funkcjonalnej kwalifikacji dowolnego akumulatora pokładowego. Jego wymagania dotyczące certyfikacji, kontroli, dokumentacji i ponownej certyfikacji są właściwym przykładem nadzoru energii ciśnieniowej tylko w tym zakresie.nasa-ground-pvs
Kwalifikacja ma wykazać margines projektu w reprezentatywnym środowisku; odbiór ma wykryć wady wykonania bez niepotrzebnego zużycia życia; nadzór ma wykryć dryf między tymi zdarzeniami. Po naprawie, zmianie materiału, procesu, prowadzenia albo objętości martwej wykonuje się ocenę wpływu. Jeśli zmiana dotyka $C_h$, $L_h$, $R_h$, sztywności, tarcia, szczelności lub granicy ciśnieniowej, odpowiednie próby regresyjne muszą wrócić.
Studium A: syntetyczny kompensator pompowanej pętli cieplnej
Wszystkie liczby w tym studium są syntetycznymi, bezpiecznymi danymi dydaktycznymi, nie opisują istniejącego wyrobu i nie są powiązane z uzbrojeniem. Medium nazwano „cieczą dydaktyczną”, a gaz „gazem obojętnym typu N₂”; podane własności są założeniami rachunku, nie kartą materiałową.
Wymaganie, granica i stan początkowy
Rozpatrujemy zamkniętą pompowaną pętlę laboratoryjną od $263{,}15$ do $323{,}15 \mathrm K$ (od −10 do 50°C). Kompensator ma utrzymać absolutne ciśnienie węzła w przedziale wyznaczonym od dołu przez ostrzejsze z dwóch kryteriów — antykawitacyjne na wlocie pompy i przeciw wrzeniu w najgorętszym punkcie — oraz od góry przez syntetyczny limit funkcjonalny $235 \mathrm{kPa(abs)}$. Limit nie jest dopuszczalnym ciśnieniem żadnego realnego zbiornika.
Kryteria dolne są:
$$p_{node}\ge p_v(T_{in})+\Delta p_{node\to in} +\rho g(NPSH_R+M_N),$$
$$p_{node}\ge p_v(T_{hot})+\Delta p_{node\to hot}+M_B.$$
Syntetycznie przyjęto $NPSH_R=4{,}0 \mathrm m$, margines $M_N=2{,}0 \mathrm m$, stratę do wlotu $18 \mathrm{kPa}$, stratę do gorącego punktu $24 \mathrm{kPa}$ i margines wrzenia $M_B=25 \mathrm{kPa}$. W analizie rzeczywistej $NPSH_R$ pochodziłoby z kwalifikowanej charakterystyki pompy, a $p_v(T)$ z danych konkretnego płynu.
Stan odniesienia w $T_0=293{,}15 \mathrm K$:
| Wielkość | Syntetyczna wartość | Znaczenie |
|---|---|---|
| całkowita objętość cieczy $V_{L0}$ | $4{,}300 \mathrm L$ | cała masa cieczy przeliczona na objętość w $T_0$ |
| pojemność zewnętrznej części pętli $V_{e0}$ | $4{,}000 \mathrm L$ | bez komory kompensatora |
| objętość komory $V_{c0}$ | $1{,}000 \mathrm L$ | gaz + ciecz w kompensatorze |
| ciecz w kompensatorze $V_{l0}$ | $0{,}300 \mathrm L$ | $V_{L0}-V_{e0}$ |
| gaz w kompensatorze $V_{g0}$ | $0{,}700 \mathrm L$ | $V_{c0}-V_{l0}$ |
| ciśnienie gazu $p_{g0}$ | $180{,}0 \mathrm{kPa(abs)}$ | po ustaleniu cieplnym |
| różnica separatora $p_l-p_g$ | $3{,}0 \mathrm{kPa}$ | syntetyczny model quasi-statyczny |
| ciśnienie otoczenia do ilustracji | $101{,}3 \mathrm{kPa(abs)}$ | tylko do przeliczenia wskazania manometrycznego |
Objętości są jawnie rozłączne: $V_{L0}$ dotyczy całej cieczy, $V_{e0}$ przestrzeni poza kompensatorem, $V_{l0}$ cieczy w nim, a $V_{g0}$ gazu. Na tej podstawie operator może odtworzyć każdy wiersz, zamiast ufać nieopisanej „pojemności układu”.
Do audytu położenia przyjmujemy dodatkowo jawną, syntetyczną kalibrację geometryczną. Oś $x$ jest zgodna z konwencją ogólną: rośnie od dolnego ogranicznika w kierunku zwiększania $V_l$. Stała efektywna powierzchnia $A_{eff}=50{,}0 \mathrm{cm^2}$ odpowiada przeliczeniu $5{,}00 \mathrm{mL/mm}$. Przy dolnym ograniczniku $x_{min}=0$ objętość cieczy wynosi $V_{l,min}=0{,}2200 \mathrm L$, a górny ogranicznik leży w $x_{max}=30{,}00 \mathrm{mm}$, czyli przy $V_{l,max}=0{,}3700 \mathrm L$. Dla swobodnego ruchu między ogranicznikami:
$$x=\frac{V_l-V_{l,min}}{A_{eff}} =\frac{V_l-0{,}2200 \mathrm L}{5{,}00 \mathrm{mL/mm}}.$$
To równoważne przeliczenie objętości na położenie służy wyłącznie temu syntetycznemu przykładowi; nie opisuje geometrii istniejącego wyrobu.
Przyjmujemy syntetyczne objętościowe współczynniki rozszerzalności: cieczy $\beta_l=4{,}00\cdot10^{-4} \mathrm{K^{-1}}$, części zewnętrznej i komory $\beta_e=\beta_c=5{,}00\cdot10^{-5} \mathrm{K^{-1}}$. Dla $\Delta T=T-T_0$:
$$V_L(T)=V_{L0}(1+\beta_l\Delta T),$$
$$V_e(T)=V_{e0}(1+\beta_e\Delta T), \qquad V_c(T)=V_{c0}(1+\beta_c\Delta T),$$
$$V_l(T)=V_L(T)-V_e(T), \qquad V_g(T)=V_c(T)-V_l(T).$$
Proces cieplny jest powolny, zatem w modelu przyjęto gaz idealny pozostający w temperaturze $T$ i stałą masę:
$$p_g(T)=p_{g0}\frac{T}{T_0}\frac{V_{g0}}{V_g(T)}, \qquad p_l(T)=p_g(T)+3{,}0 \mathrm{kPa}.$$
Dla syntetycznego $R_s=296{,}8 \mathrm{J/(kg\,K)}$ masa gazu w stanie odniesienia wynosi
$$m_g=\frac{180000\cdot0{,}000700}{296{,}8\cdot293{,}15} =0{,}001448 \mathrm{kg}=1{,}448 \mathrm g.$$
Ta kontrola $P$–$V$–$T$ jest niezależna od bilansu objętości cieczy: oba rachunki spotykają się dopiero w $V_g(T)$.
Przewidywana sygnatura i rachunek nominalny
| $T$ | $V_L$ | $V_e$ | $V_c$ | $V_l$ | $V_g$ | $p_g$ abs. | $p_l$ abs. | $p_l$ manom. przy $p_{amb}=101{,}3$ kPa |
|---|---|---|---|---|---|---|---|---|
| $263{,}15 \mathrm K$ | 4,2484 L | 3,9940 L | 0,9985 L | 0,2544 L | 0,7441 L | 152,00 kPa | 155,00 kPa | 53,70 kPa |
| $293{,}15 \mathrm K$ | 4,3000 L | 4,0000 L | 1,0000 L | 0,3000 L | 0,7000 L | 180,00 kPa | 183,00 kPa | 81,70 kPa |
| $323{,}15 \mathrm K$ | 4,3516 L | 4,0060 L | 1,0015 L | 0,3456 L | 0,6559 L | 211,76 kPa | 214,76 kPa | 113,46 kPa |
Syntetyczny budżet bezpośredniego toru położenia obejmuje standardowe niepewności: wzorcowania $0{,}30 \mathrm{mm}$, powtarzalności $0{,}25 \mathrm{mm}$ oraz wpływu temperatury i orientacji $0{,}20 \mathrm{mm}$. Daje to $u_c(x)=0{,}439 \mathrm{mm}$ oraz, dla $k=2$, $U(x)=0{,}878 \mathrm{mm}$, konserwatywnie zaokrąglone do $0{,}90 \mathrm{mm}$. Położenie każdego ogranicznika ma osobną rozszerzoną niepewność $U_{stop}=0{,}20 \mathrm{mm}$. Reguła przedziałowa wymaga więc zawarcia całego przedziału $[x-U(x),x+U(x)]$ w konserwatywnym zakresie $[0{,}20;29{,}80] \mathrm{mm}$.
| Stan | $V_l$ | $x$ | nominalna rezerwa do dolnego / górnego ogranicznika | przedział $x\pm U(x)$ | konserwatywna rezerwa dolna / górna |
|---|---|---|---|---|---|
| zimny | 0,2544 L | 6,88 mm | 6,88 / 23,12 mm | [5,98; 7,78] mm | 5,78 / 22,02 mm |
| odniesienia | 0,3000 L | 16,00 mm | 16,00 / 14,00 mm | [15,10; 16,90] mm | 14,90 / 12,90 mm |
| gorący | 0,3456 L | 25,12 mm | 25,12 / 4,88 mm | [24,22; 26,02] mm | 24,02 / 3,78 mm |
Wszystkie konserwatywne rezerwy są dodatnie, zatem trzy stany nominalne spełniają kryterium położenia.
Kontrola przykładowego gorącego wiersza:
$$V_g=1{,}0015-[4{,}3516-4{,}0060]=0{,}6559 \mathrm L,$$
$$p_g=180{,}0\frac{323{,}15}{293{,}15}\frac{0{,}7000}{0{,}6559} =211{,}76 \mathrm{kPa(abs)}.$$
Przewidywaną sygnaturą zdrowego kompensatora jest monotoniczny wzrost $V_l$ i $p_l$ z temperaturą, bez nieciągłości, z pozostałym gazem po gorącej stronie i pozostałą cieczą po zimnej. Nieciągłość położenia przy gładkim $p(T)$ wskazywałaby na pomiar położenia; nieciągłość $p$ i $x$ przy zmianie kierunku — na tarcie lub przeskok separatora; ciśnienie rosnące szybciej niż bilans — na nadmiar cieczy, mniejszą objętość gazu albo błąd termiczny.
Do sprawdzenia granicy antykawitacyjnej przyjęto syntetyczne $p_v$: $1 \mathrm{kPa}$ na zimno, $3 \mathrm{kPa}$ w stanie odniesienia i $15 \mathrm{kPa}$ na gorąco, a gęstości odpowiednio 1000, 995 i $985 \mathrm{kg/m^3}$. Nominalne $NPSH_A$ po odjęciu $18 \mathrm{kPa}$ straty wynosi:
| Stan | $p_{in,abs}=p_l-18$ kPa | $p_v$ | $NPSH_A$ | Wymagane $NPSH_R+M_N$ |
|---|---|---|---|---|
| zimny | 137,00 kPa | 1 kPa | 13,87 m | 6,00 m |
| odniesienia | 165,00 kPa | 3 kPa | 16,60 m | 6,00 m |
| gorący | 196,76 kPa | 15 kPa | 18,82 m | 6,00 m |
W gorącym punkcie $p_{hot}=214{,}76-24=190{,}76 \mathrm{kPa(abs)}$, czyli nominalna różnica ponad $p_v$ wynosi $175{,}76 \mathrm{kPa}$. Dwie dolne granice w gorącym stanie są równe około $91 \mathrm{kPa(abs)}$ z kryterium pompy oraz $64 \mathrm{kPa(abs)}$ z kryterium wrzenia. O sterowaniu dolnym progiem decyduje tu pompa. Nie wolno jednak z tego uczynić reguły: inna temperatura gorącego punktu albo strata przepływu może odwrócić kolejność.
Pomiar, niepewność, przypadek graniczny i decyzja
Kanały próby to: dwa absolutne przetworniki ciśnienia po obu stronach separatora; ciśnienie bezpośrednio przed pompą; temperatura gazu, cieczy przy pompie, gorącego punktu i korpusu; masa napełnienia; objętość dozowana wzorcem; położenie separatora; strumień oraz zapis stanu zaworów. Przed odczytem quasi-statycznym wymaga się ustalonego gradientu temperatur, a nie tylko upływu arbitralnego czasu. Bilans masy i odczyt położenia są dwiema niezależnymi obserwacjami objętości.
Syntetyczny budżet standardowy przyjmuje: $u(p_{g0})=0{,}5 \mathrm{kPa}$, $u(V_{g0})=2 \mathrm{mL}$, $u[V_g(T)]=4 \mathrm{mL}$ w skrajnych temperaturach, $u(T)=0{,}25 \mathrm K$, $u(p_l-p_g)=0{,}5 \mathrm{kPa}$ i składnik modelu cieplnego $u_{model}=1{,}5 \mathrm{kPa}$. Dla
$$p_l=p_{g0}\frac{T}{T_0}\frac{V_{g0}}{V_g}+\Delta p_s$$
użyto pochodnych analitycznych i sumy kwadratowej. Otrzymane niepewności rozszerzone dla $k=2$, po konserwatywnym zaokrągleniu w górę, wynoszą $U=4 \mathrm{kPa}$ na zimno i w stanie odniesienia oraz $U=5 \mathrm{kPa}$ na gorąco. Osobno propagowane niepewności strat i $p_v$ nie zbliżają dolnego progu do przedziału ciśnienia kompensatora.
W stanie odniesienia oba wystąpienia $V_{g0}$ w ilorazie $V_{g0}/V_g(T_0)=V_{g0}/V_{g0}$ oznaczają tę samą wielkość i są doskonale skorelowane, a nie dwa niezależne pomiary. Iloraz skraca się algebraicznie do jedności, więc jego pochodne łączne względem $V_{g0}$ znoszą się i składnika $u(V_{g0})$ nie dodaje się dwukrotnie. W przyjętym modelu $p_l(T_0)=p_{g0}+\Delta p_s$; podane $U=4 \mathrm{kPa}$ wynika z pozostałych składników budżetu i konserwatywnego zaokrąglenia.
Reguła decyzji brzmi: cały przedział $[p_l-U,p_l+U]$ ma pozostać powyżej obu dolnych granic i nie przekraczać $235 \mathrm{kPa(abs)}$, a przedział położenia nie może przecinać żadnego ogranicznika. Dla nominalnego gorącego stanu górny kraniec wynosi $214{,}76+5=219{,}76 \mathrm{kPa}$ — wynik zgodny.
Przypadek graniczny celowo żąda dodania syntetycznego nadmiaru $50 \mathrm{mL}$ cieczy, bez zmiany komory. Swobodny model objętości wymagałby na gorąco $V_l=0{,}3956 \mathrm L$ oraz $x_{free}=35{,}12 \mathrm{mm}$, czyli przekroczenia górnego ogranicznika o $5{,}12 \mathrm{mm}$, równoważnego $25{,}6 \mathrm{mL}$. Także cały kontrfaktyczny przedział żądanego położenia $[34{,}22;36{,}02] \mathrm{mm}$ leży powyżej konserwatywnej granicy $29{,}80 \mathrm{mm}$. Kontakt nastąpiłby już po dodaniu $0{,}3700-0{,}3456=0{,}0244 \mathrm L$, a więc po $24{,}4 \mathrm{mL}$. Bezpośrednio przed kontaktem $V_g=1{,}0015-0{,}3700=0{,}6315 \mathrm L$, co w modelu bez kontaktu odpowiada około $222{,}94 \mathrm{kPa(abs)}$ po stronie cieczy.
Poniższy rachunek z $V_g=0{,}6059 \mathrm L$ jest zatem świadomie kontrfaktyczną ekstrapolacją: pokazuje, jaki wynik dałby model przy założeniu dalszego swobodnego ruchu separatora, którego rzeczywisty ogranicznik nie dopuszcza:
$$p_l=180{,}0\frac{323{,}15}{293{,}15}\frac{0{,}7000}{0{,}6059}+3{,}0 =232{,}24 \mathrm{kPa(abs)}.$$
Kontrfaktyczny przedział po uwzględnieniu zwiększonej niepewności $U=6 \mathrm{kPa}$ wynosi $[226{,}24;238{,}24] \mathrm{kPa}$ i przecina limit. Nie wolno użyć go jako prognozy rzeczywistego ciśnienia po kontakcie: wtedy $F_{contact}$ przejmuje dalsze obciążenie, znika założenie stałej różnicy $p_l-p_g=3 \mathrm{kPa}$, a ciśnienie zależy od niepodanych podatności cieczy, ścian i elementów pętli. Decyzja jest mimo to jednoznacznie niezgodna/odrzucić stan napełnienia już dlatego, że żądany przedział położenia przecina ogranicznik; kontrfaktyczny rachunek pokazuje dodatkowo, że nawet model bez ogranicznika nie zamykałby górnego kryterium ciśnienia. To przykład, dlaczego niepewności i bilansu nie dopisuje się po decyzji.
Po zaplanowanej liczbie syntetycznych cykli cieplnych regresja powtarza trzy punkty $P$–$V$–$T$, krzywą pośrednią przy ogrzewaniu i chłodzeniu, próbę szczelności oraz próg ruszenia. Kryterium wymaga pozostania w pierwotnych korytarzach $p(T)$ i $x(T)$, braku wzrostu histerezy ponad ustalony budżet oraz ponownego spełnienia reguły przedziałowej. Tak zbudowana próba odpowiada metodzie JPL: bilansuje temperaturę i objętość pętli, aby utrzymać ciśnienie pomiędzy granicami kawitacji, wrzenia i konstrukcji, zamiast dobierać kompensator tylko z różnicy skrajnych objętości.jpl-loop-balance
Studium B: dochodzenie zakleszczenia mieszka łazika MSL
To studium jest historycznym, cywilnym przypadkiem opisanym publicznie przez zespół JPL/NASA. Nie uzupełniamy raportu domysłami ani liczbami z innych konstrukcji. Jego wartością jest przebieg falsyfikacji i zmiana modelu przyczynowego.
Wymaganie i obserwacja
Mars Science Laboratory miał dwie mechanicznie pompowane pętle cieplne. Akumulator gaz–ciecz z metalowym mieszkiem musiał przyjąć zmianę objętości płynu w przewidzianym zakresie temperatur i zachować wymagany skok. Podczas końcowego badania zapasowego zintegrowanego zespołu pompy mieszek nie osiągnął wymaganej długości: ruch kończył się przedwcześnie, dostępna pojemność była mniejsza niż wymagana, a zjawisko powtarzało się. Obrazowanie rentgenowskie wykazało przechylenie zespołu końca mieszka/prowadzenia i zakleszczenie o obudowę.msl-observation
Łańcuch rozpoczął się prawidłowo od funkcji i mierzalnej sygnatury: wymaganiem był pełny użyteczny skok, a nie samo utrzymanie ciśnienia czy szczelności. Stanem było urządzenie po końcowej integracji. Pomiary obejmowały objętość/skok oraz obrazowanie wnętrza. Wynik „ciśnienie istnieje” nie mógł przesłonić utraty funkcjonalnej pojemności.
Hipoteza kleju i próba negatywna
Po rozcięciu obudowy znaleziono pozostałość kleju taśmy w miejscu korelującym z początkiem zakleszczenia. Była to wiarygodna hipoteza inicjatora: lepka pozostałość mogła zwiększyć lokalne tarcie. Zespół usunął pozostałość i powtórzył badanie. Zachowanie oraz zmierzone objętości nie zmieniły się. Wynik negatywny odrzucił wyjaśnienie „bieżąca pozostałość kleju jest wystarczającą przyczyną powtarzalnego zakleszczenia”. Nie dowiódł natomiast, że pozostałość nie mogła zainicjować pierwszego przechylenia ani utrwalić niekorzystnego ułożenia.msl-adhesive-test
To rozróżnienie jest kluczowe: korelacja w oględzinach utworzyła hipotezę, ale dopiero usunięcie czynnika bez ustąpienia objawu wykazało, że model był niekompletny. Czyszczenie nie mogło być zaakceptowaną korektą tylko dlatego, że znaleziono zanieczyszczenie.
Mechanizm samowzmacniający
Raport opisuje mechanizm:
tarcie osiowe → moment przechylający → kontakt z obudową
↑ ↓
└──── większe tarcie ← większa siła normalnaTarcie w prowadzeniu wytwarzało moment przechylający. Przechylony mieszek wyginał się bocznie i stykał z obudową. W miarę przemieszczania obszarów kontaktu potrzebna siła normalna rosła; wraz z nią rosło tarcie, które dalej zwiększało moment. Powstało samowzmacniające zakleszczenie, którego większa siła osiowa nie musiała zwolnić. Ten model wyjaśniał zarazem przedwczesny koniec skoku, przechylenie na obrazie i powtarzalność po usunięciu pozostałości.msl-mechanism
Zmiana skali konstrukcji miała znaczenie. W porównaniu z wcześniejszymi rozwiązaniami wymagano większej pojemności i znacznie dłuższego życia cyklicznego. Aby ograniczyć naprężenia zmęczeniowe, zastosowano bardziej podatny mieszek, co zmniejszyło sztywność osiową i boczną oraz zwiększyło podatność na przechylenie. Raport wskazał też, że istniejące tolerancje nie tworzyły wymiarowych granic gwarantujących działanie. Nie dało się więc zamknąć niepewności stwierdzeniem, że każda część mieści się na rysunku.msl-stiffness-limits
Decyzja, prowadzenie i dowód regresyjny
Ograniczenia programu wykluczały prostą wymianę całego mieszka, wydłużenie obudowy lub zmniejszenie wymaganego skoku. Korekta skupiła się na prowadzeniu końca mieszka i ograniczeniu możliwości przechylenia. Zespół badał warianty prowadnika, obciążenia promieniowe i momentowe oraz pełny skok, a wybraną zmianę poddał powtarzanym próbom. Decyzja wynikała z nowego mechanizmu — należało przerwać pętlę wzrostu siły normalnej — a nie z kosmetycznego usunięcia śladu.
Dowód regresyjny obejmował próbę ruchu po modyfikacji przy niekorzystnych obciążeniach łączonych. Raport ujawnia zarazem granice dowodu: geometria stanowiska nie odtworzyła w pełni założonego obciążenia promieniowego, część oceny miała charakter inżynierskiego osądu, a nie każdą próbę wykonano w całkowicie zamkniętej i ciśnieniowo identycznej konfiguracji. Wynik wspierał konkretną korektę w granicach badanego urządzenia; nie ustanowił uniwersalnych luzów, materiałów ani geometrii dla wszystkich mieszków.msl-regression
Uogólnienie jest metodologiczne, nie wymiarowe. Wymaganie pełnego skoku musi mieć pomiar objętości lub położenia. Anomalia kontaktowa wymaga modelu sił normalnych i momentów, nie tylko ciśnienia osiowego. Usunięcie podejrzanego zanieczyszczenia wymaga powtórzenia identycznej metryki. A korekta prowadzenia musi być sprawdzona po integracji i w reprezentatywnej kombinacji obciążeń.
Studium C: syntetyczny laboratoryjny tłumik pulsacji
Wszystkie liczby w tym studium są syntetyczne, neutralne i dydaktyczne. Opis dotyczy zamkniętej wodnej sieci laboratoryjnej małej mocy; nie przedstawia parametrów żadnego sprzętu wojskowego.
Wymaganie, stan i model R–L–C
Źródło narzuca sinusoidalną amplitudę szczytową strumienia $|\hat Q_s|=10{,}0 \mathrm{mL/s}$ w paśmie 15–40 Hz. Bez tłumika pozostała droga ma w dydaktycznym modelu rezystancyjną impedancję $Z_N=1{,}00\cdot10^9 \mathrm{Pa\,s/m^3}$, zatem amplituda ciśnienia odniesienia wynosi $|\hat p_0|=10{,}0 \mathrm{kPa}$. Wymagania są dwa:
- $|H|=|\hat p/\hat p_0|\le0{,}55$ w całym paśmie 15–40 Hz, łącznie z niepewnością rozszerzoną;
- szczytowy strumień gałęzi nie może przekroczyć $8{,}0 \mathrm{mL/s}$, aby zachować syntetyczny limit prędkości portu i skoku separatora.
Punkt pracy gazu: $p_{g0}=200 \mathrm{kPa(abs)}$, $T_0=293{,}15 \mathrm K$, $V_{g0}=2{,}40 \mathrm{mL}$, politropowy $n=1{,}20$. Stąd
$$C_h=\frac{2{,}40\cdot10^{-6}}{1{,}20\cdot200000} =1{,}00\cdot10^{-11} \mathrm{m^3/Pa}.$$
Syntetyczne przyłącze i dławik po identyfikacji mają $L_h=4{,}05\cdot10^6 \mathrm{Pa\,s^2/m^3}$ oraz $R_h=4{,}00\cdot10^8 \mathrm{Pa\,s/m^3}$. Otrzymujemy
$$f_0=\frac{1}{2\pi\sqrt{L_hC_h}}=25{,}01 \mathrm{Hz}, \qquad \zeta=\frac{R_h}{2}\sqrt{\frac{C_h}{L_h}}=0{,}314.$$
Model ma postać $Z_A=R_h+j(\omega L_h-1/\omega C_h)$ oraz $H=Z_A/(Z_N+Z_A)$. Jego przewidywania obliczono bez zaokrąglania wartości pośrednich; tabela zaokrągla moduł do trzech cyfr, fazę do $0{,}1°$, ciśnienie i przepływ do $0{,}001$ jednostki, a moc do $0{,}01 \mathrm{mW}$.
| $f$ | $|H|$ | $\arg H$ | $|\hat p|$ | $|\hat Q_A|$ | $\overline P_A=\tfrac12\mathrm{Re}(\hat p\hat Q_A^*)$ |
|---|---|---|---|---|---|
| 5 Hz | 0,917 | −17,2° | 9,169 kPa | 2,975 mL/s | 1,77 mW |
| 15 Hz | 0,507 | −33,6° | 5,066 kPa | 6,426 mL/s | 8,26 mW |
| 25 Hz | 0,286 | −0,1° | 2,857 kPa | 7,143 mL/s | 10,20 mW |
| 40 Hz | 0,482 | +33,3° | 4,819 kPa | 6,531 mL/s | 8,53 mW |
| 80 Hz | 0,814 | +25,0° | 8,140 kPa | 4,330 mL/s | 3,75 mW |
Poniżej $f_0$ dominuje reaktancja podatnościowa i faza $H$ jest ujemna; powyżej — inertancyjna i faza zmienia znak. Przy 25 Hz gałąź jest blisko rezonansu szeregowego, ma najmniejszą impedancję i tworzy minimum transmitancji węzłowej. Szczytowy skok objętości w tym punkcie wynosi
$$|\hat V_l|=\frac{|\hat Q_A|}{2\pi f} =\frac{7{,}143\cdot10^{-6}}{2\pi\cdot25} =4{,}55\cdot10^{-8} \mathrm{m^3}=0{,}0455 \mathrm{mL}.$$
Moc $10{,}20 \mathrm{mW}$ nie wynika z samego przesunięcia fazowego, lecz z iloczynu amplitud i części rzeczywistej. Kontrola energetyczna przy rezonansie daje $\hat p=2{,}857 \mathrm{kPa}$, $\hat Q_A=7{,}143 \mathrm{mL/s}$ niemal w fazie, a więc $\overline P_A=0{,}5\cdot2857\cdot7{,}143\cdot10^{-6}=0{,}01020 \mathrm W$.
Doświadczenia NASA z pierścieniowym akumulatorem pokazały praktyczną identyfikację podatności, inertancji i oporu z czasowych przebiegów ciśnienia w wielu punktach oraz odchylenie od prostego modelu kryzy zależne od średniej prędkości. To bezpośredni powód, aby parametrów modelu nie traktować jako niezmiennych danych geometrycznych.nasa-annular-test
Wariant nieudany i korekta
Pierwszy syntetyczny wariant miał ten sam $L_h$ i $C_h$, lecz $R_h=5{,}00\cdot10^7 \mathrm{Pa\,s/m^3}$, czyli $\zeta=0{,}039$. Przy 25 Hz dawał bardzo głębokie minimum $|H|=0{,}0476$, ale szczytowy strumień gałęzi $9{,}524 \mathrm{mL/s}$. Spełniał widowiskową metrykę ciśnienia, lecz przekraczał limit 8,0 mL/s, zawężał użyteczne tłumienie i zwiększał wrażliwość na odstrojenie.
Decyzja: wariant odrzucono. Korekta zwiększyła kontrolowaną dyssypację do $R_h=4{,}00\cdot10^8$, obniżając przepływ rezonansowy do 7,143 mL/s kosztem płytszego, lecz szerszego minimum ciśnienia. To nie jest zasada „więcej oporu jest lepiej”: dalszy wzrost $R_h$ odłączyłby podatność i podniósł $|H|$. Optimum wynika równocześnie z pasma, przepływu gałęzi, mocy i niepewności.
Pomiar, niepewność i decyzja
Stanowisko mierzy synchronicznie: ciśnienie węzła, ciśnienie komory za portem, strumień źródła i strumień gałęzi; dodatkowo temperaturę gazu, cieczy i korpusu oraz średnie ciśnienie absolutne. Przetworniki ciśnienia i przepływu mają zweryfikowane płaskie pasmo co najmniej 0–200 Hz, wspólne próbkowanie 2 kHz, analogowy filtr antyaliasingowy 800 Hz i błąd synchronizacji poniżej $50 \mathrm{\mu s}$. Przy 40 Hz odpowiada to najwyżej $0{,}72°$ błędu czasu. Każdy punkt utrzymuje się 60 s; fazory wyznacza się z widma wzajemnego względem strumienia źródła, odrzucając odcinek ustalania.
Syntetyczny budżet modelu zawiera $u(C_h)/C_h=2\%$, $u(L_h)/L_h=1{,}5\%$, $u(R_h)/R_h=2{,}5\%$ i $u(Z_N)/Z_N=1{,}5\%$ jako niepewności standardowe. Budżet pomiaru obejmuje wzorcowanie amplitudy, szum, przeciek widmowy, synchronizację, powtarzalność stanu gazu i niepewność korekcji fazy. Po propagacji i zaokrągleniu przyjęto rozszerzone $U(|H|)=0{,}025$, $U(\arg H)=2{,}0°$ oraz względne $U(\overline P_A)=9\%$ dla $k=2$.
Dla amplitudy $|\hat Q_A|$ syntetyczne standardowe składniki wynoszą: wzorcowanie $0{,}08 \mathrm{mL/s}$, szum i przeciek widmowy $0{,}07 \mathrm{mL/s}$, powtarzalność $0{,}08 \mathrm{mL/s}$ oraz korekcja stanu i pasma $0{,}04 \mathrm{mL/s}$. Przy braku korelacji $u_c(Q_A)=\sqrt{0{,}08^2+0{,}07^2+0{,}08^2+0{,}04^2}=0{,}139 \mathrm{mL/s}$, stąd $U(Q_A)=2u_c=0{,}278 \mathrm{mL/s}$, po zaokrągleniu $0{,}28 \mathrm{mL/s}$. Reguła decyzji dla górnych granic jest stosowana w każdym z 26 punktów $f=15,16,\ldots,40 \mathrm{Hz}$: $|H|_{meas}+U(|H|)\le0{,}55$ oraz $|\hat Q_A|_{meas}+U(Q_A)\le8{,}0 \mathrm{mL/s}$.
Syntetyczny zapis próby kontrolnej jest następujący:
| $f$ | model $|H|\angle\phi$ | pomiar syntetyczny $|H|\angle\phi$ | zmierzone syntetyczne $|\hat Q_A|$ | model $\overline P_A$ | pomiar syntetyczny | Decyzja punktowa |
|---|---|---|---|---|---|---|
| 15 Hz | $0{,}507\angle-33{,}6°$ | $0{,}515\angle-34{,}2°$ | 6,48 mL/s | 8,26 mW | 8,4 mW | zgodny z modelem i limitem |
| 25 Hz | $0{,}286\angle-0{,}1°$ | $0{,}293\angle+0{,}6°$ | 7,25 mL/s | 10,20 mW | 10,5 mW | zgodny; minimum transmitancji |
| 40 Hz | $0{,}482\angle+33{,}3°$ | $0{,}489\angle+34{,}0°$ | 6,60 mL/s | 8,53 mW | 8,7 mW | zgodny z modelem i limitem |
Niezależne obliczenie modelu bez zaokrąglania wartości pośrednich oraz jawnie syntetyczny pomiar co $1 \mathrm{Hz}$ dały następujące ekstrema pełnego przemiatania:
| Wielkość w 26 punktach 15–40 Hz | Maksimum | Częstotliwość |
|---|---|---|
| modelowe $|H|$ | 0,506612 | 15 Hz |
| zmierzone syntetyczne $|H|$ | 0,515 | 15 Hz |
| modelowe $|\hat Q_A|$ | 7,142857 mL/s | 25 Hz |
| zmierzone syntetyczne $|\hat Q_A|$ | 7,25 mL/s | 25 Hz |
Najgorszy przedział transmitancji występuje przy 15 Hz: $[0{,}515-0{,}025;0{,}515+0{,}025]=[0{,}490;0{,}540]$, więc jego górny kraniec pozostaje o $0{,}010$ poniżej limitu 0,55. Najgorszy przedział szczytowego przepływu występuje przy 25 Hz: $[7{,}25-0{,}28;7{,}25+0{,}28]=[6{,}97;7{,}53] \mathrm{mL/s}$, a górny kraniec ma $0{,}47 \mathrm{mL/s}$ rezerwy do limitu $8{,}0 \mathrm{mL/s}$. Wszystkie pozostałe punkty mają mniejsze górne krańce. Decyzja: poprawiony tłumik spełnia oba syntetyczne wymagania w konfiguracji stanowiska.
Nie jest to jeszcze decyzja integracyjna. $Z_N$ założono rezystancyjne, podczas gdy reprezentatywny przewód wniesie propagację, odbicia i dodatkowe węzły. Po cyklach temperaturowych i mechanicznych powtarza się identyfikację $C_h,L_h,R_h$, krzywą $P$–$V$, próg ruszenia i przeciek. Następnie montuje się reprezentatywny przewód i ponownie mierzy $H$, $Q_A$ oraz moc na obu końcach pasma i w minimum. Kryteria pozostają te same; przesunięcie minimum nie może być „skorygowane” zmianą osi wykresu.
Regresja kończy się dopiero, gdy korytarz $|H|+U\le0{,}55$, limit przepływu, znakowa zmiana fazy, moc oraz brak nowych harmonicznych są zachowane po temperaturze, cyklach i integracji. To praktyczna realizacja hierarchii: model skupiony odpowiada najpierw na pytanie o mechanizm, model rozłożony na pytanie o sieć, a próba na pytanie o zgodność. Studium numeryczne bez tej ostatniej pętli pozostaje hipotezą; metodologia walidacji modeli i porównywania wariantów jest szeroko rozwinięta w pracy Tarnawskiego, choć jej przedmiotem jest inna klasa układów przepływowych.local-tarnawski-validation
Wnioski: jedna podatność nie oznacza jednej funkcji
Akumulator magazynujący energię rozlicza się z pracy po różnicy ciśnień i objętości oddanej. Kompensator rozlicza się z bilansu całej pętli, skrajnych temperatur i zapasu skoku. Tłumik pulsacji rozlicza się z zespolonej transmitancji, przepływu gałęzi i mocy w paśmie. Ogranicznik udaru rozlicza się z czasu propagacji, szczytu, impulsu oraz minimalnego ciśnienia absolutnego. Wspólny gaz i separator nie pozwalają zastąpić tych dowodów jednym testem.
Teza o zależnej od stanu impedancji bocznej porządkuje również awarie. Utrata gazu zmienia krzywą $P$–$V$ i $C_h$. Zablokowany port zmienia $L_h$ i $R_h$. Tarcie separatora zmienia równowagę sił, histerezę i moc. Ogranicznik usuwa użyteczną podatność. Przewód rozłożony zmienia $Z_N$ i położenie węzłów. Każde zjawisko ma przewidywaną sygnaturę, a każda sygnatura wymaga kanału pomiarowego oraz budżetu niepewności.
Najkrótszy poprawny proces nie polega więc na pominięciu ogniw, tylko na dobraniu ich rozdzielczości do ryzyka. Najpierw mierzalne wymaganie; potem zapis stanu, granicy i znaków; model wystarczający dla skali czasu; przewidywanie przed próbą; zsynchronizowany pomiar; jawna niepewność; przedziałowe kryterium; decyzja bez przesuwania bramki; na końcu regresja po zmianie i środowisku. Dopiero cały ten łańcuch odróżnia działający element od urządzenia, które jedynie ma właściwą nazwę i ciśnienie na manometrze.