Termiczna dyfuzja neutronów

Streszczenie

W reaktorze neutron nie leci po jednej prostej od rozszczepienia do następnego rozszczepienia. Zderza się wielokrotnie z jądrami moderatora, zmienia kierunek, traci energię i tylko statystycznie przesuwa się z miejsc, gdzie jego gęstość jest większa, do miejsc, gdzie jest mniejsza. Ten właśnie zbiorowy efekt opisuje dyfuzja neutronów. Bez niego nie da się zrozumieć, dlaczego rdzeń musi mieć określone rozmiary, po co stosuje się reflektor i czemu jedne moderatory pozwalają na pracę na uranie naturalnym, a inne nie.¹^,²

Termiczna dyfuzja neutronów dotyczy już neutronów spowolnionych, czyli takich, które zeszły z zakresu energii rozszczepieniowych do energii cieplnych. To właśnie one są najważniejsze dla większości klasycznych reaktorów. Ich rozkład przestrzenny decyduje o tym, gdzie reakcja łańcuchowa jest podtrzymywana najskuteczniej, a gdzie neutrony są tracone przez ucieczkę albo pochłonięcie.¹^,²^,³

Rozszerzenie tematu

Dyfuzja neutronów jest przybliżonym opisem transportu cząstek w ośrodku, w którym neutrony wielokrotnie rozpraszają się na jądrach i poruszają się chaotycznie. Słowiński podkreśla, że stochastyczny charakter migracji neutronów sprawia, iż lokalna koncentracja neutronów i lokalny strumień mają różne wartości w różnych miejscach rdzenia.¹ To właśnie z tego wynika makroskopowe zjawisko przepływu neutronów z obszarów o większej gęstości strumienia do obszarów o mniejszej.

W najprostszym przybliżeniu prąd dyfuzyjny neutronów opisuje prawo Ficka:

$$\vec{J}(\vec{r}) = -D(\vec{r}) \nabla \Phi(\vec{r})$$

gdzie $\vec{J}$ jest prądem neutronów, $\Phi$ skalarną gęstością strumienia, a $D$ współczynnikiem dyfuzji zależnym od właściwości materiału.¹ Znak minus oznacza, że neutrony „spływają” w dół gradientu strumienia, czyli z miejsc bardziej nasyconych do słabiej nasyconych.

Ważne jest jednak, by odróżnić dwa etapy drogi neutronu. Najpierw neutron rozszczepieniowy musi zostać spowolniony, zwykle z energii rzędu 1-2 MeV do energii cieplnych. Dopiero potem można sensownie mówić o jego termicznej dyfuzji w moderatorze i w rdzeniu. Słowiński rozróżnia tu parametry odpowiadające spowalnianiu oraz parametry odpowiadające późniejszej migracji neutronów cieplnych.¹ W praktyce projektowej oba etapy są ze sobą splecione, ale nie są tym samym zjawiskiem.

Jednym z najważniejszych parametrów jest długość dyfuzji neutronów cieplnych $L$, zdefiniowana przez zależność:

$$L^2 = \frac{D_T}{\Sigma_{a,T}}$$

gdzie $D_T$ oznacza współczynnik dyfuzji neutronów cieplnych, a $\Sigma_{a,T}$ makroskopowy przekrój czynny na ich pochłanianie.¹ Długość dyfuzji mówi w przybliżeniu, na jakiej skali strumień neutronów termicznych maleje wykładniczo wskutek dyfuzji i pochłaniania. Im większe $L$, tym dalej neutrony cieplne mogą „dojść”, zanim zostaną utracone.

W praktyce bardzo ważna jest też długość migracji $M$, która łączy w jednym parametrze etap spowalniania i późniejszej dyfuzji.¹ To ona pozwala bardziej realistycznie oszacować, jak daleko od miejsca narodzin neutron może dotrzeć, zanim albo spowoduje rozszczepienie, albo zostanie pochłonięty, albo ucieknie poza granice rdzenia.

Słowiński podaje konkretne wartości dla typowych moderatorów. Dla wody lekkiej długość dyfuzji neutronów cieplnych wynosi około 0,03 m, dla ciężkiej wody około 1,6 m, a dla grafitu około 0,52 m.¹ Te liczby pokazują od razu, jak różne własności materiałowe przekładają się na skalę geometryczną reaktora. Woda lekka moderuje świetnie, ale stosunkowo silniej pochłania neutrony. Ciężka woda pochłania ich bardzo mało, dzięki czemu daje bardzo duże długości dyfuzji i bardzo korzystny bilans neutronowy. Grafit wypada pośrodku: jest słabszym moderatorem objętościowo niż woda, ale znacznie mniej „zjada” neutrony niż zwykła woda.¹^,³

To właśnie z tego biorą się duże konsekwencje technologiczne. Reaktor moderowany ciężką wodą może pracować na uranie naturalnym, bo neutrony są spowalniane przy relatywnie małych stratach absorpcyjnych.³ Reaktor grafitowy też może dojść do korzystnego bilansu, ale wymaga odpowiednio czystego moderatora i odpowiednio dużych rozmiarów geometrycznych, aby ograniczyć straty powierzchniowe. Dlatego jakość grafitu reaktorowego była historycznie tak ważna.

Dyfuzja termiczna ma bezpośredni związek z problemem ucieczki neutronów z rdzenia. Słowiński opisuje cykl neutronowy tak, że tylko część neutronów po spowolnieniu dociera do zakresu cieplnego, a spośród tych, które już tam dotarły, tylko pewna część wywołuje dalsze rozszczepienia; pozostałe są pochłaniane przez moderator albo uciekają poza zewnętrzną granicę układu.² Z punktu widzenia projektu reaktora oznacza to, że geometria nie jest dodatkiem do fizyki materiałowej, lecz jej integralną częścią.

Stąd bierze się rola reflektora neutronów. Jeśli wokół strefy aktywnej znajduje się warstwa materiału silnie rozpraszającego i słabo pochłaniającego, część neutronów zmierzających ku granicy może zostać zawrócona do środka.⁴ W języku teorii dyfuzji oznacza to poprawę warunków brzegowych i zmniejszenie efektywnych strat powierzchniowych. Reflektor nie tworzy neutronów, ale sprawia, że większa część już istniejących neutronów pozostaje użyteczna dla układu.

Warto też pamiętać, że teoria dyfuzji jest przybliżeniem. Zakłada ona brak silnej anizotropii pola neutronowego i sensownie działa przede wszystkim tam, gdzie neutrony wykonały już wiele zderzeń i ich rozkład kierunkowy jest stosunkowo wygładzony. Dlatego jest znacznie lepsza dla neutronów termicznych w dużym moderatorze niż dla neutronów bardzo prędkich albo w małych, silnie niejednorodnych układach. Mimo to pozostaje jednym z podstawowych narzędzi opisu rdzenia reaktora, bo daje użyteczny związek między materiałem, geometrią i rozkładem strumienia.

To również wyjaśnia, dlaczego temat dyfuzji termicznej nie ma takiego samego znaczenia dla ładunku jądrowego jak dla reaktora. W bombie liczy się bardzo szybkie mnożenie na neutronach prędkich, a czas dostępny na moderację i wielokrotne rozproszenia jest zbyt mały. W reaktorze sytuacja jest odwrotna: duża część projektu polega na takim prowadzeniu neutronów przez moderator, aby jak najwięcej z nich wróciło do paliwa już jako neutrony termiczne, a nie zostało straconych przez ucieczkę albo niepożądany wychwyt.²^,³

Najkrótsze podsumowanie jest takie: termiczna dyfuzja neutronów opisuje przestrzenny ruch już spowolnionych neutronów w moderatorze i w rdzeniu. Parametry takie jak $D$, $L$ i $M$ nie są abstrakcyjnymi stałymi z podręcznika, tylko wprost przekładają się na wybór moderatora, wielkość rdzenia, sens stosowania reflektora i możliwość utrzymania reakcji łańcuchowej w rzeczywistym reaktorze.¹^,²^,⁴

Historia teorii dyfuzji neutronów

Teoria dyfuzji neutronów rozwinęła się niemal równocześnie z odkryciem samego neutronu przez Jamesa Chadwicka w 1932 roku. Już rok później Enrico Fermi w Rzymie zaczął badać, jak różne materiały wpływają na spowolnienie neutronów i ich pochłanianie. Słynna obserwacja z 1934 roku — że neutrony opóźnione przez hydrogeniczny blok parafiny powodują silniejszą aktywację srebrnej folii niż neutrony prędkie — zapoczątkowała systematyczne badania moderatorów i spowalniania neutronów.

Formalny aparat matematyczny dyfuzji neutronów wywodzi się z teorii kinetycznej gazów i teorii transportu ciepła. Fundamentalne równanie transportu Boltzmanna było znane od XIX wieku; jego zastosowanie do neutronów w reaktorze stało się możliwe dzięki pracom m.in. Placzka i Wignera w latach trzydziestych. Przybliżenie dyfuzyjne — uproszczenie równania Boltzmanna do równania Ficka z dodatkowymi wyrazami źródłowymi i pochłaniającymi — wprowadzono jako narzędzie projektowe reaktorów w ramach Projektu Manhattan (1942–1945).⁵

Na Uniwersytecie Chicago, a następnie w Oak Ridge i w Hanford, fizycy Projektu Manhattan (m.in. Fermi, Wigner, Wheeler, Teller, Seaborg) opracowywali coraz precyzyjniejsze modele rdzenia reaktora. Pierwsze obliczenia projektowe Chicago Pile-1 (CP-1) były prowadzone właśnie przy użyciu uproszczonej teorii dyfuzji neutronów — Fermi i jego współpracownicy obliczyli, ile bloków grafitu i ile naturale-uranowych wkładek potrzeba, by osiągnąć stan krytyczny. Okazało się, że obliczenia były bardzo bliskie rzeczywistości. CP-1 osiągnął krytyczność 2 grudnia 1942 roku przy dokładnie wyliczonej masie paliwa.

Po wojnie teoria dyfuzji neutronów stała się standardową częścią inżynierii reaktorowej na całym świecie. Podręczniki Glasstone'a i Edlunda (The Elements of Nuclear Reactor Theory, 1952) oraz Weinberga i Wignera (The Physical Theory of Neutron Chain Reactors, 1958) skodyfikowały aparat matematyczny na kilka dekad. Do dziś studenci fizyki reaktorowej uczą się dyfuzji neutronów z tych samych fundamentów.

Pełne równanie dyfuzji neutronów

Równanie Ficka w postaci $\vec{J} = -D\nabla\Phi$ jest tylko jednym składnikiem. Pełne równanie bilansu neutronów cieplnych w stanie ustalonym (steady-state) wygląda następująco:

$$-D\nabla^2\Phi(\vec{r}) + \Sigma_a\,\Phi(\vec{r}) = S(\vec{r})$$

gdzie:

$D$ — współczynnik dyfuzji [m]
$\Phi(\vec{r})$ — strumień neutronów cieplnych [$\mathrm{n \cdot m^{-2} \cdot s^{-1}}$]
$\Sigma_a$ — makroskopowy przekrój czynny na pochłanianie neutronów cieplnych [$\mathrm{m^{-1}}$]
$S(\vec{r})$ — źródło neutronów cieplnych (z procesu spowalniania neutronów prędkich) [$\mathrm{n \cdot m^{-3} \cdot s^{-1}}$]

Pierwszy wyraz opisuje straty dyfuzyjne (ucieczkę), drugi wyraz straty na pochłanianie, a prawa strona to źródło. Przekształcając:

$$\nabla^2\Phi - \frac{1}{L^2}\Phi = -\frac{S}{D}$$

gdzie $L^2 = D/\Sigma_a$ jest kwadratem długości dyfuzji.

Dla reaktora krytycznego w przybliżeniu jednogrupoawym źródło neutronów pochodzi z rozszczepień, a strumień neutronów cieplnych musi być samopodtrzymujący. Prowadzi to do równania własnego:

$$\nabla^2\Phi + B^2\Phi = 0$$

gdzie $B^2$ jest kwadratu całkowitego beklingu (buckling) geometrycznego. Dla kuli o promieniu $R$ rozwiązanie ma postać:

$$\Phi(r) = A\,\frac{\sin(B r)}{r}$$

z warunkiem brzegowym $\Phi(\tilde{R}) = 0$ przy ekstrapolowanym promieniu $\tilde{R} = R + d$ (gdzie $d \approx 2D$ jest długością ekstrapolacji). To prowadzi do warunku krytyczności: $B = \pi/\tilde{R}$, stąd minimalna masa krytyczna zależy od $L$ i $M$ materiałów.⁵

W przypadku geometrii cylindrycznej rdzeń reaktora (niski cylinder) daje odrębne rozwiązanie przez produkty funkcji Bessela i trygonometrycznych. Obliczenia te stały się standardowymi ćwiczeniami projektowymi reaktorów.

Spowalnianie neutronów — wiek Fermiego

Nim neutron dotrze do zakresu cieplnego, musi przejść przez etap spowalniania (thermalisation, moderation). Formalnym opisem tego procesu jest teoria wieku Fermiego (Fermi age theory), w której definiuje się wiek $\tau$ jako miarę dystansu, jaki neutron pokonuje podczas spowalniania:

$$\tau = \int_{E_T}^{E_0} \frac{D(E)}{\xi\,\Sigma_s(E)\,E}\,dE$$

gdzie $\xi$ jest logarytmicznym dekrementem energii na jedno zderzenie (dla danego pierwiastka moderatora), $\Sigma_s$ makroskopowym przekrojem czynnym na rozpraszanie.

Dla różnych moderatorów wartości $\tau$ różnią się dramatycznie:

Moderator	$\tau$ [$\mathrm{cm^2}$]	$L^2$ [$\mathrm{cm^2}$]	$M^2 = L^2 + \tau$
Woda lekka H₂O	27	8	35
Ciężka woda D₂O	131	20400	20531
Grafit (czysty)	364	3500	3864
Beryl BeO	172	450	622

Woda lekka ma bardzo mały wiek Fermiego ($\tau = 27\,\mathrm{cm^2}$), bo wodór znakomicie spowalnia neutrony w małej liczbie zderzeń. Ale ma też mały $L^2$, bo H pochłania neutrony termiczne silniej niż D lub węgiel-12. Ciężka woda ma duże $\tau$ i ogromne $L^2$ — neutrony przemierzają przez nią olbrzymie dystanse, zanim zostaną pochłonięte, co daje znakomity bilans neutronowy.

Długość migracji jest zdefiniowana przez:

$$M^2 = L^2 + \tau$$

i reprezentuje uogólnioną miarę całkowitego dystansu od narodzin neutronu rozszczepialnego do jego pochłonięcia lub ucieczki. $M$ jest tym parametrem, który bezpośrednio wiąże się z geometrycznym rozmiarem krytycznym rdzenia.⁵

Porównawcze właściwości moderatorów

Dobry moderator musi spełniać kilka warunków jednocześnie: efektywnie spowalniać neutrony (duże $\xi$, małe $A$), minimalnie je pochłaniać (mały $\sigma_a$) i być trwały mechanicznie oraz chemicznie w warunkach promieniowania.

Wodór i woda lekka (H₂O)

Wodór jest najefektywniejszym spowalniaczem pojedynczych neutronów — w każdym zderzeniu centrycznym (frontalnym) neutron może oddać nawet całą energię kinetyczną. Logarytmiczny dekrement energii $\xi_H = 1{,}000$ jest najwyższy ze wszystkich pierwiastków. Średnio potrzeba tylko ok. 18 zderzeń, by spowolnić neutron z 2 MeV do energii termicznych.

Wadą wodoru jest stosunkowo wysoki przekrój czynny na wychwyt neutronów (~0,33 barn dla termicznych). To sprawia, że w reaktorach H₂O jako moderator wymaga wzbogaconego uranu (3–5% U-235) — strat neutronowych jest zbyt wiele, by utrzymać krytyczność na uranie naturalnym.

Woda lekka jest najtańszym moderatorem i chłodziwem, dominującym w większości reaktorów komercyjnych (PWR, BWR). Mały rozmiar geometryczny rdzenia (mały $M^2$) jest zaletą inżynierską i ekonomiczną.⁵

Ciężka woda (D₂O)

Deuter ma $\xi_D = 0{,}725$ — gorszy od wodoru jako spowalniacz — ale za to przekrój czynny na wychwyt neutronów jest ponad 600 razy mniejszy ($\sigma_a \approx 0{,}0005\,\mathrm{barn}$). Neutron może przebywać w D₂O bardzo długo i przebyć olbrzymi dystans, zanim zostanie pochłonięty. Stąd $L^2 = 20\,400\,\mathrm{cm^2}$, co przekłada się na $L \approx 143\,\mathrm{cm}$ — i to jest właśnie powód, dla którego reaktory CANDU pracują na paliwie niezbogaconym.

Wadą D₂O jest cena — woda ciężka kosztuje ok. 350–400 USD/kg (cena historycznie wahała się bardzo). Produkcja D₂O to skomplikowany przemysłowy proces — metodą GS (Girdler–Sulfide) lub elektrolityczną — opisany w odrębnych artykułach. Dla porównania: woda lekka jako moderator kosztuje niemal nic.⁵

Grafit

Węgiel-12 ma $\xi_C = 0{,}158$ (potrzeba ok. 115 zderzeń dla moderacji do termicznych). Przekrój czynny na wychwyt jest mały, ale nie tak mały jak deuteru: $\sigma_a \approx 0{,}0034\,\mathrm{barn}$. Grafit może pracować jako moderator na uranie naturalnym pod warunkiem, że jest odpowiednio czysty — nawet kilka ppm boru całkowicie niszczy bilans neutronowy, bo bor ma gigantyczny przekrój czynny na wychwyt neutronów termicznych (~3840 barn).

Reaktory grafitowe dominowały historycznie w ZSRR (seria RBMK), Wielkiej Brytanii (Magnox, AGR) i USA (reaktory Hanford, X-10 w Oak Ridge). Wymagają jednak dużych rozmiarów geometrycznych ze względu na gorsze właściwości spowalniające. Reaktor CP-1 z 1942 roku ważył ponad 400 ton.

Beryl i BeO

Beryl ($^9$Be) i tlenek berylu (BeO) są interesującymi moderatorami dla specjalnych zastosowań. $\xi_{Be} = 0{,}209$, $\sigma_a \approx 0{,}0092\,\mathrm{barn}$. BeO ma lepszy bilans neutronowy niż grafit i jest stosowany w reaktorach badawczych o wysokim strumieniu oraz w niektórych projektach reaktorów kosmicznych (SNAP). Beryl emituje też neutrony pod wpływem napromieniowania gamma (fotoneumtrony), co może być zaletą lub wadą zależnie od zastosowania.

Wadą berylu jest toksyczność — pary BeO i pyły berylu są silnie rakotwórcze, co utrudnia obróbkę i nakłada surowe wymagania BHP.

Wartości liczbowe i konsekwencje projektowe

Dla projektanta reaktora kluczowym pytaniem jest, jaki minimalny rozmiar rdzenia jest konieczny, by osiągnąć krytyczność, i jak ten rozmiar zależy od moderatora i reflektora.

Warunek krytyczności bez reflektora dla jednogrupaowej teorii dyfuzji:

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

gdzie $k_\infty$ jest mnożnikiem neutronów w nieskończonym układzie, a $B^2$ jest beklingu geometrycznym (dla kuli $B^2 = (\pi/R)^2$). Stąd minimalny promień krytyczny kuli bez reflektora:

$$R = \frac{\pi}{B_c},\quad B_c^2 = \frac{k_\infty - 1}{M^2}$$

Większy $M^2$ (np. ciężka woda) wymaga mniejszego $B^2$, czyli mniejszego $B$, czyli większego rdzenia geometrycznie. Ale jednocześnie większy $M^2$ oznacza korzystniejszy bilans neutronowy ($k_\infty$ wyższe), więc mianownik ($k_\infty - 1$) jest większy i ostateczny rdzeń wcale nie musi być ogromny.

Konkretnie: reaktor D₂O moderowany może pracować na uranie naturalnym (0,71% U-235), bo straty neutronowe są znikome. Reaktor H₂O moderowany wymaga ~3% U-235, ale geometrycznie jest znacznie mniejszy. Reaktor grafitowy wymaga ok. 1–1,5% U-235 lub bardzo dużej masy czystego grafitu, stąd historyczne olbrzymy w stylu Hanford.

Reflektor zmienia warunki brzegowe i pozwala na efektywne zmniejszenie masy krytycznej. Gruba warstwa D₂O jako reflektora dla rdzenia uranowego zmniejsza minimalną masę krytyczną o 40–60%. To właśnie dlatego reaktory i — w mniejszym stopniu — ładunki jądrowe z tamperem mają mniejsze masy krytyczne niż układy nagie.⁵

Teoria transportu vs przybliżenie dyfuzyjne

Dyfuzja neutronów to przybliżenie. Pełna teoria jest opisana przez równanie transportu Boltzmanna dla neutronów:

$$\frac{1}{v}\frac{\partial\Psi}{\partial t} + \hat{\Omega}\cdot\nabla\Psi + \Sigma_t\Psi = \int\!\!\int \Sigma_s(\hat{\Omega}'\to\hat{\Omega},\,E'\to E)\Psi\,dE'\,d\hat{\Omega}' + S$$

gdzie $\Psi(\vec{r}, \hat{\Omega}, E, t)$ jest kątową fluencją neutronów — funkcją pozycji, kierunku, energii i czasu. Rozwiązanie tej równania (niezależnie od przybliżeń) daje pełny opis pola neutronowego we wszystkich geometriach.

Przybliżenie dyfuzyjne (P1 lub B1) powstaje przez rozwinięcie $\Psi$ w szereg sferycznych harmonicznych i zatrzymanie na pierwszym niebanalnym wyrazie. Jest ono poprawne tam, gdzie pole neutronowe jest niemal izotropowe — co zachodzi daleko od granic, silnych absorbentów i interfejsów materiałowych. Przy granicy rdzeń-chłodziwo lub rdzeń-reflektor przybliżenie dyfuzyjne jest mniej precyzyjne.

Dla precyzyjnych obliczeń projektowych reaktorów komercyjnych używa się dziś:

Deterministyczne kody transportu: DORT, PARTISN, DENOVO — rozwiązują dyskretyzowane równanie transportu Boltzmanna metodą $S_N$ (discrete ordinates).
Monte Carlo: MCNP (Los Alamos), OpenMC (MIT), SERPENT (Finlandia) — symulują historia każdego neutronu statystycznie, nie wymagają żadnych przybliżeń transportowych.
Kody dyfuzyjne: CITATION, PARCS, SIMULATE — używane do cyklicznych obliczeń rdzenia, bo są szybkie; mniej dokładne przy zagadnieniach związanych z granicami.

W Polsce obliczenia kinetyki reaktorowej prowadzone są m.in. w Instytucie Energii Atomowej NCBJ w Świerku, który posiada reaktor badawczy MARIA — moderowany lekką wodą, z berylowym reflektorem. Obliczenia MARIA prowadzone są z użyciem kodów transportu, uzupełnianych pomiarami aktywacyjnymi i spektrometrią.

Reaktory ciężkowodne — CANDU i NRX

Reaktor CANDU (Canada Deuterium Uranium), opracowany przez Atomic Energy of Canada Limited (AECL) w latach 1950–1971, jest najlepszym przykładem praktycznego zastosowania wielkiej długości dyfuzji w D₂O.

W CANDU moderator (D₂O) jest fizycznie oddzielony od chłodziwa (D₂O lub H₂O w kanalach) i utrzymywany pod niskim ciśnieniem w dużej, cylindrycznej kadzi kalandrii (calandria). Umożliwia to wymianę paliwa „na gorąco" bez wyłączania reaktora, co jest unikalną zaletą. Paliwo jest naturalne UO₂ w pakietach (fuel bundles) o długości 49,5 cm.

Skuteczność CANDU wynika bezpośrednio z fizyki dyfuzji: $L_{D2O} \approx 143\,\mathrm{cm}$ oznacza, że nawet przy stosunkowo rzadkim ułożeniu prętów paliwowych (co kilka centymetrów) neutrony termiczne swobodnie migrują przez moderator do sąsiednich prętów. Bilans jest niemal zerowy przy pochłanianiu przez D₂O.

Poprzednikiem CANDU był reaktor NRX (National Research Experimental) w Chalk River (Kanada), uruchomiony w 1947 roku — jeden z pierwszych reaktorów ciężkowodnych. W 1952 roku NRX przeżył poważny wypadek — zapomniany przez szeroką publiczność incydent, który przez wiele dekad pozostawał jednym z poważniejszych wypadków jądrowych na świecie. Rdzeń reaktora uległ częściowemu uszkodzeniu termicznemu, a kontaminacja radioaktywna wymagała długotrwałego oczyszczania. Spośród studentów przybyłych do pomocy przy oczyszczaniu był Jimmy Carter — późniejszy prezydent USA, jeden z nielicznych amerykańskich prezydentów ze znajomością inżynierii jądrowej.⁵

Długość dyfuzji a bezpieczeństwo i proliferacja

Parametr $L$ ma też znaczenie z punktu widzenia bezpieczeństwa i nieproliferacji. Ciężka woda jako moderator jest towarem kontrolowanym przez Grupę Dostawców Jądrowych (NSG — Nuclear Suppliers Group), bo jej zastosowanie w reaktorach CANDU pozwala na użycie naturalnego uranu, co de facto obniża barierę technologiczną wejścia do programów proliferacyjnych.

Irak pod rządami Saddama Husajna próbował w latach osiemdziesiątych zakupić D₂O od różnych dostawców (m.in. przez Niemcy i Jugosławię), by zasilić reaktor badawczy Osiraq. Izraelski atak na Osiraq w 1981 roku zatrzymał iracki program, zanim reaktor uruchomiono. Indie posiadają dużą flotę reaktorów CANDU (PHWR), a ich program jądrowy oryginalnie budował na kanadyjskiej technologii D₂O do produkcji plutonu dla testów jądrowych.

Pakistan zbudował reaktor ciężkowodny Khushab wbrew embargom, korzystając z D₂O wyprodukowanej w Pakistanie samodzielnie. Reaktor ten jest uznawany za kluczowy element programu produkcji plutonu dla broni jądrowej. Objętość moderatora i jego czytelne parametry dyfuzji są przedmiotem analiz wywiadu, bo pozwalają oszacować moc reaktora i produkowaną ilość plutonu.

Grafit jako moderator jest mniej objęty kontrolami, ale NSG ogranicza eksport wysokiej czystości grafitu reaktorowego (limit zanieczyszczeń borem 5 ppm), bo reaktory grafitowe mogą także produkować pluton (jak to miało miejsce w radzieckich reaktorach w Tomsku i Krasnojarsku do lat 90.).

Moderacja a przekrój czynny na wychwyt rezonansowy

Istotnym efektem, który komplikuje proste równania dyfuzji, jest wychwyt rezonansowy uranu-238. Podczas spowalniania neutron z 2 MeV musi przejść przez zakres energii 1 eV – 1 keV, w którym U-238 ma ogromne rezonanse absorpcyjne. Szansa, że neutron zostanie pochłonięty przez U-238 zanim zostanie spowolniony do energii termicznych, jest opisana przez prawdopodobieństwo uniknięcia rezonansu $p$ (resonance escape probability).

Parametr $p$ jest silnie zależny od stopnia jednorodności czy heterogeniczności paliwa i moderatora. W jednorodnym ośrodku (paliwo wymieszane z moderatorem) neutron spowalnia przy ciągłej bliskości do U-238 — $p$ jest małe. W układzie heterogenicznym (pręty paliwowe zanurzone w moderatorze) neutron spowalnia głównie w moderatorze, gdzie jest mało U-238, i wraca do paliwa dopiero jako neutron termiczny — $p$ jest duże.

To właśnie heterogeniczność rdzenia reaktora (pręty paliwowe + moderator) jest jedną z kluczowych innowacji Fermi-Szilarda opracowaną dla CP-1. Jednorodna mieszanina grafitu z uranem nigdy by nie osiągnęła krytyczności — straty przez rezonansowy wychwyt U-238 byłyby zbyt duże. Układ heterogeniczny (kostki uranu przesunięte względem bloków grafitu) poprawił $p$ na tyle, że krytyczność stała się możliwa.⁵

Termalizacja i efekty chemiczne wiązań

W zakresie energii poniżej ~1 eV (energia termiczna neutronów) proste klasyczne rozpraszanie nie jest już wystarczającym opisem. Neutrony o energii porównywalnej z energiami drgań sieciowych materiału wchodzą w oddziaływania z fonowymi modami drgań siatki. W wodzie atomów wodoru poruszają się w wyniku drgań molekularnych i rotacji — o energiach rzędu meV.

Efekt termalizacji chemicznej oznacza, że efektywny przekrój czynny na rozpraszanie neutronów przez wodę silnie zależy od temperatury i struktury molekularnej wody. Dane do tych obliczeń zawarte są w bibliotekach S(α,β) (thermal scattering matrices), zawartych w bazach danych ENDF/B, JEFF i JENDL.

W praktyce ważny efekt ma temperatura moderatora. Wzrost temperatury wody przesuwa rozkład Maxwella neutronów termicznych do wyższych energii — co zmienia efektywną wartość przekroju czynnego na pochłanianie przez U-235. W reaktorach H₂O wzrost temperatury moderatora powoduje też zmniejszenie gęstości i odpowiednio zmniejszenie mocy moderacji — to jest mechanizm ujemnego współczynnika temperatury moderatora, który jest ważną cechą bezpieczeństwa reaktorów PWR i BWR.

Dyfuzja w ładunkach jądrowych — tamper i reflektor

Choć dyfuzja termiczna neutronów jest zjawiskiem reaktorowym, analogiczne mechanizmy transportu i odbicia neutronów są ważne w bombach jądrowych w postaci tampera (odbijacza neutronów). Tamper w ładunku jądrowym nie jest moderatorem — nie służy do spowalniania neutronów. Ale pełni analogiczną funkcję do reaktorowego reflektora: zawraca neutrony prędkie emitowane przez rozszczepienie z powrotem do materiału rozszczepialnego, zmniejszając efektywną masę krytyczną.

Materiały tamperowe to wolfram, uran naturalny (U-238), beryl lub stal. U-238 jako tamper jest szczególnie efektywny, bo łączy dwie funkcje: efektywnie odbija neutrony prędkie (duży przekrój na rozpraszanie sprężyste) i sam ulega rozszczepieniu szybkimi neutronami (próg rozszczepienia U-238 ~1 MeV), co dodatkowo zwiększa uzysk energetyczny ładunku.

Analiza ballistic properties tampera z U-238 w bombie jądrowej wymagała tych samych narzędzi matematycznych co obliczenia reaktorowe — równań dyfuzji / transportu, ale dla neutronów prędkich i z odpowiednimi sekcjami czynymi dla U-238. Prace te były prowadzone przez tych samych fizyków, którzy projektowali reaktory. Różnica leży wyłącznie w zakresie energii neutronów i skali czasowej zjawiska.

Pomiary parametrów dyfuzji — eksperymentalne metody

Wartości $D$, $L$ i $\tau$ dla różnych moderatorów mierzone są przez szereg technik eksperymentalnych:

Metoda impulsowa (pulsed neutron method): do bloku moderatora wprowadza się krótki impuls neutronów z generatora DD lub DT, a następnie obserwuje się zanik populacji neutronów termicznych w czasie. Tempo zaniku daje $\Sigma_a$ i $D$ oddzielnie.

Metoda mnożnika neutronów (exponential assembly): pomiar geometrycznego rozkładu strumienia neutronów w exponential assembly (uran + moderator) pozwala wyznaczyć $k_\infty$ i $L$ przez dopasowanie krzywej eksponencjalnej.

Aktywacja próbek (foil activation): wkładanie małych próbek aktywowanych materiałów (Au, In, Co) do moderatora i mierzenie ich aktywności po napromieniowaniu daje rozkład strumienia termicznego w przestrzeni.

Wszystkie te metody wymagają czystych, dużych próbek moderatora i starannej korekcji na skończone efekty geometryczne. Precyzja parametrów dla H₂O i D₂O jest dziś na poziomie 0,5–1%, dla grafitu — silnie zależna od czystości próbki.

Historycznie szczególnie ważne były pomiary w zespołach krytycznych (critical assemblies) — takich jak TCA (Tank Critical Assembly) w Oak Ridge czy OECD HECTOR w Niemczech. Zanim komputery pozwoliły na numeryczne rozwiązanie równań transportu, fizycy reaktorowi opierali się na empirycznych pomiarach „eksponencjalnych": stackach uranu i moderatora o różnych proporcjach, z których mierzono przestrzenny zanik strumienia neutronów. Parametry $k_\infty$, $L^2$ i $\tau$ wyznaczano drogą dopasowania krzywych eksponencjalnych do profili aktywacyjnych, ustawiając wkładki aktywacyjne (foils) w regularnej siatce i mierząc ich aktywację po napromieniowaniu określoną fluencją. Technika ta, mimo prostoty konceptualnej, wymagała bardzo dokładnych kalibracji detektorów i korekcji na efekty krawędziowe — w szczególności wtedy, gdy moderator nie był idealnie jednorodny lub gdy obecne były interfejsy między różnymi materiałami. Wyniki tych historycznych pomiarów z lat 1942–1970 nadal stanowią część baz walidacyjnych (benchmark libraries) używanych do certyfikacji współczesnych kodów obliczeniowych takich jak MCNP czy OpenMC.

Zanieczyszczenia w moderatorach — krytyczne znaczenie

Żaden moderator reaktorowy nie jest idealnie czysty. Praktyczne ograniczenia czystości są jednym z najważniejszych czynników projektowych:

Dla grafitu: bor (~5830 barn dla termicznych), wodór (wnikanie wilgoci → H₂O jako zanieczyszczenie), siarka. Już 1–2 ppm boru skutecznie pochłania neutony termiczne i niszczy bilans. Projekt CP-1 wymagał od Fermiego zamówienia specjalnie czystego grafitu od Union Carbide — standardowy grafit elektryczny był zbyt skażony borem.

Dla D₂O: H₂O jako zanieczyszczenie. Woda lekka ma $\sigma_a(H) \approx 0{,}33\,\mathrm{barn}$ wobec $\sigma_a(D) \approx 0{,}0005\,\mathrm{barn}$. Nawet 0,5% H₂O w moderatorze D₂O reaktora CANDU pogarsza bilans neutronowy. Reaktory CANDU wymagają D₂O o czystości >99,75% i regularnego izotopowego oczyszczania (upgrading) moderatora przez elektrolizę lub destylację.

Dla H₂O: chlor (z chlorowania wody), bor (popiół lotny), gad. Gadolin-155 i gadolin-157 mają przekroje na wychwyt ~60.000 barn — są stosowane jako burnable absorbers w reaktorach, ale jako zanieczyszczenie byłyby katastrofalne.

Dyfuzja neutronów w reaktorach badawczych wysokiego strumienia

Reaktory badawcze o wysokim strumieniu neutronów (High Flux Reactors — HFR) to specjalna kategoria urządzeń, gdzie termiczna dyfuzja neutronów jest projektowana w sposób szczególnie świadomy, bo cel jest odwrotny niż w reaktorze energetycznym: zamiast minimalizować ucieczkę neutronów do wyprowadzenia energii elektrycznej, maksymalizuje się ich gęstość w wybranym miejscu eksperymentalnym.

Reaktor ILL (Institut Laue-Langevin) w Grenoble we Francji osiąga strumień neutronów cieplnych $\Phi_{th} \approx 1{,}5 \times 10^{15}\,\mathrm{n \cdot cm^{-2} \cdot s^{-1}}$ — jeden z najwyższych na świecie. Rdzeń reaktora (20 MW, paliwo wzbogacone >90% U-235) jest bardzo mały (cylinder ~20 cm średnicy), otoczony 2,5-metrową kadź D₂O moderatora. Całe urządzenie zaprojektowane jest tak, aby neutony termiczne migrujące przez D₂O tworzyły olbrzymi, jednorodny rezerwuar, z którego wyprowadzane są przez kanały do kilkudziesięciu jednoczesnych eksperymentów neutronograficznych, dyfraktometrycznych i małokątowego rozpraszania.

Długość dyfuzji $L_{D2O} \approx 143\,\mathrm{cm}$ oznacza, że rdzeń ILL jest małym „punktem" w morzu neutronów dyfundujących przez dużą kadź D₂O — efektywna przestrzeń eksperymentalna jest znacznie szersza niż sam rdzeń. To jest forma technologiczna, wynikająca wprost z dużego $L$ ciężkiej wody.⁵

Reaktor MARIA w Świerku pod Warszawą (30 MW, uruchomiony 1974) działa na zasadzie nieco innej: jest basenowy, moderowany lekką wodą, z berylowym reflektorem. Strumień termiczny $\Phi_{th} \approx 3{,}5 \times 10^{13}\,\mathrm{n \cdot cm^{-2} \cdot s^{-1}}$. Przez kilkanaście kanałów napromieniowania prowadzone są eksperymenty z aktywacją próbek, produkcją radioizotopów (m.in. Mo-99 do celów medycznych) i badaniem struktury materiałów metodą dyfrakcji neutronów.

Neutronografia i dyfrakcja neutronów

Jedną z najważniejszych aplikacji termicznych neutronów jest neutronografia — technika analogiczna do rentgenografii, lecz oparta na rozpraszaniu neutronów przez jądra atomowe zamiast na rozpraszaniu elektronów przez rentgeny.

Ponieważ neutrony mają masę i skok kwantowy (de Broglie $\lambda = h/(mv)$), termiczne neutrony z reaktora mają długości fali w zakresie 1–3 Å — porównywalnej z odległościami między atomami w kryształach. W tzw. dyfrakcji neutronowej (neutron powder diffraction lub single-crystal neutron diffraction) wiązka neutronów pada na próbkę i ulega ugięciu Bragga, dając wzorzec dyfrakcyjny, z którego odczytuje się strukturę krystaliczną materiału.

Kluczową zaletą neutronów w porównaniu do rentgenów jest to, że neutrony rozpraszają się od jąder atomowych, a nie elektronów. Oznacza to, że:

Neutrony są czułe na lekkie atomy (H, Li, B), które są niemal niewidoczne dla rentgenów.
Neutrony rozróżniają izotopy — np. H i D mają różne długości rozpraszania neutronowego.
Neutrony mają moment magnetyczny i mogą rozpraszać się na momentach magnetycznych elektronów, pozwalając na badanie struktury magnetycznej materiałów.

Reaktory badawcze z wysokim strumieniem neutronów termicznych są niezbędne do neutronografii. Synchrotrony mogą zastąpić reaktory do rentgenografii, ale do dyfrakcji neutronów potrzebny jest reaktor lub spallacyjne źródło neutronów (jak ISIS w Rutherford Appleton Laboratory, Wielka Brytania, lub SNS w Oak Ridge).

Pośrednio zatem fizyka dyfuzji termicznej neutronów i parametry $L$, $M$ są bezpośrednio powiązane z naukami materiałowymi, farmacją (struktura białek), inżynierią i fizyką kondensacji materii.

Podsumowanie

Termiczna dyfuzja neutronów — opisana przez prawo Ficka, długość dyfuzji $L$, wiek Fermiego $\tau$ i długość migracji $M$ — jest fundamentem rozumienia fizyki reaktorowej. Pozwala ilościowo powiązać właściwości materiałowe moderatora z wymaganymi rozmiarami rdzenia, skutecznością reflektora i warunkiem krytyczności.

Różnica między wodą lekką ($L \approx 3\,\mathrm{cm}$) a ciężką ($L \approx 143\,\mathrm{cm}$) nie jest akademickim drobiazgiem — to ona rozstrzyga, czy reaktor może pracować na uranie naturalnym, ile kosztuje paliwo, jak duże musi być urządzenie i jakie ryzyko proliferacji stwarza technologia. Historia reaktorów jądrowych od CP-1 (1942) do CANDU (1971) do nowoczesnych SMR jest w dużym stopniu historią inżynierów szukających optymalnego kompromisu między współczynnikiem spowalniania, współczynnikiem pochłaniania i kosztami.

Fizyka dyfuzji ma też wymiar czysto naukowy: parametry $L$, $\tau$, $D$ mierzone w laboratoriach to nie tylko stałe inżynierskie, ale wkład do tablic danych jądrowych (ENDF/B, JEFF, JENDL), które służą zarówno do projektowania reaktorów i oceny bezpieczeństwa, jak i do badań fundamentalnych w astrofizyce nuklearnej i syntezy ciężkich pierwiastków. Bez precyzyjnych pomiarów parametrów dyfuzji neutronów termicznych nie byłoby możliwe ani zaprojektowanie CP-1 przez Fermiego, ani zaprojektowanie współczesnych reaktorów PWR czy reaktorów CANDU na uranie naturalnym, ani wreszcie obliczenie minimalnej masy krytycznej dla celów obrony lub kontroli zbrojeń. To jest jedna z tych fizycznych wielkości, która łączy naukę podstawową z inżynierią i polityką nieproliferacji w jeden nierozerwalny węzeł.¹^,²^,⁴^,⁵

Dodatkowe materiały multimedialne

Do tego artykułu nie dodano jeszcze materiałów wideo. Warto wrócić do tej sekcji dopiero wtedy, gdy uda się znaleźć materiał dobrze pokazujący różnicę między samym spowalnianiem neutronów a ich późniejszą dyfuzją termiczną w rdzeniu.

Najlepiej zestawić ten tekst z artykułami o przekroju czynnym na rozszczepienie, ciężkiej wodzie jako moderatorze oraz grafticie reaktorowym i zanieczyszczeniach borem. Tam widać, jak równania dyfuzji przekładają się na rzeczywisty wybór moderatora i na sukces albo porażkę programu reaktorowego.

Powiązane kalkulatory i narzędzia

Dyfuzja i spowalnianie neutronów — porównuje moderatory przez długość spowalniania, dyfuzję i nieucieczkę neutronów.
k_eff — pokazuje, jak geometria, moderator i straty neutronów wpływają na krytyczność układu.
NKE — daje kartę nuklidu z danymi jądrowymi, energiami separacji i porównaniem lokalnych baz.
Bilans cyklu paliwowego — łączy energię, burnup, uran naturalny, ogony i SWU w jednym bilansie materiałowym.

Ćwiczenia praktyczne

Ćwiczenie laboratoryjne powinno polegać na zbudowaniu prostego modelu jednowymiarowej dyfuzji neutronów cieplnych w moderatorze. W wariancie podstawowym należy:

przyjąć wartości $D_T$ i $\Sigma_{a,T}$ dla kilku moderatorów,
policzyć długość dyfuzji $L$ dla wody lekkiej, ciężkiej wody i grafitu,
narysować przebieg $\Phi(r) \sim e^{-r/L}$ dla tych materiałów,
porównać, jak zmienia się zasięg użytecznego strumienia neutronów,
wyciągnąć wniosek, jakie rozmiary geometryczne faworyzuje każdy z moderatorów.

Celem ćwiczenia nie jest pełny projekt rdzenia, lecz uchwycenie związku między własnościami materiału a skalą przestrzenną pola neutronowego. W drugiej części należy dodać warunek brzegowy odpowiadający obecności reflektora neutronów i porównać, jak zmienia się efektywna utrata neutronów na granicy.

Drugie ćwiczenie, teoretyczno-przemysłowe, powinno polegać na wyborze moderatora dla prostego reaktora na uran naturalny lub lekko wzbogacony. Należy:

porównać zyski ze skutecznego spowalniania z kosztami pochłaniania neutronów,
ocenić wpływ jakości grafitu reaktorowego i domieszek absorbujących,
wskazać, kiedy opłaca się użyć ciężkiej wody,
opisać rolę reflektora w zmniejszaniu strat powierzchniowych,
powiązać wynik z warunkami osiągnięcia krytyczności w rzeczywistym rdzeniu.

To ćwiczenie ma pokazać, że teoria dyfuzji jest narzędziem projektowym, a nie tylko eleganckim równaniem. Pozwala bezpośrednio przejść od danych materiałowych do decyzji o wymiarach i architekturze układu.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

Przypisy

Bronisław Słowiński, Podstawy fizyczne energetyki jądrowej. Zakres wykorzystania: definicja dyfuzyjnego prądu neutronów, współczynnika dyfuzji, długości dyfuzji i długości migracji oraz wartości tych parametrów dla wody lekkiej, ciężkiej wody i grafitu.
↩
Bronisław Słowiński, Podstawy fizyczne energetyki jądrowej. Zakres wykorzystania: opis cyklu neutronowego w reaktorze, rola moderatora, ucieczki neutronów i wychwytu rezonansowego w bilansie strat.
↩
Materiały jądrowe. Zakres wykorzystania: praktyczne znaczenie ciężkiej wody i grafitu jako moderatorów oraz ich konsekwencje dla pracy reaktora.
↩
Jeremi Gdenarz, Jak się bać broni jądrowej. Zakres wykorzystania: rola reflektora neutronów jako warstwy zawracającej część uciekających neutronów z powrotem do obszaru aktywnego.
↩
Samuel Glasstone i Milton C. Edlund, The Elements of Nuclear Reactor Theory, Van Nostrand 1952; Alvin M. Weinberg i Eugene P. Wigner, The Physical Theory of Neutron Chain Reactors, University of Chicago Press 1958. Zakres wykorzystania: pełne równanie dyfuzji, wiek Fermiego, tabela parametrów moderatorów (H₂O, D₂O, grafit, beryl), heterogeniczność rdzenia i prawdopodobieństwo ucieczki rezonansowej, historia CP-1, reaktory HFR, bazy danych benchmark.
↩