Masa krytyczna i jej parametry geometryczne

Streszczenie

Najłatwiej zrozumieć masę krytyczną, wyobrażając sobie mgłę deszczu spadającą przez las. Jeżeli drzewa stoją rzadko, wiele kropel przeleci między nimi i zniknie poza lasem. Jeżeli jednak zagęścić las albo otoczyć go ścianą, która zawraca część kropel do środka, coraz większa część ruchu pozostanie wewnątrz układu. Z neutronami jest podobnie: masa krytyczna nie jest magiczną liczbą wpisaną raz na zawsze w dany izotop, lecz granicą, przy której układ zaczyna zatrzymywać dość neutronów, by reakcja łańcuchowa mogła się sama podtrzymywać.^1,2

Dlatego pytanie „ile potrzeba materiału?” jest zawsze niepełne. Trzeba od razu dopytać: w jakim kształcie, o jakiej gęstości, z jakim reflektorem i dla jakiego tempa reakcji. Kula ma przewagę nad cylindrem, bo przy tej samej objętości ma najmniejszą powierzchnię, a więc najmniej „miejsca ucieczki” dla neutronów. Z kolei kompresja obniża masę krytyczną tak silnie, że w konstrukcjach implozyjnych kilka kilogramów materiału może po sprężeniu zachowywać się jak wielokrotność masy krytycznej.^1,3

Rozszerzenie tematu

Masa krytyczna to minimalna masa bryły materiału rozszczepialnego, dla której straty neutronów są na tyle małe, że reakcja łańcuchowa może się samopodtrzymywać.¹ W języku artykułu poprzedniego oznacza to po prostu stan, w którym efektywny współczynnik mnożenia osiąga wartość bliską $k_{ef} = 1$. Poniżej tej granicy układ jest podkrytyczny i gaśnie, powyżej niej staje się nadkrytyczny. Sama idea jest prosta, ale liczba odpowiadająca masie krytycznej zależy od kilku parametrów jednocześnie.

Pierwszy z nich to geometria. Przy ustalonej objętości bryła kulista ma najmniejszą powierzchnię, a więc najmniejsze prawdopodobieństwo ucieczki neutronów przez granicę układu. To dlatego wartości masy krytycznej podaje się zazwyczaj dla kuli, o ile nie zaznaczono inaczej.¹ Zależność ta nie wynika z żadnej „preferencji natury” dla sfery, lecz z relacji powierzchni do objętości. Im większa powierzchnia przypada na jednostkę objętości, tym łatwiej neutron opuści układ, zanim jeszcze spowoduje kolejne rozszczepienie.

Ten punkt dobrze uzupełnia podręcznik Słowińskiego, który obok samej masy podaje też promień krytyczny dla jednorodnej sferycznej strefy aktywnej. Dla czystych izotopów rozszczepialnych przytacza on wartości rzędu około 8,5 cm i 48 kg dla uranu-235, około 6 cm i 17 kg dla plutonu-239 oraz około 6 cm i 16 kg dla uranu-233.² Te liczby trzeba czytać ostrożnie: są to wartości modelowe dla czystych, kulistych układów, a nie uniwersalne stałe fizyczne. Mimo to dobrze pokazują skalę zjawiska i tłumaczą, dlaczego geometria jest tak ważna.

Drugi parametr to gęstość. Gdenarz wyprowadza to bardzo jasno: dla danego kształtu prawdopodobieństwo ucieczki neutronu zależy od rozmiaru bryły wyrażonego w średnich drogach swobodnych, a ponieważ średnia droga swobodna skaluje się odwrotnie z gęstością, masa krytyczna jest w przybliżeniu proporcjonalna do $1/\rho^2$.¹ To jedna z najważniejszych zależności w całej fizyce pierwszych bomb:

$$m_{kr} \sim \frac{1}{\rho^2}$$

Nie chodzi tu o dokładny wzór projektowy dla konkretnej konstrukcji, lecz o prawo skali. Jeżeli podwoić gęstość materiału, masa krytyczna spada w przybliżeniu czterokrotnie. Dlatego w bombie implozyjnej nie trzeba zaczynać od wielkiej, „gołej” masy krytycznej. Wystarczy przygotować układ lekko podkrytyczny, a następnie bardzo szybko zwiększyć jego gęstość przez sprężenie. Podwojenie gęstości rdzenia, który początkowo jest lekko podkrytyczny, czyni z niego niemal cztery masy krytyczne i daje wystarczające wprowadzenie reaktywności dla bomby.³

To prowadzi bezpośrednio do pytania, dlaczego metoda implozyjna była tak wielkim przełomem. W konstrukcji działowej zwiększenie reaktywności osiąga się przede wszystkim przez szybkie połączenie dwóch części materiału rozszczepialnego. W konstrukcji implozyjnej osiąga się to inaczej: przez gwałtowne podniesienie gęstości tego samego rdzenia. Taki układ jest trudniejszy inżynieryjnie, ale znacznie lepiej wykorzystuje materiał i pozwala zejść do mniejszych mas fizycznych przy dużo wyższej nadkrytyczności.³

Trzeci ważny parametr to otoczenie rdzenia, czyli reflektor neutronów. Jeżeli wokół materiału rozszczepialnego znajduje się warstwa substancji, w której neutrony łatwo się rozpraszają, ale nie są chętnie pochłaniane, część neutronów wróci do rdzenia zamiast zostać utracona. Takie otoczenie zmniejsza wymaganą masę krytyczną.¹ Gdenarz pokazuje to na przykładzie kulistej bryły uranu-235: wzrost grubości reflektora z berylu, naturalnego uranu lub żelaza istotnie obniża krytyczną masę układu.¹

W praktyce reflektor nie jest jedynie dodatkiem „neutronowym”. W urządzeniach militarnych zwykle łączy się go z rolą tampra, czyli ciężkiej warstwy bezwładnej, która przez bardzo krótki czas opóźnia rozbieganie się rozgrzanego rdzenia. Znaczenie reflektora bywa przy tym przeceniane, bo neutron, który ucieknie z rdzenia i wróci dopiero po dłuższej drodze, wnosi mniej niż neutron pozostający od początku wewnątrz masy krytycznej. Mimo to dobry reflektor nadal wyraźnie poprawia bilans neutronów i obniża wymaganą masę materiału.³

Zależność od reflektora tłumaczy też, dlaczego eksperymenty krytycznościowe były tak niebezpieczne. Wypadki związane z Demon Core nie polegały na „dodawaniu nowego paliwa”, lecz na takim ustawieniu otaczających rdzeń elementów odbijających neutrony, że układ zbliżał się do krytyczności albo ją przekraczał. To bardzo dobry historyczny przykład na to, że masa krytyczna nie jest własnością samej porcji plutonu, ale całej konfiguracji geometrycznej i materiałowej wokół niego.⁴

Warto również odróżniać masę krytyczną reaktorową od „użytecznej” masy krytycznej dla bomby. W reaktorze można dopuścić moderator i wolniejszy przebieg reakcji, bo celem jest stabilne $k_{ef} = 1$. W ładunku jądrowym moderacja jest niepożądana, bo spowalnia rozwój reakcji. Gdenarz pisze wprost, że choć da się uzyskać krytyczność dzięki neutronom termicznym i moderatorowi, to w bombie taki zabieg byłby szkodliwy, ponieważ znacząco zwolniłby przebieg rozszczepienia.¹ Z militarnego punktu widzenia liczy się więc nie dowolna masa krytyczna, lecz taka konfiguracja, która daje bardzo szybki przyrost mocy.

Trzeba też pamiętać, że wartości liczbowe podawane w podręcznikach są najczęściej obliczeniami modelowymi. Słowiński wyraźnie zaznacza, że proste analityczne podejście ma przede wszystkim wartość poglądową, a praktyczne liczby użyteczne projektowo uzyskuje się dopiero z rozbudowanych modeli numerycznych.² To ważna uwaga metodologiczna. W projekcie takim jak ten można i trzeba podawać orientacyjne liczby, ale nie wolno ich przedstawiać tak, jakby były bezwarunkowymi stałymi przyrody.

Z tego wynika dobra intuicja końcowa. Masa krytyczna nie jest „progiem ilościowym” w prostym sensie: że poniżej nic się nie dzieje, a powyżej wszystko działa. Jest raczej punktem, w którym geometria, gęstość, skład izotopowy i otoczenie razem zaczynają dawać dodatni bilans neutronowy. Dlatego dwa układy z identyczną masą tego samego materiału mogą zachowywać się zupełnie inaczej. Jeden pozostanie bezpiecznie podkrytyczny, a drugi po dodaniu reflektora albo po kompresji przejdzie do stanu silnie nadkrytycznego.

To właśnie dlatego masa krytyczna jest jednym z najważniejszych pojęć łączących fizykę z inżynierią. Bez niej nie da się rozumieć ani pracy reaktora, ani sensu soczewek wybuchowych, ani tego, dlaczego rdzeń Fat Mana mógł mieć fizycznie około 6,2 kg, a mimo to po kompresji wejść w obszar wielokrotnej nadkrytyczności.³

Matematyczne podejście do masy krytycznej

Formalne wyprowadzenie masy krytycznej opiera się na teorii dyfuzji neutronów i warunku krytyczności, który jest spełniony, gdy kubatura $B^2$ spełnia:

$$k_\infty = 1 + M^2 B^2$$

gdzie $k_\infty$ jest mnożnikiem nieskończonego układu, $M^2$ kwadratem długości migracji i $B^2$ beklinem geometrycznym. Dla sfery napromieniowanej nago (bez reflektora):

$$B^2 = \left(\frac{\pi}{R + d}\right)^2$$

gdzie $R$ to promień sfery, a $d \approx 2D$ to ekstrapolacja dyfuzyjna (~2 cm dla U-235). Warunek krytyczności $k_\infty = 1 + M^2 B^2$ daje:

$$R_c = \frac{\pi M}{\sqrt{k_\infty - 1}} - d$$

a masa krytyczna:

$$m_c = \frac{4}{3}\pi R_c^3 \cdot \rho$$

Kluczowym parametrem jest $M^2 = L^2 + \tau$, gdzie $L$ to długość dyfuzji neutronów termicznych, a $\tau$ wiek Fermiego. Dla czystego $^{235}$U w gęstości normalnej ($\rho = 18{,}7\,\mathrm{g/cm^3}$) wartości przybliżone:

• $k_\infty \approx 2{,}5$ (dla 100% U-235 z neutronami prędkimi)
• $M \approx 3{,}5\,\mathrm{cm}$
• $R_c \approx 8{,}5\,\mathrm{cm}$
• $m_c \approx 48\,\mathrm{kg}$

Dla $^{239}$Pu ($\rho = 19{,}6\,\mathrm{g/cm^3}$): $k_\infty \approx 2{,}9$, $R_c \approx 6{,}2\,\mathrm{cm}$, $m_c \approx 10{-}17\,\mathrm{kg}$ (zależnie od modelu).⁵

Zależność $m_c \sim 1/\rho^2$ wynika z prostej analizy skalowania: jeśli $\rho \to \lambda\rho$, to średnia droga swobodna $\lambda_{mfp} \to \lambda_{mfp}/\lambda$, a więc $M \to M/\lambda$ i $R_c \to R_c/\lambda$. Masa $m_c \propto R_c^3 \cdot \rho \propto (1/\lambda)^3 \cdot \lambda = 1/\lambda^2$. Stąd dwukrotne zwiększenie gęstości daje czterokrotnie mniejszą masę krytyczną.

Dane dla różnych izotopów rozszczepialnych

Poniższa tabela zbiera przybliżone masy krytyczne dla nagich kul (bez reflektora) z różnych materiałów, przy gęstości nominalnej i normalnej frakcji rozszczepialnej:

Izotop	Gęstość [g/cm³]	$m_c$ naga kula [kg]	$R_c$ [cm]	Uwagi
$^{235}$U (100%)	18,7	48	8,7
$^{235}$U (90% HEU)	~18	~55	~9	typowy HEU
$^{239}$Pu (100%)	19,6	10–17	6,2	różne wg modelu
$^{233}$U (100%)	18,6	~16	~6,5
$^{241}$Pu	19,8	~12	~6,5
$^{252}$Cf	15,1	~2,7 kg	~4	choć Cf jest rzadki
$^{237}$Np	20,2	~60	~8,7	z reflektorem woda

Warto zwrócić uwagę na kilka obserwacji:

Pluton-239 ma masę krytyczną ok. 3–4 razy mniejszą niż uran-235. To właśnie jeden z powodów (obok efektywności produkcji plutonu w reaktorach) dla których pluton stał się materiałem preferowanym w nowszych konstrukcjach broni.

Uran-233 (z cyklu torowego) ma masę krytyczną bliską plutonu-239. Indie rozważają u-233 jako materiał dla swojego programu jądrowego opartego na torze.

Californium-252 ma bardzo małą masę krytyczną (~2,7 kg), ale produkowane są tylko gramy rocznie (reaktory badawcze). Stosowany jako kompaktowe źródło neutronów w miernikach grubości i aktywacji analitycznej.⁵

Wpływ reflektora — dane ilościowe

Reflektor neutronów dramatycznie redukuje masę krytyczną. Oto porównanie dla różnych grubości i materiałów reflektora dla kuli $^{239}$Pu:

Reflektor	Grubość [cm]	$m_c$ [kg]	Redukcja
Brak (naga kula)	0	~10	—
Naturalny uran	5	~5,6	~44%
Naturalny uran	10	~4,5	~55%
Naturalny uran	20	~4,0	~60%
Berill Be	10	~4,5	~55%
Woda H₂O	10	~8	~20%

Naturalny uran (U-238) jest doskonałym reflektorem: duże gęstość, duży przekrój czynny na rozpraszanie sprężyste neutronów prędkich, mały pochłanianie. Dlatego u-238 był preferowanym materiaem tamperowym w pierwszych bombach implozyjnych. Beryl jest jeszcze lepszym reflektorem neutronów prędkich (mały A, małe pochłanianie), ale jest drogie i toksyczne.

Gruba warstwa reflektora pozwala zredukować masę plutonu do ~4 kg — wartości bliskiej historycznym rdzeniom implozyjnym. „Demon Core" — plutoniowa kula używana w wypadkach Los Alamos w 1945–1946 — ważyła 6,2 kg i była wystarczająco blisko krytyczności, że reflektory berylu lub elementy z karbidu wolframu zbliżały ją do, a czasem przekraczały, granicę krytyczną.⁵

Eksperymenty krytycznościowe — Godiva i Blue Water

Historia wyznaczania parametrów krytyczności jest historią niebezpiecznych i starannie kontrolowanych eksperymentów krytycznych, prowadzonych przez dziesięciolecia w Los Alamos i innych laboratoriach jądrowych.

Godiva I, II, III (1955–1967): Seria gołych kul uranu-235 o wysokim wzbogaceniu (~93% HEU), używanych do pulsowych pomiarów krytyczności w Los Alamos. Nazwa pochodzi od Lady Godiva — nagiej (gołej) kuli bez reflektora. Godiva umożliwiała precyzyjne pomiary geometrycznego progu krytyczności, pomiary przekrojów czynnych neutronów prędkich przez transmisję i aktywację.

Jezebel: Gołe kule plutonu-239 używane w LASL (Los Alamos Scientific Laboratory) do pomiarów krytyczności. Jezebel i Godiva razem stanowiły podstawowe narzędzia doświadczalnego wyznaczania stałych jądrowych niezbędnych do projektowania broni.

VENUS (B4C Control Rod Experiment): europejski reaktor krytycznościowy w SCK·CEN w Belgii, służący do weryfikacji danych bibliotek nuklearnych.

JMTR, REBUS, Godiva-IV: Eksperymenty krytycznościowe prowadzone przez OECD/NEA i IAEA w ramach programu International Criticality Safety Benchmark Evaluation Program (ICSBEP), który zbiera i weryfikuje dane z historycznych eksperymentów. Baza danych ICSBEP (dostępna publicznie) zawiera tysiące zestawów danych dla różnych geometrii, materiałów i konfiguracji reflektorów.⁵

Wypadki krytyczne poza Demon Core

Historia jądrowa obfituje w wypadki krytyczne, z których każdy ilustruje inny aspekt fizyki masy krytycznej:

Wypadek w Wood River Junction (1964): Pracownik zakładu przeróbki uranu w Rhode Island wlał roztwór uranu wzbogaconego do nieprawidłowego naczynia (okrągłe naczynie o geometrii sprzyjającej krytyczności, z wodą jako moderatorem). Nastąpił błysk gamma i impuls neutronowy — pracownik zginął w ciągu 49 godzin po pochłonięciu dawki ~4400 rad. Wypadek ujawnił, że krytyczność może nastąpić w roztworach wodnych przy znacznie mniejszych masach niż w układach stałych.

Wypadek Sarov (Arzamas-16, 1997): Rosyjski technik pracujący z kulą HEU popełnił błąd konfiguracyjny, zbliżając dwie części uranu. Impuls krytyczny trwał ułamek sekundy. Technik zginął po ~2 tygodniach na skutek ostrego poparzenia radiacyjnego. Był to jeden z ostatnich znanych wypadków krytycznych w badanym arsenale jądrowym.

Wypadek SL-1 (1961): Opisany szczegółowo w artykule o neutronach opóźnionych — rdzeń reaktora SL-1 w Idaho osiągnął promptową nadkrytyczność po ręcznym wyciągnięciu pręta kontrolnego, co spowodowało eksplozję parową i trzy ofiary śmiertelne.

Wypadek Tokaimura (1999): Nadkrytyczność w roztworze azotanu uranu — impuls nie był gwałtowny (reakcja trwała godziny), ale dawki śmiertelne pochłonęli dwaj pracownicy.

Wspólny mianownik wszystkich wypadków: brak zrozumienia lub niedostateczne uwzględnienie wszystkich parametrów geometrycznych, moderatorowych i gęstościowych tworzących masę krytyczną. Fizyka jest prosta; przestrzeganie jej wymaga systemu bezpieczeństwa.⁵

Bezpieczeństwo krytycznościowe w przemyśle

Bezpieczeństwo krytycznościowe (Criticality Safety) jest odrębną subdyscypliną inżynierii jądrowej. Jej celem jest zapewnienie, że żaden układ materiałów jądrowych w zakładzie przemysłowym (wzbogacalnia, fabryka paliwa, zakład przeróbki) nigdy nie osiągnie stanu krytycznego w sposób niezamierzony.

Główne parametry bezpieczeństwa krytycznościowego:

• Masa podkrytyczna — ograniczenie masy w pojedynczym naczyniu lub obszarze
• Geometria podkrytyczna — użycie cienkich, wydłużonych naczyń (cylindrycznych), gdzie geometria jest niekorzystna dla krytyczności
• Rozcieńczenie — mieszanie wzbogaconego uranu z materiałami niepodtrzymującymi krytyczności
• Brak moderatora — operacje z suchymi, niezawodnionymi materiałami
• Kontrola interakcji — oddzielenie przestrzenne między kontenerami zawierającymi materiały rozszczepalne

Norma ANSI/ANS-8.1 (Criticality Safety Criteria) i przepisy NRC (Nuclear Regulatory Commission) wymagają co najmniej dwóch niezależnych marginesów podkrytyczności. Zakłady paliwowe (np. Global Nuclear Fuel w USA) muszą posiadać licencje krytycznościowe obejmujące analizy bezpieczeństwa dla każdego etapu procesu technologicznego.

Ważnym narzędziem jest pojęcie pojedynczej jednostki podkrytycznej (sub-critical unit) i interakcji między jednostkami (array criticality). Dwie paczki HEU, z których każda jest indywidualnie podkrytyczna, mogą razem stanowić układ krytyczny — bo neutrony z jednej paczki wchodzą do drugiej.⁵

Masa krytyczna a proliferacja jądrowa

Masa krytyczna jest centralnym pojęciem w analizie proliferacyjnej. IAEA definiuje istotną ilość (Significant Quantity, SQ) jako minimalną masę materiału jądrowego, z której — w teorii — można zbudować urządzenie jądrowe:

• 8 kg plutonu (90% Pu-239)
• 25 kg HEU (90% U-235)
• 75 kg LEU (20% U-235)
• 20 kg U-233

Wartości SQ są celowo konserwatywne (wyższe niż faktyczne masy krytyczne z reflektorem), by zapewnić margines bezpieczeństwa dla inspekcji. IAEA wymaga raportowania wszelkich przekroczeń lub zbliżeń do progu SQ.

Z perspektywy proliferacyjnej kluczowe jest:

Po pierwsze, rosnące ryzyko minuaturyzacji z postępem technologicznym. Dokładność obliczeniowa i lepsze materiały reflektorowe mogą obniżyć faktyczną masę używaną w urządzeniach. Dlatego Stockpile Stewardship (program bezpieczeństwa arsenału USA) nie jest tylko technologią wojskową — dane tam generowane pozwalają na lepsze szacunki minimów krytycznych, co jest istotne dla weryfikacji traktatów.

Po drugie, bezpieczeństwo fizyczne HEU i plutonu jest przedmiotem Konwencji o Ochronie Fizycznej Materiałów Jądrowych (CPPNM, 1980) i jej zmian. Summit bezpieczeństwa jądrowego (Nuclear Security Summits 2010–2016) przyczyniły się do konwersji ~93 reaktorów badawczych z HEU na LEU oraz do oczyszczenia wielu stron z niewielkich ilości HEU.⁵

Historia wyznaczania masy krytycznej — Los Alamos 1943–1945

Historia eksperymentalnego wyznaczania masy krytycznej w Los Alamos jest jednym z najmroczniejszych rozdziałów Projektu Manhattan. Fizycy pracowali z substancjami, których krytyczne parametry były znane tylko w przybliżeniu, korzystając z metod, które dziś byłyby niedopuszczalne.

Robert Serber, autor Los Alamos Primer (wewnętrznego podręcznika dla nowych pracowników Los Alamos z 1943 roku, odtajnionego w 1992), opisał fizyczną logikę masy krytycznej w sposób, który nadal jest wzorcową prezentacją: zależy ona od przekroju czynnego na rozszczepienie, od przekroju na rozpraszanie i od geometrii układu. Serber wyprowadził szacunkowe formuły, ale podkreślał, że dokładne wyniki wymagają eksperymentów.

Eksperymenty te prowadzono przede wszystkim w ramach „wydziału krytycznych zestawów" (Critical Assemblies) kierowanego przez Otto Frisch. Frisch (jeden z odkrywców mechanizmu rozszczepienia) używał techniki „szybkiego zbliżania" (fast-join experiments): dwie części HEU zbliżano do siebie ręcznie lub mechanicznie i mierzono wzrost aktywności neutronów, czerpiąc z chwilowych wartości $k_\mathrm{eff}$.

W sierpniu 1945 roku, kilka dni po próbach Trinity i bombardowaniu Japonii, Harry Daglian pracował samotnie w nocy z tym samym rdzeniem plutonowym (Demon Core). Popełnił błąd, omówiony w artykule Demon Core. W maju 1946 roku Louis Slotin zginął w wyniku analogicznego wypadku podczas demonstracji.

Po tych wypadkach całkowicie zmieniono metodologię eksperymentów krytycznych. Zbliżanie ręczne zostało zakazane; zastąpiono je zdalnymi systemami mechanicznymi z automatycznym rozdzieleniem przy wykryciu nadkrytyczności. Dziś eksperymenty krytyczne w placówkach takich jak CRITICALITY SAFETY EXPERIMENT FACILITY (CSEF) w ORNL lub NCERC (National Criticality Experiments Research Center) w Nevada National Security Site (NNSS) są prowadzone wyłącznie zdalnie, za osłonami, z wieloma redundantnymi systemami bezpieczeństwa.⁵

Krytyczność pulsacyjna i reaktory TRIGA

Szczególną kategorią urządzeń operujących z masą krytyczną są reaktory impulsacyjne i urządzenia do szybkich pomiarów krytyczności. Reaktory TRIGA (Training, Research, Isotopes, General Atomics), projektowane przez Edwarda Tellersa i Freemana Dysona pod koniec lat 50., są znane ze zdolności do kontrolowanych impulsów mocy.

W reaktorze TRIGA paliwo zawiera domieszek wodorku cyrkonu (ZrH), który zapewnia silnie ujemny współczynnik temperatury paliwa. Gdy temperatura rdzenia wzrasta (podczas impulsu), wodór w ZrH ogrzewa się, co obniża $k$ — efektywnie wyłączając reaktor termicznie bez potrzeby sterowania prętami.

Dzięki temu TRIGA można pulsować do chwilowej mocy ~1000–2000 MW przez czas ~10–20 ms, a reaktor samoistnie powraca do niskich mocy bez żadnej interwencji operatora. Fizycznie rzecz biorąc, rdzeń osiąga nadkrytyczność (wyjście powyżej masy krytycznej dla chwilowych warunków), ale ujemny współczynnik temperatury działa jak własne sprzężenie zwrotne.

To klasyczny przykład, jak parametry masy krytycznej (gęstość, skład, temperatura) można dobierać tak, by uzyskać pożądaną kinematykę — w tym przypadku niemal niesterowalną samoistną stabilizację przez ujemny współczynnik temperatury.

Reaktory TRIGA były budowane w ponad 40 krajach na całym świecie, w tym w Polsce: TRIGA Mark II na Politechnice Warszawskiej (reaktor MARIA w Świerku to oddzielny typ). Dziś wymagania IAEA dotyczące konwersji reaktorów badawczych z HEU na LEU dotyczyły m.in. TRIGA — paliwo HEU w wielu reaktorach TRIGA jest stopniowo zastępowane nowym paliwem LEU o wyższym wzbogaceniu masowym (~20%), które przy odpowiedniej geometrii pozwala utrzymać krytyczność na neutronach termicznych.⁵

Masa krytyczna w reaktorach energetycznych

W kontekście reaktorów energetycznych koncepcja masy krytycznej ulega modyfikacji: zamiast masy krytycznej pojedynczego izotopu rozszczepialnego, mówi się o krytyczności układu złożonego z paliwa, moderatora, chłodziwa, reflektora i materiałów strukturalnych.

Dla reaktora PWR z paliwem UO₂ (~3-5% U-235):

• Masa całego rdzenia: ~100–120 ton UO₂
• Masa U-235: ~2,5–4,5 tony
• Masa Pu-239 po wypale: ~200–300 kg
• Całkowity $k_\infty$ dla świeżego paliwa: ~1,5–1,8
• $k_\mathrm{eff}$ pracy: ~1,000 (z regulacją przez pręty)

Rdzeń reaktora jest projektowany tak, by $k_\mathrm{eff} = 1$ przez cały cykl paliwowy, mimo że $k_\infty$ zmienia się znacznie. Regulacja odbywa się przez:

• Pręty kontrolne (B₄C lub Ag-In-Cd): pochłaniacze neutronów wchodzące do rdzenia
• Bor rozpuszczony w chłodziwie (PWR): regulacja stężenia kwasu borowego
• Palniki pochłaniające (BP, burnable poisons): Gd₂O₃ lub B₄C wbudowane w paliwo

Masa krytyczna w reaktorach nie jest zatem stałą, lecz parametrem, którym operatorzy aktywnie zarządzają przez cały cykl paliwowy. Różnica między reaktorem a bronią polega m.in. na celowej pracy w zakresie opóźnionej krytyczności i przy $k_\mathrm{eff}$ utrzymanym bliskim jedności — zamiast prędkiego przejścia przez wielokrotną nadkrytyczność, jak w bombie.

Wyjątkiem jest rozruch reaktora i ładowanie paliwa. W trakcie wprowadzania nowych palikwodnych prętów z paliwa do rdzenia, $k_\mathrm{eff}$ zmienia się i operatorzy muszą precyzyjnie śledzić stany krytyczne. Procedury „kontrolowanej krytyczności" podczas rozruchu (approach to criticality) są wśród najdokładniej regulowanych etapów eksploatacji reaktora.⁵

Masa krytyczna dla roztworów wodnych

Szczególnie niebezpieczną kategorią są układy krytyczne w roztworach wodnych (aqueous solutions). Woda jako moderator dramatycznie obniża masę krytyczną uranu wzbogaconego:

Wzbogacenie	Masa krytyczna w H₂O [kg U]	Masa jednostki geometrycznej
90% HEU	0,8–1,0	mała butelka (cylindryczna)
20%	~8	kilka litrów
5%	>100	duża objętość

Obecność wody jako moderatora obniża masę krytyczną HEU do poziomu poniżej 1 kg uranu całkowitego — co czyni nawet małe rozlewy lub kontaminację wodą potencjalnie niebezpiecznymi.

To tłumaczy, dlaczego zakłady wzbogacania uranu i fabryki elementów paliwowych muszą:

1. Operować z surowcami w postaci ciał stałych lub w bezwodnych rozpuszczalnikach organicznych
2. Utrzymywać geometrię naczyń cylindryczną (niekorzystna dla krytyczności) zamiast kulistej
3. Skrupulatnie kontrolować stężenia roztworów i ich ilości
4. Zapobiegać gromadzeniu się wody w nieplanowanych miejscach

Wypadek Wood River Junction (1964) i Tokaimura (1999) obydwa dotyczyły nadkrytyczności w roztworach wodnych. To właśnie dlatego bezpieczeństwo krytycznościowe dla roztworów jest osobną, szczególnie rygorystyczną subdyscypliną.⁵

Masa krytyczna a koncepcja bomb

Różne koncepcje bomb jądrowych mają radykalnie różne wymagania dotyczące masy krytycznej:

Bomba działowa (Gun-type, Little Boy): Dwie podkrytyczne masy U-235 (jedna cylindryczna, jedna tarcza) łączone wybuchem konwencjonalnym. Łączna masa 64 kg U-235, z czego ~0,85 g uległo rozszczepieniu w wybuchu. Sprawność zaledwie ~1,4%. Nie nadaje się dla plutonu (zbyt duże ryzyko preinicjacji przez Pu-240).

Bomba implozyjna (Fat Man): Rdzeń z ~6,2 kg Pu-239 + podkrytyczna konfiguracja U-238/Pu, skompresowany do ~2× gęstości normalnej. Sprawność ~17%. Dzięki kompresji osiągnięta masa efektywna (przeliczona na gęstość krytyczną) to odpowiednik ~40–50 mas krytycznych kuli Pu przy gęstości normalnej.

Boosted fission: Rdzeń Pu-239 z rdzeniem D-T, mniejszy niż w Fat Man. Kompresja i produkcja neutronów z fuzji pozwalają na zmniejszenie masy fizycznej przy zachowaniu lub zwiększeniu sprawności.

Warhead termojądrowy (staged weapon, Teller-Ulam): Rdzeń rozszczepialny (primary) jest miniaturowy — może ważyć 1–3 kg Pu przy wysokim stopniu kompresji. Jego zadaniem jest wytworzenie impulsu rentgenowskiego do implozyji drugiego etapu (secondary, fuzja).⁵

Podsumowanie i znaczenie interdyscyplinarne

Masa krytyczna jest pojęciem na przecięciu fizyki jądrowej, hydrodynamiki, inżynierii materiałowej i fizyki transportu neutronów. Jej wartość dla danego układu zależy od co najmniej pięciu niezależnych parametrów: materiału (izotopu), gęstości, geometrii, reflektora i czystości materiału. Zmiana któregokolwiek z tych parametrów może przenieść układ z bezpiecznej podkrytyczności do niebezpiecznej nadkrytyczności.

To właśnie ta złożoność czyni wypadki krytycznościowe tak groźnymi: operatorzy nieświadomi jednego parametru — na przykład że naczynie okrągłe jest bardziej krytyczne niż wydłużone przy tej samej objętości, lub że zwykła woda pitna jest skutecznym moderatorem — mogą nieświadomie stworzyć konfigurację spełniającą wszystkie warunki krytyczności i tym samym wywołać niekontrolowany impuls neutronów.

Jednocześnie ta sama złożoność daje inżynierom ogromne możliwości kontroli: przez zmianę geometrii, gęstości, reflektora i moderatora można zarówno projektować reaktory, jak i budować urządzenia, w których krytyczność jest osiągana na żądanie i pod kontrolą. Masa krytyczna jest więc nie tylko granicą, za którą kryje się destrukcja — jest też parametrem projektowym, który decyduje o bezpieczeństwie i efektywności całej energetyki jądrowej.^1,2,3,5

Masa krytyczna w kontekście małych reaktorów modułowych (SMR)

W ostatnich latach koncepcja małych reaktorów modułowych (SMR — Small Modular Reactors) wymusza ponowne przemyślenie zagadnień masy krytycznej. Reaktory o mocy 50–300 MWe muszą mieć znacznie mniejsze rdzenie niż konwencjonalne PWR (>1000 MWe), co wymaga wyższego wzbogacenia paliwa lub zastosowania specjalnych materiałów, aby utrzymać krytyczność.

Kilka kluczowych implikacji:

HALEU (High-Assay Low-Enriched Uranium): Paliwo o wzbogaceniu 5–20% U-235, potrzebne dla wielu projektów SMR (np. reaktor BWRX-300 GE-Hitachi, NuScale, Kairos FHR). Przy wyższym wzbogaceniu masa krytyczna maleje, ale wzrasta wrażliwość układu na błędy operacyjne i zagrożenia proliferacyjne. HALEU jest objęty szczególnym reżimem kontrolnym w systemie IAEA.

Reaktory na szybkich neutronach SMR: Projekty jak TerraPower Natrium (SFR — Sodium-cooled Fast Reactor) lub ARC-100 używają paliwa metalicznego lub tlenków z wysokim wzbogaceniem. Masa krytyczna dla neutronów prędkich jest wyższa (bo nie ma moderacji), ale czas pokolenia neutronów jest krótszy, co stawia inne wymagania na systemy bezpieczeństwa.

Reaktory kosmiczne: Reaktory dla misji kosmicznych (np. Kilopower NASA/DoE, ~1–10 kWe) używają uranu HEU lub HALEU w bardzo kompaktowych rdzeniach. Masa rdzenia to kilka kilogramów, a minimalna masa krytyczna jest osiągana przez geometrię i reflektor berylowy. Bezpieczeństwo w razie awarii wymagało zaprojektowania geometrii, w której wodna woda morska nie może prowadzić do krytyczności w razie lądowania awaryjnego na oceanie.⁵

Zestawienie modeli i kodów do wyznaczania masy krytycznej

Precyzyjne wyznaczenie masy krytycznej w skomplikowanych geometriach jest dziś domeną kodów komputerowych, a nie obliczeń analitycznych.

MCNP (Monte Carlo N-Particle): Kod opracowany w Los Alamos, dostępny dla uprawnionych użytkowników przez RSICC (Radiation Safety Information Computational Center). Monte Carlo symuluje historię każdego neutronu indywidualnie, z losowaniem kolejnych zdarzeń (rozpraszanie, pochłanianie, rozszczepienie, ucieczka) zgodnie z przekrojami czynnymi z baz danych ENDF/B. Dokładność jest limitowana tylko dokładnością przekrojów czynnych i liczbą historii neutronów; błędy statystyczne można zmniejszać przez wydłużenie obliczeń.

KENO (w zestawie SCALE): Kod Monte Carlo opracowany przez ORNL, szczególnie popularny dla obliczeń bezpieczeństwa krytycznościowego. Dostępny dla uprawnionych użytkowników przez RSICC.

OpenMC: Open-source Monte Carlo code opracowany na MIT, dostępny publicznie przez GitHub. Jest używany w badaniach akademickich i w niepoufnych obliczeniach projektów reaktorowych.

Kody deterministyczne: DORT, PARTISN, DENOVO (metoda $S_N$ discrete ordinates) — rozwiązują przestrzennie-kątowe równanie transportu neutronów bez metody Monte Carlo. Szybkie, ale wymagają wielogrupowego spektrum energetycznego i dokładnej geometrii.

Każdy z tych kodów wymaga biblioteki danych jądrowych (ENDF/B, JEFF, JENDL) zawierającej przekroje czynne dla każdego izotopu i każdego procesu (elastyczne, nieelastyczne rozpraszanie, wychwyt, rozszczepienie) w szerokim zakresie energii. Jakość danych jest kluczowa: błędy rzędu 1–2% w przekrojach czynnych przekładają się na porównywalne błędy w masie krytycznej. Najlepsze kody (OpenMC, MCNP6) wyznaczają masę krytyczną z błędem statistycznym poniżej 0,1% przy milionach historii neutronów, pod warunkiem dobrej biblioteki danych jądrowych i poprawnie zbudowanej geometrii wejściowej. Walidacja kodów przez benchmarki ICSBEP jest obowiązkowa przy ich certyfikacji dla zastosowań przemysłowych i licencjonowania reaktorów.⁵

Wpływ zanieczyszczeń i izotopów ubocznych

W praktycznej fizyce masy krytycznej istotną rolę odgrywają zanieczyszczenia i izotopy uboczne obecne w materiałach rozszczepialnych:

Pluton-240 w plutonie reaktorowym (WGPu vs RGPu): Pluton z reaktorów energetycznych (Reactor Grade Plutonium, RGPu) zawiera 20–25% Pu-240 (wobec <6% w Weapons Grade Plutonium, WGPu). Pu-240 nie jest bezpośrednio szkodliwy dla masy krytycznej, ale ma wysoką spontaniczną szybkość rozszczepienia (~1060 n/s/g), która powoduje ryzyko preinicjacji w bombach działowych. To właśnie jeden z powodów, dla których RGPu jest uważany za niepraktyczny dla budowy broni działowej.

Uran-234 w HEU: Naturalne paliwo zawiera ~0,0054% U-234. W wysoko wzbogaconym uranie (>90% U-235) zawartość U-234 wzrasta proporcjonalnie do wzbogacenia, osiągając ~0,5–1%. U-234 ma bardzo duży przekrój czynny na wychyt neutronów termicznych (~100 barn), co zwiększa straty w układach z moderatorem. Dla broni (neutony prędkie) to mniejszy problem.

Kontaminacja borem lub innym absorbentem: Nawet kilka ppm boru (B-10, $\sigma_a = 3840\,\mathrm{barn}$) lub gadolinu ($\sigma_a(Gd-157) \approx 60{,}000\,\mathrm{barn}$) w materiale rozszczepialnym może podnieść masę krytyczną powyżej wartości dla czystego izotopu. To nakłada wysokie wymagania na czystość izotopową w produkcji materiałów jądrowych.⁵

Analogy klasyczne i popularnonaukowe

Intuicja masy krytycznej daje się opisać kilkoma użytecznymi analogiami dydaktycznymi:

Analogia deszczu w lesie (opisana w streszczeniu artykułu): Rzadki las (mała gęstość materiału) pozwala wielu „kroplom-neutronom" ulec ucieczce. Gęsty las (duża gęstość lub duża geometria) zatrzymuje więcej kropel. Ściana wokół lasu (reflektor) zawraca część do środka. Analogia jest jakościowo trafna, choć ilościowo neutronowa fizyka jest bardziej złożona.

Analogia lawiny: Masa krytyczna jest jak nachylenie górskie, po którym zaczyna się samopodtrzymująca lawina. Poniżej progu każdy kamień zatrzymuje się sam. Powyżej progu każdy następny kamień uderza w kolejne, tworząc reakcję łańcuchową. Dodanie reflektora to jak zwiększenie nachylenia stosu przy tej samej ilości śniegu.

Analogia sieci komórkowej: Reaktywność $k_\mathrm{eff}$ to jak R₀ w epidemiologii (liczba reprodukcji). Gdy R₀ > 1, epidemia rośnie; gdy < 1 gaśnie. Masa krytyczna odpowiada wartości populacji, poniżej której gęstość kontaktów jest za mała, by utrzymać R₀ > 1. Analogia jest szczególnie pouczająca dla osób z biologicznym backgroundem.

Analogie te są przydatne w popularyzacji, ale mają swoje granice. Najważniejszy aspekt, który analogie pomijają, to wykładniczy, lawinowy charakter wzrostu w fazie nadkrytycznej — różnica między $k=1{,}01$ a $k=1{,}5$ jest dramatyczna w skali czasu mikroskopowej, ale w analogiach liniowych nie jest uchwytna.

Dla osób zainteresowanych głębszym rozumieniem: najdokładniejsze publiczne ujęcie masy krytycznej dla różnych materiałów i konfiguracji zawiera baza danych ICSBEP (International Criticality Safety Benchmark Evaluation Program), dostępna przez OECD/NEA. Zawiera ona ponad 5000 zestawów benchmarkowych (experimental critical assemblies) z lat 1942–2020, skrupulatnie zweryfikowanych i zrekonstruowanych przez ekspertów. Każdy zestaw opisuje geometrię, skład materiałowy i mierzony $k_\mathrm{eff}$ wraz z niepewnościami, i służy jako dane walidacyjne dla kodów obliczeniowych. To najbardziej kompleksowe jawne źródło danych o krytyczności w historii fizyki jądrowej, a zarazem główne narzędzie do certyfikacji kodów bezpieczeństwa krytycznościowego na całym świecie.^1,2,3,5

Dodatkowe materiały multimedialne

Powiązane kalkulatory i narzędzia

• Masa krytyczna — porównuje wpływ materiału, gęstości, reflektora i geometrii na masę krytyczną.
• k_eff — pokazuje, jak geometria, moderator i straty neutronów wpływają na krytyczność układu.
• Wizualizacja: Układ rozszczepieniowy — Koncepcyjny przekrój pokazuje złożenie typu gun, implozję i ogólny wzrost k_eff bez parametrów konstrukcyjnych.
• Model 3D: Geometrie krytyczne — 6 modeli — sfera/walec/płyta/sześcian.

Ćwiczenia praktyczne

Ćwiczenie laboratoryjne powinno polegać na zbudowaniu komputerowego modelu estymującego masę krytyczną dla uproszczonych geometrii. Najprostsza wersja może być oparta na symulacji Monte Carlo ruchu neutronów w sferze i cylindrze, z losowaniem kolejnych zdarzeń: rozproszenia, pochłonięcia, rozszczepienia i ucieczki przez granicę układu. Celem nie jest wierne odwzorowanie pełnego transportu neutronów, lecz porównanie, jak zmienia się estymowane $k_{eff}$ przy tej samej masie materiału, ale innym kształcie i promieniu. Wynik powinien pokazać, że kula osiąga krytyczność przy mniejszej masie niż układ o mniej korzystnym stosunku powierzchni do objętości.

Wariant rozszerzony powinien dodać dwa parametry: zmianę gęstości oraz zewnętrzny reflektor. W pierwszej części wystarczy przyjąć, że średnia droga swobodna skaluje się odwrotnie z gęstością i sprawdzić, jak estymowany próg krytyczności przesuwa się po zwiększeniu $\rho$. W drugiej części można otoczyć rdzeń warstwą materiału reflektorowego o zadanym prawdopodobieństwie rozpraszania i małym prawdopodobieństwie pochłaniania. Raport z laboratorium powinien zawierać:

1. zestaw przyjętych uproszczeń fizycznych,
2. porównanie sfery i cylindra,
3. wpływ wzrostu gęstości na estymowaną masę krytyczną,
4. wpływ reflektora na redukcję strat neutronów,
5. ocenę, które uproszczenia najbardziej zniekształcają wynik.

Drugie ćwiczenie, teoretyczno-przemysłowe, powinno polegać na przeliczeniu skali oszczędności materiału rozszczepialnego przy zmianie gęstości oraz przy zastosowaniu reflektora. Punktem wyjścia jest prawo skali:

$$m_{kr} \sim \frac{1}{\rho^2}$$

Należy więc:

1. przyjąć modelową wartość masy krytycznej dla sfery o gęstości normalnej,
2. policzyć, jak zmieni się ta masa po wzroście gęstości o 25%, 50% i 100%,
3. oszacować dodatkowy efekt reflektora jako redukcję strat neutronów w przyjętym modelu,
4. porównać, która metoda daje większy zysk: sama geometria, sama kompresja czy kompresja z reflektorem,
5. opisać, jakie ograniczenia technologiczne pojawiają się przy przejściu od modelu teoretycznego do rzeczywistego urządzenia.

To ćwiczenie ma pokazać rzecz praktyczną: w rzeczywistym projekcie przemysłowym nie „kupuje się po prostu większej ilości materiału”, lecz optymalizuje cały układ tak, by jak najmniejsza masa fizyczna zachowywała się jak układ nadkrytyczny w wymaganym przedziale czasu.

Ten temat najlepiej spina się z tekstami o reakcji łańcuchowej i współczynniku mnożenia k, przekroju czynnym na rozszczepienie oraz o reflektorze i tamperze uranowym. Razem pokazują one, że masa krytyczna nie jest tylko problemem ilości materiału, lecz także geometrii, zatrzymywania neutronów i szybkości osiągania nadkrytyczności.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Reakcja łańcuchowa i współczynnik mnożenia k
• Metoda implozyjna - mechanizm Fat Man
• Demon Core - tragiczne wypadki w Los Alamos
• One-point safety: dlaczego przypadkowe odpalenie jednego detonatora nie może dawać istotnego uzysku jądrowego