Mach stem i wysokość wybuchu: pomiary mikrofonowe, szybka kamera i symulacja

Streszczenie

Mach stem jest pionowym, wzmocnionym frontem fali uderzeniowej powstającym wtedy, gdy fala padająca i odbita od powierzchni przestają być dwiema osobnymi falami. Dla wybuchu powietrznego oznacza to, że geometria wysokości, odległości od punktu zero i kąta padania może zmienić rozkład nadciśnienia przy gruncie silniej niż sugerowałaby sama energia w kt TNT. Ten artykuł pokazuje Mach stem jako problem pomiarowy: trzeba mierzyć czas przybycia, przebieg ciśnienia, prędkość czoła, wysokość punktu potrójnego i zgodność z modelem, a nie tylko narysować większy promień skutków.^1,2

Tekst rozwija wcześniejszy artykuł o wysokości detonacji a fali podmuchowej, ale kładzie nacisk na aparaturę i modele: rząd mikrofonów lub czujników ciśnienia, szybką kamerę, ślady dymne lub metody optyczne, rekonstrukcję trajektorii punktu potrójnego oraz porównanie z uproszczonym kalkulatorem Height-of-burst i Mach stem. W ujęciu dydaktycznym nie chodzi o operacyjne dobieranie wysokości wybuchu, lecz o zrozumienie, dlaczego model HOB bez pomiarowej walidacji jest tylko szkicem zjawiska.

Rozszerzenie tematu

Zwykłe odbicie fali od ściany albo podłoża można sobie wyobrazić jak piłkę odbijającą się od podłogi: coś nadchodzi, coś wraca, a punkt kontaktu leży na powierzchni. Dla fali uderzeniowej ten obraz działa tylko przez chwilę. Jeżeli fala padająca dociera do sztywnego podłoża pod odpowiednim kątem i z odpowiednią siłą, fala odbita może dogonić falę padającą, a ich układ przechodzi w odbicie nieregularne. Wtedy między gruntem i punktem spotkania fal powstaje trzeci front, czyli Mach stem.¹

Geometria ma tu trzy główne elementy. Pierwszy to fala padająca idąca od punktu wybuchu. Drugi to fala odbita od podłoża. Trzeci to punkt potrójny, w którym spotykają się fala padająca, odbita i pionowy front Macha. Gdy ten punkt porusza się w górę i na zewnątrz, wysokość Mach stem rośnie, a pole ciśnienia przy powierzchni nie jest już prostą wersją fali swobodnej.²

Dlatego Mach stem nie jest tylko ciekawym rysunkiem z podręcznika fal uderzeniowych. W analizie skutków wybuchu powietrznego decyduje o tym, czy przy gruncie widzimy lokalnie słabszą falę swobodną, falę odbitą, czy połączony front o innym rozkładzie ciśnienia. To bezpośrednio łączy się z równoważnikiem trotylowym, nadciśnieniem i impulsem: sama energia nie wystarcza, jeżeli nie znamy geometrii odbicia.

Podstawy fizyczne: warunki Rankine'a-Hugoniota i zachowanie fali

Fala uderzeniowa to nieciągłość: po jednej stronie frontu gaz ma gęstość ρ₁, ciśnienie P₁ i prędkość u₁, po drugiej ρ₂, P₂ i u₂. Przejście przez front podlega prawom zachowania masy, pędu i energii, znanych jako warunki Rankine'a-Hugoniota. Dla fali uderzeniowej w powietrzu przy założeniu gazu doskonałego:

ρ₂/ρ₁ = (γ+1)·M² / [(γ-1)·M² + 2]

P₂/P₁ = [2γ·M² - (γ-1)] / (γ+1)

gdzie M to liczba Macha fali (prędkość frontu do prędkości dźwięku za falą), a γ ≈ 1,4 dla powietrza. Dla słabych fal M → 1 dostajemy P₂ ≈ P₁ — fala dźwiękowa. Dla silnych M >> 1 stosunek gęstości osiąga granicę (γ+1)/(γ-1) = 6 dla powietrza. Silna fala uderzeniowa nie może zatem zagęścić powietrza bardziej niż sześciokrotnie, niezależnie od nadciśnienia.

Nadciśnienie ΔP = P₂ - P₁ jest tym parametrem, który najczęściej cytują modele blast. W artykule o Kingery-Bulmash, Brode i CONWEP tabele K-B dostarczają ΔP(Z) jako funkcję skalowanej odległości Z = R/W^(1/3), gdzie R to odległość w metrach, a W to masa TNT w kilogramach. To skalowanie Hopkinsona-Cranza wynika właśnie z warunków Rankine'a-Hugoniota: dla geometrycznie podobnych wybuchów w powietrzu wszystkie parametry blast zależą tylko od Z.^1,5

W odbiciu od sztywnego podłoża warunki Rankine'a-Hugoniota muszą być spełnione dwa razy: raz dla fali padającej, raz dla fali odbitej. Stąd wynika ciśnienie odbite Pr, które jest kilkakrotnie większe od ciśnienia padającego Pso. Dla słabych fal współczynnik odbicia Cr = Pr/Pso zbliża się do 2 (odbicie akustyczne), ale dla silnych fal może sięgać 8 lub więcej. To wyjaśnia, dlaczego ściana prostopadła do kierunku propagacji fali dostaje znacznie większe obciążenie niż ściana równoległa.

Odbicie regularne i nieregularne

W odbiciu regularnym fala padająca i odbita spotykają się na powierzchni. Można wtedy mówić o kącie padania, kącie odbicia i lokalnym wzroście ciśnienia przy podłożu. To nadal jest problem gazodynamiczny, bo za frontem fali zmieniają się prędkość, temperatura, gęstość i ciśnienie, ale geometria pozostaje względnie prosta.

W odbiciu nieregularnym układ ma już trzy fale. Front Macha porusza się prawie równolegle do powierzchni, a punkt potrójny unosi się nad gruntem. Nad i pod tym punktem gaz przeszedł przez różne sekwencje frontów, dlatego w obrazie mogą pojawić się też nieciągłości kontaktowe i poślizgowe. Dla studentów najważniejsze jest to, że ciśnienie wzdłuż wysokości Mach stem nie musi być jednorodne.²

Właśnie ta niejednorodność sprawia, że pojedynczy czujnik na gruncie może dać zbyt ubogi obraz. Pomiar na powierzchni mówi, kiedy fala przyszła i jaki był przebieg ciśnienia w tym punkcie. Nie mówi automatycznie, co działo się pół metra, metr albo kilkadziesiąt metrów wyżej wzdłuż frontu. Dlatego w badaniach nad Mach stem potrzebne są szeregi czujników i dane optyczne, a nie jeden punkt pomiarowy.^1,2

Dlaczego wysokość wybuchu zmienia wynik

Dla wybuchu powietrznego przy powierzchni gruntu liczą się dwie odległości: wysokość nad gruntem i odległość pozioma od punktu zero. Fala swobodna idzie po odległości skośnej, a fala odbita tworzy dodatkową ścieżkę i zmienia kąt, pod którym energia trafia w powierzchnię. Gdy odległość pozioma rośnie, kąt padania staje się coraz bardziej płaski, a warunki przejścia do odbicia Macha zmieniają się razem z osłabianiem fali.

To prowadzi do pozornie paradoksalnego wniosku: wybuch nie może być po prostu „jak najniżej” albo „jak najwyżej”. Przy zbyt niskim wybuchu rośnie rola lokalnego kraterowania, gruntu i silnego sprzężenia z powierzchnią. Przy zbyt wysokim wybuchu fala przy gruncie może być bliższa osłabionej fali swobodnej. Pomiędzy tymi skrajnościami istnieją zakresy, w których odbicie i Mach stem zwiększają promień danego progu nadciśnienia.³

W tym miejscu pojawia się pojęcie height of burst, czyli HOB. W analizie dydaktycznej nie należy go traktować jako jednej uniwersalnej „optymalnej wysokości”. Optimum zależy od progu nadciśnienia, uzysku, celu, wysokości obiektu, topografii i tego, czy interesuje nas grunt, ściana, dach, wieża, samolot czy czujnik w pionowej osi pomiarowej. Ten sam HOB, który zwiększa obszar przy jednym progu, nie musi być optymalny dla innego progu ani innej wysokości celu.^2,3

Co mierzy czujnik na gruncie

W starszych danych empirycznych najczęściej mierzono maksymalne nadciśnienie i czas przybycia fali. Jeżeli czujniki ciśnienia ustawiono w szeregu, można było wyznaczyć także średnią prędkość czoła. Znacznie rzadziej rejestrowano pełną historię pressure-time, bo aparatura była trudniejsza, a fazę ujemną często uznawano za mniej istotną dla typowego projektowania ochronnego.¹

W nowszych badaniach nadal dominują czujniki ciśnienia, ale coraz większą rolę mają metody optyczne. Omang i współautorzy opisują użycie technik fotograficznych, śladów dymnych oraz szybkiego filmowania do wyznaczania ruchu cząstek i gęstości. Wspominają też o pełnoskalowym obrazowaniu typu Schlieren, chociaż podkreślają, że takie metody są praktycznie najbardziej użyteczne w wizualizacji i precyzyjnej pracy laboratoryjnej.¹

Z punktu widzenia serwisu najważniejszy jest komplet danych, który pozwala porównać eksperyment z modelem. Minimalny zestaw obejmuje:

• Pmax, czyli maksymalne nadciśnienie,
• t_toa, czyli czas przybycia czoła fali,
• t+, czyli czas dodatniej fazy,
• I, czyli impuls dodatniej fazy,
• położenie czoła fali w kolejnych chwilach,
• informację, czy w danym punkcie mamy odbicie regularne czy Mach stem.

Sama wartość Pmax jest zbyt mało informacyjna. W badaniach porównawczych Height of burst explosions zgodność z danymi empirycznymi była dobra dla nadciśnienia i czasu przybycia, ale dla czasu dodatniej fazy i impulsu różnice między bazami danych były już widoczne. To dobra lekcja metrologiczna: parametry całkowe i czasowe są trudniejsze niż pik.¹

Mikrofony i czujniki ciśnienia

W ćwiczeniu dydaktycznym słowo „mikrofon” trzeba rozumieć ostrożnie. Zwykły mikrofon audio dobrze pokazuje czas przybycia impulsu akustycznego i odbicia, ale nie jest automatycznie czujnikiem nadciśnienia fali uderzeniowej. Do ilościowego pomiaru blastu używa się czujników ciśnienia o znanym paśmie, czasie narastania, zakresie dynamicznym i odporności mechanicznej. W bezpiecznym ćwiczeniu akademickim można jednak użyć mikrofonów do zrozumienia geometrii czasów przybycia bez generowania rzeczywistej fali uderzeniowej.

Układ pomiarowy dla analogii akustycznej może mieć rząd mikrofonów położonych wzdłuż linii od źródła impulsu oraz drugi rząd na małych wysokościach nad powierzchnią. Z różnic czasu przybycia można wyznaczyć prędkość czoła, wykryć odbicia od podłoża i sprawdzić, gdzie sygnały zaczynają wyglądać jak złożenie fali bezpośredniej i odbitej. To nie jest pełny Mach stem, ale jest dobrym wstępem do myślenia o trajektorii czoła i o tym, dlaczego jeden czujnik nie wystarcza.

W realnym pomiarze blastowym najważniejsze problemy aparaturowe to pasmo, czas narastania, kalibracja i montaż. Czujnik wystający ponad powierzchnię może sam zaburzać falę. Czujnik wpuszczony w płytę mierzy inny przypadek geometryczny niż sonda w wolnym polu. Długi przewód może zniekształcić impuls, a zbyt wolna akwizycja wygładza szczyt. Te same problemy pojawiają się w artykułach o oscyloskopie w impulsach nanosekundowych i jitterze synchronizacji, tylko tutaj sygnałem jest ciśnienie zamiast napięcia.

Szybka kamera i obrazowanie frontu

Czujnik ciśnienia daje wykres w punkcie. Kamera daje geometrię. Przy odpowiednim kontraście można śledzić front fali, cień zaburzenia gęstości, ślady dymne albo ruch lekkich markerów. Jeżeli obraz jest zsynchronizowany z akwizycją ciśnienia, można przypisać pikom w sygnale konkretne położenie frontu, odbicia i punktu potrójnego.¹

Szybka kamera nie rozwiązuje jednak problemu sama. Dla słabnącej fali kontrast spada, a dokładne położenie punktu potrójnego może być trudne do wyznaczenia. Obraz 2D może mieszać zjawiska spoza płaszczyzny obserwacji. Tło, perspektywa, ostrość i częstotliwość klatek mogą ograniczać dokładność. Dlatego dane optyczne najlepiej traktować razem z czujnikami ciśnienia: kamera podpowiada geometrię, a czujniki dają lokalny przebieg P(t).

W badaniach pełnoskalowych szczególnie kuszące są metody typu Schlieren i shadowgraph, bo bezpośrednio pokazują gradienty gęstości. Są jednak trudne technicznie na dużą skalę. W dydaktyce warto zacząć od prostszego poziomu: symulacja numeryczna frontów, analiza nagrania demonstracyjnego albo optyczny eksperyment z falą na płytkiej wodzie jako analogią geometrii odbicia. Taka analogia nie odtwarza gazodynamiki, ale pozwala zrozumieć, czym jest punkt potrójny.

Punkt potrójny jako główny obiekt pomiaru

W modelu Mach stem najcenniejszym obiektem nie jest pojedyncze maksimum ciśnienia, lecz trajektoria punktu potrójnego. Jeżeli wiemy, gdzie w kolejnych chwilach znajduje się punkt przecięcia trzech frontów, możemy rozdzielić strefę odbicia regularnego, strefę Macha i obszar, w którym profil pionowy ciśnienia staje się niejednorodny.

Praca Mach Stem Modeling with Spherical Shock Waves zwraca uwagę, że w sferycznym szoku ciśnienie w Mach stem może najpierw rosnąć, a potem maleć względem fali padającej, a rozkład wzdłuż wysokości może mieć wyraźne „wybrzuszenie” nadciśnienia. Autor podkreśla też, że czas pomiaru jest krytyczny: bardzo wcześnie, pośrednio i późno pionowy profil ciśnienia może wyglądać inaczej.²

To prowadzi do ważnej konsekwencji dla wysokich obiektów. Maksymalizacja nadciśnienia na gruncie nie jest tym samym, co maksymalizacja obciążenia na określonej wysokości nad gruntem. Budynek, maszt, wieża, czujnik na wysokości albo samolot w pobliżu frontu nie „widzą” tylko wartości z mapy powierzchniowej. Jeżeli interesuje nas obiekt przestrzenny, trzeba znać wysokość punktu potrójnego i pionowy profil ciśnienia.²

Jak działa lokalny kalkulator HOB

Lokalny kalkulator Height-of-burst i Mach stem jest celowo opisowy. Jego model robi trzy rzeczy: liczy falę padającą po odległości skośnej modelem Kingery-Bulmash, nakłada prosty współczynnik odbicia zależny od kąta padania i dodaje gładki mnożnik Mach stem dopiero wtedy, gdy odległość skalowana jest wystarczająco duża względem wysokości skalowanej.⁴

To nie jest pełna tabela CONWEP ani kompletny model przejścia regularne-nieregularne. Kod mówi to wprost: celem jest dydaktyczne pokazanie, dlaczego dla danego progu nadciśnienia może istnieć optimum wysokości, a nie projektowa predykcja. Dodatkowo kalkulator pokazuje wskaźnik słabego szoku oparty na bezwymiarowym parametrze a, który ma charakter diagnostyczny, nie normowy.⁴

Wynik kalkulatora należy więc czytać w trzech warstwach:

• jako geometrię: wysokość, odległość skośna, odległość pozioma, kąt padania,
• jako porównanie reżimów: fala swobodna, powierzchniowa, odbicie regularne, Mach stem,
• jako ostrzeżenie metrologiczne: pełny wynik wymagałby tablic, czujników i walidacji z eksperymentem.

Takie podejście jest lepsze dydaktycznie niż jedna liczba „optymalnej wysokości”. Student widzi, że zmiana progu nadciśnienia zmienia skan wysokości, a zmiana odległości zmienia reżim odbicia.

Symulacja: od RSPH do prostego modelu strony

Omang i współautorzy pokazują, że symulacje RSPH z prostym modelem wybuchu mogą dobrze odtwarzać ciśnienie i czas przybycia dla testów wolnego pola oraz HOB. Jednocześnie porównanie baz empirycznych ujawnia różnice w czasie dodatniej fazy i impulsie. To znaczy, że „zgodność z eksperymentem” trzeba rozbijać na parametry: model może dobrze przewidywać Pmax, a słabiej przewidywać całkę po czasie.¹

W praktyce mamy więc hierarchię narzędzi. Najniżej jest model jakościowy, taki jak lokalny kalkulator. Wyżej są empiryczne procedury K-B/CONWEP i Brode, które dają szybkie wyniki dla typowych geometrii. Jeszcze wyżej są symulacje CFD, ALE, RSPH albo sprzężone modele FEM, które potrafią śledzić geometrię frontu i odpowiedź konstrukcji, ale wymagają większej liczby założeń i walidacji.^1,4

Najgorszym błędem jest mieszanie poziomów. Jeżeli prosty kalkulator ma charakter dydaktyczny, nie wolno z niego robić normowego narzędzia projektowego. Jeżeli symulacja daje piękny obraz frontu, nie wolno zapominać, że jej parametry źródła, równanie stanu, siatka, warunki brzegowe i model turbulentny też wnoszą niepewność. W obu przypadkach potrzebna jest kontrola względem pomiaru.

Jak projektować bezpieczne ćwiczenie

Ćwiczenie akademickie nie powinno polegać na wytwarzaniu fali uderzeniowej z materiału wybuchowego. Bezpieczny wariant ma trzy części: akustyczny analog czasów przybycia, symulację geometrii odbicia oraz analizę danych syntetycznych z czujników. Każda część uczy innego fragmentu problemu, a razem dają zrozumienie Mach stem bez niebezpiecznego protokołu laboratoryjnego.

W części akustycznej można użyć impulsu dźwiękowego, rzędu mikrofonów i płaskiej powierzchni odbijającej. Celem nie jest uzyskanie nieliniowego szoku, tylko rekonstrukcja fali bezpośredniej i odbitej z czasów przybycia. Studenci mierzą opóźnienia, porównują je z geometrią i widzą, jak łatwo pomylić odbicie z osobnym zdarzeniem, jeżeli nie ma kilku kanałów.

W części symulacyjnej można użyć prostego skryptu rysującego kręgi fali padającej i odbitej metodą obrazu zwierciadlanego. Następnie dodaje się warunek przejścia do frontu Macha jako regułę dydaktyczną, a nie model fizyczny. Studenci śledzą punkt przecięcia frontów i zapisują jego wysokość jako funkcję odległości.

W części danych syntetycznych otrzymują tabelę: położenie czujnika, czas przybycia, Pmax, t+, impuls i etykietę reżimu. Ich zadaniem jest odtworzyć, gdzie zaczęło się odbicie nieregularne, oraz wskazać, które parametry są najbardziej niepewne. To ćwiczenie dobrze łączy się z artykułem o Kingery-Bulmash, Brode i CONWEP, bo pokazuje, jakie dane w ogóle trafiają do takich modeli.

Historia badań eksperymentalnych: od WWII do prób jądrowych

Badania Mach stem i wysokości wybuchu mają bezpośrednie korzenie w programach broni z okresu II wojny światowej i pierwszej dekady zimnowojennej.

Operacja Trinity (1945). Pierwsza detonacja jądrowa wygenerowała falę blastową, którą mierzono sieciami czujników ciśnienia i seimofoniczno-akustycznych rozmieszczonych w promieniu kilku mil. Dane te, zebrane przez zespół Los Alamos i Oppenheimera, pozwoliły na wstępne potwierdzenie, że sferyczna fala blast zachowuje się zgodnie z prawem sześciennym skalowania. Czas przybycia fali zmierzono mikrofonami kondensatorowymi i fotokomórkami — mniej dokładnie niż późniejsze czujniki piezoelektryczne, ale wystarczająco do rekonstrukcji prędkości frontu.⁶

Testy na Bikini i Enewetak (1946–1952). Operacja Crossroads i kolejne testy dostarczyły pierwszych danych o wybuchu powietrznym nad morzem i lądem. Mierzono przede wszystkim efekty na okrętach, budynkach i manekiny do oceny skutków biologicznych. Opracowania z tych testów (dostępne dziś przez DTRA i DOE OSTI) zawierają szczegółowe wykresy HOB i ciśnień.

Castle Bravo (1954). Wybuch 15 Mt przy stosunkowo niskiej wysokości wytworzył ogromny krater i wielkie skażenie Bikini. Z punktu widzenia blast physics był to przypadek wybuchu prawie-powierzchniowego, gdzie fala Macha tworzy się natychmiast i nie ma „wolnego pola" dla czystego sferycznego szoku. Analiza Castle Bravo wpłynęła na decyzje o wyższych HOB dla ładunków termonuklearnych w następnych testach.

Serie Nevada (NTS, lata 50.–60.). Setki detonacji od kilkutonowych do kilkusetnukilotonowych ze szczegółowymi układami czujników. Kingery i Bulmash zebrali te dane i skompilowali w tabele wielomianowe — to właśnie podstawa artykułu o CONWEP i modelach Kingery-Bulmash. Mapy HOB z NTS są do dziś podstawowym odniesieniem empirycznym.

Próby radzieckie (Semipalatynsk, Nowa Ziemia). Radzieckie bazy danych blast są mniej dostępne, ale analiza porównawcza z lat 90. (po otwarciu części archiwów) wykazała dobrą zgodność z modelami K-B w szerokim zakresie kalibrów. Różnice wynikały głównie z różnic w topografii i warunkach atmosferycznych testów, a nie z odmiennej fizyki fali.

Prawa skalowania: Hopkinson-Cranz i ich granice

Prawo skalowania Hopkinsona-Cranza mówi, że dwa wybuchy o różnych energiach dają identyczne rozkłady ciśnienia, jeżeli porówna się je w skalowanych odległościach Z = R/W^(1/3). Skalowanie działa dlatego, że dla gazu doskonałego w geometrii sferycznej wszystkie długości skalują się sześciennym pierwiastkiem energii.

Konsekwencje dla Mach stem:

• Wysokość HOB skaluje się tak samo: optymalna h/W^(1/3) jest taka sama dla 1 kt i dla 1 Mt.
• Czas przybycia skaluje się jak W^(1/3) / c₀, gdzie c₀ to prędkość dźwięku.
• Impuls ciśnienia skaluje się jak W^(1/3).

Granice prawa skalowania: atmosfera nie jest gęstością-bezwymiarową. Przy porównaniu małych i dużych ładunków zmienia się prędkość dźwięku (zależna od temperatury i wysokości), gęstość powietrza (zależna od wysokości) i wpływ grawitacji. Dla ładunków termonuklearnych o Mt kulę ognistą musi się też skalować z energią — i tu pojawiają się efekty termiczne i radiacyjne wpływające na kształt fazy wstępnej. Prawo sześcienne jest dobrą aproksymacją dla konwencjonalnych zakresów, ale należy je traktować jako przybliżenie.

Innym ograniczeniem jest ciśnienie atmosferyczne na wysokości. Badania HOB w górach lub na dużych wysokościach (np. wybuchy w stratosferze dla efektów EMP) wymagają korekty gęstości atmosferycznej, bo skalowanie odnosi się do warunków standardowych.⁵

HOB dla celów jądrowych: różne kryteria, różne optimum

W analizie efektów taktycznych i strategicznych bomby jądrowej wybór HOB zależy od tego, jakie kryterium chce się maksymalizować:

Kryterium	Konsekwencja dla HOB
Maksymalna powierzchnia przy progu 10 psi	Wyższy HOB — Mach stem rozszerza zasięg poziomy
Maksymalne ciśnienie w punkcie zero	Eksplozja przy/na powierzchni
Minimalizacja krateru i skażenia gruntowego	Wyższy HOB — unikamy kontaktu z gruntem
Skutki EMP (impuls elektromagnetyczny)	Bardzo wysoki HOB (lub kosmiczny)
Zniszczenie twardych bunkrów	Penetrator o niskim HOB lub uderzenie w ziemię

Ta tabela pokazuje, że nie istnieje jeden „optymalny" HOB — każde kryterium efektu ma swoje optimum. W wojskowych doktrynach jądrowych historia doboru HOB dla konkretnych rodzajów celów jest szczegółowo udokumentowana w odtajnionych dokumentach NWS (Nuclear Weapon School) i SNL (Sandia National Laboratories), dostępnych przez DTRA RDAT&E.

Wzmocnienie Mach stem w terenie zabudowanym

Teren zabudowany zmienia warunki odbicia zasadniczo. Zamiast jednej płaskiej powierzchni gruntu fala blast trafia na trójwymiarowy labirynt ścian, dachów, uliczek i dziedzińców. Konsekwencje:

1. Wielokrotne odbicia. Każda ściana generuje własną falę odbitą. Interferencja tych fal może tworzyć lokalne wzmocnienia ciśnienia lub osłabienia.

2. Kanalizacja ciśnienia. Wąskie uliczki prowadzą falę jak rury — attenuation jest mniejsze niż w otwartym terenie, bo geometria bloków budynków ogranicza boczną ekspansję fali.

3. Cień blastowy. Za masywnym budynkiem ciśnienie może być znacząco mniejsze niż przed nim. Zjawisko cienia blastowego było dobrze widoczne w Hiroszimie, gdzie niektóre budynki bezpośrednio za większymi strukturami (np. bankiem) doznały mniejszych uszkodzeń.

4. Mach stem w wąwozie ulic. Fala Macha między budynkami może osiągać inną wysokość punktu potrójnego niż w otwartym terenie, bo ściany tworzą dodatkowe powierzchnie odbicia.

Symulacje numeryczne terenu zabudowanego (np. SHAMRC, AUTODYN, Ansys LS-DYNA) wymagają modelowania geometrii budynków, co jest komputerowo drogie, ale coraz bardziej dostępne. Badania m.in. US Army ERDC (Engineering Research and Development Center) i UK DSTL obejmują szczegółowe modele blast w warunkach urban.¹

Rodzaje czujników i ich charakterystyki techniczne

Praktyczny pomiar blastowy wymaga wyboru czujników zoptymalizowanych do konkretnych parametrów. Przegląd techniczny:

Czujniki piezoelektryczne (PCB Piezotronics, Kistler). Najszerzej stosowane do pomiaru ciśnienia fali uderzeniowej. Czas narastania < 1 µs, zakres do kilkudziesięciu MPa (kilkadziesiąt tys. psi). Problem: czułość na drgania mechaniczne (artefakty podczas przejścia frontu przez podstawę czujnika). Wymagają kondycjonowania sygnału i separacji galwanicznej.

Czujniki piezorezystancyjne (Endevco 8530C). Szeroko stosowane w rządowych laboratoriach USA. Liniowe w szerokim zakresie. Czas narastania < 2 µs. Wymagają separacji galwanicznej. Rejestracja fazy ujemnej możliwa (w odróżnieniu od niektórych piezoelektrycznych).

Czujniki pojemnościowe (condenser microphone). Dobre dla ciśnień bliskich akustycznym i niskich eksplozji. W blastach wysokoenergetycznych szybko osiągają nasycenie. Stosowane w historycznych testach nuklearnych (lata 40.–50.).

Rejestratory optyczne (interferometria VISAR, PDV). Mierzą prędkość powierzchni lub cząstek przez zmianę fazy wiązki laserowej. Stosowane m.in. do pomiaru prędkości detonacji i prędkości Macha na powierzchni tarczy.

Porównanie z aparaturą radioelektroniczną: czas narastania 1 µs odpowiada pasmu ok. 350 kHz — znacznie niżej niż MHz-y typowe dla impulsów nanosekundowych. Blast jest wolniejszy niż impulsy elektroniczne, ale trwa dłużej: kilkadziesiąt do kilkuset milisekund fazy pozytywnej i fazy negatywnej.¹

Faza ujemna (suction phase) i jej znaczenie dla konstrukcji

Model blast często skupia się na fazie dodatniej: gwałtownym wzroście ciśnienia nad atmosferyczne, który trwa czas t+. Ale po fazie dodatniej następuje faza ujemna: ciśnienie spada poniżej atmosferycznego, tworząc podciśnienie (suction). Faza ujemna jest słabsza, ale trwa dłużej.

Dla konstrukcji faza ujemna oznacza odwrócenie kierunku obciążenia. Ściana, która chwilę wcześniej była ściskana do środka, jest teraz ciągniona na zewnątrz. Może to prowadzić do:

• wyrwania elementów połączeniowych (śrub, nitów) zaprojektowanych tylko na ściskanie,
• odłamywania fragmentów betonu po stronie zewnętrznej od fazy dodatniej (spalling), po stronie wewnętrznej od fazy ujemnej,
• drgań rezonansowych, jeżeli częstotliwość fazowa pasuje do własnej budynku.

W modelu pressure-impulse opisanym w artykule o SDOF faza ujemna pojawia się jako odcinek trzeciej fazy pełnego profilu blast. Dla struktur podatnych (ogrodzenia, witryny, zadaszenia) faza ujemna może być decydująca o zniszczeniu, choć nadciśnienie fazy dodatniej jest małe.

Przy Mach stem profil temporalny ciśnienia przy gruncie może różnić się od profilu w swobodnym polu lub od profilu odbicia regularnego. Modelowanie fazy ujemnej jest trudniejsze, bo wymaga pełnego równania stanu powietrza, a proste modele K-B mają tu największe błędy.^1,5

Symulacje numeryczne: ALE, SPH, FEM i ich porównanie

Współczesne narzędzia numeryczne do modelowania blastów i Mach stem:

ALE (Arbitrary Lagrangian-Eulerian). Siatka może poruszać się z materiałem (Lagrangian) lub pozostawać nieruchoma (Eulerian) lub jedno i drugie. Kody: AUTODYN, LS-DYNA z opcją ALE. Dobry do symulacji interakcji fala-konstrukcja (FSI), ale wymaga drobnej siatki przy froncie fali.

SPH (Smoothed Particle Hydrodynamics). Bezsiatkowa metoda: materiał jest reprezentowany przez cząstki. Kod Omanga i współautorów RSPH to sferyczna wersja SPH. Zalety: brak problemów z siatką przy dużych deformacjach. Wady: trudniejsza dokładność przy słabych falach, wyższy koszt obliczeniowy.

FEM (Finite Element Method). Dobry dla konstrukcji (stalowe ramy, płyty betonowe). Kod ABAQUS, LS-DYNA opcja Lagrangian. Kłopot z dużymi przemieszczeniami elementów.

Sprzężone ALE+FEM. Najpopularniejsze podejście: powietrze i produkty wybuchu modelowane ALE, konstrukcja modelowana FEM, sprzężenie przez warunek brzegowy ciśnienie=obciążenie. Wymagane dla dokładnej oceny skutków konstrukcyjnych.

Wartość kalibracyjna: każdy kod przed użyciem do projektowania powinien odtworzyć dane z testów DIPRO, CONWEP/K-B lub ICSBEP-Blast, bo „piękna wizualizacja" bez walidacji nie jest wynikiem inżynierskim.¹

Modelowanie w Polsce: ośrodki i narzędzia badawcze

Badania nad falami uderzeniowymi i blast physics w Polsce są prowadzone w kilku instytucjach:

• WAT (Wojskowa Akademia Techniczna) — Instytut Techniki Uzbrojenia zajmuje się balistyką wewnętrzną i zewnętrzną, eksplozjami, odłamkami i falami uderzeniowymi. Posiadają komorę detonacyjną i aparaturę pomiarową.
• Instytut Techniczny Wojsk Lotniczych (ITWL) — badania odporności konstrukcji lotniczych na uderzenia i obciążenia impulsowe.
• Politechniki — AGH (Kraków), PW (Warszawa), PG (Gdańsk) prowadzą badania z zakresu mechaniki wybuchu, ochrony konstrukcji i numeryki.
• OBRUM (Ośrodek Badawczo-Rozwojowy Urządzeń Mechanicznych) — pojazdy opancerzone, ochrona przed IED, blast protection.

Polska nie prowadzi własnych testów blast jądrowego, ale korzysta z danych NATO i STANAG-ów odnoszących się do chemicznych materiałów wybuchowych. W zakresie konwencjonalnych wybuchów ważne są normy ochrony budynków rządowych i infrastruktury krytycznej na fale blast, wdrożone po zamachach bombowych (np. standardy NATO ATP-55 DOZER dla budynków NATO).⁵

Pojęcia pokrewne: Prandtl-Meyer i kąt Macha

Przy opisie odbicia fali uderzeniowej i warunków pojawienia się Mach stem przydatne są dwa pojęcia z gazodynamiki supersonicznej.

Kąt Macha (μ) to kąt, jaki tworzy czoło małych zaburzeń z kierunkiem przepływu supersonicznego:

sin(μ) = 1/M

Dla M = 1 kąt Macha wynosi 90° — czoło fali jest prostopadłe do kierunku ruchu. Dla M = 2 kąt wynosi 30°, dla M = ∞ maleje do zera. Mach stem pojawia się w geometriach, gdzie lokalne warunki przepływu za falą padającą nie pozwalają na spełnienie warunków odbicia regularnego przy danym kącie padania i danej liczbie Macha fali.

Rozwiązanie Mach triple-point. W punkcie potrójnym trzy fale: padająca (incident), odbita (reflected) i Mach stem spotykają się z warstewką kontaktu (slip line / contact discontinuity). Gaz który przeszedł przez front padający jest inaczej przyspieszony niż gaz który przeszedł przez Mach stem. Warstwa ślizgowa (ang. shear layer) separuje te dwa stany. W warunkach idealnych jest to nieciągłość, w warunkach realnych tworzy się turbulentna warstwa mieszania, widoczna na obrazach Schlieren.

Przejście z odbicia regularnego (RR) do nieregularnego odbicia Mach (MR) zachodzi, gdy kąt padania fali překroczy kąt graniczny. Dokładna wartość tego kąta zależy od liczby Macha fali i modelu gazu. Tzw. „paradoks von Neumanna" polega na tym, że warunki RR→MR z teorii i z eksperymentu nie pokrywają się dokładnie w pewnym zakresie kątów: teoria przewiduje MR, a eksperyment pokazuje jeszcze RR lub odwrotnie. Wyjaśnienie paradoksu wymagało uwzględnienia efektów trójwymiarowych, lepkości i grubości skończonej frontu.

Analogie z mechaniką falową: fale na wodzie

Dla studentów nieznających gazodynamiki użyteczna jest analogia z falami na płytkiej wodzie. W falowym tunelu wodnym można ustawić przeszkodę klinową i obserwować Mach stem na powierzchni wody — tak samo jak w tłumionym gazie. Prędkość rozchodzenia się zaburzeń odpowiada prędkości dźwięku. Punkt potrójny jest widoczny gołym okiem. Wartość edukacyjna jest duża, bo student może ręcznie zmieniać kąt klina i obserwować przejście odbicia regularnego w nieregularne.

W badaniach laboratoryjnych takie tanie, skalowalne analogie często poprzedzają drogie eksperymenty z materiałami wybuchowymi. Stosuje ich wiele europejskich i amerykańskich laboratoriów blast physics jako pierwszy etap walidacji geometrycznej dla kodu symulacyjnego przed testem z pełnowymiarową eksplozją.

Analogia ma granice: fale na wodzie są nieliniowe tylko przy dużej amplitudzie, a gaz doskonały przy wybuchu jest znacznie bardziej nieliniowy. Ponadto płytka woda nie odtwarza fazy ujemnej blast (podciśnienia) ani efektów termicznych kuli ogniowej. Mimo to jako pomoc dydaktyczna do zrozumienia geometrii punktu potrójnego jest niezastąpiona.

Implikacje dla projektowania schronów i ochrony infrastruktury

Znajomość Mach stem ma bezpośrednie przełożenie na projektowanie schronów i wzmocnień infrastruktury. Kilka ważnych reguł:

1. Lokalizacja podziemna. Schron zagłębiony w ziemi jest chroniony przed nadciśnieniem i fragmentami, ale musi wytrzymać fale sejsmiczne wywołane przez falę blast. Fala Macha przy powierzchni ma mniejsze znaczenie dla schronu głębokiego, bo grunt tłumi i dysperguje falę szybko.

2. Geometria ściany frontowej. Ściana prostopadła do źródła wybuchu dostaje ciśnienie odbite Pr (kilkakrotnie większe niż Pso). Ściana nachylona pod kątem zmniejsza efektywne obciążenie. Schron z pochyłymi ścianami czołowymi jest lepszy niż z prostopadłymi, mimo że wygląda na mniej masywny.

3. Krytyczny HOB dla schronu naziemnego. Jeżeli schron jest naziemny, przejście w strefę Macha zwiększa obciążenie nie tylko przez wyższe ciśnienie, ale przez zmieniony profil czasowy. Projekt musi uwzględniać oba parametry: pik i impuls.

4. Polska norma PN-EN 1337. W zakresie ochrony budynków przed wybuchem istnieje seria norm europejskich (EN 13123, EN 13124 dla okien, EN 1337 dla urządzeń łożyskowych) i normy NATO. WAT i COBIK (Centrum Ochrony Biernej i Kryptografii) współpracują przy wdrażaniu standardów ochrony infrastruktury krytycznej w Polsce, uwzględniających parametry blast.

Typowe błędy interpretacyjne

Pierwszy błąd to utożsamienie Mach stem z każdą falą odbitą. Odbicie regularne nie jest jeszcze frontem Macha. Potrzebny jest układ trzech fal i punkt potrójny nad powierzchnią.

Drugi błąd to traktowanie wysokości wybuchu jako jednej liczby oderwanej od progu. HOB ma sens tylko razem z pytaniem: jaki próg ciśnienia, jaka wysokość celu i jaki typ skutku? Zmiana kryterium zmienia odpowiedź.

Trzeci błąd to czytanie Pmax bez impulsu. Jeżeli dwa punkty mają podobny pik, ale różny czas dodatniej fazy, konstrukcja może zareagować inaczej. Przy Mach stem taka różnica jest szczególnie ważna, bo pionowy profil ciśnienia i czas pomiaru są powiązane.

Czwarty błąd to brak rozdzielenia geometrii i modelu. To, że kalkulator rysuje maksimum promienia dla progu, nie znaczy, że zna lokalną topografię, zabudowę, atmosferę ani wielokrotne odbicia. Model HOB jest szkieletem fizycznym, nie mapą realnego miasta.

Piąty błąd to ignorowanie prawa skalowania przy porównywaniu ładunków o różnej masie. Zestawienie HOB w metrach dla 1 kg i dla 1 t TNT jest bez sensu bez normalizacji przez W^(1/3). Dopiero skalowana wysokość h/W^(1/3) jest sensownym parametrem porównawczym między różnymi kalibrami.

Szósty błąd to mylenie akustycznego analogu (demonstracja z mikrofonem i impulsem dźwiękowym) z rzeczywistym pomiarem blast. Analogia jest użyteczna dydaktycznie do odczytywania geometrii, ale nic nie mówi o nieliniowości, ciśnieniu odbitym ani fazie ujemnej. Oba zadania — nauczanie geometrii odbicia i pomiar nadciśnienia — wymagają różnych narzędzi i nie wolno ich mieszać. Student, który zrozumie tę granicę, jest gotowy do pracy z prawdziwymi danymi pomiarowymi z testów empirycznych.

Dodatkowe materiały multimedialne

Powiązane kalkulatory i narzędzia

• Height-of-burst i Mach stem - porównuje falę swobodną, powierzchniową i dydaktyczny model HOB/Mach-stem.
• Podmuch - dostarcza parametrów fali K-B/CONWEP potrzebnych jako wejście do prostego modelu HOB.
• Pressure-Impulse - pokazuje, co dzieje się dalej, gdy fala odbita obciąża konstrukcję.
• Odporność płyty - liczy modalną odpowiedź płyty na zadany przebieg ciśnienia p(t).

Ćwiczenia praktyczne

Pierwsze ćwiczenie powinno polegać na rekonstrukcji geometrii z czasów przybycia. Należy:

1. ustawić w symulacji rząd czujników na gruncie i drugi rząd na kilku wysokościach,
2. wygenerować syntetyczne czasy przybycia fali bezpośredniej i odbitej,
3. narysować krzywe równego czasu,
4. wyznaczyć miejsce, w którym dwa fronty zaczynają zlewać się w front Macha,
5. opisać, dlaczego pojedynczy czujnik nie wystarcza do rozpoznania reżimu odbicia.

Drugie ćwiczenie powinno dotyczyć danych ciśnieniowych. Należy:

1. otrzymać tabelę Pmax, t+ i I dla kilku czujników,
2. porównać punkty na gruncie i punkty na wysokości,
3. wskazać, gdzie pionowy profil ciśnienia ma maksimum poza gruntem,
4. wyjaśnić, dlaczego maksymalizacja nadciśnienia na powierzchni nie jest tym samym, co maksymalizacja obciążenia wyższego obiektu,
5. zapisać niepewności związane z czasem próbkowania i pasmem czujników.

Trzecie ćwiczenie powinno porównać model dydaktyczny z danymi syntetycznymi. Należy:

1. użyć kalkulatora HOB/Mach stem dla kilku wysokości,
2. odczytać reżim, mnożnik i promień dla progu,
3. porównać wynik z przygotowaną tabelą „pomiarową”,
4. wskazać, czy rozbieżność dotyczy piku, czasu przybycia czy impulsu,
5. przygotować krótką notatkę o tym, czego lokalny model nie obejmuje.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Wysokość detonacji a fala podmuchowa
• Od nadciśnienia do impulsu: praktyczne modele Kingery-Bulmash, Brode i CONWEP
• Pressure-Impulse i model SDOF: dlaczego konstrukcje nie reagują tylko na pik ciśnienia
• Fala odbita i obciążenie konstrukcji
• Oscyloskop w pomiarach impulsów nanosekundowych