Reakcja łańcuchowa i współczynnik mnożenia k

Streszczenie

Jeżeli w szeregu kostek domina jedna przewraca dokładnie jedną następną, cały układ pozostaje w równowadze: proces trwa, ale nie przyspiesza. Jeżeli jedna kostka przewraca średnio mniej niż jedną kolejną, ciąg zaraz się urywa. Jeżeli zaś każda uruchamia więcej niż jedną następną, proces narasta lawinowo. Reakcja łańcuchowa rozszczepienia działa według tej samej logiki, tylko zamiast kostek domina mamy jądra atomowe, a zamiast impulsu mechanicznego neutrony. To właśnie dlatego fizycy tak pilnują jednej liczby: współczynnika mnożenia $k$. Od niego zależy, czy układ z materiałem rozszczepialnym jest podkrytyczny, krytyczny czy nadkrytyczny.^1,2

W praktyce ten sam język opisuje zarówno reaktor, jak i bombę, ale znaczenie fizyczne nie jest identyczne. W reaktorze dąży się do $k = 1$, aby moc była stała i możliwa do sterowania. W ładunku jądrowym trzeba na bardzo krótki czas uzyskać $k > 1$, i to wyraźnie większe od jedności, bo tylko wtedy liczba rozszczepień rośnie dostatecznie szybko, zanim rozgrzany materiał zacznie się rozszerzać i sam obniży reaktywność układu. Zrozumienie tej różnicy jest ważniejsze niż zapamiętanie samej definicji.^1,3

Rozszerzenie tematu

Najprostsza definicja reakcji łańcuchowej rozszczepienia jest taka: neutron wywołuje rozszczepienie jądra, a to rozszczepienie emituje kolejne neutrony, z których część powoduje następne rozszczepienia. Wymuszone rozszczepienie uranu-235 daje średnio około 2,52 neutronu, a plutonu-239 około 2,95 neutronu. Sama liczba wyemitowanych neutronów nie wystarcza jednak do powstania samopodtrzymującego procesu. Część neutronów ucieka z układu, część zostaje pochłonięta bez rozszczepienia, a tylko pewna frakcja trafia do następnego jądra w taki sposób, że rzeczywiście uruchamia kolejny akt rozszczepienia.¹

Właśnie tę skuteczną część opisuje współczynnik mnożenia $k$. W podręcznikowym zapisie Słowińskiego jest to stosunek liczby neutronów w kolejnym pokoleniu do liczby neutronów w pokoleniu poprzednim:

$$k = \frac{n_i}{n_{i-1}}$$

Jeżeli $k < 1$, układ jest podkrytyczny i reakcja wygasa. Jeżeli $k = 1$, układ jest krytyczny i populacja neutronów utrzymuje się średnio na stałym poziomie. Jeżeli $k > 1$, układ jest nadkrytyczny i liczba neutronów, a więc także tempo wydzielania energii, rośnie z pokolenia na pokolenie.² To brzmi prosto, ale za tą prostotą stoi cała fizyka bilansu neutronów.

W układzie idealnym, nieskończenie dużym, można mówić o $k_{\infty}$, czyli współczynniku mnożenia w objętości na tyle dużej, że ucieczkę neutronów można pominąć. Taki parametr jest wygodny, bo pozwala najpierw opisać samą jakość środowiska neutronowego: materiał paliwa, moderator, prawdopodobieństwo pochłaniania i rozszczepienia. W realnym układzie neutrony jednak uciekają przez powierzchnię, dlatego dla rzeczywistego reaktora używa się $k_{ef}$ lub $k_{eff}$, czyli efektywnego współczynnika mnożenia. Słowiński zapisuje tę zależność jako:

$$k_{ef} = k_{\infty} \cdot p_u$$

gdzie $p_u$ jest prawdopodobieństwem uniknięcia ucieczki neutronów ze strefy aktywnej.² Z tego od razu widać, dlaczego geometria układu jest tak ważna: nawet dobry materiał rozszczepialny nie utrzyma reakcji, jeśli bryła będzie zbyt mała lub źle ukształtowana. To zagadnienie rozwija osobny artykuł o masie krytycznej i parametrach geometrycznych.

W reaktorze termicznym pełny cykl neutronowy wygodnie opisuje tak zwana formuła czteroczynnikowa:

$$k_{\infty} = \epsilon \cdot p \cdot f \cdot \eta_T$$

Każdy z tych czynników odpowiada innemu etapowi losu neutronu. $\epsilon$ uwzględnia dodatkowe rozszczepienia na neutronach prędkich, $p$ oznacza prawdopodobieństwo uniknięcia wychwytu rezonansowego podczas spowalniania, $f$ mówi, jaka część neutronów cieplnych zostanie pochłonięta tam, gdzie chcemy, czyli w paliwie, a nie w materiale konstrukcyjnym lub moderatorze, natomiast $\eta_T$ określa, ile neutronów rozszczepieniowych powstaje średnio na jeden neutron pochłonięty w paliwie.² Ta formuła nie jest tylko rachunkową ciekawostką. To właśnie ona pozwala zrozumieć, dlaczego o powodzeniu programu jądrowego decydowały czasem detale materiałowe.

Dobrym przykładem jest kwestia moderatora grafitowego. Jeżeli moderator albo materiały konstrukcyjne pochłaniają zbyt wiele neutronów, bilans neutronowy od razu się pogarsza i układ może nie dojść do $k_{ef} = 1$, mimo że sam materiał paliwowy jest poprawnie dobrany. Wczesne eksperymenty Fermiego pokazują to dobrze: kolejne stosy eksperymentalne wymagały materiałów o dostatecznej czystości, tak aby projektowane $k$ przestało być zbyt małe i zaczęło zbliżać się do krytyczności.^2,4 Właśnie tu widać, że reakcja łańcuchowa jest zjawiskiem fizycznym, ale jej praktyczna realizacja jest zawsze problemem materiałowym i inżynierskim.

Warto też odróżnić opis reaktorowy od opisu eksplozji jądrowej. W reaktorze znaczenie mają zarówno neutrony natychmiastowe, jak i opóźnione, ponieważ to właśnie neutrony opóźnione dają operatorowi czas na sterowanie i utrzymywanie mocy w pobliżu stanu krytycznego. W ładunku jądrowym skala czasu jest jednak tak mała, że w praktyce liczą się tylko neutrony natychmiastowe. Gdenarz podkreśla, że przy opisie eksplozji przez $k_{eff}$ uwzględnia się właśnie neutrony natychmiastowe, a nie opóźnione, bo cały proces musi rozwinąć się, zanim materiał zdąży się rozprężyć.³

Samo $k > 1$ jeszcze nie mówi, jak szybko rośnie liczba rozszczepień. Do tego potrzebny jest czas życia jednej generacji neutronów. Gdenarz zapisuje tempo wzrostu populacji przez parametr:

$$\alpha = \frac{k_{eff} - 1}{\tau}$$

gdzie $\tau$ jest średnim czasem trwania pokolenia neutronów.³ Jeżeli $\alpha > 0$, populacja neutronów rośnie wykładniczo. W reaktorze ten wzrost ma być praktycznie zatrzymany i sprowadzony do stabilnego punktu pracy. W bombie przeciwnie: trzeba uzyskać bardzo duże dodatnie $\alpha$, zanim wzrost temperatury i ciśnienia rozerwie układ. Dlatego z punktu widzenia konstrukcji ładunku jądrowego samo osiągnięcie krytyczności nie wystarcza. Potrzebna jest jeszcze odpowiednio szybka droga do silnej nadkrytyczności.

Jest jeszcze jeden istotny efekt: podkrytyczne namnażanie neutronów. Jeżeli masa jest bliska krytyczności, ale jeszcze jej nie osiągnęła, pojedynczy neutron może mimo wszystko wywołać całą krótką, gasnącą kaskadę rozszczepień. Dla $k = 0{,}9$ taki łańcuch trwa średnio około 10 generacji, a dla $k = 0{,}99$ już około 100 generacji.¹ To ma duże znaczenie historyczne. Pokazuje bowiem, dlaczego proste mechaniczne składanie materiału rozszczepialnego było tak trudne: układ nie może zbyt długo przebywać w pobliżu krytyczności, bo przypadkowy neutron uruchomi reakcję za wcześnie. Ten problem prowadzi wprost do zagadnień omawianych w artykułach o predetonacji, metodzie działowej i metodzie implozyjnej.

Zależność między $k$, geometrią i gęstością tłumaczy też, dlaczego w konstrukcjach militarnych tak ważne są reflektory i tampers. Reflektor może zawrócić część neutronów, które inaczej opuściłyby rdzeń, a tamper opóźnia hydrodynamiczne rozbieganie się materiału po rozpoczęciu reakcji. W obu przypadkach chodzi o to samo: przez bardzo krótki czas poprawić gospodarkę neutronową układu i utrzymać $k$ powyżej jedności wystarczająco długo, by rozszczepienie zdążyło uwolnić znaczną energię.^1,3

Z praktycznego punktu widzenia najważniejsza lekcja jest więc następująca. Reakcja łańcuchowa nie jest po prostu serią rozszczepień. Jest bilansem neutronów. Każdy proces, który poprawia ten bilans, podnosi $k$: większa masa, lepszy kształt, większa gęstość, mniejsza liczba pochłaniaczy, skuteczniejszy reflektor. Każdy proces, który go pogarsza, obniża $k$: mała bryła, niekorzystna geometria, zanieczyszczenia, ucieczka neutronów, słabe dobranie moderatora. Dlatego współczynnik mnożenia jest jednocześnie pojęciem fizycznym, inżynierskim i historycznym. Na nim spotykają się teoria reaktora, konstrukcja pierwszych bomb i pytanie, dlaczego jedne programy jądrowe osiągnęły krytyczność, a inne zatrzymały się przed nią.

Podsumowując: $k$ nie jest tylko symbolem z podręcznika. To najkrótszy możliwy opis odpowiedzi na pytanie, czy neutrony w danym układzie zarabiają na siebie. Gdy zarabiają dokładnie na jednego następcę, mamy stabilną krytyczność. Gdy nie zarabiają nawet na jednego, reakcja gaśnie. Gdy zarabiają na więcej niż jednego i robią to wystarczająco szybko, układ wchodzi w obszar gwałtownego wzrostu mocy.

Historia: Szilard, Fermi i droga do CP-1

Pojęcie reakcji łańcuchowej ma zaskakująco konkretną datę narodzin jako idea naukowa. Leo Szilard — węgierski fizyk, który uciekł z nazistowskich Niemiec do Anglii — doznał olśnienia podczas spaceru w Londynie w 1933 roku, rok po odkryciu neutronu przez Chadwicka. Wyobraził sobie, że istnieje pierwiastek, który po trafieniu neutronem emituje dwa neutrony, co pozwoliłoby na łańcuchowy wzrost reakcji. Opatentował tę ideę w 1934 roku i przekazał patent Admiralicji Brytyjskiej, starając się utajnić koncepcję.⁴

Szilard myślał początkowo o berze, indu i innych pierwiastkach — nie o uranie, bo rozszczepienia jądrowego jeszcze nie odkryto. Gdy w 1939 roku Meitner i Frisch wyjaśnili fizykę rozszczepienia U-235, Szilard natychmiast połączył fakty: oto pierwiastek, który po trafieniu neutronem emituje ~2,5 neutronu. Razem z Fermim złożył petycję do prezydenta Roosevelta (słynny list Einsteina-Szilarda z 2 sierpnia 1939) ostrzegającą przed możliwością zbudowania bomby przez Niemcy.

Program Manhattan wyłonił się z tej inicjatywy. Kluczowym krokiem była weryfikacja eksperymentalna: czy rozszczepienie U-235 z neutronami termicznymi może dać $k > 1$ w realnym układzie? Fermi i jego zespół budowali kolejne stosy eksperymentalne (Chicago Pile 0, Chicago Pile 1), mierząc współczynnik mnożenia. Każda iteracja była ulepszeniem: czystszy grafit (mniej pochłaniania neutronów), czystszy UO₂ i metal uranowy, lepsza geometria.

Chicago Pile-1 osiągnął krytyczność 2 grudnia 1942 roku o godzinie 15:25 — George Weil ostrożnie wyciągał pręt kontrolny z grafitu, a Fermi śledził licznik neutronów. Przy pewnym położeniu pręta współczynnik mnożenia osiągnął $k_{eff} = 1{,}0006$ i moc zaczęła rosnąć. Po 28 minutach Fermi wydał rozkaz wbicia pręta z powrotem. Herbert Anderson nadał do Waszyngtonu zaszyfrowaną wiadomość: „Włoski żeglarz dopłynął do Nowego Świata i rdzenni mieszkańcy są przyjacielscy."⁴

CP-1 miał moc zaledwie 200 W — celowo tak mała, bo nie miał układu chłodzenia. Ale zademonstrował kluczowy fakt: reakcja łańcuchowa oparta na U-235 i neutronach termicznych z grafitem jako moderatorem jest możliwa. To przełom dający zielone światło dla reaktorów produkcyjnych w Hanford i dla programu broni.

Formuła sześcioczynnikowa

Uproszczona formuła czteroczynnikowa opisuje tylko reaktory termiczne i pomija kilka efektów. Pełna formuła sześcioczynnikowa ujmuje całość:

$$k_{eff} = \epsilon \cdot p \cdot f \cdot \eta_T \cdot P_{NL}^T \cdot P_{NL}^F$$

Nowe czynniki to:

• $P_{NL}^T$ — prawdopodobieństwo nieuciekania termicznych neutronów (thermal non-leakage probability)
• $P_{NL}^F$ — prawdopodobieństwo nieuciekania prędkich neutronów (fast non-leakage probability)

Iloczyn $P_{NL}^T \cdot P_{NL}^F$ odpowiada $p_u$ ze wzoru Słowińskiego. Jednak rozbicie na dwa składniki — dla neutronów termicznych i prędkich — jest ważne, bo mają one różne długości dyfuzji i różne migracje:

$$P_{NL}^F = \frac{1}{1 + B^2 \tau} \qquad P_{NL}^T = \frac{1}{1 + L^2 B^2}$$

gdzie $B^2$ to geometryczny buckling (parametr kształtu), $\tau$ to wiek Fermiego moderatora, a $L$ to długość dyfuzji neutronów termicznych. Te wzory wywodzą się z teorii dyfuzji neutronów, opisanej szerzej w artykule o termicznej dyfuzji neutronów.

Dla reaktora sferycznego o promieniu $R$:

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

Zależność $k_{eff}$ od rozmiaru przez $B^2$ wyjaśnia, dlaczego mała bryła materiału rozszczepialnego jest podkrytyczna: małe $R$ daje duże $B^2$, co obniża oba czynniki nieuciekania i w konsekwencji $k_{eff} < 1$. To właśnie matematyczna postać warunku masy krytycznej.²

Kinetyka punktowa i równania kinetyki neutronowej

W reaktorze potrzebna jest nie tylko znajomość stanu stacjonarnego ($k_{eff} = 1$), ale i dynamiki zmian mocy w czasie. Opis taki dają równania kinetyki punktowej (Point Kinetics Equations, PKE):

$$\frac{dn(t)}{dt} = \frac{\rho - \beta_{eff}}{\Lambda} n(t) + \sum_{i=1}^{6} \lambda_i C_i(t)$$

$$\frac{dC_i(t)}{dt} = \frac{\beta_i}{\Lambda} n(t) - \lambda_i C_i(t), \quad i = 1, \ldots, 6$$

gdzie:

• $n(t)$ — populacja neutronów (lub moc reaktora)
• $\rho$ — reaktywność $= (k_{eff}-1)/k_{eff}$
• $\beta_{eff}$ — efektywna frakcja neutronów opóźnionych
• $\Lambda$ — sredni czas generacji neutronów (ang. prompt neutron lifetime)
• $C_i(t)$ — stężenie i-tego prekursora neutronów opóźnionych
• $\lambda_i$, $\beta_i$ — stała rozpadu i frakcja i-tej grupy prekursorów

Dla U-235: $\beta_{eff} \approx 0{,}0065$, $\Lambda \approx 50 \mu\mathrm{s}$ (reaktory LWR z moderatorem). Tak długo jak $\rho < \beta_{eff}$ (poniżej prompt criticality), czas podwajania mocy jest rzędu sekund — bo neutrony opóźnione wydłużają efektywny czas generacji do $\sim \beta_{eff} \cdot \Lambda / |\beta_{eff} - \rho| \sim\$sekund. To daje czas reakcji operatorowi lub układom automatyki.⁶

Gdy $\rho > \beta_{eff}$ — stan prompt supercritical — moc rośnie z czasem charakterystycznym $\Lambda \approx 50 \mu\mathrm{s}$, czyli milion razy szybciej. To jest stan, który nastąpił w Czarnobylu 26 kwietnia 1986 roku: po testie turbiny moc reaktora RBMK wzrosła do ~30 000 MW (nominalnie 3200 MW) w czasie ~3 sekund.

Czynniki reaktywności i współczynniki temperatury

W rzeczywistym reaktorze $k_{eff}$ zmienia się wraz z temperaturą, ciśnieniem i składem paliwa. Te zmiany opisują współczynniki reaktywności:

Współczynnik temperatury moderatora ($\alpha_M = d\rho/dT_M$): jeśli temperatura wody rośnie, jej gęstość spada → mniej moderacji → neutrony termiczne mają wyższe energie → mniejszy $\sigma_f$ U-235 → mniejszy $\eta_T$ → mniejszy $k_{eff}$. Dla LWR: $\alpha_M \approx -30 \mathrm{pcm/K}$ (pcm = per cent mille, $10^{-5}$). Ujemna wartość to kluczowy mechanizm bezpieczeństwa.

Współczynnik temperatury paliwa (efekt Dopplera) ($\alpha_D = d\rho/dT_P$): wzrost temperatury paliwa poszerza rezonanse U-238 przez efekt Dopplera → więcej pochłaniania w rezonansach → mniejszy $p$ → mniejszy $k_{eff}$. Dla LWR: $\alpha_D \approx -3 \mathrm{pcm/K}$. Działa natychmiastowo (bo paliwo reaguje szybciej niż moderator) — to pierwsza linia samoochrony reaktora.⁶

Współczynnik pustek ($\alpha_V$): dla LWR (gdzie moderatorem jest woda) tworzenie się pęcherzyków pary obniża moderację → mniejszy $k_{eff}$. Dla RBMK (grafit jako moderator, woda jako chłodziwo) tworzenie się pary w kanałach chłodziwa zmniejsza pochłanianie neutronów przez wodę → więcej neutronów dociera do grafitu → większy $k_{eff}$. To jest tzw. dodatni współczynnik pustek RBMK — główna przyczyna niestabilności Czarnobyla.

Efekt zatruwania ksenonem: po zmniejszeniu mocy Xe-135 narasta przez 6–8 h (przejściowo $k_{eff}$ spada). Jeśli operator próbuje utrzymać moc przez wyciąganie prętów kontrolnych, układ staje się trudny do kontroli. Ten właśnie scenariusz rozegrał się w Czarnobylu w nocy z 25 na 26 kwietnia 1986.⁶

Reaktory prędkie i krytyczność bez moderatora

Reaktory termiczne wymagają moderatora (H₂O, D₂O, grafit), bo neutrony prędkie ($E \sim 1 \mathrm{MeV}$) mają zbyt niski $\eta$ i zbyt małe $\sigma_f$ dla U-235, by utrzymać $k > 1$ bez spowolnienia. Jednak reaktory prędkie (SFR — Sodium-cooled Fast Reactor, GFR — Gas-cooled Fast Reactor) pracują właśnie z neutronami prędkimi — cielęciem chłodziwa jest ciekły sód (lub hel, ołów), który nie moderuje neutronów.

W środowisku prędkim bilans neutronowy jest inny:

• Przekroje czynne są mniejsze (brak 1/v zależności dla neutronów termicznych)
• Ale $\eta$ (neutrony na pochłoniecie) jest wyższy dla Pu-239: $\eta_{fast} \approx 3{,}0$ vs $\eta_{thermal} \approx 2{,}9$
• U-238 ma znaczący $\sigma_f$ przy $E > 1 \mathrm{MeV}$, więc może uczestniczyć w łańcuchu
• Współczynnik $\epsilon$ (rozszczepienia prędkie) jest bardzo duży — właściwie prawie każde rozszczepienie jest prędkie

Reaktory prędkie mogą być „breeders" (rozkładacze paliwa) — produkujące więcej Pu-239 niż spalają U-235. Warunkiem jest $\eta_{fast}(Pu-239) - 1 > 1$, czyli jeden neutron na podtrzymanie reakcji i jeden na zamianę U-238 w Pu-239. Historyczne reaktory prędkie: EBR-I (Idaho, 1951 — pierwszy reaktor generujący prąd), EBR-II, BN-800 (Rosja, 2016), Superphénix (Francja, 1985–1998). Planowane: BN-1200 (Rosja), ASTRID (Francja) — anulowany, PFBR (Indie, opóźniony).⁶

Ładunek jądrowy: superkrytyczność i czas eksplozji

W ładunku jądrowym cel jest diametralnie inny niż w reaktorze. Musi nastąpić gwałtowny wzrost mocy od zera do petawatów w mikrosekundy. Wymaga to:

1. Szybkiego przejścia do głębokiej superkrytyczności: $k_{eff}$ musi wzrosnąć od ~0,9 (stan montowania) do ~2 lub więcej (stan eksplozji) w czasie ~1 μs. Osiąga się to przez implozję (szybkie ściskanie materiału, wzrost gęstości → wzrost $k_{eff}$) lub metodę działową (szybkie składanie masy w konfigurację nadkrytyczną).

2. Inicjacji w odpowiednim momencie: neutron inicjujący musi pojawić się w momencie szczytowej superkrytyczności. Za wcześnie → predetonacja (fizzle). Za późno → brak inicjacji.

3. Podtrzymania łańcucha przez co najmniej kilka czasów generacji: typowy czas generacji neutronów prędkich w HEU $\approx 10 \mathrm{ns}$. Efektywna eksplozja wymaga kilkudziesięciu generacji przed rozerwaniem ładunku przez rozprężanie plazmowe.

Parametr $\alpha$ z równania $dn/dt = \alpha n$ wyznacza czas podwajania populacji:

$$t_{2\times} = \frac{\ln 2}{\alpha} = \frac{\tau \ln 2}{k_{eff} - 1}$$

Dla $k_{eff} = 2$ i $\tau = 10 \mathrm{ns}$: $t_{2\times} = 6{,}9 \mathrm{ns}$. Po 60 ns (zaledwie 9 podwojeń) moc wzrasta 512-krotnie. Ten wykładniczy wzrost trwa do momentu, gdy grawitacja materiału i ciśnienie promieniowania X zaczną rozsadzać układ — dopiero wtedy $k_{eff}$ zaczyna spadać.³

Wypadki krytyczne: historia i lekcje bezpieczeństwa

Historia fizyki jądrowej zawiera kilkanaście nieumyślnych zdarzeń krytycznych — sytuacji, gdy $k_{eff}$ niespodziewanie przekroczył jedność w sposób niekontrolowany. Każdy z tych wypadków wzbogacił rozumienie reakcji łańcuchowej o nowe lekcje:

Demon's Core (1945–1946, Los Alamos): Sfera Pu-239 (~6,2 kg) była używana do krytycznych eksperymentów. Harry K. Daghlian Jr. upuścił klocek reflektora na sferę (21 sierpnia 1945) — błysk promieniowania Czerenkowa i impuls neutronów. Daghlian zmarł 25 dni później. Rok później (21 maja 1946) Louis Slotin podczas pokazywania eksperymentu omyłkowo opuścił reflektor berylowy na kulę plutonu — ponowny impuls. Slotin zasłonił sferę własnym ciałem chroniąc obserwatorów i zmarł 9 dni później.

Obydwa wypadki wynikały z tej samej przyczyny: ręczne kontrolowanie konfiguracji bliskiej krytycznej, bez automatycznych zabezpieczeń. Po 1946 roku wszystkie eksperymenty krytyczne w USA były prowadzone przez zdalne sterowanie.⁴

SL-1 (3 stycznia 1961, Idaho): Eksperymentalny reaktor SL-1 (Stationary Low-Power Reactor-1) uległ katastrofalnemu wybuchowi parowo-jądrowemu w wyniku ręcznego wyjęcia pręta kontrolnego na pełną długość (normalne wyjęcie wynosiło 10 cm, operater wyciągnął ~50 cm). Moc wzrosła od ~3 kW do ~20 000 MW w ułamku sekundy. Trzy osoby zginęły — jeden operator został przygnieciony do sufitu wyrzutym prętem.

Analiza wykazała, że $k_{eff}$ wzrósł powyżej 1+$\beta$ w ciągu ~3 ms — wejście w prompt supercriticality. Energia wyzwolona (~1,3 GJ) zamieniła wodę moderatora w parę, co rozsadziło reaktor.⁶

Tokaimura (30 września 1999, Japonia): Wypadek podczas przetwarzania paliwa uranowego w zakładzie JCO w Tokaimura. Trzech pracowników wlało do zbiornika 16 kg wzbogaconego uranu (wzbogacenie ~19%) zamiast standardowych porcji — po przekroczeniu masy krytycznej (roztworu, nie metalicznego uranu) układ przeszedł do oscylacyjnej reakcji łańcuchowej. Trwała ona 20 godzin. Dwóch pracowników przyjęło dawki śmiertelne (10 Gy i 6 Gy). Poza terenem zakładu 660 osób zostało narażonych na promieniowanie. Wypadek był wynikiem złamania procedur bezpieczeństwa i braku szkolenia personelu.

Lekcja wspólna: Każdy z tych wypadków był możliwy dlatego, że operatorzy nie mieli wystarczającego wyczucia tego, jak blisko krytyczności jest układ, albo celowo łamali procedury bezpieczeństwa. Bezpieczeństwo krytycznościowe (criticality safety) jest dziś sformalizowane w normie ANSI/ANS-8.1 i analogicznych normach ISO — definiując dozwolone geometrie, masy, stężenia i procedury dla każdego etapu obchodzenia się z materiałami rozszczepialnymi.⁴

Kody numeryczne i obliczenia krytyczności

Współczesne obliczenia krytyczności są prowadzone kodami Monte Carlo transportu neutronów, które śledzą każdy neutron od narodzin (rozszczepienie) do śmierci (pochłonięcie lub ucieczka z układu). Główne kody:

MCNP (Monte Carlo N-Particle, LANL): najpowszechniejszy na świecie, rozwijany od 1963 roku. Bieżąca wersja MCNP6 obsługuje neutronowe, fotonowe i elektronowe transporty. Oblicza $k_{eff}$ metodą fission source iteration (pętla: generuj neutrony z rozszczepień, propaguj, zapisz nowe lokalizacje rozszczepień, powtarzaj). Wymagana dokładność $\Delta k_{eff} < 0{,}001$ dla aplikacji bezpieczeństwa.

KENO (SCALE/ORNL): kod Monte Carlo specjalnie zoptymalizowany dla obliczeń krytyczności. Część pakietu SCALE (Standardized Computer Analyses for Licensing Evaluation), używanego przez UE i NRC do certyfikacji projektów reaktorów.

OpenMC: otwartoźródłowy kod Monte Carlo (MIT), rozwijany od 2011 roku. Obsługuje ciągłe przekroje czynne (ACE format) i jest coraz szerzej używany w badaniach akademickich.

SERPENT2: fińsko-europejski kod Monte Carlo (VTT, Finlandia), specjalnie zoptymalizowany dla obliczeń wypalania paliwa i generowania stałych grupowych dla kodów deterministycznych.

Do walidacji kodów używa się bazy ICSBEP (International Criticality Safety Benchmark Evaluation Project), zawierającej >5000 zestawów benchmarkowych. Standard wymagany to zgodność z eksperymentem w granicach $\pm 0{,}003$ w $k_{eff}$.⁶

Monitoring i systemy ostrzegania o krytyczności

Zakłady jądrowe, w których obchodzi się materiałami rozszczepialnymi (wzbogacarnie, zakłady przetwarzania paliwa, laboratoria badawcze), wyposażone są w systemy monitoringu krytyczności. Ich zadaniem jest wykrycie gwałtownego wzrostu promieniowania (wskazującego na nadkrytyczność) i uruchomienie alarmów ewakuacyjnych zanim dawki będą śmiertelne.

Typowe systemy: CADS (Criticality Accident Detection System) — detektory He-3 lub scyntylacyjne rozmieszczone w zakładzie; przy wzroście tempa zliczeń o >10% nadtła w czasie <1 s wyzwalają alarm dźwiękowy i sygnał ewakuacji. Standard: DOE-STD-1039-2016 (USA), PN-EN 15237 (UE).

W Polsce program jądrowy jest nadzorowany przez Państwową Agencję Atomistyki (PAA), która nadzoruje stosowanie norm bezpieczeństwa krytycznościowego w ośrodkach badawczych (NCBJ Świerk, reaktor MARIA) i w przyszłości ma nadzorować planowane elektrownie jądrowe. PAA stosuje normę PN-EN ISO 817 i PN-EN 15237 jako odpowiedniki ANS-8.1.⁶

Reakcja łańcuchowa w broni termojądrowej

W ładunkach termojądrowych (układ Teller-Ulam) reakcja łańcuchowa rozszczepienia odgrywa dwie role:

1. Inicjacja: pierwotny ładunek rozszczepieniowy (tzw. primary lub fissile pit) jest inicjatorem całości. Musi osiągnąć $k_{eff} \gg 1$ i wytworzyć wystarczającą temperaturę (~100 milionów kelwinów) i strumień fotonów X do ściśnięcia i zapłonu paliwa fuzyjnego (secondary).

2. Wzmocnienie rozszczepienne (boosting): neutrony z fuzji D+T (14,1 MeV) wracają do paliwa rozszczepieniowego i wywołują dalsze rozszczepienia. Daje to efektywne wzmocnienie materiału rozszczepialnego, bo te wysokoenergetyczne neutrony mają duże $\sigma_f$ dla U-238 — co pozwala na rozszczepienie „koszulki" U-238 (jacketu) w bombach Mk/B61 i podobnych.

Całkowita wydajność termojądrowego ładunku (~50% z rozszczepienia, ~50% z fuzji i rozszczepienia koszulki) wynika z tego bilansu neutronowego, który jest precyzyjną wersją tego samego $k_{eff}$, tylko że teraz neutronami prędkimi z fuzji.^3,4

Zastosowania poza bronią i reaktorami

Reakcja łańcuchowa i współczynnik mnożenia mają zastosowania daleko poza reaktorami i bronią:

Pomiary pasywne i aktywne MAEA: weryfikacja deklaracji ilości materiałów rozszczepialnych przez MAEA często polega na pomiarach neutronów spontanicznych z rozszczepienia (pasywne systemy neutronowe, PNMM) lub na aktywnej weryfikacji $k_{eff}$ przez pomiary reaktywności (Differential Die-Away Analysis, DDA). Systemem takim można stwierdzić, czy masa deklarowanego materiału jest zgodna z mierzoną — bo $k_{eff}$ jest czuły na masę i geometrię.

Transmutacja odpadów: jądrowe odpady długożyciowe (Am-241, Am-243, Np-237, Cm-244) mogą być „spalane" w reaktorach prędkich przez rozszczepienie. Zamiana długożyciowych aktynowców w krótkotrwałe produkty rozszczepienia redukuje czas koniecznego składowania odpadów od 300 000 lat do ~300 lat. Program europejski MYRRHA (Multi-purpose hYbrid Research Reactor for High-tech Applications, Belgia) testuje podkrytyczny reaktor napędzany akceleratorem (ADS — Accelerator Driven System), gdzie źródło neutronów zewnętrznych (spallacyjnych) utrzymuje $k = 0{,}97$, a odpad jest transmutowany nawet bez osiągania krytyczności.

Systemy ADS i energia jądrowa IV generacji: zasada jest prosta — reaktor podkrytyczny ($k < 1$) jest bezpieczniejszy, bo wygaśnie samoistnie po wyłączeniu akceleratora. Produkcja neutronów spallacyjnych przez uderzanie protonami w grube tarcze Pb lub Bi-Pb jest efektywna: ~20–30 neutronów na proton przy energiach >600 MeV. To otwiera możliwość budowy systemów, które nie mogą „uciec" w stan nadkrytyczny z definicji geometrycznej.⁶

Chicago Pile-1 i pierwsze reaktory produkcyjne

Po sukcesie CP-1 (grudzień 1942) Fermi i jego zespół szybko zbudowali CP-2 i CP-3 — reaktory badawcze na potrzeby Projektu Manhattan. Ale kluczowe były reaktory produkcyjne w Hanford, WA, których celem była produkcja plutonu-239 dla bomby.

Reaktory B, D i F w Hanford (budowane w 1943–1944) były pierwszymi reaktorami przemysłowymi na świecie. Każdy o mocy ~250 MW termicznych, moderowane grafitem, chłodzone wodą z rzeki Kolumbia. Paliwo: niezbogacony uran metaliczny. Cel: napromieniowanie U-238 → Np-239 → Pu-239 przez wychwyt neutronu i dwa β⁻ rozpady. Zbiorczy wynik: produkcja ~10 kg Pu-239 miesięcznie z jednego reaktora — wystarczające do wyprodukowania jednej bomby co 2–3 miesiące.

Projekt reaktorów napotykał na nieoczekiwany problem: po uruchomieniu reaktora B moc zaczęła nagle spadać po kilku godzinach pracy. Fermi zidentyfikował przyczynę jako zatrucie ksenonem-135 — produkt rozszczepienia, który pochłania neutrony z ogromnym przekrojem czynnym ($\sigma = 2{,}65 \times 10^6 \mathrm{b}$). Efekt był przewidziany przez Johna Wheelera, który zalecił rozbudowanie rdzeni i dodanie zapasowych prętów paliwowych — „na wszelki wypadek". Ten margines okazał się kluczowy: reactor był wystarczająco duży, by skompensować efekt ksenonu i kontynuować pracę.⁴

Ta historia jest ważna pedagogicznie: nawet przy perfekcyjnym rozumieniu fizyki reakcji łańcuchowej (Fermi wiedział o ksenonie-135), niespodzianki mogą się zdarzyć. Reaktor Hanford-B był pierwszą demonstracją, że kontrola $k_{eff}$ w warunkach przemysłowych jest złożonym problemem inżynierskim, wykraczającym poza fizyczne obliczenia.

{{image-full:https://upload.wikimedia.org/wikipedia/commons/thumb/f/f3/Chicago_Pile-1_-_GPN-2003-00013.jpg/640px-Chicago_Pile-1_-_GPN-2003-00013.jpg:Chicago Pile-1 (CP-1) — pierwszy reaktor jądrowy na świecie, zbudowany przez Enrico Fermiego i jego zespół pod boiskiem squash na Uniwersytecie Chicago. Osiągnął krytyczność 2 grudnia 1942 roku. Sterta waży 360 ton grafitu z 5 ton metalicznego uranu i 46 ton dwutlenku uranu.}}

Reaktywność i pcm — jednostka praktyczna

W reaktorach jądrowych reaktywność $\rho$ jest wyraźną miarą odchylenia od stanu krytycznego:

$$\rho = \frac{k_{eff} - 1}{k_{eff}} \approx k_{eff} - 1 \quad \text{dla } k_{eff} \approx 1$$

Jednostką praktyczną jest pcm (per cent mille, $10^{-5}$), β (jako wielokrotność $\beta_{eff}$), lub nile/milli-nile (historyczne). Dla U-235: $\beta_{eff} = 650 \mathrm{pcm}$.

Typowe wartości w PWR:

• Świeże paliwo, zimny reaktor: $\rho \approx +2000 \mathrm{pcm}$ (kompensowane przez bor w wodzie)
• Stan pracy (gorący, z ksenonem, wypalony do 50%): $\rho \approx 0 \mathrm{pcm}$ (stan krytyczny)
• Wyłączony przez pręty kontrolne: $\rho \approx -3000 \ldots -5000 \mathrm{pcm}$ (margines zamknięcia)

Margines zamknięcia (shutdown margin) musi być wystarczający, by reaktor pozostał podkrytyczny nawet przy jednym zaciętym pręcie kontrolnym (single rod stuck) — wymóg licencyjny w każdym kraju posiadającym elektrownie jądrowe.⁶

Historia pojęcia: od Szilarda do Fermi Prize

Pojęcie reakcji łańcuchowej przeszło długą drogę od patentu Szilarda do precyzyjnego sformalizowania przez Fermiego. Kluczowe daty i nazwiska:

• 1933: Szilard wymyśla ideę, patent złożony 1934
• 1939: Frisch, Meitner, Bohr, Wheeler — rozszczepienie i jego fizyka
• 1940: Frisch, Peierls — memorandum o bombowej bombie uranowej (HEU), wyznaczenie masy krytycznej ~1 kg (błędne, ale w prawidłowym rzędzie wielkości)
• 1941: Glenn Seaborg odkrywa Pu-239 jako alternatywny materiał rozszczepialny
• 1942: CP-1 — doświadczalne potwierdzenie $k > 1$ w cywilnym układzie
• 1942–1945: reaktory Hanford (B, D, F) — przemysłowa produkcja Pu-239
• 1943: reaktor X-10 w Oak Ridge (graphite pile) — pierwsza produkcja plutonu w skali laboratoryjnej
• 1944–1945: Y-12 (wzbogacalnia elektromagnetyczna U-235) i K-25 (dyfuzyjna) — HEU dla Little Boy
• 1945: testy Trinity (Pu-239), Little Boy (HEU), Fat Man (Pu-239)

Fermi w 1954 roku jako jeden z pierwszych laureatów otrzymał nagrodę Komisji Energii Atomowej (późn. Fermi Award), przyznawan przez Departament Energii USA naukowcom wybitnie zasłużonym w dziedzinie energii jądrowej. Szilard otrzymał ją w 1959 roku, Teller w 1962, Bethe w 1993.⁴

Współczynnik mnożenia w zastosowaniach medycznych

Koncepcja współczynnika mnożenia ma niespodziewane zastosowanie w radioterapii neutronowej: BNCT (Boron Neutron Capture Therapy). W tej metodzie leczenia nowotworów pacjent otrzymuje związek boru-10 (np. BPA — boronofenyloalaninę), który preferencyjnie gromadzi się w komórkach nowotworowych. Następnie napromieniowuje się obszar nowotworu neutronami termicznymi. Reakcja $^{10}B(n,\alpha)^7Li$ wyzwala cząstki alfa i jony litu — o krótkim zasięgu (~10 μm), dokładnie zabijające komórkę nowotworową.

Związek z $k_{eff}$: reaktory badawcze dostarczające neutronów do BNCT muszą pracować stabilnie przy $k_{eff} = 1$. Pierwsze kliniczne zastosowania BNCT prowadzono przy reaktorze JRR-4 w Japonii (zamknięty 2006) i MIT Research Reactor (zamknięty 2018). Współcześnie opracowywane są kompaktowe akceleratorowe źródła neutronów (ABNS), eliminujące potrzebę reaktora — co czyni BNCT dostępnym w szpitalach. Instalacje kliniczne: STELLA (Japonia), IAEA-wspierane projekty w Argentynie i Chinach.⁶

Znaczenie historyczne i polityczne

Reakcja łańcuchowa jest jednym z tych odkryć naukowych, które bezpośrednio ukształtowały historię polityczną świata. Bez zrozumienia $k_{eff}$ Projekt Manhattan nie mógłby zaistnieć — bo nie byłoby wiadomo, ile materiału potrzeba, jak go ułożyć i kiedy nastąpi eksplozja.

Ale dokładnie te same równania, które opisują eksplozję atomową, opisują działanie reaktorów energetycznych. Podział jest pojęciowy, nie fizyczny: ten sam $k_{eff}$ może oznaczać stabilną pracę elektrowni (gdy $= 1$ przez lata) lub gwałtowną eksplozję (gdy $\gg 1$ przez mikrosekundy). Różnica tkwi w czasie i w kontroli.

To dlatego opanowanie reakcji łańcuchowej nie jest wyłącznie militarnym faktem historycznym. Jest fundamentem wszystkich 440 działających reaktorów energetycznych na świecie, które w 2023 roku dostarczyły ~10% globalnej elektryczności. Jest też podstawą przyszłych reaktorów fuzyjnych (choć tam łańcuchem jest fuzja, a nie rozszczepienie), transmutacji odpadów i medycznych zastosowań radioizotopów.

Polska planuje budowę swoich pierwszych elektrowni jądrowych (Choczewo, Pątnów) z reaktorami AP1000 (Westinghouse) lub APR1400 (Korea). Rozumienie $k_{eff}$, kinetyki reaktora i bezpieczeństwa krytycznościowego jest warunkiem wstępnym kształcenia kadry dla tych instalacji — i jest bezpośrednim powodem istnienia tej wiki.

Podsumowanie: bilans neutronów jako fundament

Cała fizyka jądrowa aplikacyjna — reaktory energetyczne, bronie jądrowe, bezpieczeństwo krytycznościowe, systemy ADS, ochrona przed proliferacją, pomiary MAEA, transmutacja odpadów i radioterapia BNCT — sprowadza się do jednego pytania: czy w danym układzie jeden neutron wywołuje więcej czy mniej niż jedno następne rozszczepienie?

Jeśli mniej ($k_{eff} < 1$): reakcja gaśnie.
Jeśli dokładnie jeden ($k_{eff} = 1$): stabilna praca reaktora.
Jeśli więcej ($k_{eff} > 1$): wykładniczy wzrost mocy — od kontrolowanego rozruchu reaktora aż po eksplozję atomową, zależnie od tego, o ile więcej i na jak długo.

Czynniki wpływające na $k_{eff}$ — masa i geometria, gęstość, skład materiałowy, obecność moderatora i reflektora, temperatura, stężenie trucizn neutronowych — są tematami powiązanych artykułów tej wiki. Zrozumienie każdego z nich oznacza głębsze zrozumienie jednego parametru: $k$.

Artykuły powiązane pozwalają zbudować pełny obraz: masa krytyczna opisuje geometrię i gęstość; dyfuzja neutronów opisuje ich przestrzenne zachowanie; neutrony opóźnione tłumaczą, dlaczego reaktorami można w ogóle sterować; przekrój czynny łączy fizykę jądrową z bilansem neutronów. Współczynnik mnożenia jest syntezą tych wszystkich zjawisk w jeden parametr operacyjny.

Dodatkowe materiały multimedialne

W praktyce ten temat spina się bezpośrednio z masą krytyczną i jej parametrami geometrycznymi, neutronami natychmiastowymi i opóźnionymi oraz Chicago Pile-1. Te trzy teksty pokazują odpowiednio geometryczną, czasową i historyczno-reaktorową stronę tego samego bilansu neutronowego.

Powiązane kalkulatory i narzędzia

• k_eff — pokazuje, jak geometria, moderator i straty neutronów wpływają na krytyczność układu.
• Wizualizacja: Wizualizator reakcji łańcuchowej — Siatka jąder rozszczepialnych, neutrony aktywne, rozszczepienia i przejścia między stanem podkrytycznym, krytycznym oraz nadkrytycznym.

Ćwiczenia praktyczne

Ćwiczenie laboratoryjne powinno mieć postać komputerowej symulacji Monte Carlo rozgałęzionej reakcji łańcuchowej. Minimalny wariant można wykonać we własnym kodzie w Pythonie, Julia lub MATLAB-ie, a wariant bardziej zaawansowany oprzeć na gotowym środowisku do modelowania stochastycznego. Kluczowe jest nie samo rysowanie wykresów, lecz rozdzielenie trzech procesów: rozszczepienia, pochłaniania bez rozszczepienia i ucieczki neutronów z układu. Dla każdego neutronu w pokoleniu należy więc losować jego dalszy los na podstawie zadanych prawdopodobieństw, a następnie liczyć populację neutronów, liczbę rozszczepień i estymator $k_{eff}$ w kolejnych generacjach. W części obowiązkowej trzeba porównać co najmniej trzy konfiguracje: podkrytyczną, krytyczną i nadkrytyczną, a następnie wykonać analizę czułości na zmianę prawdopodobieństwa ucieczki i pochłaniania.

Jeżeli laboratorium ma mocniejsze zaplecze obliczeniowe, sensownym rozszerzeniem jest model przestrzenny: sferyczny układ z prostym opisem długości drogi swobodnej neutronów i warunkiem ucieczki po przekroczeniu granicy geometrycznej. Taka wersja nie odtwarza pełnego transportu neutronów, ale już wnosi coś naukowego: pozwala zmierzyć, jak sam promień układu i przyjęta długość drogi swobodnej wpływają na estymowane $k_{eff}$. Wynikiem nie powinien być tylko kod, lecz także krótki raport z tabelą parametrów wejściowych, wykresami zmian populacji oraz dyskusją, które uproszczenia modelu najbardziej zniekształcają fizykę problemu.

Drugie ćwiczenie, teoretyczno-przemysłowe, powinno polegać na przeliczeniu, jak duży zapas reaktywności trzeba przewidzieć w układzie, który w stanie roboczym ma utrzymywać $k_{ef} = 1$, mimo że w czasie pracy pogarsza się bilans neutronowy przez wypalanie paliwa i gromadzenie produktów rozszczepienia. Punktem wyjścia jest wzór na reaktywność:

$$\rho = \frac{k_{ef} - 1}{k_{ef}}$$

opisany przez Słowińskiego.⁵ Należy więc:

1. założyć początkowe $k_{ef}$, na przykład $1{,}05$,
2. przeliczyć odpowiadającą mu reaktywność,
3. założyć spadek $k_{ef}$ w kolejnych etapach pracy,
4. wyznaczyć, w którym momencie układ schodzi do $k_{ef} = 1$,
5. policzyć, jaki margines sterowania musi zapewnić moderator, geometria i układ prętów pochłaniających.

Na skalę przemysłową to ćwiczenie nie prowadzi do budowy reaktora, lecz do oszacowania, jak precyzyjnie trzeba zarządzać bilansem neutronów, aby duży układ nie przechodził samorzutnie ani w stan zbyt głęboko podkrytyczny, ani w niebezpieczny dodatni zapas reaktywności. To właśnie tutaj widać, że pojęcie współczynnika mnożenia nie jest abstrakcją, ale parametrem projektowym, który musi być stale mierzony, modelowany i korygowany.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Masa krytyczna i jej parametry geometryczne
• Chicago Pile-1 (CP-1)
• Neutrony natychmiastowe a opóźnione
• Termiczna dyfuzja neutronów