skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: THE APPLICATION OF MAMMOTH FOR A DETAILED TIGHTLY COUPLED FUEL PIN SIMULATION WITH A STATION BLACKOUT

Abstract

Accurate calculation of desired quantities to predict fuel behavior requires the solution of interlinked equations representing different physics. Traditional fuels performance codes often rely on internal empirical models for the pin power density and a simplified boundary condition on the cladding edge. These simplifications are performed because of the difficulty of coupling applications or codes on differing domains and mapping the required data. To demonstrate an approach closer to first principles, the neutronics application Rattlesnake and the thermal hydraulics application RELAP-7 were coupled to the fuels performance application BISON under the master application MAMMOTH. A single fuel pin was modeled based on the dimensions of a Westinghouse 17x17 fuel rod. The simulation consisted of a depletion period of 1343 days, roughly equal to three full operating cycles, followed by a station blackout (SBO) event. The fuel rod was depleted for 1343 days for a near constant total power loading of 65.81 kW. After 1343 days the fission power was reduced to zero (simulating a reactor shut-down). Decay heat calculations provided the time-varying energy source after this time. For this problem, Rattlesnake, BISON, and RELAP-7 are coupled under MAMMOTH in a split operator approach. Each system solves its physics on amore » separate mesh and, for RELAP-7 and BISON, on only a subset of the full problem domain. Rattlesnake solves the neutronics over the whole domain that includes the fuel, cladding, gaps, water, and top and bottom rod holders. Here BISON is applied to the fuel and cladding with a 2D axi-symmetric domain, and RELAP-7 is applied to the flow of the circular outer water channel with a set of 1D flow equations. The mesh on the Rattlesnake side can either be 3D (for low order transport) or 2D (for diffusion). BISON has a matching ring structure mesh for the fuel so both the power density and local burn up are copied accurately from Rattlesnake. At each depletion time step, Rattlesnake calculates a power density, fission density rate, burn-up distribution and fast flux based on the current water density and fuel temperature. These are then mapped to the BISON mesh for a fuels performance solve. BISON calculates the fuel temperature and cladding surface temperature based upon the current power density and bulk fluid temperature. RELAP-7 then calculates the fluid temperature, water density fraction and water phase velocity based upon the cladding surface temperature. The fuel temperature and the fluid density are then passed back to Rattlesnake for another neutronics calculation. Six Picard or fixed-point style iterations are preformed in this manner to obtain consistent tightly coupled and stable results. For this paper a set of results from the detailed calculation are provided for both during depletion and the SBO event. We demonstrate that a detailed calculation closer to first principles can be done under MAMMOTH between different applications on differing domains.« less

Authors:
; ; ; ; ; ; ; ; ; ; ;
Publication Date:
Research Org.:
Idaho National Lab. (INL), Idaho Falls, ID (United States)
Sponsoring Org.:
USDOE Office of Nuclear Energy (NE)
OSTI Identifier:
1363897
Report Number(s):
INL/CON-16-37597
DOE Contract Number:  
DE-AC07-05ID14517
Resource Type:
Conference
Resource Relation:
Conference: Top Fuel 2016, Boise, Idaho, September 11–16, 2016
Country of Publication:
United States
Language:
English
Subject:
36 MATERIALS SCIENCE; 73 NUCLEAR PHYSICS AND RADIATION PHYSICS; Fuels Performance; Multiphysics; Neutronics; Thermal Hydraulics

Citation Formats

Gleicher, Frederick, Ortensi, Javier, DeHart, Mark, Wang, Yaqi, Schunert, Sebastian, Novascone, Stephen, Hales, Jason, Williamson, Rich, Slaughter, Andrew, Permann, Cody, Andrs, David, and Martineau, Richard. THE APPLICATION OF MAMMOTH FOR A DETAILED TIGHTLY COUPLED FUEL PIN SIMULATION WITH A STATION BLACKOUT. United States: N. p., 2016. Web.
Gleicher, Frederick, Ortensi, Javier, DeHart, Mark, Wang, Yaqi, Schunert, Sebastian, Novascone, Stephen, Hales, Jason, Williamson, Rich, Slaughter, Andrew, Permann, Cody, Andrs, David, & Martineau, Richard. THE APPLICATION OF MAMMOTH FOR A DETAILED TIGHTLY COUPLED FUEL PIN SIMULATION WITH A STATION BLACKOUT. United States.
Gleicher, Frederick, Ortensi, Javier, DeHart, Mark, Wang, Yaqi, Schunert, Sebastian, Novascone, Stephen, Hales, Jason, Williamson, Rich, Slaughter, Andrew, Permann, Cody, Andrs, David, and Martineau, Richard. Thu . "THE APPLICATION OF MAMMOTH FOR A DETAILED TIGHTLY COUPLED FUEL PIN SIMULATION WITH A STATION BLACKOUT". United States. https://www.osti.gov/servlets/purl/1363897.
@article{osti_1363897,
title = {THE APPLICATION OF MAMMOTH FOR A DETAILED TIGHTLY COUPLED FUEL PIN SIMULATION WITH A STATION BLACKOUT},
author = {Gleicher, Frederick and Ortensi, Javier and DeHart, Mark and Wang, Yaqi and Schunert, Sebastian and Novascone, Stephen and Hales, Jason and Williamson, Rich and Slaughter, Andrew and Permann, Cody and Andrs, David and Martineau, Richard},
abstractNote = {Accurate calculation of desired quantities to predict fuel behavior requires the solution of interlinked equations representing different physics. Traditional fuels performance codes often rely on internal empirical models for the pin power density and a simplified boundary condition on the cladding edge. These simplifications are performed because of the difficulty of coupling applications or codes on differing domains and mapping the required data. To demonstrate an approach closer to first principles, the neutronics application Rattlesnake and the thermal hydraulics application RELAP-7 were coupled to the fuels performance application BISON under the master application MAMMOTH. A single fuel pin was modeled based on the dimensions of a Westinghouse 17x17 fuel rod. The simulation consisted of a depletion period of 1343 days, roughly equal to three full operating cycles, followed by a station blackout (SBO) event. The fuel rod was depleted for 1343 days for a near constant total power loading of 65.81 kW. After 1343 days the fission power was reduced to zero (simulating a reactor shut-down). Decay heat calculations provided the time-varying energy source after this time. For this problem, Rattlesnake, BISON, and RELAP-7 are coupled under MAMMOTH in a split operator approach. Each system solves its physics on a separate mesh and, for RELAP-7 and BISON, on only a subset of the full problem domain. Rattlesnake solves the neutronics over the whole domain that includes the fuel, cladding, gaps, water, and top and bottom rod holders. Here BISON is applied to the fuel and cladding with a 2D axi-symmetric domain, and RELAP-7 is applied to the flow of the circular outer water channel with a set of 1D flow equations. The mesh on the Rattlesnake side can either be 3D (for low order transport) or 2D (for diffusion). BISON has a matching ring structure mesh for the fuel so both the power density and local burn up are copied accurately from Rattlesnake. At each depletion time step, Rattlesnake calculates a power density, fission density rate, burn-up distribution and fast flux based on the current water density and fuel temperature. These are then mapped to the BISON mesh for a fuels performance solve. BISON calculates the fuel temperature and cladding surface temperature based upon the current power density and bulk fluid temperature. RELAP-7 then calculates the fluid temperature, water density fraction and water phase velocity based upon the cladding surface temperature. The fuel temperature and the fluid density are then passed back to Rattlesnake for another neutronics calculation. Six Picard or fixed-point style iterations are preformed in this manner to obtain consistent tightly coupled and stable results. For this paper a set of results from the detailed calculation are provided for both during depletion and the SBO event. We demonstrate that a detailed calculation closer to first principles can be done under MAMMOTH between different applications on differing domains.},
doi = {},
journal = {},
number = ,
volume = ,
place = {United States},
year = {2016},
month = {9}
}

Conference:
Other availability
Please see Document Availability for additional information on obtaining the full-text document. Library patrons may search WorldCat to identify libraries that hold this conference proceeding.

Save / Share: