Optimizing the RMC Code Using the Decay Chain Method for Large-Scale Decay Calculations
Decay calculations play an important role in reactor physics simulations. In particular, the isotopic decay of the burned fuel during refueling is important for predicting the startup reactivity of the following burn-up cycle. In addition, there is a growing interest in high-fidelity simulations whe...
| Main Authors: | , , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Frontiers Media S.A.
2021-07-01
|
| Series: | Frontiers in Energy Research |
| Subjects: | |
| Online Access: | https://www.frontiersin.org/articles/10.3389/fenrg.2021.676889/full |
| _version_ | 1818865517138542592 |
|---|---|
| author | Hao Li Ganglin Yu Shanfang Huang Kan Wang |
| author_facet | Hao Li Ganglin Yu Shanfang Huang Kan Wang |
| author_sort | Hao Li |
| collection | DOAJ |
| description | Decay calculations play an important role in reactor physics simulations. In particular, the isotopic decay of the burned fuel during refueling is important for predicting the startup reactivity of the following burn-up cycle. In addition, there is a growing interest in high-fidelity simulations where the mesh in the burn-up region can involve millions of regions. However, existing models repeatedly solve the same Bateman equations for each region, which is a waste of calculational resources. RMC is a Monte Carlo neutron transport code developed for advanced reactor physics analysis including criticality calculations and burn-up calculations. This paper presents a decay calculation method named the Decay Chain Method (DCM) to optimize the RMC code for large-scale decay calculations. Unlike traditional methods, the Decay Chain Method solves the Bateman equations one decay chain at a time rather than one region at a time. The decay calculation in the burn-up mode then treats the decay steps as zero power burn-up steps with some optimized calculational methods to further reduce the calculational time. These methods were evaluated for a single pin example and for a Virtual Environment for Reactor Applications (VERA) full-core example. The calculations for the single pin example verify the accuracy of the decay step treatment in the burn-up mode and show the improved efficiency. The single pin is divided into 1–1,000,000 decay regions to study the efficiency differences between the Transmutation Trajectory Analysis (TTA) and DCM methods. Both methods have a linear complexity with respect to the number of regions but DCM costs just one-sixtieth of the TTA time. In the simplified VERA full core example, the DCM method reduces the decay calculation time to 0.32 min from 75.26 min while the accuracy remains unchanged. |
| first_indexed | 2024-12-19T10:48:48Z |
| format | Article |
| id | doaj.art-67142489d4f7445bba191efe021291ab |
| institution | Directory Open Access Journal |
| issn | 2296-598X |
| language | English |
| last_indexed | 2024-12-19T10:48:48Z |
| publishDate | 2021-07-01 |
| publisher | Frontiers Media S.A. |
| record_format | Article |
| series | Frontiers in Energy Research |
| spelling | doaj.art-67142489d4f7445bba191efe021291ab2022-12-21T20:25:08ZengFrontiers Media S.A.Frontiers in Energy Research2296-598X2021-07-01910.3389/fenrg.2021.676889676889Optimizing the RMC Code Using the Decay Chain Method for Large-Scale Decay CalculationsHao LiGanglin YuShanfang HuangKan WangDecay calculations play an important role in reactor physics simulations. In particular, the isotopic decay of the burned fuel during refueling is important for predicting the startup reactivity of the following burn-up cycle. In addition, there is a growing interest in high-fidelity simulations where the mesh in the burn-up region can involve millions of regions. However, existing models repeatedly solve the same Bateman equations for each region, which is a waste of calculational resources. RMC is a Monte Carlo neutron transport code developed for advanced reactor physics analysis including criticality calculations and burn-up calculations. This paper presents a decay calculation method named the Decay Chain Method (DCM) to optimize the RMC code for large-scale decay calculations. Unlike traditional methods, the Decay Chain Method solves the Bateman equations one decay chain at a time rather than one region at a time. The decay calculation in the burn-up mode then treats the decay steps as zero power burn-up steps with some optimized calculational methods to further reduce the calculational time. These methods were evaluated for a single pin example and for a Virtual Environment for Reactor Applications (VERA) full-core example. The calculations for the single pin example verify the accuracy of the decay step treatment in the burn-up mode and show the improved efficiency. The single pin is divided into 1–1,000,000 decay regions to study the efficiency differences between the Transmutation Trajectory Analysis (TTA) and DCM methods. Both methods have a linear complexity with respect to the number of regions but DCM costs just one-sixtieth of the TTA time. In the simplified VERA full core example, the DCM method reduces the decay calculation time to 0.32 min from 75.26 min while the accuracy remains unchanged.https://www.frontiersin.org/articles/10.3389/fenrg.2021.676889/fulldecay chain methodtransmutation trajectory analysislarge-scale burnup calculationvirtual environment for reactor applications (VERA)reactor Monte Carlo code (RMC) |
| spellingShingle | Hao Li Ganglin Yu Shanfang Huang Kan Wang Optimizing the RMC Code Using the Decay Chain Method for Large-Scale Decay Calculations Frontiers in Energy Research decay chain method transmutation trajectory analysis large-scale burnup calculation virtual environment for reactor applications (VERA) reactor Monte Carlo code (RMC) |
| title | Optimizing the RMC Code Using the Decay Chain Method for Large-Scale Decay Calculations |
| title_full | Optimizing the RMC Code Using the Decay Chain Method for Large-Scale Decay Calculations |
| title_fullStr | Optimizing the RMC Code Using the Decay Chain Method for Large-Scale Decay Calculations |
| title_full_unstemmed | Optimizing the RMC Code Using the Decay Chain Method for Large-Scale Decay Calculations |
| title_short | Optimizing the RMC Code Using the Decay Chain Method for Large-Scale Decay Calculations |
| title_sort | optimizing the rmc code using the decay chain method for large scale decay calculations |
| topic | decay chain method transmutation trajectory analysis large-scale burnup calculation virtual environment for reactor applications (VERA) reactor Monte Carlo code (RMC) |
| url | https://www.frontiersin.org/articles/10.3389/fenrg.2021.676889/full |
| work_keys_str_mv | AT haoli optimizingthermccodeusingthedecaychainmethodforlargescaledecaycalculations AT ganglinyu optimizingthermccodeusingthedecaychainmethodforlargescaledecaycalculations AT shanfanghuang optimizingthermccodeusingthedecaychainmethodforlargescaledecaycalculations AT kanwang optimizingthermccodeusingthedecaychainmethodforlargescaledecaycalculations |