Przejdź do narzędzia

ChainSolver rozwiązuje układy równań różniczkowych zwyczajnych opisujące ewolucję stężeń izotopów w czasie — zarówno podczas napromieniowania w reaktorze, jak i w trakcie spontanicznego chłodzenia. Obliczenia uwzględniają: rozpad radioaktywny (α, β⁻, EC, IT, SF), wychwyt neutronów (n,γ), reakcje rozszczepienia (n,f) oraz opcjonalnie depresję strumienia termicznego i rezonansowe samoochronienie.

Wynikiem jest pełna historia stężeń wszystkich nuklidów łańcucha transmutacji — od materiału wyjściowego do produktów końcowych.

Typowe zastosowania: obliczanie czystości radiochemicznej produkowanych izotopów medycznych i przemysłowych, planowanie harmonogramów napromieniowania, prognozowanie aktywności tarcz po wyjęciu z reaktora, analiza zanieczyszczeń reaktorowych (Xe-135, Sm-149, Pu-240) oraz ocena bilansu materiałów paliwowych w pełnym cyklu paliwowym.1

Wersja webowa a pełny program

Wersja webowa pozwala przeprowadzić kompletną symulację bezpośrednio w przeglądarce — bez instalacji, konfiguracji środowiska Delphi ani kompilowania bibliotek Fortran. Dzięki temu można szybko sprawdzić, czy interesujący nuklid pojawia się w łańcuchu, ocenić rząd wielkości stężeń, porównać kilka harmonogramów napromieniowania i zobaczyć, które parametry mają największy wpływ.

Gdy dzięki wersji webowej znajdziemy interesujący przypadek i zależy nam na najwyższej precyzji (np. do przygotowania wniosku na produkcję izotopu lub do dokumentacji technicznej), warto sięgnąć po pełny program ORIP_XXI ChainSolver — korzysta on z zaawansowanych solverów Fortran (VODE, LSODA, RADAU, MEBDF) z adaptacyjną kontrolą błędu i obsługuje pliki danych OOB z pełnymi przekrojami.1

Audyt modelu i czytanie wyniku

Strona kalkulatora zawiera sekcję Audyt modelu, która ma ułatwić rozróżnienie trzech poziomów obliczenia:

  • równanie transmutacji — czyli macierzowy bilans rozpadów, wychwytów i rozszczepień liczony w czasie;
  • jakość danych jądrowych — okresy półtrwania, przekroje czynne, branching ratio i wydajności rozszczepienia odczytywane z lokalnych baz;
  • warunki fizyczne próbki — widmo neutronów, geometria, samoochrona, historia strumienia, temperatura i skład materiału.

Wynik ChainSolvera należy więc czytać jako ilościową symulację dla zadanych założeń, a nie jako automatyczną odpowiedź na pytanie, jaki będzie skład dowolnej rzeczywistej próbki. Ten sam nuklid startowy może dać inny inwentarz końcowy przy innym widmie neutronów, innym strumieniu, dłuższym chłodzeniu albo przy silnej samoochronie rezonansowej. Audyt modelu pokazuje, które z tych elementów są liczone bezpośrednio, które są ujęte uproszczonym współczynnikiem, a które trzeba świadomie przenieść do osobnego scenariusza wejściowego.

W praktyce edukacyjnej najważniejsza jest kolejność pracy: ChainFinder pokazuje, jakie drogi są formalnie dostępne, ChainSolver sprawdza, które z nich mają istotny wkład liczbowy, a NKE pozwala wrócić do danych pojedynczego nuklidu.

Formularz zawiera presety scenariuszy, które ustawiają typowe punkty startowe i zakresy symulacji: rozpad Co-60, chłodzenie Cs-137, aktywację stali, ciąg U-238 → Pu-239, produkcję Po-210 oraz rozpad trytu. Presety nie są „prawdą technologiczną”; są skrótem dydaktycznym pozwalającym szybciej przejść do pytania, jak czytać wynik.

Wynik zawiera tabelę dominujących wkładów. Osobno wskazywane są nuklidy dominujące masą, aktywnością, ciepłem rozpadu i prostym wskaźnikiem gamma. To ważna korekta interpretacyjna: nuklid o największej masie często nie jest nuklidem o największej aktywności, a nuklid dominujący w cieple nie musi dominować w widmie gamma.

Układ równań różniczkowych (ODE)

Równanie bilansu koncentracji Ni nuklidu i w łańcuchu transmutacji:

dNi/dt = Σj [λj · Nj · bji] — napływ przez rozpad j→i
    + Σj [φ · σc,j · Nj · Yji] — napływ przez wychwyt/rozszczepienie j→i
    − (λi + φ·σc,i + φ·σf,i) · Ni — odpływ: własny rozpad i reakcje

Gdzie:

  • Ni — stężenie nuklidu i [at/cm³]
  • λi = ln(2)/T½ — stała rozpadu [s⁻¹]
  • bji — branching ratio (prawdopodobieństwo przejścia j→i)
  • φ — strumień neutronów termicznych [n/(cm²·s)]
  • σc,i — przekrój czynny wychwytu neutronów [cm²]
  • σf,i — przekrój czynny rozszczepienia [cm²]
  • Yji — wydajność rozszczepienia (produkcja i przez rozszczepienie j)

Układ jest liniowy dla stałego φ w danym etapie — macierz transmutacji A pozostaje stała, co pozwala na bezpośrednie obliczenie eksponentu macierzowego e^(At)·N(0). Z tego korzystają metody Cranka-Nicolson i Padé implementowane w wersji webowej.

Metody numeryczne — wersja webowa

Crank-Nicolson (domyślna, zalecana)

Implementuje aproksymację Padé [1/1] eksponentu macierzowego:

N(t+h) = (I − h/2·A)⁻¹ · (I + h/2·A) · N(t)

Metoda jest bezwarunkowo stabilna — krok h może być dowolnie duży bez ryzyka rozbieżności numerycznej. Jest to szczególnie ważne przy jednoczesnej obecności nuklidów o T½ od nanosekund (izomery jądrowe) do miliardów lat (aktynidy) w tym samym łańcuchu. Rozkład LU macierzy (I − h/2·A) wykonywany jest raz na etap napromieniowania, co czyni każdy krok obliczeniowo tani (samo podstawianie wsteczne). Rząd dokładności: 2 (błąd ∝ h²).

Padé [2/2] ze skalowaniem i potęgowaniem

Aproksymacja wyższego rzędu z algorytmem scaling & squaring:

exp(hA) ≈ (I − hA/2 + h²A²/12)⁻¹ · (I + hA/2 + h²A²/12)

1. Skaluj: B = hA / 2^s tak, żeby ‖B‖∞ ≤ 0,5
2. Oblicz Padé[2/2](B)
3. Potęguj: exp(hA) = (Padé[2/2](B))^(2^s)

Rząd dokładności: 4. Wolniejsza od CN dla dużych układów (wymaga obliczenia B²), ale dokładniejsza dla małej liczby kroków.

RK4 — klasyczna metoda Runge-Kutty

Dostępna jako opcja dla scenariuszy niesztywnych. Metoda jawna rzędu 4 z stałym krokiem. Dla układów z dużą rozpiętością stałych czasowych (sztywnych) wymaga bardzo małego kroku — niezalecana do symulacji obejmujących jednocześnie izomery i długożyjące aktynidy.

RKF45 — adaptacyjna metoda Runge-Kutty-Fehlberga

Rozszerza RK4 o automatyczną kontrolę kroku. Używa tablicy Butchera Fehlberga z 6 etapami i oblicza jednocześnie dwa rozwiązania — rzędu 4 i rzędu 5. Różnica między nimi daje szacunek błędu lokalnego:

  • błąd ≤ tolerancja → akceptuje krok, ewentualnie zwiększa h (wzór:   hnowy = h · 0,9 · (1/błąd)^0.2)
  • błąd > tolerancja → odrzuca krok, zmniejsza h, powtarza obliczenie

Nadal metoda jawna — dla silnie sztywnych układów (jednoczesna obecność izomerów i aktynidów) wymaga małych kroków i będzie wolna. Właściwa dla scenariuszy bez dużej rozpiętości stałych czasowych, gdy ważna jest dokładność przy minimalnej liczbie ewaluacji funkcji.

Solvery oryginalne (pełny program Fortran)

Pełny program ChainSolver korzysta z czterech zaawansowanych solverów wywoływanych jako biblioteki DLL:

LSODA (ODEPACK, LLNL) — automatycznie przełącza się między metodą Adamsa (niesztywna, rzędy 1–12) a BDF (sztywna, rzędy 1–5) na podstawie analizy wartości własnych Jacobiana. Obsługuje rzadkie macierze (sparse Jacobian) i adaptacyjny krok z lokalną kontrolą błędu.

VODE (Brown, Hindmarsh, Byrne) — solver BDF/Adams zmienno-współczynnikowy w postaci Nordsieck. Oferuje bardziej niezawodne sterowanie krokiem niż LSODE i jest szczególnie polecany do układów transmutacji z szeroką rozpiętością stałych czasowych.

RADAU (Hairer & Wanner, Universität Genf) — niejawna metoda Runge-Kutty Radau IIA rzędu 5/9/13 ze zmiennym rzędem. Obsługuje układy DAE (różniczkowo-algebraiczne, M·y'=f). Wyjątkowo dokładna dla problemów o silnej sztywności.

MEBDF (Cash, CWI) — zmodyfikowane rozszerzone BDF. Stabilność A(α) przy wyższych rzędach niż klasyczne BDF.

Wieloetapowe napromieniowanie

Wersja webowa obsługuje dowolną liczbę etapów napromieniowania, każdy z własnym strumieniem φ:

  • φ = 0 → czysty rozpad (chłodzenie po wyjęciu z reaktora)
  • φ > 0 → aktywne napromieniowanie
  • Etapy mogą się przeplatać: napromieniowanie → chłodzenie → ponowne napromieniowanie

Stężenia z końca etapu i stają się warunkami początkowymi etapu i+1. Nowa macierz transmutacji A jest budowana dla każdego etapu (różne φ mogą zmienić względną wagę wychwytu neutronów i czystego rozpadu).

Przykładowy harmonogram 3-etapowy dla Pu-239:

  1. 2 lata napromieniowania, φ = 10¹⁴ n/(cm²·s)
  2. 1 rok chłodzenia, φ = 0
  3. 1 rok ponownego napromieniowania, φ = 5·10¹³ n/(cm²·s)

Korekcje fizyczne

Depresja strumienia termicznego

W gęstych próbkach silnych pochłaniaczy (Gd, Cd, B) strumień neutronów termicznych maleje wzdłuż drogi przez próbkę. Czynnik depresji D obliczany metodą cięciwy:

D = (1 − exp(−Σa · l̄)) / (Σa · l̄)

Gdzie l̄ to średnia cięciwa geometryczna próbki (dla kuli r: l̄ = 4r/3), a Σa = Σi Ni·σa,i to makroskopowy przekrój pochłaniania. Efektywny strumień φ' = D·φ.

Samoochrona rezonansowa (SSK)

Nuklidy z dużymi rezonansami (U-238, Xe-135, Sm-149, Gd-157) pochłaniają neutrony o energiach rezonansowych tak skutecznie, że wewnątrz próbki ich „widzialny" strumień jest zmniejszony. Efekt opisuje czynnik samoochrony SSK:

SSK(N) = σeff(N) / σ0 ≤ 1

SSK maleje wraz z rosnącą koncentracją N — gęstsza próbka = większa samoochrona.

Implementacja korzysta z metody Goldsteina oraz tabel prawdopodobieństwa z biblioteki danych rezonansowych ENDF/B-VIII.1 (234 nuklidy, 33410 rezonansów).2

Odkrywanie łańcucha transmutacji

Przed każdą symulacją ChainFinder (algorytm BFS) automatycznie buduje sieć nuklidów osiągalnych z zadanego punktu startowego przez:

  • Rozpad radioaktywny (α, β⁻, EC, IT, SF) — przy dowolnym φ
  • Wychwyt neutronów (n,γ) → A+1 — tylko gdy φ > 0
  • Rozszczepienie (n,f) → fragmenty — tylko gdy φ > 0 i σ_f > 0

Głębokość przeszukiwania (1–8) kontroluje rozmiar układu ODE — typowo 20–50 nuklidów. Tryby rozpadu dla nuklidów nieobecnych w bazie σa wyznaczane są heurystycznie z reguł fizyki jądrowej (linia stabilności, reguła ekscesu neutronów).

Przykład: Pu-239 po 2 latach w reaktorze

Ewolucja stężeń nuklidów — Pu-239 napromieniowany przez 2 lata przy φ = 10¹⁴ n/(cm²·s)
Ewolucja stężeń nuklidów — Pu-239 napromieniowany przez 2 lata przy φ = 10¹⁴ n/(cm²·s)

Symulacja: 1 g Pu-239, φ = 10¹⁴ n/(cm²·s), 2 lata, metoda Crank-Nicolson, łańcuch 39 nuklidów. Uruchom tę symulację

Nuklid Z A T½ Stężenie końcowe Frakcja masowa
Pu-240 94 240 6,56 ky 7,0442% 85,98%
Pu-241 94 241 14,3 lat 0,6726% 8,24%
Pu-239 94 239 24,1 ky 0,4732% 5,75%
U-236 92 236 23,4 My 0,0013% 0,016%
U-237 92 237 6,75 dni 0,0009% 0,011%

Pu-240 („trucizna") dominuje po intensywnym napromieniowaniu — po 2 latach stanowi ok. 86% frakcji masowej. Potwierdza to, że produkcja plutonu broniowego wymaga krótkich cykli i dużego φ (reaktory Hanford: φ ≈ 3·10¹⁴, cykl ~100 dni, Pu-240 < 1%).

Dane wejściowe i wyjściowe — wersja webowa

Parametry wejściowe:

Parametr Zakres / jednostki Opis
Nuklid startowy Symbol-A (np. U-238, Am-241m1) Punkt startowy sieci transmutacji
Głębokość łańcucha 1–8 generacji Limit BFS przy budowaniu sieci nuklidów
Metoda ODE CN / Padé [2/2] / RK4 Solver układu równań transmutacji
Etapy napromieniowania dowolna liczba; t [s–My], φ [n/cm²/s] Harmonogram: czas trwania + strumień neutronów
Cięciwa geometryczna [cm] Do obliczenia depresji strumienia (opcja)
Koncentracja N₀ [at/cm³] Do obliczenia samoochrony SSK (opcja)
Temperatura [K], domyślnie 300 K Wpływ na poszerzenie Dopplera rezonansów

Wyniki wersji webowej:

  • Wykres — względne stężenia Ni(t)/N0 wszystkich nuklidów (Chart.js, interaktywny)
  • Tabela stanu końcowego — stężenie procentowe i frakcja masowa każdego nuklidu, posortowane malejąco
  • Grafy połączeń (ChainFinder) — sieć przejść reaktywnych, klikalny SVG

Wyniki dostępne w oryginalnej wersji (Fortran)

  • Liczba jąder, masy [g], aktywności α/β/EC/IT/całkowita, aktywność właściwa na gram
  • Udziały izotopów, wydajność masowa na gram materiału wyjściowego
  • Depresja strumienia (jeden czynnik D na każdy interwał czasowy)
  • SSK: wartości dokładne i przybliżone
  • Depozyt energii rozszczepień [MeV] (przy założeniu 200 MeV/rozszczepienie)
  • Ścieżki syntezy dominujących izotopów (analiza pathways)

Instalacja oryginalnego programu ORIP_XXI

Wymagania:

  • Windows (32-bit lub 64-bit z kompatybilnością Win32)
  • Git for Windows
  • Delphi 5 Professional z biblioteką RX 2.75
  • solvery ODE jako DLL (VODE, LSODA, RADAU, MEBDF) — dostarczone z paczką, skompilowane GCC 2.95.2/MinGW.

Pobieranie:

Instalacja krok po kroku:

  1. Umieść plik ORIP_XXI.oob w podkatalogu bin/. Plik OOB zawiera dane nuklearne dla >2800 izotopów.
  2. Otwórz projekt ChainSolver.dpr w Delphi 5 i skompiluj, lub uruchom gotowy bin/ChainSolver.exe.
  3. Upewnij się, że pliki DLL solverów (LSODA.dll, VODE.dll, RADAU.dll, MEBDF.dll) są w tym samym katalogu co EXE.
  4. W zakładce „Conditions" wybierz solver, ustaw strumień neutronów i harmonogram. W zakładce „Chain" wybierz nuklid startowy i głębokość łańcucha. Kliknij „Calculate".

Wersja webowa korzysta z tego samego pliku danych OOB (skonwertowanego do SQLite), co pozwala porównywać wyniki między obiema wersjami.1

Tak powinny wyglądać okna wszystkich trzech programów w oryginalnych wersjach dla Windows:

NKE — okno programu
ChainSolver — okno programu
ChainFinder — okno programu

Streszczenie

ChainSolver rozwiązuje układy równań różniczkowych opisujące ewolucję stężeń izotopów w łańcuchach transmutacji — zarówno podczas napromieniowania neutronami, jak i spontanicznego rozpadu. Narzędzie służy do obliczania czystości izotopów medycznych, prognozowania aktywności tarcz, analizy skażenia paliwa jądrowego i badania zachowania zanieczyszczeń reaktorowych.

Rozszerzenie tematu

ChainSolver jest końcowym ogniwem w trzypoziomowej hierarchii kalkulatorów ORIP_XXI. Podczas gdy NKE dostarcza danych o pojedynczym nuklidzie, a ChainFinder mapuje topologię łańcucha transmutacji, ChainSolver odpowiada na pytanie ilościowe: ile danego nuklidu powstanie po określonym czasie napromieniowania i chłodzenia? Pozornie proste pytanie kryje znaczną złożoność matematyczną i fizyczną — a jej zrozumienie jest kluczem do świadomego korzystania z narzędzia i interpretacji wyników.

Matematyczne podstawy transmutacji — równania Batemana

Historia. W 1910 roku Harry Bateman, angielski matematyk, opublikował analityczne rozwiązanie układu równań różniczkowych opisujących szeregowy łańcuch rozpadów radioaktywnych. Rozwiązanie Batemana zakłada: stały strumień zewnętrzny (lub jego brak), brak sprzężeń zwrotnych, liniową zależność tempa produkcji od stężeń nuklidów "rodziców". Dla prostych łańcuchów szeregowych (A→B→C→D) istnieją zamknięte wzory analityczne — np. aktywność córki jako funkcja czasu.

Ograniczenia rozwiązania analitycznego. Analityczne rozwiązanie Batemana jest możliwe dla prostych łańcuchów liniowych. Dla rzeczywistych sieci transmutacji z rozgałęzieniami, cyklami (nuklid X może rozpaść się do Y, ale Y może wychwycić neutron i dać Z, które z kolei daje X) i setkami nuklidów — konieczne jest podejście numeryczne. Stąd potrzeba solverów numerycznych, takich jak implementowane w ChainSolverze.

Macierzowe sformułowanie. Układ N równań dN_i/dt = Σ A_ij · N_j można zapisać jako:
dN/dt = A·N
gdzie A jest macierzą transmutacji (N×N), a N — wektorem stężeń. Elementy macierzy A kodują: szybkości rozpadów (człony diagonalne < 0), przepływy z nuklidów "rodziców" (człony pozadiagonalne > 0), wychwyty neutronów (człony skalowane przez φ).

Formalne rozwiązanie: N(t) = exp(At) · N(0), gdzie exp(At) to eksponent macierzowy. Obliczenie eksponentu macierzowego efektywnie (bez jawnego obliczania wszystkich wartości własnych) jest sednem różnych metod numerycznych implementowanych w ChainSolverze.

Eksponent macierzowy — dlaczego to trudny problem

Obliczenie exp(At) dla macierzy A transmutacji jądrowej jest numerycznie trudne z kilku powodów:

Sztywność (Stiffness). Układy równań transmutacji są ekstremalnie "sztywne": łańcuchy zawierają zarówno izomery jądrowe (T½ < 1 ns) jak i aktynidy (T½ > 10⁸ lat). Stosunek najkrótszej do najdłuższej stałej czasowej wynosi często 10²⁰–10²⁵. Metody jawne (jak standardowy RK4) dla takich układów wymagają kroków czasowych rzędu nanosekund — nawet gdy interesujące zjawiska zachodzą na skali lat. To czyni je praktycznie bezużytecznymi dla realnych obliczeń.

Skalowanie numeryczne. Elementy macierzy A rozciągają się od λ ≈ 10¹⁰ s⁻¹ (izomery jądrowe) do λ ≈ 10⁻¹⁸ s⁻¹ (Th-232). Bezpośrednie operacje macierzowe mogą tracić precyzję numeryczną (underflow/overflow). Stąd metody skalowania (scaling & squaring w Padé [2/2]) lub użycie arytmetyki o rozszerzonej precyzji.

Nilpotentność. Macierze transmutacji dla izotopów stabilnych (λ=0, σ=0) mają wiersze złożone z samych zer — macierz jest osobliwa. Metody oparte na rozkładzie LU muszą obsłużyć takie przypadki.

Wybór metody. Crank-Nicolson (domyślna w wersji webowej): bezwarunkowo stabilna, rząd 2, dobra równowaga między dokładnością i efektywnością. LSODA i VODE (pełny program): adaptacyjne, automatycznie przełączają się między metodami Adams-BDF zależnie od lokalnej sztywności. RADAU: idealny dla najcięższych problemów. Brak jednej "najlepszej" metody — wybór zależy od charakterystyki konkretnego układu.

Fizyka wychwytu neutronów i jej wpływ na obliczenia

ChainSolver obsługuje nie tylko rozpad radioaktywny, ale i transmutację w polu neutronów. To wymaga poprawnego ujęcia fizyki wychwytu:

Strumień neutronów i przekroje czynne. Parametr φ (strumień neutronów termicznych) determinuje tempo wychwytu. Dla nuklidów z dużymi przekrojami czynnymi (np. Xe-135: σ_a ≈ 2,6×10⁶ barn!), nawet umiarkowany φ może radykalnie zmienić ewolucję łańcucha. Xe-135 "zatruje" reaktor — pochłaniając neutrony tak skutecznie, że łańcuchowa reakcja jądrowa może się zatrzymać. ChainSolver pozwala modelować ten efekt.

Rezonansowe samoochronienie (SSK). W grubych próbkach nuklidy z dużymi rezonansami pochłaniają neutrony z energii rezonansowej tak skutecznie, że w głębi próbki brakuje tych neutronów. Efektywny przekrój czynny jest mniejszy niż przekrój z tabel (σ₀). Współczynnik SSK < 1 koryguje to efekt. Ignorowanie SSK daje zawyżone tempo wychwytu w środku próbki.

Widmo neutronów. Reaktory badawcze (basenowe, jak MARIA) mają widmo termiczno-epitermalne. Reaktory energetyczne — termiczne lub prędkie. Izotopy produkowane w różnych reaktorach mogą mieć różną czystość ze względu na różne wartości σ_a dla różnych energii neutronów. ChainSolver zakłada jedno widmo na etap — dla precyzyjnych obliczeń konieczna jest dobra charakteryzacja widma.

Depresja strumienia. W próbkach silnie pochłaniających (Gd, Cd) strumień termiczny maleje wzdłuż drogi w próbce. Dla cienkich próbek efekt jest pomijalny; dla grubych — kluczowy. Współczynnik depresji D < 1 redukuje efektywny strumień. Ignorowanie depresji daje zawyżone tempo pochłaniania dla grubych próbek.

Zastosowania medyczne — produkcja izotopów terapeutycznych i diagnostycznych

Jednym z kluczowych zastosowań ChainSolvera jest planowanie i weryfikacja produkcji izotopów medycznych:

Technet-99m (⁹⁹ᵐTc). Najczerzej stosowany izotop w medycynie nuklearnej (~80% procedur diagnostycznych na świecie). Produkowany z Mo-99, który z kolei jest produktem rozszczepienia U-235 lub produktem napromieniowania Mo-98 (n,γ). ChainSolver może modelować narastanie ⁹⁹Mo w tarczy Mo-98 w reaktorze i późniejszy wzrost ⁹⁹ᵐTc. Kluczowy parametr: "czas szczytowy" aktywności Tc-99m — moment, w którym wyodrębnienie Tc daje maksymalną aktywność.

Lutet-177 (¹⁷⁷Lu). Izotop terapeutyczny (emiter β⁻ i gamma) do celowanej terapii guzów neuroendokrynnych i prostaty (177Lu-DOTATATE, 177Lu-PSMA). Produkcja dwoma drogami:

  • Bezpośrednia: napromieniowanie ¹⁷⁶Lu (n,γ) → ¹⁷⁷Lu. Daje Lu-177 "carrier-free" (bez stabilnych izotopów Lu — ważne dla znakowania)
  • Pośrednia: napromieniowanie ¹⁷⁶Yb → ¹⁷⁷Yb → ¹⁷⁷Lu (beta-rozpad Yb-177). Wymaga chemicznego rozdzielenia Yb od Lu.
    ChainSolver modeluje obie drogi i ich zanieczyszczenia (Lu-177m, Lu-178) jako funkcję czasu napromieniowania.

Rad-223 (²²³Ra). Emiter alfa do leczenia przerzutów kostnych raka prostaty (lek Xofigo). Produkowany przez rozkład Ac-227 → ²²³Ra → ²¹⁹Rn → ... → ²⁰⁷Pb. ChainSolver modeluje narastanie Ra-223 w generatorze Ac-227 (analogia do generatora Mo/Tc) i optymalizuje harmonogram "dojenienia" generatora.

Aktualny problem braku Mo-99. IEA POLATOM Świerk jest jednym z europejskich dostawców ⁹⁹Mo produkowanego w reaktorze MARIA. Globalne braki Mo-99 w 2009–2010 i 2019 roku pokazały podatność rynku izotopów medycznych na przerwy w reaktorach. Modelowanie alternatywnych harmonogramów napromieniowania i czystości produktów przez ChainSolver jest praktycznie użyteczne dla planowania dostaw.

Zastosowania w fizyce reaktorów

Poza produkcją izotopów, ChainSolver ma zastosowania w analizie reaktorów:

Zatrucie Xenonowe. Po rozruchu lub zmianie mocy reaktora: ¹³⁵I (T½ = 6,6h) rozkłada się do ¹³⁵Xe (T½ = 9,2h, σ_a = 2,6×10⁶ barn). Xe-135 pochłania neutrony, zmniejszając reaktywność. Po awaryjnym wyłączeniu reaktora (SCRAM): produkcja Xe-135 chwilowo wzrasta (bo zanika pochłanianie przez sam reaktor), a dopiero po kilku godzinach zaczyna maleć. To "xenon pit" — reaktor przez kilka godzin nie może być ponownie uruchomiony (reaktywność za niska). Awaria w Czarnobylu (1986) była częściowo spowodowana próbą uruchomienia reaktora właśnie w trakcie xenon pit. ChainSolver modeluje ewolucję Xe-135 jako funkcję czasu po SCRAM.

Zatrucie Samarowe. Sm-149 (stable, σ_a = 41 000 barn) produkowany z Pm-149 (T½ = 53h, produkt rozszczepienia). W trakcie pracy reaktora osiąga stan stacjonarny. Po SCRAM: produkcja Pm ustaje, ale Sm-149 nie znika (jest stabilny, nie rozpada się). Reaktywność stopniowo spada przez akumulację Sm-149 przez kilka tygodni.

Bilans materiałów paliwowych. W typowym reaktorze PWR po jednym cyklu paliwowym (18 miesięcy): paliwo zawiera U-235, U-236, Pu-239, Pu-240, Pu-241, Am-241, Cm-242 i setki innych aktynidów i produktów rozszczepienia. ChainSolver (z odpowiednimi danymi wejściowymi) może modelować tę ewolucję na poziomie pojedynczego elementu paliwowego. Precyzja jest tutaj niższa niż w dedykowanych kodach (ORIGEN2, SERPENT) ze względu na uproszczone widmo neutronów.

Aktywacja materiałów strukturalnych. Stal nierdzewna reaktora po wieloletniej pracy zawiera aktywowane izotopy: Fe-55, Fe-59, Co-58, Co-60, Ni-63. ChainSolver modeluje narastanie tych izotopów podczas pracy i ich zanikanie po wygaszeniu reaktora. Wyniki są podstawą do planowania demontażu i zarządzania odpadami.

Interpretacja wyników — pułapki i najczęstsze błędy

Mylenie aktywności z masą. Nuklid o największej frakcji masowej ≠ nuklid o największej aktywności. Przykład z symulacji Pu-239: Pu-240 dominuje masą (86%), ale jego aktywność właściwa (Ci/g) jest niższa niż krótszożyjących zanieczyszczeń. ChainSolver wyraźnie rozdziela kolumny "stężenie" i "aktywność" — ważne, by używać właściwej przy projektowaniu tarcz.

Problem czystości radiochemicznej. W produkcji izotopów medycznych ważna jest czystość radiochemiczna: stosunek aktywności pożądanego izotopu do całkowitej aktywności próbki. ChainSolver daje tę wartość bezpośrednio w tabeli wyników. Typowy wymóg regulacyjny: czystość > 99% (lub wyższa dla terapeutyków).

Zbyt mała głębokość łańcucha. Wybranie głębokości 1–2 może pominąć ważne nuklidy pojawiające się w dalszych generacjach. Przykład: napromieniowanie Gd-157 przy głębokości 2 może pominąć Tb-158 i Tb-160 produkowane przez wychwyt przez produkty wychwytu. Dla nieznanych układów warto uruchamiać symulacje z rosnącą głębokością i obserwować, kiedy wyniki stabilizują się.

Różne wyniki dla różnych metod. Jeśli CN i Padé [2/2] dają znacznie różne wyniki — to sygnał problemów (zbyt duży krok, źle dobrane parametry). W normalnych warunkach obie metody powinny zgadzać się do 4+ cyfr znaczących.

Ekstrapolacja poza zakres ważności. ChainSolver zakłada stały strumień neutronów w każdym etapie. W rzeczywistym reaktorze strumień zmienia się z wypaleniem paliwa. Dla krótkich napromieniowań (kilka dni) to nie problem; dla symulacji kilkuletnich cykli — konieczna jest dyskretyzacja na wiele etapów z różnymi φ.

ChainSolver w kontekście certyfikacji GMP izotopów medycznych

Produkcja izotopów medycznych dla celów diagnostycznych i terapeutycznych podlega regulacjom GMP (Good Manufacturing Practice):

Europejskie Farmakopeae i wymagania jakościowe. Monografie farmakopei europejskiej (Ph. Eur.) dla izotopów medycznych precyzują: wymagania czystości radionuklidowej (purity), limit zanieczyszczeń chemicznych, kontrolę czasu produkcji, dokumentację kalibracji aktywności. ChainSolver (lub jego Fortran odpowiednik) jest narzędziem do prognozowania składu produktu przed napromieniowaniem (pre-check) i weryfikacji po napromieniowaniu (post-check).

Walidacja oprogramowania. W środowisku GMP: oprogramowanie do obliczeń musi być zwalidowane (validated per FDA 21 CFR Part 11 lub EU Annex 11). Oznacza to: dokumentację metod obliczeniowych (protokół walidacyjny), weryfikację wyników na znanych przypadkach, kontrolę wersji oprogramowania. Wersja webowa kalkulatora nie jest zwalidowanym narzędziem GMP — służy edukacji; do decyzji produkcyjnych używa się pełnego ORIP_XXI z odpowiednią dokumentacją.

Kalibracja systemu danych jądrowych. Wyniki ChainSolvera zależą od jakości danych jądrowych (T½, σ_a, branching ratios). Dla certyfikowanej produkcji izotopów: używane biblioteki danych muszą być udokumentowane (ENDF/B, JEFF) i regularnie aktualizowane gdy ukazują się nowe wersje bibliotek.

8 Otwarte pytania badawcze

  1. Jak dokładnie różne biblioteki danych jądrowych (ENDF/B-VIII, JEFF-3.3, JENDL-5) wpływają na wyniki symulacji produkcji ¹⁷⁷Lu? Która biblioteka lepiej odtwarza wyniki eksperymentalne?

  2. Dla jakich kategorii problemów transmutacyjnych metoda Crank-Nicolson daje nieakceptowalnie duże błędy — i jak zidentyfikować takie przypadki bez a priori wiedzy o układzie?

  3. Jak modelować zmienność strumienia neutronów w trakcie cyklu reaktora dla precyzyjnych obliczeń bilansu paliwa? Ile etapów ChainSolvera jest optymalnych?

  4. Jakie są granice stosowalności modelu jednoenergetycznego w ChainSolverze dla izotopów z dużymi przekrojami czynnymi w obszarze epitermicznym?

  5. Jak dobrać głębokość łańcucha ChainFindera dla nieznanego układu, żeby zagwarantować uwzględnienie wszystkich istotnych nuklidów bez nadmiernego wzrostu rozmiaru układu ODE?

  6. Jaki jest potencjał ChainSolvera do modelowania transmutacji odpadów jądrowych (Minor Actinides — Np, Am, Cm) w reaktorach prędkich lub ADS (Accelerator-Driven Systems)?

  7. Jak zmieniają się wyniki ChainSolvera dla standardowych scenariuszy produkcji Lu-177 przy przejściu z reaktora badawczego (widmo termiczne) na reaktor SMR z widmem epitermicznym?

  8. Czy możliwa jest automatyczna walidacja wyników ChainSolvera przez porównanie z eksperymentalnymi danymi aktywacji z reaktora MARIA w Świerku?

Słownik pojęć kluczowych

Bateman, równania — układ liniowych równań różniczkowych opisujących ewolucję stężeń nuklidów w łańcuchu radioaktywnym; sformułowany przez Harry'ego Batemana w 1910 roku.

BDF (Backward Differentiation Formula) — rodzina niejawnych metod numerycznych do rozwiązywania sztywnych układów ODE; rzędy 1–6; stosowane w solverach LSODA i VODE.

Branching ratio — prawdopodobieństwo danego kanału rozpadowego nuklidu; np. Cs-137 rozpada się 94,6% przez β⁻ do metastabilnego Ba-137m, a 5,4% bezpośrednio do podstawowego Ba-137.

Crank-Nicolson — niejawna metoda różnicowa rzędu 2, bezwarunkowo stabilna; domyślna metoda wersji webowej ChainSolvera.

Czystość radionuklidowa — stosunek aktywności pożądanego izotopu do całkowitej aktywności próbki; krytyczny parametr w produkcji izotopów medycznych.

Depresja strumienia — redukcja strumienia neutronów termicznych wewnątrz grubej próbki silnie pochłaniającej; modelowana współczynnikiem D ≤ 1.

Eksponent macierzowy — operator exp(At) rozwiązujący układ dN/dt = AN; jego efektywne obliczenie jest sednem metod numerycznych transmutacji.

Samoochrona rezonansowa (SSK) — efekt redukcji efektywnego przekroju czynnego nuklidów z silnymi rezonansami w grubych próbkach; modelowany współczynnikiem SSK ≤ 1.

Sztywność (stiffness) — właściwość układu ODE z bardzo różnymi stałymi czasowymi; wymaga specjalnych niejawnych metod numerycznych.

Xenon pit — stan po awaryjnym wyłączeniu reaktora, gdy akumulacja ¹³⁵Xe uniemożliwia ponowne uruchomienie przez kilka godzin; czynnik istotny w awarii w Czarnobylu.

8 Podsumowanie dydaktyczne

  1. ChainSolver odpowiada na pytanie ilościowe: "ile danego nuklidu będzie po X godzinach napromieniowania i Y godzinach chłodzenia?" — uzupełniając topologiczną mapę ChainFindera.

  2. Matematycznym rdzeniem obliczeń jest eksponent macierzowy układu równań Batemana; jego efektywne i stabilne obliczenie jest kluczowym wyzwaniem implementacyjnym.

  3. Sztywność układu (stosunek stałych czasowych izomery/aktynidy ≈ 10²⁵) wymaga niejawnych metod numerycznych (CN, BDF, RADAU) — metody jawne (RK4) są praktycznie bezużyteczne dla realnych łańcuchów.

  4. Korekcje fizyczne (depresja strumienia, SSK) są kluczowe dla próbek grubych lub o dużej koncentracji silnych pochłaniaczy; ignorowanie ich może dać wyniki błędne o kilkadziesiąt procent.

  5. Wynik ChainSolvera jest ilościową symulacją dla zadanych założeń — nie "automatyczną prawdą" o składzie dowolnej próbki; zmiana widma neutronów, strumienia lub geometrii zmienia wyniki.

  6. Aplikacje medyczne (Lu-177, Tc-99m, Ra-223) wymagają uwzględnienia nie tylko dominującego nuklidu, ale i zanieczyszczeń promieniujących (Lu-177m, Mo-99) dla oceny czystości radionuklidowej.

  7. Fizyka reaktorów (xenon pit, zatrucie samarowe, bilans paliwa) jest naturalnym obszarem zastosowań ChainSolvera — z ważnym zastrzeżeniem, że dla pełnych cykli paliwowych dedykowane kody (ORIGEN, SERPENT) dają wyższą precyzję.

  8. Ten serwis udostępnia wersję webową ChainSolvera bezpośrednio z poziomu przeglądarki — bez instalacji Delphi ani kompilacji Fortranu — co obniża próg wejścia dla studentów i pozwala skupić się na fizyce zamiast na konfiguracji środowiska.

Porównanie ChainSolvera z innymi kodami transmutacji jądrowej

Zrozumienie miejsca ChainSolvera w szerszym ekosystemie oprogramowania jądrowego pomaga właściwie ocenić jego możliwości i ograniczenia:

ORIGEN2 i ORIGEN-S (Oak Ridge National Laboratory). Klasyczny kod transmutacji jądrowej, opracowany w ORNL w latach 1970. ORIGEN2 oblicza inwentarze radionuklidów w paliwie jądrowym jako funkcję czasu. Używa zintegrowanych bibliotek przekrojów czynnych dla standardowych typów reaktorów (PWR, BWR, AGR). Wbudowane biblioteki dla ~860 aktynidów, produktów rozszczepienia i produktów aktywacji. Ograniczenie: biblioteki są "zamrożone" dla konkretnych typów reaktorów; dla niestandardowych warunków napromieniowania precyzja spada. ORIGEN-S (wersja nowsza): może być sprzężony z kodami transportu neutronów (SCALE package).

FISPACT-II (UK Atomic Energy Authority). Następca FISPACT, rozwijany przez UKAEA. Specjalizuje się w aktywacji materiałów i produkcji izotopów. Obsługuje biblioteki EAF-2010, TENDL, JEFF. Popularny w Europie do obliczeń dla ITER/JET (fuzja) i reaktorów. Silna strona: automatyczne generowanie raportów o aktywności, ciepłu rozpadu i toksyczności dla materiałów aktywowanych.

SERPENT (VTT Finland). Kod Monte Carlo do fizyki reaktorów z wbudowanymi obliczeniami transmutacji (burnup). Łączy precyzję obliczeń transportowych z ewolucją składu paliwa. Używany intensywnie do analizy cyklu paliwowego i projektowania reaktorów SMR. Wymaga dużych zasobów obliczeniowych (Monte Carlo jest wolny). Nie nadaje się do "szybkich" obliczeń produkcji izotopu medycznego.

MCNP6 + CINDER90 (Los Alamos). MCNP: standard Monte Carlo do transportu neutronów i gamma. CINDER90: biblioteka transmutacji sprzężona z MCNP do obliczeń burnup. Używany w analizach bezpieczeństwa jądrowego (NRC) i projektach reaktorowych.

Porównanie tabelaryczne:

Cecha ChainSolver (ORIP_XXI) ORIGEN2 FISPACT-II SERPENT
Transport neutronów Nie (φ stały) Nie (φ stały) Nie (φ stały) Monte Carlo
Biblioteki danych ENDF/B, JEFF Wbudowane EAF, TENDL, JEFF JEFF, ENDF
Zakres nuklidów ~2800 ~860 ~2800+ wszystkie
Prędkość Bardzo szybka Szybka Szybka Wolna
Interface GUI + Web CLI CLI CLI/Python
Dostępność Open source Licencja RSICC Licencja Open source
Zastosowania Izotopy med., edukacja Paliwo jądrowe Aktywacja mat. Projektowanie reaktorów

Nisza ChainSolvera. ChainSolver wypełnia lukę między prostymi (ale nieprecyzyjnymi) kalkulatorami online a ciężkimi kodami przemysłowymi. Jest optymalny dla: szybkiej analizy konkretnego łańcucha transmutacji, produkcji izotopów medycznych/przemysłowych, scenariuszy niestandardowych (reaktory badawcze, własne harmonogramy napromieniowania), edukacji (zrozumienie numerycznych aspektów transmutacji).

Rozpad radioaktywny — przypomnienie podstaw dla kontekstu

Zrozumienie ChainSolvera wymaga solidnych podstaw z fizyki jądrowej. Poniżej zwięzłe przypomnienie kluczowych koncepcji:

Prawo rozpadu. Aktywność A(t) nuklidu o stałej rozpadu λ:
A(t) = A₀ · exp(−λt), gdzie λ = ln2/T½

Stałą λ wyznacza się eksperymentalnie lub oblicza z modeli jądrowych. T½ (okres półtrwania) rozciąga się od 10⁻²³ s (rezonanse jądrowe) do ~10²² lat (²⁰⁸Pb — praktycznie stabilny).

Kanały rozpadu. Każdy nuklid może rozpadać się przez kilka kanałów jednocześnie:

  • α (alfa): emisja ²He⁴, Z−2, A−4; typowe dla ciężkich nuklidów (Z > 82)
  • β⁻: emisja elektronu i antyneutrino; nuklidy bogatoneutronowe
  • β⁺ (EC): emission pozytonu lub wychwyt elektronu; nuklidy ubogoneutronowe
  • IT (izomeryczny): przejście z metastabilnego stanu wzbudzonego do stanu podstawowego (np. ⁹⁹ᵐTc → ⁹⁹Tc)
  • SF (spontaniczne rozszczepienie): dla ciężkich aktynidów (Cf-252, Pu-240)
  • Inne: n-emisja, p-emisja (rzadkie, dla egzotycznych nuklidów daleko od stabilności)

Serie promieniotwórcze. Naturalne szeregi radioaktywne:

  • Szereg uranowy: ²³⁸U → ... → ²⁰⁶Pb (14 ogniw)
  • Szereg aktynowy: ²³⁵U → ... → ²⁰⁷Pb (11 ogniw)
  • Szereg torowy: ²³²Th → ... → ²⁰⁸Pb (10 ogniw)
  • Szereg neptunowy: ²³⁷Np → ... → ²⁰⁹Bi (nie istnieje w naturze)

Każdy szereg jest typowym zastosowaniem ChainSolvera — modelowanie osiągania równowagi wiekowej (secular equilibrium) lub jej zaburzenia przez czynniki zewnętrzne.

Równowaga wiekowa. W długim łańcuchu rozpadów, gdy T½(rodzic) >> T½(córki): aktywności wszystkich ogniw wyrównują się do aktywności rodzica. To "równowaga wiekowa" (secular equilibrium). Przykład: ²²⁶Ra (T½ = 1600 lat) → ²²²Rn → ... → ²¹⁰Pb → ²¹⁰Bi → ²¹⁰Po → ²⁰⁶Pb. W skale zawierającej Ra i rówieśne Ra próbki gleby: aktywności wszystkich ogniw są równe aktywności Ra.

Produkcja izotopów dla przemysłu — zastosowania niemedyczne

Izotopy promieniotwórcze mają szerokie zastosowania poza medycyną. ChainSolver jest użyteczny do modelowania ich produkcji:

Irydu-192 (¹⁹²Ir) do radiografii przemysłowej. Iryd-192 (T½ = 73,8 dni) jest standardowym źródłem gamma do kontroli spawów, rur i odlewów w przemyśle naftowym, lotniczym i budowlanym. Produkcja: napromieniowanie ¹⁹¹Ir (n,γ) → ¹⁹²Ir. Zanieczyszczenie ¹⁹⁴Ir (T½ = 19,3h) pochodzi z ¹⁹³Ir (n,γ) → ¹⁹⁴Ir. ChainSolver modeluje optymalny czas napromieniowania i czas chłodzenia dla minimalizacji ¹⁹⁴Ir.

Kobaltu-60 (⁶⁰Co) do radioterapii i radiosterylizacji. Ko-60 (T½ = 5,27 roku) produkowany z Ko-59 (n,γ) → Ko-60. Jedna z najważniejszych aplikacji produkcji izotopów — źródła Ko-60 napędzają urządzenia do radioterapii (Gamma Knife, teleterapia), sterylizacji sprzętu medycznego i napromieniowania żywności. Reaktor MARIA Świerk produkuje Ko-60 na eksport. Optymalizacja: czas napromieniowania (dłuższy → wyższa aktywność Ko-60) vs. zanieczyszczenie Ko-58 (z Ni pochłaniającego neutrony).

Tulu-170 (¹⁷⁰Tm) do przenośnych defektoskopów. Tu-170 (T½ = 128,6 dni) jako źródło gamma/beta do kompaktowych przenośnych defektoskopów. Mniejszy i prostszy w obsłudze niż urządzenia z rentgenem, ale krótszy czas życia źródła. Produkcja: napromieniowanie Er-170 lub Tm-169 (n,γ) → Tm-170.

Plutonu-238 (²³⁸Pu) do RTG kosmicznych. Pu-238 (T½ = 87,7 roku) jest paliwem do Radioizotopowych Generatorów Termoelektrycznych (RTG) stosowanych w misjach NASA: Cassini, New Horizons, Perseverance. Produkowany przez napromieniowanie Np-237 → Np-238 (β⁻) → Pu-238. USA zaprzestały produkcji Pu-238 w 1988 roku i używały rosyjskich zapasów; produkcja wznowiona przez DoE w 2013 roku. ChainSolver może modelować ten łańcuch, choć produkcja Pu-238 w reaktorach badawczych jest nieopłacalna w skali — wymaga reaktorów mocy.

Promieniowanie jonizujące w naukach przyrodniczych poza fizyką

ChainSolver, jako narzędzie do modelowania ewolucji nuklidów, ma zastosowania wykraczające poza techniczne reaktory i produkcję izotopów:

Datowanie radiowęglowe (¹⁴C). ¹⁴C (T½ = 5730 lat) produkowany przez kosmiczne promieniowanie neutronowe: ¹⁴N (n,p) → ¹⁴C. Proporcja ¹⁴C/¹²C w atmosferze jest stała (z korektą kalibracyjną). Po śmierci organizmu: ¹⁴C zanika bez odnawiania. Wiek próbki wyznacza się z stosunku ¹⁴C/¹²C. ChainSolver może modelować zanikanie ¹⁴C jako prosty przykład jednoetapowego rozpadu — dobry punkt startowy dla ćwiczeń.

Geologia izotopowa. Systemy Rb-Sr, U-Pb, K-Ar, Sm-Nd do datowania skał i minerałów. ChainSolver może modelować narastanie izotopu radiogenicznego (np. ⁸⁷Sr z rozpadów ⁸⁷Rb) przez geologiczny czas. Dla systemów takich jak Pb-Pb (dwa kadle uranowe konwergujące do Pb) użycie ChainSolvera z wieloetapową historią termiczną ilustruje złożoność geochronologii.

Ochrona środowiska. Skażenie środowiska nuklidami sztucznymi (Cs-137, Sr-90, Am-241 z testów broni jądrowej) lub naturalnymi (Ra-226 w wodach kopalnianych). Modelowanie losu radionuklidów w środowisku wymaga rozumienia ich ewolucji w czasie. ChainSolver dostarcza odpowiedzi na pytanie "jakie nuklidy będą obecne za 50 lat" dla zadanego składu początkowego.

Astrofizyka nuklearna. Procesy r i s w gwiazdach tworzą ciężkie pierwiastki przez wychwyt neutronów. Łańcuchy transmutacji w ChainSolverze odpowiadają fizycznie tym samym procesom, choć parametry (strumień neutronów, czas, temperatura) są radykalnie różne od ziemskich reaktorów. Modelowanie procesu r (rapid neutron capture) w wyniku zderzenia gwiazd neutronowych (kilonova) to jedno z aktywnych pól badań astrofizyki nuklearnej.

Instrukcja krok po kroku — ćwiczenie z Co-60

Poniżej instrukcja szczegółowa dla typowej sesji edukacyjnej z ChainSolverem:

Cel: Wyznaczenie optymalnego czasu napromieniowania próbki Co-59 dla uzyskania maksymalnej aktywności Ko-60 z minimalnymi zanieczyszczeniami.

Krok 1: Uruchom ChainFinder dla Co-59. Wpisz "Co-59" jako nuklid startowy. Ustaw strumień φ = 10¹³ n/(cm²·s) (typowy dla reaktora badawczego). Wybierz głębokość 3. Zidentyfikuj kluczowe nuklidy: Co-60 (pożądany), Ni-60, Ni-61, Fe-56 (zanieczyszczenia od normalnych metali towarzyszących).

Krok 2: Uruchom ChainSolver — krótkie napromieniowanie. Ustaw: nuklid startowy Co-59, jeden etap, t = 7 dni, φ = 10¹³. Odczytaj z tabeli: ile % Co-60 w porównaniu do całości? Jaka jest aktywność właściwa Co-60?

Krok 3: Wydłuż czas. Zmień etap na t = 30 dni, potem 90 dni, 180 dni. Dla każdego odczytaj: frakcję masową Co-60, aktywność Co-60, zanieczyszczenia (Co-57, Co-58 z powolnych reakcji). Narysuj wykres aktywności Co-60 vs. czas napromieniowania.

Krok 4: Dodaj etap chłodzenia. Po 30 dniach napromieniowania dodaj 7 dni chłodzenia (φ=0). Jak zmienia się skład? Które krótkotrwałe zanieczyszczenia zanikają podczas chłodzenia? Jak zmienia się czystość radiochemiczna?

Krok 5: Zinterpretuj wyniki. Dlaczego aktywność Co-60 asymptotycznie dąży do nasycenia zamiast rosnąć nieograniczenie? (Odpowiedź: pochłanianie Co-60 przez Co-60+n → Co-61, oraz rozpad Co-60, który skraca czas życia izotopu w wiązce). Co to oznacza dla optymalnej długości napromieniowania reaktorowego?

ChainSolver a kodeksy ochrony radiologicznej

Obliczenia ChainSolvera mają bezpośrednie znaczenie dla planowania ochrony radiologicznej w laboratoriach produkujących izotopy:

Klasyfikacja odpadów radioaktywnych. Wyniki ChainSolvera (stężenia nuklidów, aktywności, czas zaniku) są podstawą do klasyfikacji odpadów radiochemicznych. Odpady o aktywności poniżej limitów zwolnienia (clearance levels — IAEA RS-G-1.7) mogą być traktowane jak konwencjonalne odpady. ChainSolver pomaga odpowiedzieć: "po ilu dniach aktywność próbki spadnie poniżej limitu zwolnienia?" — co minimalizuje koszty składowania.

Dawki dla personelu — obliczenia konserwatywne. Planowanie napromieniowania tarcz w reaktorze: nie tylko interesuje nas produkcja pożądanego izotopu, ale i promieniowanie gamma z zanieczyszczeń (które narażają personel manipulujący tarczą po wyjęciu z reaktora). ChainSolver daje aktywności wszystkich nuklidów tuż po wyjęciu z reaktora i w kolejnych dniach — podstawa do obliczenia dawki według modelu: D ≈ Σ_i A_i · Γ_i / r², gdzie Γ_i to czynnik dawkowy izotopu i.

Krytyczność i bezpieczeństwo jądrowe. ChainSolver modeluje narastanie Pu-239 z U-238 w reaktorze. Rozumienie tej dynamiki jest istotne dla oceny ryzyka krytyczności w instalacjach przetwarzania paliwa — choć do właściwych obliczeń krytyczności używa się dedykowanych kodów (MCNP, KENO).

Monitoring środowiskowy. Efluenty z zakładów produkcji izotopów zawierają śladowe ilości radionuklidów. Rozumienie ich ewolucji w czasie (jakie nuklidy będą obecne 24h, 7 dni, rok po emisji) pomaga w projektowaniu systemów monitoringu i interpretacji pomiarów.

NCBJ Świerk — kontekst polski dla obliczeń transmutacji

Polska ma unikalną w Europie infrastrukturę do produkcji izotopów jądrowych, związaną z reaktorem MARIA w Narodowym Centrum Badań Jądrowych (NCBJ) w Świerku:

Reaktor MARIA. Jedyny działający reaktor badawczy w Polsce. Moc: 20 MW (cieplnych). Typ: basenowy, wysokokontrastowy (HEU, a od 1999 LEU). Strumień neutronów termicznych w centrum: do ~5×10¹⁴ n/(cm²·s). Reaktor pracuje w cyklach roboczych z przerwami na ładowanie paliwa i konserwacje. Planowe modernizacje — MARIA ma pracować do co najmniej 2030 roku.

Produkcja izotopów w MARIA. NCBJ produkuje i eksportuje do Europy izotopy medyczne i przemysłowe: Mo-99/Tc-99m (prekursor generatora technetowego), I-131 (terapia tarczycy), Co-60 (radioterapia, sterylizacja), Lu-177 (terapia nowotworów), Ir-192 (defektoskopia przemysłowa), Se-75 (defektoskopia), Ho-166 (terapia celowana). Każda produkcja wymaga planowania harmonogramu napromieniowania — dokładnie to robi ChainSolver.

Centrum Radiofarmaceutyczne POLATOM. IEA POLATOM (teraz część NCBJ) jest wytwórcą i dystrybutorem radiofarmaceutyków. Produkuje zestawy do znakowania technetem (DMSA, DTPA, HMDP, MAG3) dla diagnostyki medycznej (scyntygrafia nerek, kości, perfuzji). Produkuje Lu-177 chlorek do radiofarmaceutyków terapeutycznych.

MARIA a obliczenia ChainSolvera. Wiedząc, że MARIA produkuje Ko-60 z pełną mocą przy φ ≈ 3×10¹⁴ n/(cm²·s), student może sprawdzić w ChainSolverze: jaka jest aktywność Ko-60 po 20-dniowym napromieniowaniu? Jak zmienia się czystość? Jakie są czasy transportu i jak wpływają na aktywność produktu dostarczanego do klienta? To konkretne, praktyczne ćwiczenie łączące lokalny kontekst z obliczeniami numerycznymi.

Połączenia między kalkulatorami ORIP_XXI

Typowy scenariusz pracy. Student lub badacz pracuje sekwencyjnie z trzema kalkulatorami:

  1. NKE: "Jakie mam dane dla Lu-177?" → sprawdza T½ (6,647 dni), kanały rozpadu (β⁻ 100%, γ 208 keV, 113 keV), przekroje czynne (σ_a Lu-176 = 2100 barn dla n termicznych). To dane wejściowe dla kolejnych kroków.

  2. ChainFinder: "Jakie nuklidy wchodzą w grę przy napromieniowaniu Lu-176?" → mapa łańcucha: Lu-176 → Lu-177 (główne), Lu-176 + n → Lu-177, Lu-177 + n → Lu-178 (zanieczyszczenie), Lu-177m (T½ = 160,4 dni, problem!). Wizualizacja grafu pokazuje, że Lu-177m może być kłopotliwym zanieczyszczeniem.

  3. ChainSolver: "Ile Lu-177 vs Lu-177m powstaje przy φ = 10¹⁴ przez 5 dni?" → ilościowa odpowiedź: frakcja masowa, aktywność, czystość radionuklidowa. Porównanie z wymaganiami Farmakopei Europejskiej.

Pętle informacyjne. Wyniki ChainSolvera często generują nowe pytania do ChainFindera i NKE: "Skąd pochodzi nieoczekiwane zanieczyszczenie Lu-173? Który nuklid jest jego rodzicem?" → powrót do ChainFindera z głębszym przeszukiwaniem.

Zewnętrzne bazy danych jako weryfikacja. Po obliczeniu ChainSolverem warto porównać wyniki z eksperymentalnymi danymi z baz:

  • NuDat 3 (Brookhaven National Laboratory): baza danych nuklidów z eksperymentalnymi danymi rozpadowymi, mierzalnymi przekrojami czynnymi i widmami gamma. Dostępna pod nndc.bnl.gov/nudat3.
  • Livechart of Nuclides (IAEA, Wiedeń): interaktywna mapa nuklidów z danymi zaczerpniętymi z ENSDF. Dostępna pod nds.iaea.org/relnsd/vcharthtml/VChartHTML.html.
  • OECD/NEA JANIS: Java-based Nuclear Data Information System — przeglądarka bibliotek danych jądrowych (ENDF/B, JEFF, JENDL). Pozwala porównać przekroje czynne z różnych bibliotek.

Weryfikacja wybranych wyników obliczeń w eksperymentalnych bazach danych buduje intuicję co do dokładności ChainSolvera i wiarygodności używanych bibliotek danych jądrowych.

Znaczenie dydaktyczne sekwencji NKE → ChainFinder → ChainSolver. Ta sekwencja odwzorowuje naturalny tok myślenia badacza: od danych atomowych, przez topologię sieci transmutacji, do ilościowej symulacji. Uczy podejścia "bottom-up" do złożonych problemów jądrowych — zamiast traktowania wyników jako "magicznej czarnej skrzynki", student rozumie każdy krok obliczeń i może weryfikować wyniki cząstkowe na każdym poziomie.

Dla polskiego studenta fizyki lub inżynierii jądrowej — szczególnie wobec planowanej budowy elektrowni jądrowej i rozwijającego się sektora medycyny nuklearnej — opanowanie tych narzędzi jest inwestycją w konkretne, poszukiwane kompetencje zawodowe. NCBJ w Świerku, instytucje regulacyjne (PAA), instytuty badawcze i przyszłe elektrownie jądrowe potrzebują specjalistów rozumiejących transmutację jądrową na poziomie obliczeniowym — nie tylko opisowym. Kalkulator webowy ChainSolvera dostępny w tym serwisie pozwala zacząć tę przygodę już dziś — bez instalacji specjalistycznego oprogramowania i z wbudowanymi presetsami dydaktycznymi obejmującymi typowe scenariusze produkcji izotopów.

Dodatkowe materiały multimedialne

  • Dokumentacja ORIGEN2: ORNL/CCC-371, Oak Ridge National Laboratory
  • Dokumentacja ODEPACK/LSODA: Hindmarsh, A.C., "ODEPACK, A Systematized Collection of ODE Solvers", 1983
  • ENDF/B-VIII.1 — baza danych jądrowych: nndc.bnl.gov/endf-releases
  • IAEA Nuclear Data Section: nds.iaea.org
  • European Pharmacopoeia monographs on radiopharmaceuticals: edqm.eu

Narzędzia interaktywne

Ćwiczenia praktyczne