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 summary report provides an overview of all LDRD projects at Argonne that concluded in Fiscal Year 2021. Many projects are funded for multiple years, the initial fiscal year for each project is indicated by the first four digits of the LDRD project number.
The Advanced Dimensional Depletion for Engineering of Reactors (ADDER) software is being developed in the Research and Test Reactor (RTR) Program at Argonne National Laboratory to meet the reactor design and analysis needs of the Conversion Program. ADDER is a flexible tool that (1) provides a depletion capability through coupling external neutronics codes with a built-in CRAM solver or external depletion code and (2) provides a user-friendly interface to perform fuel management and criticality search operations. The ADDER software is a Python 3 application written using modern software development practices subject to a compliant implementation of NQA-1 and applicable Department of Energy software quality assurance standards. This report is the user guide for the software release referred to as ADDER v1.0.1. The motivation for a software to have flexible capabilities that ADDER possesses is the need to support a wide variety of geometries that are commonly required in analysis of research and test reactors. These reactors can have complex fuel, experiment, or control material shuffling patterns that persist over several years with many fuel management and partial refueling intervals. The scale of fuel management analysis can require tracking of an inventory that is multiple times the core loading. Many reactors, both power and non-power reactors of various types, will find the features of ADDER useful to facilitate key tasks that a fuel or core design engineer must perform with the convenience of concise input and validated functionality.
This report describes improvements to the Feasible Actuator Range Modifier (FARM) component of the RAVEN-based HYBRID framework for analysis of Integrated Energy Systems (IES). FARM supports the HERON plug-in that solves the power dispatch problem. The solution to the dispatch problem involves economically optimal dispatches that satisfy limits on production variables and their rates of variation (explicit constraints) as well as process variables tied to the service life of equipment (implicit constraints). FARM serves to validate or confirm that a HERON solution for explicit constraints also satisfies the implicit constraints. FARM-alpha was released by Argonne National Laboratory in January 2021 followed by FARM-Beta in January 2022 with the latter providing increased flexibility for the user. In this report, FARM-Gamma, the latest version of the code, is described. The major improvement is the implementation of a system identification algorithm based on the Dynamic Mode Decomposition with Control (DMDc) coupled with a “Rolling Window” scheme that allows obtaining linear time-varying state-space models. This feature equips FARM with the most accurate approximation of system dynamics, and it relieves the user from the burden of performing an exhaustive off-line characterization of the dynamics. FARM-Gamma capabilities are assessed by solving the power dispatch problem for a representative IES unit. The simulation times corresponding to the different releases are estimated and compared. These values capture the increasing computational burden of the successively higher fidelity state-space models adopted by FARM-Alpha, FARM-Beta and FARM-Gamma. The code implementation provides significant flexibility, i.e., the user can always select the most suitable version of FARM according to the problem to be solved and the available computational resources. It is anticipated that FARM will play a role in addressing several future IES applications. We outline how it can support the coordinated management and safe operation of a nuclear plant coupled to industrial processes to produce hydrogen and synfuels.
The System Analysis Module (SAM) is under development at Argonne National Laboratory as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. It utilizes the object-oriented application framework MOOSE to leverage the modern software environment and advanced numerical methods. The capabilities of SAM are being extended to enable the transient modeling, analysis, and design of various advanced nuclear reactor systems. The molten-salt-cooled pebble-bed reactor, or pebble-bed FHR (PB-FHR) is a promising candidate among advanced nuclear reactor concepts with its improved passive safety characteristics and high thermal efficiency. To support the development and utilization of the SAM code for PB-FHR safety analysis, activities on SAM code enhancements, reference plant model developments, and code validations have been performed in the past a few years to support near-term industry and NRC needs. This report summarizes recent progress under DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program in SAM code development and demonstration for transient safety analysis of Fluoride-salt-cooled High-temperature Reactors. SAM capabilities has been significantly enhanced over the years to add FHR specific modeling features, including salt freezing and thawing, spherical core channel and pebble bed core modeling, solid-fluid thermal radiation, tritium transport and general species transport in fluids and solids, and the general code enhancements on solver schemes of point kinetics module and reactivity feedback models. A reference PB-FHR model is developed, based on publicly available information from Kairos Power’s generic FHR design and the University of California, Berkeley (UCB) Mk1 design. A reference reactor model is foundational to the methodologies employed by NRC to verify the adequacy of computer codes and evaluation models. The reference FHR model was utilized for a number of selected FHR design basis accidents, including station blackout, loss of heat sink, loss of flow, transient overpower, and overcooling events.
The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code.
The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code.
PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.
This report documents the Argonne-Exelon effort to develop and utilize an agent-based model (ATEAM) of charging demand and infrastructure expansion applicable to the Washington, DC–Baltimore, MD consolidated metropolitan area. This study extends the ATEAM model time horizon to 10 years (from 2020 to 2030), expands agent behavior modeling capabilities, incorporates more granular and extensive empirical data on charging behavior, and analyzes charging needs for a much larger population of PEVs, in keeping with regional goals for significant adoption of ZEVs. With given targets for annual BEV adoption, five scenarios were developed to examine public infrastructure needs and resulting charging load, considering different home charging availabilities, as well as different PEV consumer profiles and public charging infrastructure deployment strategies. Scenario results show that if new chargers (both L2 and DCFC) are spread more widely (as with ubiquitous deployment strategies), there will be less variation in the number of chargers added to each census tract in the study area. More importantly, widespread public charging infrastructure with ubiquitous deployment strategies reduces unmet charging demand and improves charging success, even with heavy reliance on public charging. About 80 percent of BEV drivers can charge on their first attempt in scenarios with ubiquitous deployment strategies. Moreover, widespread public charging infrastructure better meets the demand for more charging, and in return, increases BEV adoption. Low home charging availability produces higher charging loads in public locations, especially during the early morning (around 8:00 a.m.) and late afternoon (around 6:00 p.m.). The evening peak load indicates that drivers are taking advantage of public charging before heading home. Study results also indicate that even with 20 percent home charging availability in 2030, just 20 percent of drivers attempt to charge on a given day. With their relatively high electric range (200+ miles), the BEVs expected to be on the road in 2030 can handle daily commutes without re-charging for a couple of days.
An NNSA-Argonne National Laboratory workshop to define NNSA and external partner mission needs and opportunities at the Advanced Photon Source was held February 15-17, 2022. This workshop summary provides a fresh look at NNSA’s evolving mission needs, partner agency and university interests, x-ray photon characteristics, and technical requirements to meet those needs. The workshop was organized according to four breakout groups: Coupling dynamic drivers to a synchrotron, High explosives and chemical sciences, Process-structure-performance: X-ray experiments to understand materials behaviors, and Process-structure-performance: manufacturing and aging. This summary report describes the workshop outcomes, and the preliminary conceptual design for a dedicated Defense Materials Science Sector (DMSS) at the Advanced Photon Source.
This report summarizes spreading and heat transfer results calculated using MELTSPREAD for a scenario of molten FLiNaK spilling onto a flat stainless steel substrate that was inspired by the “maximum credible accident” scenario of the Molten Salt Reactor Experiment (MSRE). MELTSPREAD was developed at Argonne to model the one-dimensional flowing and freezing of molten corium and was applied to model molten salts for the first time during FY21. The molten salt MELTSPREAD model that was developed during FY21 was updated as part of this work using a corrected heat of fusion value of FLiNaK and to accommodate a larger salt spill volume. The model was run with and without the inclusion of contributions from decay heat and a sensitivity analysis of initial spill conditions was performed to determine the importance of those factors on model outcome. The results provide insight into the expected spreading and heat transfer behavior of simulated and irradiated fuel salt that has been spilled and indicate aspects of the model to be revised in future versions of the MELTSPREAD model for applications to molten salt.
The DIF3D software is a key part of the legacy fast reactor analysis software package developed and maintained at Argonne National Laboratory. The DIF3D-VARIANT capability has been pushed to the forefront of usage in the last decade at ANL because of improved computing resources and the focus on reactors with significant leakage. The purpose of this document is to discuss when the DIF3D-VARIANT solver provides a solid benefit to analysis work and discuss the resource requirements and explain why they are so much greater than the other solver options in DIF3D. Additional discussion will be provided on solution accuracy and any deficiencies relative to other solvers available in the literature.
The main objective of this CRADA is to provide Argonne National Laboratory (ANL)’s modeling and simulation capabilities via ARPA-E's MEITNER (Modeling-Enhanced Innovations Trailblazing Nuclear Energy Reinvigoration) Resource Team (RT) arrangement to support the demonstration of the viability of HolosGen’s Holos-Quad reactor design. The Holos-Quad reactor design is an advanced reactor concept that incorporates many new design features, such as the neutron-coupled Subcritical Power Modules (SPMs) in its core design, elimination of balance of plant (BOP) by direct integration of a helium Brayton cycle power conversion system with each SPM, among many other innovative features. The demonstration of the viability of such an innovative reactor design warrants iterations of modeling and simulation and testing. The purpose of this project is to utilize ANL’s modeling and simulation capabilities in nuclear reactor analysis and power conversion system analysis to inform HolosGen and the Design Team (DT) in the design and optimization of the Holos-Quad concept. The major work scopes of this project include: to investigate conceptual designs and materials for radiation shielding to protect personnel and internal components such as turbomachinery; Assess the the helium Brayton cycle power conversion system performance in both nominal and load following conditions; Investigate power conversion components and overall system performance; Identify control strategies to enable load following; Perform simulations of HolosGen’s subscale simulator and using available test data for code validation/benchmark purpose; Perform core thermal-hydraulics, safety analysis, and structural analysis.
NREL, Argonne National Laboratory (ANL) and Oak Ridge National Laboratory (ORNL), which are also referred to collectively herein as the "Contractors" will perform electrode ink development and characterization, focusing on 1) inks with appropriate properties for R2R (gravure and slot-die) coating onto PFSA membrane and 2) catalysts suitable for electrolysis. Participant will develop acceptance criteria and characterization methods for direct coating on the membrane. The purpose of this project is to show feasibility for cost reduction of the electrode through increased throughput and improved material utilization. Successful completion of this project will result in more affordable hydrogen production via proton exchange membrane (PEM) electrolysis by substantially reducing the anode precious metal content and catalyst-coated membrane processing costs. Assists laboratory in achieving programmatic scope, adds new capability to the laboratory's core competencies. This project uses and enhances the laboratory's core competencies, and enhances U.S. competitiveness by utilizing DOE developed intellectual property and/or capabilities. This project will benefit the Fuel Cell subprogram within the EERE Fuel Cell Technologies Office (FCTO) by actively developing materials, formulations, components, and processes for PEM electrolyzer stacks for overcoming critical technical barriers to widespread hydrogen production for zero-emission transportation applications. Through this JWS, NREL, ANL and ORNL are assisting a key domestic industry partner, Proton Energy Systems, Inc., in the area of low-cost electrolyzer electrode formulation and coating technology. After completion of this project, NREL, ANL and ORNL will be able to verify whether the direct deposit coating processes investigated are suitable for producing CCMs for PEM electrolyzers, which is of great interest and benefit to the DOE EERE Fuel Cell Technologies Office (FCTO) and H2@Scale missions.
The optimization of the fuel-engine interface and associated performance is a significant goal for CoOptima and the industry. The strategic driver for fuel-engine optimization is clear: diesel optimization can significantly reduce emissions including NOx, particulate matter, and CO2 emissions. Previous evaluations of ethanol and diesel blends raised significant concerns with fuel injector erosion using ethanol blends. While Sylvatex’s technology can use ethanol, it also utilizes a surfactant and water to compatibilize and modify ethanol blends, thereby controlling its impact in diesel engines. The results from the project will (1) provide important insights into complex fuels and erosion, which will be used to design additional, improved formulations; and (2) expand capabilities of Argonne’s cavitation-induced erosion risk assessment (CIERA) tool to model multi-component systems, furthering Co-Optima objectives.
In high-performance Monte Carlo (MC) radiation transport codes, cross section lookups often account for the majority of computational expens. In response to this observation, a number of nuclear data mini-apps have been developed to profile and optimize the lookup process. A review of these mini-apps can be found in Ref. In contrast to a fully fledged MC code, nuclear data mini-apps do not feature particle tracking; instead, the cross section lookup process is considered in isolation, a simplification that allows for easier implementation and evaluation of various lookup schemes and techniques. GIDI+, as described on its public repository1, is ...a collection of C++ libraries for accessing evaluated and processed nuclear data stored in the Generalized Nuclear Database Structure (GNDS). In addition to reading GNDS files, GIDI+ has functions to sum and collapse multi-group data as needed by deterministic transport codes, and to sample GNDS data as needed by Monte Carlo transport codes. Certain modes of Mercury, a LLNL-developed MC radiation transport code, utilize GIDI+ for cross section lookups. Specifically, Mercury uses the MCGIDI library contained within GIDI+ to perform lookups on GPUs. A mini-app exercising MCGIDI’s cross section lookup capabilities would prove useful for optimization, which in turn could improve the performance of client codes such as Mercury. Such a mini-app could also be used to assess the performance of MCGIDI compared to other nuclear data mini-apps. In this document, we introduce a nuclear data mini-app built with the MCGIDI library. The MCGIDI mini-app (MCGIDI-MA) was developed to closely resemble the methodology present in XSBench, a nuclear data mini-app developed at Argonne National Laboratory. The remainder of this section provides an overview on XSBench. Section 2 describes MCGIDI-MA and its usage. Section 3 contains performance results of MCGIDI-MA on LLNL’s Lassen and Quartz compute platforms, as well as comparisons to XSBench performance where appropriate. In Section 4, we present a summary of our findings.
Municipal solid waste infrastructure in the United States contains complex supply chains to transport, sort, and sequester waste. This project developed models to estimate the environmental and economic costs of waste infrastructure in the U.S. using data from several municipalities. Life cycle costing (LCC) models were developed for waste management relevant supply chain steps, including waste collection, sorting, landfilling, waste-to-energy, ash processing, and recycling. A waste-to-materials and energy (WTM&E) life cycle analysis (LCA) and LCC dashboard was developed in collaboration with Argonne National Laboratory (ANL) and an LCC tool was developed for additional analysis.
This report presents the load following analysis of the “Holos-Quad” 10MWe micro-reactor design developed by HolosGen LLC. This analysis is focused on the Holos-Quad micro-reactor ability to match the changing grid demand at 10%/min rate. The control mechanisms for the plant are identified, simulated, and compared. Based on the comparison results, control strategy for load following of the Holos-Quad micro-reactor is developed. The control strategy and load following capabilities are demonstrated in a full-range down-and-up reactor power transient from 100% to 0% and back to 100% load at 10%/min rate. All calculations are carried out with Argonne’s Plant Dynamics Code.
The Subcritical Power Module Sub-scale Simulator (SPM-SS) has been designed and constructed by HolosGen LLC under the ARPA-E MEITNER program to simulate the thermal-hydraulic and heat transfer behavior of the full-scale Holos-Quad Subcritical Power Modules (SPMs). Four coupled SPMs, each rated at 5.5MW, form the Holos-Quad gas-cooled microreactor design. The SPM-SS represents a substantially scaled-down system with a power rating less than 40 kW, equipped with an electrically heated fuel cartridge heat exchanger, an electrically heated compressor heat exchanger, and a valve actuated turbine heat exchanger, in addition to a recuperator and a cooler heat exchanger. The fuel cartridge represents a portion of the full-scale SPM core, the compressor heat exchanger mimics the temperature changes resulting from the compressor’s turbomachinery inefficiencies, the turbine heat exchanger mimics the expansion process normally occurring through the turbine, while the recuperator and cooler heat exchangers complete the subscale simulator loop. The heaters equipping the fuel cartridge and the compressor heat exchangers are electronically controlled to simulate normal and off-normal SPM operating conditions. The full-scale Holos-Quad SPM design eliminates the traditional balance of plant and executes thermal-to-electric energy conversion by means of an intercooled Brayton cycle with decoupled compressor-turbine turbomachinery. The Holos-Quad full-scale design is equipped with a multi-stage axial Low- and High-Pressure compressor, and a multistage axial turbine. The SPM-SS is designed for testing and validation of selected components which are instead coupled by a traditional balance of plant. The SPMSS is not equipped with turbo-machinery (compressor and turbine) as the development of these components were excluded from the scope of work under the ARPA-E MEITNER funding program. The SPM-SS balance of plant enables modifications, replacement and testing of individual components with different working fluids and is designed to include the turbo-machinery components that will be developed in future research . The SPM-SS can be operated with different gases, variable mass-flow-rates, pressures, and temperatures to obtain test data for selected components, whose performance can be scaled to validate the computer model of the full-scale SPM at various conditions (e.g., start-up, transients conditions). The SPM-SS can operate at the maximum Holos-Quad design pressure of 7 MPa, and a maximum temperature limited to 650 °C by the electrical heaters. Several SPM-SS tests have been conducted and analyzed with the Plant Dynamics Code (PDC) developed at the Argonne National Laboratory (ANL). These tests aimed at validating the PDC modeled predictions of the full-scale Holos-Quad design with data from selected SPM-SS components. In order to address SPM-SS specific characteristics, such as components heat losses and absence of turbomachinery components, some modifications to the PDC have been implemented to factor the design differences from the full-scale Holos-Quad SPM to the SPM-SS. As the PDC offers capabilities to analyze systems with different working fluids, air, nitrogen, and helium were utilized as the SPM-SS working fluids. Air was utilized to fine-tune the SPM-SS Systems Structures and Components (SSCs), nitrogen was utilized to pressure test the SPM-SS loop at the SPMs design maximum pressure of 7MPa. Helium was utilized as the working fluid circulating through the SPM-SS SSCs for specific tests to validate the PDC predictions of the fuel cartridge heat exchanger. SPM-SS tests data were also analyzed with both the steady-state and transient analysis capabilities offered by the PDC. This report describes the PDC analysis of the SPM-SS tests data, including the necessary code modifications and comparison of the code results with the experimental data. Based on the results, a discussion is presented on how the analysis supports design and transient calculations of the full-scale Holos-Quad microreactor. Also based on the results of this work, recommendations are made for future optimizations of the SPM-SS components and PDC model development needs.