Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Fuel Performance Modeling”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Validation of MFUEL Metal Fuel Performance Models of SAS4A/SASSYS-1

The fuel characterization models of SAS4A/SASSYS-1 (SAS) have recently been extended to include a new U-Pu-Zr metal fuel model, MFUEL. MFUEL is equipped with mechanistic physics based models to predict the pre-transient characterization and transient response of metal fuel, with emphasis on fuel melting, cladding failure, and the metal fuel’s impact on core reactivity. The MFUEL model will be available in the full version of SAS4A/SASSYS-1 5.7, which is scheduled to be released in June 2023. Fast reactor fuel pins that operated in EBR-II and FFTF with low smear density U-Zr and U-Pu-Zr metal fuels and irradiation resistant ferritic-martensitic cladding showed significant advantages in achieving high burnups and assuring inherent safety characteristics during anticipated transients, design basis events, and beyond design basis events. For safety analysis, a fuel performance model must be able to predict (1) Fuel pin mechanics and compositional and dimensional changes, (2) Clad failure, and (3) Fuel pin thermal resistance. Achieving these high level goals accurately is strongly related to the model performance of individual physical processes taking place within a fuel pin during its lifetime. Metal fuels typically operate above the mid-point temperature of melting during normal, as well as off-normal, conditions. At these elevated temperatures, the availability of thermal activation provides a driving force for various diffusional processes leading to complex phase transformations, micro-structure evolution, significant amounts of fuel swelling, interconnected porosity formation, excessive amounts of fission gas release, and fuel clad chemical interactions. Clad failure in fast reactors primarily occurs as a result of creep rupture augmented by clad wastage formation. The reaction is driven by thermal creep induced dislocation motion, grain boundary cavity nucleation, growth and breakup of grain boundaries. The high level complexity and limited available data requires introducing physics-based modeling approaches to gain extrapolation ability and sensitivity with respect to various conditions. The objective of this report is to perform validation of MFUEL using the experimental data for (1) Normal operation EBR-II fuel behavior, (2) Normal operation PHENIX fuel behavior, (3) HT9 Pressure tube ramp-and-hold creep rupture tests, (4) Whole Pin Furnace (WPF) creep strain, creep rupture and eutectic tests, (5) Fuel Behavior Test Apparatus (FBTA) eutectic tests, and (6) TREAT M5-7 OverPower tests up to clad failure. Section-2 includes a brief description of the MFUEL models. A detailed description of the MFUEL physics-based, semi-empirical models will be presented in the SAS V 5.7 theory manual. Section-3, 4, and 5 describes the validation effort for the pre-transient irradiation, furnace transients, and TREAT M-Series transients, respectively.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

MELCOR Input Model for Spent Fuel Transportation and Storage Canister

This report documents the progress and current results of the MELCOR spent fuel cask input model. The MELCOR model is being developed to investigate aerosol transport and deposition given the aerosol physical phenomena models within MELCOR. To perform the analyses, a general portrayal of the MAGNASTOR® cask system has been employed; however, this system was selected based on available information to provide a reasonable representation of a spent fuel cask. The analytical results are not intended to characterize the performance of the MAGNASTOR® cask. Instead, the provided results are intended to enhance our general understanding of the aerosol behavior within casks and the validity of current models. The current model efforts are being performed to investigate hypothetical UO 2 release from failed fuel pins within a spent fuel cask. The existing MELCOR model of the MAGNASTOR® cask system has been adapted to permit future comparative analyses with the GOTHIC representation of the MAGNASTOR® cask. To support this comparison, the PNNL model characteristics that are unrelated to the aerosol modeling were applied to the MELCOR model. These characteristics included improved comparability of the axial fidelity, total spent fuel power, fuel pin axial power profile, and heat losses from cannister. The thermal-hydraulic solutions are improved within the capability of the MELCOR code and will permit better overall agreement with the GOTHIC results. Detailed results are presented on the thermal-hydraulic analysis of the MELCOR cask as well as characterization of UO 2 aerosol dispersion and deposition within the cask.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

TRITIUM TRANSPORT MODELING AT THE PORE SCALE IN CERAMIC BEEDER MATERIALS USING TMAP8

Fusion reactors depend on the blanket material to breed and release tritium at the same rate or faster than it is consumed by the fusion reaction. Cellular ceramic breeders are dense materials that can maintain a high tritium breeding ratio while promoting tritium release thanks to highly connected pores. Assessing the tritium breeding capabilities of these materials requires a combination of extensive experimental and modeling efforts. In this work, we develop and calibrate a multiphysics model of tritium transport. This novel model accounts for ceramic and pore diffusion, trapping and detrapping, and several surface reactions at the pore surface. We perform a sensitivity analysis and calibrate the model by comparing its predictions against experimental measurements of deuterium absorption. The calibrated model is then used to model tritium absorption in samples with different pore microstructures to investigate the effect of pore interconnectivity on tritium absorption. The model is part of the development of the multiscale, multiphysics framework for tritium transport (i.e., the Tritium Migration Analysis Program [TMAP8]), which is itself built on top of the finite element multiphysics framework MOOSE. This study demonstrates some of TMAP8’s capabilities, and is the first step toward assessing the tritium breeding capabilities of ceramic breeder material designs.

08 HYDROGEN↗
DOE OSTI2021

Uranium Bed Design Parameters for Tritium Plants Supporting Fusion Reactors

Uranium hydride is commonly used to store hydrogen or its isotopes in a solid state. The Self-Assaying Tritium Accountancy and Containment Unit for ITER (STACI) is a 5.2-kg bed of depleted uranium (dU) capable of holding up to 33 mol of hydrogen or its isotopes. Here we summarize the data analysis of past experimental campaigns with STACI, with the aim of describing the kinetics and thermodynamics of the hydriding process. Computed tomography imaging was performed on STACI both before and after its experimental campaign, and a high degree of swelling was observed in the dU. Literature on studies in regard to the swelling of large (multikilogram) quantities of uranium hydride for storage applications was not identified during this study. Data from the experimental campaign, as well as data on the formation reaction, are presented. The authors hope to create an analytical model of STACI based on these data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Synthetic Data Generation Using Machine Learning

Robust machine learning techniques for image analysis require a substantial amount of data to yield confident results. In the nuclear domain, data scarcity is a substantial challenge because there are so few facilities worldwide. This research focuses on being able alleviate the data scarcity problem by generating synthetic data to bridge the gap between large and small datasets. This work achieves that goal using a Generative Adversarial Network (GAN) architectural approach, by training a model on real-world data and expands that small dataset through synthetic data amendments. Model performance is impacted by the size of the real-world dataset and the number of training epochs utilized. This means that 1) It is important to develop your GAN to be optimized with the specific data type, and 2) approaches taken when training the GAN should be specialized to encompass important aspects of the dataset that it is generating. By taking a step to improve dataset sizes in this way, the gap between models trained by parties with significant amount of data and those without access to large data, closes, allowing for robust analyses of satellite imagery for nuclear domain applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Assessment of MOOSE-Based Tools for Calculating Radial Core Expansion

Radial core expansion in liquid-metal cooled fast reactor systems is a well-known phenomenon that produces strong reactivity feedback effects. An inherently safe reactor design takes advantage of negative reactivity feedback in accident scenarios by utilizing a core restraint system which produces a bowed shape that allows for radial expansion of the fuel regions. Detailed modeling and simulation of core radial expansion itself as well as subsequent reactivity feedback is a challenging task involving contact of many fuel assembly elements and physics feedback from neutronics, thermal hydraulics, and thermal mechanical response. A variety of physics codes have been developed to model aspects of radial core expansion but in general invoke geometrical or physics approximations. No code system currently exists which tightly and robustly couples these physics with enough detail to fully resolve the complex core radial expansion reactivity feedback effects. The future availability of such a code system is of vital importance to fully understanding the reactivity feedback effects that occur due to radial expansion, and consequently to optimizing the design of the core restraint system. A high-fidelity code will also be used to benchmark existing lower fidelity, faster running models to understand their benefits, limitations, and range of applications. A code development path using MOOSE-based tools is proposed in order to leverage the detailed geometry capabilities and natural tight coupling and robustness of MOOSE-based applications for modeling this complex phenomena. While simulation of the full phenomenon involves several physics, an assessment has been initiated on the capabilities and readiness of the currently available Tensor Mechanics module within MOOSE for calculation of the structural mechanical responses which occur within the reactor core. This report focuses on modeling the force-deformation response which mimics the physics of a duct contact deformation, as well as differential thermal expansion which produces a thermal bowed shaped for the fuel assemblies in a core. Simple examples were initially performed such as simple supported beam bending under load. The complexity of examples was progressively increased to better mimic the duct behavior by including a differential thermal example and inclusion of hexagonal cross-sections in the geometry. Further assessment of the structural mechanical response simulation capability is still required for modeling duct contact interactions and irradiation creep and swelling. Companion thermal hydraulic and neutronics assessments will also be required; these activities are planned for future years. Finally, integration of the multiple physics components through MOOSE is required to predict the core radial expansion and subsequent feedback effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

On the effect of scalar flux weighing of linearly anisotropic scattering matrices in few-group transport calculations

The majority of the codes available for homogenized group constant generation for deterministic transport calculations apply the approximation of scalar flux weighting during energy group condensation of higher-order anisotropic scattering matrices. In this paper, we discuss the effect of scalar flux weighting of the linearly anisotropic scattering matrices in the frame of S P{sub 3} and S{sub 12} calculations performed for a two-dimensional VVER-440 reactor benchmark. To compare group constants generated for two neutron energy groups, an infinite pin cell was homogenized with Serpent 2 and ERANOS ECCO. Serpent 2 applies scalar flux, while ERANOS ECCO performs current weighting of the linearly anisotropic scattering matrices during energy group condensation. For analyzing the effect of the various weighting options, three simple reactor models were built assuming different core sizes using standard rectangular assemblies with 15*15 fuel pins. Diffusion, SP{sub 3} and S{sub 12} calculations were performed for the 3 models using group constants generated with Serpent 2 and ERANOS ECCO. The effect of scalar flux weighting of linearly anisotropic scattering matrices is shown by comparing the decrease in reactivity due to the decreased reactor size, as well as the assembly power distribution to reference results obtained with Serpent 2 Monte Carlo calculations. Neglecting higher than linearly anisotropic scattering and indirect application of diffusion coefficients in higher-order transport calculations is advised if angular flux-moment spectra weighted higher-order scattering matrices cannot be generated. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Summary of GNU Octave Lady Godiva Simulator

A pseudo real-time interactive reactor simulator for the Lady Godiva reactor has been developed using the GNU Octave software package. The simulator tracks the time-dependent power and fuel temperature behavior of the Lady Godiva reactor during sub-prompt critical transients. Users interact with the simulator by controlling reactivity input. This input occurs through users repositioning a control rod as desired during simulator operation. Various reactor parameters are displayed graphically in near “real time,” facilitating user assessment and feedback input. The reactor model consists of a lumped-fuel core, a single control rod, and a simple heat transfer mechanism that allows heat to be rejected from the core. The goal of the simulator is to demonstrate in real-time the system response to reactivity perturbations, the effect of reactivity feedback, and the effect of temperature-dependent equilibrium conditions. The dynamic behavior of the reactor is primarily driven by the relative characteristic time-constants between the delayed neutrons and the overall heat transfer parameter. The simulator will allow users to gain experience performing standard reactor operations, such as establishing steady-state at different fuel temperatures and making power level changes

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Methodology for Sensitivity Analysis of Homogenized Cross-Sections to Instantaneous and Historical Lattice Conditions with Application to AP1000® PWR Lattice

In the two-step method for nuclear reactor simulation, lattice physics calculations are performed to compute homogenized cross-sections for a variety of burnups and lattice configurations. A nodal code is then used to perform full-core analysis using the pre-calculated homogenized cross-sections. One source of uncertainty introduced in this method is that the lattice configuration or depletion conditions typically do not match a pre-calculated one from the lattice physics simulations. Therefore, some interpolation model must be used to estimate the homogenized cross-sections in the nodal code. This current study provides a methodology for sensitivity analysis to quantify the impact of state variables on the homogenized cross-sections. This methodology also allows for analyses of the historical effect that the state variables have on homogenized cross-sections. An application of this methodology on a lattice for the Westinghouse AP1000® reactor is presented where coolant density, fuel temperature, soluble boron concentration, and control rod insertion are the state variables of interest. The effects of considering the instantaneous values of the state variables, historical values of the state variables, and burnup-averaged values of the state variables are analyzed. Using these methods, it was found that a linear model that only considers the instantaneous and burnup-averaged values of state variables can fail to capture some variations in the homogenized cross-sections.

Price, Dean (ORCID:0000000309990111)↗
DOE OSTI2024

Final Report on Predictive Analyses of PRD as Function of Anomalies

This report summarizes the work completed in FY-2024 to analyze the power reactivity decrement (PRD) concepts of the ARC-100 core. The PRD has been traditionally defined by the reactivity change from a hot zero power (HZP) to a particular power. Consequently, the PRD accounts for the core reactivity changes due to increase coolant temperature gradient axially and radially across the core, and increased fuel temperature. The coolant temperature gradient leads to sodium and structure density changes, to radial core expansion from assembly flowering and bowing (due to axial and radial temperature gradients within the assemblies), and to control rod driveline thermal expansion. The fuel temperature increase associated with coolant temperature and power increases leads to Doppler effect and axial thermal expansion. In this work, the normal operating Hot Full Power (HFP) state is the only power level of interest, so analyses focus on the PRD calculated from HZP to HFP. The PRD has been used to assess the reactor safety features asymptotically in unprotected accident scenarios, including the loss of heat sink (LOHS), loss of flow (LOF), and transient overpower (TOP) without scram. The PRD concept relies on the “global” reactivity coefficients A, B, and C that are estimated based on “individual” reactivity effects (Doppler, sodium density, etc.). The objectives of this work are: 1) to improve and verify the methodology used to calculate the ABC coefficients used in the PRD, 2) to assess if the PRD can be used to reliably identify abnormal events. This report fulfills the FY-2024 scope of WBS#1.15.8.3 activity, “ANL0120 – Predictive Analyses of PRD as function of deformation”. The PRD concept is described in Sections 2. Additional effort in refining the methodology for core bowing modeling is performed in Sections 3. Then, two verification exercises are proposed in Section 4 to benchmark these coefficients based on direct neutronic-only calculations and on dynamic core transient simulations. Finally, the PRD approach is assessed for the detection of several unexpected events, such as primary flow perturbation or improper fuel loading, in Section 5.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Microstructurally validated stable and predictable swelling in low-enriched uranium monolithic U-10Mo fuel mini-plates

Qualification of the low-enriched uranium (LEU) monolithic U-10 wt%Mo (U-10Mo) plate-type fuel system requires a demonstration of a stable and predictable fuel swelling behavior over the anticipated operating conditions of the United States high-performance research reactors (USHPRRs) selected for conversion to LEU operation. This will allow each reactor to develop appropriate safety margins that will retain fuel element lifetime coolability. Additionally, the fuel system must maintain performance attributes when fabricated at a commercial scale. The Mini-plate 1 experiment represents the first irradiation test of commercially fabricated miniaturized monolithic LEU U-10Mo fuel plates. Here, the swelling behavior within this experiment was compared against that of historical fuel developmental tests to reveal that the commercially fabricated fuel performed within the current recommended U-10Mo swelling model's predictions. Additionally, the fuel microstructural evolution was evaluated to link initial conditions to subtle variations detected in the swelling response, providing validation and confidence that the fuel system is robust.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology for pin lattice models is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. The user can alternatively provide their own geometry and power profile instead of inheriting it from DIF3D. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and a sufficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

High-Fidelity Modeling of Fuel-To-Coolant Thermomechanical Transport Behaviors Under Transient Conditions

This report summarizes the work completed under NEUP project number 21-24006. The objectives of this project are to advance the high-fidelity modeling capabilities and important phenomena that is important for high-burnup UO 2 and accident tolerant fuels (ATF) during transient conditions. Accurate modeling of the time-dependent phenomena that impact material performance must be used to determine the figures of merit and safety margin. Phenomena such as fuel fragmentation, cladding oxidation, pellet-clad interaction, clad ballooning, and clad rupture are examples that pose challenges to modeling during these transients. This project focused on the development of high-fidelity tightly coupled multiphysics models that can capture the time-dependent material response and associated thermal hydraulic conditions during these events. These models can then be validated against existing separate effects tests and in-pile integral experiments and will be used to model Transient Reactor Test facility (TREAT) loss-of-coolant accidents (LOCA) experiments. To achieve the project objective, we used a combination of NEAMS and NRC codes to model various LOCA test sets for the separate effects and in-pile integral experiments. BlueCRAB tool set, which can accurately predict material response at a sub-fuel pin level, as well as modeling the entire reactor system response to these events. Fuel performance was modeled using BISON (various versions) and FAST (version 1.2.1). BISON and FAST can model on a sub-fuel pin level the fuel performance under transient conditions. Both have simplified thermal hydraulic models that are capable of providing basic coolant boundary conditions. To better capture the thermomechanical interaction between the fuel, clad, and coolant, more sophisticated thermal hydraulic models are necessary.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

SNF Interim Storage Canister Corrosion and Surface Environment Investigations (FY2020 Status Report)

This progress report describes work performed during FY20 at Sandia National Laboratories (SNL) to assess the localized corrosion performance of container/cask materials used in the interim storage of spent nuclear fuel (SNF). Of particular concern is stress corrosion cracking (SCC), by which a through-wall crack could potentially form in a canister outer wall over time intervals that are shorter than possible dry storage times. Work in FY20 further defined our understanding of the potential chemical and physical environment present on canister surfaces, evaluated the relationship between the environment and the resultant corrosion that occurs, and initiated crack growth rate testing under relevant environmental conditions. In FY20, work to define dry storage canister surface environments included several tasks. First, collection of dust deposition specimens from independent spent fuel storage installation (ISFSI) site locations helped to establish a more complete understanding of the potential chemical environment formed on the canister. Second, the predicted evolution of canister surface relative humidity RH) values was estimated using ISFSI site weather data and the horizontal canister thermal model used by the SNL probabilistic SCC model. These calculations determined that for typical ISFSI weather conditions, seasalt deliquescence to produce MgCl2-rich brines could occur in less than 20 years at the coolest locations on the canister surface, and, even after nearly 300 years, conditions for NaCl deliquescence (75% RH) are not reached. This work illustrates the importance of understanding the stability of MgCl2-rich brines on the heated canister surface, and the potential impact of brine composition on corrosion processes, including pitting and stress corrosion cracking. In an additional study, the description of the canister surface environment was refined in order to define more realistic corrosion testing environments including diurnal cycles, soluble salt chemistries, and inert mineral particles. The potential impacts of these phenomena on canister corrosion are being evaluated experimentally. Finally, work over the past few years to evaluate the stability of magnesium chloride brines continued in FY20. MgCl2 degassing experiments were carried out, confirming that MgCl2 brines slowly degas HCl on heated surfaces, converting to less deliquescent magnesium hydroxychloride phases and potentially leading to brine dryout.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Validation of Thermohydraulic Simulations using RELAP for Critical Dual Purpose Canisters

Simulations using the RELAP5-3D systems analysis code and STAR-CCM+, a high-fidelity CFD code, were performed on a model of a real-scale MPC-32 dual purpose canister undergoing a criticality event. These simulations were performed to provide code-to-code verification of the natural convection capabilities in RELAP5. The STAR-CCM+ simulations revealed various natural convection flow patterns that aided in the cooling of fission heat produced in the fuel rods in the canister. Negligible spatial temperature variations were observed across the canister in the STAR-CCM+ results, supporting the use of RELAP5 in predicting the sub-cooled regime. The surface-averaged rod temperature predicted by RELAP5 compared well with STAR-CCM+, verifying the predictive capability of RELAP5 for sub-cooled conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Do not forget the electrons: Extending moderately-sized nuclear networks for multidimensional hydrodynamic codes

Context.Nuclear networks are widely used coupled with hydrodynamical simulations of explosive scenarios to account for the change of nuclear species and energy generation rate due to nuclear reactions. In this way, there is a feedback mechanism between the hydrodynamical state and the nuclear processes. Unfortunately, the timescale of nuclear reactions is orders of magnitude smaller than the dynamical timescale that drives hydrodynamical simulations. Therefore, these nuclear networks are usually very small, reduced in most cases to a dozen elements, especially when simulations are carried out in more than one dimension. Aims.We present here an extended nuclear network, with 90 species, designed for being coupled with hydrodynamic simulations, which includes neutrons, protons, electrons, positrons, and the corresponding neutrino and anti-neutrino emission. This network is also coupled with temperature, making it extremely robust and, together with its size, unique of its kind. The inclusion of electron captures on free protons makes the network very appropriate for multidimensional studies of Type Ia supernova explosions, especially when the exploding object is a massive white dwarf. Methods.We perform several tests that are relevant to simulate explosive scenarios, such as Type Ia supernovae and core-collapse supernovae. We compare the results of the 90 nuclei network with a standardα-chain network with 14 elements to evaluate the differences in the energy generation rate. We also evaluate the relevance of including the electrons in the network in terms of generated yields and how it affects the pressure of a degenerate fluid such as that of white dwarfs. The results obtained with the 90-nuclei network have been verified with a much larger 2000-nuclei network built from REACLIB (WinNet), in terms of nuclear energy generation rate, pressure, and produced yields. Results.The results obtained with the proposed medium-sized network compare fairly well, to a few percent, with those computed withWinNetin scenarios reproducing the gross physical conditions of current Type Ia supernova explosion models. In those cases where the carbon and oxygen fuel ignites at high density, the high-temperature plateau typical of the nuclear statistical equilibrium regime is well defined and stable, allowing large integration time steps. We show that the inclusion of electron captures on free protons substantially improves the estimation of the electron fraction of the mixture. Therefore, the pressure is better determined than in networks where electron captures are excluded, which will ultimately lead to more reliable hydrodynamic models. Explosive combustion of helium at low density, occurring near the surface layer of a white dwarf, is also better described with the proposed network, which gives nuclear energy generation rates much closer toWinNetthan typical reduced alpha networks. Conclusions.A nuclear network withN= 90 species, including electrons, aimed at multidimensional calculations of supernova explosions is described and verified. The proposed network is suitable for the study of Type Ia supernova explosions because it provides better values of pressure and electron abundance than other existing networks with smaller or even a similar size but without including electron capture processes.

Astronomy & Astrophysics↗
DOE OSTI2024

Uranium measurements in the field using high-resolution cadmium zinc telluride detectors

A new generation of cadmium zinc telluride (CZT) detectors has become available and is being evaluated by the International Atomic Energy Agency (IAEA) for safeguards verifications in the field. The new CZT detector, model M400, is a room temperature spectrometer manufactured by H3D, Inc. The M400 demonstrates superior energy resolution, effective isotope identification capabilities, and convenient usability features when tested in a controlled laboratory environment. These characteristics define the M400 as a potential platform for IAEA field detection applications, which could become suitable for nuclear material characterization (e.g., enrichment verification) and nuclear safeguards missions. The capabilities of gamma spectrometry codes including Fixed energy, Response function Analysis with Multiple efficiencies (FRAM) from Los Alamos National Laboratory, CZT for Uranium (CZTU) from Lawrence Livermore National Laboratory, and Gamma Detector Response and Analysis Software (GADRAS) from Sandia National Laboratories were adapted for M400 spectra, and the performance of the codes has been validated. This was reported in a prior work. To further validate the performance of the high-energy resolution CZT detector and the isotopic analysis codes, a field measurement campaign consisting of uranium hexafluoride (UF 6 ) cylinder measurements was conducted at a fuel fabrication facility. A total of 34 Type 30B cylinders containing UF6 were measured using three different M400 CZT detectors. Each detector was outfitted with a custom rectangular collimator and shield made out T-Flex®, a tungsten-impregnated polymer. Measurements were performed at three different locations of the cylinder, ensuring that the measurement geometry satisfied the infinite thickness criterion. The spectra from the M400 CZT were analyzed using the code General Enrichment Meter (GEM). For analyzing the gamma-ray spectra from UF 6 cylinder, the GEM code is the appropriate tool since it relies only on the gamma-ray emissions from 235 U and not from other isotopes. Results from the spectral analysis were compared with the known abundance of 235 U in the cylinders, as well as with the International Target Values 2020 (ITV2020). The suitability of the different underlying techniques used by the various codes for UF 6 analysis is discussed. The challenges of measuring UF 6 contained in cylinders and mitigation strategies are highlighted.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Construction of a Code Verification Matrix for Heat Conduction With Finite Element Code Applications

When establishing the pedigree of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model. During this process, numerical results on various meshes are systematically compared to a reference analytic solution. The selection of analytic solutions can be a laborious process, as it is difficult to establish adequate code confidence without performing redundant work. In this report we address this issue by applying a physics-based process that establishes a set of reference problems. In this process, code simulation options are categorized and systematically tested, which ensures that gaps in testing are easily identified and addressed. The resulting problems are primarily intended for code verification analysis but may also be useful for comparison to other simulation codes, troubleshooting activities, or training exercises. The process is used to select fifteen code verification problems relevant for the one-dimensional steady-state heat conduction equation. These problems are applicable to a wide variety of simulation tools, but, in this work, a demonstration is performed using the finite element-based nuclear fuel performance code BISON. Convergence to the analytic solution at the theoretical rate is quantified for a selection of the problems, which establishes a baseline pedigree for the code. Not only can this standard set of conduction solutions be used for verification of other codes, but also the physics-based process for selecting problems can be utilized to quantify and expand testing for any simulation tool.

42 ENGINEERING↗