Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “volume distribution”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site

The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2021

Spectral Induced Polarization-Biogeochemical Relationships for Remediation Amendment Monitoring

Geophysical tools such as electrical resistivity (ER) can indirectly monitor subsurface changes in response to remedial injections. These methods exhibit relatively low spatial resolution compared to sediment core characterization but are advantageous due to the ability to collect measurements non-intrusively over time across large volumes of the subsurface. Moreover, along with confirmatory groundwater or core sampling, geophysical methods can be used during active biogeochemical remedies to monitor short-term contaminant transformations and mobility, as well as part of an overall strategy for long-term monitoring of subsurface contaminated sites. The use of alternating current spectral induced polarization (SIP) provides significantly more information than conventional geophysical methods like direct current ER. SIP allows for monitoring changes in both solution and surface conductivity by separation of real and imaginary conductivity, respectively, as well as surface capacitance. In principle, SIP can measure indicators of remedy progression such as precipitation reactions that sequester contaminants, potentially providing a better indication of amendment delivery and reactivity as compared to conventional ER methods. However, multiple processes and material properties have overlapping (interacting) electrical responses across a range of frequencies. Hence, the purpose of this scoping study was to evaluate if SIP can be used to monitor (a) amendment delivery and (b) precipitation and reactivity of amendments under consideration at Hanford. SIP measurements were collected in fully saturated columns packed with sand and Hanford formation sediments containing (a) amendments that were highly conductive with significant capacitance (zero valent iron – ZVI, sulfur modified iron – SMI) and (b) amendments that exhibited low electrical conductivity with a small capacitance (calcite, apatite, bismuth). The sand was a quartz material with homogenous particle size that exhibited a relatively low surface conductivity. It was used as a control for comparison with the sediments from the Hanford Site, which have a greater surface conductivity because of their complex mineralogy and heterogeneous particle size distribution and may have complex interactions with amendments. The amendment mass fraction was varied to represent the different stages and subsurface locations associated with the delivery of a remedy. The primary objective was to identify the solution and solid surface changes associated with the delivery amendments and their secondary reactions within the subsurface. The figure below summarizes results for the amendments tested through this project. The SIP phase shift or imaginary conductivity change for the high conductivity amendments was more than 10 times that of the low conductivity amendments, highlighting the relative ease of detection of ZVI and SMI independent of the background signal from sand or Hanford sediments. The ZVI phase shift and imaginary conductivity changes occur primarily at high frequency (> 100 Hz) whereas SMI changes were at low frequency (0.01 to 10 Hz). SIP signals of SMI also increased over time and the maximum shifted to lower frequencies. The low conductivity amendments (calcite, abiotic and biotic apatite, bismuth subnitrate) exhibited relatively small phase and imaginary conductivity changes when added to sediments. The changes were above minimum detection limits (0.5 mrad for phase shift, 0.03 µS/cm for imaginary conductivity) for the highest concentration except for the commercial bismuth material. The lowest amendment concentration that can be detected is likely sediment specific, as minerals in sediments (clays, magnetite, Fe-oxides) have some capacitance and, therefore, exhibit a variable phase shift.

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

Opportunities to Implement Solutions to Achieve Remedial Action Objectives in Consideration of Stakeholder Interests - 20365

A case study in the implementation of a soils excavation and removal project conducted at a Formerly Utilized Sites Remedial Action Program (FUSRAP) Maywood Superfund Site (FMSS) vicinity property; an active commercial business, in a densely populated area, with significant operational, technical and, logistical constraints that were expected to limit areas that could be remediated. The United States Army Corps of Engineers-led (USACE) team successfully navigated complex overlapping stakeholder interests to achieve a more efficient and complete removal of contaminated soils and debris while minimizing unnecessary excess costs to the Government. This paper centers on FUSRAP activities at the FMSS vicinity property located at 149-151 Maywood Avenue, Maywood, Bergen County, New Jersey. For most of its FUSRAP history, this ∼109,000 square meter (27-acre) property housed a now-demolished ∼26,000 square meter (6.5-acre) warehouse operated by Sears Logistics Services, a unit of the Sears Roebuck Company (Figure 2). Sears ended its property lease and vacated the warehouse in December 2016. Given that long history and for the purposes of this paper, the property will be referred to as the 'Sears property' or simply 'the property.' The Sears property is approximately 109,000 square meters (27-acres) and is currently zoned for commercial use. The property is bound to the north and northwest by 100 West Hunter Avenue (the Stepan Company), to the northeast by 205 Maywood Avenue, to the east by Maywood Avenue, to the south by 23 West Howcroft Road, and to the west by businesses on NJ Route 17 and the NJ Route 17 roadway. Until December 2016, the property housed a warehouse and distribution center operated by Sears Logistical Services. On-site structures were demolished by the property owner in 2017. The warehouse covered the north section of the property along the Stepan Company property line. A railroad spur from the adjacent property now known as the Maywood Interim Storage Site (MISS) ended at the northeast corner of the warehouse. Wetlands were located east of the warehouse; the rest of the property was covered with paved parking lots and grassed areas. During remediation at the property, opportunities frequently arose where the project team was able to coordinate effectively with stakeholders to sequence remediation activities to facilitate removal of otherwise inaccessible soils and identify opportunities for materials reuse when supported by residual radiological and chemical levels. Some specific opportunities include: - remediating individual truck loading dock bays to maintain tenant operations; - working along a busy highway and in a utility corridor containing a 30-inch high pressure gas main; - remediating Lodi Brook and associated wetlands requiring bypass pumping and other diversionary structures to maintain local community stormwater drainage; - protecting during remediation and supporting access to facilitate tenant/owner maintenance of critical fire protection and water supply system service lines buried in contaminated soils and necessary for safe warehouse operations; and - coordinating with stakeholders to ensure non-radiological contaminants of concern for the site (unaffiliated with FUSRAP) were addressed by the responsible party in a manner that maximized benefit to all parties while supporting an efficient overall site remediation program. Once the property was vacant and the warehouse structure and radiologically non-impacted above-grade structures were demolished, the project identified an opportunity to significantly reduce the volume of waste associated with the foundation of the former warehouse, a foundation suspected of being partially constructed in radiologically contaminated soils. The project developed and implemented supplemental radiological verification survey and sampling strategies based on radiological cross-contamination risk potential with process flow-charts and screening-level based decision points; a program that classified saw-cut sections of concrete based on observed residual surface radioactivity conditions using a combination of gamma sensitive sodium-iodide scans and beta sensitive Geiger-Muller direct measurements. The proposed approach was reviewed with stakeholders with feedback, including consideration of applicable State of New Jersey Site Remediation Program criteria and guidance before implementation. Once classified, additional sampling was performed at frequencies driven by screening. The sampling methods, stockpile sampling frequencies, criteria and, approach to data evaluation were developed in consideration of existing site cleanup standards, regional background ranges, State of New Jersey radiological remediation program guidance, and specific stakeholder input With the property remediation and survey efforts nearly completed, it is relevant to examine retrospectively the objectives, assumptions and limitations in the remedial design (i.e., what was planned) versus what was able to be accomplished through effective teaming. Benefits to other environmental remediation site programs include better understanding of how to work effectively with stakeholders to achieve win-win outcomes, and approaches to site remediation, waste minimization and materials reuse that were successful in their overall outcomes. (authors)

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

WHONDRS River Corridor Sediment and Water Geochemistry and In Situ Sensor Data from Machine-Learning-Informed Sites across the Contiguous United States (v6)

This dataset supports a broader study examining hyporheic zone respiration rates to improve predictive models at a contiguous United States (CONUS) scale. The CONUS-Scale Model-Sample Study (CM) was designed following ICON (integrated, coordinated, open, and networked) principles to facilitate a model-experiment (ModEx) iteration approach, leveraging crowdsourced sampling across the CONUS. New machine learning models were created every month to guide sampling locations. Data from the resulting samples were used to test and rebuild the machine learning models for the next round of sampling guidance. Sampling began in April 2022 and ended in October 2023. In addition to the widely distributed CONUS sites, a more spatially focused sampling occurred in the Yakima River Basin, WA in summer 2022. Data from this more spatially intensive sampling occurred under the label “Second Spatial Study (SSS)” and were also included in the machine learning models. Other data types collected from SSS that were not part of CM were published in a separate data package (https://data.ess-dive.lbl.gov/view/doi:10.15485/1969566). This data package was originally published in February 2023. It was updated in June 2023 (v2; new and modified files); December 2023 (v3; new and modified files); June 2024 (v4; new and modified files); April 2024 (v5; new and modified files); and September 2025 (v6; modified files). See the change history section in the readme for more details. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. This dataset is comprised of two folders of field photos and videos, one folder of raw Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) data and one main data folder containing (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) readme; (5) international generic sample number (IGSN) mapping file; (6) field protocols; (7) a subfolder with sample data; and (8) a subfolder with sensor data. The sample data subfolder contains (1) surface water and sediment dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC) data and averages; (2) surface water and sediment total nitrogen data and averages; (3) surface water major cations and anions and averages; (4) sediment grain size data; (5) sediment iron (II) data and averages; (6) wet sediment mass, dry sediment mass, water mass, and wet sediment volume in incubation and sediment ICR vials; (7) sediment incubation respiration rate data and averages; (8) normalized respiration rate data and averages; (9) methods codes; (10) sediment specific surface area; (11) sediment percent carbon and nitrogen; (12) sediment gravimetric moisture and averages; (15) sediment X-ray diffraction (XRD) data; (16) sediment adenosine triphosphate (ATP) and averages; (17) a subfolder with sediment incubation respiration data, scripts, and plots; (18) surface water and sediment FTICR methods; and (19) a subfolder of 9.4 Tesla (9.4T) FTICR-MS data. This folder contains five subfolders, one containing the sediment .xml data files, one containing the water .xml files, one containing the sediment CoreMS output files, one containing the water CoreMS output files, and the other containing instructions and scripts for processing the files in CoreMS (https://github.com/EMSL-Computing/CoreMS).The sensor data subfolder contains (1) a subfolder with miniDOT dissolved oxygen and temperature data and plots; (2) miniDOT dissolved oxygen and temperature summary data; and (3) miniDOT installation methods. All files are .csv, .pdf, .R, .xml, .d, .html, .Rmd, .py, .cal, .json, .jpg, .jpeg, .png, .mov, or .mp4. CORRECTION: Carbon and nitrogen content are reported as percentages. The current column headers "01395_C_percent_per_mg" and "01397_N_percent_per_mg" are incorrect. These should read "01395_C_percent" and "01397_N_percent" and will be corrected in the next version of this data package. We thank the United States Forest Service, Washington Department of Fish and Wildlife, Washington Department of Natural Resources, Cowiche Canyon Conservatory, Washington State Parks and Recreation Commission (Scientific Research Permit #210901), and the Confederated Tribes and Bands of the Yakama Nation for access to field locations where the samples labeled “SSS” were collected. We also thank the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview. WHONDRS consortium members were asked to provide any acknowledgments for the collection of samples labeled “CM” and the following is a list of acknowledgments that were submitted with their corresponding Site IDs: (MART) Research activities were conducted in part on the Wind River Experimental Forest within the Gifford Pinchot National Forest; (MP- 100379) Philadelphia is part of Lenapehoking, the ancestral homelands of the Lenape peoples; (MP-102398) Land surveyed is the ancestral homelands of the Nookhose'iinenno (Arapaho), Tsis tsis'tas (Cheyenne), and Nuuchu (Ute); (MP-100749 and MP- 100747) Georgia Coastal Ecosystem LTER, OCE-1832178; (SP-70 and SP-72) Eastern Shoshone, Shoshone-Bannock; (MP- 102944) Funded by Oregon Watershed Enhancement Board. On the traditional lands of the Confederated Tribes of the Siletz, Confederated Tribes of the Grand Rhonde, and the Clatsop-Nehalem Confederated Tribe; (MP- 100607) Holiday Creek is located on the traditional territory of the Monacan Indian Nation; (SP-45) Lafayette Blue Springs State Park; (MP-102420) NSF DEB-2016749; (MP-100019) New Hampshire Agriculture Experiment Station; (SP-35) Rayonier (land owner; https://www.rayonier.com/); (MP- 101276) US Department of Energy, Office of Science, Biological and Environmental Research, Subsurface Biogeochemical Research, Watershed Dynamics and Evolution SFA at ORNL; (MP- 103224) Watershed Dynamics and Evolution SFA at ORNL; (MP- 101584) Traditional lands of the Oceti Sakowin (Dakota, Lakota, Nakoda) and Anishinaabe Peoples.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Geomechanical in situ testing of fault reactivation in argillite repositories

Abstract. Pressurization of natural faults as a result of repository-induced effects can lead to their reactivation and permeability generation in case such features are present near disposal tunnels. Potential driving forces for such pressurization are the temperature increase caused by heat-producing high-level radioactive waste and the generation of hydrogen and other gases due to corrosion of engineered materials. We are primarily concerned about pressurization and fault activation in host rocks and faults that have very low natural permeability for pore pressure increases to dissipate, such as the argillite rocks currently investigated in Switzerland, France and other countries. This presentation discusses a series of in situ experiments of fault activation by fluid injection conducted in the argillite rock (Opalinus clay) at the Mont Terri underground research laboratory in Switzerland. A multi-model monitoring test bed was installed at Mont Terri that includes distributed fiber optics for strain, temperature, and acoustics; local fault pore pressure and three-dimensional displacement sensors; active seismic imaging; and passive seismic monitoring. The fault experiments (and their subsequent analysis via hydromechanical modeling) provide a new fundamental understanding of the coupling between pore pressure, fault deformation and permeability generation as a function of time and allow exploring how seismic and aseismic events may impact the integrity of faulted argillite host rock. Our experimental observations demonstrate that significant flow (and transport) can occur along the initially impermeable argillite fault when rupture is activated. However, the fault permeability decreases to almost its pre-activation value when fluid injection ceases and fluid pressure drops. Dilatant slip on the fault plane alone does not explain the observed increase in fault permeability; pressure-induced fault opening also plays a role, favored by the softness of the shale along with the fact that the structure of the fault zone prevents fluids from diffusing into the adjacent damage zone. Rupture initiation and permeability generation is initially aseismic, associated only with an increase in noise level and emerging tremors. Micro-earthquakes are initiated later in the experiments and typically occur away from the fluid-pressurized area. Hydromechanical models show that stress transferred from the initial aseismic deformation can build up to stress criticality and later induce seismic rupture. After presenting the experimental results, we will close the presentation with an outlook to future experimental campaigns using the fault test bed at Mont Terri. We are currently planning a controlled thermal stimulation of the fault, where instead of fluid injection as a trigger mechanism we will heat up the nearby rock volume and measure potential effects on fault stability. Lessons learned from our past and future experiments will help inform the safety assessment of geologic disposal in argillite host rock.

Birkholzer, Jens T.↗
DOE OSTI2022

Application of a Machine Learning Algorithm in Generating an Evapotranspiration Data Product From Coupled Thermal Infrared and Microwave Satellite Observations

Land surface evapotranspiration (ET) is one of the main energy sources for atmospheric dynamics and a critical component of the local, regional, and global water cycles. Consequently, accurate measurement or estimation of ET is one of the most active topics in hydro-climatology research. With massive and spatially distributed observational data sets of land surface properties and environmental conditions being collected from the ground, airborne or space-borne platforms daily over the past few decades, many research teams have started to use big data science to advance the ET estimation methods. The Geostationary satellite Evapotranspiration and Drought (GET-D) product system was developed at the National Oceanic and Atmospheric Administration (NOAA) in 2016 to generate daily ET and drought maps operationally. The primary inputs of the current GET-D system are the thermal infrared (TIR) observations from NOAA GOES satellite series. Because of the cloud contamination to the TIR observations, the spatial coverage of the daily GET-D ET product has been severely impacted. Based on the most recent advances, we have tested a machine learning algorithm to estimate all-weather land surface temperature (LST) from TIR and microwave (MW) combined satellite observations. With the regression tree machine learning approach, we can combine the high accuracy and high spatial resolution of GOES TIR data with the better spatial coverage of passive microwave observations and LST simulations from a land surface model (LSM). The regression tree model combines the three LST data sources for both clear and cloudy days, which enables the GET-D system to derive an all-weather ET product. This paper reports how the all-weather LST and ET are generated in the upgraded GET-D system and provides an evaluation of these LST and ET estimates with ground measurements. The results demonstrate that the regression tree machine learning method is feasible and effective for generating daily ET under all weather conditions with satisfactory accuracy from the big volume of satellite observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

SolarAPP+ Performance Review (2022 Data)

The Solar Automated Permit Processing Plus (SolarAPP+) platform is an online portal to facilitate and expedite rooftop solar photovoltaic (PV) permitting processes. SolarAPP+ allows PV contractors to upload system specifications, have those specifications automatically reviewed for code compliance, and receive instant approval for code-compliant systems. SolarAPP+ also provides inspection checklists to verify installation practices and adherence to approved designs. SolarAPP+ is available to authorities having jurisdiction (AHJs) at no cost. This report is part of an ongoing series of reviews of SolarAPP+ performance. Consistent with previous performance reviews, we summarize SolarAPP+ adoption trends to date and compare various metrics for PV systems permitted through SolarAPP+ versus systems permitted through conventional AHJ permitting processes. As of the end of 2022, the National Renewable Energy Laboratory (NREL) had contacted over 1,500 AHJs with significant solar permitting volume regarding SolarAPP+. Of those, 607 AHJs had at least expressed interest in the platform. 16 AHJs had begun piloting the platform and 15 of these had publicly launched the platform by the end of 2022. In 2022, 206 installers submitted 11,092 permits through the SolarAPP+ platform, including 708 permits for solar+storage systems. SolarAPP+ permits accounted for around 37% of all permits issued in participating AHJs. We compare permitting timelines through SolarAPP+ to traditional AHJ permitting processes to assess the platform's performance. Consistent with previous SolarAPP+ performance reviews, we find that permitting timelines are significantly shorter for SolarAPP+ projects. Based on median timelines, a typical SolarAPP+ project is permitted and inspected 8 business days sooner than traditional projects. We estimate that automatic SolarAPP+ permitting saved between 3,500 and 13,900 hours of AHJ staff time in 2022. Finally, we find evidence that SolarAPP+ may improve inspection outcomes, with SolarAPP+ projects failing inspections about 28% less frequently than traditional projects.

14 SOLAR ENERGY↗
DOE OSTI2020

Comparison of active interrogation methods for source location in a scattering and absorbing medium, consisting of PGNAA, and an AmBe quasi-forward biased directional source

Source location of Special Nuclear Material (SNM) encompassing 95% 235 U and 239 Pu is identified by utilizing a directional source from Patent No: US20190013109A1 using Prompt Gamma Neutron Activation Analysis (PGNAA) and neutron spectroscopy simulated with Monte-Carlo N-Particle transport 6.2 (MCNP). BC-408, HPGe, LaBr 3 detector arrays were used to identify the location of the SNM using total counts incident on each detector, and PGNAA photopeaks from HPGe and LaBr 3 detector arrays in a polyethylene shield. The conducted simulations varied the volume and location of the SNM in the MCNP input files to observe how the source location method behaved. PGNAA photopeaks used for source identification include 61 keV from fission, 2.223 MeV prompt gamma from hydrogen, 511 keV annihilation, and a single and double escape peaks from the prompt gamma interaction from hydrogen. The capabilities of each detector systems to acquire well resolved photopeaks with a 1% relative error or less, and total relative error for F4 and F8 tallies were less than 0.015% relative error. Source predictions of the SNM with uneven amounts of polyethylene shielding between the source and detectors was observed to overpredict and give invalid source location predictions. Source locations of the SNM with even amounts of polyethylene material between the source and each detector were found to be valid. With a 1 Ci 241 Am source activity, it was determined that 1630 s were needed to obtain the results for each detector system with the quasi-forward directional AmBe source. Coupling source and material identification together would increase acquisition time but would only require one system to determine.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
DOE OSTI2025

High-Speed and High-Quality Field Welding Repair Based on Advanced Non-Destructive Evaluation and Numerical Modeling

Creep strength-enhanced ferritic (CSEF) steels such as Grade 91 (9Cr-1Mo-V) and Grade 92 (Fe-9Cr-2W-0.5Mo) steels are widely used in the fossil-fuel-fired and nuclear power plants. The weld integrity of these steels is crucial for power plants' safe and reliable operations. Due to harsh service conditions, the steel weld can become susceptible to environmental degradation. Field welding repair is used to restore the degraded weld’s performance where a controlled temper-bead welding technique is commonly used to temper the freshly formed martensite during welding. However, knowledge of weld repairability is limited and experimental trial and error optimization to achieve desired microstructure and joint properties is expensive and time-consuming. Many existing computational models, e.g., finite element models, are limited to solving heat conduction equation and ignoring convective heat transfer due to molten metal flow. These models can result in over-prediction of peak temperatures of weld pool and heat-affected zone (HAZ), which in turn can affect the accuracy of tempering prediction. Moreover, these finite element models require an input of the deposit profiles in advance and thus limits the usability of these models. Here, a molten pool-based, multi-pass multi-layer model has been developed based on computational fluid dynamics (CFD) approach with the Volume of Fluid (VOF) method. The model calculates the bead formation, thereby eliminating the need for pre-determined bead profiles required by finite element models. For computational efficiency, a coordinate system attached to the moving heat source is utilized. A subroutine is developed to convert the temperature profiles in the reference frame stationary to the heat source to that stationary to the workpiece. The converted thermal cycles are then imported into a microstructure model to compute the tempering kinetics and resultant hardness using a Johnson-Mehl-Avrami-Kolmogorov (JMAK), and modified Grange-Baughman parameter. The modeling approach is first developed and validated on single- and multi-pass deposition of stainless steel filler metal onto a SA-533 high strength steel substrate. The models are then applied to a multi-pass V-groove repair weld of Grade 91 steel plate as well as directed energy deposition of Grade 92 steel. Non-destructive characterization of microstructures was performed on Grade 91 and 92 steel welds. Two welding processes, cold metal transfer (CMT) and flux-cored arc welding (FCAW), were investigated for the Grade 91 steel weld samples. For the Grade 92 weld samples, three different heat inputs (low, medium, and high) of gas tungsten arc welding (GTAW) were utilized to replicate traditional field welding processes. The non-destructive evaluation (NDE) method used for this research was immersion ultrasonic testing (UT) using a micro-resolution ultrasonic imaging methodology specifically designed to operate in the through-transmission configuration operating at 20 MHz of frequency. The system used a focused ultrasonic beam spot size diameter between 250-300 μm, and a 6 μm laser vibrometer spot size for detection, to produce highly defined images with longitudinal and mode-converted shear waves. From the micro-resolution ultrasonic C-scan images, three microstructural regions, i.e., weld metal (WM), HAZ, and base metal (BM), were clearly identifiable. Various levels of ultrasonic amplitudes distributed over the three regions were correlated with electron beam backscattered diffraction (EBSD) images using grain size, grain boundaries, and dislocation densities. The results showed that areas with relatively higher ultrasonic amplitude levels were associated with smaller grains and higher dislocation densities, while areas with lower amplitude levels were associated with larger grains and lower dislocation densities. In addition, ultrasonic velocity data obtained across the three different weld microstructural regions of Grade 91 test samples were correlated with optical metallographic images and hardness measurements. The results showed distinctive decreases in ultrasonic velocity and hardness over the HAZ region, where weld failures often occur during service.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Marine Hydrogen Demonstration

This report summarizes Phase 1 of a project involving the design of a Floating Hydrogen Production and Dispensing Barge destined for the Port of San Francisco (SF). The H 2 Barge is designed to produce renewable H 2 at the rate of ~ 530 kg/day, storing 512 kg of hydrogen at 517-bar, allowing fast refueling of hydrogen fuel cell vessels and land-side hydrogen delivery trailers for distribution into the nascent SF hydrogen ecosystem. The broader considerations that impacted the H2 Barge design are also described. An account is given of a new review process formulated by the United States Coast Guard (USCG) to review this first-of-its-kind maritime implementation of hydrogen technology. The immediate goals of the H 2 Barge Project are to 1) demonstrate the feasibility, viability and methods of hydrogen production, storage and fueling in a maritime context, 2) help shape (where needed) and navigate the required local, state and federal regulatory gauntlet and 3) catalyze a “green hydrogen ecosystem” (both marine and landside) with locally produced renewable hydrogen at the San Francisco waterfront. A summary is also given of the modeling and experimental activity of Phase 1 directed to the development of science-based refueling protocols for large marine Type IV 250-bar hydrogen tanks. Combined modeling and experimental studies are reported of the filling of large (28 kg) 250-bar Type IV hydrogen tanks of the type being deployed on early hydrogen ferries, such as the MV Sea Change. The primary question was to determine how such tanks can be successfully filled (state of charge greater than 97%) within 45 minutes without exceeding the 82 °C temperature limit historically set for such tanks. The studies show that a gas injector is needed avoid thermal stratification during filling which can result in potential hot spots. Pre-cooling of the hydrogen was found to be essential in most cases, as ambient conditions greatly affect the need for a pre-cooling to achieve the 45-minute fill time. Pre-cooling cannot be supplied by nearby water, such as that found in nature (bays, lakes, rivers, etc.) because pre-cooling cooling below 0 °C was found to be necessary to avoid excessive compression heating. The experimental results afforded a calibration of the engineering model (SOFIL) for these large 250-bar tanks, which now enables using SOFIL to predict volume-averaged hydrogen filling temperatures to an accuracy of +/- 2.7°C for these tanks. The model can therefore be used to evaluate potential scenarios for development of a standardized fueling methodology for ferries utilizing large Type-IV tanks of the type examined here.

08 HYDROGEN↗
DOE OSTI2021

Understanding the Impacts of Stellar Companions on Planet Formation and Evolution: A Survey of Stellar and Planetary Companions within 25 pc

We explore the impact of outer stellar companions on the occurrence rate of giant planets detected with radial velocities. We searched for stellar and planetary companions to a volume-limited sample of solar-type stars within 25 pc. Using adaptive optics imaging observations from the Lick 3 m and Palomar 200″ Telescopes, we characterized the multiplicity of our sample stars, down to the bottom of the main sequence. With these data, we confirm field star multiplicity statistics from previous surveys. We additionally combined three decades of radial velocity (RV) data from the California Planet Search with newly collected RV data from Keck/HIRES and the Automated Planet Finder/Levy Spectrometer to search for planetary companions in these same systems. Using an updated catalog of both stellar and planetary companions, as well as detailed injection/recovery tests to determine our sensitivity and completeness, we measured the occurrence rate of planets among the single- and multiple-star systems. We found that planets with masses in the range of 0.1–10 M {sub J} and with semimajor axes of 0.1–10 au have an occurrence rate of 0.18{sub −0.03}{sup +0.04} planets per star when they orbit single stars and an occurrence rate of 0.12 ± 0.04 planets per star when they orbit a star in a binary system. Breaking the sample down by the binary separation, we found that only one planet-hosting binary system had a binary separation <100 au, and none had a separation <50 au. These numbers yielded planet occurrence rates of 0.20{sub −0.06}{sup +0.07} planets per star for binaries with separation a {sub B} > 100 au and 0.04{sub −0.02}{sup +0.04} planets per star for binaries with separation a {sub B} < 100 au. The similarity in the planet occurrence rate around single stars and wide primaries implies that wide binary systems should actually host more planets than single-star systems, since they have more potential host stars. We estimated a system-wide planet occurrence rate of 0.3 planets per wide binary system for binaries with separations a {sub B} > 100 au. Finally, we found evidence that giant planets in binary systems have a different semimajor-axis distribution than their counterparts in single-star systems. The planets in the single-star sample had a significantly higher occurrence rate outside of 1 au than inside 1 au by nearly 4σ, in line with expectations that giant planets are most common near the snow line. However, the planets in the wide binary systems did not follow this distribution, but rather had equivalent occurrence rates interior and exterior to 1 au. This may point to binary-mediated planet migration acting on our sample, even in binaries wider than 100 au.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Controlling Factors Analyses of Sediment Dispersal From Source-to-Sink Systems of Lacustrine Rift Basin: A Case Study From Paleogene Shahejie-Dongying Formations of Nanpu Sag, Eastern China

Tectonic activity not only shapes the basic stratigraphic framework of rift basins, but also profoundly affects the sediment dispersal in rift basins. In this study, analyses of heavy mineral assemblages in different periods demonstrate that there are three obvious tectono-sedimentary evolutionary stages (E s 3 –E s 2 , E s 1 –E d 2 , and E d 1 , respectively) in the Paleogene provenance area of Nanpu Sag, and the volume of sand bodies increases from the bottom of the Paleogene Shahejie (E s ) Formation to the top of the Dongying (E d ) Formation in Nanpu Sag. Besides, this study comprehensively utilize the analyses of seismic interpretation, palynology, heavy mineral, and borehole core samples to investigate the controlling factors of sediment dispersal in the rift basin. The assemblages of heavy minerals in different periods reflect the rock composition and catchment area of different provenance areas, and their vertical differences reflect the evolution process of the provenance area and reflect the uplift-denudation process of the provenance area. The results reveal that the synergy of the evolution of tectonic activity and the adjustment of topographical evolution are the main controlling factors of sediment dispersal in Nanpu Sag, while climate change is not the main controlling factor. We conclude that an increased sediment supply rate in the long term reflects the control of tectonic activity on basin topography, rather than climate fluctuations. The differences in morphological modification result in differential sediment dispersal, which is principally related to the differential extrusion of the fault system. The catchment area and provenance distance adjustment is evidenced by the vertical changes of heavy mineral characteristics of single-well and interaction and linkage of boundary faults, and the adjustment of topography evolution. A consideration is that the interaction and linkage of boundary faults and complex subsidence history are multi-directional, and differential evolution of provenance area is universal in lacustrine rift basins, all of this highlights the adjustment of sediment pathways generated by this characteristic of rift basins and emphasizes the importance of controlling factors analyses in understanding differential sediment dispersal that presents in the rift basins. Besides, four sets of sediment dispersal patterns were delineated based on different developmental regions in the rift basin, which are fault segmental point and multi-stage fault terrace, single-stage fault terrace and axial fault valley, axial fault terrace, and paleo-terrace and axial fault valley, respectively. This study has a certain guiding significance for the prediction of the spatial distribution of sand bodies in the rift basin and the exploration of potential oil and gas targets in the rift basin.

Yu, Jifeng↗
DOE OSTI2020

InSAR Monitoring to Evaluate Surface Changes with CO 2 Storage in a Depleted Oil Field in Northern Michigan (Volume II.E)

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. This report describes the effort under MRCP to evaluate the potential use of Interferometric Synthetic Aperture Radar (InSAR) for monitoring the outcome of injected CO 2 at the Dover 33 reef near Gaylord, Michigan. InSAR is a satellite-based technology that provides high-precision information on the movement of ground surface in areas with high radar coherence (e.g., roads, buildings, bare soils). Depending on the setting, this technique may provide a useful tool for characterizing reservoirs by measuring surface deformations from activities such as brine water disposal; production of water, oil, and/or gas; and carbon capture, utilization, and storage (CCUS).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Hexagonal Distributed Embedded Energy Converters (HexDEECs)

Distributed Embedded Energy Converter Technologies (DEEC-Tec) is a new domain for marine renewable energy research that utilizes a conglomeration of small distributed embedded energy converters (DEECs) that, in aggregate, form larger metamaterial frameworks. These resulting DEEC-Tec metamaterials can then, in turn, be used to construct flexible ocean wave energy converters called flexWECs. DEEC-Tec enables flexWECs: (i) to be inherently broad-banded ocean wave frequency energy converters and (ii) to have an inherent lack of highly loaded rigid bodies. The DEEC-Tec domain also benefits the marine renewable energy domain by inherently availing ways that marine energy can be harvested and converted that heretofore has not yet been considered possible: real-time execution of transforming topologies (e.g., actively changing a flexWEC's shape and form) and morphologies (e.g., actively changing a flexWEC's stiffness and damping throughout its entire structure). Presented, is one specific type of DEEC, a HexDEEC, that shows promise in aiding the adoption and further development of the DEEC-Tec domain - it is a small energy transducer being developed by the United States National Renewable Energy Laboratory. The HexDEEC is a small (characteristic length approximating a centimeter) energy transducer that converts the dynamic deformations of an elastomer into electricity through a charging-discharging cycle of a capacitor whose capacitance is varied by those elastic deformations. The HexDEEC is composed of a hyperelastic hexagonal housing (nominally silicon rubber) with six electrodes on its inner faces. The upper three electrodes share the same charge while the lower three electrodes oppose the upper electrode charges. Externally, the HexDEEC has two arms extending away from the middle vertices of the hexagon. Via principles governing the relationship between electrical capacitance and electrical potential (voltage and charge), electricity is generated when the HexDEEC's arms are dynamically pulled or released under tensile loading as doing so causes the distance between the upper and lower sets of electrodes to change - varying the energy converter's overall capacitance. Analytical and numerical modeling have already been used to estimate the electrical energy produced by a HexDEEC. The cursory models approximate the HexDEEC as a parallel plate variable capacitor - simplifying from six to two opposing plates with a constant dielectric volume between those two plates. To account for the elastic HexDEEC material properties, software such as SolidWorks and STAR-CCM+ have been used to generate hyperelastic models; notably, Mooney-Rivlin based models. Individual HexDEECs have been fabricated by drawing uncured liquid silicon rubber into molds via vacuum pressure. To simplify manufacturing, HexDEEC sub-components - e.g., electrodes, wires - can be placed within those molds such that they are directly embedded into the hexagonal housing during the curing process. Furthermore, DEEC-Tec metamaterials can be created by interweaving or sequentially layering multiple HexDEEC strands together. The HexDEEC based metamaterial could then generate electricity through its gross deformations. Ultimately, HexDEECs represent a specific type of energy transducer that can be leveraged, by the DEEC-Tec domain, to create metamaterials used to construct novel flexWECs.

DEEC-Tec↗
DOE OSTI2022

Fabrication of Advanced Nanocarbon-Metal Composites for Improved Energy Efficiency

We fabricated nanocarbon metal composites (NCMC) of Al alloys with a process called “electrocharging assisted process” (EAP). This method consists of the application of a high current to a mixture of liquid metal and carbon particles. We investigated aluminum alloys (6061 and 1350) and used activated carbon and graphite powder as the source of carbon. The high current induces the formation of carbon chains and ribbons in the liquid metal. Upon solidification of the metal an epitaxial relation between the carbon nanostructures and the metal lattice is produced. The purpose of the project was to find the parameters during the reaction that would give rise to a high density of nanoribbons extending throughout the metal such that the electrical conductivity and the mechanical strength of the composite increased and that the method could be extended to large scale manufacturing. For this purpose, we designed two reactors for the incorporation of nanocarbon ribbons in aluminum metal. The first reactor was designed to have more control of the region with the high current density. However, there were problems with getting good mixing of the carbon in small volumes. A second reactor was designed with a stirrer that allowed for better mixing of the carbon by the introduction of argon gas through the shaft of the stirrer. NCMC were fabricated with a series of parameters to understand the role of current, time of applied current, type and shape of the cathode electrode and stirring speed of the mixture. We analyzed the composites by Raman scattering to gain information on the crystallite size of the nanocarbon, XRD to obtain the crystal structure of the composite, SEM and TEM to characterize the grain size of the aluminum grains and the quality of the crystal structure. We also measured the electrical conductivity and mechanical properties of selected samples. The nanocomposites showed increase in crystallite size of the nanocarbon with a linear dependence of the crystallite size on the duration of the applied current. There is a minimum current density necessary for the crystallite size to increase compared to the ~ 10 nm size of the activated carbon source used in the fabrication. The electrical conductivity of the nanocomposites increased with crystallite size of the nanocarbon and with the concentration of converted carbon. The maximum increase in electrical conductivity was 5.7% above the baseline for samples with ~ 4 wt % nanocarbon with crystallite size larger than 30 nm. These samples presented a hardness ~ 8% higher than the baseline samples with no carbon.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

PowerAmerica (Final Technical Report)

The U.S. Department of Energy’s Advanced Manufacturing Office (predecessor to AMMTO) established PowerAmerica in December 2014 to develop and accelerate the adoption of wide bandgap (WBG) semiconductor chips and power electronics in manufacturing. The objective was to spark early commercialization of energy efficient products; educate and train the workforce; create high-tech jobs; and nurture the growth of the U.S. WBG semiconductor manufacturing industry. PowerAmerica isled by North Carolina State University by way of a five-year, $140M cooperative agreement with DoE. This private-public partnership with DoE includes member companies ranging from startups and small-medium enterprises to large system integrators, world-class universities, and national labs. The WBG power electronic ecosystem formed by our diverse membership is focused on using advanced manufacturing to 1) lower the cost of silicon carbide and gallium nitride semiconductor devices to be comparable to silicon devices; 2) demonstrate the system benefits and energy efficiency advantages of WBG semiconductor power electronics through system demonstrations that validate their effectiveness; and 3) build an education pipeline for a skilled workforce to meet the future demand for emerging WBG semiconductor markets — and enhance U.S. economic competitiveness globally. PowerAmerica was initially funded in six budget periods, each lasting 12 to 18 months. By the end of Budget Period 6 (August 2023), PowerAmerica had achieved its major objectives in technology development, semiconductor device cost reduction, ecosystem growth and engagement, and education and workforce development. Through strategic foundry investments, we helped to establish the first U.S. SiC foundry (X-FAB) and supported a lab (Microchip) to start SiC volume production. We’ve helped bring together companies from different parts of the supply chain, resulting in several successful new partnerships and business relationships. Through projects with some of the largest manufacturers of energy-intensive equipment in the U.S. — John Deere, GE, United Technologies, Raytheon, Carrier, Toshiba, and others — we have successfully demonstrated the energy and system benefits of WBG technology. We have also harnessed the unique capabilities and facilities of national labs — the Naval Research Laboratory, National Renewable Energy Laboratory, and Argonne National Laboratory — to help solve challenging technical problems for industry. Thanks to our many webinars, tutorials, annual events, and presence at major energy and electronics conferences around the world, the PowerAmerica name has become synonymous with WBG technology advancement. The commercial interest in WBG technology is higher than ever; companies and governments around the world have announced hundreds of millions of dollarsin future investment to build capacity — and compete with silicon in many markets and applications. We have trained hundreds of engineering students, from universities across the U.S., through hands-on projects and WBG coursework. Working professionals have also benefited through the years from our many targeted short courses, tutorials, and technical webinars. In short, PowerAmerica has made great strides in each of our key objectives, and the organization has been operating without government or NC State

14 SOLAR ENERGY↗
DOE OSTI2024

Segmentation of RDX and TNT in X‐Ray Computed Tomography Reconstructions of Melt‐Cast Explosives

ABSTRACT Three‐dimensional mesoscale characterization of heterogeneous melt‐cast high explosives is challenging because of the difficulty differentiating binder from explosive crystals: two functionally different materials which are typically similar in density by design. Here, we report an algorithm which can differentiate hexahydro‐1,3,5‐trinitro‐1,3,5‐triazine (RDX) from 2,4,6‐trinitrotoluene (TNT) in x‐ray computed tomography (CT) volumes with tens of microns resolution. This method allows us to quantify RDX/TNT content, porosity, and RDX domain size. We calibrated the segmentation algorithm using simulated x‐ray CT volumes containing object models of RDX crystals within a TNT matrix. We then segmented and analyzed CT data for Composition B (Comp B), a 60/40 RDX/TNT mixture, and Cyclotol, a 75/25 RDX/TNT mixture. We examined melt‐cast samples fabricated with 100% theoretical maximum density (TMD) and 85% TMD. For the 100% TMD Comp B and Cyclotol samples, the RDX content values calculated by segmentation were 3% and 9% lower, respectively, than the values measured by high‐performance liquid chromatography on material from the same synthesis lots. This result is consistent with the expected underreporting of RDX content resulting from x‐ray CT resolution limits on RDX particles with diameters smaller than 25 µm. The 85% TMD samples were less accurately segmented with our algorithm due to the confounding presence of voids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗