An iterative fission matrix scheme for calculating steady-state power and critical control rod position in a TRIGA reactor

Tyler J. Topham, Adam Rau, William J. Walters

Research output: Contribution to journalArticle

Abstract

Although Monte Carlo simulations are the most accurate representation of nuclear reactor cores, they can be computationally expensive. This computational expense can prohibit the use of Monte Carlo in situations that require sequences of neutronics calculations, such as accounting for temperature feedback. In the present work, we demonstrate a fission matrix methodology to calculate the fission source in the Penn State Breazeale Reactor while accounting for control rod movement and fuel temperature feedback. A fission matrix A, represents the rate of fission neutron production in a reactor core, where the aij entry of the matrix is equal to the neutrons produced in cell i per neutron born from cell j. The fission source distribution and multiplication factor keff are the principle eigenvector and principle eigenvalue, respectively, of the fission matrix. Given a varied core temperature distribution, a new fission matrix can be obtained by interpolating from a database of pre-calculated fission matrices. This approach shows promise in accurately and quickly calculating the reactivity and fission source distributions for a system with an arbitrary temperature distribution. In the present work, the largest resulting relative error of the calculated power distribution per pin was 1.52% compared to a Serpent criticality simulation, and all calculated eigenvalues were within 50 pcm, which corresponds to 0.1–0.4 cm of rod insertion depending on whether rods are inserted in the center or edge of the core. In contrast, the critical rod position varies by 13.49 cm for the cases considered (0–650 kW). The present database is composed of data from 78 Monte Carlo calculations, but enables subsequent neutronics calculations to be done in a fraction of a second. Future work would be to implement correction factors for neighboring temperature differences in three-dimensional models, incorporate xenon feedback, and investigate ways to decrease database size.

Original languageEnglish (US)
Article number106984
JournalAnnals of Nuclear Energy
Volume135
DOIs
StatePublished - Jan 1 2020

Fingerprint

Control rods
Neutrons
Reactor cores
Feedback
Nuclear fuel accounting
Temperature distribution
Xenon
Eigenvalues and eigenfunctions
Temperature

All Science Journal Classification (ASJC) codes

  • Nuclear Energy and Engineering

Cite this

@article{d29b5bb0400a4a16a1414bc6b2c491f3,
title = "An iterative fission matrix scheme for calculating steady-state power and critical control rod position in a TRIGA reactor",
abstract = "Although Monte Carlo simulations are the most accurate representation of nuclear reactor cores, they can be computationally expensive. This computational expense can prohibit the use of Monte Carlo in situations that require sequences of neutronics calculations, such as accounting for temperature feedback. In the present work, we demonstrate a fission matrix methodology to calculate the fission source in the Penn State Breazeale Reactor while accounting for control rod movement and fuel temperature feedback. A fission matrix A, represents the rate of fission neutron production in a reactor core, where the aij entry of the matrix is equal to the neutrons produced in cell i per neutron born from cell j. The fission source distribution and multiplication factor keff are the principle eigenvector and principle eigenvalue, respectively, of the fission matrix. Given a varied core temperature distribution, a new fission matrix can be obtained by interpolating from a database of pre-calculated fission matrices. This approach shows promise in accurately and quickly calculating the reactivity and fission source distributions for a system with an arbitrary temperature distribution. In the present work, the largest resulting relative error of the calculated power distribution per pin was 1.52{\%} compared to a Serpent criticality simulation, and all calculated eigenvalues were within 50 pcm, which corresponds to 0.1–0.4 cm of rod insertion depending on whether rods are inserted in the center or edge of the core. In contrast, the critical rod position varies by 13.49 cm for the cases considered (0–650 kW). The present database is composed of data from 78 Monte Carlo calculations, but enables subsequent neutronics calculations to be done in a fraction of a second. Future work would be to implement correction factors for neighboring temperature differences in three-dimensional models, incorporate xenon feedback, and investigate ways to decrease database size.",
author = "Topham, {Tyler J.} and Adam Rau and Walters, {William J.}",
year = "2020",
month = "1",
day = "1",
doi = "10.1016/j.anucene.2019.106984",
language = "English (US)",
volume = "135",
journal = "Annals of Nuclear Energy",
issn = "0306-4549",
publisher = "Elsevier Limited",

}

An iterative fission matrix scheme for calculating steady-state power and critical control rod position in a TRIGA reactor. / Topham, Tyler J.; Rau, Adam; Walters, William J.

In: Annals of Nuclear Energy, Vol. 135, 106984, 01.01.2020.

Research output: Contribution to journalArticle

TY - JOUR

T1 - An iterative fission matrix scheme for calculating steady-state power and critical control rod position in a TRIGA reactor

AU - Topham, Tyler J.

AU - Rau, Adam

AU - Walters, William J.

PY - 2020/1/1

Y1 - 2020/1/1

N2 - Although Monte Carlo simulations are the most accurate representation of nuclear reactor cores, they can be computationally expensive. This computational expense can prohibit the use of Monte Carlo in situations that require sequences of neutronics calculations, such as accounting for temperature feedback. In the present work, we demonstrate a fission matrix methodology to calculate the fission source in the Penn State Breazeale Reactor while accounting for control rod movement and fuel temperature feedback. A fission matrix A, represents the rate of fission neutron production in a reactor core, where the aij entry of the matrix is equal to the neutrons produced in cell i per neutron born from cell j. The fission source distribution and multiplication factor keff are the principle eigenvector and principle eigenvalue, respectively, of the fission matrix. Given a varied core temperature distribution, a new fission matrix can be obtained by interpolating from a database of pre-calculated fission matrices. This approach shows promise in accurately and quickly calculating the reactivity and fission source distributions for a system with an arbitrary temperature distribution. In the present work, the largest resulting relative error of the calculated power distribution per pin was 1.52% compared to a Serpent criticality simulation, and all calculated eigenvalues were within 50 pcm, which corresponds to 0.1–0.4 cm of rod insertion depending on whether rods are inserted in the center or edge of the core. In contrast, the critical rod position varies by 13.49 cm for the cases considered (0–650 kW). The present database is composed of data from 78 Monte Carlo calculations, but enables subsequent neutronics calculations to be done in a fraction of a second. Future work would be to implement correction factors for neighboring temperature differences in three-dimensional models, incorporate xenon feedback, and investigate ways to decrease database size.

AB - Although Monte Carlo simulations are the most accurate representation of nuclear reactor cores, they can be computationally expensive. This computational expense can prohibit the use of Monte Carlo in situations that require sequences of neutronics calculations, such as accounting for temperature feedback. In the present work, we demonstrate a fission matrix methodology to calculate the fission source in the Penn State Breazeale Reactor while accounting for control rod movement and fuel temperature feedback. A fission matrix A, represents the rate of fission neutron production in a reactor core, where the aij entry of the matrix is equal to the neutrons produced in cell i per neutron born from cell j. The fission source distribution and multiplication factor keff are the principle eigenvector and principle eigenvalue, respectively, of the fission matrix. Given a varied core temperature distribution, a new fission matrix can be obtained by interpolating from a database of pre-calculated fission matrices. This approach shows promise in accurately and quickly calculating the reactivity and fission source distributions for a system with an arbitrary temperature distribution. In the present work, the largest resulting relative error of the calculated power distribution per pin was 1.52% compared to a Serpent criticality simulation, and all calculated eigenvalues were within 50 pcm, which corresponds to 0.1–0.4 cm of rod insertion depending on whether rods are inserted in the center or edge of the core. In contrast, the critical rod position varies by 13.49 cm for the cases considered (0–650 kW). The present database is composed of data from 78 Monte Carlo calculations, but enables subsequent neutronics calculations to be done in a fraction of a second. Future work would be to implement correction factors for neighboring temperature differences in three-dimensional models, incorporate xenon feedback, and investigate ways to decrease database size.

UR - http://www.scopus.com/inward/record.url?scp=85070687972&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85070687972&partnerID=8YFLogxK

U2 - 10.1016/j.anucene.2019.106984

DO - 10.1016/j.anucene.2019.106984

M3 - Article

AN - SCOPUS:85070687972

VL - 135

JO - Annals of Nuclear Energy

JF - Annals of Nuclear Energy

SN - 0306-4549

M1 - 106984

ER -