Streszczenie

Albedo neutronów mówi, jaka część strumienia padającego na materiał wraca z powrotem zamiast zniknąć w pochłanianiu albo przejść dalej. To nie jest lustrzane odbicie pojedynczej cząstki, lecz uśredniony efekt wielu rozproszeń, absorpcji i dyfuzji w materiale.1

Pojęcie jest przydatne w metrologii neutronowej, osłonach, reflektorach reaktorowych i interpretacji eksperymentów z moderatorami. Dobrze pokazuje, że materiał może być „dobry neutronowo” nie dlatego, że odbija wszystko, lecz dlatego, że częściej rozprasza niż pochłania i ma korzystną długość dyfuzji.

Rozszerzenie tematu

Odbicie, które nie jest lustrem

Słowo „albedo” kojarzy się z optyką: jasna powierzchnia odbija dużo światła, ciemna mało. Dla neutronów intuicja jest podobna tylko na bardzo ogólnym poziomie. Neutron wpadający do materiału nie odbija się od powierzchni jak promień od lustra. Wnika do środka, zderza się z jądrami, zmienia kierunek i energię, a czasem zostaje pochłonięty. Jeżeli po serii takich zdarzeń wraca przez powierzchnię, traktujemy go jako część strumienia odbitego.

Dziunikowski i Kalita definiują dyfuzyjny współczynnik odbicia neutronów termicznych jako liczbę mówiącą, jaka część strumienia neutronów padających na powierzchnię ośrodka zostaje od tej powierzchni odbita.1 Jeżeli Psi_p oznacza strumień padający, a Psi_o strumień odbity, to:

Psi_o = alpha * Psi_p,

czyli:

alpha = Psi_o / Psi_p.

Współczynnik alpha jest więc liczbą bezwymiarową. alpha = 0 oznacza brak odbicia w rozważanym modelu, a alpha = 1 idealny powrót całego strumienia. Rzeczywiste materiały leżą pomiędzy tymi granicami.

Dlaczego neutron wraca

Neutron wraca nie dlatego, że powierzchnia jest twarda, lecz dlatego, że rozpraszanie w materiale jest losowe. Po wejściu do ośrodka neutron wykonuje coś w rodzaju trójwymiarowego błądzenia przypadkowego: raz leci głębiej, raz skośnie, raz z powrotem ku powierzchni. Jeżeli zanim zostanie pochłonięty, trafi z powrotem do granicy, opuszcza materiał po stronie, z której przyszedł.

Im większe prawdopodobieństwo rozpraszania i im mniejsze prawdopodobieństwo absorpcji, tym większa szansa powrotu. Dlatego dobry reflektor neutronów powinien mieć duży przekrój rozpraszania i mały przekrój pochłaniania. To samo kryterium pojawia się w artykule o spowalnianiu neutronów i miernikach wodoru, tylko tam interesowała nas przede wszystkim utrata energii, a tutaj powrót strumienia przez granicę.

Rozproszenie i absorpcja

W opisie dyfuzyjnym trzeba rozróżnić średnią drogę swobodną na rozproszenie i średnią drogę do absorpcji. Dziunikowski i Kalita zapisują je przez liczbę jąder w jednostce objętości oraz odpowiednie przekroje czynne:

lambda_r = 1 / (n0 * sigma_r),

lambda_a = 1 / (n0 * sigma_a).1

Pierwsza wielkość mówi, jak często neutron zmienia kierunek po rozproszeniu. Druga mówi, jak długa jest średnia zygzakowata droga neutronu do pochłonięcia. W praktyce wygodniejsza bywa długość dyfuzji L, czyli miara średniej odległości w linii prostej od miejsca powstania albo wejścia neutronu do miejsca jego pochłonięcia.1

Jeżeli lambda_a jest mała, neutrony szybko znikają i albedo jest niskie. Jeżeli rozpraszanie jest częste, a absorpcja mała, neutron może długo błądzić i ma większą szansę wrócić. Z tego powodu materiały bogate w silne pochłaniacze, takie jak bor albo kadm w zakresie termicznym, mogą radykalnie obniżać odbicie neutronów termicznych, nawet jeśli geometrycznie wyglądają jak „osłona”.

Długość dyfuzji i grubość materiału

Albedo rośnie z grubością materiału odbijającego tylko do pewnego momentu. Cienka warstwa może nie dać neutronom wielu okazji do rozproszenia i powrotu. Warstwa grubsza działa lepiej, ale po kilku długościach dyfuzji dokładanie kolejnych centymetrów daje coraz mniejszy zysk, bo neutrony, które weszły bardzo głęboko, mają małą szansę wrócić na właściwą stronę przed absorpcją.

To jest praktycznie ważne. Reflektor neutronowy nie staje się proporcjonalnie lepszy z każdą kolejną warstwą materiału. W pewnym zakresie grubości przyrost jest duży, potem następuje nasycenie. Dziunikowski i Kalita omawiają ten efekt w modelu płyty o skończonej grubości, używając gęstości strumienia, współczynnika dyfuzji i długości dyfuzji.1

W dydaktycznym modelu wystarczy funkcja nasycająca:

alpha(d) = alpha_inf * (1 - exp(-d / L_eff)),

gdzie d jest grubością warstwy, alpha_inf albedo dla warstwy bardzo grubej, a L_eff umowną skalą dochodzenia do nasycenia. To nie jest pełna teoria transportu neutronów, ale dobrze pokazuje sens: cienka warstwa działa słabo, warstwa rzędu kilku długości dyfuzji działa prawie jak półprzestrzeń.

Albedo termiczne i szybkie

Nie ma jednego albedo materiału „w ogóle”. Współczynnik zależy od energii neutronów, kąta padania, geometrii i składu materiału. Ten sam materiał może dobrze odbijać neutrony termiczne, a inaczej zachowywać się wobec neutronów szybkich. Neutron szybki może najpierw tracić energię w zderzeniach, a dopiero potem zachowywać się jak neutron termiczny. Jeżeli w trakcie spowalniania zostanie pochłonięty albo ucieknie w inną stronę, wkład do albedo będzie inny.

Dlatego trzeba odróżniać trzy pytania:

  • ile neutronów wraca,
  • z jaką energią wracają,
  • po jakim czasie wracają.

W reaktorze termicznym powrót neutronu po moderacji może być korzystny. W bardzo szybkim procesie, takim jak rozwój reakcji w małym układzie nadkrytycznym, neutron wracający późno i z mniejszą energią może mieć inną wartość niż neutron, który nigdy nie uciekł z rdzenia. Dlatego albedo jest pojęciem potrzebnym, ale niewystarczającym do pełnej oceny reflektora.

Pomiar albedo jako problem metrologiczny

Źródłowe ćwiczenie Dziunikowskiego i Kality pokazuje dwa klasyczne sposoby wyznaczania albedo: przez wskaźnik aktywacyjny oraz przez liczniki neutronów termicznych.1 Dla serwisu ważniejsza od programu ćwiczenia jest logika pomiaru. Trzeba porównać sygnał związany ze strumieniem padającym i odbitym, a następnie skorygować wynik o geometrię, absorpcję w samym wskaźniku, tło i wydajność detektora.

Model aktywacyjny można wyjaśnić bez instrukcji laboratoryjnej. Jeżeli wskaźnik rejestruje neutrony przechodzące przez pewną powierzchnię, a neutrony mogą wracać po kolejnych rozproszeniach, sygnał może zawierać sumę geometryczną:

Psi + alpha Psi + alpha^2 Psi + ... = Psi / (1 - alpha).

Taka suma dobrze pokazuje, dlaczego nawet umiarkowane albedo może istotnie zwiększyć liczbę przejść neutronów przez obszar pomiarowy. W realnym pomiarze ta intuicja wymaga poprawek, ale dydaktycznie jest bardzo użyteczna.

Związek z reflektorami neutronów

Reflektor neutronów jest materiałem, którego zadaniem jest zwiększyć prawdopodobieństwo powrotu neutronów do obszaru, z którego uciekły. W reaktorze oznacza to zmniejszenie strat neutronów przez powierzchnię rdzenia. W fizyce krytyczności oznacza to zmianę efektywnego bilansu neutronów. W metrologii oznacza to, że obiekty wokół stanowiska mogą zmieniać pole neutronowe przez odbicie i moderację.

W artykule o reflektorze i tamprze uranowym albedo pojawia się jako element wyjaśnienia, dlaczego warstwa otaczająca rdzeń może zawracać część neutronów. Tutaj patrzymy na to samo pojęcie od strony metrologicznej: nie projektujemy układu, lecz pytamy, jak materiał zmienia strumień na granicy.

To rozróżnienie jest istotne dla bezpieczeństwa treści. Albedo jako pojęcie dyfuzyjne można omawiać przez bezpieczne modele, syntetyczne dane i porównanie materiałów. Nie trzeba podawać procedur z materiałami rozszczepialnymi ani instrukcji eksperymentów krytyczności.

Albedo w osłonach i pomieszczeniach

W pomiarach neutronowych otoczenie nigdy nie jest całkowicie neutralne. Ściana, podłoga, stół, pojemnik, moderator, woda w materiale i elementy metalowe mogą rozpraszać neutrony z powrotem do detektora. W spektrometrii neutronowej mówi się czasem o room return, czyli wkładzie neutronów odbitych i rozproszonych przez pomieszczenie. To jeden z powodów, dla których funkcje odpowiedzi i kalibracja detektorów neutronowych są trudne.2

Albedo pomaga nazwać ten problem. Jeżeli przyrząd jest wzorcowany w jednej geometrii, a używany w innej, odbicia od otoczenia mogą zmienić sygnał. Dla promieniowania gamma często myślimy o rozproszeniu w ścianach jako o tle. Dla neutronów dodatkowo dochodzi moderacja: neutron odbity może wrócić z inną energią niż miał pierwotnie.

W dobrym opisie pomiaru neutronowego warto zapisać:

  • materiał i odległość dużych powierzchni w pobliżu,
  • obecność moderatorów i pochłaniaczy,
  • wysokość detektora nad podłożem,
  • czy pomiar odbywał się w geometrii otwartej, czy w pomieszczeniu,
  • czy stosowano korektę tła i odbić.

Materiał odbijający a materiał moderujący

Materiał może dobrze spowalniać neutrony, ale niekoniecznie najlepiej odbijać je w użytecznym sensie. Wodór bardzo skutecznie odbiera energię neutronom, ale zwykła woda również pochłania część neutronów termicznych. Ciężka woda ma mniejsze pochłanianie, więc w reaktorach jest wyjątkowo cenna. Grafit ma małe pochłanianie przy wysokiej czystości, ale spowalnia wolniej. Beryl może być dobrym lekkim reflektorem neutronów szybkich, lecz jego wartość zależy od czasu powrotu i widma energii.3

Albedo porządkuje te kompromisy. Wysokie albedo wymaga nie tylko rozpraszania, ale też przeżycia neutronu do momentu powrotu przez powierzchnię. Dlatego silny moderator z dużą absorpcją może mieć gorszy bilans niż materiał, który spowalnia wolniej, ale mniej pochłania.

Prosty model Monte Carlo

Albedo można dobrze pokazać w miniaturowej symulacji losowej. Użytkownik nie potrzebuje rzeczywistych neutronów ani źródeł promieniowania. Wystarczy model jednowymiarowy:

  • neutron startuje przy powierzchni i idzie w głąb materiału,
  • w każdym kroku może zostać pochłonięty z prawdopodobieństwem p_abs,
  • jeżeli nie zostanie pochłonięty, zmienia kierunek losowo,
  • jeżeli wróci przez powierzchnię wejściową, liczymy go jako odbity,
  • jeżeli przejdzie przez tylną granicę cienkiej warstwy, liczymy go jako przepuszczony.

Po zasymulowaniu wielu neutronów albedo jest po prostu:

alpha = N_return / N_in.

Taki model jest bardzo uproszczony, ale dobrze uczy trzech rzeczy. Po pierwsze, albedo jest wielkością statystyczną. Po drugie, absorpcja i grubość silnie zmieniają wynik. Po trzecie, zwiększanie grubości daje nasycenie, bo neutrony bardzo głęboko rozproszone rzadko wracają.

Czego nie wolno wyczytać z samego albedo

Albedo nie mówi wszystkiego o materiale. Dwa materiały mogą mieć podobny ułamek neutronów odbitych, ale różne widmo energii neutronów wracających. Mogą też mieć różny czas opóźnienia powrotu albo inną zależność od kąta padania. Dla detektora dozymetrycznego, reaktora i szybkiego procesu jądrowego te różnice mogą mieć bardzo różne znaczenie.

Nie można też traktować albedo jako stałej tablicowej bez geometrii. Współczynnik zależy od tego, czy materiał jest półprzestrzenią, płytą, kulą, cienką warstwą, mieszaniną, proszkiem czy układem z pustkami. Zależy od energii i kierunkowości strumienia. Zależy od temperatury tylko pośrednio, przez widmo termiczne i własności materiału, ale w precyzyjnych pomiarach również to może mieć znaczenie.

Najkrótsze podsumowanie: albedo neutronów jest praktyczną miarą powrotu strumienia z materiału. Jego wartość wynika z losowego transportu neutronów, a nie z lustrzanego odbicia, dlatego trzeba je interpretować razem z widmem, geometrią, absorpcją i długością dyfuzji.

Teoria dyfuzji neutronów: równanie i długość dyfuzji

Formalna teoria albedo opiera się na teorii dyfuzji neutronów. Strumień neutronów Φ(r) w jednorodnym ośrodku opisuje równanie dyfuzji Ficka:

$$D \nabla^2 \Phi - \Sigma_a \Phi + S = 0$$

gdzie:

  • D = 1/(3Σ_tr) — współczynnik dyfuzji [cm], Σ_tr = makroskopowy przekrój transportu
  • Σ_a = n_0 σ_a — makroskopowy przekrój pochłaniania [cm⁻¹]
  • S — gęstość źródła neutronów [n cm⁻³ s⁻¹]

Dla układu bez źródeł (S=0) rozwiązanie w prostej geometrii (nieskończona płyta) ma postać:

$$\Phi(x) = A e^{-x/L} + B e^{+x/L}$$

gdzie długość dyfuzji L jest kluczowym parametrem:

$$L = \sqrt{D/\Sigma_a} = \sqrt{1/(3\Sigma_{tr}\Sigma_a)}$$

Fizyczny sens L: L² jest proporcjonalne do średniokwadratowej odległości między miejscem powstania a miejscem pochłonięcia neutronu termicznego. Dlatego L jest miarą „zasięgu" neutronu w materiale.

Tabela długości dyfuzji dla typowych materiałów (neutony termiczne, 25,3 meV):

Materiał D [cm] Σ_a [cm⁻¹] L [cm] Albedo α_∞
Woda H₂O 0,160 0,0197 2,85 ~0,82
Ciężka woda D₂O 0,870 0,000280 171 ~0,998
Grafit (czysty) 0,916 0,000318 54,0 ~0,996
Beryl Be 0,502 0,00120 20,4 ~0,990
Żelazo Fe — — 1,01 ~0,67
Ołów Pb — — 1,28 ~0,70

Dane: NIST, Glasstone & Edlund „The Elements of Nuclear Reactor Theory". Albedo dla półprzestrzeni (α_∞) obliczone z modelu dyfuzyjnego.

Długość dyfuzji D₂O = 171 cm (ponad 60× więcej niż H₂O) wyjaśnia wyjątkowość ciężkiej wody jako moderatora — neutrony przemierzają ogromne odległości przed pochłonięciem, a granica ciężkiej wody działa jak doskonałe lustro neutronowe.

Albedo dla neutronów szybkich: teoria transportu

Dla neutronów szybkich (>1 keV) teoria dyfuzji przestaje być dobrym przybliżeniem — neutrony szybkie mają małą długość dyfuzji, ale duże przekroje (n,2n), inelastyczne rozpraszanie i silne anizotropię rozpraszania. Potrzebne jest pełniejsze równanie transportu Boltzmanna.

Albedo dla neutronów szybkich zależy silnie od energii. Ołów jest wyjątkiem — ma małe pochłanianie i duże przekroje (n,2n), więc jest skutecznym reflektorem neutronów szybkich (stąd jego użycie jako tamper w niektórych urządzeniach). Beton, mimo gęstości podobnej do żelaza, jest efektywnym reflektorem neutronów szybkich ze względu na wodór — spowalnia neutrony, a spowolnione mają wyższy przekrój (n,γ) na azotowo-silikonowej matrycy.

Symulacje Monte Carlo (MCNP, OpenMC, Geant4) są standardowym narzędziem do obliczania albedo dla złożonych geometrii i szerokich widm. Wynik jest funkcją α(E_in, θ_in, E_out, θ_out, material) — 4D obiekt, którego 1D projekcja (całkowanie po kątach i energiach) daje skalarne albedo używane w modelowaniu reaktorów.

Albedo w reaktorze jądrowym: wpływ na k_eff

W fizyce reaktorowej albedo reflektora determinuje, jaka część neutronów uciekających z rdzenia wraca i uczestniczy w kolejnych rozszczepieniach. Wpływ na mnożnik efektywny k_eff:

Bez reflektora (rdzeń otoczony próżnią, albedo = 0):

  • Znaczna część neutronów ucieka przez powierzchnię
  • Straty nieliniowe wzrastają dla małych rdzeni (duży stosunek powierzchnia/objętość)

Z reflektorem (albedo α):

  • Neutrony uciekające i wracające redukują wymagane rozmiary krytyczne
  • Grubość reflektora BU (reflector savings): dla D₂O wynosi kilkadziesiąt centymetrów, dla grafitu ~30 cm, dla H₂O ~10 cm
  • Reflektor obniża wymaganą masę krytyczną i zmniejsza koszty paliwa

W reaktorze MARIA (NCBJ Świerk): basen z wodą otaczający rdzeń pełni rolę moderatora, reflektora (albedo ~0,82) i osłony biologicznej. Bez wody bieżącej reaktor automatycznie wyłącza się (bezpieczeństwo pasywne bazujące na albedo wody).

W reaktorach ciężkowodnych (CANDU, PHWR): beczka moderatora D₂O otacza paliwo, a albedo D₂O ≈ 0,998 pozwala używać niezbogaconego paliwa — neutronowy bilans jest tak dobry, że nawet duże straty na pochłanianiu U-238 mogą być tolerowane.

Porównanie materiałów pod kątem albedo: kompromisy projektowe

Przy wyborze reflektora lub moderatora w reaktorze albo przy projektowaniu osłony neutronowej ważny jest nie tylko stosunek σ_scatter/σ_abs, ale cały zestaw własności materiału:

Materiał L [cm] Gęstość [g/cm³] α_∞ (termiczne) Uwagi
D₂O 171 1,10 ~0,998 Najlepszy moderator, drogi
Grafit (czysty) 54 1,60 ~0,995 Dobry, wrażliwy na B-zanieczyszczenie
Beryl Be 20,4 1,85 ~0,990 Dobry, toksyczny i drogi
H₂O 2,85 1,00 ~0,82 Tani, dostępny, dobre chłodzenie
Polietylen CH₂ 2,7 0,95 ~0,85 Dobry dla termicznych, degraduje pod neutronami
Parafina 2,5 0,87 ~0,83 Jak polietylen, tańszy
Żelazo 1,01 7,87 ~0,67 Złe albedo, ale dobra osłona gamma
Ołów 1,28 11,35 ~0,70 Dobre na neutrony szybkie przez reakcję (n,n')

Beryl ma wyjątkową własność: reakcja (n, 2n) — pochłonięcie neutronu termicznego może dać dwa neutrony z jądra Be-9 przy Q = -1,67 MeV. Przy energiach > 1,67 MeV ten kanał zwiększa „efektywne albedo" berylu powyżej jedności (!). To jest jeden z powodów, dla których beryl (i niektóre reakcje w ołowiu) jest stosowany jako mnożnik neutronów w blanketach fuzyjnych.

Bor w graficie jest klasycznym przykładem niszczenia albedo przez zanieczyszczenie: grafit z zawartością 1 ppm wagowych boru (1 µg B / g C) obniża L z 54 cm do ok. 38 cm, a α_∞ z 0,995 do ok. 0,980. Reaktory z moderatorem grafitowym (Windscale, RBMK) wymagały grafitu o czystości <0,5 ppm B — kontrola ta była strategicznie ważna dla powodzenia programu reaktorowego.

Eksperyment Fermiego z „Chicago Pile-1" (2 grudnia 1941, pierwsza kontrolowana reakcja łańcuchowa) wykazał, jak krytyczna jest czystość grafitu. Fermi i Szilard wielokrotnie mierzyli absorpcję neutronów w próbkach grafitu, szukając dostawcy o najniższej zawartości zanieczyszczeń, i ostatecznie wybrali grafit od National Carbon Company o zawartości boru ~0,4 ppm. Bez tego wyboru masa grafitu w stosie byłaby zbyt duża i projekt mógłby się nie udać.

Albedo w osłonach biologicznych od neutronów

Osłona od neutronów jest bardziej złożona niż osłona od promieniowania gamma, bo wymaga zarówno spowolnienia (moderacji) neutronów jak i ich pochłonięcia. Albedo odgrywa tu nieoczekiwaną rolę.

Niezamierzone albedo od ścian: w laboratorium z silnym źródłem neutronów (np. 252Cf, neutrony z cyklotronu, reaktor) neutrony odbijają się od ścian betonowych i mogą wracać do obszaru pracy jako „tło neutronowe" o mieszanym widmie. Ten „room scatter" może znacząco przewyższyć strumień bezpośredni w obszarach z boku i za osłoną.

Żelazo w osłonach: żelazo ma złe albedo termiczne (~0,67), ale jest dobrym ekranem gamma. W osłonach mieszanych (neutrony + gamma) często stosuje się żelazo + warstwę polietylenu lub boru. Żelazo pochłania gamma, polietylen moderuje i pochłania neutrony.

Streaming przez szczeliny: najważniejszy problem projektowy. Neutron może wielokrotnie odbijać się od ścian wąskiego kanału (streaming), przemierzając długą drogę bez pochłaniania. Albedo ściany kanału determinuje efektywność bariery. Dla betonu (albedo ~0,6) neutron przy każdym odbiciu traci 40% prawdopodobieństwa powrotu — kanał długości 5 × HVL_eff daje redukcję analogiczną do ołowiu dla gamma.

Beton o zwiększonej zawartości wodoru: beton baritowy (BaSO₄ + portlandyt) jest stosowany w osłonach reaktorów i akceleratorów ze względu na dobry kompromis: baru zapewnia pochłanianie gamma, a wodór w portlandycie (Ca(OH)₂) zapewnia moderację neutronów. Albedo betonu baritowego dla termicznych: ~0,72.

Osłony reaktorów wodnych: w reaktorach ciśnieniowych (PWR/VVER) rdzeń jest otoczony wodą pod ciśnieniem, a dalej stalowym dzwonem ciśnieniowym i betonową osłoną biologiczną. Woda jako pierwszy reflektor (albedo ~0,82) zwraca część neutronów do rdzenia — to korzystne dla ekonomii reaktora. Zewnętrzna osłona betonowa pochłania neutrony i gamma — tu albedo wody wewnętrznej jest elementem projektowym bilansu neutronowego.

Monte Carlo i symulacje albedo: MCNP, OpenMC, Geant4

Dla złożonych geometrii i szerokich widm albedo neutronów oblicza się numerycznie metodą Monte Carlo. Standardowe kody to:

MCNP6/MCNPX (Los Alamos National Laboratory): standardowy kod w ochronie radiologicznej, reaktorach i fyzyce cząstek. Bazy danych przekrojów czynnych: ENDF/B-VIII.0, JEFF-3.3 (europejska). Użycie: definiuje się geometrię (materiał, grubości, granice), importuje widmo neutronów padających i liczy transmisję i albedo w trybie F2 tally (strumień przez powierzchnię).

OpenMC (MIT): otwartoźródłowy kod Monte Carlo z podobną funkcjonalnością co MCNP. Coraz szerzej stosowany w akademii (dostępny bezpłatnie). Polska: grupy na AGH i Politechnice Warszawskiej używają go w badaniach reaktorowych.

Geant4 (CERN): stosowany głównie w fizyce cząstek i diagnostyce medycznej. Ma moduł G4HadronPhysics dla neutronów do 20 MeV.

Typowa procedura obliczenia albedo w MCNP:

  1. Definicja geometrii: nieskończona płyta z materiałem badanym, próżnia po obu stronach
  2. Źródło neutronów: punkt na powierzchni z kierunkiem prostopadłym do płyty, widmo monoenergetyczne lub z pliku
  3. Licznik F2 tally na powierzchni wejściowej (mierzy strumień powrotny), normalizowany do strumienia padającego
  4. Wynik: α = F2_return / F2_incident

Zalety MC nad teorią dyfuzji:

  • Pełna geometria 3D, niejednorodne materiały, nieregularne kształty
  • Pełne widmo energii, anizotropowe rozpraszanie
  • Promieniowanie wtórne (gamma z (n,γ), fotony charakterystyczne)
  • Brak aproksymacji dyfuzyjnej — walidowany eksperymentalnie

Wady: czasochłonne obliczenia dla złożonych geometrii (godziny na serwerze klastrowym) i konieczność dobrej znajomości kodu.

Weryfikacja Monte Carlo z eksperymentem jest obowiązkowa dla zastosowań inżynierskich. Dla albedo neutronów termicznych od płyt grafitowych i H₂O obliczenia MCNP zgadzają się z pomiarami aktywacyjnymi w granicach 1–3% (Maerker & Muckenthaler, ORNL, 1965). Dla neutronów szybkich od żelaza błędy mogą być 5–15% ze względu na niepewność przekrojów czynnych w bazie danych. Właśnie dlatego prowadzone są programy benchmarking, w których kody MC walidowane są na zestandaryzowanych eksperymentach (baza SINBAD — Shielding Integral Benchmark Archive and Database, NEA).

Dozymetry albedo: pomiar dawki od neutronów

Zastosowanie albedo w dozymetrii osobistej jest eleganckie i praktyczne. Dozymetr albedo (albedo dosimeter) składa się z detektora neutronowego umieszczonego na ciele człowieka lub na plecach. Ciało człowieka (głównie woda i tkanki miękkie) działa jako moderator i reflektor — neutrony termiczne „odbijają się" od ciała i powracają do detektora od strony tylnej.

Typowa konfiguracja dozymetru albedo:

  • Detektor TLD (LiF:Mg,Ti, TLD-600 bogatszy w ⁶Li lub TLD-700 uboższy)
  • Para detektorów: jeden od strony napromieniowania (mierzy bezpośrednie neutrony + albedo), drugi od strony ciała (mierzy głównie albedo)
  • Okienko kadmowe (Cd, 0,5 mm) pochłania termiczne neutrony bezpośrednie, pozostawiając fast neutrons i albedo

Zasada działania: przy napromieniowaniu szybkimi neutronami od przodu ciało moderuje je, termiczne albedo powraca do tylnego detektora. Stosunek sygnałów detektor-tylny/detektor-przedni jest funkcją energii neutronów — pozwala szacować energię i dawkę ekwiwalentną H*(10).

Norma ISO 9397:1996 określa kalibrację dozymetrów albedo. W Polsce dozymetria albedo jest stosowana w NCBJ, IFJ PAN i kilku szpitalach z cyklotronami (produkcja izotopów medycznych). Pracownicy w polach neutronowych (reaktor MARIA, cyklotrony C16 i C18 w Krakowie) stosują dozymetry albedo jako obowiązkowe elementy nadzoru indywidualnego.

Wyznaczanie albedo metodą aktywacyjną: szczegóły

Metoda aktywacyjna jest klasycznym podejściem do wyznaczania albedo termicznego. Schemat eksperymentu:

  1. Wskaźnik aktywacyjny (np. drut Au, In lub folia Dy) umieszczony między źródłem neutronów a materiałem badanym
  2. Aktywacja A_padające wskaźnika bez materiału reflektora (tylko strumień padający)
  3. Aktywacja A_total wskaźnika z materiałem reflektora (strumień padający + odbity)
  4. Albedo: α = (A_total - A_padające) / A_padające = ΔA/A_padające

Problemy praktyczne:

  • Wskaźnik sam pochłania neutrony i zmienia pole. Korekta: cienki wskaźnik (<1 mg/cm² Aur)
  • Granica geometryczna: jak daleko od powierzchni mierzyć? Konwencja: bezpośrednio przy granicy
  • Rozproszenie od ścian pomieszczenia: pomiar w warunkach o-room (w dużym helu lub wodzie)
  • Aktywacja folii zawiera wkład zarówno od neutronów termicznych jak i rezonansowych. Korekty przez pomiar z filtrem Cd (szybkie neutrony) i bez (termiczne + epi)

Wynik typowy: α(woda, termiczne) = 0,78–0,84 w zależności od geometrii; α(beton, termiczne) ≈ 0,60–0,75; α(polietylen) ≈ 0,85–0,92.

Albedo w polu mieszanym: interpretacja i pułapki pomiarowe

Wiele realnych środowisk pracy z neutronami to pola mieszane: neutrony + gamma z różnych źródeł, o różnych energiach, przychodzące z różnych kierunków. Interpretacja albedo w takim polu jest trudna.

Problem izotropowego pola: teoria dyfuzji i prosta definicja α = Φ_return/Φ_incident zakłada jednoznaczny kierunek padania. W pomieszczeniu, gdzie neutrony przychodzą ze wszystkich stron, „strumień padający" i „odbity" są nieodróżnialne. Detektory kierunkowe (komory szczelinowe, kształtki „shadow cone") są potrzebne do wydzielenia składowych.

Shadow cone method: stosuje się pojemnik (shadow cone) pochłaniający neutrony z określonego kierunku. Pomiar z pojemnikiem i bez niego pozwala wydzielić składową bezpośrednią (strumień padający) i składową scattered (albedo + rozproszenie od ścian). Metoda stosowana przy kalibracji przenośnych monitorów neutronowych.

Energie wracające vs padające: neutron wpadający jako szybki może wrócić jako termiczny po moderacji w grubej ścianie betonowej. Detektor termicznych neutronów przy ścianie mierzy inny albedo niż detektor szerokopasmowy. Dlatego albedo jest własnością układu detektor–materiał–geometria, a nie samego materiału.

Kalibracja detektorów w obecności odbić: przy kalibracji monitorów neutronowych (kule Bonnera, FNT) w przestrzeni laboratoryjnej odbicia od ścian zmieniają widmo przy detektorze. Norma ISO 8529-2 zaleca kalibracje w geometrii „na zewnątrz" lub w dużych pomieszczeniach. Dla neutronów termicznych minimalna odległość od ściany betonowej: >3 m dla pomiarów o niepewności <5% (oszacowanie z NCRP 112).

Kolekcja energii: przy pomiarze dawki ekwiwalentnej H*(10) od neutronów pole mieszane wymaga znajomości widma energii. Albedo zmienia widmo przez preferencyjny powrót składowych termicznych i spowolnionych. Detektor mierzący tylko neutrony termiczne w pobliżu betonowej ściany może zawyżyć mierzoną H*(10) wskutek albedo ściany — podczas gdy neutrony termiczne mają niższy współczynnik jakości Q niż szybkie.

Trzy przykłady obliczeniowe

Przykład 1. Albedo z modelu dyfuzyjnego dla nieskończonej półprzestrzeni

W modelu dyfuzyjnym Milne'a albedo dla wiązki padającej prostopadle na nieskończoną półprzestrzeń wynosi:

α = (1 - 2D/L) / (1 + 2D/L)

gdzie stosujemy warunek brzegowy znikania strumienia w ekstrapolowanej granicy.

Dla wody (D = 0,160 cm, L = 2,85 cm):
2D/L = 2 × 0,160 / 2,85 = 0,112
α = (1 - 0,112) / (1 + 0,112) = 0,888/1,112 = **0,799**

Dla grafitu (D = 0,916 cm, L = 54 cm):
2D/L = 2 × 0,916 / 54 = 0,0339
α = (1 - 0,0339) / (1 + 0,0339) = 0,966/1,034 = **0,935**

Dla D₂O (D = 0,870 cm, L = 171 cm):
2D/L = 2 × 0,870 / 171 = 0,0102
α = (1 - 0,0102) / (1 + 0,0102) = 0,990/1,010 = **0,980**

Wyniki: α(H₂O) ≈ 0,80, α(C) ≈ 0,94, α(D₂O) ≈ 0,98 — zgodne z porządkiem wartości tabel nuklearnych.

Przykład 2. Efektywne zwiększenie strumienia przez wielokrotne odbicia

W układzie: źródło emituje 10⁶ n/s, detektor otoczony od tyłu reflektorem o albedo α = 0,85. Suma geometryczna odbitych neutronów:

Strumień pierwotny + odbity + ponownie odbity + ...:
Φ_eff = Φ₀ × (1 + α + α² + ...) = Φ₀ / (1 - α) = 10⁶ / (1 - 0,85) = 10⁶ / 0,15 = **6,67 × 10⁶ n/s**

Reflektor o albedo 0,85 zwiększa efektywny strumień przy detektorze 6,67× ponad strumień bez reflektora! To wyjaśnia, dlaczego reflektory w reaktorach mogą dramatycznie zmniejszać wymaganą masę paliwa.

Przy albedo 0,95:
Φ_eff = 10⁶ / 0,05 = 20 × 10⁶ — 20-krotny wzrost!

Przykład 3. Obliczenie grubości reflektora dla 90% albedo

Używając modelu α(d) = α_∞ × (1 - exp(-d/L)):

Dla wody, α_∞ = 0,80, L_eff ≈ 2L = 5,7 cm (efektywna skala w modelu z podwójną długością dyfuzji):

Chcemy α(d) = 0,90 × α_∞ = 0,72:
0,72 = 0,80 × (1 - exp(-d/5,7))
1 - exp(-d/5,7) = 0,72/0,80 = 0,90
exp(-d/5,7) = 0,10
d/5,7 = -ln(0,10) = 2,30
d = 2,30 × 5,7 = **13,1 cm** wody

Innymi słowy: warstwa wody 13 cm daje 90% maksymalnego albedo. Kolejne 13 cm dodaje tylko 9% (bo 0,90 × 0,80 = 0,72 → do 0,79). Efekt nasycenia jest oczywisty.

Pytania otwarte dla badaczy i studentów

  • Model dyfuzyjny zakłada, że neutron traci całą pamięć o kierunku po każdym rozproszeniu (izotropowe rozpraszanie). Jak zmodyfikuje się obliczone albedo, gdy uwzględnimy anizotropię rozpraszania (neutron po kolizji częściej leci „do przodu")? Jak parametr asinmetrii rozpraszania ⟨cos θ⟩ wchodzi do D i L?
  • Dozymetr albedo działa poprawnie, gdy ciało ludzkie jest dobrym moderatorem dla neutronów szybkich. Jaka jest optymalna energia neutronów, przy której metoda albedo jest najbardziej czuła? Jak zmienia się wydajność dozymetru albedo dla neutronów o energii 1 keV, 100 keV i 10 MeV z punktu widzenia moderacji w ciele i aktywacji TLD?
  • Reflektor reaktora CANDU jest wykonany z D₂O, a RBMK używa grafitu jako moderatora (i częściowo reflektora). Porównaj albedo tych materiałów i pokaż, jak różnica w albedo tłumaczy, dlaczego CANDU może pracować na niezbogaconym uranie, a RBMK wymagał niskozboganego.
  • W układzie z wieloma reflektorami (rdzeń + reflektor + osłona) albedo na każdej granicy może być inne. Jak formalnie przelicza się strumień przez serię warstw, każda z własnym albedo? Czy jest to analogia do szeregowego układu oporników?
  • Monte Carlo i teoria dyfuzji mogą dawać różne albedo dla tej samej geometrii. Jakie są fizyczne przyczyny rozbieżności i kiedy teoria dyfuzji jest wystarczającym przybliżeniem? Podaj przykłady sytuacji, w których błąd teorii dyfuzji przekracza 10%.
  • Albedo neutronów w berylu (α_∞ ≈ 0,99) jest bliskie jednostce — skąd ta wartość prawie idealna i dlaczego beryl jest tak wartościowy w reaktorach fuzyjnych (jako mnożnik neutronów w blanketach)? Jaka jest rola reakcji ⁹Be + n → 2⁴He + 2n i jak wpływa na efektywne albedo berylowego reflektora?
  • Jaka jest interpretacja fizyczna współczynnika dyfuzji D = λ_tr/3 = 1/(3Σ_tr) i dlaczego stosunek D/Σ_a daje kwadrat długości dyfuzji? Przeanalizuj wymiary: [D] = cm, [Σ_a] = cm⁻¹, [L²] = cm² — pokaż że wzór L = √(D/Σ_a) jest dimensionally consistent.
  • W pomiarach albedo metodą aktywacyjną używa się fotek Au, In lub Dy. Dlaczego te materiały, a nie np. Na lub Fe? Co decyduje o wyborze wskaźnika aktywacyjnego dla pomiaru strumienia neutronów termicznych w obecności neutronów epi- i szybkich?
  • Promieniowanie (n,γ) zachodzące w materiale reflektora produkuje fotony gamma, które są osobnym wkładem do dozy za reflektorem. Jak albedo neutronowe powiązane jest z wtórnym polem gamma i jak uwzględnia się ten efekt przy projektowaniu osłon reaktorowych z grafitu lub D₂O?

Podsumowanie dydaktyczne

  1. Albedo nie jest lustrzanym odbiciem: neutron wraca z materiału przez serię losowych rozpraszań i jest zliczany jako odbity, jeśli wróci przez powierzchnię wejściową — to efekt statystyczny, nie mechanistyczny.

  2. Długość dyfuzji L jest kluczowym parametrem: L = √(D/Σ_a) — większa L oznacza, że neutron może dłużej błądzić przed pochłonięciem, co zwiększa szansę powrotu i daje wyższe albedo.

  3. Albedo nasyca się z grubością materiału: w modelu α(d) = α_∞(1 - exp(-d/L)) po kilku długościach dyfuzji dodawanie materiału przestaje znacząco poprawiać odbicie — istnieje grubość efektywna.

  4. Tabela albedo jest nieliniowo wrażliwa na Σ_a: przejście od H₂O (α ≈ 0,80) do D₂O (α ≈ 0,98) wynika ze zmniejszenia Σ_a o ponad 70×; to nie jest proporcjonalny efekt.

  5. W reaktorze albedo reflektora zmniejsza wymaganą masę krytyczną: efektywne albedo 0,85 może zwiększyć efektywny strumień 6,7× przez wielokrotne odbicia — wyjaśnia to, dlaczego reflektory D₂O i grafitowe pozwalają na ekonomiczne reaktory na naturalne paliwo.

  6. Dozymetry albedo mierzą dawkę od neutronów pośrednio: ciało ludzkie jako moderator + TLD jako detektor to elegancki system detekcji, działający poprawnie w szerokim zakresie energii neutronów (1 keV – 10 MeV).

  7. Pomiar albedo metodą aktywacyjną wymaga korekcji tła i geometrii: aktywacja wskaźnika zależy od sumy strumienia bezpośredniego i odbitego — poprawne wyznaczenie albedo wymaga ostrożnego odjęcia wkładu strumienia bezpośredniego.

  8. Room return jest praktycznym problemem metrologii neutronowej: odbicia od ścian, podłogi i sprzętu laboratoryjnego zmieniają widmo i strumień neutronów, co powoduje błędy przy kalibracji detektorów w małych pomieszczeniach — dlatego kalibracje wykonuje się w dużych halach lub na otwartym powietrzu.

Kontekst fizyki jądrowej: albedo neutronów jest pojęciem, które elegancko łączy trzy obszary — fizykę reaktorów (projekcja reflektora, masa krytyczna), metrologie neutronową (kalibracja detektorów, korekcja tła) i ochronę radiologiczną (projektowanie osłon, dozymetry personalne). Student rozumiejący albedo rozumie, dlaczego D₂O jest tak cenne w reaktorach CANDU, dlaczego grafit RBMK musiał być ekstremalnie czysty, dlaczego kalibracje neutronowe wykonuje się w dużych przestrzeniach, i dlaczego dozymetr na plecach może wskazywać inne wartości niż dozymetr z przodu. To nie jest sucha teoria transportu — to inżynierski i metrologiczny fundament całej fizyki neutronowej.

Dodatkowe materiały multimedialne

Warto przygotować wizualizację random walk: neutrony wchodzą do warstwy, losowo się rozpraszają, część wraca, część zostaje pochłonięta, część przechodzi dalej. Suwaki p_abs, grubość warstwy i średnia droga kroku powinny pokazywać zmianę albedo.

Druga miniaplikacja może porównywać funkcję alpha(d) = alpha_inf * (1 - exp(-d / L_eff)) dla kilku materiałów modelowych. Użytkownik zobaczy, że po kilku długościach dyfuzji przyrost odbicia staje się mały.

Gotowe uzupełnienia interaktywne: kalkulator albedo neutronów, wizualizacja dyfuzji i albedo neutronów oraz model 3D sfer Bonnera.

Ćwiczenia praktyczne

Pierwsze ćwiczenie: dla danych syntetycznych Psi_p = 10000 i Psi_o = 4200 obliczyć alpha. Następnie policzyć strumień odbity dla trzech innych wartości strumienia padającego.

Drugie ćwiczenie: użyć modelu alpha(d) = alpha_inf * (1 - exp(-d / L_eff)). Dla alpha_inf = 0,8 i L_eff = 4 cm obliczyć albedo dla d = 1, 2, 4, 8, 12, 20 cm i narysować wykres nasycenia.

Trzecie ćwiczenie: zasymulować jednowymiarowy random walk dla 10000 neutronów. Porównać trzy przypadki: mała absorpcja, duża absorpcja i cienka warstwa. Wynik zapisać jako N_return, N_abs, N_trans oraz alpha.

Czwarte ćwiczenie: porównać dwa materiały modelowe. Materiał A silnie rozprasza, ale też silnie pochłania; materiał B rozprasza słabiej, ale prawie nie pochłania. Student ma wskazać, który może mieć większe albedo dla grubej warstwy i dlaczego.

Piąte ćwiczenie: opisać, jak albedo otoczenia mogłoby zniekształcić pomiar detektorem neutronów w małym pomieszczeniu. Należy wskazać co najmniej trzy elementy otoczenia, które mogą działać jako moderator, reflektor albo pochłaniacz.

Szóste ćwiczenie: zaprojektować kartę wyników dla pomiaru albedo na danych syntetycznych. Karta powinna zawierać definicję alpha, dane wejściowe, korektę tła, niepewność statystyczną, opis geometrii modelowej i komentarz, czy wynik dotyczy neutronów termicznych, czy ogólnego pola neutronowego.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły