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 523 records · Page 29

Integrating Advanced Modeling and Accelerated Testing for a Modernized Fuel Qualification Paradigm

With the increasing interest in sodium fast reactor technology, as seen by applications to the U.S. Nuclear Regulatory Commission for the OKLO Aurora plant, fuel testing for the TerraPower Traveling Wave Reactor, and the impending construction and startup of the versatile test reactor (VTR), a modernized, accelerated approach to fuel qualification is needed. To guide this effort, a Phenomena Identification Ranking Table–styled analysis was performed for a U-Pu-Zr sodium-free annular fuel system. This analysis evaluated a series of fuel design properties and parameters against their contributions to key fuel performance phenomena. The resulting priority parameters were then reviewed against existing modeling and experimental capabilities to support investigation of the highest-priority parameters. A pathway for qualification was then established using highthroughput, high-volume experiments from MiniFuel and FAST in parallel with advanced physics-based model development. This effort outlines how the first stages of qualification can be reduced from the typical 20+-year development cycle to 5 to 7 years by deploying accelerated irradiation testing platforms. As with any accelerated test, these methods are prototypic in some aspects and less so in others; however, by coupling with advanced fuel performance modeling and simulation capabilities, the larger space of irradiation parameters and material response provided offers advantages for the validation of physics-based models supporting the deployment of novel fuel designs. As a test case, this paper utilizes a proposed Mark II fuel system for the upcoming VTR. Thus, an accelerated qualification method can be tested for the development of MARK II driver fuel so that by the time of VTR startup, lead test assemblies for a Mark II fuel can be initiated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Making sense of uncertain nuclear data

At its core, and in an imperfect world, the scientific method relies heavily on the concept of uncertain knowledge. Understanding the nature and magnitude of our lack of knowledge on nuclear data has important ramifications on a wide range of nuclear technologies that encompass defense, energy, medicine, astrophysics, among others. Over the span of a long and fruitful carrier, Dr. Massimo Salvatores has made significant contributions to the fields of nuclear reactor physics, nuclear data, and uncertainty quantification and propagation. Advances in machine learning techniques and computational capabilities are enabling a new level of confidence in our ability to predict key nuclear metrics as well as their associated uncertainties across many applications. Here this short article reviews some of those concepts with an eye towards developing a new paradigm to provide nuclear data libraries of benefit to various applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development and Implementation of CFD-Informed Models for the Advanced Subchannel Code CTF

The research described in this PhD thesis contributes to the development of efficient methods for utilization of high-fidelity models and codes to inform low-fidelity models and codes in the area of nuclear reactor core thermal-hydraulics. The objective is to increase the accuracy of predictions of quantities of interests using high-fidelity CFD models while preserving the efficiency of low-fidelity subchannel core calculations. An original methodology named Physics- based Approach for High-to-Low Model Information has been further developed and tested. The overall physical phenomena and corresponding localized effects, which are introduced by the presence of spacer grids in light water reactor (LWR) cores, are dissected in corresponding four building basic processes, and corresponding models are informed using high-fidelity CFD codes. These models are a spacer grid-directed cross-flow model, a grid-enhanced turbulent mixing model, a heat transfer enhancement model, and a spacer grid pressure loss model. The localized CFD-models are developed and tested using the CFD code STAR-CCM+, and the corresponding global model development and testing in sub-channel formulation is performed in the thermal- hydraulic subchannel code CTF. The improved CTF simulations utilize data-files derived from CFD STAR-CCM+ simulation results covering the spacer grid design desired for inclusion in the CTF calculation. The current implementation of these models is examined and possibilities for improvement and further development are suggested. The validation experimental database is extended by including the OECD/NRC PSBT benchmark data. The outcome is an enhanced accuracy of CTF predictions while preserving the computational efficiency of a low-fidelity subchannel code.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Thermal transport at (001) twist grain boundaries in UO{sub 2}

The thermal conductivity of the nuclear fuels critically affects performance and safety of the nuclear reactor and thus continues to attract interest of engineers and researches. In spite of the large volume of experimental data on the thermal transport properties in UO{sub 2}, new discoveries in recent years include anisotropic thermal conductivity, and unexpectedly large contribution of the optical phonons to thermal transport. On the other hand, advances in computer power allow to simulate directly from atomic level the effect of structural defects such as grain boundaries, dislocations and pores on the thermal conductivity in UO{sub 2}. Unfortunately, such simulations do not provide detailed information about phonon scattering processes, which are important for the complete understanding of the thermal transport and generation of the physics-based thermal transport models in nuclear performance codes. In this work, we apply phonons wave-packet (PWD) molecular dynamics simulation technique in order to elucidate the individual phonons contribution to the thermal transport through the grain boundaries (GB) in UO{sub 2}. In this work we elucidated the effect of the grain boundary atomic structure on the details of the thermal conductance through the interface. Namely, we demonstrated that transmission coefficient as a function of the wave vector is a strong function of the disorder at the interface: for ideally structured GBs it has sharp peaks and pronounced valleys. Such structure disappears upon introduction of the disorder near the interface, and as a result transmission coefficients behave very similarly between GBs of different misorientation, but the same degree of disorder. Kapitza conductance through these GBs is also nearly identical. It is unlikely for experimentally observed GBs to be atomically sharp, and thus Kapitza conductance in polycrystalline UO{sub 2} can be considered misorientation independent with good accuracy, and it does not to be accounted for in physics-based thermal transport model for fuel performance codes. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Reducibility of Reactor Physics Depletion and Decay Calculation

Depletion and decay calculations are integral part of reactor physics calculations, necessary to estimate fuel characteristics following a given exposure to radiation. Given the complexity of reactor design and the unprecedented sophistication of the new high fidelity modeling and simulation tools, attempting to calculate a very accurate estimation of fuel characteristics on a pin-by-pin level, there is an ever-increasing need to perform these calculations in a computationally efficient manner. This summary presents a new enabling methodology to achieve higher computational gains while addressing one of the fundamental challenges of existing methods, that is how to preserve the accuracy of the original high fidelity models. This is difficult with current approaches relying on look-up tables or the elimination of non-influential isotopes. To achieve that we employ a reduced order modeling approach which reduces the effective dimensionality of the solution space, represented by the number of tracked nuclides. This approach allows the analyst to devise a user-defined level of reduction supported by rigorous error bounds that quantify the accuracy of the reduced model. This summary explores this new methodology by calculating the effective dimensionality for representative BWR depletion problems that are required to meet user-defined accuracy for the nuclides concentrations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Recent Testing of ENDF/B-VIII.0 at ORNL [Slides]

This presentation provides an overview of comparisons of ENDF/B-VII.1 and ENDF/B-VIII.0 and uses VALID suite to test >600 cases with ENDF/B-VII.1 and ENDF/B-VIII.0 in KENO V.a/KENO-VI, although a limited number of cases have been tested with TSUNAMI. It provides a results summary of KENO V.a/KENO-VI with a pooled results summary. The presentation also discusses changes for LCT experiments, MCT experiments, MST experiments, and PST experiments. A few words about graphite are also provided and the presentation states that graphite performance investigated across ENDF/B-VII.0, ENDF/B-VII.1, and ENDF/B-VIII.0 for MSR and HTGR models. A brief summary of the results is presented to make sure the Validation Committee is aware of the work. Also presented is an Eigenvalue comparison with a graphite moderated system, and comparison between an HTR-10 pebble model ENDF/B-VII.1 vs. ENDF/B-VIII.0, which provides different graphite data. In conclusion, ORNL has recently completed a comparison of ENDF/B-VIII.0 and ENDF/B-VII.1 using the VALID library. Work is also being performed for advanced reactor applications. Other efforts on-going and they are also trying to investigate more intermediate spectrum cases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

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↗

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↗

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↗

Nonlinear rescaling of control values simplifies fuzzy control

Traditional control theory is well-developed mainly for linear control situations. In non-linear cases there is no general method of generating a good control, so we have to rely on the ability of the experts (operators) to control them. If we want to automate their control, we must acquire their knowledge and translate it into a precise control strategy. The experts' knowledge is usually represented in non-numeric terms, namely, in terms of uncertain statements of the type 'if the obstacle is straight ahead, the distance to it is small, and the velocity of the car is medium, press the brakes hard'. Fuzzy control is a methodology that translates such statements into precise formulas for control. The necessary first step of this strategy consists of assigning membership functions to all the terms that the expert uses in his rules (in our sample phrase these words are 'small', 'medium', and 'hard'). The appropriate choice of a membership function can drastically improve the quality of a fuzzy control. In the simplest cases, we can take the functions whose domains have equally spaced endpoints. Because of that, many software packages for fuzzy control are based on this choice of membership functions. This choice is not very efficient in more complicated cases. Therefore, methods have been developed that use neural networks or generic algorithms to 'tune' membership functions. But this tuning takes lots of time (for example, several thousands iterations are typical for neural networks). In some cases there are evident physical reasons why equally space domains do not work: e.g., if the control variable u is always positive (i.e., if we control temperature in a reactor), then negative values (that are generated by equal spacing) simply make no sense. In this case it sounds reasonable to choose another scale u' = f(u) to represent u, so that equal spacing will work fine for u'. In the present paper we formulate the problem of finding the best rescaling function, solve this problem, and show (on a real-life example) that after an optimal rescaling, the un-tuned fuzzy control can be as good as the best state-of-art traditional non-linear controls.

Vanlangingham, H.↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗