Heat Source Characterization In A TREAT Fuel Particle Using Coupled Neutronics Binary Collision MonteCarlo Calculations
Abstract
This work presents a multiphysics, multiscale approach to modeling the Transient Test Reactor (TREAT) currently prepared for restart at the Idaho National Laboratory. TREAT fuel is made up of microscopic fuel grains (r ˜ 20µm) dispersed in a graphite matrix. The novelty of this work is in coupling a binary collision MonteCarlo (BCMC) model to the Finite Element based code Moose for solving a microsopic heatconduction problem whose driving source is provided by the BCMC model tracking fission fragment energy deposition. This microscopic model is driven by a transient, engineering scale neutronics model coupled to an adiabatic heating model. The macroscopic model provides local power densities and neutron energy spectra to the microscpic model. Currently, no feedback from the microscopic to the macroscopic model is considered. TREAT transient 15 is used to exemplify the capabilities of the multiphysics, multiscale model, and it is found that the average fuel grain temperature differs from the average graphite temperature by 80 K despite the lowpower transient. The large temperature difference has strong implications on the Doppler feedback a potential LEU TREAT core would see, and it underpins the need for multiphysics, multiscale modeling of a TREAT LEU core.
 Authors:
 Publication Date:
 Research Org.:
 Idaho National Lab. (INL), Idaho Falls, ID (United States)
 Sponsoring Org.:
 USDOE Office of Nuclear Energy (NE)
 OSTI Identifier:
 1369423
 Report Number(s):
 INL/CON1640179
 DOE Contract Number:
 DEAC0705ID14517
 Resource Type:
 Conference
 Resource Relation:
 Conference: M&C 2017 International Conference on Mathematics & Computational Methods Applied to Nuclear Science & Engineering, Jeju, Republic of Korea (South Korea), April 16–20, 2017
 Country of Publication:
 United States
 Language:
 English
 Subject:
 97 MATHEMATICS AND COMPUTING; 36 MATERIALS SCIENCE; Binary Collision MonteCarlo; Multiphysics; TREAT
Citation Formats
Schunert, Sebastian, Schwen, Daniel, Ghassemi, Pedram, Baker, Benjamin, Zabriskie, Adam, Ortensi, Javier, Wang, Yaqi, Gleicher, Frederick, DeHart, Mark, and Martineau, Richard. Heat Source Characterization In A TREAT Fuel Particle Using Coupled Neutronics Binary Collision MonteCarlo Calculations. United States: N. p., 2017.
Web.
Schunert, Sebastian, Schwen, Daniel, Ghassemi, Pedram, Baker, Benjamin, Zabriskie, Adam, Ortensi, Javier, Wang, Yaqi, Gleicher, Frederick, DeHart, Mark, & Martineau, Richard. Heat Source Characterization In A TREAT Fuel Particle Using Coupled Neutronics Binary Collision MonteCarlo Calculations. United States.
Schunert, Sebastian, Schwen, Daniel, Ghassemi, Pedram, Baker, Benjamin, Zabriskie, Adam, Ortensi, Javier, Wang, Yaqi, Gleicher, Frederick, DeHart, Mark, and Martineau, Richard. Sat .
"Heat Source Characterization In A TREAT Fuel Particle Using Coupled Neutronics Binary Collision MonteCarlo Calculations". United States.
doi:. https://www.osti.gov/servlets/purl/1369423.
@article{osti_1369423,
title = {Heat Source Characterization In A TREAT Fuel Particle Using Coupled Neutronics Binary Collision MonteCarlo Calculations},
author = {Schunert, Sebastian and Schwen, Daniel and Ghassemi, Pedram and Baker, Benjamin and Zabriskie, Adam and Ortensi, Javier and Wang, Yaqi and Gleicher, Frederick and DeHart, Mark and Martineau, Richard},
abstractNote = {This work presents a multiphysics, multiscale approach to modeling the Transient Test Reactor (TREAT) currently prepared for restart at the Idaho National Laboratory. TREAT fuel is made up of microscopic fuel grains (r ˜ 20µm) dispersed in a graphite matrix. The novelty of this work is in coupling a binary collision MonteCarlo (BCMC) model to the Finite Element based code Moose for solving a microsopic heatconduction problem whose driving source is provided by the BCMC model tracking fission fragment energy deposition. This microscopic model is driven by a transient, engineering scale neutronics model coupled to an adiabatic heating model. The macroscopic model provides local power densities and neutron energy spectra to the microscpic model. Currently, no feedback from the microscopic to the macroscopic model is considered. TREAT transient 15 is used to exemplify the capabilities of the multiphysics, multiscale model, and it is found that the average fuel grain temperature differs from the average graphite temperature by 80 K despite the lowpower transient. The large temperature difference has strong implications on the Doppler feedback a potential LEU TREAT core would see, and it underpins the need for multiphysics, multiscale modeling of a TREAT LEU core.},
doi = {},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Sat Apr 01 00:00:00 EDT 2017},
month = {Sat Apr 01 00:00:00 EDT 2017}
}

At the present time a Monte Carlo transport computer code is being designed and implemented at Lawrence Livermore National Laboratory to include the transport of: neutrons, photons, electrons and light charged particles as well as the coupling between all species of particles, e.g., photon induced electron emission. Since this code is being designed to handle all particles this approach is called the ''All Particle Method''. The code is designed as a test bed code to include as many different methods as possible (e.g., electron single or multiple scattering) and will be data driven to minimize the number of methods andmore »

A Comparison of Monte Carlo Particle Transport Algorithms for an Interior Source Binary Stochastic Medium Benchmark Suite
Particle transport through binary stochastic mixtures has received considerable research attention in the last two decades. Zimmerman and Adams proposed a Monte Carlo algorithm (Algorithm A) that solves the LevermorePomraning equations and another Monte Carlo algorithm (Algorithm B) that should be more accurate as a result of improved local material realization modeling. Zimmerman and Adams numerically confirmed these aspects of the Monte Carlo algorithms by comparing the reflection and transmission values computed using these algorithms to a standard suite of planar geometry binary stochastic mixture benchmark transport solutions. The benchmark transport problems are driven by an isotropic angular flux incidentmore » 
Coupled proton/neutron transport calculations using the S sub N and Monte Carlo methods
Coupled charged/neutral article transport calculations are most often carried out using the Monte Carol technique. For example, the ITS, EGS, and MCNP (Version 4) codes are used extensively for electron/photon transport calculations while HETC models the transport of protons, neutrons and heavy ions. In recent years there has been considerable progress in deterministic models of electron transport, and many of these models are applicable to protons. However, even with these new models (and the well established models for neutron transport) deterministic coupled neutron/proton transport calculations have not been feasible for most problems of interest, due to a lack of coupledmore »