Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “USeS”

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 739 records · Page 41

Electrification Futures Study: Methodological Approaches for Assessing Long-Term Power System Impacts of End-Use Electrification

By its nature, electrification enhances the coupling between the electric sector and end-use sectors. Assessing the impacts of electrification requires both an examination of the complex interactions between sectors and a broader assessment of multiple parts of the energy system. The Electrification Futures Study (EFS) uses several complementary modeling and analysis tools to analyze the impacts of electrification on the U.S. energy system. In particular, the EFS relies on an overarching scenario analysis approach, but through the use of separate modeling approaches designed to assess various electricity demand- and supply-side futures. The primary model employed to generate the supply-side scenarios is the Regional Energy Deployment System (ReEDS) model, which is a capacity expansion model for the U.S. electricity system through 2050. Traditionally, the model has been primarily exercised in scenario analysis that implicitly assumed limited electrification. Because of this assumption, resulting future load profiles are approximated by historical ones, load growth is driven primarily by population and economic growth only, and changes in direct end-use natural gas consumption do not effect natural gas costs for electricity generation. In this report, we (1) reflect the potential for resource sharing between regions given these changes in demand, (2) represent how changes in natural gas consumption in end-use sectors could impact the economics of natural gas-fired generation, and (3) document a new model representation of demand-side flexibility used for the EFS. These improvements to ReEDS are employed for the EFS supply-wide analysis, which is summarized in a companion EFS report titled Electrification Futures Study: Scenarios of Power System Evolution and Infrastructure Development for the United States (Murphy et al. 2019). The data and methods documented in this report could also be adapted for other models with similar scope and limitations as ReEDS, and these data and methods could be used to assess future electric system scenarios.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

First-principles Modeling and Design of Solid-State Interfaces for the Protection and Use of Lithium Metal Anodes

Li-ion batteries are one of the most advanced energy storage technologies in use today. Li-ion batteries are used in a multitude of applications ranging from consumer electronics, medical devices, sensors and grid storage. However, improving the capacity and energy density delivered by current Li-ion technology requires advanced materials research into novel chemical systems. In this project we have focused particularly in the use of solid-state electrolytes with lithium metal electrodes. Research into all solid-state batteries (ASSB) with Li metal electrodes has significantly expanded in recent years, however most studies reported experimental findings, which left substantial room for theoretical and modeling work as a tool to understand and determine design principles allowing reliable and safe use of ASSBs with Li metal. Among the remaining obstacles preventing reliable use of ASSBs with a Li metal electrode, the stability of the interface between the solid electrolyte and Li metal, and the propagation/dendrite formation of Li metal and resulting mechanical degradation of the electrolyte are key phenomenon that are yet to be fully understood. In the current project we have addressed these two coupled phenomena using first principles calculations and mesoscale continuum modeling. We have obtained chemical and electrochemical stability windows for several solid electrolyte materials. Additionally, from mathematical and numerical modeling of Li protrusion and dendrite initiation during plating and stripping we have determined design criteria in terms of chemical, electrochemical, and mechanical properties and operating conditions for which stable deposition can occur. We also considered the effects of mixed electronic-ionic conduction in solid electrolytes, which has more recently been suggested as another important mechanism involved in ASSB failure. Throughout our work we have successfully addressed important questions necessary for the use of ASSB’s. We have determined guiding principles for materials properties and operating conditions necessary to operate ASSB’s. And have proposed novel solid electrolyte materials with predicted chemical stability an ionic conductivity. Although this represents significant progress in our understanding, open questions remain in order to fully develop reliable and safely operate ASSBs with Li metal. Future work, building on this project will require further experimental, theoretical and simulation efforts to address remaining questions.

25 ENERGY STORAGE↗

Blind Modeling Validation Exercises Using the Horizontal Dry Cask Simulator

The U.S. Department of Energy (DOE) established a need to understand the thermal-hydraulic properties of dry storage systems for commercial spent nuclear fuel (SNF) in response to a shift towards the storage of high-burnup (HBU) fuel (> 45 gigawatt days per metric ton of uranium, or GWd/MTU). This shift raises concerns regarding cladding integrity, which faces increased risk at the higher temperatures within spent fuel assemblies present within HBU fuel compared to low-burnup fuel (≤ 45 GWd/MTU). A dry cask simulator (DCS) was built at Sandia National Laboratories (SNL) in Albuquerque, New Mexico to produce validation-quality data that can be used to test the accuracy of the modeling used to predict cladding temperatures. These temperatures are critical to evaluating cladding integrity throughout the storage cycle of commercial spent nuclear fuel. A model validation exercise was previously carried out for the DCS in a vertical configuration. Lessons learned during the previous validation exercise have been applied to a new, blind study using a horizontal dry cask simulator (HDCS). Three modeling institutions – the Nuclear Regulatory Commission (NRC), Pacific Northwest National Laboratory (PNNL), and Empresa Nacional del Uranio, S.A., S.M.E. (ENUSA) – were granted access to the input parameters from the DCS Handbook, SAND2017-13058R, and results from a limited data set from the horizontal BWR dry cask simulator tests reported in the HDCS update report, SAND2019-11688R. With this information, each institution was tasked to calculate peak cladding temperatures and air mass flow rates for ten HDCS test cases. Axial as well as vertical and horizontal transverse temperature profiles were also calculated. These calculations were done using modeling codes (ANSYS/Fluent, STAR-CCM+, or COBRA-SFS), each with their own unique combination of modeling assumptions and boundary conditions. For this validation study, the ten test cases of the horizontal dry cask simulator were defined by three independent variables – fuel assembly decay heat (0.5 kW, 1 kW, 2.5 W, and 5 kW), internal backfill pressure (100 kPa and 800 kPa), and backfill gas (helium and air). The plots provided in Chapter 3 of this report show the axial, vertical, and horizontal temperature profiles obtained from the dry cask simulator experiments in the horizontal configuration and the corresponding models used to describe the thermal-hydraulic behavior of this system. The tables provided in Chapter 3 illustrate the closeness of fit of the model data to the experiment data through root mean square (RMS) calculations of the error in peak cladding temperatures (PCTs), PCT axial locations, axial temperature profiles, vertical and horizontal temperature profiles at two different axial locations, and air mass flow rates for the ten test cases, normalized by the experimental results. The model results are assigned arbitrary model numbers to retain anonymity. Due to the relatively flat axial temperature profiles, small temperature gradients resulted in large deviations of all models’ PCT axial location from the experimental PCT axial location. When the PCT axial location error is excluded in the calculation of the combined RMS of the normalized errors that considers PCT, the temperature profiles, and the air mass flow rates, the model data fits the experimental data to within 5%. When the vault information is excluded, the model data fits the experimental data to within 2.5%. An error analysis was developed further for one model, using the model and experimental uncertainties in each validation parameter to calculate validation uncertainties. The uncertainties for each parameter were used to define quantifiable validation criteria. For this analysis, the model was considered validated for a given comparison metric if the normalized error in that metric divided by the validation uncertainty was less than or equal to 1. When considering the combined RMS of the normalized errors of all metrics divided by their validation uncertainties, the model was found to have satisfied the criterion for model validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Design of Efficient Molecular Electrocatalysts for Water and Carbon Dioxide Reduction Using Predictive Models of Thermodynamic Properties

Objectives: 1) To design and synthesize water soluble transition metal complexes that function as electrocatalysts for the reduction of water to hydrogen or carbon dioxide to formate 2) To experimentally measure thermochemical properties, such as hydride donor ability, and use these quantities to optimize catalyst performance in solutions of varying pH 3) To incorporate results from detailed mechanistic and kinetic studies to direct improvements in catalyst design 4) To use thermochemical data from objective 2 to generate predictive models for the properties of new complexes to apply to the development of catalysts for other reductive reactions, with an emphasis on carbon dioxide reduction to methanol Description: Wide-spread implementation of renewable but intermittent energy sources, such as solar, requires the development of efficient methods for energy storage. The high energy density of chemical bonds makes chemical fuels an ideal solution for energy storage. Hydrogen and reduced carbon compounds have been proposed as ideal candidates for chemical fuels. However, the generation of chemical fuels from electricity requires competent electrocatalysts. The proposed research focuses on developing electrocatalysts for the reduction of water to hydrogen, and carbon dioxide to formate. Both products can be used directly as an energy carrier in fuel cells, or as a reductant for more saturated chemical fuels. Hydrogen is the most common fuel used in current commercial fuel cells. Additionally, there is interest in formic acid as a liquid carrier for hydrogen because of its increased storage density and ease of dehydrogenation. Formate is also an intermediate in the sequential reduction of CO 2 to methanol, another potential chemical fuel with high energy storage density. Despite the utility of formate as a chemical fuel or precursor, there are very few examples of electrocatalysts for its production from CO 2 , and even fewer demonstrate high product selectivity. In heterogeneous catalysis, the Sabatier principle is used to generate volcano curves that describe optimal thermochemical properties for key intermediates that result in peak catalytic activity. In most cases, such as hydrogen production and oxidation, the most favorable metals (Pt, Re, Rh, and Ir) are rare and expensive. Molecular inorganic complexes provide an opportunity to use electronic and steric ligand effects to tune the critical thermodynamic parameters of abundant metals to values comparable to key surface intermediates on precious metals. This principle will be applied to the design of aqueous homogeneous catalysts for the reduction of H 2 O and CO 2 optimized to function at specific pH ranges. The critical intermediate in both of these reactions is a metal hydride. The strength of this bond, or hydricity (ΔGH - ) dictates the overall thermodynamics of the reduction of H + to H 2 and the sequential reduction of C1 substrates, such as CO 2 , CO, and H 2 CO (the latter to CH 3 OH). ΔGH - will be systematically measured for a series of first row metal complexes to form predictive models for metal and ligand electronic effects. Benefits and Outcomes: This approach addresses thermodynamic requirements essential for energy-efficient catalyst design and transition state barrier considerations for faster catalysis. The advancements made from this research will include the discovery of efficient, stable, and fast catalysts for hydrogen and formate generation from electricity. The ability to store energy efficiently is critical to the widespread use of renewable energy schemes. This provides a viable solution to produce high energy density fuels for stationary and mobile energy applications. The fundamental research conducted is also directly applicable to the design of catalysts for the generation of further reduced products such as methanol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of High-Performance Ni-base Alloys for Gas Turbine Wheels Using a Coprecipitation Approach (Final Technical Report)

In this project, a synergistic approach involving experiments informed by modelling was used to develop γ'/γ'' coprecipitation-strengthened Ni-base superalloys for land-based gas turbine wheels. In thick sections of such turbine wheels, the slow cooling rate normally observed, as well as variations in cooling rate could potentially result in drastically inhomogeneous and coarsened microstructure. Thus, ten alloy compositions were selected based on coprecipitation and sluggish γ’ growth and coarsening during three iterations of alloy development, assisted by extensive CALPHAD modeling. Phase equilibrium calculations, informed by the characterization of early alloys, were successful in predicting suitable aging treatments to obtain bimodal coprecipitate distribution in several alloys. These aging treatments were applied to the subsequent round of alloys and the resulting microstructures were characterized in detail. Two different cooling rates were used during each of these heat treatments to simulate the range of cooling rates observed in a thick section forging. The microstructures were used to calibrate a multi-phase field model (MPF), which was used to study the formation mechanisms of various coprecipitate morphologies observed in the experiments. Following this, the most promising alloys in each round were selected and studied using high temperature tensile and creep tests. Considering that a slower cooling rate would lead to a strength debit, while a faster cooling rate might lead to a decline in creep rupture life of a full-scale wheel, tensile tests were carried out for slower cooled samples, while creep tests were carried out for faster cooled samples. The deformation microstructures were then characterized using advanced characterization techniques such as Diffraction Contrast-Scanning Transmission Electron Microscopy (DC-STEM) and Electron Channeling Contrast Imaging (ECCI). Further, a mean-field creep model was developed and calibrated using these experimental findings to predict long-term creep behavior of these coprecipitation-strengthened alloys.

36 MATERIALS SCIENCE↗

Perform Life Cycle Energy and GHG Emission Analysis, Select Candidate Refrigerant(s), Developing Heat Pump by Using Low-GWP Refrigerants and 5 mm Diameter Tubes

Reducing environmental impacts from the energy conversion systems has been an important research topic due to recent severe global climate changes. In residential buildings, the space heating and cooling are main energy consumers, which are relying on vapor compression based heat pumps. To mitigate the environmental impacts of residential systems, the heat pumps using low-GWP refrigerants and 5mm-tube heat exchangers are proposed to replace the conventional 9 mm or 7 mm tube heat exchangers. This research uses a simulation-based design method to optimize the refrigerant circuitries of a 5-ton low-GWP refrigerant heat pumps whose indoor and outdoor heat exchangers adopt 5 mm tubes. Heat pump systems using four commercially available low-GWP refrigerants (R454A, R454C, ARM20A and ARM20B) are optimized and compared with the baseline R410A system. The results indicate that directly replacing 9 mm tubes into 5mm tubes in the baseline systems induces heat pump performance degradation; by optimizing the 5 mm tube heat pumps with added tube rows, the system refrigerant charge is reduced up to 33% compared with the baseline; the optimized low-GWP 5 mm tube heat pumps offer the same or higher performance under both cooling and heating mode; that under cooling mode, the EER improvement is up to 10.7% and under heating mode the COP improvement is up to 3.4%; and a smart flow control device to maintain counter-flow circuitry pattern on both operating modes can improve the heating performance of 5 mm tube optimized system by up to 3.0%. In an effort to better understand the environmental impacts of these systems, life cycle climate performance (LCCP) evaluation method is used to assess the baseline and optimized systems. LCCP considers the direct and indirect emissions of the system over the course of its lifetime from manufacturing to disposal. Our analyses indicate that the optimized systems using 5mm tubes and low-GWP refrigerants can reduce life-cycle CO 2 emission significantly. And improving system efficiency is the most effective way to reduce the emission of low-GWP refrigerant systems. This report shows that using 5 mm diameter tubes in heat exchangers is a feasible approach to develop heat pumps with low-GWP refrigerants. Using the proposed design approach can yield higher system efficiency with lower environmental impact.

42 ENGINEERING↗

Development of Digital Twin Predictive Model for PWR Components: Updates on Multi Times Series Temperature Prediction Using Recurrent Neural Network, DMW Fatigue Tests, System Level Thermal-Mechanical-Stress Analysis

The long-term operation (LTO) of nuclear power plant (NPP) beyond their original design life of 40 years, can lead to more material damage associated with cyclic fatigue under thermal-mechanical loading cycles and associated long-term exposure of reactor material to the deleterious reactor-coolant environments. However, under this LTO condition the reactor components can still safely operate but may require more frequent Nondestructive Evaluation (NDE) of reactor components. Frequent NDE requirement may lead to frequent shutdown of the NPP. This in turn can lead to power outage and additional NDE-inspection-cost related economic loss. The economic loss can be minimized by reducing uncertainty in life estimation of safety-critical pressure boundary components and by implementing more digital approach such as by using upcoming digital-twin (DT) technology for predicting the structural states (e.g., time and location dependent inside/outside thickness temperature, stress, strain, plastic deformation, etc.) and associated fatigue life of a component in real time. Towards this goal Argonne National Laboratory (ANL) with the sponsorship of DOE Light Water Reactor Sustainability (LWRS) program is working on the development of a DT framework that can be used for real time environmental fatigue prediction of reactor components. The DT framework is based on limited experiment-data, Artificial-intelligence (AI) – Machine-Learning (ML) - Deep-Learning (DL) based techniques and Multiphysics-computational-mechanics such as finite element (FE) based modeling tools. Towards this overall goal, following are some of the major contributions made during the FY21: 1) Multiple 82/182 dissimilar metal weld (DMW) specimens (both solid-weld and joint-weld representing the actual reactor multi-metal nozzles) were fatigue tested. The resulting fatigue lives were compared to the NUREG-6909 based best-fit and design fatigue curves. Additionally, the results of 52/152 DMW fatigue specimens (which were recently tested at Republic of Korea under the sponsorship of International Nuclear Energy Research Initiative - INERI program) were compared to the NUREG-6909 based best-fit and design fatigue curves. From the comparison of 82/182 and 52/152 DMW test data with NUREG-6909 best-fit curve, most of the reported test data fall way away from the NUREG-6909 suggested best-fit or mean curve. The NUREG-6909 suggested best-fit curve is the best-fit curve of austenitic stainless steel and due to lack of enough data on Nickel-based welds, this is currently being used for predicting the life of Nickel-alloy-based welded components. However, the above observation may require higher scaling factor (e.g., ASME suggested factor of 20 on cycles rather than the current NUREG-6909 suggested factor of 12 on cycles) for scaling the austenitic-stainless-steel best-fit-curve for estimating the design or safe-life of a welded component. Accordingly, for example, if a DMW component experience a strain amplitude of 0.6% the PWR-water life of the component would be 52 cycles instead of 85 cycles. However, more DMW tests are required to further ascertain the above-mentioned observations. 2) A system level CAD and finite element model were developed which consists of reactor pressure vessel (RPV), part of steam generator (SG), part of pressurizer (PRZ), hot leg (HL), and surge line (SL). This is with detailed nozzle geometry and thermal-mechanical material properties of different metals to simulate realistic thermal-mechanical stress under connected system global thermal-mechanical boundary conditions. 3) Different system level heat transfer analyses were performed with estimation of relevant heat transfer coefficients. The resulting data were used in subsequent system level thermal-mechanical stress analysis and for generating spatial-temporal training and validation data for a system level digital-twin based temperature predictor. Transient heat transfer analyses were performed considering thermal boundary condition under design-basis (DB) loading and EDF (Électricité de France) data-based grid-load-following (EDF-GLF) loading cycles. 4) System level thermal-mechanical stress analysis was performed for identifying damage-prone hotspots and for future extension of the model for cyclic state prediction. From the system-level model simulation under DB loading cycle it is found that HL and the SL nozzle that connect to the HL can experience significant stress and strain and could be one of the weakest links in the overall reactor coolant system (RCS). 5) An AI/ML based DT model was developed for multi-time-series temperature prediction at any inside/outside thickness locations of PWR pressure boundary components. This is by using Recurrent-neural-network (RNN) and keras machine learning libraries. The RNN model was validated against two laboratory test-based data sets with one obtained through ANL’s in-air fatigue test system and other through PWR-water test loop. The experimentally validated DT model further validated against FE model results to predict thermal scarification related spatialtemporal temperatures at random locations of a component. The well validated DT model was then used for demonstrating spatial-temporal temperature prediction under 100+ years of reactor operation subjected to combined DB, EDF-GLF and randomized grid-load-following (RANDOMGLF) loading Cycles. The expert-elicitation DT model framework was developed assuming field/input/process measurements can be available from a few existing plant sensors and can readily be used by the NPP operators. The above temperature prediction model will feed to the next-step stress analysis model based on which the life of a component can be predicted in realtime, which is one of our future works.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Concentrating Rare Earth Elements in Coal mine drainage Using Coal Combustion Products through Abandoned Mine Land Reclamation

Rare earth elements (REEs) (including scandium, yttrium and a group of 15 lanthanides) are often considered to be critical components in the productions of renewable energy hardware, electric vehicles, health care and military equipment, and consumer electronic products. The demand of REEs has been projected to be growing at an annual rate of 5-9% in the next 25 years. In 2011, the global demand of total rare earth oxides (REOs) was estimated to be approximately 105,000 tons, which is expected to grow to 210,000 tons by 2025. China overwhelmingly dominates the current worldwide rare earth productions but has strategically restricted its exports, causing significant instability for the global market. In response to the increasing demand for REEs and the supply dominance of China, identifying alternative sources of REEs has become a critical issue for the United States and other countries. Coal, coal ash, and coal mine drainage (CMD) are considered to be the alternative sources of REEs. In the U.S., high REE concentrations have been reported to be closely associated with coal deposits, including the Appalachian Basins. When surface and/or groundwater come in contact with geologic strata containing sulfide minerals exposed by coal mining, the accelerated oxidation of sulfide minerals in the presence of ferric iron and/or oxygen can produce sulfuric acid. The process promotes the weathering of REE-bearing rocks and minerals in the host geologic strata. Compared to average river water and seawater, the concentrations of REEs can be orders of magnitude higher in CMD. In this study, we demonstrated a trap-extract-precipitate (TEP) process that can effectively recover REEs from CMD. The three-stage TEP process uses alkaline industrial by-products to capture REEs from CMD and then applies an extraction/precipitation procedure to produce a feedstock that can be economically processed to produce marketable rare earth oxides. The alkaline industrial by-products tested in this study include the residual from a water softening process (DRWP sludge) and two types of stabilized flue gas desulfurization materials (sFGDs). sFGD material is a mixture of lime (CaO) and two coal combustion by-products, calcium sulfite FGD by-product and fly ash. The objectives of this study are to (1) validate the effectiveness and feasibility; (2) determine mechanisms controlling the rare earth recovery, (3) quantify the associated economic and environmental benefits, and (4) evaluate the full-scale application. To achieve these objectives, tasks to be carried out in this proposed project are organized into three phases. In the first phase, the research team collaborated with Ohio Department of Natural Resources, American Electric Power, The Wilds (a nonprofit wildlife conservation organization), and a private landowner to carry out field investigations aimed to screen and evaluate the seasonal changes of rare earths in the CMD discharges that have high recovery potentials. Next, the recovery of REEs from CMD was tested using a series of lab-scale column and batch tests under, respectively, percolation and completely mixed conditions. Results obtained from these lab-scale studies show that all three tested solids are very effective in retaining REEs. Over 98% of the CMD REEs that contacted the solids were captured before the solids exhausted their neutralization capacities. We also determined an extraction process using a non-acid, organic ligand extraction solution that can effectively remobilize the retained REEs from the spent solids (over 90%). The REE concentrate (>7.5 wt. % of total REEs) is then formed in an aeration process. The TEP process uses environmentally benign industrial by-products and a naturally-occurring organic ligand to mitigate CMD and recover REEs. Techno-economic analysis (TEA) and life-cycle assessment (LCA) was carried out in the third phase. The engineering-economic costs and net energy, net CO2 emissions, and water and other requirements were investigated to understand the economic and environmental implications of this process. This work uses mass and energy balances from laboratory-scale experiments to estimate the economic costs and environmental impacts. The results suggest that passive treatment systems that use DRWP sludge are preferred over those that use sFGD material, because of lower economic costs ($89,300/yr with a unit cost of $86/gT-REE vs. $89,800/yr, or $278/gT-REE) and improved environmental performance across all indicators from two different impact assessment methods. These differences are largely attributable to the larger capacity of DRWP sludge in the passive treatment application. We envision this TEP process can be integrated with abandoned mine land (AML) reclamation to create an approach that can add economic incentives for AML reclamation, remediate CMD discharge, and eliminate public safety hazards and threats to local environment and ecological systems posed by AMLs. It can restore lands and communities that are adversely impacted by legacy mining.

01 COAL, LIGNITE, AND PEAT↗

Assessment of the Griffin Reactor Multiphysics Application Using the Empire Micro Reactor Design Concept

In late 2019, INL and ANL agreed to jointly develop the reactor physics code named Griffin based on the integration of the two code suites, MAMMOTH/Rattlesnake (INL) and MC2 - 3/PROTEUS (ANL). Griffin is being developed based on the MOOSE framework and MOOSE quality assurance procedures. This decision was made to be able to allow DOE-NE to efficiently invest funding to this area and to provide effective and timely support for existing and potential users; the latter includes industry and government organizations who are developing various types of advanced reactors in the near and long term. Since MAMMOTH/Rattlesnake has been developed based on the MOOSE framework, the INL/ANL Griffin development team agreed to build Griffin beginning with a merger of MAMMOTH and Rattlesnake into a single code and moving forward by implementing capabilities from the PROTEUS suite into Griffin. Moving forward, both ANL and INL efforts are equally invested in the Griffin project, with management support, to provide an advanced reactor multiphysics tool to assist in reactor design, optimization, and safety analysis. Much work remains in moving Griffin forward to migrate PROTEUS capabilities and to optimize performance to meet user needs. The main objective of this work is to assess the current status of Griffin capabilities in terms of performance and accuracy, to determine priorities for PROTEUS migration, and to identify capabilities and features to improve for supporting the code integration effort. For this assessment, the Empire micro reactor problem that was developed in the ARPA-E MEITNER program was selected as an advance reactor concept of interest to the technical community. The Empire reactor problem was expanded from its original incomplete specification to be a small heat-pipe-cooled micro reactor core with ~113 cm radius and 70 cm in height, composed of 18 fuel assemblies, 12 control drums, and beryllium radial and axial reflectors. In the current model, using 5 cm axial reflectors specified in the original Empire assembly model, more than 10% of neutrons leak axially and through the empty center safety hole, as well as through heat pipe channels in fuel assembly elements that extend through the top reflector region. Several calculation models of the core were defined for systematic assessment, including 2-D and 3-D fuel assemblies and whole cores with cylindrical boundaries. Cross sections were generated using Serpent 2, and meshes were produced using the Argonne mesh tool or the INL neutronics meshing tools combined with CUBIT. Cross sections and meshes were converted to the ISOXML and Exodus formats, respectively, so that Griffin and PROTEUS could use consistent data for solving the reactor problems. With the prepared cross sections and meshes, PROTEUS was run first to ensure that all input data were correctly generated and input options in terms of angle, mesh, and energy group were accurately determined. Comparisons against Serpent 2 solutions were made in terms of eigenvalue and pin power. The same calculations and comparisons were then conducted using Griffin. For the fuel assembly and whole core problems, the PROTEUS eigenvalues agreed well with reference Serpent 2 solutions within 100 and 30 pcm, respectively, and pin power differences relative to Serpent 2 were overall less than 2.2% and RMS 0.8% for the whole core models. This indicated that all input data were properly prepared. Using the same data, Griffin was run selecting the SAAF-CFEM SN solver with Legendre-Gaussian quadrature and NDA and DSA for acceleration. It was found that the SAAF-CFEM solver of Griffin required finer meshes to achieve eigenvalue and pin power solutions in good agreement with Serpent 2, consequently requiring more memory requirement and longer computation time. On the other hand, the SPH-Diffusion 2-D core calculations performed using Griffin were able to recover the exact eigenvalue from the reference Serpent 2 solutions, resulting in a pin-power distribution with an RMS of 0.6% and maximum absolute difference of less than 1.4%. The runtimes for SPH-Diffusion for the 2-D core were less than 3 minutes on 40 cores. During this evolution of this evaluation, many updates were made in Griffin by the Griffin development team of INL (focusing on software updates) and ANL (reviewing and supporting software updates) to complete this assessment. Observations from the code assessment are presented in the conclusion section of this report, followed by a discussion of recommendations for future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Non-Intrusive Measurements Using Femtosecond Laser in Atmospheric Pressure Plasma Jet

Femto-second (fs) lasers offer state-of-the art non-intrusive diagnostic capabilities in plasma flow environments. The diagnostics in plasma flows using a fs-laser that are relevant to this proposal include Two-photon Absorption Laser Induced Fluorescence (TALIF) for measuring the number density of an atomic species, and Femto-second Laser Electronic Excitation Tagging (FLEET) for performing velocity measurements in plasma flows. Various aspects of these diagnostics using fs-lasers are still an active domain of research. The TALIF technique has been used to measure the number density of atomic species in various plasma flow environments like atmospheric pressure plasma jet (APPJ) and arc-jet flows. This technique utilizes a fs-laser to excite an atomic species and the subsequent relaxation process through fluorescence is recorded to make measurements. A narrow band-pass filter is used to isolate the fluorescence signal that is in the infra-red regime. However, broad-band background emissions can potentially pass through the filter and induce noise in the TALIF signal. Particularly, arc-jet plasma flows inherently contain background emissions. In such plasma flow environments, in addition to the filter, the background-subtraction technique will play an important role in reducing the noise that is induced in the number density measurements using TALIF. The proposed research focuses on evaluating the effectiveness of a background-subtraction technique by measuring atomic oxygen using a femtosecond (fs) TALIF in APPJ with induced background emissions and comparing measurements obtained without inducing the emissions. Therefore, by quantifying the changes in the signal-to-noise ratio caused by a background-subtraction technique applied for APPJ, it will enable a better understanding of the uncertainties in the background subtracted fs-TALIF measurements for arc-jet plasma flows. In addition to this, FLEET will be used to obtain velocity in the APPJ, and the impact of background emissions will be analyzed. Under this effort, the experimenters and research personnel will gain the technical know-how in the operational and theoretical aspects of the advanced fs-laser-based diagnostic techniques. The knowledge gained in this process will then also directly be applied to a separate ONR funded project at the Aerodynamics Research Center (ARC), UT Arlington, related to arc-jet flow characterization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Involute Working Group – FSI Analysis of Fuel Plates Using Finite Volume and Finite Element Methods

The three involute plate research reactors RHF, HFIR, and FRM II have expressed an interest in using computational software to carry their steady-state safety analysis. Since these tools represent a significant departure from the methods used currently (one-dimensional), the acceptability of the new approach by regulators requires thorough verification and validation of these tools. Therefore, Argonne National Laboratory and the three involute-plate reactors formed an informal group called the Involute Working Group aiming at qualifying computational tools to perform steady-state safety analysis. The present report focuses on a comparison of finite volume and finite element methods to model solids in fluid-structure interaction problems with the goal to estimate the coolant flow-induced fuel plate deflections obtained with the two methods. The finite volume method will be obsoleted in STARCCM+ by the end of 2021, nevertheless, this evaluation is important because the method was used by ANL researchers to model the response of the fuel plates, despite its drawbacks, which are discussed in the report. It was essential to check how those estimates compare to the results obtained with the finite element method that is considered superior for structural analysis. Various geometries, i.e., flat, cylindrical and circle-involute fuel plates, as well as coolant flow speed, were considered. The comparison shows that, independently of the plate geometry, the finite volume method significantly underestimates the deflection as compared to finite element method for coarser meshes. When the discretization is developed as a result of a mesh sensitivity study using finite element method, the result obtained using finite volume method can be a few times smaller than the corresponding finite element method solution. A code-to-code comparison , between STAR-CCM+ and LS-DYNA was included in the analysis. Within the LS-DYNA models, two types of finite element formulations were used: solid and shell finite elements. Mesh sensitivity study showed that both approaches converge to a similar value that was obtained with STAR-CCM+ finite element solver. The evaluation of the computational solvers was extended by adding two benchmark cases from the STAR-CCM+ Verification Suite and presented in the Appendix A. The selected cases are: (1) bending of a cantilever beam under external load, and (2) cylindrical shell deformation analysis, known in the literature as ‘Scordelis-Lo roof’. The problems were solved with finite volume, and finite element methods, and the results confirmed the previously discussed findings. The analysis shows that the finite element solver is superior to the finite volume solver in terms of representation of model geometry and estimating the structural behavior of fuel plates. Depending on the ratio of the load to the flexibility of the plate, the finite volume solver can greatly under- or overestimate the structural response if a very carefully selected mesh is not used.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Minimizing Fraud in the Carbon Offset Market Using Blockchain Technologies

Fraud in the Environmental Benefit Credit (EBC) markets is pervasive. To make matters worse, the cost of creating EBCs is often higher than the market price. Consequently, a method to create, validate, and verify EBCs and their relevance is needed to mitigate fraud. The EBC market has focused on geologic (fossil fuel) CO 2 sequestration projects that are often over budget and behind schedule and has failed to capture the "lowest hanging fruit" EBCs - terrestrial sequestration via the agricultural industry. This project reviews a methodology to attain possibly the least costly EBCs by tracking the reduction of inputs required to grow crops. The use of bio- stimulant products, such as humate, allows a farmer to use less nitrogen without adversely affecting crop yield. Using less nitrogen qualifies for EBCs by reducing nitrous oxide emissions and nitrate runoff from a farmer's field. A blockchain that tracks the bio-stimulant material from source to application provides a link between a tangible (bio-stimulant commodity) and the associated intangible (EBCs) assets. Covert insertion of taggants in the bio-stimulant products creates a unique barcode that allows a product to be digitally tracked from beginning to end. This process (blockchain technology) is so robust, logical, and transparent that it will enhance the value of the associated EBCs by mitigating fraud. It provides a real time method for monetizing the benefits of the material. Substantial amounts of energy are required to produce, transport, and distribute agricultural inputs including fertilizer and water. Intelligent optimization of the use of agricultural inputs can drive meaningful cost savings. Tagging and verification of product application provides a valuable understanding of the dynamics in the water/food energy nexus, a major food security and sustainability issue. As technology in agriculture evolves so to must methods to verify the Enterprise Resource Planning (ERP) potential of innovative solutions. The technology reviewed provides the ability to combine blockchain and taggants ("taggant blockchains") as the engine by which to (1) mitigate fraudulent carbon credits; (2) improve food chain security, and (3) monitor and manage sustainability. The verification of product quality and application is a requirement to validate benefits. Recent upgrades to humic and fulvic quality protocols known as ISO CD 19822 TC134 offers an analytical procedure. This work has been assisted by the Humic Products Trade Association and International Humic Substance Society. In addition, providing proof of application of these products and verification of the correct application of prescriptive humic and bio-stimulant products is required. Individual sources of humate have unique and verifiable characteristics. Additionally, methods for prescription of site- specific agricultural inputs in agricultural fields are available. (See US Patents 734867B2, US 90658633B2.) Finally, a method to assure application rate is required through the use of taggants. Sensors using organic solid to liquid phase change nanoparticles of various types and melting temperatures added to the naturally occurring materials provide a barcode. Over 100 types of nanoparticles exist ensuring numerous possible barcodes to reduce industry fraud. Taggant materials can be collected from soil samples of plant material to validate a blockchain of humic, fulvic and other soil amendment products. Other non-organic materials are also available as taggants; however, the organic tags are biodegradable and safe in the environment allowing for use during differing application timeliness.

54 ENVIRONMENTAL SCIENCES↗

Thermal Process Intensification: Transforming the Way Industry Uses Thermal Process Energy (Workshop Report)

The US Department of Energy’s (DOE’s) Advanced Manufacturing Office held the virtual workshop entitled “Thermal Process Intensification: Transforming the Way Industry Uses Thermal Process Energy” in November and December 2020. The workshop brought together participants from universities/laboratories, industries, equipment manufacturers, technology vendors, nongovernmental organizations, and subject-matter experts to discuss transformative technologies and strategies to substantially improve the performance (e.g., energy productivity, thermal efficiency, reduced greenhouse gas [GHG] emissions, reduced number of process steps) of thermal processing systems in the industrial sector. The US industrial sector accounts for 32% of the nation’s primary energy use (including feedstocks), and refining, chemicals, pulp and paper, iron and steel, and food products represent the top energy-consuming sectors. Thermal processing or process heating represents the largest energy use category; it accounts for 63% of all energy use in manufacturing. Additionally, thermal processing is the largest contributor of carbon dioxide (CO2) generation, resulting from combustion of fuels and process related chemical reactions, such as in the case of cement and lime production. The challenge of achieving net-zero industrial GHG emissions is colossal considering the established industrial base that depends mainly on carbon-based processes and energy sources; the time frame and cost to replace carbon-based energy sources and feedstocks; and the long-term outlook for development with large-scale adaptation of alternative non carbon–based technologies. The goals of the DOE Advanced Manufacturing Office Thermal Process Intensification Workshop were as follows: (1) Identify R&D gaps and opportunities to facilitate transformative improvement in industrial thermal processes beyond current technologies and allow for entirely new methods for processing materials; (2) Gain insight into new and innovative approaches to thermally intensify processes, reduce heat demand, harness waste heat, and use fuels and hydrocarbon feedstocks more efficiently; (3) Identify the R&D pathways to thermal process intensification (TPI) with the highest potential for impact and adoption by the industrial sector; (4) Define areas of research, development, and demonstration (RD&D) activities to accelerate development and application of emerging and transformative technologies to intensify thermal processes in industry. The scope and focus of the workshop were defined to meet these goals. Based on the available data for energy use and GHG emissions, the industries that collectively use more than 80% of the total process heating energy consumption were selected as primary focus areas. The chosen industries were combined into the following four groups based on similarities in their thermal processes: high-temperature metal processing (iron and steel industry, alumina-aluminum industry); high-temperature nonmetal and mineral processing (cement and glass industry); medium- to low-temperature thermal processing (food processing and pulp and paper industry as part of forest products sector); and hydrocarbon processing (petroleum refining and chemical industry). Furthermore, all potential TPI technologies associated with the processes defined were considered as part of the workshop. Different types of TPI technologies possible in industrial in the document.

42 ENGINEERING↗

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

We present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. We inverted waveforms from 192 Global Centroid Moment Tensor earthquakes (MW 5.5-7.0) recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region. Spatial coverage of the available data is highly uneven due to earthquakes clustered along plate boundaries and sparse coverage of open seismic networks in the region. We considered three possible starting models: the SPiRaL global model (Simmons et al., 2021); MEC-1 (Kaviani et al., 2020); and CSEM2.0 (Noe et al., 2023). Because the SPiRaL model provides good fits to the observed waveforms measured by the time-bandwidth product of selected windows in several period bands, provides all the necessary parameters and covers the entire domain we used it for the starting model with the period band 50-100 seconds. Inversion iterations proceeded using time-frequency phase misfits in six stages and 54 total iterations reducing the minimum period to 30 seconds. Our final model, MESWA, provides improved waveform fits compared to the starting model for both the data used in the inversion and an independent validation data set of 66 events. Two metrics of waveform fit (the time-frequency phase misfit used in the optimization and normalized L2 misfit) were both reduced by nearly 60% for both data sets and MESWA provides significantly larger misfit reductions relative to the SPiRaL model than the MEC-1 or CSEM models. We also find that MESWA provides a larger time-bandwidth product of selected windows indicating that more information content of the observed waveforms is explained by MESWA than the other models. Our new model reveals tectonic features imaged by other studies and methods but in a new holistic model of shear and compressional wavespeeds (v S and v P , respectively) with anisotropy covering the crust and uppermost mantle of a larger domain. MESWA has smaller scale-length features and tends to sharpen some features relative to the SPiRaL starting model. Examples include: low crustal v S in the TurkishIranian Plateau, Zagros Mountains, Afghan Central Blocks and Sulaiman Fold Belt; low mantle vSfollowing divergent (Gulf of Aden, Red Sea) and transform (Dead Sea Fault) margins of the Arabian Plate; low and high v S in the mantle beneath the Arabian Shield and Platform, respectively. Low vS is imaged below Cenozoic volcanic centers of the Arabian Peninsula, the so-called Mecca-Madina-Nafud (MMN) Line. Positive anisotropy (v SH > v SV ) is inferred for asthenospheric depths across the region except where up/downwelling may influence fabric alignment (e.g. Afar, Red Sea, Arabian Shield). Elevated vS tracks Makran subduction under southeast Iran. MESWA resembles the SPiRaL model in its long-wavelength structure, but enhances shorter wavelengths features on the order of 200 km and smaller. The resulting model could be used for as a starting model for further improvements, say using waveforms from in-country seismic networks that are not openly available or smaller-scale studies targeting shorter period waveforms. The model also could be used for source characterization and moment tensor inversion to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗

End-Use Savings Shapes Measure Documentation: LED Lighting

Building on the successfully completed effort to calibrate and validate the U.S. Department of Energy's ResStock and ComStock models over the past three years, the objective of this work is to produce national data sets that empower analysts working for federal, state, utility, city, and manufacturer stakeholders to answer a broad range of analysis questions. The goal of this work is to develop energy efficiency, electrification, and demand flexibility end-use load shapes (electricity, gas, propane, or fuel oil) that cover a majority of the high-impact, market-ready (or nearly market-ready) measures. "Measures" refers to energy efficiency variables that can be applied to buildings during modeling. An end-use savings shape is the difference in energy consumption between a baseline building and a building with an energy efficiency, electrification, or demand flexibility measure applied. It results in a timeseries profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each timestep. ComStock is a highly granular, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual subhourly energy consumption of the commercial building stock across the United States. The baseline model intends to represent the U.S. commercial building stock as it existed in 2018. The methodology and results of the baseline model are discussed in the final technical report of the End-Use Load Profiles project. This documentation focuses on a single end-use savings shape measure - LED lighting. The LED lighting measure replaces interior lights with LEDs where applicable. The measure resulted in 1.25 TBtu energy savings across the building stock, with the majority of savings occurring in retail and warehouse buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Exploring the Nature of f-Element Soft Donor Interactions Using Electronically Tunable Azolate Ionic Liquids

This project was undertaken to advance the understanding of how f-elements interact with moderately soft donors, a heavily investigated yet open question which is of prime importance in spent nuclear fuel processing and fundamental inorganic chemistry. During the course of the project, based on exciting results, a stretch goal was developed to study the hydrolysis products of transuranic actinide metals, a somewhat understudied field even with its significance in nuclear fuel cycle and impacts in environmental chemistry. The stretch goal was to take our serendipitous discovery of an easy route to isolation of crystalline multinuclear ƒ-element hydrolysis products, and apply it to gaining a mechanistic understanding of Pu(III/IV) colloid formation. The simplicity of our techniques should lend themselves to the remote handling required for study of many transuranic elements. We developed several methodologies using azolium azolate chemistry to overcome ƒ-element hydrolysis problems that prohibit the isolation of ƒ-element soft donor complexes and to build a crystallographic library of ƒ-element N-donor complexes as a means to understand the fundamental differences between actinide and lanthanide interactions with moderately soft donor ligands. Our next major endeavor will be to transfer this chemistry from 4ƒ elements to transuranic elements, particularly in the study of hydrolysis of Pu(III/IV). While our work is fundamental in nature, applications of the knowledge we are generating should be felt in such diverse fields as catalysis, separations in general, nuclear waste remediation specifically, and many other applications in f-element magnetic and luminescent properties. The potential ramifications of the consistent and robust formation of hydrolysis controlled hexanuclear lanthanide structures are enormous, with future uses being catalyst formation, higher-nuclearity structure synthesis using our hexanuclear motif as a fundamental building block, Pu waste remediation, separations, and many other potential applications resulting from characteristic magnetic and luminescent properties of lanthanide polynuclear structures. Three synthetic methodologies (direct mixing with variable stoichiometries, use of volatile solvent, metathesis) were developed starting with 7 acidic and 6 basic azoles to obtain ionic liquids suitable for f-element coordination. Proton transfer by acidic/basic azole combination led to suitable low melting salts and two cocrystals. Acid/base reaction of azoles with soft-donor permanent cations of ([X 4444 ][OH] (where [X 4444 ] + = tetrabutylammonium [N 4444 ] + or tetrabutylphosphonium [P 4444 ] + ) with weakly acidic azoles including imidazole, 1,2,3-triazole, 1,2,4-triazole, 5-aminotetrazole, 4,5-dicyanoimidazole, and 2-amino-4,5-dicyanoimidazole) revealed several suitable low-melting salts. Metathesis reactions of Na(azolate) were conducted by first using weakly acidic azoles including 4,5-dicyanoimidazole, 2-amino-4,5-dicyanoimidazole, 5-aminotetrazole, and 1,2,4-triazole to form sodium or lithium salts using group(I) hydroxides in methanolic solutions. The best results were obtained by reacting the basic and acidic azoles directly in 1:1 or 3:1 ratios at elevated temperatures. Twenty-two azole mixtures were identified which are either low melting solids or room temperature liquids. Each of the low melting solids was confirmed to be a new solid phase, each of which is being further characterized. The liquids and solids are anticipated to be ILs, eutectics, or partially ionized systems, all of which will be suitable for the dissolution of f-element salts. Five new synthetic methodologies were developed to finding suitable crystallization conditions for f-element complexation with the goal of finding simple one pot reaction syntheses and crystallization strategies that could be used under the demanding conditions of transuranic chemistries. These synthetic methods yield many new crystalline phases which were studied by single crystal X-ray diffraction.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modeling graphene sheet growth and dynamical matrix calculations using molecular dynamics

Molecular dynamics (MD) has been an incredibly useful tool to model physical processes that were synthesized experimentally but not fully understood. MD, through the use of semi-empirical inter-atomic potentials, has allowed understanding of different physical processes in materials science. Yet as well as providing useful insights into materials science, molecular dynamics has a wider range of usability. In this report, I will be detailing how MD can be used to study graphene formation from a carbon liquid which requires high temperatures and pressures. Beyond this, I will describe the usefulness of MD for understanding the physics for phonon transport quantum sensors. To do this, MD was employed to determine the dynamical matrix by treating atoms as coupled oscillators. An accurate understanding of the dynamical matrix of a system is required to calculate the non-equilibrium Green’s function used to describe the phonon transport within phonon wave-guides. I found that, across multiple pressures and temperatures, randomly placed carbon atoms will show evidence of pent-first formation with semi-empirical models. Density functional theory (DFT), on the other hand, was too computationally expensive to use for full scale MD simulations, but we have the possibility of training a machine learned interatomic potential to approximate DFT for carbon in the environments being studied for pent-first graphene sheet formation.

36 MATERIALS SCIENCE↗

Nanostructured Alumina Forming Austenitic Alloy (NAFA) Production Using Mechanical Alloying and High-Temperature Consolidation

Alumina-forming austenitic (AFA) alloys are a promising class of nuclear materials because of their high-temperature oxidation/corrosion resistance and mechanical properties. Unfortunately, these alloys are limited in use for core material applications, and they are specifically limited for use as nuclear fuel cladding because of their high Ni-transmutation and helium generation rate in-service. Improving these alloys through a fine dispersion of oxide precipitates and thus increasing the effective irradiation sink strength of the alloy system may mitigate many of the degradation phenomena expected during alloy deployment. These phenomena include high-temperature helium embrittlement and cavity swelling for lead-cooled fast reactor applications. This work effort uses combination of conventional and advanced manufacturing approaches are being used to fabricate nanostructured AFA (NAFA) materials. As the first objective of this initiative, a conventional AFA was modified using mechanical alloying and extrusion to alter the precipitation characteristics to include a fine dispersion of nanoscale oxides intended to serve as traps for irradiation-induced point defects and transmuted He within the lattice. This analysis compared the efficacy of the conventional mechanical alloying and extrusion approach with the unalloyed AFA consolidated approach using hot isostatic pressing (HIP). It was found that, although the mechanical alloying approach is successful in producing a fine distribution of oxides within the first nanostructured AFA (NAFA-1), the distribution is heterogeneous because of the mild milling parameters used to prevent cold welding of powder to the spherical milling media. The additional dispersion of oxides, coupled with a higher volume fraction of other secondary phases in the NAFA-1, produces higher alloy strengths that range up to 600°C in comparison to the unalloyed HIP AFA, thus exceeding the operating temperature of lead-cooled fast reactors. However, the strength of the NAFA-1 is lower than that of the HIP AFA at 800°C, which is presumed to be a function of increased secondary phases from the nonoptimized AFA chemistry. Future work is planned on a newly procured NAFA-2 chemistry that is more suitable for advanced reactor applications exploring new manufacturing routes.

36 MATERIALS SCIENCE↗