Generation of fission yield covariance matrices and its application in uncertainty analysis of decay heat

Dit artikel beschrijft de generatie van covariantiematrices voor splijtingsopbrengsten met behulp van de GLS-updating-methode, wat leidt tot een aanzienlijke vermindering van de onzekerheid in de berekende vervalwarmte voor thermische neutroneninduced splijting van $^{235}$U.

De kern

Stel je voor dat een kernreactor een enorme, complexe machine is die elektriciteit maakt door atomen te "splijten" (verpletteren). Wanneer je de reactor stopt, is het alsof je de motor uitschakelt, maar de machine blijft nog even heet van de restwarmte. Deze warmte noemen we vervalwarmte. Als je deze warmte niet goed kunt voorspellen en beheersen, kan het gevaarlijk worden voor de veiligheid van de centrale.

Om deze warmte te berekenen, moeten wetenschappers een soort "recept" volgen. Ze kijken naar duizenden verschillende stukjes atoom (splijtingsproducten) die ontstaan bij het splijten. Elk van deze stukjes heeft een eigen levensduur en geeft een beetje warmte af terwijl het vervalt.

Het probleem waar dit papier over gaat, is dat de "receptboeken" (de databanken met atoomgegevens) tot nu toe een groot gat hadden: ze wisten precies hoeveel warmte elk stukje atoom zou geven, maar ze wisten niet hoe onzeker die gegevens waren en ze wisten niet of de fouten bij de ene atoomsoort samenhangen met de andere.

Hier is een simpele uitleg van wat de auteurs hebben gedaan, met behulp van een paar creatieve vergelijkingen:

1. Het probleem: De onbekende vrienden

Stel je voor dat je een orkest hebt met 100 muzikanten. Je weet dat elke muzikant een beetje vals kan spelen (onzekerheid). Maar in de oude databanken werd aangenomen dat als de fluitist een noot mist, de trompettist dat niet doet. Ze dachten dat alle fouten los van elkaar stonden.

In werkelijkheid is dat niet zo. Als de dirigent (de natuurwetten) een fout maakt in de snelheid, spelen alle muzikanten tegelijk een beetje te snel. De fouten zijn gekoppeld. Als je dit niet meeneemt in je berekening, denk je dat het orkest veel meer kan variëren dan dat het eigenlijk kan. Je krijgt dus een veel te pessimistisch beeld van de onzekerheid.

De auteurs van dit papier wilden die "gekoppelde fouten" (de covariantie) berekenen, omdat die cruciaal zijn voor een eerlijke inschatting van de veiligheid.

2. De oplossing: De "Super-Detective" methode (GLS)

De auteurs gebruikten een wiskundige techniek genaamd Generalized Least Squares (GLS). Je kunt dit zien als een slimme detective die drie verschillende receptboeken (ENDF/B, JENDL en JEFF) naast elkaar legt.

De detective heeft een paar harde regels (de "wetten van de natuur"):

• De balans: Als je een atoom splijt, moet de totale hoeveelheid materie en lading precies hetzelfde blijven als voor de splijting. (Net als bij een weegschaal: wat erin gaat, moet eruit komen).
• De ketting: De atomen die ontstaan, vervallen in een kettingreactie naar stabiele atomen. De detective kijkt of de eindresultaten van die kettingen kloppen met wat we in het echt hebben gemeten.

De detective neemt de oude, onvolmaakte gegevens en past ze een klein beetje aan zodat ze perfect voldoen aan deze natuurwetten. Tegelijkertijd berekent hij: "Als deze atoomsoort een beetje meer warmte geeft, betekent dat dan dat die andere atoomsoort er iets minder van geeft?"

Dit proces creëert een correlatiematrix. In plaats van een lijstje met losse foutmarges, krijgen we een web van verbanden. Het is alsof we van een lijstje met losse losse schroeven gaan naar een ingewikkeld, maar perfect op elkaar afgestemd mechanisme.

3. Het resultaat: Van paniek naar rust

Wat gebeurde er toen ze deze nieuwe, slimme berekeningen toepasten op de warmteberekening?

• Vroeger (zonder de nieuwe methode): De onzekerheid was enorm. Het leek alsof de warmte op elk moment 10% of meer kon afwijken. Dit kwam vooral omdat ze dachten dat elke atoomsoort een eigen, losse fout had.
• Nu (met de nieuwe methode): Omdat ze zagen dat de fouten elkaar vaak opheffen (als de ene te hoog is, is de andere te laag), viel de totale onzekerheid drastisch.
  • Bij korte tijd (0,1 seconde na stoppen) daalde de onzekerheid van ongeveer 10% naar slechts 5%.
  • Bij langere tijd (na 100.000 seconden) was de onzekerheid zelfs gedaald tot ongeveer 1%.

Het is alsof je eerst dacht dat je een brug misschien wel 10% te zwaar zou zijn, maar toen je ontdekte dat de balken elkaar steunden, realiseerde je je dat de brug eigenlijk 99% veilig is.

4. Waarom is dit belangrijk?

Voor de veiligheid van kerncentrales is het cruciaal om te weten hoeveel warmte er nog overblijft.

• Als je denkt dat de onzekerheid groot is, moet je je koelsystemen enorm overdimensioneren (te duur en inefficiënt).
• Als je de onzekerheid juist goed kunt berekenen (zoals deze auteurs deden), kun je de systemen veiliger en efficiënter ontwerpen.

Kortom:
De auteurs hebben een slimme wiskundige "rekenmachine" gebouwd die de oude atoomgegevens heeft "opgepoetst" door ze te laten voldoen aan de wetten van de natuur. Hierdoor konden ze laten zien dat de onzekerheid in de warmteberekening veel kleiner is dan men dacht. Ze hebben de "losse schroeven" van de oude data omgezet in een "perfect op elkaar afgestemd orkest", waardoor we veel zekerder zijn over de veiligheid van kernenergie.

Technische samenvatting

Titel: Generatie van covariantiematrices voor splijtingsopbrengst en toepassing in onzekerheidsanalyse van vervalwarmte

Auteurs: Wendi Chen, Tao Ye, et al. (Instituut voor Toegepaste Fysica en Computatiele Wiskunde, Peking, China)

---

1. Het Probleem

Onzekerheidsanalyse is van cruciaal belang in nucleair engineering, met name voor de berekening van vervalwarmte (decay heat) van splijtingsproducten. Deze warmte is essentieel voor reactorontwerp, veiligheidsanalyses en het beheer van nucleair afval.

• Huidige tekortkoming: Bestaande geëvalueerde nucleaire databanken (zoals ENDF/B-VIII.0, JENDL-5 en JEFF-3.3) bevatten wel de beste schattingen en onzekerheden (standaardafwijkingen) voor splijtingsopbrengsten (fission yields), maar geen covariantiematrices.
• Gevolg: Zonder covarianties (correlaties tussen verschillende isotopen) worden de onzekerheden in de totale vervalwarmte vaak overschat. De correlaties tussen splijtingsproducten met dezelfde massagetal of lading zijn fysiek noodzakelijk om de totale onzekerheid correct te kwantificeren.

2. Methodologie

De auteurs hebben een geavanceerde aanpak ontwikkeld om covariantiematrices te genereren en deze toe te passen in onzekerheidsanalyses.

A. Generatie van Covariantiematrices (GLS-benadering)
Er werd gebruikgemaakt van de Generalized Least Squares (GLS) update-methode om de onafhankelijke splijtingsopbrengsten (IFY) en hun covarianties te evalueren.

• Invoerdata: Gebruik van bestaande data uit ENDF/B-VIII.0, JENDL-5 en JEFF-3.3.
• Fysieke Beperkingen (Constraints): De GLS-update werd uitgevoerd onder strikte fysieke behoudswetten:
  1. Behoud van massagetal en lading: De som van de opbrengsten moet overeenkomen met het oorspronkelijke kernmassa en -lading (correctie voor prompte neutronen).
  2. Normalisatie: De som van de opbrengsten (exclusief lichte geladen deeltjes) moet gelijk zijn aan 2.
  3. Symmetrie van lading: Bij lage-energie neutroneninduced splijting moet de ladingverdeling symmetrisch zijn rond de kernlading.
  4. Ketenopbrengst (Chain Yield) data: Gebruik van experimentele ketenopbrengstdata (van England & Rider en Nichols et al.) als sterke beperking. Dit zorgt voor een sterke negatieve correlatie tussen isotopen met hetzelfde massagetal.
• Resultaat: De methode genereert niet alleen de covariantiematrix, maar past ook licht de oorspronkelijke IFY-waarden aan om te voldoen aan de beperkingen.

B. Sommatieberekening en Generalized Perturbation Theory

• Vervalwarmteberekening: Er werd een "Fission Pulse Decay Heat" (FPDH) berekening uitgevoerd voor thermische neutroneninduced splijting van $^{235}$U. Dit omvatte zowel de warmte door lichte deeltjes ($\beta$-straling) als elektromagnetische straling ($\gamma$-straling).
• Onzekerheidspropagatie: De onzekerheid in de vervalwarmte werd berekend met behulp van de Generalized Perturbation Theory (Sandwich-formule). Hierbij werden onzekerheden meegenomen van:
  • Splijtingsopbrengst (IFY)
  • Vervalenergie
  • Vervalconstante
  • Vertakkingsverhoudingen (branching ratios)

3. Belangrijkste Bijdragen

1. Eerste generatie van volledige covariantiematrices: Voor de eerste keer werden covariantiematrices gegenereerd voor IFY's in de drie belangrijkste databanken (ENDF, JENDL, JEFF) met behulp van de GLS-methode en fysieke beperkingen.
2. Kwantificering van correlaties: De studie toont aan dat de correlaties voornamelijk negatief zijn tussen isotopen met naburige massagetallen, vooral in gebieden waar de ketenopbrengst afwijkt van de onafhankelijke opbrengst door kruisende massaketens.
3. Validatie van onzekerheidsreductie: Er werd aangetoond dat het gebruik van deze covariantiematrices de totale onzekerheid in vervalwarmteberekeningen drastisch verlaagt.

4. Resultaten

De analyse werd uitgevoerd voor koeltijden variërend van 0,1 seconden tot $10^5$ seconden.

• Invloed van de GLS-update op de data:

  • De aanpassing van de IFY-waarden zelf was gematigd voor ENDF/B-VIII.0 en JENDL-5, en verwaarloosbaar voor JEFF-3.3.
  • De grootste verandering was de reductie van de onzekerheid (variatie) door de introductie van correlaties.

• Onzekerheid in Vervalwarmte:

  • Zonder covarianties (Set-A): De ongecorreleerde opbrengstdata leverde een constante onzekerheid van ongeveer 4% op, wat dominant was voor koeltijden langer dan 100 seconden.
  • Met covarianties (Set-B): De bijdrage van de splijtingsopbrengst aan de totale onzekerheid nam sterk af.
    • Bij 0,1 s koeltijd: De totale onzekerheid daalde naar ~10% voor ENDF/B-VIII.0 en JEFF-3.3, en ~5% voor JENDL-5.
    • Bij $10^5$ s koeltijd: De onzekerheid daalde naar ongeveer 1% voor alle drie de bibliotheken.
  • Dominante bron: Na de introductie van covarianties werd de vervalenergie de belangrijkste bron van onzekerheid in de berekende vervalwarmte, in plaats van de splijtingsopbrengst.

• Vergelijking van bibliotheken:

  • Voor koeltijden langer dan 1000 s waren de resultaten van alle drie de bibliotheken in goede overeenstemming met experimentele data (CoNDERC database).
  • Voor korte koeltijden (< 1000 s) waren er discrepanties. JEFF-3.3 onderschatte de elektromagnetische vervalwarmte significant tussen 5 en 50 seconden.
  • JENDL-5 leverde over het algemeen kleinere onzekerheden op voor korte koeltijden omdat het vervalenergie-onzekerheden bevatte voor meer kortlevende radionucliden.

• Sensitiviteitscoëfficiënten:

  • De aanpassing van de IFY's had weinig invloed op de bijdrage van belangrijke splijtingsproducten voor ENDF/B-VIII.0 en JEFF-3.3.
  • Voor JENDL-5 waren er echter merkbare veranderingen in de bijdragen en sensitiviteiten van specifieke nucliden (zoals $^{97}$Y, $^{135}$Te, $^{137,138}$I en $^{89}$Br) na de GLS-update.

5. Betekenis en Conclusie

Dit werk biedt een robuuste methodologie om covariantiematrices voor splijtingsopbrengsten te genereren waar deze in bestaande databanken ontbreken.

• Technische impact: Het toont aan dat het negeren van correlaties leidt tot een overschatting van de onzekerheid in vervalwarmte. Door covarianties in te bouwen, kan de onzekerheid met een factor 2 tot 4 worden verlaagd (afhankelijk van de koeltijd).
• Praktische toepassing: De gegenereerde data en matrices zijn direct toepasbaar in veiligheidsanalyses en reactorontwerp, waardoor nauwkeurigere en betrouwbaardere berekeningen van vervalwarmte mogelijk zijn.
• Toekomst: De auteurs bevelen aan dat deze bijgewerkte data en covariantiematrices worden geïntegreerd in toekomstige edities van nucleaire databanken om de nauwkeurigheid van nucleaire veiligheidsanalyses wereldwijd te verbeteren.