Modeling Decay Heat with a Simplified Depletion Chain in OpenMC

In dit werk wordt de CASL-verarmingsketen voor OpenMC gemodificeerd door het toevoegen van pseudo-nucliden en "delay nuclides" om de nauwkeurigheid van de schattingen voor vervalwarmte aanzienlijk te verbeteren terwijl de rekenefficiëntie behouden blijft.

De kern

Hoe je een complexe kernreactor-berekening simpeler en slimmer maakt: Een verhaal over "Pseudo-atomen" en "Vertraging"

Stel je voor dat een kernreactor een enorme, complexe stad is. In deze stad wonen duizenden verschillende soorten "inwoners" (atoomsoorten of nucliden). Sommige inwoners zijn heel actief, anderen slapen, en weer anderen veranderen langzaam in iets anders.

Wanneer de reactor stopt met werken (de "schakelaar" wordt uitgezet), stoppen de nieuwe inwoners niet met produceren. De oude, onstabiele inwoners blijven nog langzaam veranderen en geven daarbij warmte af. Dit noemen we vervalwarmte (decay heat). Zelfs als de reactor uit is, moet deze warmte worden afgevoerd, anders kan het net als in Fukushima tot een ramp leiden.

Het Probleem: De Te Grote Lijst

Om te weten hoeveel warmte er nog vrijkomt, moeten wetenschappers een lijst bijhouden van al deze inwoners en hun veranderingen.

• De "Grote Lijst" (ENDF): Dit is de perfecte, maar enorme lijst met meer dan 3.800 inwoners. Het is accuraat, maar het is alsof je probeert een heel land te besturen met een pen en papier. Het kost enorm veel tijd en rekenkracht.
• De "Korte Lijst" (CASL): Om het sneller te maken, hebben wetenschappers een vereenvoudigde lijst gemaakt met slechts 228 inwoners. Dit werkt heel goed om te weten hoe de reactor draait, maar het is een ramp voor het voorspellen van de warmte na het uitschakelen. Het is alsof je alleen de burgemeesters en politieagenten op je lijst zet, maar vergeet de duizenden burgers die toch nog warmte produceren. De berekende warmte is dan veel te laag.

De Oplossing 1: De "Pseudo-Inwoner" (Pseudo-nuclides)

De auteurs van dit onderzoek (Tanmay Gupta en Benoit Forget) dachten: "Laten we de korte lijst niet volledig vervangen, maar laten we hem slim maken."

Ze bedachten een truc: Pseudo-inwoners.
In plaats van elke enkele inwoner apart te volgen, groepeerden ze honderden kleine inwoners in één "fictieve" inwoner.

• De Analogie: Stel je voor dat je in een stad 100 verschillende soorten bloemen hebt die allemaal dezelfde hoeveelheid water nodig hebben. In plaats van 100 aparte regelaars te hebben, maak je één "Super-regelaar" die namens al die 100 bloemen werkt.
• In de computercode creëerden ze 10 van deze "Super-inwoners" (pseudo-nucliden). Elke Super-inwoner vertegenwoordigt een groep echte atomen met een vergelijkbare snelheid van verandering.
• Het Resultaat: Door deze 10 extra figuren toe te voegen aan de korte lijst, werd de warmteberekening plotseling veel nauwkeuriger, zonder dat ze de hele lange lijst hoefden te gebruiken. Het was alsof ze de stad weer goed konden besturen met een klein team.

De Oplossing 2: De "Vertraging" (Delay Nuclides)

Maar er was nog een probleem.
Stel je voor dat een inwoner (A) verandert in inwoner (B), en die verandert weer in inwoner (C). In de echte wereld duurt dit even. In de computercode met de "Super-inwoners" deden ze dit echter te snel: ze lieten A direct veranderen in C, alsof er geen tijd tussen zat.

• Het Effect: Als je de reactor uitschakelt, blijft de echte inwoner A nog even doorgaan om B en C te maken. Maar de computer dacht: "A is weg, dus C is ook weg." Dit leidde tot fouten in de warmteberekening net na het uitschakelen.
• De Oplossing: Ze introduceerden "Vertraging-inwoners".
  • In plaats van dat A direct C maakt, maakt A eerst een "Vertraging-inwoner" (een tussenstap). Deze vertraging-inwoner wacht even (met de juiste snelheid) voordat hij verandert in C.
  • De Analogie: Het is alsof je een pakketje niet direct bezorgt, maar eerst in een postkantoor legt dat even wacht voordat het wordt verzonden. Zo blijft de stroom van "warmte-producerende" inwoners kloppen, zelfs als de bron stopt.

Wat is het eindresultaat?

Door deze twee slimme trucjes te combineren (Super-inwoners + Vertraging-inwoners) hebben de onderzoekers een systeem gemaakt dat:

1. Snel is: Het gebruikt nog steeds de korte lijst, dus het rekent 50% sneller dan de enorme lijst.
2. Accuraat is: De voorspelling van de warmte na het uitschakelen is nu bijna perfect (minder dan 5% fout, en zelfs minder dan 0,3% voor de meest precieze versie).

Samenvattend:
Ze hebben de "korte lijst" van de reactor niet zomaar geknipt, maar hem verrijkt met slimme vertegenwoordigers en tijdbalken. Hierdoor kunnen ingenieurs veiligere en snellere berekeningen maken voor kernreactoren, zonder de complexiteit van duizenden losse atomen hoeven te volgen. Het is een mooi voorbeeld van hoe je een complex probleem oplost door niet alles apart te tellen, maar de juiste groepen te maken.

Technische samenvatting

Titel: Modellering van Vervalwarmte met een Vereenvoudigde Uitputtingsketen in OpenMC

1. Probleemstelling

In nucleaire reactoren is vervalwarmte (decay heat) een kritieke parameter voor veiligheidsanalyses, vooral na het afsluiten van de reactor. Hoewel vervalwarmte tijdens bedrijf ongeveer 6-7% van het totale vermogen uitmaakt, wordt het de dominante warmtebron na shutdown. Accurate modellering is essentieel voor het ontwerp van koelsystemen om incidenten zoals Fukushima te voorkomen.

OpenMC, een open-source Monte Carlo-transportcode, kan uitputtingsberekeningen (depletion) uitvoeren om vervalwarmte te schatten. Hiervoor is een "depletion chain" (uitputtingsketen) nodig die de transmutatiepaden en vervalroutes van nucliden beschrijft.

• De Benchmark: De volledige ENDF/B-VII.1 keten bevat data voor meer dan 3.800 nucliden en dient als referentiewaarde, maar is computationally zeer duur (hoge geheugenvraag en rekentijd).
• De Uitdaging: De vereenvoudigde CASL-keten (ontwikkeld voor het Consortium for Advanced Simulation of Light Water Reactors) bevat slechts 228 nucliden. Hoewel deze keten zeer efficiënt is voor het berekenen van de effectieve vermenigvuldigingsfactor ($k_{eff}$) en reactiviteit tijdens bedrijf, underschat deze de vervalwarmte drastisch. Dit komt doordat de keten honderden specifieke nucliden mist die weliswaar een klein aandeel hebben in de reactiviteit, maar gezamenlijk een groot deel van de vervalwarmte leveren.

2. Methodologie

De auteurs hebben een methode ontwikkeld om de nauwkeurigheid van de CASL-keten te verbeteren zonder de rekenkosten significant te verhogen. Dit werd gedaan door twee concepten te introduceren: Pseudo-nucliden (PNs) en Vertraging-nucliden (Delay Nuclides).

A. Pseudo-nucliden (Pseudo-Nuclides - PNs)
Om de bijdrage van de honderden ontbrekende ENDF-exclusieve nucliden te vangen zonder ze individueel te tracken, worden PNs gebruikt.

• Groepering: De ontbrekende nucliden worden gegroepeerd op basis van hun vervalconstante ($\lambda$). De auteurs hebben 10 groepen gedefinieerd op een logaritmische schaal.
• Representatie: Elke groep wordt vertegenwoordigd door één fictief nuclide (PN) met een representatieve vervalconstante ($\lambda_{PN}$), berekend als het gemiddelde van de groep.
• Productie: In plaats van de nuclide-dichtheid te tracken, trackt OpenMC de som van de producten van de vervalenergie ($Q_i$) en de dichtheid ($N_i$) voor elke groep. De productiesnelheid van een PN wordt afgeleid uit de productieroutes in de volledige ENDF-keten, geschaald met de vervalenergie ($Q_i$).
• Voordeel: Dit benadert de totale vervalwarmte van honderden nucliden met slechts 10 extra nucliden.

B. Het Probleem met Directe Vervalroutes (Startup/Shutdown Fouten)
Bij het testen van de CASL+10PNs keten bleek dat er aanzienlijke fouten optreden tijdens de overgangsfasen (opstarten en afsluiten).

• Oorzaak: De implementatie van PNs veronderstelde vaak "instantane verval" (instantaneous decay). In de werkelijkheid hebben sommige vervalroutes tussen een CASL-voorouder en een PN een tussenliggende nuclide met een relatief lange halfwaardetijd. Door dit te negeren, wordt de productie van de PN te snel (tijdens opstarten) of stopt deze te vroeg (tijdens afsluiten), wat leidt tot overschatting of onderschatting van de warmte.

C. Introductie van Vertraging-nucliden (Delay Nuclides)
Om de timing van de productie te corrigeren, introduceerden de auteurs "delay nuclides".

• Functie: Een delay nuclide fungeert als een tijdelijke opslag in de uitputtingsketen die de structuur van de vervalroute nabootst. Het vertraagt de productie van de PN totdat de fysieke halfwaardetijd van de tussenliggende nuclide is verstreken.
• Implementatie: Voor vervalroutes met een maximale halfwaardetijd groter dan een bepaalde drempel (bijv. 1 minuut), wordt een delay nuclide toegevoegd. Dit nuclide verval met zijn eigen constante ($\lambda_{delay}$) naar de juiste PN.
• Optimalisatie: Voor de vereenvoudigde keten (CASL+10PNs) werden de delay nuclides ook gegroepeerd (10 delay nucliden) om het aantal te tracken nucliden laag te houden.

3. Belangrijkste Resultaten

De auteurs testten de methoden op twee modellen: een Pressurized Water Reactor (PWR) en een Sodium-cooled Fast Reactor (SFR).

• Verbetering Vervalwarmte:

  • De originele CASL-keten onderschatte de vervalwarmte met tot wel 70% na shutdown.
  • De toevoeging van 10 PNs (CASL+10PNs) bracht de fout tijdens bedrijf terug tot <5%, maar liet nog steeds pieken in fouten zien bij opstarten en afsluiten.
  • De toevoeging van delay nuclides (CASL+10PNs v2) elimineerde deze pieken. De fout in de tijd-geïntegreerde vervalwarmte bleef onder de 5% voor alle tijdsintervallen na shutdown.
  • De "All PNs" variant met delay nuclides (CASL+All PNs v2) bereikte een foutmarge van <0,3%, vergelijkbaar met de volledige ENDF-benchmark.

• Rekenprestaties:

  • Hoewel de volledige ENDF-keten de meest accurate is, is deze zeer duur.
  • De vereenvoudigde ketens (CASL+10PNs v2) waren ongeveer 50% sneller in totale rekentijd (inclusief transportoplossing) dan de ENDF-keten.
  • De tijd voor de uitputtingsoplosser (depletion solver) was 70% sneller dan bij de ENDF-keten.
  • De geheugeneisen bleven laag doordat het aantal te tracken nucliden beperkt bleef tot ongeveer 240 (228 CASL + 10 PNs + 10 delay).

• Generaliseerbaarheid:

  • De methode bleek ook effectief voor het SFR-model (snelle neutronenspectrum), wat aantoont dat de aanpak niet specifiek is voor thermische reactoren, maar algemeen toepasbaar is.

4. Bijdragen en Betekenis

De belangrijkste bijdragen van dit werk zijn:

1. Oplossing voor een fundamenteel tekort: Het oplossen van het probleem waarbij vereenvoudigde uitputtingsketens (zoals CASL) nuttig zijn voor reactiviteitsberekeningen, maar onbruikbaar voor veiligheidsanalyses van vervalwarmte.
2. Innovatieve methodiek: De combinatie van pseudo-nucliden (voor het samenvatten van energiebijdragen) en delay nuclides (voor het corrigeren van tijdsafhankelijke productiefouten) is een krachtige nieuwe aanpak.
3. Efficiëntie-Nauwkeurigheid Trade-off: Het demonstreert dat het mogelijk is om de nauwkeurigheid van een volledige keten (ENDF) te benaderen met slechts een fractie van de rekenkosten door slimme fysieke benaderingen toe te passen in de ketenstructuur.
4. Praktische Toepassing: De methode maakt het mogelijk om OpenMC te gebruiken voor gedetailleerde veiligheidsanalyses van vervalwarmte in complexe scenario's (zoals langdurige shutdowns) zonder de onhoudbare rekenlast van duizenden nucliden.

Conclusie:
Dit onderzoek biedt een robuuste, schaalbare oplossing voor het modelleren van vervalwarmte in OpenMC. Door de introductie van delay nuclides worden de tijdsafhankelijke fouten die inherent waren aan eerdere benaderingen met pseudo-nucliden grotendeels opgelost, waardoor vereenvoudigde ketens nu veilig en accuraat kunnen worden gebruikt voor zowel operationele als post-shutdown veiligheidsanalyses.