SEARCH · Engineering Papers
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.
This presentation covers the background, Theory, Simulation Results, Experimental Attempts, Conclusions, and Future Work for a model for simultaneous evaluation of thermal neutron scattering & resolved resonance parameters.
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.
Abstract We simulate a black hole accretion disk system with full-transport general relativistic neutrino radiation magnetohydrodynamics for 1.2 s. This system is likely to form after the merger of two compact objects and is thought to be a robust site of r -process nucleosynthesis. We consider the case of a black hole accretion disk arising from the merger of two neutron stars. Our simulation time coincides with the nucleosynthesis timescale of the r -process (∼1 s). Because these simulations are time-consuming, it is common practice to run for a “short” duration of approximately 0.1–0.3 s. We analyze the nucleosynthetic outflow from this system and compare the results of stopping at 0.12 and 1.2 s. We find that the addition of mass ejected in the longer simulation as well as more favorable thermodynamic conditions from emergent viscous ejecta greatly impacts the nucleosynthetic outcome. We quantify the error in nucleosynthetic outcomes between short and long cuts.
This paper reviews five critical algorithms for BISON’s simulation of nuclear fuel: the burnup model, the coolant channel model, fuel cracking models, mechanical contact, and thermal contact.
Transient simulations of nuclear systems face the computational challenge of resolving both space and time during reactivity changes. A common strategy for tackling this issue is to split the neutron flux into shape and amplitude functions. This split can be solved with high- order/low-order methods. While this multi-fidelity approach has traditionally been reserved for deterministic methods, it is also possible to implement in Monte Carlo as an efficient alternative to Dynamic Monte Carlo. This work implements the frequency transform method with thermal feedback in high-order/low-order Monte Carlo by blending static coupling methods such as single-batch Monte Carlo, with a simple thermal-fluids calculation. While previous work focused solely on prescribed transients, the addition of time-dependent thermal-fluids allows transients to be self-propagating. Tests were run in a fluids-initiated transient to showcase the basic functionalities of this methodology. Preliminary results behave as expected, paving the way for studying more sophisticated thermal-fluids coupling methods. (authors)
The dispositioning of nuclear waste generated at facilities across the country is an ongoing battle that affects us all. National laboratories and research centers dealing in medical research, clean energy, and other nuclear activities such as the Department of Energy (DOE) facilities face the need to properly manage and dispose of nuclear waste. A dynamic modeling solution would enable the DOE and others to make decisions on waste disposal and technological options. In doing so, this research explores modeling techniques using available data to address these situations. The focus being on developing an initial robust and adaptable discrete event model using the ExtendSim tool. This modeling effort will target the dispositioning of transuranic waste at the Savannah River National Laboratory (SRNL) which can be expanded to represent the current state of disposition process for waste generated at other DOE facilities. The model aims to assess resource allocation and waste processing options to stabilize productivity and cut the backlog of nuclear waste. By assessing the results of different scenarios, this research aims to provide actionable insights for the DOE. This approach has the potential to significantly improve the management of radioactive waste, offering the capability of evaluating options for optimizing the process for nuclear waste disposal. The findings of this study can serve as a valuable resource for decision-makers and other national laboratories, or research entities engaged in nuclear operations by enabling them to make more informed choices.
Understanding how irradiation affects the thermo-physical and mechanical properties of nuclear materials, such as thermal conductivity degradation in fuels and embrittlement of structural components, is critical to the safety and efficiency of nuclear reactors. These effects are largely governed by the formation and evolution of atomic-scale point defects and defect clusters. Due to their small sizes, however, these defects are invisible under high-resolution scanning transmission electron microscopy. This project aims to fill this experimental knowledge gap by integrating density functional theory (DFT), machine learning interatomic potential (MLIP), and kinetic Monte Carlo (KMC) techniques to predict longtime evolution of irradiation-induced defects in nuclear fuels.
Application packaging and distribution are the final steps for delivering software to end-users; both are frequently neglected when creating scientific software. Commercial businesses rely on electronic distribution systems that have rendered disk drives obsolete. Still, national laboratories continue to rely heavily on removable media to distribute and limit access to controlled applications. With increasing concerns of unauthorized copying of sensitive applications, a modern distribution system that utilizes cryptographically secure communication and authentication protocols has been developed. This new distribution system will secure the chain of custody for nuclear software while simultaneously simplifying access to these tools. This report summarizes four primary advancements made toward the secure distribution of Nuclear Energy Advanced Modeling and Simulation (NEAMS)-developed, Multiphysics Object Oriented Simulation Environment (MOOSE)-based applications: application installation, package distribution, automated package building, and distribution of documentation. NEAMS is currently developing more than ten separate applications based on the open-source MOOSE Framework. Distribution of these applications has primarily been accomplished by distributing source code, with end-users compiling the applications themselves. This work created a mechanism where MOOSE applications can be installed in a similar way to any other software. This allows both administrators and end-users simplified access to runnable executables. With this new installation capability, it was then possible to rethink distribution. A new, secure capability for delivering MOOSE-based applications over the internet has been created. This system requires unique cryptographic tokens for authentication, greatly securing the custody chain for software. Once granted access, installation of any NEAMS code can be accomplished with these terminal commands: "conda install ncrc" "ncrc install ncrc-bison." After these two commands (and authenticating) the BISON application will be securely down- loaded from Idaho National Laboratory (INL)’s servers, installed, and ready to use. To enable this new distribution capability to be successful, the open-source Continuous Integration, Verification, Enhancement, and Testing (CIVET) Continuous Integration (CI) capability was augmented to add Continuous Delivery (CD). CD enables the automated building and packaging of MOOSE-based applications as they are modified by development teams, ensuring that our customers can obtain up-to-date versions of the software at any time. The need for instruction on how to use these applications was addressed through modifications to the MOOSE documentation system. The MooseDocs capability, which enables robust documentation of MOOSE-based applications, has been extended to allow both for the installation of documentation and the packaging of documentation with installed applications. Together, these enhancements form the core of a new, secure distribution mechanism for nuclear simulation tools. In concert with the Nuclear Computational Resource Center (NCRC), NEAMS- developed applications will now be straightforward to obtain securely.
The direct Urca process of rapid neutrino emission can occur in nonuniform nuclear pasta phases that are expected in the inner crusts of neutron stars. Here, the periodic potential for a nucleon in nuclear pasta allows momentum conservation to be satisfied for direct Urca reactions. We improve on earlier work by modeling a rich variety of pasta phases (gnocchi, waffle, lasagna, and anti-spaghetti) with large-scale molecular dynamics simulations. We find that the neutrino luminosity due to direct Urca reactions in nuclear pasta can be 3 to 4 orders of magnitude larger than that from the modified Urca process in the NS core. Thus neutrino radiation from pasta could dominate radiation from the core and this could significantly impact the cooling of neutron stars.
The overall aim of the CyclusJS project has been to design novel approaches that support fuel cycle simulations and empower decision-makers and researchers with better tools to compare and contrast results from multiple simulations to understand the landscape of potential outcomes better. In particular, nuclear fuel cycle simulations focus on modeling the nuclear industry and ecosystem at a macroscopic level, and the analysis tools connect such simulations with decision support systems. As a demonstration of the novelty and effectiveness of the proposed approach, we study scenarios for transitioning from one technology, Light Water Reactors (LWR), to a newer Sodium-cooled Fast breeder Reactor (SFR) technology.
Hanford Site nuclear waste is to be vitrified at the Waste Treatment and Immobilization Plant (WTP), which is a part of the safe and efficient retrieval, treatment, and disposal mission of the U.S. Department of Energy Office of River Protection. Hanford tank 241-AP-105 (referred to herein as AP-105) is one of the initial Hanford radioactive tank wastes planned to be processed and vitrified. A portion of AP-105 waste was retrieved by Washington River Protection Solutions, LLC (WRPS) and transferred to Pacific Northwest National Laboratory (PNNL). The waste went through dilution by Columbia River water to reach a target sodium (Na) concentration of 7 M, solids filtration, and cesium removal by ion exchange. A glass composition was calculated from the Kim et al. glass models to satisfy the WTP baseline requirements based on the as-received sample and the target dilution to 7 M, from which a simulant was calculated and glass forming chemical (GFC) additions were determined to form a liquid/solids mixture called melter feed. To prepare for the processing of the 7 M Na AP-105 waste melter feed and learn about the production expectations, the melter feed simulant of 7 M Na AP-105 waste was processed in a non-radioactive, continuous laboratory-scale melter (CLSM) system.
Nuclear tank waste at the Hanford site is slated for vitrification in large scale refractory-lined melters to transform it into a stable borosilicate waste form suitable for long-term storage or disposal. Molten glass corrodes the refractory lining over time at a rate correlated to the velocity of molten glass against the refractory-lined wall. To properly design a melter and plan for maintenance, it is necessary to accurately model the glass flow and find the correlation between flow velocity and corrosion rate. This modeling was started in STAR-CCM+ CFD software and is being continued in Nek5000 for its fast-running and quick turnaround code. This poster explains the basic process of reconstructing the geometry, fluid properties, and heating of test melters in Nek5000. The geometry is imported from a mesh file and boundary conditions assigned based on surface IDs. The fluid properties are set in the .par and .usr case files. The heating of the actual crucibles in joule heating with electrodes but simulated in Nek5000 with a volumetric heat source term. This volumetric heat source is fit to the actual heating profile with piecewise polynomials and exponentials. Basic results and verification methods are explained in the poster, as well as further work that must be done to complete modeling.
The nuclear industry intends to increase burnup past the regulatory limit of 62 GWd/MTU, and to be successful, the dispersal behavior of fragmented high-burnup fuel during loss-of-coolant accidents must be better understood. Part of this behavior is fuel fragmentation, relocation, and dispersal. This work reports on testing concerning the dispersal of fuel from rods following burst and pressure blowdown through the burst. The Fuel Relocation Induced by Accident Recovery (FRIAR) system was used to examine surrogate material and pre-burst cladding specimens to evaluate the effects of vibrational loading, mixture composition, and burst size on post-burst dispersal. Vibrational loading resulted in maximum dispersal around 15 Hz when the loading was performed with peak-to-peak oscillation amplitudes at 0.5 and 1 mm. Here, under this aggressive vibrational loading, dispersal was found to become significant as burst widths increased to approximately 5 mm and larger. However, below 4 mm widths, little material dispersed out of the opening, even at 1 mm amplitudes and for more dispersive mixtures (smaller average fragments). At and above 7 mm widths, dispersal was rapid for all mixtures. Finally, dispersal in the threshold region (4–5 mm wide bursts) strongly depended on the mixture. For mixtures with large fragments, dispersal was minimized to a few fragments, as the large fragments proved to be effective at minimizing dispersal. Conversely, mixtures with smaller fragments showed increased dispersal in this region.
Nuclear energy is a promising low-carbon energy candidate to meet the increased demand for green energy, where the integration of fuel recycling can have significant benefits for material usage and waste reduction. Utilizing in situ monitoring tools can provide ample opportunities to better control and safeguard nuclear material recycle processes while also offering knowledge and insight into real-time solution properties. The simultaneous measurement of analytical targets in multiple process locations can enable real-time mass balance and material accountancy calculations. This is demonstrated here with a mass balance study of Nd 3+ on countercurrent aqueous/organic metal extraction within a single centrifugal contactor. The Nd 3+ concentration was simultaneously monitored at the inlets and outlets of both aqueous and organic phases using a visible absorbance detector that allowed for the simultaneous measurement of up to six locations. The Nd 3+ concentration was calculated by using chemical data science algorithms, where model training sets were collected on a single track of the detector. The discussion includes addressing the challenges of using a model collected on a single track and applying it as a model across the other tracks on the detector. Each track of the detector corresponds to one measurement location on the contactor. The difference in the integrated moles of Nd 3+ between the inlet and outlet at the end of the experiment was near zero, indicating that the mass balance of this experiment was maintained. Overall, the online spectroscopic monitoring was able to follow changing solution conditions and accurately measure the concentration of Nd 3+ in different locations within the contactor system.
Hanford Site nuclear waste is to be vitrified at the Waste Treatment and Immobilization Plant (WTP), which is a part of the safe and efficient retrieval, treatment, and disposal mission of the U.S. Department of Energy Office of River Protection. Hanford tank 241-AN-107 (referred to herein as AN-107) is one of the initial Hanford radioactive tank wastes planned to be processed and vitrified. A portion of AN-107 waste was retrieved by Washington River Protection Solutions, LLC (WRPS) and transferred to Pacific Northwest National Laboratory (PNNL). Compared to previously received and vitrified wastes (AP-107, AP-105, and AP-105), the concentration of organics in AN-107 was greater by an order of magnitude, while the activity of radionuclides was multiple orders of magnitude greater.
Hanford Site nuclear waste is to be vitrified at the Waste Treatment and Immobilization Plant (WTP), which is a part of the safe and efficient retrieval, treatment, and disposal mission of the U.S. Department of Energy - Hanford Field Office. A portion of Hanford tank 241-AW-105 (referred to herein as AW-105) waste was retrieved by Hanford Tank Waste Operations and Closure (H2C) and transferred to Pacific Northwest National Laboratory (PNNL). Compared to previously received and vitrified wastes (AP-107, AP-101, AN-107 and AP-105), the concentration of potassium in AW-105 was greater by an order of magnitude.
As fuel vendors and designers pursue the development of advanced reactors or the increase in burnup limits for existing reactors, advanced computational tools are necessary to understand the fuel performance. BISON, a fuel performance code developed primarily at Idaho National Laboratory, which is built upon the Multiphysics Object-Oriented Simulation Environment (MOOSE), provides capabilities to analyze multiple nuclear fuel forms in a wide variety of dimensions. Since its inception, BISON has been used to investigate the performance of light-water reactor fuel rods, accident tolerant fuel concepts, metallic and mixed-oxide (MOX) fuels for fast reactors, plate fuels for research reactors, and tri-structural isotropic (TRISO) fuel particles. This talk will provide a history of BISON, highlights of major development milestones, tributes to key contributors, and applications of its use to various fuel forms in one-, two-, and three-dimensions.
Abstract The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, f V , was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p . The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of f V vs. p did not change significantly. The percolation threshold p c was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.