Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “function unit testing”

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 937 records · Page 52

Field Validation of a Grid-Interactive Efficient Building Software Solution

The U.S. General Services Administration's (GSA's) Green Proving Ground (GPG) program, in partnership with the National Laboratory of the Rockies (NLR), completed a field study of a Grid-Interactive Efficient Buildings (GEB) software solution. The study focused on a single testbed facility to test the GEB functionality of the software solution, along with other features. The testbed facility - a courthouse - is a common building type in GSA's vast building portfolio, offering potentially impactful findings on a scalable level. The study evaluated Prescriptive Data's technology, Nantum OS, a connected building operating system ("GEB Solution") which aggregates multiple sources of previously siloed building data and combines that data with external sources, such as weather information or utility signals, into a single integrated platform. A GEB Solution is a type of Energy Management Information System (EMIS). EMIS is defined as a system of devices, data services, and software applications that communicates with any building system or third-party data source to aggregate and transform data into new capabilities to aid in the optimization of energy use at the building, campus, or agency level. This specific GEB Solution is an EMIS with ASO, automated system optimization, offering supervisory control of certain aspects of the Building Automation System (BAS). Multiple features were evaluated including, but not limited to, Continuous Demand Management to avoid setting new monthly kilowatt (kW) peaks, energy efficiency for reduction of kilowatt hours (kWh) and natural gas consumption, and automated demand response (ADR) for purposes of lowering demand during a utility called Demand Response (DR) event. The testbed facility was the Foley Federal Building and US Courthouse ("Foley Federal Building") located in Las Vegas, NV. This is a 209,496 sq. ft. building constructed in the 1960s with major renovations in 2004. The facility was a good candidate due to the large prevalence of office and courthouse spaces in the GSA portfolio of buildings. It also has many features which allow integration into and control of the building and a strong facilities team to assist with the study. Quantitative and qualitative performance objectives were developed using GSA's GPG GEB project template along with input from the vendor and building facility staff; these are outlined in Table 1. The quantitative performance objectives focused on continuous demand management, energy efficiency, and automated demand response. The qualitative performance objectives focused on the ease of installation and commissioning as well as the operability of the GEB solution. Other performance metrics that are reported on include carbon reduction, cost effectiveness, and occupant acceptance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Enhancement of Optical Efficiency of CSP Mirrors for Reducing O&M Cost via Near-Continuous Operation of Self-Cleaning Electrodynamic Screens (EDS). Final Report

Over the past decade, techno-economical advancements in solar energy systems, particularly the improvement in the conversion efficiency of PV modules made with mono-crystalline silicon solar cells from 12% to 20%, as well as the cost reduction in manufacturing by a factor of about 10%, have made it possible to achieve a levelized cost of electricity (LCOE) in PV plants that is comparable to, or less than, the cost of deriving electricity from fossil fuels. Operating PV plants in mid-latitude sun-belt regions, where the solar irradiance level is highest, provides a high annual energy-yield (kWh/kWp) due to two factors: (1) the availability of predictable high solar irradiance throughout the year with the fewest interruptions in solar flux from clouds and rain, and (2) the increased conversion efficiency of crystalline solar cells, as recombination loss has decreased with increased intensity of the sunlight that illuminates the silicon solar cells. Semi-arid and desert regions, however, are plagued by high atmospheric dust concentrations and frequent sand storms. The deposition of a layer of dust on the optical surfaces of solar collectors such as PV modules and concentrating mirrors reduces the transmission efficiency of sunlight that actually reaches the solar cells or receivers, resulting in high energy-yield soiling loss. There are two major cost components to operating a solar plant: (1) installation costs, and (2) operation-and-maintenance (O&M) costs. There is no fuel cost; hence operating a solar plant in a semi-arid or desert region provides high returns on investment if soiling losses are mitigated via efficient cleaning methods and optimized cleaning frequency. If solar collectors are not cleaned, the accumulation of dust layers on solar collectors may cause the operation of such plants in arid regions to become economically unviable. Washing solar collectors with water and detergent, as is most commonly done now, is an efficient method for cleaning. The conventional approach in utility-scale solar plants is to use a large truck with a water tank and pump system for spraying deionized water on the surface of the solar collectors. Robotic cleaning with brushes, used for many solar plants, requires lesser water for cleaning. The water consumed using semi-automated cleaning of PV modules in utility-scale solar plants is approximately 2 liters/m 2 per cleaning cycle. The total optical surface area of the solar collectors in 1 TW-scale solar installation will be more than 3 × 109 m 2 ; hence an enormous amount of water be needed for cleaning. There simply is not enough fresh water in the sun-belt areas of the world for predicted cleaning needs. In solar power plants, the estimated cost of cleaning solar collectors includes expenses related to the equipment used, labor, cost of transportation of water, the energy required for cleaning, plus ancillary costs whereas the cost of water is not considered. The water used is obtained from sources located close to plant sites, unmindful of the environmental and societal impact as the power plants are oftentimes located in regions that face severe drought. This practice is very similar to the cost calculations in deriving the levelized cost of electricity (LCOE) in conventional power plants based on burning fossil fuel such as coal and gas, while disregarding the cost of climate change and health effects. Unless a water-free or low-water cleaning method is established, the expansion of solar plants may lose public support in areas suffering long intervals of drought. The goal of this research project has been the development and application of the Electrodynamic Screen (EDS) as a means for a water-free, scalable cleaning process applicable to solar-power installations, including rooftop applications. We describe here the development of an EDS film-based cleaning process as an emerging method for use on PV modules, parabolic troughs, and heliostats. This report aims to show the feasibility of integrating or retrofitting EDS films onto the optical surfaces of solar collectors (both PV and CSP) while maintaining high transmission or reflection efficiency. The cleaning action provided by the EDS film is an active method to remove dust deposits by electrodynamic force. Current lab-scale prototype EDS films, retrofitted onto solar panels and mirrors, have shown to be capable of maintaining optical transmission or specular-reflection efficiencies higher than 90% of initial values under clean conditions. The optical surfaces of solar collectors laminated with EDS films can remove more than 90% of deposited dust when the EDS is activated for less than two minutes. As an electrodynamic dust removal process, the EDS film-based method is designed primarily for the removal of dust in solar installations located in semi-arid and desert areas, where the atmosphere is often dry and dusty and rainfall is infrequent. While EDS film application minimizes water consumption and facilitates cleaning as frequently as needed, it has limitations in removing contaminants such as soot, organic pollutants deposited as fine films on the surface, and bird droppings. The dust removal efficiency of the EDS is maximum at relative humidity RH is 40 to 50% and decreases at when RH > 65%. Many solar plant sites undergo diurnal and seasonal cycles of high ambient, early afternoon temperatures and an RH that reaches the dew point early in the morning. These variations in atmospheric conditions do not limit the operation of the EDS film. This report presents a brief review of the progress and the potential of EDS film technology for mitigating the impact of dust on solar collectors via water-free cleaning, as well as current technical challenges regarding efficiency and durability. Our experimental data on the performance of EDS films show that: (1) the dust-removal efficiency (DRE) can reach levels higher than 90%, (2) the specular reflectivity (SR) of EDS film-laminated second-surface mirrors reach levels in excess of 90%, (3) the specular reflectivity restoration (SRR) can exceed 90%, (4) the output-power restoration (OPR) of PV modules can exceed 95%, and (5) the optical transmission efficiency (TE) of the EDS films can be greater than 90%. Working with Sandia National Laboratories (NM), Corning Research and Development Corporation (NY), Eastman Kodak (NY), Tomark-Worthen Industries (NH), and EDS Chile SPA (Chile), we have produced EDS film-laminated PV modules and demonstrated their self-cleaning functions without requiring water. We have demonstrated that the operational range of EDS films will cover the expected ambient temperature of solar fields at RH cycling varying from 20 to 95% as long as the EDS films are activated in the RH range 20 to 50%. (Typical solar-field climates in deserts and semi-arid lands often reach near dew point in coastal areas.) Our experiments on the application of hydrophobic-fluorinated nanoparticle coatings on EDS film surfaces show that the EDS operational range can be extended to higher RH levels that approach the dew point. At Eastman Kodak, as one of our industrial partners, we were able to establish a process for manufacturing EDS films using flexographic printing of the electrodes onto transparent polymer films. This process utilizes an existing manufacturing line at Eastman Kodak that allows fabrication of medium-scale EDS films (26 cm × 30 cm). The manufacturing process has the capacity to produce EDS films at high production speeds. The medium-scale EDS films that have been printed at Kodak were evaluated in the lab at Boston University. The EDS films produced at Kodak were laminated at Tomark-Worthen using an industrial scale vacuum laminator to produce EDS film stacks, which can be affixed onto the optical surfaces of PV modules or concentrating mirrors. The EDS film stack consists of the EDS films that have Willow® Glass (WG) which has a thickness of 100 μm as the front surface. The WG sheets obtained from Corning Research and Development Corporation are customized in size and shape to cover the active area of the EDS films. The back surface of the EDS film stack is integrated onto the optical surface of the solar panel or mirror using optically clear adhesive (OCA) films or silicone adhesives. The OCA films (thickness 25 μm) are produced by 3M. The EDS film stacks have the architecture: WG/OCA/EDS Film/OCA/ over PV module or solar mirror. The power-supply units needed for activating the electrodes of the EDS film were designed and produced at Boston University. These power supply units provide three-phase, 1.2 kV voltage pulses at a very low current (micro-ampere) level and at a low frequency (≈ 5 Hz). The voltage pulses are applied to the electrodes in a sequence such that the train of pulses resembles a unidirectional traveling wave of electrical field on the surface of the EDS film. The dust particles on the surface become charged electrostatically and are levitated by the Coulomb force. The lateral sweeping action of the traveling electric field created by the three-phase voltages pulses then sweeps the dust off the surface. The energy consumed by the EDS electrodes is less than 0.2 Wh/m 2 /cleaning cycle, enabling energy-efficient restoration of output power (OPR) of PV module or specular reflectivity restoration (SRR) for solar mirrors. The EDS system consists of (1) an EDS film stack laminated onto the solar collectors, (2) connection of the EDS film to its power supply unit and (3) Interconnection of the power supply to PV modules or solar mirrors. Design, construction and assessment of field-testing units that have EDS stack laminated PV modules for evaluating the performance of the EDS films in solar fields is being carried out at BU. Our progress under this project, aimed at the advancement of EDS film technology, has reached DOE Technology Readiness Level (TRL) 6. Based on the extensive laboratory evaluations and limited field trials, as well as contacts with potential users, we believe that the technology has reached its commercial stage. We are conducting a cost analysis using the National Renewable Energy Laboratory (NREL) System Advisory Model (SAM) and are preparing for field trials of EDS films in different solar fields in the US, Chile, India and in the Middle East. A brief description of EDS film performance, construction and testing of the field-test unit, autonomous operation of the field-test unit for evaluating EDS performance in increasing energy yield, associated revenue savings, and water conservation is presented.

14 SOLAR ENERGY↗

FUN-BioCROP model with litter decomposition parameters derived from the LIDET dataset

This repository contains the code and data necessary to run the FUN-BioCROP (Fixation and Uptake of Nitrogen-Bioenergy Carbon, Rhizosphere, Organisms, and Protection) model with litter decomposition parameters derived from a modified Monte Carlo simulation that used the Long-term Intersite Decomposition Experiment Team dataset. Related publication:Juice, S.M., Ridgeway, J.R., Hartman, M.D., Parton, W.J., Berardi, D.M., Sulman, B.N., Allen, K.E., & Brzostek, E.R. Reparameterizing litter decomposition using a simplified Monte Carlo method improves litter decay simulated by a microbial model and alters bioenergy soil carbon estimates. Description of Files: FUNBioCROP_LIDET Study.Rmd R code with FUN-BioCROP model that can be run with 10 different sets of parameters for litter decomposition (Baseline, LIDET, or eight other best parameter sets identified in the modified Monte Carlo simulation. CORPSE Functions_Bioenergy_V2.R Code with CORPSE model functions, called by FUNBioCROP_LIDET Study.Rmd Model Input Data: bulk.csv, bulk_till.csv, rhizo.csv, rhizo_till.csv, litter.csv Initial C and N (kg C or N/m2) pool values for each soil compartment, final values from spin up. All five files have the same columns: (Column - Description - Units) uFastC - Unprotected fast decomposing carbon - kg carbon/m2 uSlowC - Unprotected slow decomposing carbon - kg carbon/m2 uNecroC - Unprotected necromass carbon - kg carbon/m2 pFastC - Protected fast decomposing carbon - kg carbon/m2 pSlowC - Protected slow decomposing carbon - kg carbon/m2 pNecroC - Protected necromass carbon - kg carbon/m2 livingMicrobeC - Carbon in living microbial biomass - kg carbon/m2 uFastN - Unprotected fast decomposing nitrogen - kg nitrogen/m2 uSlowN - Unprotected slow decomposing nitrogen - kg nitrogen/m2 uNecroN - Unprotected necromass nitrogen - kg nitrogen/m2 pFastN - Protected fast decomposing nitrogen - kg nitrogen/m2 pSlowN - Protected slow decomposing nitrogen - kg nitrogen/m2 pNecroN - Protected necromass nitrogen - kg nitrogen/m2 inorganicN - Inorganic nitrogen - kg nitrogen/m2 CO2 - Carbon in carbon dioxide - kg carbon/m2 livingMicrobeN - Nitrogen in living microbial biomass - kg nitrogen/m2 Model Input Data: FluxTower_AvgSoilT.csv: Average daily soil temperature (oC) at 10 cm depth at University of Illinois Urbana-Champaign (UIUC) Energy Farm flux tower from 7/2008-3/2016. (One year of averaged data) Model Input Data: FluxTower_AvgSoilVWC.csv: Average daily soil volumetric water content (VWC) at 10 cm depth at UIUC Energy Farm flux tower from 7/2008-3/2016. (One year of averaged data) Model Input Data: input_CCS_LIDET Study.csv: This file has daily data to run FUN-BioCROP (Column - Description - Units): yr - calendar year - year doy - day of year (1 to 365) (no leap year) - day anpp - aboveground NPP (DayCent) - kg C/m2/day bnpp - belowground NPP (DayCent) - kg C/m2/day aglivc - live aboveground biomass carbon (DayCent) - kg C/m2 bglivcj - live juvenile fine root biomass carbon (DayCent) - kg C/m2 bglivcm - live mature fine root biomass carbon (DayCent) - kg C/m2 aglivn - live aboveground biomass nitrogen (DayCent) - kg N/m2 bglivnj - live juvenile fine root biomass nitrogen (DayCent) - kg N/m2 bglivnm - live mature fine root biomass nitrogen (DayCent) - kg N/m2 nyr - simulation year - year cult - indicates a cultivation event (0 or 1) crop - indicates a new crop (0 or 1) fert - indicates a fertilizer event (0 or 1) frst - indicates the first day of the growing season (0 or 1) harv - indicates a harvest event (0 or 1) last - indicates the end of the growing season (0 or 1) croptype - crop type (0=none; 1=alfalfa; 2=corn; 3=grass clover pasture; 4=soybean; 5=wheat) cropsrl - crop specific root length - mm/g root cultrhizmix - fraction of rhizosphere mixed with bulk soil during cultivation (0.0-1.0) - fraction cultlitmix - fraction of litter mixed with bulk soil during cultivation (0.0-1.0) - fraction harvremov - fraction of above ground biomass removed during harvest (0.0-1.0) - fraction fertamt - fertilization amount - g N/m2 lifehist - plant life history (0 = annual, 1 = perennial) froot_turnover_c - amount of C in fine root turnover - kg C/m2 froot_turnover_n - amount of N in fine root turnover - kg N/m2 agrd_turnover_c - amount of C in aboveground biomass turnover - kg C/m2 agrd_turnover_n - amount of N in aboveground biomass turnover - kg N/m2 leaf_litter_fastfrac - Fast decomposing fraction of leaf litter (0.0-1.0) - fraction root_litter_fastfrac - Fast decomposing fraction of root litter (0.0-1.0) - fraction root_diameter - root diameter - mm root_length - root length - mm root/m2 rhizo_frac - fraction of total soil volume that is rhizosphere (0.0 - 1.0) - fraction date - date in format YYYY-MM-DD Instructions: Save the model code ("FUN-BioCROP_LIDET Study.Rmd") and accompanything files (data streams and CORPSE function code) in the same folder. In model code "Chunk 3: Load CORPSE Data Streams" set the working directory (setwd) to the folder with the files saved in step #1. In "Chunk 5: Define LIDET parameter sets" select the litter decomposition parameter set to be used in the run, and comment out all other sets. If changing any parameter values, edit them in "Chunk 6: Load parameters." Run all chunks up to and including "Chunk 10: Prepare Data for Export." In "Chunk 11: Export Output Data" edit data frames for export and filenames, as necessary. "Chunk 12: Graph Total Soil C" makes a figure of C remaining over the model run period. Description of each model chunk (in file FUN-BioCROP_LIDET Study.Rmd): Chunk 1: Remove all functions, clear memory. Removes all functions from R environment, clears the memory. Chunk 2: Load Packages. Loads packages necessary to run the code. Chunk 3: Load CORPSE Data Streams. Sets the working directory and loads the data files necessary to run CORPSE. Chunk 4: Load CORPSE Functions. Loads the R script with CORPSE functions from the working directory, "CORPSE Functions_Bioenergy_V2.R". Chunk 5: Define LIDET parameter sets. Has ten different parameter sets for litter decomposition tested in this study: Baseline parameters, LIDET parameters, and the other 8 best performing parameter sets identified in the modified Monte Carlo. To run the model, all but one parameter set must be commented out. Chunk 6: Load Parameters. Loads all fixed parameters to run the model. Data frame with definitions of parameters is in the CORPSE function script "CORPSE Functions_Bioenergy_V2.R" Chunk 7: Prepare Data Streams. Takes data streams loaded in Chunk 3 and puts them in the format necessary to run the model. The model is coded to run at least two sites at a time, so if only one site is being run it must be run in duplicate. Individual data tables of daily values are created in this chunk from the input data file. Chunk 8: Set Initial Conditions. Creates data tables of soil C and N pools for each soil compartment (rhizo_till, rhizo, bulk_till, bulk, litter) and loads initial values into the data tables. Creates lists for each soil compartment to hold model output. Chunk 9: Load FUN Data and Set Up Matrices. Uses DayCent data to calculate FUN input data: root and leaf N demand, total N demand, plant CN, leaf N available for retranslocation, and litter production. Creates matrices for FUN model outputs. Chunk 10: Run Model. Runs the model. Chunk 11: Prepare Data for Export. Combines data from each day saved as lists into data frames for each soil compartment. Adds values from all soil compartments together to calculate total soil values, creates separate data frames for each soil C and N pool (e.g., protected slow C) for the total soil value. Adds different C and N pools together to calculate total soil C and N for all layers. Creates data frame of ratio of protected to unprotected SOC. Organizes FUN data for export. Chunk 12: Export Results. Exports CSV files of model results to the working directory. Chunk 13: Graph Total Soil C. Makes figure of C remaining over time. Related Links: Original FUN-BioCROP model: https://github.com/BrzostekEcologyLab/FUN-BioCROP LIDET dataset: https://andlter.forestry.oregonstate.edu/data/abstract.aspx?dbcode=TD023

Juice, Stephanie↗

Pilot Testing of a Highly Efficient Pre-combustion Sorbent-based Carbon Capture System

TDA developed and demonstrated a highly efficient pre-combustion carbon capture system. The overall objective of this work was to develop a new sorbent-based pre-combustion carbon capture technology for Integrated Gasification Combined Cycle (IGCC) power plants. In this project our goal was to demonstrate the techno-economic viability of the new technology by 1) demonstrating it in large-scale slipstream tests, and 2) carrying out a high fidelity engineering and cost analysis. TDA’s process used an advanced physical adsorbent that selectively removes CO 2 from coal-derived synthesis gas above the dew point of the gas at temperatures as high as 300°C. The sorbent consists of a mesoporous carbon whose surface was grafted with functional groups that remove CO 2 via a well-known acid-base interaction. As documented in bench-scale experiments and field tests with actual coal gas, the sorbent achieved a very high capacity for CO 2 at temperatures as high as 300°C. The sorbent bound CO 2 more strongly than common physical adsorbents, providing the chemical potential needed for the high temperature operation. However, because CO 2 does not form a true covalent bond with the surface sites (as is the case with chemical absorbents), the sorbent regeneration could be carried out with only a very small energy input. The heat input to regenerate our sorbent was only 4.9 kcal per mol of CO 2 , which is much lower than that for chemical absorbents (e.g., 29.9 kcal/mol CO 2 for sodium carbonate) and was similar to the requirements of physical solvents (e.g., 4 kcal/mol CO 2 for Selexol TM ). Because the sorbent operates above the dew point of the synthesis gas (unlike the Selexol TM process), a higher power cycle efficiency can be achieved. With previous DOE/NETL funding (Contract No. DE-FE-0000469), we demonstrated the techno-economic viability of the technology in bench-scale tests and slipstream demonstrations at the National Carbon Capture Center (NCCC), Wilsonville, Alabama and Wabash River IGCC plant in Terra Haute, Indiana. We demonstrated a stable working CO 2 capacity for over 11,650 cycles with simulated synthesis gas. We also evaluated its performance with actual synthesis gas in two test campaigns at the Wabash River IGCC Plant, Terre Haute, IN and the National Carbon Capture Center (NCCC), Wilsonville, AL. The slipstream tests clearly showed that the actual coal gas constituents and potential contaminants (e.g. trace metals, halides, tars) had no effect on the sorbent’s ability to remove CO 2 (the same sorbent beds were used in both field tests with no sign of deactivation for 2,000 cycles with over 26,750 SCF of gas treated). As expected, due to the high temperature CO 2 removal capability and low energy needed to regenerate the sorbent, the power cycle efficiency with our process was greater than 34% on a higher heating value (HHV) basis; in comparison, the same IGCC plant equipped with the Selexol TM solvent for carbon capture can only achieve 31.4% HHV efficiency. The capital cost for an IGCC system with TDA’s process is estimated as $2,417/kW e , which is 12% lower than that of the IGCC/ Selexol TM process. The levelized cost for electricity including the transport, storage and monitoring (TS&M) cost for CO 2 was calculated as $\$ $92.9/MWh (lowest reported to our knowledge), which is much better than the $105.2/MWh estimated for the IGCC/ Selexol TM process. In this project (DE-FE0013105), TDA Research, in collaboration with our partners Gas Technology Institute (GTI), Illinois Clean Coal Institute (ICCI), University of California, Irvine (UCI), University of Alberta (UOA), Siemens, NCCC and Sinopec advanced the technical maturity of the technology; scaling it up by a factor of 100. We optimized the reactor design using computational fluid dynamics (CFD); using adsorption modeling we improved the pressure swing adsorption (PSA) cycle sequence. We carried out two field test campaigns with a fully-equipped 0.1 MW e prototype unit (for a total of 844 hours) using actual synthesis gas to prove the viability of the new technology. A successful 30 day (707 hrs) evaluation was completed at NCCC under air blown gasification conditions. We demonstrated 97.3% carbon capture at 1,500 SLPM, 93% carbon capture at 1,800 SLPM, and 90% carbon capture at 2,100 SLPM in the NCCC tests. We also demonstrated the system for 137 hours at a Sinopec petrochemical plant under oxygen blown gasification, demonstrating 86% carbon capture at 2,660 SLPM. In collaboration with University of California, Irvine (UCI), we completed a techno-economic analysis (TEA) for TDA’s warm gas cleanup technology integrated to IGCC power plant. The net plant efficiencies (on a coal HHV basis) for the warm gas cleanup cases were estimated to be 34.0% for E-GasTM gasifier, 34.4% for GE gasifier, 33.4 for the Shell gasifier and 34.2 for the TRIG TM gasifier (Cases 2, 4, 6 and 8 in this study) with a catalytic combustor for CO 2 purification, which are significantly higher than those for the Cold Gas Case, or an increase of as much as 12% in the heat rate for Case 2, 6% for Case 4, 9% for Case 6, and 9% for Case 8. The 1st year cost of electricity with the transport, storage and monitoring (TS&M) costs for the CO 2 included was $\$ $129.2/MWh for the E-GasTM gasifier Warm Gas Cleanup Case, $\$ $131.9/MWh for the GE gasifier Warm Gas Cleanup Case, $\$ $146.8/MWh for the Shell Gasifier Warm Gas Cleanup Case, and $\$ $129.9/MWh for the TRIG TM gasifier Warm Gas Cleanup Case. For comparison, the costs for the baseline Cold Gas CO 2 removal with Selexol for the different gasifiers were: $\$ $146.6/MWh for the E-Gas TM gasifier, $\$ $142.2/MWh for the GE gasifier, $\$ $159.0/MWh for the Shell gasifier and $\$ $144.3/MWh for the TRIG TM gasifier. In summary, the costs for our system were 7 to 12% lower than the corresponding Cold Gas Cleanup cases. The results of this techno-economic analysis suggested that TDA’s high temperature PSA-based Warm Gas Clean-up Technology can make a substantial improvement in the IGCC plant thermal performance for achieving near zero CO 2 emissions for E-Gas TM , GE, Shell and TRIG TM gasifier based IGCC power plants. The capital expenses were estimated to be lower than that of Selexol’s™. Taken together, the higher net plant efficiency and lower capital and operating costs resulted in substantial reduction in the cost of carbon capture for the IGCC plant equipped with TDA’s high temperature PSA-based carbon capture system. Finally, in collaboration with Gas Technology Institute (GTI) we completed the environmental health and safety assessment for TDA’s warm gas carbon capture technology.

01 COAL, LIGNITE, AND PEAT↗

NRAP-Open-IAM: NRAP Open Source Integrated Assessment Model

Note: This is the last version (a2.6.1) of NRAP-Open-IAM released during NRAP Phase II in 2022. The latest version of NRAP-Open-IAM is available here: https://edx.netl.doe.gov/dataset/phase-iii-nrap-open-iam NRAP-Open-IAM is an open-source software product that enables quantification of containment effectiveness and leakage risk at storage sites in the context of system uncertainties and variability. NRAP-Open-IAM represents the next-generation in a line of systems-based computational models developed for quantitative geological carbon storage (GCS) risk assessment. The model comprises a set of reduced-order and analytical models of various components of the GCS system, potential leakage pathways, receptors of concern including impact to groundwater resources and the atmosphere, a framework to support stochastic simulation, time stepping, uncertainty quantification, other analytical functionality for scenario and risk-performance evaluation, and a basic graphical user interface to support scenario development, data input simulation definition, and basic post-processing and results display. As the NRAP Open-IAM functionality continues to evolve, we continue to add to its capability to develop quantitative, probabilistic, and time-dependent profiles of the evolution of risk at a GCS site and evaluate the influence of uncertain parameters on uncertainty in predicted risk. It can be used to quantify the dynamics of reservoir saturation plume and pressure-affected area, for evaluation of the area of potential groundwater impact (i.e., Area of Review) and monitoring requirements to support cost and regulatory analysis, and for consideration of different post-injection site care and closure scenarios. This submission contains the current version of NRAP-Open-IAM available for evaluation and testing. To use the NRAP-Open-IAM, download the source code (https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/4c24a3da-3b40-4ffe-9892-c807ae9f8760) then open the NRAP-Open-IAM user's guide (https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/8b27335a-343c-4836-b8a3-3ad0bdc9e669) to read more about the tool. Installation instructions for Windows, Mac, and Linux can be found in the "installers" folder of the extracted NRAP-Open-IAM folder and describe setup of environment (e.g., Python libraries) needed for proper work of the tool. Test of installation can be done by running "python openiam_setup_tests.py" in the "setup" folder. The installation test also runs a test suite to see if the NRAP-Open-IAM has been installed correctly. To run the test suite separately, run "python iam_test.py" in the "test" folder. User's guide: https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/8b27335a-343c-4836-b8a3-3ad0bdc9e669 Developer's guide: https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/3bc6ee7d-609d-4eb6-80ba-fa6130ee0313 Reservoir simulation data used in some examples distributed with NRAP-Open-IAM: - Kimberlina: https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/eb62cece-61b2-4037-9b6d-32407dde2ab8 - Kimberlina (compartmentalized): https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/366f9530-3b32-4b84-affe-ab2df1d9a8b5 - FutureGen 2.0: https://edx.netl.doe.gov/dataset/futuregen-2-0-1008-simulation-reservoir-lookup-table NRAP-Open-IAM GitLab repository: https://gitlab.com/NRAP/OpenIAM Related publications: - Bacon, D., Yonkofski, C., Brown, C., Demirkanli, D. and Whiting, J., 2019. Risk-based post injection site care and monitoring for commercial-scale carbon storage: Reevaluation of the FutureGen 2.0 site using NRAP-Open-IAM and DREAM. International Journal of Greenhouse Gas Control 90: 102784. - Bacon, D. Demirkanli, D., and White, S., 2020. Probabilistic risk-based Area of Review (AoR) determination for a deep-saline carbon storage site. International Journal of Greenhouse Gas Control 102: 103153. - Harp, D., Oldenburg, C., and Pawar, R., 2019. A metric for evaluating conformance robustness during geologic CO2 sequestration operations. International Journal of Greenhouse Gas Control 85: 100-108. - Lackey, G., Vasylkivska, V., Huerta, N., King, S., and Dilmore, R., 2019. Managing well leakage risks at a geologic carbon storage site with many wells, International Journal of Greenhouse Gas Control, 88 :182-194. - Vasylkivska, V., Dilmore, R., Lackey, G., Zhang, Y., King, S., Bacon, D., Chen, B., Mansoor, K., and Harp, D., 2021. NRAP-Open-IAM: A flexible open-source integrated assessment model for geologic carbon storage risk assessment and management, Environmental Modelling & Software, 143: 105114. Presentations: - Chen, B., Harp, D., and Pawar, R., A data assimilation approach (ES-MDA) coupling with NRAP-Open-IAM for quantifying uncertainty reduction in geological CO2 sequestration. AGUFM 2019: T44A-02. - Chen, B., and Harp, D., Improving risk analysis precision for geologic CO2 sequestration by quantifying the uncertainty reduction before and after acquiring monitoring data. 14th Greenhouse Gas Control Technologies Conference, Melbourne, Australia, 2018, pp. 21-26. - Harp, D., National Risk Assessment Partnership Task 2: Containment Assurance. No. LA-UR-19-28654, Los Alamos National Laboratory (LANL), Los Alamos, NM (United States), 2019. - Vasylkivska, V., King, S., Bacon, D., Harp, D., Chen, B., Mansoor, K., Onishi, T., Yang, Y., Zhang, Y., and Keating, E., NRAP-Open-IAM: An open-source integrated assessment model, poster, Mastering the Subsurface Through Technology Innovation, Partnerships and Collaboration: Carbon Storage and Oil and Natural Gas Technologies Review Meeting, Pittsburgh, PA, August 13-16, 2018. - Vasylkivska, V., Lackey, G., King, S., Wentworth, A., Huerta, N., Creason, C., DiGiulio, J., Yang, Y., and Dilmore, R., Long-term risk analysis of a geologic CO2 storage project during the post-injection period, SIAM Conference on Computational Science and Engineering, Spokane, WA, February 25-March 1, 2019. - Vasylkivska, V., Overview of the NRAP-Open-IAM tool for carbon storage (beta release), 2019 Annual NRAP Tool Users Meeting, Pittsburgh, PA, August 27, 2019. - Vasylkivska, V., Bacon, D., Chen, B., Dilmore, R., Harp, D., King, S., Lackey, G., Lindner, E., Liu, G., Mansoor, K. and Zhang, Y., NRAP-Open-IAM: A new, open-source code for integrated assessment of geologic carbon storage containment effectiveness and leakage risk, poster, American Geophysical Union Fall Meeting 2020 (virtual meeting), December 2020. - Vasylkivska, V., NRAP open-source integrated assessment model and relevant application, oral presentation, NRAP workshop "NRAP Tools for Geologic Carbon Storage Risk-Based Decision Making" held in conjunction with Groundwater Protection Council (GWPC) 2021 Annual Forum (virtual meeting), Salt Lake City, UT, September 2021. - Vasylkivska, V., NRAP-Open-IAM: open-source integrated assessment model, digital poster/demonstration, software demonstration session, 2022 Carbon Management Project Review Meeting, August 16, 2022

AoR↗

Selective Shielding Bone Marrow: An Approach to Create Personal Protective Equipment for Gamma Radiation

Many casualties of the Hiroshima and Nagasaki atomic bombs, and many of the firefighters who first responded to the Chernobyl nuclear power plant accident, became ill with acute radiation syndrome (ARS). Most of the people who do not recover from ARS will die within a few weeks to a few months after exposure, with the primary cause of death being the destruction of the body's bone marrow (BM). This is exemplified by the Chernobyl firefighters, the majority of whom received whole body doses between 0.8 and 10 Gy, and the United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR) final report on the disaster concluded that the underlying cause of death of all 28 firefighters who succumbed to ARS was from bone marrow failure. BM is comprised of hematopoietic stem cells (HSCs) which are responsible for the constant renewal of blood cells. Due to their high rate of proliferation, HSCs are especially vulnerable to ionizing radiation, but are endowed with remarkable regenerative potential. Owing to their central role in blood production, lethal irradiation of HSCs leads to death from severe anemia, infection and internal bleeding. This relationship between high doses of radiation and HSC apoptosis has led to the use of HSC transplantation as a life-saving intervention in cases of acute exposure. This is demonstrable in medical practice, where thousands of individuals have undergone supra-lethal Total Body Irradiation (TBI) for purposes of cancer therapy, and were rescued by well-matched BM transplantation. Indeed, life-threatening damage may be reversed by BM transplantation in individuals receiving doses of radiation of up to 10 Gy. At doses > 10 Gy, damage to gastrointestinal (GI) tissue may become a limiting factor in survival. Remarkably, in the common procedure of BM transplantation, the number of HSCs extracted from a single active BM site containing less than 5% of the donor's BM tissue is sufficient to support the complete reconstitution of the HSC compartment in a lethally irradiated recipient. This capacity of BM to expand and replenish is due to the high regenerative potential of the HSCs it harbors. In contrast to mature blood cells, which are dispersed throughout the body, HSCs are confined to the bones, allowing for effective targeted shielding. This has been confirmed in several animal models, where sparing the HSC-rich area of a subject receiving otherwise full-body irradiation is sufficient to support hematopoietic functions and allow survival. Approximately 48% of the human body's active BM is contained within the lumbar vertebrae and the pelvic region. This high concentration of BM makes the pelvic region an ideal area to shield for preventing the serious hematopoietic effects of radiation exposure. Stem-cell rich tissues like active BM are also more sensitive to the stochastic effects (i.e. cancer) of radiation exposure. The high amount of BM in the pelvic region, combined with the presence of the sensitive gastrointestinal system and female gonads in this same area of the body allows for the possibility of significantly reducing stochastic effects of radiation by limiting radiation exposure to this vital area. Thus, shielding this region holds great promise for both acute and protracted exposures. Thus, a specialized radiation shielding device for the protection of active BM concentrations is presented here. This belt-like selective shielding device focuses on the protection of BM that is present in the pelvis and hip bones, protecting the medullary volume from which BM is commonly extracted for transplantation (i.e. the iliac bones) while allowing relatively unhindered movement of the wearer. In order to optimize the use of shielding materials towards the protection of active BM, the shielding uniquely brings into account the natural shielding properties of human tissue, with thickness being inversely related to the thickness and radiation attenuation of the underlying tissue at each point surrounding the area being protected. The findings resulting from both experimental and simulation testing of this selective shielding device are presented here.

61 RADIATION PROTECTION AND DOSIMETRY↗

Multimetallic Layered Composites (MMLCs) for Rapid, Economical Advanced Reactor Deployment (Final Report)

This project focused on the development of multi-metallic layered composites (MMLCs) for advanced fission reactor technologies. There are many instances where one alloy or material simply cannot meet all the demands thrown at it by a reactor system, or cannot allow it to perform as strongly as one would like. Instead of focusing all our effort on developing one perfect alloy, we seek to leverage the design principle of “separation of functionality,” used in many other arenas in design, to boost performance beyond single alloys alone. One illustrative example shows the power of this approach for molten salt-cooled reactors: A three meter tall, three meter diameter reactor vessel made of Incoloy 800 was quoted at $\$$500k in 2018. A Hastelloy N vessel was quoted at $\$$5M. An MMLC vessel, in which a layer of Hastelloy N would be weld-overlaid onto Incoloy 800, was quoted at $\$$700k, and it would achieve the same performance. The potential economic gains of leveraging this approach are therefore substantial. At a minimum, each MMLC would contain one core structural layer and one coolant-facing corrosion-resistant layer. Sometimes, MMLCs required buffer layers, as the structural and corrosion-resistant layers were metallurgically incompatible. In other words, they didn’t always play nice, thus separating layers compatible with both functioned as intermediaries to keep the composite together. However, in doing so we inevitably produce new interfaces, where new issues can arise. Therefore, this project focused on what happens at these interfaces from a combination of high temperatures, irradiation, corrosion, and time. After all, a reactor makes money when it is operating, and outages of any kind erode its economic viability. First, we set out to experimentally prove that MMLCs for at least two advanced reactor systems can be made, today, in US domestic facilities. In this respect we were successful – one MMLC (a Ni-201/Incoloy 800H composite) was successfully made and drawn into two-inch coolant piping. Others were attempted, though new issues relating to cracking in vanadium layers for one and radiation damage performance of the corrosion-resistant layer in another prevented us from moving further in those specific arenas – these are engineering problems which deserve continued focus after this project. Additional experimental work focused on long-term corrosion testing of the outermost layers of the salt-cooled and liquid lead-cooled MMLC concepts, which would then be fed into predictions of how long the MMLCs could last. Next, computational (thermodynamics and atomistic) simulation studies studied how much we expect the interfaces to “blend,” due to the mixing action of neutron irradiation. This eats into both the margin for the structural layer of each MMLC, as dilution from the corrosion-resistant layer into the structural layer would decrease the total load-bearing capacity of an MMLC of finite size. On the other hand, dilution of the corrosion-resistant layer into the structural layer further reduced the margin of corrodible material, reducing the lifetime of the MMLC or necessitating extra thickness to be imparted to the MMLC to meet its functional requirements. Work here focused on irradiation-induced segregation to predict new phases which may embrittle the MMLCs, as well as quantifying irradiation-induced mixing at each interface. The results showed that mixing is expected, but it is both steady and therefore predictable, and not lifetime-limiting for most MMLC concepts – it simply has to be accounted for in calculations of reactor performance when utilizing an MMLC. Then, full-core simulations using the experimentally-derived corrosion data, the computationally discovered irradiation-induced mixing data (partially validated by experiment), and existing, benchmarked core designs for large and small sized reactor concepts (one salt-cooled, one lead-cooled) were conducted to quantify any expansion of reactor operating envelopes achieved by utilizing these MMLCs. This new framework, called REX (Reactor Envelope Expansion), incorporates a combination of core neutronics, thermal hydraulics, and the material performance data derived from this project to see how using an MMLC expands advanced fission reactor operating envelopes. It was discovered that in some cases, MMLC utilization does indeed increase the maximum operating temperatures and cycle lengths of reactor concepts, while in other cases it does not. Finally, our tech-to-market (T2M) strategy was not necessarily to create specific embodiments of MMLCs for immediate sale (because getting into the nuclear market is incredibly slow and laden with regulation, this is a long-term goal), but rather immediate stimulation of US industry using the design approach of MMLCs derived from this project. In this respect we were successful, as one of the PhD students funded on this project co-founded Allium Engineering, Inc., which created a stainless steel / low-alloy steel MMLC to function as chloride corrosion-resistant rebar for embedding into concrete structures. Allium Engineering continues to be successful, having recently opened their first factory as of this writing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Investigations into the Ternary NaF-KF-UF4 Salt System – Phase A

A knowledge gap exists in the data and understanding of fresh fuel salt and irradiated multicomponent fuel salt systems thermophysical properties. Quantifying these properties is necessary for the design and construction of test reactors, as well as the licensing of future commercial molten-salt reactors. To facilitate thermal property determination on a proposed fuel salt composition for Seaborg Technologies, several samples containing depleted uranium tetrafluoride (UF4), sodium fluoride (NaF), and potassium fluoride (KF) were blended, and a melt temperature analysis was performed. From the melting temperature analysis, it was determined that sample Seaborg-7, a ternary salt composition of 26.4UF4-24.7KF-48.9NaF (mol%), was very near a ternary eutectic point. Therefore, thermal properties such as melting temperature, salt stability, density, heat capacity, thermal diffusivity, and viscosity were experimentally determined on the Seaborg-7 salt. These measurements document the baseline properties of fresh fuel salt as a function of temperature, where future experiments on irradiated fuel salt will provide a holistic perspective on the change of thermophysical properties during reactor operations. Several precision instruments were used to collect property data, and instrument calibrations and data collection were performed and documented in a standardized and reproducible manner with meticulous detail. This process ensured that the measurement procedures and resulting data can readily be duplicated elsewhere. The Seaborg-7 salt was shown to be stable at temperatures up to 900°C, as no mass change was observed upon repeated heating and cooling. The peak melting temperature was determined to be 547°C (557°C endset). The enthalpy of fusion (??H?_fus^o) was determined to be 167.5 ± 2.7 J/g while the enthalpy of crystallization (??H?_c^o) was determined to be -147.8 ± 13.3 J/g. In addition to the eutectic melting peak, upon heating, several pre eutectic peaks were observed, occurring at 470°C (onset) and 499°C (peak). Specific heat capacity measurements showed a slightly increasing trend with respect to temperature in the solid phase, while the liquid-specific heat capacity showed a somewhat flat trend with an average value of 106.1 ± 1.24 J/mol·K between 600 to 800°C. Three independent trials using the Seaborg-7 salt determined the density to be ?(T) = 4.908 – 0.000363·T(°C), validated between 32 to 200°C, and ?(T) = 4.808 – 0.00113·T(°C), validated between ~575 to 850°C. Thermal diffusivity was determined for the liquid state and is represented by the linear equation y = 0.1581 + 0.000207·T(°C) between 550 to 850°C. The viscosity was determined from 600 to 800°C and is represented by the exponential fit equation, ? (mPa·s) = 736.58e^(-0.006·T(°C)). This report documents the conclusion of fuel salt thermophysical property measurements for the Seaborg SPP, Phase A project.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Rapid Bayesian High Entropy Alloy Designs Fabricated via Wire Arc Additive Manufacturing

Purpose: This project seeks to demonstrate a new high-throughput (rapid) alloy design technique applied to creating new high entropy alloys (HEAs) for extreme environments. High entropy alloys shift the design paradigm from being focused on a single principal element (e.g. nickel-based alloys) to target alloys that include high atomic fractions (X >10%) of multiple elements. These HEA materials can exhibit sluggish diffusion and enhanced corrosion resistance, ideal for potential applications in advanced ultra supercritical (A-USC) steam cycles for power generation. Scope: The addition of multiple elements in high atomic fractions creates an enormous design space that cannot easily be investigated by traditional material design strategies such as designed of experiments (DOE). This project utilizes a Bayesian machine learning algorithm that has been modified to work with calculation of phase diagrams (CALPHAD) software. This Bayesian algorithm reduces manual inputs and increase the likelihood of achieving an optimal solution. Compositional inputs to this algorithm will be assessed using existing material property models for high temperature strength and corrosion resistance. The target for alloy performance will be a 15% (~100 ⁰C) increase in allowable service temperature beyond heat-resistant stainless steels while maintaining or improving alloy cost and corrosion resistance. Haynes 230 was selected as a baseline, which is 57 wt% Ni with 22 wt% Cr 14 wt% W, and 2 wt% Mo as solid solution strengtheners. In addition to rapid design via Bayesian machine learning, the alloys were rapidly fabricated using a multi-wire arc additive manufacturing (mWAAM) technique which allows for precise control of alloy composition and assessing of alloy design “windows” to study composition effects. Build speeds for wire-arc additive processes are among the highest for additive technologies enabling rapid and reliable sample fabrication when compared to conventional methods such as arc button melting. The mWAAM samples will be rapidly characterized via instrumented indentation for room temperature modulus and strength and for elevated temperature strength via hot hardness tests. After being screened with hardness testing, potential alloys will be further evaluated with conventional microscopy techniques including scanning electron microscopy (SEM) and transmission electron microscopy (TEM) to assess agreement with modeling results. The most promising compositions will also be evaluated by printing full sized tensile specimens for mechanical behavior tests at elevated temperatures. Results: Bayesian machine learning of a single performance function was initially used to optimize five performance metrics: 1) single phase stability, 2) yield strength, 3) creep resistance (low diffusion coefficient), 4) freezing range (weldability), and 5) material cost. The single performance function was suboptimal as assumptions had to be made about the results while formulating the optimization. A goal-oriented Bayesian optimization strategy (Hanaoka, 2021) was implemented with CALPHAD for use with the five metrics above. This multi-objective Bayesian optimization (MOBO) enabled the design of NiCrCoFe alloys with V and W additions. A base composition of NiCoCr was selected as Ni provides a stable FCC matrix, Cr aids corrosion/oxidation resistance, and Co is a solid-solutions strengthener that also improves creep by increasing the activation energy. Fe helps reduce diffusion coefficients and cost. Finally, V and W were selected for their reasonable solubility and high atomic misfit to aid in solid solution strengthening. Cracking of the mWAAM specimens was an early issue, and the Easton solidification cracking model (Easton et al., 2014a) was selected for addition to the MOBO function. High performing alloys fabricated by mWAAM included Ni 28 Cr 25 Co 26 Fe 15 V 8 and Ni 62 Cr 18 Co 1 Fe 3 W 15 . It was observed that even after adapting the mWAAM process for W, the W did not fully dissolve. To fully evaluate the Ni 62 Cr 18 Co 1 Fe 3 W 15 composition, a cored wire (80-20 NiCr sheath/powder core) was manufactured and printed via WAAM, and HIP’ing was utilized to homogenize and densify the printed alloy. The V and W alloys produced met metrics 1 (solid solution), 4 (solidification cracking), and 5 (cost). However, an unmodeled mechanism of thermal stress cracking was identified in the WAAM produced materials, perhaps exacerbated by the lack of grain boundary strengthening elements (B, C). Conclusions & Recommendations: A high-throughput (rapid) alloy design technique was applied to designing and manufacturing new high entropy alloys (HEAs) for extreme environments utilizing MOBO and mWAAM. The developed process was rapid and effective in addressing the mechanisms included in the model. The lack of grain boundary strengthening element additions (e.g., B, C) was a simplification that likely produced thermal stress cracking that turned into a large part of the investigation. Additions on the order of 0.005 wt% B and 0.05 wt% C likely would have minimized thermal stress grain boundary cracking. Overall, the high throughput design strategy is promising for rapid design of metrics-driven alloys for advanced ultra supercritical (A-USC) steam cycles for power generation. The MOBO and mWAAM process could be commercialized to accelerate metrics-driven alloy design. In addition, the cored-wire process utilized for scale-up is a promising high-volume process for WAAM alloy development and scale-up.

36 MATERIALS SCIENCE↗

Development and Demonstration of a Fuel-Efficient, Class 8 Tractor & Trailer Engine System (SuperTruck II)

Navistar presents the SuperTruck II (ST II) Final Report to the Unites States Department of Energy (US DOE), which covers the five Budget Periods (BPs) from 10-1-2016 through 6-30-2022. For ST II, Navistar built on the achievements of the SuperTruck I (ST I) Program as a catalyst to continue critical research, design and development, testing, and operations to reach the ambitious goals of the ST II project. This approach allowed Navistar to continue contributing to the essential needs of our nation for safe, efficient, and cost-effective delivery of goods and services, as we reduced negative environmental effects and improved operational productivity. This document contains information specified in DOE F 4600.2, Final Scientific/Technical Report DOE F 241.3, B. SCIENTIFIC/TECHNICAL REPORTS, explaining how we met and exceeded program requirements. Throughout this Final Report, Navistar extracted information from documents prepared during the project that represent our management, design and development, building, and testing efforts to meet and exceed SuperTruck II project goals. Navistar followed Plan requirements to achieve / exceed Project Objectives: a) >100% improvement in vehicle freight efficiency (FE) (on ton-MPG basis) relative to 2009 baseline with stretch goal of 140% improvement [actual: 170%); b) >55% engine brake thermal efficiency (BTE) demonstrated in operational engine at a 65-mph cruise point on a dynamometer – ≥31% increase from 2009 baseline [actual: 55.20% of combined BTE) ; and c) development and implementation of commercially cost effective technologies (in terms of a simple payback). Technology selection / development path focused on developing technologies applicable for production within 3-year approach, while ensuring technology readiness and cost of ownership for end users. The Program was organized into five budget periods: Requirements / Technology Assessment and Initial Hardware Testing; Technology Development and Concept Readiness Demonstration; Technology Finalization and Validation Tractor / Trailer Fabrication, Integration and Commissioning Demonstration; and Fuel Economy (FE) and Brake Thermal Efficiency (BTE) and Program Completion. Leadership was provided by DOE, with tasks performed by laboratories (Argonne National Laboratory, Lawrence Livermore National Laboratory); partners at Bosch, TPI, Dana, and J.B. Hunt; , and support from University of Michigan and Clemson University. Navistar lead this team with Principal Investigator / Contracting Officer; Project Manager (PM); Vehicle, Engine, and Aftertreatment Engineers; Finance Manager, Technical Program Leads, and Legal/IP; and other key personnel. Work also included personnel in risk management; funding / budget / finance. Work involved analysis, development, testing, and down selection of individual/system engine, aftertreatment, and vehicle technologies, with integration of selected technologies into a prototype vehicle for demonstration of fuel-efficiency gain. Work also included component/integrated system level development of truck and trailer aerodynamics, base engine efficiency, advanced aftertreatment, combustion efficiency, waste heat recovery, hybrid powertrain, reduced rolling resistance, weight reduction, idle reduction, and driver feedback. As ST II progressed, Navistar performed computer-based modeling / simulations of technologies focused on the primary operational areas: Engine, Aftertreatment, and Vehicle. During the ST II Program, the COVID Virus outbreak unexpectedly challenged by the effects of, which affected staffing, scheduling, design, supplies, availability of materials, production procedures, and testing. The DOE responded by extending the program by three quarters to ensure that project tasks were completed for this vital project. Focus continued on analyzing, developing, testing, and down selecting individual-/system-level engine and vehicle technologies for integration of the final selected technologies into a prototype vehicle that would demonstrate fuel-efficiency gains made possible through these technologies. This included component/integrated system-level development of truck and trailer aerodynamics, base engine efficiency, advanced aftertreatment, combustion efficiency, waste heat recovery, solar power, distributed and intelligent vehicle power, hybrid powertrain, reduced rolling resistance, weight reduction, idle reduction, and driver feedback. Throughout the program, function, reliability, and performance at all levels were ensured through testing. Proof of this approach was demonstrated in multiple, on-road demonstrations: Scenario A (Flatland) Fuel Economy, Scenario B (Hilly) Fuel Economy, and City Cycle Tests. Other benefits derived from ST II included new/improved products, publications, patents, and next-step capabilities related to electric/hydrogen vehicles and autonomous driving.

Zukouski, Russ↗

Timelike Compton Scattering from a longitudinally polarised target with CLAS12 at Jefferson Lab

Explorations into the internal dynamics of hadrons are constantly evolving, and the requirement for experimental results to verify theoretical models of hadron structure is paramount. A key area in this field is the study of Generalised Parton Distributions (GPDs), which are functions used to model the momenta of quarks and gluons within hadrons, and the methods to access GPDs experimentally. One such scattering process that allows access to these is Timelike Compton Scattering. TCS complements existing Deeply Virtual Compton Scattering experiments and allows investigation into the universality of GPDs through access to the real and imaginary parts of the parton helicity independent GPD Hq via beam spin asymmetries (BSA), and it provides novel access to the real and imaginary parts of the parton helicity dependent GPD ˜Hq through target polarisation asymmetries (TSA). This thesis work presents a comparative study with the first published BSA for TCS at the Thomas Jefferson National Accelerator Facility (JLab), alongside a first time extraction of a Target Spin Asymmetry with the Summer 2022 data taking run. JLab hosts the Continuous Electron Beam Accelerator Facility (CEBAF) which provides a 12 GeV electron beam to four experimental halls. Hall-B contains the CEBAF Large Acceptance Spectrometer, which took data across three run periods on a longitudinally polarised NH3 and ND3 fixed target from 2022-2023, to extract measurements of electron-proton scattering, from which a TCS signal could be extracted. The thesis discusses work done to understand and eliminate contributions from the non-/low-polarised nuclear background, testing pre-established cuts to eliminate pion background from a dilepton (e+e-) final state and modifying them as needed for the new experimental run, and attempts to hone in on a clean TCS signal from which to extract the two asymmetry observables. A comparison with existing BSA results was performed; however, the statistical errors are too large to draw a significant conclusion as to whether there is agreement across each bin. More data is needed for a multidimensionally binned extraction. A proof of principle was achieved in the TSA measurements, with two out of four kinematic bins showing preliminary agreement in shape with theoretical values. Again, the errors are significant due to the contributions from the nuclear background. To support these conclusions, a further study was done, which takes into account an estimate of the asymmetries with the full available dataset (this thesis is based only on data taken in the summer set; at the time of writing processing was still being conducted for the final two datasets), as well as an estimate including additional future experiment days that were awarded in July 2024. Additional work was done on a secondary project exploring the feasibility of measuring TCS at the upcoming Electron Ion Collider, supporting the design proposal for the detector for the first interaction region and giving a positive outlook for the future of these types of measurements beyond JLab.

Gates, Kayleigh [Univ. of Glasgow, Scotland (Unite↗

Pauling’s rules for oxide-based minerals: A re-examination based on quantum mechanical constraints and modern applications of bond-valence theory to Earth materials

Since their introduction in 1929, Pauling’s five rules have been used by scientists from many disciplines to rationalize and predict stable arrangements of atoms and coordination polyhedra in crystalline solids; amorphous materials such as silicate glasses and melts; nanomaterials, poorly crystalline solids; aqueous cation and anion complexes; and sorption complexes at mineral-aqueous solution interfaces. The predictive power of these simple yet powerful rules was challenged recently by George et al. (2020), who performed a statistical analysis of the performance of Pauling’s five rules for about 5000 oxide crystal structures. They concluded that only 13% of the oxides satisfy the last four rules simultaneously and that the second rule has the most exceptions. They also found that Pauling’s first rule is satisfied for only 66% of the coordination environments tested and concluded that no simple rule linking ionic radius to coordination environment will be predictive due to the variable quality of univalent radii. We address these concerns and discuss quantum mechanical calculations that complement Pauling’s rules, particularly his first (radius sum and radius ratio rule) and second (electrostatic valence rule) rules. We also present a more realistic view of the bonded radii of atoms, derived by determining the local minimum in the electron density distribution measured along trajectories between bonded atoms known as bond paths, i.e., the bond critical point (rc). Electron density at the bond critical point is a quantum mechanical observable that correlates well with Pauling bond strength. Moreover, a metal atom in a polyhedron has as many bonded radii as it has bonded interactions, resulting in metal and O atoms that may not be spherical. O atoms, for example, are not spherical in many oxide-based crystal structures. Instead, the electron density of a bonded oxygen is often highly distorted or polarized, with its bonded radius decreasing systematically from ~1.38 Å when bonded to highly electropositive atoms like sodium to 0.64 Å when bonded to highly electronegative atoms like nitrogen. Bonded radii determined for metal atoms match the Shannon (1976) radii for more electropositive atoms, but the match decreases systematically as the electronegativities of the M atoms increase. As a result, significant departures from the radius ratio rule in the analysis by George et al. (2020) is not surprising. We offer a modified, more fundamental version of Pauling’s first rule and demonstrate that the second rule has a one-to-one connection between the electron density accumulated between the bonded atoms at the bond critical point and the Pauling bond strength of the bonded interaction. Pauling’s second rule implicitly assumes that bond strength is invariant with bond length for a given pair of bonded atoms. Many studies have since shown that this is not the case, and Brown and Shannon (1973) developed an equation and a set of parameters to describe the relation between bond length and bond strength, now redefined as bond valence to avoid confusion with Pauling bond-strength. Brown (1980) used the valence-sum rule, together with the path rule and the valence-matching principle, as the three axioms of bond-valence theory (BVT), a powerful method for understanding many otherwise elusive aspects of crystals and also their participation in dynamic processes. We show how a priori bond-valence calculations can predict unstrained bond-lengths and how bond-valence mapping can locate low-Z atoms in a crystal structure (e.g., Li) or examine possible diffusion pathways for atoms through crystal structures. In addition, we briefly discuss Pauling’s third, fourth, and fifth rules, the first two of which concern the sharing of polyhedron elements (edges and faces) and the common instability associated with structures in which a polyhedron shares an edge or face with another polyhedron and contains high-valence cations. The olivine [α-(MgxFe1–x)2SiO4] crystal structure is used to illustrate the distortions from hexagonal close-packing of O atoms caused by metal-metal repulsion across shared polyhedron edges. We conclude by discussing several applications of BVT to Earth materials, including the use of BVT to: (1) locate H+ ions in crystal structures, including the location of protons in the crystal structures of nominally anhydrous minerals in Earth’s mantle; (2) determine how strongly bonded (usually anionic) structural units interact with weakly bonded (usually cationic) interstitial complexes in complex uranyl-oxide and uranyl-oxysalt minerals using the valence-matching principle; (3) calculate Lewis acid strengths of cations and Lewis base strengths of anions; (4) determine how (H2O) groups can function as bond-valence transformers by dividing one bond into two bonds of half the bond valence; (5) help characterize products of sorption reactions of aqueous cations (e.g., Co2+ and Pb2+) and oxyanions [e.g., selenate (Se6+O4)2- and selenite (Se4+O3)2-] at mineral-aqueous solution interfaces and the important role of protons in these reactions; and (6) help characterize the local coordination environments of highly charged cations (e.g., Zr4+, Ti4+, U4+, U5+, and U6+) in silicate glasses and melts.

Geochemistry & Geophysics↗

First simultaneous measurement of the γ-ray and neutron emission probabilities in inverse kinematics at a heavy-ion storage ring

The probabilities for γ-ray and particle emission as a function of the excitation energy of a decaying nucleus are valuable observables for constraining the ingredients of the models that describe the deexcitation of nuclei near the particle emission threshold. These models are essential in nuclear astrophysics and applications. In this paper, we have for the first time simultaneously measured the γ-ray and neutron emission probabilities of 208 Pb. The measurement was performed in inverse kinematics at the Experimental Storage Ring (ESR) of the GSI/FAIR facility, where a 208 Pb beam interacted through the 208 Pb (p,p') reaction with a hydrogen gas jet target. Instead of detecting the γ rays and neutrons emitted by 208 Pb, we detected the heavy beamlike residues produced after γ and neutron emission. These heavy residues were fully separated by a dipole magnet of the ESR and were detected with outstanding efficiencies. The comparison of the measured probabilities with model calculations has allowed us to test and select different descriptions of the γ-ray strength function and the nuclear level density available in the literature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Engineered Barrier System R&D and International Collaborations – LANL (FY20): Spent Fuel and Waste Disposition

The United States Department of Energy’s Spent Fuel and Waste Disposition (SFWD) program is investigating the design and safety function of generic nuclear geologic repositories in a variety of geologic settings (salt, argillite, and crystalline rock). Different configurations and loadings of spent nuclear fuel and waste within disposal canisters are also being investigated, some of which have the potential to generate repository temperatures higher than previously considered (i.e., temperatures >100ºC) by foreign and domestic concepts. This report expands on engineered barrier material stability in a high temperature crystalline rock repository through high temperature hydrothermal experiments. Experiments were designed to develop engineered barrier system (EBS) concepts in a hightemperature crystalline environment in 1) bentonite-Grimsel Granodiorite interactions, 2) bentonite-cement reactions, and 3) interaction between waste canister materials and bentonite. Experiment results are applied to understanding long-term repository performance in terms of radionuclide isolation. One hydrothermal experiment was completed in the rocking autoclaves at LANL in FY-20: IEBS-6 (Grimsel Granodiorite + Wyoming bentonite + cured ordinary Portland cement + Grimsel Granodiorite synthetic groundwater, 250ºC/150 bar, 8 weeks). Several other experiments were planned but were delayed due to the COVID-19 pause in laboratory work. In addition, some characterization of the reaction products of experiments conducted in FY-20 was prevented by COVID-19. Mineral phase chemistry was not measured via electron microprobe analyses for IEBS-6. Quantitative X-ray diffraction results from IEBS-6 and HBT-1 were not completed. The missing analyses and discussion of the results will be included in next year’s report. New characterization that was completed includes scanning electron microscopy of reaction products from IEBS-6, quantitative X-ray diffraction results are presented from IEBS-1 through IBES-5, preliminary scanning electron microscope images and chemical analyses for IEBS-6, X-ray diffraction of the clay fraction from all experiments, measurement/imaging of mineral growth on the surface of steel coupons. Major observations pertaining to bentonite stability in a Grimsel Granodiorite environment include the stability of Na-montmorillonite at 250°C, the formation of trace CSH phases, and the formation of bentonite colloids on experiment cooling. The addition of a cured chip of Portland cement to the bentonite-Grimsel system results in slightly higher pH values and the formation of diverse secondary mineral phases that were not observed in the previous experiments (e.g., analcime, garronite, CSH phases). The new characterization efforts related to the interaction of stainless-steel coupons and bentonite clay focused on thickness and mineralogy of phases that formed at the steel surface. In Wyoming bentonite + Grimsel Granodiorite systems, newly formed minerals at the bentonitesteel coupon interface included alteration of the outermost steel edge to Fe,Cr-oxide phases, followed by Fe-rich phyllosilicates (Fe-saponite, chlorite) and interbedded Fe,Cr,Ni-sulfide phases (pentlandite). Hydrothermal experiments were completed to assess uranium-carbonate complexation at conditions relevant to high-temperature disposal. Autoclave solubility experiments were conducted at 150 to 250ºC with a range of carbonate and uranium concentrations. The experiment results were characterized via situ UV-Visible spectroscopy and synchrotron-based in situ XAS techniques. Results show a significant decrease in the stability of uranyl-carbonate complexes at temperatures above 100ºC along with the precipitation of uranium oxides. Further, at T > 200ºC, results show that uranyl-hydroxyl complexes control solubility of uranium instead of the previously predicted uranyl-carbonate species. These results are significant for understanding the mobility of uranium in the EBS, which will likely contain carbonate-rich fluids. International research efforts focused on three main areas: 1) participation in international conferences, 2) building collaborations with foreign repository programs, and 3) the initiation of an experimental program to complement the full-scale HotBENT test at the Grimsel test site. This experiment included Wyoming bentonite + low carbon steel + Grimsel Granodiorite synthetic groundwater and was run at the planned maximum temperature of the HotBENT test (200°C). Complete characterization of reaction products was hindered by disruptions to laboratory work but will be reported in the next FY. The experimental results obtained in FY-20 continue to document the wide-ranging effects of bulk composition and pressure-temperature conditions in the mineralogical and geochemical evolution of a high-temperature repository environment. Concepts developed will be used to inform models of long-term material stability in a generic crystalline rock-hosted repository.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of an ab initio learned model of electron deposition range in deuterium-tritium plasmas through time-dependent density functional theory calculations and machine learning

Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Plutonium Oxidation State Distribution in the Presence of WIPP-Relevant Organics and Iron Corrosion Products

The oxidation state and solubility of plutonium (Pu) in high ionic strength synthetic WIPP (Waste Isolation Pilot Plant) brines as a function of pC H+ in the presence and absence of WIPP-relevant organic ligands (EDTA [Ethylenediaminetetraacetic acid], oxalate, citrate, acetate) and iron corrosion products (magnetite and metallic iron) at 𝑇 = 23 ± 2 ∘C was thoroughly studied by long-term batch solubility experiments (between approximately 800-1,100 days) from an undersaturation approach. The oxidation state of Pu in the WIPP environment has been a topic of interest since the initial Compliance Certification Application (CCA). This study aims to investigate the solubility of Pu under the expected WIPP conditions and to determine the oxidation state of the solid phase that will control the solubility. One of the most important results of this study is that the Pu oxidation state was analyzed both from the surface area of the corrosion products and in the precipitated solid. The analysis of the Pu oxidation state on the surface of the iron mineral from the ongoing undersaturated experiments is more relevant to the performance assessment of nuclear waste disposal than short term batch (plutonium-iron phase) experiments. The X-ray Absorption Near-Edge Spectroscopy (XANES) analysis showed that Pu oxidation state is different on the metal surface and in the precipitated solid. Pu(III) is the dominant oxidation state in the metallic iron (Fe 0 ) system in the presence and absence of organics. Pu(IV) is the dominant oxidation state in the magnetite system in the presence and absence of organics. Also, organics stabilize Pu(IV) in the magnetite system. Analysis of Pu in the precipitated solid by Extended X-ray Absorption Fine Structure (EXAFS) analysis showed that Pu formed an inner-sphere complex with iron with minor amounts of PuO 2 present. X-ray diffraction (XRD) results indicate that metallic iron and magnetite did not oxidize in three years in the alkaline and high ionic strength system. Under these conditions (8 < pC H+ < 10 at T = (22 ± 2) °C under nitrogen atmosphere, the solubility of Pu changes by up to three orders of magnitude (10 -5 and 10 -8 M). The spread in solubility is highest at pC H+ = 9. Pu(III) and Pu(IV) showed different solubility behavior in the synthetic WIPP brine. A summary of the data collected in this report will be submitted to Sandia National Laboratories as part of a parameter update report which will outline the changes to the OXSTAT parameter for the 2026 Compliance Recertification Application (CRA-2026). The experiments performed were done according to the U.S. Department of Energy (DOE) approved Test Plan entitled “Effects of Radiolysis, Organic Complexation, and Redox Conditions on the Speciation and Oxidation State Distribution of Pu(III/IV)” (LCO-ACP-25). All data reported were obtained under the Los Alamos National Laboratory-Carlsbad Office (LANL-CO) Quality Assurance Program, which is compliant with the DOE Carlsbad Field Office, Quality Assurance Program Document (CBFO/QAPD).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Effect of Strength of Soluble and Solid Acid Catalysts on Selectivity in Aldol Reactions

Aldol reactions can follow the classical condensation path and produce an α,β-unsaturated carbonyl compound and water, however in some cases, aldol chemistry may be steered towards a fission that gives an olefin and a carboxylic acid. This investigation examines how the strength of Brønsted acid sites—both in homogeneous and heterogeneous catalysis—controls these two pathways and their kinetics. The cross-aldol reaction between benzaldehyde and 3-pentanone served as test case. Batch reactions, conducted in toluene as solvent at a temperature of 140 °C under autogenous pressure, were analyzed by GC and in situ ATR-FTIR spectroscopy to determine product distributions and rate constants for condensation and fission pathways. A series of soluble acids including a family of sulfonic acids mostly favored the condensation pathway, with formic acid as the weakest in the series by its pKa being inactive for aldol chemistry. Significant amounts of fission products were rare except for the known selectivity of phosphoric acid, and the circumstances were difficult to delineate. For the sulfonic acid family, the logarithm of the first-order condensation rate constants scaled only roughly with pKa (water) values, whereas a good correlation was obtained with calculated deprotonation Gibbs energies in toluene. A series of H-forms of isomorphously substituted beta zeolites, HESiBEA with E = Al, Ga, Fe, or B, favored the fission pathway. Site density and site strength were characterized by calorimetric measurements of the heats of adsorption of isopropylamine, which decreased in the order Al > Ga, Fe > B. The logarithm of the site-normalized first-order fission rate constants scaled roughly with the heats of adsorption and correlated excellently with reported deprotonation energies. In conclusion, acid strength mainly affects activity and can be seen as a pre-requisite for either aldol condensation or fission chemistry, whereas additional, yet to be full clarified catalyst properties and reaction conditions are required to steer aldol chemistry towards fission selectivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗