DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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 » it achieves accuracy close to that of the full, high-fidelity depletion chain with respect to the studied responses.« less

Authors:
 [1]; ORCiD logo [2];  [3];  [4]
  1. Argonne National Lab. (ANL), Argonne, IL (United States)
  2. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
  3. Georgia Inst. of Technology, Atlanta, GA (United States)
  4. 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


A simple implicit coupling scheme for Monte Carlo neutronics and isotopic depletion
journal, June 2020


The stochastic implicit Euler method – A stable coupling scheme for Monte Carlo burnup calculations
journal, October 2013


Fourth- and sixth-order commutator-free Magnus integrators for linear and non-linear dynamical systems
journal, December 2006


Initial MCNP6 Release Overview
journal, December 2012

  • Goorley, T.; James, M.; Booth, T.
  • Nuclear Technology, Vol. 180, Issue 3
  • DOI: 10.13182/NT11-135

Commutator-free Lie group methods
journal, April 2003


RMC – A Monte Carlo code for reactor core analysis
journal, August 2015


Nuclide depletion capabilities in the Shift Monte Carlo code
journal, April 2018


Extended development of a Monte Carlo code OpenMC for fuel cycle simulation of molten salt reactor
journal, January 2020


FISPACT-II: An Advanced Simulation System for Activation, Transmutation and Material Modelling
journal, January 2017


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
  • DOI: 10.1115/ICONE25-66584

Modeling and simulation functional needs for molten salt reactor licensing
journal, December 2019


Substep methods for burnup calculations with Bateman solutions
journal, November 2011


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
  • DOI: 10.13182/NSE10-99

ENDF/B-VII.1 Nuclear Data for Science and Technology: Cross Sections, Covariances, Fission Product Yields and Decay Data
journal, December 2011


Testing actinide fission yield treatment in CINDER90 for use in MCNP6 burnup calculations
journal, November 2015


Windowed multipole for cross section Doppler broadening
journal, February 2016


Stochastic semi-implicit substep method for coupled depletion Monte-Carlo codes
journal, June 2016


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
  • DOI: 10.13182/NT08-A4003

Rational Approximations to the Matrix Exponential in Burnup Calculations
journal, October 2011

  • Pusa, Maria
  • Nuclear Science and Engineering, Vol. 169, Issue 2
  • DOI: 10.13182/NSE10-81

Numerically stable Monte Carlo-burnup-thermal hydraulic coupling schemes
journal, January 2014


Development of a versatile depletion code AMAC
journal, August 2020


Computing the Matrix Exponential in Burnup Calculations
journal, February 2010

  • Pusa, Maria; Leppänen, Jaakko
  • Nuclear Science and Engineering, Vol. 164, Issue 2
  • DOI: 10.13182/NSE09-14

Development and verification of code IMPC-Depletion for nuclide depletion calculation
journal, July 2020


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
  • DOI: 10.1063/1.2720459

Investigation of Inconsistent ENDF/B-VII.1 Independent and Cumulative Fission Product Yields with Proposed Revisions
journal, January 2015


The Magnus expansion and some of its applications
journal, January 2009


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
  • DOI: 10.13182/NT11-3

Higher order methods for burnup calculations with Bateman solutions
journal, September 2011


MC21 v.6.0 – A continuous-energy Monte Carlo particle transport code with integrated reactor feedback capabilities
journal, August 2015


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
  • DOI: 10.1007/s41365-019-0560-z

The Serpent Monte Carlo code: Status, development and applications in 2013
journal, August 2015


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
  • DOI: 10.13182/NSE11-36

Development and verification of LOOP: A Linkage of ORIGEN2.2 and OpenMC
journal, January 2017


OpenMC: A state-of-the-art Monte Carlo code for research and development
journal, August 2015


A Stochastic Approximation Method
journal, September 1951

  • Robbins, Herbert; Monro, Sutton
  • The Annals of Mathematical Statistics, Vol. 22, Issue 3
  • DOI: 10.1214/aoms/1177729586