An elementary method to determine the critical mass of a sphere of fissile material based on a separation of neutron transport and nuclear reaction processes

Dit pedagogische artikel presenteert een vereenvoudigde methode om de kritieke massa van een bolvormige splijtstof te berekenen door neutronentransport en kernreacties te scheiden, wat resulteert in een nauwkeurige schatting zonder de diffusievergelijking op te hoeven lossen.

De kern

De Kern van het verhaal: Hoe groot moet een kogel zijn om te ontploffen?

Stel je voor dat je een kogel van splijtstof (zoals uranium of plutonium) hebt. De vraag is: Hoe groot moet die kogel zijn voordat hij zichzelf in een kettingreactie kan houden en ontploft? Dit heet de "kritieke massa".

Als de kogel te klein is, ontsnappen er te veel deeltjes (neutronen) naar buiten en stopt de reactie. Als hij groot genoeg is, blijven er genoeg deeltjes binnen om nieuwe ontploffingen te veroorzaken.

Lamoreaux heeft een nieuwe, simpele manier bedacht om dit te berekenen. In plaats van ingewikkelde wiskunde te gebruiken die alleen natuurkundigen begrijpen, gebruikt hij een simpele "wiskundige schatting" die gebaseerd is op twee simpele ideeën:

1. De reactie: Hoeveel nieuwe deeltjes worden er gemaakt?
2. Het verlies: Hoeveel deeltjes ontsnappen er?

---

De Analogie: De "Neutronen-Feestzaal"

Om dit te begrijpen, laten we een feestzaal voorstellen.

1. De Deeltjes (De Gasten)

De neutronen zijn de gasten op het feest. Ze rennen rond in de zaal (de splijtstofkogel).

• Soms botsen ze tegen een muur en worden ze gevangen (geabsorbeerd).
• Soms botsen ze tegen een speciale machine die, als er een gast tegenaan botst, nieuwe gasten produceert (dit is kernsplijting).
• Soms lopen ze gewoon naar de uitgang en verlaten de zaal (lekken).

2. Het Doel: Een evenwicht

Voor een "kritieke" kogel moet het aantal gasten in de zaal constant blijven.

• Als er 100 gasten zijn, moeten er precies 100 nieuwe gasten worden gemaakt (of de oude blijven) om degenen die weglopen of gevangen worden, te vervangen.
• Als er minder dan 100 nieuwe gasten zijn, sterft het feest uit (subkritisch).
• Als er meer zijn, groeit het feest explosief (supercritisch = ontploffing).

---

De Twee Delen van de Oplossing

Lamoreaux splitst het probleem op in twee simpele stappen, net als het oplossen van een raadsel.

Stap 1: De "Reis" (De Wiskunde van de Reactie)

Eerst kijken we alleen naar wat er gebeurt als een gast door de zaal loopt.

• Hoeveel kans is er dat een gast een nieuwe gast maakt voordat hij de zaal verlaat of gevangen wordt?
• Dit hangt af van hoe "dicht" de zaal is gevuld en hoe vaak de machines werken.
• De auteur berekent een minimale reistaf (een afstand) die een gast moet afleggen voordat hij genoeg nieuwe gasten heeft gegenereerd om het verlies te compenseren.

Vergelijking: Stel je voor dat je in een bos loopt. Je moet een bepaalde afstand lopen om genoeg bomen te vinden om een vuurtje te stoken. Als je te kort loopt, heb je niet genoeg hout. De auteur berekent precies hoe lang die wandeling moet zijn.

Stap 2: De "Dwarslopers" (De Wiskunde van het Bewegen)

Nu komt het leuke deel. De gasten (neutronen) lopen niet in een rechte lijn naar de uitgang. Ze botsen constant tegen muren en andere gasten aan. Ze slenteren of dwalen rond. Dit heet een "random walk" (willekeurige wandeling).

• Omdat ze dwalen, is de afstand die ze afleggen veel langer dan de rechte lijn naar de uitgang.
• Door te dwalen, blijven ze langer in de zaal. Dit is goed! Het geeft ze meer kans om nieuwe gasten te maken.
• De auteur gebruikt simpele statistiek om te berekenen: "Als een gast een bepaalde afstand moet dwalen om genoeg nieuwe gasten te maken, hoe groot moet de zaal dan zijn zodat hij er niet uitvalt?"

Vergelijking: Stel je voor dat je een bal in een kamer met honderden pilaartjes gooit. De bal stuitert overal tegenaan. Als de kamer klein is, valt de bal snel uit een raam. Als de kamer groot is, botst de bal zo vaak tegen pilaartjes dat hij de kamer niet uitkomt voordat hij genoeg "werk" heeft gedaan. De auteur berekent de perfecte grootte van die kamer.

---

Het Resultaat: Een Simpele Formule

Het mooie aan deze methode is dat je geen ingewikkelde differentiaalvergelijkingen (de zware wiskunde die normaal gesproken nodig is) hoeft op te lossen.

De auteur komt uit op een simpele formule die zegt:

De straal van de kogel is evenredig met de afstand die een deeltje moet dwalen, vermenigvuldigd met een correctiefactor.

Toen hij deze simpele formule toepaste op echte materialen (Plutonium-239 en Uranium-235), bleek het resultaat verbazingwekkend nauwkeurig te zijn.

• Voor Plutonium-239: De berekende kritieke massa was ongeveer 10 kg.
• De super-computers (die duizenden jaren aan berekeningen doen) zeggen: 10,2 kg.

Dat is een verschil van minder dan 3%! Dat is alsof je de afstand van Amsterdam naar Rotterdam schat met je pasjes en uitkomt op 57 km, terwijl de GPS 58,5 km aangeeft.

Waarom is dit belangrijk?

1. Educatie: Het helpt studenten (en ons) om te begrijpen waarom iets ontploft, zonder dat ze eerst 10 jaar wiskunde moeten studeren. Het maakt de fysica "voelbaar".
2. Controle: Het kan gebruikt worden om te checken of ingewikkelde computerprogramma's (zoals MCNP) wel goed werken. Als je simpele schatting en de supercomputer niet overeenkomen, weet je dat er iets mis is.
3. Veiligheid: Het helpt bij het begrijpen van hoe verrijking (zuiverheid van het materiaal) de kritieke massa beïnvloedt. Bijvoorbeeld: als je uranium minder zuiver maakt, moet de kogel veel groter zijn om te ontploffen.

Conclusie

Lamoreaux heeft laten zien dat je de gevaarlijke en complexe wereld van kernontploffingen kunt begrijpen met simpele logica en een beetje statistiek. Het is als het oplossen van een puzzel: als je weet hoeveel stukjes je nodig hebt om het plaatje compleet te maken, en je weet hoe de stukjes bewegen, kun je de grootte van het plaatje berekenen zonder de hele doos te openen.

Het bewijst dat soms de simpelste manier van kijken naar een probleem de meest accurate antwoorden geeft.

Technische samenvatting

Titel: Een elementaire methode om de kritieke massa van een bol splijtbaar materiaal te bepalen, gebaseerd op een scheiding van neutronentransport en kernreactieprocessen

1. Probleemstelling

Het bepalen van de kritieke massa van een bol van splijtbaar materiaal (zoals $^{235}\text{U}$ of $^{239}\text{Pu}$) vereist doorgaans complexe numerieke methoden, zoals het oplossen van de diffusievergelijking of het gebruik van geavanceerde Monte Carlo-simulaties (bijv. MCNP). Hoewel deze methoden zeer nauwkeurig zijn, vormen ze een hoge drempel voor educatieve doeleinden en voor snelle, "uit het hoofd"-schattingen (back-of-the-envelope estimates).

De auteur stelt een vereenvoudigde, pedagogische methode voor die de fysica van het probleem ontleedt zonder geavanceerde wiskunde te vereisen. Het doel is om de kritieke straal af te leiden door het probleem te scheiden in twee delen:

1. De kernreactieprocessen (absorptie en splijting).
2. De geometrie en mechanica van het neutronentransport (diffusie/robuust wandelen).

2. Methodologie

De methode maakt gebruik van elementaire calculus en statistische argumenten. De kern van de benadering is het definiëren van de kritieke toestand als het moment waarop het aantal vrije neutronen ($N$) in de bol constant blijft in de tijd.

A. Scheiding van Processen:

• Kernreacties: De auteur definieert een kritieke wandelafstand $\ell_c$ die een neutron moet afleggen voordat het wordt geabsorbeerd of een splijting veroorzaakt. Deze afstand wordt bepaald door de totale absorptiewaarschijnlijkheid (splijting + niet-splijtingsabsorptie) en het aantal vrijgekomen neutronen per splijting ($\nu$).
• Transport: De relatie tussen deze lineaire wandelafstand $\ell_c$ en de fysieke straal van de bol ($R_c$) wordt bepaald door een willekeurig wandelmodel (random walk) gebaseerd op verstrooiing.

B. Afleiding van de Kritieke Straal:

1. Kritieke Afstand ($\ell_c$):
   Op basis van de cross-sections ($\sigma_f$ voor splijting, $\sigma_a$ voor absorptie, en $\sigma_0 = \sigma_f + \sigma_a$) en het aantal neutronen per splijting ($\nu$), wordt de voorwaarde voor kritikaliteit afgeleid. Een neutron moet een afstand $\ell$ afleggen zodat de productie van nieuwe neutronen de verliezen door absorptie en lekken compenseert. Dit leidt tot de vergelijking:
   $$\ell_c = -\frac{1}{n\sigma_0} \ln\left(1 - \frac{\sigma_0}{\nu \sigma_f - \sigma_0}\right)$$
   waarbij $n$ de nucleaire dichtheid is.

2. Random Walk en Geometrie:
   Neutronen ondergaan een willekeurig wandelproces met een gemiddelde vrije weglengte tussen verstrooiingen $\ell_s = 1/(n\sigma_{tr})$, waarbij $\sigma_{tr}$ de transportverstrooiingscross-section is.
   Volgens de theorie van willekeurige wandelingen is het aantal stappen $N_s$ nodig om een geometrische dimensie (de straal $R$) te bereiken gerelateerd aan $\ell_s$. De auteur introduceert een correctiefactor $\epsilon$ (afgeleid van de verdeling van afstanden in een bol en statistische correcties voor de wortel van het kwadraat van de verplaatsing) om de relatie tussen de totale afgelegde afstand en de straal te koppelen:
   $$R_c^2 = \epsilon^2 \ell_c \ell_s$$
   Hieruit volgt de definitieve formule voor de kritieke straal:
   $$R_c = \epsilon \sqrt{\ell_c \ell_s}$$
   De kritieke massa wordt vervolgens berekend als $m_c = \frac{4}{3}\pi \rho R_c^3$.

3. Transport Cross-section:
   Een belangrijke bijdrage is de benadering voor de transportcross-section $\sigma_{tr}$. De auteur toont aan dat rekening moet worden gehouden met de voorwaartse piek in elastische verstrooiing. De benadering die wordt gebruikt is:
   $$\sigma_{tr} \approx \sigma_e + \frac{1}{2}\sigma_i$$
   waarbij $\sigma_e$ de elastische en $\sigma_i$ de inelastische cross-section is.

3. Belangrijkste Resultaten

De methode is toegepast op $^{235}\text{U}$, $^{238}\text{U}$ en $^{239}\text{Pu}$ met gebruik van data uit ENDF/B-VIII.0.

• Nauwkeurigheid: De berekende kritieke massa's komen verrassend nauwkeurig overeen met geavanceerde MCNP6-simulaties (binnen enkele procenten).
  • $^{239}\text{Pu}$: Berekende massa is 10,0 kg (bij 1,9 MeV neutronenenergie), vergeleken met 10,2 kg voor MCNP6.
  • $^{235}\text{U}$: Berekende massa is 46,5 kg (bij 93,71% verrijking), vergeleken met 46,4 kg voor MCNP6.
• Vergelijking met Oppenheimer-Bethe: De resultaten zijn nauwkeuriger dan de historische Oppenheimer-Bethe formule (uit 1943), die bij $^{235}\text{U}$ een massa van 55 kg voorspelde (te hoog) en bij $^{239}\text{Pu}$ 8,6 kg (te laag).
• Invloed van Verrijking: De methode kan worden toegepast op verrijkte materialen. Voor $^{235}\text{U}$ met lagere verrijking (bijv. 20%) is de kritieke massa aanzienlijk groter, en de auteur benadrukt dat de effecten van neutronenmoderatie (energieverlies door verstrooiing) kleiner zijn dan vaak wordt gedacht, omdat neutronen snel worden geabsorbeerd of vervangen door nieuwe, hoog-energetische splijtingsneutronen.
• Tijdschaal: De auteur schat de tijd die nodig is om een mol $^{239}\text{Pu}$ te laten splijten op ongeveer 0,57 microseconden, wat de microseconden-schaal van kernexplosies bevestigt.

4. Bijdragen en Innovaties

• Pedagogische Toegang: De methode maakt het mogelijk om kritieke massa's te begrijpen en te berekenen zonder de complexe wiskunde van de diffusievergelijking of numerieke simulaties.
• Scheiding van Problemen: Door kernreacties en transportmechanica expliciet te scheiden, wordt de fysica duidelijker. De enige koppeling is de totale afgelegde afstand.
• Nieuwe Randvoorwaarde: In de appendices wordt een nieuwe behandeling van randvoorwaarden voor de diffusievergelijking gepresenteerd die leidt tot een kwadratische vergelijking in plaats van een transcendentale, wat de vergelijking met de elementaire methode vereenvoudigt en een sterke overeenkomst toont ($\epsilon \approx 0,89$ vs. $\epsilon' \approx 0,875$).
• Toepasbaarheid: De methode is uitbreidbaar naar onzuivere materialen, isotopische mengsels en andere geometrieën (zoals staven) voor snelle schattingen en validatie van complexe simulaties.

5. Betekenis en Conclusie

Dit artikel demonstreert dat een sterk vereenvoudigd model, gebaseerd op elementaire statistiek en een duidelijke fysieke scheiding van processen, kan leiden tot resultaten die binnen de foutmarges van geavanceerde, dure simulaties vallen.

De belangrijkste conclusie is dat de kritieke massa voornamelijk wordt bepaald door de balans tussen neutronenproductie en -verlies via absorptie en lekken, en dat de details van de neutronenspectrum-moderatie minder kritiek zijn dan vaak wordt aangenomen voor snelle neutronen in een kritieke assemblage. Deze aanpak biedt een waardevol hulpmiddel voor onderwijs, voor het begrijpen van non-proliferatiekwesties, en als snelle "sanity check" voor ingenieurs die werken met complexe neutronische codes.