Depletion capabilities in the OpenMC Monte Carlo particle transport code
Abstract
A depletion solver has been implemented in OpenMC and is described herein. The depletion solver is implemented in Python and interfaces with OpenMC’s transport solver through a C++ application programming interface, which enables an in-memory transport-depletion coupling. Multiple integration methods for advancing in time have been implemented and exhibit tradeoffs in cost, accuracy, and memory use. For all time integration methods, evaluation of the matrix exponential is performed by using the incomplete partial fraction form of the Chebyshev rational approximation method. Simulations of a pressurized water reactor (PWR) pincell and a sodium-cooled fast reactor (SFR) assembly were carried out with OpenMC and Serpent. For both problems, the use of a high-fidelity depletion chain results in predictions of keff that agree within 20–30 pcm between OpenMC and Serpent. Predicted actinide concentrations were found to agree to a fraction of a percent, and most fission product concentrations were found to agree within 1%. Here, the few cases where larger differences were observed can be attributed either to differences in how the energy dependence of fission product yields is handled or deficiencies in the nuclear data used. OpenMC simulations of the PWR and SFR problems using a simplified 228-nuclide depletion chain demonstrate thatmore »
- Authors:
-
- Argonne National Lab. (ANL), Argonne, IL (United States)
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- Georgia Inst. of Technology, Atlanta, GA (United States)
- Tsinghua Univ., Beijing (China)
- Publication Date:
- Research Org.:
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Argonne National Lab. (ANL), Argonne, IL (United States)
- Sponsoring Org.:
- USDOE National Nuclear Security Administration (NNSA); USDOE Office of Science (SC)
- OSTI Identifier:
- 1726218
- Alternate Identifier(s):
- OSTI ID: 1776818
- Report Number(s):
- LA-UR-20-25299
Journal ID: ISSN 0306-4549; TRN: US2205088
- Grant/Contract Number:
- 89233218CNA000001; 17-SC-20-SC; AC02-06CH11357; AC07-05ID14517
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Annals of Nuclear Energy (Oxford)
- Additional Journal Information:
- Journal Name: Annals of Nuclear Energy (Oxford); Journal Volume: 152; Journal Issue: C; Journal ID: ISSN 0306-4549
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 73 NUCLEAR PHYSICS AND RADIATION PHYSICS; Monte Carlo; particle transport; depletion; transmutation
Citation Formats
Romano, Paul K., Josey, Colin J., Johnson, Andrew E., and Liang, Jingang. Depletion capabilities in the OpenMC Monte Carlo particle transport code. United States: N. p., 2020.
Web. doi:10.1016/j.anucene.2020.107989.
Romano, Paul K., Josey, Colin J., Johnson, Andrew E., & Liang, Jingang. Depletion capabilities in the OpenMC Monte Carlo particle transport code. United States. https://doi.org/10.1016/j.anucene.2020.107989
Romano, Paul K., Josey, Colin J., Johnson, Andrew E., and Liang, Jingang. Sat .
"Depletion capabilities in the OpenMC Monte Carlo particle transport code". United States. https://doi.org/10.1016/j.anucene.2020.107989. https://www.osti.gov/servlets/purl/1726218.
@article{osti_1726218,
title = {Depletion capabilities in the OpenMC Monte Carlo particle transport code},
author = {Romano, Paul K. and Josey, Colin J. and Johnson, Andrew E. and Liang, Jingang},
abstractNote = {A depletion solver has been implemented in OpenMC and is described herein. The depletion solver is implemented in Python and interfaces with OpenMC’s transport solver through a C++ application programming interface, which enables an in-memory transport-depletion coupling. Multiple integration methods for advancing in time have been implemented and exhibit tradeoffs in cost, accuracy, and memory use. For all time integration methods, evaluation of the matrix exponential is performed by using the incomplete partial fraction form of the Chebyshev rational approximation method. Simulations of a pressurized water reactor (PWR) pincell and a sodium-cooled fast reactor (SFR) assembly were carried out with OpenMC and Serpent. For both problems, the use of a high-fidelity depletion chain results in predictions of keff that agree within 20–30 pcm between OpenMC and Serpent. Predicted actinide concentrations were found to agree to a fraction of a percent, and most fission product concentrations were found to agree within 1%. Here, the few cases where larger differences were observed can be attributed either to differences in how the energy dependence of fission product yields is handled or deficiencies in the nuclear data used. OpenMC simulations of the PWR and SFR problems using a simplified 228-nuclide depletion chain demonstrate that it achieves accuracy close to that of the full, high-fidelity depletion chain with respect to the studied responses.},
doi = {10.1016/j.anucene.2020.107989},
journal = {Annals of Nuclear Energy (Oxford)},
number = C,
volume = 152,
place = {United States},
year = {Sat Nov 14 00:00:00 EST 2020},
month = {Sat Nov 14 00:00:00 EST 2020}
}
Works referenced in this record:
SciPy 1.0: fundamental algorithms for scientific computing in Python
journal, February 2020
- Virtanen, Pauli; Gommers, Ralf; Oliphant, Travis E.
- Nature Methods
A simple implicit coupling scheme for Monte Carlo neutronics and isotopic depletion
journal, June 2020
- Cosgrove, P.; Shwageraus, E.; Parks, G. T.
- Annals of Nuclear Energy, Vol. 141
The stochastic implicit Euler method – A stable coupling scheme for Monte Carlo burnup calculations
journal, October 2013
- Dufek, Jan; Kotlyar, Dan; Shwageraus, Eugene
- Annals of Nuclear Energy, Vol. 60
Fourth- and sixth-order commutator-free Magnus integrators for linear and non-linear dynamical systems
journal, December 2006
- Blanes, S.; Moan, P. C.
- Applied Numerical Mathematics, Vol. 56, Issue 12
ORIGEN2: A Versatile Computer Code for Calculating the Nuclide Compositions and Characteristics of Nuclear Materials
journal, September 1983
- Croff, Allen G.
- Nuclear Technology, Vol. 62, Issue 3
Initial MCNP6 Release Overview
journal, December 2012
- Goorley, T.; James, M.; Booth, T.
- Nuclear Technology, Vol. 180, Issue 3
Commutator-free Lie group methods
journal, April 2003
- Celledoni, Elena; Marthinsen, Arne; Owren, Brynjulf
- Future Generation Computer Systems, Vol. 19, Issue 3
RMC – A Monte Carlo code for reactor core analysis
journal, August 2015
- Wang, Kan; Li, Zeguang; She, Ding
- Annals of Nuclear Energy, Vol. 82
Nuclide depletion capabilities in the Shift Monte Carlo code
journal, April 2018
- Davidson, Gregory G.; Pandya, Tara M.; Johnson, Seth R.
- Annals of Nuclear Energy, Vol. 114
Extended development of a Monte Carlo code OpenMC for fuel cycle simulation of molten salt reactor
journal, January 2020
- Zhuang, Kun; Li, Ting; Zhang, Qian
- Progress in Nuclear Energy, Vol. 118
FISPACT-II: An Advanced Simulation System for Activation, Transmutation and Material Modelling
journal, January 2017
- Sublet, J. -Ch.; Eastwood, J. W.; Morgan, J. G.
- Nuclear Data Sheets, Vol. 139
Toward an Open-Source Neutronics Code for Circulating-Fuel Reactors
conference, October 2017
- de Troullioud de Lanversin, Julien; Glaser, Alexander; Göttsche, Malte
- 2017 25th International Conference on Nuclear Engineering, Volume 5: Advanced and Next Generation Reactors, Fusion Technology; Codes, Standards, Conformity Assessment, Licensing, and Regulatory Issues
Modeling and simulation functional needs for molten salt reactor licensing
journal, December 2019
- Betzler, Benjamin R.; Heidet, Florent; Feng, Bo
- Nuclear Engineering and Design, Vol. 355
Substep methods for burnup calculations with Bateman solutions
journal, November 2011
- Isotalo, A. E.; Aarnio, P. A.
- Annals of Nuclear Energy, Vol. 38, Issue 11
A Burnup-Dependent Isomeric Production Branching Ratio Treatment
journal, May 2012
- Haeck, W.; Cochet, B.; Aguiar, L.
- Nuclear Science and Engineering, Vol. 171, Issue 1
ENDF/B-VII.1 Nuclear Data for Science and Technology: Cross Sections, Covariances, Fission Product Yields and Decay Data
journal, December 2011
- Chadwick, M. B.; Herman, M.; Obložinský, P.
- Nuclear Data Sheets, Vol. 112, Issue 12
Testing actinide fission yield treatment in CINDER90 for use in MCNP6 burnup calculations
journal, November 2015
- Fensin, Michael Lorne; Umbel, Marissa
- Progress in Nuclear Energy, Vol. 85
Windowed multipole for cross section Doppler broadening
journal, February 2016
- Josey, C.; Ducru, P.; Forget, B.
- Journal of Computational Physics, Vol. 307
Stochastic semi-implicit substep method for coupled depletion Monte-Carlo codes
journal, June 2016
- Kotlyar, D.; Shwageraus, E.
- Annals of Nuclear Energy, Vol. 92
Improved Reaction Rate Tracking and Fission Product Yield Determinations for the Monte Carlo-Linked Depletion Capability in MCNPX
journal, October 2008
- Fensin, Michael L.; Hendricks, John S.; Anghaie, Samim
- Nuclear Technology, Vol. 164, Issue 1
Rational Approximations to the Matrix Exponential in Burnup Calculations
journal, October 2011
- Pusa, Maria
- Nuclear Science and Engineering, Vol. 169, Issue 2
Numerically stable Monte Carlo-burnup-thermal hydraulic coupling schemes
journal, January 2014
- Kotlyar, D.; Shwageraus, E.
- Annals of Nuclear Energy, Vol. 63
Development of a versatile depletion code AMAC
journal, August 2020
- Zhang, Binhang; Yuan, XianBao; Zhang, Yonghong
- Annals of Nuclear Energy, Vol. 143
Computing the Matrix Exponential in Burnup Calculations
journal, February 2010
- Pusa, Maria; Leppänen, Jaakko
- Nuclear Science and Engineering, Vol. 164, Issue 2
Development and verification of code IMPC-Depletion for nuclide depletion calculation
journal, July 2020
- Zhao, Zelong; Yang, Yongwei; Gao, Qingyu
- Nuclear Engineering and Design, Vol. 363
The MCNPX Monte Carlo Radiation Transport Code
conference, January 2007
- Waters, Laurie S.; McKinney, Gregg W.; Durkee, Joe W.
- HADRONIC SHOWER SIMULATION WORKSHOP, AIP Conference Proceedings
Investigation of Inconsistent ENDF/B-VII.1 Independent and Cumulative Fission Product Yields with Proposed Revisions
journal, January 2015
- Pigni, M. T.; Francis, M. W.; Gauld, I. C.
- Nuclear Data Sheets, Vol. 123
The Magnus expansion and some of its applications
journal, January 2009
- Blanes, S.; Casas, F.; Oteo, J. A.
- Physics Reports, Vol. 470, Issue 5-6
Isotopic Depletion and Decay Methods and Analysis Capabilities in SCALE
journal, May 2011
- Gauld, Ian C.; Radulescu, Georgeta; Ilas, Germina
- Nuclear Technology, Vol. 174, Issue 2
Higher order methods for burnup calculations with Bateman solutions
journal, September 2011
- Isotalo, A. E.; Aarnio, P. A.
- Annals of Nuclear Energy, Vol. 38, Issue 9
MC21 v.6.0 – A continuous-energy Monte Carlo particle transport code with integrated reactor feedback capabilities
journal, August 2015
- Griesheimer, D. P.; Gill, D. F.; Nease, B. R.
- Annals of Nuclear Energy, Vol. 82
Development and validation of depletion code system IMPC-Burnup for ADS
journal, February 2019
- Liu, Zhao-Qing; Zhao, Ze-Long; Yang, Yong-Wei
- Nuclear Science and Techniques, Vol. 30, Issue 3
The Serpent Monte Carlo code: Status, development and applications in 2013
journal, August 2015
- Leppänen, Jaakko; Pusa, Maria; Viitanen, Tuomas
- Annals of Nuclear Energy, Vol. 82
Explicit Treatment of Thermal Motion in Continuous-Energy Monte Carlo Tracking Routines
journal, June 2012
- Viitanen, Tuomas; Leppänen, Jaakko
- Nuclear Science and Engineering, Vol. 171, Issue 2
Development and verification of LOOP: A Linkage of ORIGEN2.2 and OpenMC
journal, January 2017
- Gul, Anas; Chaudri, K. S.; Khan, R.
- Annals of Nuclear Energy, Vol. 99
OpenMC: A state-of-the-art Monte Carlo code for research and development
journal, August 2015
- Romano, Paul K.; Horelik, Nicholas E.; Herman, Bryan R.
- Annals of Nuclear Energy, Vol. 82
A Stochastic Approximation Method
journal, September 1951
- Robbins, Herbert; Monro, Sutton
- The Annals of Mathematical Statistics, Vol. 22, Issue 3