Multiphysics Simulation of the NASA SIRIUS-CAL Fuel Experiment in the Transient Test Reactor Using Griffin
Abstract
:1. Introduction
2. Overview of the TREAT Facility
3. The SIRIUS-CAL Experiment
4. Modeling Approach
- A detailed Monte Carlo transport solution is used to generate region-homogenized cross-sections and reference fluxes from full 3-D core simulations.
- These cross-sections are converted to the ISOXML format used by Griffin [15]. This cross-section file also carries the reference flux computed in the previous step.
- The ISOXML library is used in a Griffin SPH solver to generate SPH correction factors for each homogenization equivalence region and each energy group.
- Correction factors are added to the ISOXML file.
- Griffin reads cross-sections and correction factors from the ISOXML file, internally computes corrected cross-sections or diffusion coefficients, as appropriate, and performs the transport or diffusion calculation on the homogenized mesh.
4.1. Cross-Section Generation
4.2. Mesh Development
4.3. SIRIUS-CAL Steady-State Calculation
4.4. Griffin Transient Calculation
4.5. Thermal Calculation for the Specimen
- Convective heat transfer of all surfaces around the gas gap with a gas temperature that depends on time, but not on space.
- Radiative heat transfer between surfaces indicated by the black dashed line in Figure 10 with a surface of temperature that depends on time, but not on space. The radiative heat transfer coefficient is computed by assuming that all surfaces indicated in Figure 10 are surrounded by a large enclosure held at temperature .
- CERMET with a volume fraction of 55% UN, 45% W-Re-25:
- -
- -
- -
- Stainless steel 304: All material properties are taken from Ref. [48].
- Stainless steel 316: All material properties are taken from Ref. [53].
- Inconel 718: All material properties are taken from Ref. [48].
- Ti-6Al-4V: All material properties are taken from Ref. [48].
- Microsil insulation: All material properties are taken from Ref. [54].
- Mullite (crucible): Mullite properties depend on its porosity and porosity can vary somewhat. Thermal conductivities vary between 4 and 6 W/m K [55] and densities vary between and 3 for commercial mullite [55,56]. The specific heat is of the order of 250 J/kg K [56]. We adopt a thermal conductivity of 6 W/m K, a density of .
5. Input Uncertainties
5.1. Neutronics Uncertainties
- The value of the first delayed neutron precursor group in ENDF/B-VII.1 and ENDF/B-VIII.0 is , which is different from the measured value of .
- The delayed fission neutron spectrum produced by Serpent is averaged over all delayed neutron groups.
- The adjoint flux used for preparing kinetics parameters should not be SPH-corrected.
- The raw data (i.e., the evaluated nuclear data) used by Serpent as well as the method Serpent uses for computing reactor physics parameters from the raw data create uncertainties for the input data of Griffin.
5.2. Experiment Thermal-Hydraulics Uncertainties
- Cover gas and cavity temperature: The temperature at the boundary of the domain and the temperature of the cover gas are assumed to be uniform in space and equal to the temperature averaged over the surfaces of the TREAT assemblies facing the experiment. The main source of uncertainty is expected to be neglecting temperature distribution in the cover gas and the difference in TREAT and cover gas temperature. Given the relatively low temperature of the experimental peak temperature during this transient, we conservatively assume the cover gas temperature is within K, or approximately %.
- CERMET thermal conductivity: The thermal conductivity of the CERMET is computed from the thermal conductivity of the components WRe25 and UN using Maxwell’s model for the thermal conductivity [51]. The uncertainties are thus made up of the uncertainties of the thermal conductivity of the components, and the uncertainty in using the Maxwell model of thermal conductivity. The range of thermal conductivity for UN is estimated to be % about the regression [48] and the range of WRe25 thermal conductivity is estimated to be % about the regression [50]. From Ref. [51] we estimate the range of thermal conductivities computed with different models for spheres in a matrix to be 10% for the relevant range of volume fractions. The overall thermal conductivity is computed to be % of the nominal value.
- CERMET specific heat: The CERMET specific heat is computed from the specific heat of the constituent materials tungsten, rhenium and UN using the law of mixture. We assume that the law of mixtures does not introduce additional uncertainty for computing specific heat. The upper and lower bounds for CERMET specific heat are computed from the range of measured data for tungsten and UN both of which are around 5% [48] about the regression. No values for rhenium are available, so we assume that rhenium introduces a similar uncertainty as W and UN. The upper and lower bounds for CERMET specific heat are therefore estimated to be 5%.
- Specimen power: The specimen power is uncertain because (1) the TREAT reactor power is uncertain and (2) the CFs may change with specimen and TREAT power. We consider specimen power an input parameter for this study and do not break it down into more atomic sources of uncertainty. For simplicity, the specimen power is bounded between the largest and smallest power curves either computed or measured by any of the TREAT detectors.
- Natural convection heat transfer coefficient: convection heat transfer coefficients for vertical and horizontal surfaces in the experiment cavity are computed using correlations from Ref. [60] for the “isothermal vertical and horizontal surfaces of length l”, respectively. For assessing typical uncertainties in natural heat transfer coefficients, we compare Churchill and Cho’s correlation for vertical plates [61] with the corresponding correlation in Ref. [60] between Rayleigh of and and find that relative differences are between 10 and 15%. We opted to use a range of % for the heat transfer correlations because they are used for geometries that differ somewhat from a vertical plate.
- CERMET density: The percent theoretical density of the experiment specimen is between 70% and 95% according to Ref. [46]. The change in density changes thermal properties and CF values.
- CERMET emissivity: The emissivity of CERMET is assumed to be bounded by the emissivity of pure tungsten taken from Ref. [62] on the low end and the emissivity of pure uranium nitride at small wavelength. Note that the emissivity of uranium nitride is somewhere between 0.7 and 0.8, but for this study the bounding value of is assumed.
- Metal emissivities: Stainless steel emissivities are adopted from Ref. [62], the range of Inconel® emissivities are adopted from Ref. [63] and titanium grade 5 emissivities are adopted from [64]. The difference between the upper and lower bound of metal surface emissivities is mostly caused by the surface finish (polished versus oxidized).
- Mullite emissivity: The emissivity of mullite is well characterized in Ref. [65] to within .
- Cover gas conductivity: The uncertainty of the cover gas thermal conductivity is driven by not knowing the temperature of the cover gas well. The upper and lower limits are obtained by evaluating the thermal conductivity correlation Equation (2) at 300 and 1000 K.
6. Simulation Results
6.1. Neutronics Solutions
6.2. Specimen Temperature Transient
- Local: Experiment power driven by reactor power (i.e., the power coupling factor), experiment thermal properties, thermal-hydraulics conditions such as heat transfer coefficients.
- Global: Reactor power as a function of time, temperature at the outside of the containment (TREAT fans operate during a transient, drawing air from the building above the top of the reactor and blowing down through the fuel elements).
Validation and Uncertainty Quantification
- Cover gas thermal conductivity has the largest of up to 40 K. The cover gas thermal conductivity is important because the experiment is small and temperature levels are below 1000 K so that thermal conductivity through the cover gas is an important heat transfer mechanism. This uncertainty can be reduced significantly by increasing the knowledge of the cover gas temperature that drives the uncertainty in cover gas thermal conductivity. For that purpose, a thermal-hydraulics or CFD model of the cavity cover gas should be added to the SIRIUS-CAL experiment model.
- The specimen power accounts for about 35 K of uncertainty.
- The CERMET nominal density contributes up to 25 K to the uncertainty. This uncertainty stems entirely from lack of information from the experimenters and can be reduced by a more careful pre-experiment characterization of the specimen.
- The CERMET specific heat contributes up to 15 K. This uncertainty can be reduced by measurements and modeling and simulation of the material properties of the CERMET such as performed in [69].
- The heat transfer coefficient on the vertical surfaces contributes up to K. The uncertainty of heat transfer coefficients can be reduced by improving the fidelity of the cover gas thermal-hydraulics model, using a CFD model for the cover gas region, or using CFD to inform a more accurate thermal-hydraulics model with closure relationships such as heat transfer coefficients.
- The CERMET thermal conductivity contributes up to K. The same comments as for the specific heat applies.
- The cover gas temperature contributes up to K. Note that cover gas thermal conductivity and cover gas temperature are changed independently so this sensitivity is only to the temperature of the cover gas assuming that the thermal conductivity does not change. The knowledge of the cover gas temperatures can be improved by using a better thermal-hydraulics model of the cover gas region.
- The effects of emissivities and heat transfer coefficients on the horizontal surfaces are negligible.
- The specimen power is affected by the reactor power and the power coupling factor. In this work, we assume that the PCF does not change during the transient. Therefore, uncertainties are entirely due to the reactor power. The prediction of the reactor power with Griffin is an active field of work. The two main problems identified in Ref. [19] are the over-prediction of control rod worth by the existing Serpent model (similarly TREAT reactor operation’s MCNP models also overestimate control rod worth) and the discrepancy of neutron kinetics data in SPH corrected neutron diffusion calculations. The issue with effective kinetics data in the Serpent–Griffin–SPH data pipeline are discussed in Ref. [59].
- Material properties and specimen characterization are outside of the realm of neutronics/ thermal/thermal-hydraulics modeling of TREAT experiments. The uncertainties from these sources are large for the SIRIUS-CAL experiment because of unusually poor characterization of the specimen density and the relative novelty of the W-Re-25 CERMET. Future experiments, even within the SIRIUS series, will have more accurate material properties available.
- The thermal-hydraulics assumptions in the presented model lead to large uncertainties. In particular, gas temperature in the cavity and radiative transfer in the cavity (e.g., via a net radiation transfer method) would reduce the uncertainties significantly.
7. Conclusions
Author Contributions
Funding
Data Availability Statement
Acknowledgments
Conflicts of Interest
References
- Levack, D.J.; Horton, J.F.; Joyner, C.R.; Kokan, T.S.; Widman, F.; Guzek, B.J. Mars NTP Architecture Elements Using the Lunar Orbital Platform-Gateway. In Proceedings of the 2018 AIAA SPACE and Astronautics Forum and Exposition, Orlando, FL, USA, 17–19 September 2018. [Google Scholar] [CrossRef]
- DeHart, M.D.; Schunert, S.; Labouré, V.M. Nuclear Thermal Propulsion. In Nuclear Reactors; Pope, C.L., Ed.; IntechOpen: London, UK, 2022. [Google Scholar] [CrossRef]
- National Academies of Sciences, Engineering, and Medicine. Space Nuclear Propulsion for Human Mars Exploration; The National Academies Press: Washington, DC, USA, 2021. [Google Scholar] [CrossRef]
- Borowski, S.K.; McCurdy, D.R.; Packard, T.W. Nuclear Thermal Propulsion (NTP): A proven growth technology for human NEO/Mars exploration missions. In Proceedings of the 2012 IEEE Aerospace Conference, Big Sky, MT, USA, 3–10 March 2012; pp. 1–20. [Google Scholar] [CrossRef]
- Johnson, L.; Jones, J.; Kos, L.; Trausch, A.; Farris, R.; Woodcock, G. Benefits of Nuclear Electric Propulsion to Outer Planets Exploration. In Proceedings of the 38th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, Indianapolis, IN, USA, 7–10 July 2002. [Google Scholar] [CrossRef]
- Heath, B.K.; Race, C.C. TREAT Restart Project. Nucl. Technol. 2019, 205, 1369–1377. [Google Scholar] [CrossRef]
- Houts, M.G.; Kim, T.; Emrich, W.J.; Hickman, R.R.; Broadway, J.W.; Gerrish, H.P.; Doughty, G.; Belvin, A.; Borowski, S.K.; Scott, J. The Nuclear Cryogenic Propulsion Stage. In Proceedings of the Nuclear and Emerging Technologies for Space (NETS) 2014, Hancock County, IN, USA, 24–26 February 2014. American Nuclear Society Topical Meeting. [Google Scholar]
- Gustafson, J.L. Space Nuclear Propulsion Fuel and Moderator Development Plan Conceptual Testing Reference Design. Nucl. Technol. 2021, 207, 882–884. [Google Scholar] [CrossRef]
- Esselman, W.H. The NERVA Nuclear Rocket Reactor Program. Westinghouse Eng. 1965, 23, 66–75. [Google Scholar]
- Tucker, D. CERMETS for Use in Nuclear Thermal Propulsion; IntechOpen: London, UK, 2019. [Google Scholar] [CrossRef]
- O’BRIEN, R.; Burns, D.E.; Woolstenhulme, N.E.; Dempsey, D.; Hill, C.M.; Hale, C. The SIRIUS-1 Nuclear Thermal Propulsion Fuels Transient Test Series in the Idaho National Laboratory TREAT Reactor. In Proceedings of the Nuclear and Emerging Technologies for Space (NETS) 2019, Richland, WA, USA, 25–28 February 2019. American Nuclear Society Topical Meeting. [Google Scholar]
- DeHart, M.D.; Baker, B.A.; Ortensi, J. Interpretation of energy deposition data from historical operation of the transient test facility (TREAT). Nucl. Eng. Des. 2017, 322, 504–521. [Google Scholar] [CrossRef]
- Ortensi, J.; Wang, Y.; Laboure, V.M.; Prince, Z.M.; Lee, C.; Jung, Y.S.; Park, H.; Shemon, E.R. Griffin Software Development Plan; Technical Report, INL/EXT-21-63185; Idaho National Laboratory: Idaho Falls, ID, USA, 30 June 2019. [Google Scholar] [CrossRef]
- Laboure, V.M.; Ortensi, J.; Wang, Y.; Schunert, S.; Gleicher, F.N.; DeHart, M.D.; Martineau, R.C. Multiphysics Steady-state simulation of the High Temperature Test Reactor with MAMMOTH, BISON and RELAP-7. In Proceedings of the International Conference on Mathematics Computational Methods Applied to Nuclear Science and Engineering (M&C 2019), Portland, OR, USA, 25–29 August 2019. [Google Scholar]
- Wang, Y.; Schunert, S.; Ortensi, J.; Laboure, V.; DeHart, M.; Prince, Z.; Gleicher, F. Rattlesnake: A MOOSE-Based Multiphysics Multischeme Radiation Transport Application. Nucl. Technol. 2021, 207, 1047–1072. [Google Scholar] [CrossRef]
- Shemon, E.R.; Smith, M.A.; Lee, C. PROTEUS-SN User Manual Revision 1.2; Technical Report ANL/NE-14/6; Argonne National Laboratory: Argonne, IL, USA, 2014. [Google Scholar]
- Lee, C.; Yang, W.S. MC2-3: Multigroup Cross Section Generation Code for Fast Reactor Analysis; Technical Report ANL/NE-11/41 (Rev. 3); Argonne National Laboratory: Argonne, IL, USA, 2018. [Google Scholar]
- Idaho National Laboratory. MOOSE Heat Conduction Module. Available online: https://mooseframework.inl.gov/modules/index.html (accessed on 6 May 2022).
- Ortensi, J.; Baker, B.A.; Johnson, M.P.; Wang, Y.; Labouré, V.M.; Schunert, S.; Gleicher, F.N.; DeHart, M.D. Validation of the Griffin application for TREAT transient modeling and simulation. Nucl. Eng. Des. 2021, 385, 111478. [Google Scholar] [CrossRef]
- MacFarlane, D.R.; Freund, G.A.; Boland, J.F. Hazards Summary Report on the Transient Reactor Test Facility (TREAT); Technical Report ANL-5923; Argonne National Laboratory: Argonne, IL, USA, 1958. [Google Scholar]
- Bess, J.D.; DeHart, M.D. Baseline Assessment of TREAT for Modeling and Analysis Needs; Technical Report INL/EXT-15-35372; Idaho National Laboratory: Idaho Falls, ID, USA, 2015. [Google Scholar]
- Okrent, D.; Dickerman, C.E.; Gasidlo, J.; O’Shea, D.M.; Schoeberl, D.F. The Reactor Kinetics of the Transient Reactor Test Facility (TREAT); Technical Report ANL-6174; Argonne National Laboratory: Argonne, IL, USA, 1960. [Google Scholar]
- Baker, B.A.; Fielding, K.D.; Hansen, J.E.; Ellsworth, T. Enhancements to the New TREAT Automatic Reactor Control System (ARCS). Nucl. Technol. 2021, 207, 1746–1752. [Google Scholar] [CrossRef]
- Woolstenhulme, N.E. Status Report on Development of MARCH Module Designs; Idaho National Laboratory: Idaho Falls, ID, USA, 2019. [Google Scholar] [CrossRef]
- Woolum, C.; Imholte, D.D.; Fleming, A.; Kamerman, D.; Tritthart, K. Fabrication and Assembly of the First Accident Tolerant Fuel Concept for TREAT Testing. In Proceedings of the Pressure Vessels and Piping Conference, Online, 19–24 July 2020; Volume 6. Materials and Fabrication. [Google Scholar] [CrossRef]
- Idaho National Laboratory. Calibration Transient for the Sirius Experiment (SIRIUS-CAL); Technical Report TFR-1040; Idaho National Laboratory: Idaho Falls, ID, USA, 2019. [Google Scholar]
- Labouré, V.; Wang, Y.; Ortensi, J.; Schunert, S.; Gleicher, F.; DeHart, M.; Martineau, R. Hybrid super homogenization and discontinuity factor method for continuous finite element diffusion. Ann. Nucl. Energy 2019, 128, 443–454. [Google Scholar] [CrossRef]
- Ortensi, J.; Wang, Y.; Laurier, A.; Schunert, S.; Hébert, A.; DeHart, M. A Newton Solution for the Superhomogenization Method: The PJFNK-SPH. Ann. Nucl. Energy 2018, 111, 579–594. [Google Scholar] [CrossRef]
- Hébert, A. Développement de la Methode SPH: Homogénéisation de Cellules dans Un réseau non Uniforme et Calcul des Paramètres de Réflecteur; Technical Report CEA-N-2209; Commissariat à l’Energie Atomique Centre d’Etudes Nucleaires de Saclay: Gif-sur-Yvette, France, 1981. [Google Scholar]
- Prince, Z.M.; Ragusa, J.C.; Wang, Y. Implementation of the Improved Quasi-Static Method in RATTLESNAKE/MOOSE for Time-dependent Radiation Transport Modelling. In Proceedings of the PHYSOR 2016: Unifying Theory and Experiments in the 21st Century, Sun Valley, ID, USA, 1–5 May 2016. [Google Scholar]
- Leppänen, J.; Pusa, M.; Viitanen, T.; Valtavirta, V.; Kaltiaisenaho, T. The Serpent Monte Carlo code: Status, development and applications in 2013. Ann. Nucl. Energy 2012, 82, 142–150. [Google Scholar] [CrossRef]
- Brown, D.A.; Chadwick, M.B.; Capote, R.; Kahler, A.C.; Trkov, A.; Herman, M.W.; Zhu, Y. ENDF/B-VIII.0: The 8th Major Release of the Nuclear Reaction Data Library with CIELO-project Cross Sections, New Standards and Thermal Scattering Data. Nucl. Data Sheets 2018, 148, 1–142. [Google Scholar] [CrossRef]
- Ortensi, J.; Alberti, A.; Wang, Y.; DeHart, M.D.; Gleicher, F.N.; Schunert, S.; Palmer, T.S. Methodologies and Requirements for the Generation of Physics Data Inputs to MAMMOTH Transient Simulations in Support of the Transient Reactor Test Facility; Technical Report INL/LTD-15-36265; Idaho National Laboratory: Idaho Falls, ID, USA, 2015. [Google Scholar]
- Ortensi, J.; Baker, B.; Schunert, S.; Wang, Y.; Gleicher, F.N.; DeHart, M.D. Updates to the Generation of Physics Data Inputs for MAMMOTH Simulations of the Transient Reactor Test Facility—FY2016; Technical Report INL/EXT-16-39120; Idaho National Laboratory: Idaho Falls, ID, USA, 2016. [Google Scholar]
- Ortensi, J.; DeHart, M.D.; Gleicher, F.N.; Wang, Y.; Schunert, S.; Alberti, A.L.; Palmer, T.S. Full Core TREAT Kinetics Demonstration Using Rattlesnake/BISON Coupling within MAMMOTH; Technical Report INL/EXT-15-36268; Idaho National Laboratory: Idaho Falls, ID, USA, 2015. [Google Scholar]
- Ortensi, J.; Labouré, V.; Schunert, S.; Stewart, R. Neutronics Enhanced Meshing Operations (NEMO) Users Manual Version 1.0; Technical Report INL/MIS-20-00730; Idaho National Laboratory: Idaho Falls, ID, USA, 2020. [Google Scholar]
- Blacker, T.D.; Owen, S.J.; Staten, M.K.; Quadros, W.R.; Hanks, B.; Clark, B.W.; Meyers, R.J.; Ernst, C.; Merkley, K.; Morris, R.; et al. CUBIT Geometry and Mesh Generation Toolkit 15.2 User Documentation; Technical Report, SAND2016-6850R; Sandia National Laboratories: Albuquerque, NM, USA, 1 May 2016. [Google Scholar]
- Schoof, L.A.; Yarberry, V.R. EXODUS II: A Finite Element Data Model, SAND92-2137; Technical Report; Sandia National Laboratory: Albuquerque, NM, USA, 1995. [Google Scholar]
- Weems, Z.; Goluoglu, S.; DeHart, M.D. Investigations of Rod Positions for Treat M8CAL Analyses. ASME J. Nucl. Radiat. Sci. 2018, 4, 031017. [Google Scholar] [CrossRef] [Green Version]
- Goluoglu, S.; DeHart, M.D. On the effects of pre- and post-transient rod positions for TREAT temperature-limited transient powers. Anuclear Eng. Des. 2018, 331, 97–102. [Google Scholar] [CrossRef]
- Idaho National Laboratory. MOOSE Executioner System. Available online: https://mooseframework.inl.gov/syntax/Executioner/index.html (accessed on 6 May 2022).
- Prince, Z.M.; Ragusa, J.C. Multiphysics reactor-core simulations using the improved quasi-static method. Ann. Nucl. Energy 2019, 125, 186–200. [Google Scholar] [CrossRef]
- Handwerk, J.; Lied, R. The Manufacture of the Graphite-Urania Fuel Matrix for TREAT; Technical Report ANL-5963; Argonne National Laboratory: Argonne, IL, USA, 1960. [Google Scholar]
- Gaston, D.R.; Permann, C.J.; Peterson, J.W.; Slaughter, A.E.; Andrš, D.; Wang, Y.; Short, M.P.; Perez, D.M.; Tonks, M.R.; Ortensi, J.; et al. Physics-based multiscale coupling for full core nuclear reactor simulation. Ann. Nucl. Energy 2015, 84, 45–54. [Google Scholar] [CrossRef]
- X-5 Monte Carlo Team. MCNP—A General Monte Carlo N-Particle Transport Code, Version 5; Technical Report; Los Alamos National Laboratory: Los Alamos, NM, USA, 2000. [Google Scholar]
- Gross, B.J. Physics Analysis of SIRIUS-CAL Experiment Irradiated in TREAT; Technical Report ECAR-4443; Idaho National Laboratory: Idaho Falls, ID, USA, 2019. [Google Scholar]
- Williamson, R.L.; Hales, J.D.; Novascone, S.R.; Tonks, M.R.; Gaston, D.R.; Permann, C.J.; Andrs, D.; Martineau, R.C. Multidimensional multiphysics simulation of nuclear fuel behavior. J. Nucl. Mater. 2012, 423, 149–163. [Google Scholar] [CrossRef]
- Milner, J. Space Nuclear Propulsion Material Property Handbook; Technical Report SNP-HDBK-0008; Glenn Research Center: Cleveland, OH, USA.
- Rhenium. Available online: https://en.wikipedia.org/wiki/Rhenium (accessed on 25 April 2020).
- Wang, Y.; Zhao, J. Understanding the thermal conductivity of pristine W and W-Re alloys from a physics-based model. J. Nucl. Mater. 2020, 529, 151931. [Google Scholar] [CrossRef]
- Toptan, A.; Jiang, W.; Hales, J.D.; Spencer, B.W.; Casagranda, A.; Novascone, S.R. FEA-aided investigation of the effective thermal conductivity in a medium with embedded spheres. Nucl. Eng. Des. 2021, 381, 111355. [Google Scholar] [CrossRef]
- Taylor, R.; Finch, R. The specific heats and resistivities of molybdenum, tantalum, and rhenium. J. Less Common Met. 1964, 6, 283–294. [Google Scholar] [CrossRef]
- Kim, C.S. Thermophysical Properties of Stainless Steels; Argonne National Laboratory: Argonne, IL, USA, 1975. [Google Scholar] [CrossRef]
- Microsil Data Sheet. Available online: https://www.zircarceramics.com/wp-content/uploads/2017/01/MICROSIL.pdf (accessed on 25 April 2020).
- Mullite—An Introduction. Available online: https://www.azom.com/article.aspx?ArticleID=925 (accessed on 25 April 2020).
- Mindrum Cermamics. Available online: https://www.mindrum.com/ceramics/mullite/ (accessed on 25 April 2020).
- Siefken, L. SCDAP/RELAP5/MOD 3.3 Code Manual; Technical Report NUREG/CR-6150; Volume 4, Rev. 2, NRC. Available online: http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.409.5896&rep=rep1&type=pdf (accessed on 25 April 2020).
- Baker, B. Rod characterization in Serpent; unpublished.
- Takasugi, C.; Martin, N.; Laboure, V.; Ortensi, J.; Ivanov, K.; Avramova, M. Preservation of Kinetics Parameters Generated by Monte Carlo Calculations in Two-Step Deterministic Calculations. EPJ Nucl. Sci. Technol. 2022. accepted. [Google Scholar]
- Zhang, H.; Che, F.; Lin, T.; Zhao, W. 3—Thermal modeling, analysis, and design. In Modeling, Analysis, Design, and Tests for Electronics Packaging beyond Moore; Zhang, H., Che, F., Lin, T., Zhao, W., Eds.; Woodhead Publishing Series in Electronic and Optical Materials; Woodhead Publishing: Sawston, UK, 2020; pp. 59–129. [Google Scholar] [CrossRef]
- Churchill, S.W.; Chu, H.H. Correlating equations for laminar and turbulent free convection from a horizontal cylinder. Int. J. Heat Mass Transf. 1975, 18, 1049–1053. [Google Scholar] [CrossRef]
- Modest, M.F.; Mazumder, S. Radiative Heat Transfer; Academic Press: Cambridge, MA, USA, 2021. [Google Scholar]
- Keller, B.P.; Nelson, S.E.; Walton, K.L.; Ghosh, T.K.; Tompson, R.V.; Loyalka, S.K. Total hemispherical emissivity of Inconel 718. Nucl. Eng. Des. 2015, 287, 11–18. [Google Scholar] [CrossRef]
- González-Fernández, L.; Risueño, E.; Pérez-Sáez, R.; Tello, M. Infrared normal spectral emissivity of Ti–6Al–4V alloy in the 500–1150 K temperature range. J. Alloys Compd. 2012, 541, 144–149. [Google Scholar] [CrossRef]
- Optotherm Thermal Imaging. Optotherm Support, Emissivity in the Infrared. 2018. Available online: https://www.optotherm.com/emiss-table.htm (accessed on 6 May 2022).
- Carvajal Nunez, U.; Prieur, D.; Bohler, R.; Manara, D. Melting point determination of uranium nitride and uranium plutonium nitride: A laser heating study. J. Nucl. Mater. 2014, 449, 1–8. [Google Scholar] [CrossRef]
- Geelhood, K.J.; Luscher, W.G.; Senor, D.J.; Cunningham, M.E.; Lanning, D.D.; Adkins, H.E. Predictive Bias and Sensitivity in NRC Fuel Performance Codes; Technical Report NUREG/CR-7001; Pacific Northwest National Laboratory: Richland, WA, USA, 2009. [Google Scholar]
- Huntington, D.; Lyrintzis, C. Improvements to and limitations of Latin hypercube sampling. Probabilistic Eng. Mech. 1998, 13, 245–253. [Google Scholar] [CrossRef]
- Sessim, M.; Tonks, M.R. Multiscale simulations of thermal transport in W-UO2 CERMET Fuel for nuclear thermal propulsion. Nucl. Technol. 2021, 207, 1004–1014. [Google Scholar] [CrossRef]
- Spence, R.W. The Rover Nuclear Rocket Program. Science 1968, 160, 953–959. [Google Scholar] [CrossRef]
- Stewart, M.E. A Historical Review of Cermet Fuel Development and the Engine Performance Implications. In Proceedings of the Nuclear and Emerging Technologies for Space (NETS) 2015, Albuquerque, NM, USA, 23–26 February 2015. American Nuclear Society Topical Meeting. [Google Scholar]
Group | Upper Energy Bound (eV) | Group | Upper Energy Bound (eV) |
---|---|---|---|
1 | 4.00000 × 107 | 6 | 6.25000 × 10−1 |
2 | 1.15620 × 106 | 7 | 2.09610 × 10−1 |
3 | 3.48110 × 103 | 8 | 7.64970 × 10−2 |
4 | 1.32700 × 102 | 9 | 4.73020 × 10−2 |
5 | 8.10003 × 100 | 10 | 2.00100 × 10−2 |
Thermal Property | Unit | Value |
---|---|---|
Mass density | g/cm3 | 1.53 |
Heat conduction coefficient | W/cm-K | + |
Component | CF Total (Low Density) | CF Total (High Density) |
---|---|---|
Fuel specimen | 1.84 × 100 | 1.61 × 100 |
Hanger rod | 2.60 × 10 | 2.75 × 10 |
Bottom collar | 2.92 × 10 | 2.92 × 10 |
Top collar | 2.51 × 10 | 3.03 × 10 |
Flask | 3.04 × 10 | 3.13 × 10 |
Crucible | 2.52 × 10 | 2.40 × 10 |
Capsule | 2.79 × 10 | 2.78 × 10 |
BUSTER pipe | 2.52 × 10 | 2.51 × 10 |
Insulation | 1.94 × 10−3 | 1.95 × 10−3 |
Uncertainty | Reactivity Effect (pcm) | Control Rod Position | |
---|---|---|---|
Core boron | 6–7.6 ppm | 1200 | C |
Core height | 0.3172 cm | 150 | C |
CR C density | 1.4–1.6 | 100 | V |
Zr-Steel composition | - | 140 | V |
C/S rods offset | 1.27 cm | 200 | V |
Transient rods offset | 0.762 cm | 150 | V |
Parameter | Variation | Reference |
---|---|---|
Cover gas and | see text | |
cavity temperature | ||
CERMET thermal conduct | % nominal | see text |
CERMET specific heat | % nominal | see text |
Specimen power | see text | see text |
Natural convection heat | % | see text |
transfer coefficient | ||
CERMET density | 70–95% | [46] |
theoretical density | ||
CERMET emissivity | 0.4–0.8 | Low value W-Re from [62] |
High value is | ||
at small wavelength [66] | ||
SS316 & 304 emissivity | 0.26–0.66 | [62] |
Inconel® 718 emissivity | 0.2–0.53 | [63] |
Titanium grade 5 emissivity | 0.3–0.47 | [64] |
Mullite emissivity | 0.8–0.85 | [65] |
Cover gas conductivity | 0.15–0.35 W/mK | [67] |
Transient Rod Position, cm | Reactivity Insertion, %k/k | Asymptotic Period, s | |
---|---|---|---|
Measurement | 85.598 | 0.55 | 1.826 |
Griffin, without adjustment | 85.598 | 0.63 | 0.892 |
Griffin, reactivity match | 87.075 | 0.55 | 1.460 |
Griffin, period match | 87.878 | 0.51 | 1.826 |
Quantity | Reported | Griffin | Relative Difference |
---|---|---|---|
, total | 7.178 × 10−3 | 6.932 × 10−3 | −3.43% |
, s | 9.000 × 10−4 | 8.928 × 10−4 | −0.80% |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Jing, T.; Schunert, S.; Labouré, V.M.; DeHart, M.D.; Lin, C.-S.; Ortensi, J. Multiphysics Simulation of the NASA SIRIUS-CAL Fuel Experiment in the Transient Test Reactor Using Griffin. Energies 2022, 15, 6181. https://doi.org/10.3390/en15176181
Jing T, Schunert S, Labouré VM, DeHart MD, Lin C-S, Ortensi J. Multiphysics Simulation of the NASA SIRIUS-CAL Fuel Experiment in the Transient Test Reactor Using Griffin. Energies. 2022; 15(17):6181. https://doi.org/10.3390/en15176181
Chicago/Turabian StyleJing, Tian, Sebastian Schunert, Vincent M. Labouré, Mark D. DeHart, Ching-Sheng Lin, and Javier Ortensi. 2022. "Multiphysics Simulation of the NASA SIRIUS-CAL Fuel Experiment in the Transient Test Reactor Using Griffin" Energies 15, no. 17: 6181. https://doi.org/10.3390/en15176181
APA StyleJing, T., Schunert, S., Labouré, V. M., DeHart, M. D., Lin, C. -S., & Ortensi, J. (2022). Multiphysics Simulation of the NASA SIRIUS-CAL Fuel Experiment in the Transient Test Reactor Using Griffin. Energies, 15(17), 6181. https://doi.org/10.3390/en15176181