Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reactor Physics Method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Catalyst support structure, catalyst including the structure, reactor including a catalyst, and methods of forming same

Structures, catalysts, and reactors suitable for use for a variety of applications, including gas-to-liquid and coal-to-liquid processes and methods of forming the structures, catalysts, and reactors are disclosed. The catalyst material can be deposited onto an inner wall of a microtubular reactor and/or onto porous tungsten support structures using atomic layer deposition techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Application of Generalized Standard Addition Method and Ultraviolet Spectroscopy to Quantify Electrolytic Depolymerization of Lignin

Industrial lignin waste is an important byproduct of bio-refineries and the paper industry. Depolymerization of this industrial lignin could generate useful aromatic compounds. This group has focused on electrolytic decomposition of biorefinery lignin. To quantify electrolytic decomposition of the lignin in a highly caustic solution, ultraviolet (UV) spectroscopy provides a useful probe. The conversion of the neat lignin to the oxidized products achieved by the electrolytic reactor may be measured by quantifying the amount of unreacted neat lignin that remains in the effluent. Because the properties of electrolytic decomposition products are largely unknown, a useful approach to quantify decomposition of the neat lignin is to use a multivariate calibration method referred to as the generalized standard addition method (GSAM). In this approach, the electrolytic decomposition products represent a background interference and the neat lignin that remains can be quantified. Furthermore, this approach allows the conversion of the neat lignin to be calculated in a solution that is a complex mixture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Asynchronous τ-leaping

Stochastic simulation of cell signaling pathways and genetic regulatory networks has contributed to the understanding of cell function; however, investigation of larger, more complicated systems requires computationally efficient algorithms. τ-leaping methods, which improve efficiency when some molecules have high copy numbers, either use a fixed leap size, which does not adapt to changing state, or recalculate leap size at a heavy computational cost. We present a hybrid simulation method for reaction-diffusion systems which combines exact stochastic simulation and τ-leaping in a dynamic way. Putative times of events are stored in a priority queue, which reduces the cost of each step of the simulation. For every reaction and diffusion channel at each step of the simulation the more efficient of an exact stochastic event or a τ-leap is chosen. This new approach removes the inherent trade-off between speed and accuracy in stiff systems which was present in all τ-leaping methods by allowing each reaction channel to proceed at its own pace. Both directions of reversible reactions and diffusion are combined in a single event, allowing bigger leaps to be taken. This improves efficiency for systems near equilibrium where forward and backward events are approximately equally frequent. Comparison with existing algorithms and behaviour for five test cases of varying complexity shows that the new method is almost as accurate as exact stochastic simulation, scales well for large systems, and for various problems can be significantly faster than τ-leaping.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

TREAT Neutronics Analysis and Design Support, Part I: Multi-SERTTA

Experiment vehicle design is necessary in preparation for Transient Reactor Test (TREAT) facility restart and the resumption of transient testing to support Accident Tolerant Fuel (ATF) characterization and other future fuels testing requirements. Currently the most mature vehicle design is the Multi-SERTTA (Static Environments Rodlet Transient Test Apparatuses), which can accommodate up to four concurrent rodlet-sized specimens under separate environmental conditions. Robust test vehicle design requires neutronics analyses to support design development, optimization of the power coupling factor (PCF) to efficiently maximize energy generation in the test fuel rodlets, and experiment safety analyses. Calculations were performed to support analysis of a near-final design of the Multi-SERTTA vehicle, the design process for future TREAT test vehicles, and establish analytical practices for upcoming transient test experiments. Models of the Multi-SERTTA vehicle containing typical PWR-fuel rodlets were prepared and neutronics calculations were performed using MCNP6.1 with ENDF/B-VII.1 nuclear data libraries. Calculation of the PCF for reference conditions of a PWR fuel rodlet in clean water at operational temperature and pressure provided results between 1.10 and 1.74 W/g-MW depending on the location of the four Multi-SERTTA units with the stack. Basic changes to the Multi-SERTTA secondary vessel containment and support have minimal impact on PCF; using materials with less neutron absorption can improve expected PCF values, especially in the primary containment. An optimized balance is needed between structural integrity, experiment safety, and energy deposition in the experiment. Type of medium and environmental conditions within the primary vessel surrounding the fuel rodlet can also have a significant impact on resultant PCF values. The estimated reactivity insertion worth into the TREAT core is impacted more by the primary and secondary Multi-SERTTA vehicle structure with the experiment content and contained environment having a near negligible impact on overall system reactivity. Additional calculations were performed to evaluate the peak-to-average assembly powers throughout the TREAT core, as well as the nuclear heat generation for the various structural components of the Multi-SERTTA assembly. Future efforts include the evaluation of flux collars to shape the PCF for individual Multi-SERTTA units during an experiment such as to achieve uniformity in test unit environmental conditions impacted by the non-uniform axial flux/power profile of TREAT. Upon resumption of transient testing, experimental results from both the Multi-SERTTA and Multi-SERTTA-CAL will be compared against calculational results and methods for further optimization and design strategies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Design workflow of a symmetric traveling wave antenna for fast ion production on DD tokamaks

Initial computational plasma physics scoping and a finite element method antenna modeling design workflow for a symmetric center-fed high-field side high harmonic fast wave traveling wave array (TWA) antenna are reported here. The TWA is designed to generate a test population of fast deuterium ions in an existing D–D tokamak by heating neutral beam deuterium ions, accelerating them from 80 keV to several hundred keV. The resulting fast particles are tailored to mimic key reactor energetic particle parameters with regards to exciting Alfven eigenmode instabilities, allowing for a D–D tokamak like DIII-D or ASDEX-U to replicate reactor-relevant conditions experimentally. Initial scenario scoping for high single-pass absorption as well as good preferential fast ion damping relative to electron damping was completed using the ray-tracing/Fokker–Planck codes GENRAY and CQL3D. Python RF network analysis packages were used to create a custom TWA optimization tool to inform a COMSOL flat antenna design, and Petra-M was used to study cold plasma effects. The TWA produced by this workflow has several novel features when compared to previous TWA studies, including symmetric center feeding, and passive end straps for image current cancellation for reduced impurity production. We show here that the antenna design workflow can readily produce TWA antennas optimized for reflection coefficient, image current cancellation, and launched power spectrum shape; and that a population of fast ions can be generated in the correct region of parameter space, warranting future more detailed studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Enhancements to the SHARP Build System and NEK5000 Coupling

The SHARP project for the Department of Energy's Nuclear Energy Advanced Modeling and Simulation (NEAMS) program provides a multiphysics framework for coupled simulations of advanced nuclear reactor designs. It provides an overall coupling environment that utilizes custom interfaces to couple existing physics codes through a common spatial decomposition and unique solution transfer component. As of this writing, SHARP couples neutronics, thermal hydraulics, and structural mechanics using PROTEUS, Nek5000, and Diablo respectively. This report details two primary SHARP improvements regarding the Nek5000 and Diablo individual physics codes: (1) an improved Nek5000 coupling interface that lets SHARP achieve a vast increase in overall solution accuracy by manipulating the structure of the internal Nek5000 spatial mesh, and (2) the capability to seamlessly couple structural mechanics calculations into the framework through improvements to the SHARP build system. The Nek5000 coupling interface now uses a barycentric Lagrange interpolation method that takes the vertex-based power and density computed from the PROTEUS neutronics solver and maps it to the user-specified, general-order Nek5000 spectral element mesh. Before this work, SHARP handled this vertex-based solution transfer in an averaging-based manner. SHARP users can now achieve higher levels of accuracy by specifying any arbitrary Nek5000 spectral mesh order. This improvement takes the average percentage error between the PROTEUS power solution and the Nek5000 interpolated result down drastically from over 23 % to just above 2 %, and maintains the correct power profile. We have integrated Diablo into the SHARP build system to facilitate the future coupling of structural mechanics calculations into SHARP. Previously, simulations involving Diablo were done in an iterative manner, requiring a large amount manual work, and left only as a task for advanced users. This report will detail a new Diablo build system that was implemented using GNU Autotools, mirroring much of the current SHARP build system, and easing the use of structural mechanics calculations for end-users of the SHARP multiphysics framework. It lets users easily build and use Diablo as a stand-alone simulation, as well as fully couple with the other SHARP physics modules. The top-level SHARP build system was modified to allow Diablo to hook in directly. New dependency handlers were implemented to let SHARP users easily build the framework with these new simulation capabilities. The remainder of this report will describe this work in full, with a detailed discussion of the overall design philosophy of SHARP, the new solution interpolation method introduced, and the Diablo integration work. We will conclude with a discussion of possible future SHARP improvements that will serve to increase solution accuracy and framework capability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Flow-through reactor for electrocatalytic reactions

A flow-through electrolysis cell includes a hierarchical nanoporous metal cathode. A method of reducing CO 2 includes flowing the CO 2 through the hierarchical nanoporous metal cathode of the flow-through electrolysis cell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Biofuels Refining Engineering

The goal of this project is the development of novel catalysts and knowledge of reaction pathways and mechanisms for conversion of biomass-based compounds to fuels that are compatible with oil-based fuels and with acceptable or superior fuel properties. The research scope included both catalysts to convert lignocellulosic biomass-based molecules (from pyrolysis) and vegetable oil-based molecules (i.e., triglycerides and fatty acid methyl esters). This project comprised five technical tasks. Each task is briefly introduced below, and major technical accomplishments summarized. Technical accomplishments were described in greater detail in the quarterly progress reports, and in even more detail in the >50 publications acknowledging this DoE project funding (list of publications and presentations included at the end of this report). The results of this research added greatly to the knowledge base necessary for upgrading of pyrolysis oil to hydrocarbon fuels and chemicals, and for conversion of vegetable oils to fungible diesel fuel. Numerous new catalysts and catalytic reaction systems were developed for upgrading particular compounds or compound families found in the biomass-based pyrolysis oils and vegetable oils. Methods to mitigate catalyst deactivation were investigated, including novel reaction/separation systems. Performance and emission characteristics of biofuels in flames and engines were measured. Importantly, the knowledge developed from this project became the basis for a subsequent collaborative proposal led by our research group, involving researchers from the University of Wisconsin, the University of Pittsburg, and the Idaho National Lab, for the DoE Carbon, Hydrogen and Separations Efficiency (CHASE) program, which was subsequently funded (one of only four projects awarded in the CHASE program). The CHASE project examined novel catalytic processes for lignocellulosic biomass conversion as well as technoeconomic analyses for process options for maximum carbon capture and hydrogen efficiency. Our research approach combined catalyst synthesis, measurements of catalyst activity and selectivity in different reactor systems and conditions, and detailed catalyst characterization to develop fundamental understanding of reaction pathways and the capability to predict product distributions. Nearly all of the candidate catalysts were prepared in-house via standard techniques such as impregnation, co-impregnation, or chemical vapor deposition. Supports were usually purchased, but in some cases coprecipitation was used to simultaneously create the support and active component, which can be advantageous for strong active component-support interactions and for achieving high active component dispersion. In-house synthesis also allowed for studies of the effects on catalyst activity and selectivity of such factors as support porosity, calcination temperature, and reduction/activation conditions. Depending on the physical characteristics of the molecule, catalyst activity measurements were carried out in tubular flow reactors (for vapor phase reactions) or stirred tank reactors (for liquid phase reactions) over a wide range of pressures and temperatures. Reactant and product concentrations were measured using gas chromatography (both on-line and off-line, with TCD, FID, and/or mass spectrometric detection). For promising catalysts, detailed physicochemical characterization was carried out using FTIR, Raman, XPS, and XRD spectroscopies (all available in our laboratories) and TEM spectroscopy (available at OU). Additional methods included temperature programmed techniques (TPD, TPO) and surface area measurements by nitrogen adsorption techniques.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

In-Canal Assay of High Specific Activity 60 Co at the Advanced Test Reactor

Production of high specific activity (HSA) 60 Co has recently resumed at the Advanced Test Reactor (ATR) at the Idaho National Laboratory. The technical steps of performing in-canal assay of HSA 60 Co targets at the ATR are described herein. The HSA targets are assayed on a regular basis, between cycles at the ATR, in order to assess the progress of activation. The targets are also assayed at the conclusion of the irradiation in order to provide activity estimates for the distributor and to be used for safety and shipping evaluations. These target assay activities must take place in the ATR canal to provide sufficient radiological shielding. Furthermore, a specialized assay fixture is used in conjunction with custom 60 Co radiation standards to assess the irradiated target. The specific activity of each irradiated target was determined by first measuring the 60 Co standards to determine a sensitivity factor. Then, each irradiated target was measured, and the integrated measurement was multiplied by the sensitivity factor to determine the total activity of the irradiated target. Finally, a correction factor was determined to adjust the reported activity, accounting for differences in the physical geometries of the standards and targets. These methods were used to accurately assay the gross activity of nine irradiated 60 Co targets, two of which were delivered to the distributor where confirmatory hot-cell assays were performed verifying the accuracy of the in-canal assay method.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

Spent fuel measurements. passive neutron albedo reactivity (PNAR) and photon signatures

The International Atomic Energy Agency’s (IAEA) safeguards technical objective is the timely detection of a diversion of a significant quantity of nuclear material from peaceful activities to the manufacture of nuclear weapons or of other nuclear explosive devices or for purposes unknown, and deterrence of such diversion by the risk of early detection. An important IAEA task towards meeting this objective is the ability to accurately and reliably measure spent nuclear fuel (SNF) to verify reactor operating parameters and verify that the fuel has not been removed from reactors or SNF storage facilities. This dissertation analyzes a method to improve the state-of-the-art of nuclear material safeguards measurements using two combined measurement techniques: passive neutron albedo reactivity (PNAR) and passive spectral photon measurements.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2022

Perspectives on the integration between first-principles and data-driven modeling

Efficiently embedding and/or integrating mechanistic information with data-driven models is essential if it is desired to simultaneously take advantage of both engineering principles and data-science. Further the opportunity for hybridization occurs in many scenarios, such as the development of a faster model of an accurate high-fidelity computer model; the correction of a mechanistic model that does not fully-capture the physical phenomena of the system; or the integration of a data-driven component approximating an unknown correlation within a mechanistic model. At the same time, different techniques have been proposed and applied in different literatures to achieve this hybridization, such as hybrid modeling, physics-informed Machine Learning (ML) and model calibration. In this paper we review the methods, challenges, applications and algorithms of these three research areas and discuss them in the context of the different hybridization scenarios. Moreover, we provide a comprehensive comparison of the hybridization techniques with respect to their differences and similarities, as well as advantages and limitations and future perspectives. Finally, we apply and illustrate hybrid modeling, physics-informed ML and model calibration via a chemical reactor case study.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Uranium–Molybdenum Alloy Critical Experiments for the Design of the Health Physics Research Reactor

Clean critical experiments with a uranium-molybdenum alloy (average of 10.1616 wt. % Mo with a density of 17.08 g/cm 3 ) were performed at the Oak Ridge Critical Experiments Facility in 1961 to support the design of the Health Physics Research Reactor (HPRR). The HPRR was similar to the Godiva burst reactor at Los Alamos National Laboratory and was designed to produce 50 microseconds burst of 10 17 fission pulses of radiation for dosimetry measurements, initially in support of the determination of the doses from the nuclear detonations in Japan during World War II. These experiments reported here were used to verify the calculational methods used to design the HPRR. These delayed critical measurements were:1) a solid unreflected and unmoderated 8-in.-dimeter U-Mo cylinder, 2) an unmoderated and unreflected annulus with 8-in.-outside diameter, 2-in.-inside diameter cylinder with a central void, 3) an unmoderated and unreflected annulus with 8-in.-outside diameter, 2-in.-inside diameter cylinder with a central void filled with stainless steel, 4) Same as 3) but with 3-in-thick Plexiglas reflector on top with and without cadmium between the reflector and the U-Mo alloy assembly with steel in the center, and 5) an unmoderated and unreflected annulus which was a modification of the second but with the lower 5 inches of the central hole enlarged to 3.5 in. with various reflector conditions. The reflector conditions were: 1-in.-thick Plexiglas on all outer surfaces-void in the center; 1-in.-thick Plexiglas on all outer surfaces-Plexiglas in the center; 2-in.- thick Plexiglas on radial surface-void in the center; 6-in.-thick Plexiglas on the bottom only-Plexiglas in the center; and 6-in.-thick Plexiglas on bottom, 1-in.-thick on top and on the lower 8.25-cm.-section of the radial surface-void in the center. For some of these reflector conditions 0.025-cm.thick cadmium was located between the reflector and the U-Mo alloy. The uranium contained 93.17 wt. % 235 U. Reflection was a safety concern for this unmoderated and unreflected reactor and reduction of reflection effects was also investigated by insertion of neutron absorber around the U-Mo alloy. The stainless steel 304 contained 18% nickel and 8% chromium and the rest iron. The reflector material was a methacrylate plastic (Plexiglas) containing 5.8 x 10 22 atoms/cm 3 of hydrogen and 3.6 x 10 22 atoms/cm 3 of carbon with a density of 1.20 g/cm 3 . The purpose of this report is to document the experimental information for the measurements performed so that at a later date researchers could perform the required uncertainty and calculational analyses and documentation to use these data for an International Nuclear Criticality Safety Benchmark Program (ICSBEP) or a EURATON Nuclear Energy Agency (NEA) benchmark. The data from the experiments described should be acceptable for use as criticality safety benchmark experiments for the ICSBEP and the NEA nuclear criticality safety benchmark program, once the uncertainty analysis is completed. Based on previous ICSBEP benchmarks with this enriched uranium metal at ORCEF, the uncertainties in k eff could be as low as ±0.0002 for some configurations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

MATCHED-INDEX-OF-REFRACTION FLOW FACILITY FOR FUNDAMENTAL AND APPLIED RESEARCH

Significant challenges face reactor designers with regard to thermal hydraulic design and associated modeling for advanced reactor concepts. Computational thermal hydraulic codes solve only a piece of the core. There is a need for a whole core dynamics system code with local resolution to investigate and understand flow behavior with all the relevant physics and thermo-mechanics. The matched index of refraction (MIR) flow facility at Idaho National Laboratory (INL) has a unique capability to contribute to the development of validated computational fluid dynamics (CFD) codes through the use of state-of-the-art optical measurement techniques, such as Laser Doppler Velocimetry (LDV) and Particle Image Velocimetry (PIV). PIV is a non-intrusive velocity measurement technique that tracks flow by imaging the movement of small tracer particles within a fluid. At the heart of a PIV calculation is the cross correlation algorithm, which is used to estimate the displacement of particles in some small part of the image over the time span between two images. Generally, the displacement is indicated by the location of the largest peak. To quantify these measurements accurately, sophisticated processing algorithms correlate the locations of particles within the image to estimate the velocity (Ref. 1). Prior to use with reactor deign, the CFD codes have to be experimentally validated, which requires rigorous experimental measurements to produce high quality, multi-dimensional flow field data with error quantification methodologies. Computational thermal hydraulic codes solve only a piece of the core. There is a need for a whole core dynamics system code with local resolution to investigate and understand flow behavior with all the relevant physics and thermo-mechanics. Computational techniques with supporting test data may be needed to address the heat transfer from the fuel to the coolant during the transition from turbulent to laminar flow, including the possibility of an early laminarization of the flow (Refs. 2 and 3) (laminarization is caused when the coolant velocity is theoretically in the turbulent regime, but the heat transfer properties are indicative of the coolant velocity being in the laminar regime). Such studies are complicated enough that computational fluid dynamics (CFD) models may not converge to the same conclusion. Thus, experimentally scaled thermal hydraulic data with uncertainties should be developed to support modeling and simulation for verification and validation activities. The fluid/solid index of refraction matching technique allows optical access in and around geometries that would otherwise be impossible while the large test section of the INL system provides better spatial and temporal resolution than comparable facilities. Benchmark data for assessing computational fluid dynamics can be acquired for external flows, internal flows, and coupled internal/external flows for better understanding of physical phenomena of interest. The core objective of this study is to describe MIR and its capabilities, and mention current development areas for uncertainty quantification, mainly the uncertainty surface method and cross-correlation method. Using these methods, it is anticipated to establish a suitable approach to quantify PIV uncertainty for experiments performed in the MIR.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

The thermal conductivity of mixed fuel U x Pu 1-x O 2 : molecular dynamics simulations

Mixed oxides (MOX), in the context of nuclear fuels, are a mixture of the oxides of heavy actinide elements such as uranium, plutonium and thorium. The interest in the UO 2 -PuO 2 system arises from the fact that these oxides are used both in fast breeder reactors (FBRs) as well as in pressurized water reactors (PWRs). The thermal conductivity of UO 2 fuel is an important material property that affects fuel performance since it is the key parameter determining the temperature distribution in the fuel, thus governing, e.g., dimensional changes due to thermal expansion, fission gas release rates, etc. For this reason it is important to understand the thermal conductivity of MOX fuel and how it differs from UO 2 . Here, molecular dynamics (MD) simulations are carried out to determine quantitatively, the effect of mixing on the thermal conductivity of U x Pu 1-x O 2 , as a function of PuO 2 concentrations, for a range of temperatures, 300 – 1500 K. The results will be used to develop enhanced continuum thermal conductivity models for MARMOT and BISON by INL. These models express the thermal conductivity as a function of microstructure state-variables, thus enabling thermal conductivity models with closer connection to the physical state of the fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Thermal Conductivity of Metallic Uranium

This project has developed a modeling and simulation approaches to predict the thermal conductivity of metallic fuels and their alloys. We focus on two methods. The first method has been developed by the team at the University of Wisconsin Madison. They developed a practical and general modeling approach for thermal conductivity of metals and metal alloys that integrates ab-initio and semi-empirical physics-based models to maximize the strengths of both techniques. The second method has been developed by the team at Virginia Tech. This approach consists of a determining the thermal conductivity using only ab-initio methods without any fitting parameters. Both methods were complementary. The models incorporated both phonon and electron contributions. Good agreement with experimental data over a wide temperature range were found. The models also provided insight into the different physical factors that govern the thermal conductivity under different temperatures. The models were general enough to incorporate more complex effects like additional alloying species, defects, transmutation products and noble gas bubbles to predict the behavior of complex metallic alloys like U-alloy fuel systems under burnup. 3 Introduction Thermal conductivity is an important thermal physical property affecting the performance and efficiency of metallic fuels. Some experimental measurement of thermal conductivity and its correlation with composition and temperature from empirical fitting are available for U, Zr and their alloys with Pu and other minor actinides. However, as reviewed in by Kim, Cho and Sohn, due to the difficulty in doing experiments on actinide materials, thermal conductivities of metallic fuels have only been measured at limited alloy compositions and temperatures, some of them even being negative and unphysical. Furthermore, the correlations developed so far are empirical in nature and may not be accurate when used for prediction at conditions far from those used in the original fitting. Moreover, as fuels burn up in the reactor and fission products are built up, thermal conductivity is also significantly changed. Unfortunately, fundamental understanding of the effect of fission products is also currently lacking. In this project, we probe thermal conductivity of metallic fuels with ab initio calculations, a theoretical tool with the potential to yield better accuracy and predictive power than empirical fitting. This work will both complement experimental data by determining thermal conductivity in wider composition and temperature ranges than is available experimentally, and also develop mechanistic understanding to guide better design of metallic fuels in the future. So far, we focused on α-U perfect crystal, the ground-state phase of U metal. We focus on two methods. The first method has been developed by the team at the University of Wisconsin Madison. They developed a practical and general modeling approach for thermal conductivity of metals and metal alloys that integrates ab-initio and semi-empirical physics-based models to maximize the strengths of both techniques. The second method has been developed by the team at Virginia Tech. This approach consists of a determining the thermal conductivity using only ab-initio methods without any fitting parameters. Both methods were complementary and very helpful to understand the physics behind the thermal conductivity in metallic uranium and other materials with similar characteristics. In Section I, the combined model developed at UWM is explained. In Section II, the ab-initio method developed at VT is described along with the uranium pseudo-potential and its validation. Section III is devoted to the work done by Jianguo Yu at INL. Finally, we will present the performance of the project in terms of milestones, publications, and presentations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

SCALE 6.2 Lattice Physics Performance Assessment

The US Nuclear Regulatory Commission relies on the lattice physics analysis capabilities of the SCALE code system to perform confirmatory licensing analyses. Either SCALE lattice physics code—TRITON/NEWT or Polaris—can be used to generate cross section data used by the PARCS nodal core simulator for full-core neutronics calculations. This report presents an assessment of the accuracy of SCALE lattice physics codes for preparation of lattice physics data that are used to support simulator codes such as the NRC’s PARCS, for UO 2 -mixed oxide (MOX)/Zr fueled light water reactor (LWR) analyses. Due to the nature of lattice physics calculations, critical reactor experiment benchmarks cannot be modeled in explicit detail in a lattice physics code. However, this limitation does not mean that these measurement data are not usable for lattice physics studies. Therefore, either geometry approximations or axial buckling must be implemented to determine the critical water height. These modeling limitations have led to development of a three-phase assessment strategy. In the first phase, selected critical experiment benchmarks are modeled using the SCALE 3D continuous-energy (CE) Monte Carlo (MC) code KENO, which is the most rigorous neutron transport method available in SCALE, with no approximations in the spatial, angle, or energy treatments. Biases and statistical uncertainties in quantities of interest such as k eff and pin power distributions are determined by comparing CE KENO results to experimental data. This first phase in the accuracy assessment is to demonstrate and establish the use of CE KENO as a reference solution for the second phase. In the second phase, 14 numerical test suites are used to compare SCALE lattice physics calculations with CE KENO as a reference solution. In the third phase, SCALE lattice physics depletion calculations are performed, and the spent fuel isotopic results are compared with available radiochemical assay measurements. Isotopic measurement comparisons provide quantitative assessment of isotopic density distribution predictions with the depletion models in SCALE lattice physics codes. This report documents results for all test suites. The assessment was performed using standard production techniques unless otherwise noted. Both TRITON/NEWT and Polaris exhibited acceptable accuracy for most test cases. For the few test cases in which acceptable accuracy criteria were not met, further code and data development are planned. The computer codes used in this assessment are as follows: SCALE 6.2 rev19189 (pre-release of 6.2.1) was used for CE KENO, TRITON/NEWT and Polaris (PWR only) calculations; SCALE 6.2.2 was used for Polaris calculations to address several updates in support of support boiling water reactor (BWR) lattice geometries. For this work, there is no difference between SCALE 6.2.1 and SCALE 6.2.2 beyond the inclusion of the Polaris BWR analysis capability; The ENDF/B-VII.1 continuous energy and 252 group libraries deployed in SCALE 6.2 were generated with AMPX 6.2 that is distributed with SCALE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Uncertainty and Sensitivity Analysis of TREAT Transient Test no. 15

The U.S. Department of Energy (DOE) is preparing for the resumption of transient testing using the Transient Reactor Test Facility (TREAT) located at the Idaho National Laboratory (INL). In order to support the analysis of TREAT operation and experiment planning, the MAMMOTH reactor physics team has recently begun development of uncertainty and sensitivity analysis for the TREAT core simulated via MAMMOTH. MAMMOTH is a high-fidelity three-dimensional transient multi-physics code developed at INL under the DOE-NE Advanced Modeling and Simulation program (NEAMS). This research is part of the development presenting some initial uncertainty and sensitivity analysis of TREAT Transient Test no. 15. In this summary, we employ the Sobol sampling strategy implemented in RAVEN to efficiently perform the sensitivity and uncertainty analysis. The idea is using Sobol functional decomposition with generalized polynomial chaos (gPC) expansions to estimate the responses of interest. The gPC coefficients are then computed via the sparse grid spectral projection method. In addition, this strategy includes evaluation of the second or higher order sensitivity indices, allowing the study of interactions between several input variables. The principal focus of this analysis is to quantify the uncertainties in the total power and maximum fuel temperatures with respect to uncertainties in the initial transient conditions, such as rated power, fuel temperature, graphite fuel density and control rod boron concentration. In addition, both the first order and second order time-dependent sensitivity indices are computed and analyzed. The results show that the interactions between initial control rod boron concentration and the rated power will contribute at most 1% to the total uncertainty of the power peak.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗