Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scenario selection”

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 595 records · Page 33

Direct Air Capture of CO 2 and Delivery to Photobioreactors for Algal Biofuel Production (Final Report)

A mobile DAC system was designed and constructed to pair with photobioreactors growing algae for biofuel production. The DAC system was designed as a versatile research system, rather than a compact production unit. The system was constructed and mounted on a mobile skid to facilitate transportation to the algae production site. Within the DAC system, CO 2 was captured using amine-loaded monoliths that allow for high CO 2 uptake with low pressure drop. The CO 2 is collected using a Global Thermostat patented temperature/vacuum swing adsorption (TVSA) process. Amine sorbents and process conditions were optimized to produce 10 to >20 g CO 2 .h-1. The stability of the amine sorbents was also studied, with sorbent modifications made to improve stability to degradation by oxidation. An Algenol-developed Spirulina strain (Arthrospira platensis AB2293) was selected as the production cyanobacterial strain. AB2293 cultured was inoculum for outdoor production following PBR installation by Algenol. The PBR system was composed of three independent PBRs, with each PBR composed of four hanging bags internally recirculated by a liquid turnover pump. The PBRs were operated outdoors in Atlanta, GA, and integrated with the DAC system. Algae were grown with similar productivity using DAC-CO 2 as algae grown using pure CO 2 obtained commercially (Airgas). Throughout the experimental duration, no discoloration was observed, and cellular morphology was consistent between the two experimental treatments. An LCA including lifecycle greenhouse gas emissions, full life cycle inventory of the Algenol system and the DAC system and integrated DAC+PBR system was developed. Lifecycle greenhouse gas emissions were calculated for capture of carbon dioxide using input from Global Thermostat and the National Renewable Energy Laboratory. Three scenarios for energy provision were evaluated: a natural gas combined heat and power system sized to meet the electricity requirement, a natural gas combined heat and power system sized to meet the process heat requirements, and a system without on-site power that procures the electricity from the grid. In all three cases, as expected, the major contributor to the emissions is the energy consumption associated with the desorption step of the DAC process. The LCA quantified the reduced potential energy and greenhouse gas emissions of heat and mass integration of DAC and Algenol compared to unintegrated DAC and Algenol systems. A life cycle assessment of the role of sorbent productivity and lifetime was also developed. The development of more robust, oxidation resistant DAC sorbents may enable small reductions in energy requirements and in lifecycle greenhouse gas emissions and other environmental impacts. NREL performed techno-economic analysis (TEA) to identify the integration scenario most likely to achieve a 15% cost reduction target versus the baseline. Heat and mass integration of DAC and the PBR is critical to minimizing the MFSP. The baseline case utilizes no heat and mass integration, and the DAC system provides 100% of the CO 2 required by the photobioreactors (20 tonnes/hr), operating for 12 hours/day capturing 40 tonnes CO 2 /operating hour. The minimum fuel selling price (MFSP) of ethanol calculated from the baseline case was $10.68/gal ethanol. This corresponds with a targeted MFSP of $9.07/gal ethanol (or 15% reduction). This target was achieved by integration Option 2a with the greatest cost reduction of 17.8% (or $8.78/gal) and integration Option 2b with a cost reduction of 16.4% (or $8.93/gal). Reductions in MFSP are attributed to two primary process considerations: (a) CO 2 storage at night reduces the capital expenses associated with DAC (i.e., increasing on-stream time); and (b) distributed DAC scenarios (DAC-PBR integration Options 2a and 2b) make use of boiler and DAC CHP flue gas CO 2 (free). Direct air capture on-stream time was one of the largest contributors to MFSP reduction.

09 BIOMASS FUELS↗

Optimal Electrification Using Renewable Energies: Microgrid Installation Model with Combined Mixture k-Means Clustering Algorithm, Mixed Integer Linear Programming, and Onsset Method

Optimal planning and design of microgrids are priorities in the electrification of off-grid areas. Indeed, in one of the Sustainable Development Goals (SDG 7), the UN recommends universal access to electricity for all at the lowest cost. Several optimization methods with different strategies have been proposed in the literature as ways to achieve this goal. This paper proposes a microgrid installation and planning model based on a combination of several techniques. The programming language Python 3.10 was used in conjunction with machine learning techniques such as unsupervised learning based on K-means clustering and deterministic optimization methods based on mixed linear programming. These methods were complemented by the open-source spatial method for optimal electrification planning: onsset. Four levels of study were carried out. The first level consisted of simulating the model obtained with a cluster, which is considered based on the elbow and k-means clustering method as a case study. The second level involved sizing the microgrid with a capacity of 40 kW and optimizing all the resources available on site. The example of the different resources in the Togo case was considered. At the third level, the work consisted of proposing an optimal connection model for the microgrid based on voltage stability constraints and considering, above all, the capacity limit of the source substation. Finally, the fourth level involved a planning study of electrification strategies based mainly on microgrids according to the study scenario. The results of the first level of study enabled us to obtain an optimal location for the centroid of the cluster under consideration, according to the different load positions of this cluster. Then, the results of the second level of study were used to highlight the optimal resources obtained and proposed by the optimization model formulated based on the various technology costs, such as investment, maintenance, and operating costs, which were based on the technical limits of the various technologies. In these results, solar systems account for 80% of the maximum load considered, compared to 7.5% for wind systems and 12.5% for battery systems. Next, an optimal microgrid connection model was proposed based on the constraints of a voltage stability limit estimated to be 10% of the maximum voltage drop. The results obtained for the third level of study enabled us to present selective results for load nodes in relation to the source station node. Finally, the last results made it possible to plan electrification using different network technologies and systems in the short and long term. The case study of Togo was taken into account. The various results obtained from the different techniques provide the necessary leads for a feasibility study for optimal electrification of off-grid areas using microgrid systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Optimizing Selection Pressures and Pest Management to Maximize Cultivation Yield (OSPREY) (Final Technical Report)

This project was proposed in response to AOI 1, Cultivation Intensification Processes for Algae, within the FY19 Bioenergy Technologies Office Multi-Topic Funding Opportunity Announcement (FOA Number: DE-FOA-0002029). The work was designed to address a critical industry need to improve annualized productivity, stability, and quality of algal production strains for biofuels and bioproducts. The overall project goals were to generate process innovations rooted in established outdoor systems for strain selection, improvement, maintenance, and cultivation as well as pest detection and tracking. Planned advances included a 50% improvement in harvest yield based on AFDW (g m 2 d -1 ), 50% improvement in robustness based on stability metrics (e.g., high-productivity cultivation days, pond uptime), and 20% improvement in conversion yield. Individually, each of our planned process improvements (e.g., pest tracking) had the potential to increase productivity. However, to realize increases in yield at the system level, improvements to one unit’s process must be balanced against potential effects on other processes. For example, changes to strains, cultivation, and pest management developed in isolation may hurt other unit operations. Therefore, a critical success factor of the project was the integration of the pipeline components, achieved through iterative field-to- (short term) lab testing. In addition, through sustainability models, we evaluated how improvements would alter industry scenarios.

09 BIOMASS FUELS↗

Site Characterization of the Highest-Priority Geologic Formations for CO2 Storage in Wyoming

The project Site Characterization of the Highest-Priority Geologic Formations for CO2 Storage in Wyoming is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The Wyoming Carbon Underground Storage Project (WY-CUSP) consisted of CO2 storage site characterization and evaluation, focusing on Wyoming’s most promising CO2 storage reservoirs (the Pennsylvanian Weber/Tensleep Sandstone and Mississippian Madison Limestone) and premier CO2 storage site (Rock Springs Uplift). Results from the WY-CUSP project suggest the two reservoirs could store up to 17,000 million tons of CO2. The WY-CUSP team drilled a stratigraphic test well and acquired a 3-D seismic survey covering 25 square miles of the Rock Springs Uplift site. The team retrieved 916 feet of core from the 12,810-foot-deep well, along with a complete log suite, borehole images, fluid samples, and other data. Project partners (1) provided continuous visual documentation of the core, including grain size, mineralogy, facies distribution, and porosity; (2) performed continuous permeability and velocity scans of selected reservoir intervals; and (3) chemically analyzed the fluid samples. WY-CUSP scientists integrated seismic attributes with observations from log suites, a VSP survey, core, fluid samples, and laboratory analyses, including continuous permeability scans. From these integrations, researchers constructed 3-D spatial distribution volumes of reservoir and seal properties that represent geological heterogeneity at the targeted CO2 storage site. The WY-CUSP team used this data to perform new CO2 plume migration simulations. Baker Hughes, Inc., completed a series of small-scale, in-situ water injectivity measurements. A database was formed when observations, analyses, and experiments from the stratigraphic test well were integrated. Correlation of these data allowed petrophysical parameters to be extrapolated from the test well out into the storage domain (5x5 mile 3-D seismic survey volume). This resulted in an improved, realistic understanding of performance assessments for potential CO2 storage scenarios. The WY-CUSP team worked on (1) improving CO2 storage resource estimates, (2) establishing long-term integrity and permanence of confining layers, (3) designing a profitable strategy for pressure management, and (4) evaluating the utilization of stored CO2 at the Rock Spring Uplift. Finally, Baker Hughes developed a microseismic baseline for the test site using in-bore geophones to complete field operations.

3-D seismic↗

Electricity Costs and Carbon Implications for CO 2 -to-Fuels in Selected Locations in 2030

With the growing interest in converting carbon dioxide (CO 2 ) to fuels and products to reduce overall greenhouse gas emissions and extend carbon from biogenic and other sources, the development of carbon capture and utilization (CO2U) technologies and industry is crucial. The Markets, Resources, and Environmental and Energy Justice of CO 2 -to-Fuels Technologies (short title: MarkeRs & EEJ) project supports this goal by assessing the resource and market potential and infrastructure requirements for mid-term (-2030) and long-term (-2050) deployment of CO2U technologies. The overall project analyzes CO2U economic and resource requirements as well as sustainability and environmental and energy justice (EEJ) metrics. This report focuses on short-term (-2030) aspects of the project. It identifies three CO2U locations in the Midwestern United States and quantifies the costs for potential electricity resources and marginal emissions for each site in the near term (e.g., 2030). We focus on four electricity purchase options, including retail rate, physical power purchase agreement (PPA), financial PPA, and real-time pricing (RTP), for each of the three sites and consider potential policies such as production tax credits (PTCs). We also assess time-dependent hourly marginal electricity costs and the marginal emission rates using modeled scenarios from Cambium (Gagnon et al. 2021).The results of the analysis are expected to inform other projects within the U.S. Department of Energy's CO 2 -to-Fuels Consortium, especially the Economics and Sustainability of CO 2 Utilization Technologies with Techno-Economic Analysis and Life Cycle Analysis (TEA/LCA) project, led by Michael Wang of Argonne National Laboratory and Ling Tao of the National Renewable Energy Laboratory. The project results are also expected to inform the investment and technology communities and policymakers at the U.S. Department of Energy, state, and regional levels and guide investment by government and industry in research and development portfolios.

09 BIOMASS FUELS↗

Pilot-Scale Testing of an Integrated Circuit for the Extraction of Rare Earth Minerals and Elements from Coal and Coal Byproducts Using Advanced Separation Technologies

The primary objective of this project was to develop and demonstrate an integrated pilot-scale circuitry for recovering high-value rare earth elements (REEs) from coal and coal byproducts. The target performance was to produce a mixed REE product with content of at least two percent by weight on a dry mass basis in a cost-effective and environmentally benign manner. During the first nine months of the project period (Phase 2 Budget Period 2), pilot plant construction was completed including all field site startup activities such as permitting, engineering design, procurement/bidding, unit fabrication, site construction, equipment installation, module assembly, safety training, and circuit shakedown. During the remaining 21 months of project period (Phase 2 Budget Period 3), detailed field-testing activities were performed including feedstock sample collection and preparation, exploratory testing, circuit modification, detailed parametric study, and performance optimization. A detailed techno-economic analysis was performed based on the pilot plant testing findings which provided various scenarios for REE production. The project successfully accomplished the proposed target performance by producing mixed rare earth oxide (REO) with greater than 90% purity by weight in a continuous pilot scale operation from two distinctly different coarse refuse materials (i.e., West Kentucky No. 13 and Fire Clay coal seams), and at least three secondary sources (i.e., heap leach process and naturally formed acid mine drainage system). Project partners included the University of Kentucky, Virginia Tech, West Virginia University, Alliance Coal, Blackhawk Mining, Mineral Refining Company, and Mineral Separation Technologies. The pilot scale test facility was constructed at a former mining complex owned by Alliance Natural Resource Partners (Alliance Coal). The site was rehabilitated to accommodate the equipment installation, construction and fabrication, electrical power requirement, water line management and containment. The process units constructed and installed included X-ray sorting unit, crushing and grinding unit, physical separation unit, acid leaching unit, solvent extraction unit, and wastewater management unit. A rare earth mineral concentration unit was constructed as a standalone unit for flexible operation. A detailed environmental assessment and control plan was carried out to identify and quantify any potential impacts of the pilot-scale processing circuitry on the human and eco-system health and well-being. Corresponding mitigation strategies and control measures were provided. A conceptual flowsheet was developed to effectively remove thorium and uranium from high purity rare earth oxide mix or any potential radionuclide enriched stream. The two distinct feedstock materials were secured from the Blackhawk Mining Complex in eastern Kentucky where the Fire Clay (Hazard No. 4) seam is processed. The West Kentucky No. 13 (Baker) coarse refuse material was collected from an active process stream at an Alliance coal preparation plant located in western Kentucky. Characterization analysis indicated that both of feed materials generated from the two sources contained >300 ppm of TREEs on a dry whole mass basis which met the requirements for a qualified feed stock. The two feedstocks were further upgraded using a dual x-ray sorter to prepare the feed material for hydrometallurgical circuit. Thermal treatment on feed material prior to leaching was found to: 1) improve the leaching recovery of REEs, 2) increase the leaching kinetics, and 3) allow the leaching reaction to occur at lower acidity. Roasting at 600°C was selected as the pre-treatment condition for both West Kentucky No. 13 and Fire Clay coarse refuse material. Over 40% of leaching recovery was achieved by roasting West Kentucky No. 13 material having a top particle size of 3 mm in the pilot scale operation using 1.2M sulfuric acid leaching at 75OC. Initial pilot scale testing involved continuous operation of the pilot plant for 94 hours. The leaching unit was operated at solid-to-liquid ratio of 1 to 10 (w/v) using 0.5M sulfuric acid solution at a temperature of 75°C. The continuous solvent extraction circuit utilized rougher and cleaner units with DEHPA and TBP as the extractants. An innovative stripping circuit was developed to accumulate the REE concentration in the stripping solution to a level above 600 ppm. A bleed stream from the recycled strip solution was treated using oxalic acid precipitation which produced a high grade rare earth oxalate. The oxalate product was roasted to remove the oxalate which produced a rare earth oxide product having a purity greater 90%. Due to high concentrations of contaminant ions in the pregnant leach solution (PLS), a modified flowsheet was developed that involved pre-concentration of the REEs using multiple stages of precipitation and redissolution. The advantage of this process was improved removal of contamination before the downstream purification process and a significant cost reduction relative to the circuit that utilized the solvent extraction process. The modified circuitry included processes involving leaching, multistage precipitation, redissolution, and oxalate precipitation followed by roasting of the oxalate product. The circuit produced a mixed REO that was 92.96% pure from the initial test. A detailed parametric test plan was carried out which involved varying key parameters including solids feed rate, acid flowrate, acid concentration, multistage precipitation pH, redissolution pH, oxalate precipitation dosage and pH. The response variables included REE recovery, contaminant recovery, REO product grade and overall chemical consumption. Test results indicated that the acid-to-solid ratio is the key parameter to leaching efficiency as performance deteriorated with an increase in solids concentration. The optimal pH determined for REE precipitation and redissolution was 6.5 and 2.5, respectively. Additional tests were conducted to further improve the flowsheet. Recirculating a portion of the PLS to the feed of the leach tanks improved the leaching performance by lowering the pH of the leaching system and reducing the contamination recovery by shortening the residence time. Moreover, the removal of Al prior to REE precipitation significantly reduced the oxalic acid consumption in the oxalate precipitation circuit. The modified circuit produced over 90% grade REO by weight from both West Kentucky No. 13 and Fire Clay coarse refuse material in pilot scale continuous test programs. A case analysis model was developed to project the REE and major contaminants concentration in each PLS stream based on the leaching condition, pH cut point, and oxalic acid dosage. A correlation was established using empirical and semi-empirical models. Using the models, chemical consumption required for each stage was predicted based on the projected performance of the hydrometallurgy circuit. After identifying the optimum conditions, validation tests were carried out for the treatment of both West Kentucky No. 13 and Fire Clay coarse refuse materials in the pilot plant. The actual circuit performance and chemical consumptions were very close to the model predictions. Other than the two coarse refuse sources, several secondary feedstocks were also tested in the pilot plant facility. A “heap leach” system was constructed using the coarse refuse material generated from cleaning the West Kentucky No. 13 seam coal. Using the two stage SX rougher and cleaner circuit, a concentrate with a grade >90% REO was produced while recovering >97% of the REEs from the heap leach PLS. Naturally generated acid mine drainage (AMD) from West Kentucky No.13 mine was processed using the multistage precipitation circuit in the pilot plant in a test conducted for a period of 32 hours. The final grade of the mix RE oxide produced from the AMD was 90.84% with an overall circuit recovery of 64%. The primary source of REE was the selective precipitation steps involving iron and aluminum rejection. The hydrophobic-hydrophilic separation (HHS) process was proven to effectively recover coal from fine waste materials. For REM recovery, the HHS process was able to produce concentrates at grades of approximately 1.8% REE on an ash basis; however, recovery values were typically low, <10%, under the optimal conditions determined in the laboratory-scale testing. Staged testing of the pilot-scale HHS process for coal recovery and semi-continuous laboratory testing for REM testing showed that a total concentration ratio of more than 15x was observed for the REM recovery process. A circuit simulation package was developed for REE extraction and purification using a spreadsheet-based platform (Microsoft Excel). The REESim circuit simulation package is configured to track the mass and volume flows of components passing through a series of unit operations specified and configured by the user. The mass rates can then be utilized by the user to determine important performance indicators such as product mass yields, concentrate purity levels, element-by-element recoveries, and so forth. The techno-economic analysis showed that the roasting and leaching operations were the most expensive capital items, each contributing approximately 30% to the total capital cost. One notable contributor to the high production costs was the low REE recovery observed in the pilot scale trials. The product basket price was shown to have a strong influence on the economic viability of the scenarios, with the scandium price being the most significant influencer. Operating cost was shown to be extremely sensitive to REE recovery, REE feed grade, and leaching acid consumption. An analysis of ten different scenarios for a 500 t/h commercial operation revealed that three were economically favorable, producing internal rates of return varying from 27.7% to 33.1% and payback periods of 4 to 5 years. The project successfully developed and demonstrated a process to recover REEs from coal and coal byproducts in a pilot-plant operation which consistently produced over 90% grade REO mix from varies types of feedstocks. Commercialization analysis showed that the technology readiness level successfully achieved TRL 6 at the end of the project and demonstrated the need and the potential for scaling the process to further advance the technologies toward the goal of providing a domestic supply of REEs at a commercial scale.

01 COAL, LIGNITE, AND PEAT↗

Explainable AI (XAI)-driven vibration sensing scheme for surface quality monitoring in a smart surface grinding process

Local Interpretable and Model-agnostic Explanation (LIME), an explainable artificial intelligence (XAI) approach is adapted to identify the globally important time-frequency bands for predicting average surface roughness (Ra) in a smart grinding process. The smart grinding setup consisted of a Supertech CNC precision surface grinding machine, instrumented with a Dytran piezoelectric accelerometer attached to the tailstock quill along the tangential direction (Y-axis). For every grinding pass, vibration signatures were captured, and the ground truth surface roughness values were recorded using a Mahr Marsurf M300C portable surface roughness profilometer. The roughness values ranged from 0.06 to 0.14 microns over the complete set of experiments. Time-frequency domain spectrogram frames were extracted for each of the vibration signals collected during the grinding process. Convolutional Neural Networks (CNNs) were modeled to predict the surface roughness based on these spectrogram frames and their image augmentations. The best CNN model was able to predict the roughness values with an overall R2-score of 0.95, training R2-score of 0.99, and testing R2-score of 0.81 with only 80 sets of vibration signals corresponding to 4 experiments with 20 trials each. Although the data size is not large enough to guarantee such performance metrics in real-world scenarios, one can extract statistically consistent explanations underlying the relationships these complex deep learning models capture. Further, the LIME methodology was implemented on the developed surface roughness CNN model to identify the important time-frequency bands (i.e., the superpixels of a spectrogram) influencing the predictions. Based on the identified important regions on the spectrogram frames, the corresponding frequency characteristics were determined that influence the surface roughness predictions. The important frequency range based on LIME results was approximately 11.7 to 19.1 kHz. The power of XAI was demonstrated by cutting down the sampling rate from 160 kHz to 30, 20, 10, and 5 kHz based on the important frequency range and considering Nyquist criteria. Separate CNN models were developed for these ranges by only extracting time-frequency contents below their corresponding Nyquist cut-offs. A proper data acquisition strategy is proposed by comparing the model performances to argue the selection of a sufficient sampling rate to capture the grinding process successfully and robustly.

42 ENGINEERING↗

Reprint of: Explainable AI (XAI)-driven vibration sensing scheme for surface quality monitoring in a smart surface grinding process

Local Interpretable and Model-agnostic Explanation (LIME), an explainable artificial intelligence (XAI) approach is adapted to identify the globally important time-frequency bands for predicting average surface roughness (Ra) in a smart grinding process. The smart grinding setup consisted of a Supertech CNC precision surface grinding machine, instrumented with a Dytran piezoelectric accelerometer attached to the tailstock quill along the tangential direction (Y-axis). For every grinding pass, vibration signatures were captured, and the ground truth surface roughness values were recorded using a Mahr Marsurf M300C portable surface roughness profilometer. The roughness values ranged from 0.06 to 0.14 microns over the complete set of experiments. Time-frequency domain spectrogram frames were extracted for each of the vibration signals collected during the grinding process. Convolutional Neural Networks (CNNs) were modeled to predict the surface roughness based on these spectrogram frames and their image augmentations. The best CNN model was able to predict the roughness values with an overall R2-score of 0.95, training R2-score of 0.99, and testing R2-score of 0.81 with only 80 sets of vibration signals corresponding to 4 experiments with 20 trials each. Although the data size is not large enough to guarantee such performance metrics in real-world scenarios, one can extract statistically consistent explanations underlying the relationships these complex deep learning models capture. Further, the LIME methodology was implemented on the developed surface roughness CNN model to identify the important time-frequency bands (i.e., the superpixels of a spectrogram) influencing the predictions. Based on the identified important regions on the spectrogram frames, the corresponding frequency characteristics were determined that influence the surface roughness predictions. The important frequency range based on LIME results was approximately 11.7 to 19.1 kHz. The power of XAI was demonstrated by cutting down the sampling rate from 160 kHz to 30, 20, 10, and 5 kHz based on the important frequency range and considering Nyquist criteria. Separate CNN models were developed for these ranges by only extracting time-frequency contents below their corresponding Nyquist cut-offs. A proper data acquisition strategy is proposed by comparing the model performances to argue the selection of a sufficient sampling rate to capture the grinding process successfully and robustly.

47 OTHER INSTRUMENTATION↗

Enhancing nuclear data validation analysis by using machine learning

In this work, we demonstrate how machine learning (ML) techniques can be used as an auxiliary tool for nuclear data validation analysis. The output of the ML analysis can inform evaluators and validators on the quality, or lack thereof, of specific nuclear data and benchmark experiments with respect to simulating these experimental benchmark values. To this end, measured and simulated effective neutron multiplication factors, k eff , of 875 selected ICSBEP critical assemblies and the latter's sensitivities with respect to nuclear data as well as benchmarks' features (e.g., material nuclide, core geometry) are used as input for a random forest (RF) regression model. The RF is used to encode the complex inter-dependencies between thousands of nuclear data features (i.e., sensitivity profiles and aspects describing the measurements) and 875 simulated and experimental k eff values in order to understand what nuclear data features are most informative for predicting bias. The complexity of relationships and high-dimensional space is difficult-to-impossible to search fully using simply expert judgment. As a first proof-of-concept—step, it is demonstrated that this technique is able to correctly trace large discrepancies between simulated and experimental k eff back to fabricated shortcomings in nuclear data that were perturbed to simulated k eff values. In a second, real-case scenario, step, the RF algorithm is used to validate the ENDF/B-VIII.0 library in comparison to ENDF/B-VII.1 nuclear data. One case is showcased where the chosen ML algorithms highlighted nuclear data (the 19 F(n,inl) cross section from 0.4–0.9 MeV) that are shown to be problematic by comparing them to associated differential experimental data and nuclear data from other libraries. In addition to that, it is shown that the RF results point towards poor benchmark experiments and associated underestimated uncertainties (e.g., the PU-SOL-THERM-028 series). However, using the RF algorithm for validating nuclear data with respect to k eff is currently limited to pinpointing groups of questionable nuclear data due to the inherent correlations between features introduced by the nuclear data themselves and how k eff is simulated. Due to this, we recommend that the ML methods presented be used to augment—rather than replace—the expert knowledge of evaluators and validators.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Modeling Summary of ASNF in DOE Sealed Standard Canisters

A pathway for road-ready and final disposition packaging configurations for the aluminum-clad spent nuclear fuel (ANSF) fuel dictates storage within helium backfilled sealed DOE standard canisters. The typical packaging configuration for the 15-foot DOE canisters places 10 advance test reactor (ATR) elements within a Type 1a basket, and three baskets are loaded within each DOE canister. During in-reactor operations and cooling pond storage conditions, oxyhydroxide layers form on the surface of the aluminum clad fuel. These layers produce hydrogen gas over time due to the fuel’s radiation field. As part of the packing procedure, the ATR fuel should be dried to remove any residual physio-/chemi- sorbed water bound to the surface. The results of this modeling will include results at fully saturated and fully dried conditions. In addition, fuels that are currently stored at the Savannah River Site were also studied for their potential for hydrogen and pressure build up. These two additional fuels modeled were the Missouri University Research Reactor (MURR) fuel, which is packaged in the same configuration as the ATR, but with a 10-foot-tall DOE standard canister. This was selected due to its relatively high decay heat compared to other DOE-managed ASNF. The second additional fuel studied with the modeling effort was the High Flux Isotope Reactor (HFIR) fuel. This fuel is modeled as two separate DOE canisters with the inner and outer annulus split for storage. The HFIR was selected for study due to its high aluminum cladding surface area. In the associated experimental work, Task 2 - Oxyhydroxide Layer Radiolytic Gas Generation Resolution, additional tests were completed in a helium environment, and updated G-values for the radiolytic production of hydrogen from the oxyhydroxide layers were provided. These values were 2.92 ×10 -4 µmol/J at 50% relative humidity and 4.12 ×10 -4 µmol/J at 100% relative humidity. These values are lower than the value for Argon that was used in prior modeling results. In addition, prior modeling results have been completed with the G-value applied to just the mass of the corrosion layer, and this has been updated to apply the G-value for hydrogen generation to the full mass of the fuel. These two effects combine to show much smaller pressure and hydrogen build up for the sealed canister model. For a nominal scenario of stored ATR fuel, after 50 years the model results give a 1.36 atm total pressure, 7% mole percent hydrogen, for the upper decay heat, 1.51 atm total pressure, 16% mole percent hydrogen, and for upper decay heat with undried fuel 2.6 atm total pressure, 15% mole percent hydrogen. For the MURR nominal case, the model results give 1.34 atm and 6% hydrogen, for upper decay heat this gives 1.41 atm total pressure with 10.8 % hydrogen, and for upper decay heat with undried fuel, this gives 2.38 atm total pressure with 9.9% hydrogen. The nominal scenario for HFIR fuel gives 1.39 atm total pressure with 9.9% hydrogen, the upper decay heat case gives 1.43 atm with 12.1% hydrogen, and the upper decay heat with undried fuel gives 2.17 atm total pressure with 11.9% hydrogen. These results confirm the ATR scenario bounds the other intact ASNF modeled here. No case modeled yields significant oxygen, and the lower decay heat cases for all fuels modeled have hydrogen concentrations that are under the 4% flammability limit after 50 years of storage. In addition, the modeled pressures for all cases are all significantly below the 500-psi limit for the DOE standard sealed canister.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗

High-Resolution Fire Weather Index Data for the Conterminous US (1980–2099), Version 1

This dataset presents a suite of high-resolution fire weather index datasets calculated from observation (gridMet, Livneh, Daymet V4), reanalysis (AgERA5), downscaled hydro-climate projections over the conterminous United States (CONUS) based on multiple selected Global Climate Models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). Aside from the daily FWI datasets, we also include a set of FWI extreme indicators at annual, seasonal, and monthly scales, including 1) fwixx: maximum FWI; 2) fwisa: mean FWI. For annual fwisa, it refers to the season with the maximum seasonal average; 3) fwils: Length of fire season over a specified period, where fire season is defined as the days exceeding the median value of the normalized FWI during the reference period (1980-1984); 4) fwixd: Number of extreme fire weather days over a specified period, where extreme day is defined as the day with FWI > the 95th percentile of the FWI during the reference period (1980-1984). All FWI datasets cover 1980-2020 baseline and the model simulated products including the downscaled products additionally include 2021-2099 near-future periods under the high-end (SSP585) emission scenario.

54 ENVIRONMENTAL SCIENCES↗

HEAT TRANSFER ASPECTS OF TWO PISTON SPLAT QUENCHING

This thesis investigates the heat transfer aspects of the two-piston splat quench system by modifying surface roughness and by using copper and stainless steel platens. In additive manufacturing, specifically selective laser melting (SLM), metal is subjected to a very rapid heating and cooling that produces cooling rates in the range from 105-106 °C/s. This rapid solidification produces microstructures that deviate far from equilibrium. Splat quenching is an experimental method that produces rapidly solidified structures under cooling rates that are comparable to those found in additive manufacturing. The splat quenching approach allows for more efficient testing of the rapid solidification behavior of a wide variety of alloys. In this thesis, two platen characteristics, surface roughness and material thermal conductivity, were examined to determine their role in controlling the heat transfer and the subsequent cooling rate during splat quenching of austenitic stainless steel 316L. To that objective, 5 surface finishes were produced using metallographic techniques, from 7046-203 Å Rms. Copper and 316L stainless steel platen material were used to reach different levels of thermal conductance. The combination of these 2 variables created 10 scenarios, which were examined using, backscatter imaging, secondary imaging, and electron back scatter diffraction. Cooling rates were estimated using cell sizing obtained through electrolytic etching and secondary imaging. This thesis determined that both surface modification and changes in material type had a statistically significant impact on the cell sizes observed. Where material type far outweighed surface roughness in its effect on cell size.

Roze, Jonathan↗

Doubling off-axis electron cyclotron current drive efficiency via velocity space engineering

For the first time, experiments on the DIII-D tokamak have demonstrated electron cyclotron current drive with more than double the conventional efficiency by tailoring the wave–particle interactions in velocity space using a novel ‘top launch’ geometry. Steering the EC waves to propagate nearly parallel to the resonance drives current more efficiently by (1) selective damping on electrons with higher parallel velocity v ||, and (2) longer absorption path to compensate for inherently weaker absorption at higher v ||. Experiments using a fixed-injection top launch system find an optimal velocity space interaction for maximum current drive efficiency at ρ ~ 0.5 where the ease of drawing out a high v || electron tail is balanced by sufficient absorption.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spatially resolved gas and stellar kinematics in compact starburst galaxies

The kinematics of galaxies provide valuable insights into their physics and assembly history. Kinematics are governed not only by the gravitational potential, but also by merger events and stellar feedback processes such as stellar winds and supernova explosions. We aim to identify what governs the kinematics in a sample of SDSS-selected nearby starburst galaxies, by obtaining spatially resolved measurements of the gas and stellar kinematics. We obtained near-infrared integral-field K-band spectroscopy with VLT/SINFONI for 15 compact starburst galaxies. We derived the integrated as well as spatially resolved stellar and gas kinematics. The stellar kinematics were derived from the CO absorption bands, and Paα and Brγ emission lines were used for the gas kinematics. Based on the integrated spectra, we find that the majority of galaxies have gas and stellar velocity dispersion that are comparable. A spatially resolved comparison shows that the six galaxies that deviate show evidence for a bulge or stellar feedback. Two galaxies are identified as mergers based on their double-peaked emission lines. In our sample, we find a negative correlation between the ratio of the rotational velocity over the velocity dispersion (v rot /σ) and the star formation rate surface density. We propose a scenario where the global kinematics of the galaxies are determined by gravitational instabilities that affect both the stars and gas. This process could be driven by mergers or accretion events. Effects of stellar feedback on the ionised gas are more localised and detected only in the spatially resolved analysis. The mass derived from the velocity dispersion provides a reliable mass even if the galaxy cannot be spatially resolved. The technique used in this paper is applicable to galaxies at low and high redshift with the next generation of infrared-focussed telescopes (JWST and ELT).

79 ASTRONOMY AND ASTROPHYSICS↗

Integrating static PRA information with risk informed safety margin characterization (RISMC) simulation methods

The overall objective of the project was to develop a computationally feasible and user-friendly mechanized process to integrate traditional probabilistic risk assessment (PRA) and dynamic PRA (DPRA) results. Starting with the systematic identification of items in an existing PRA that need dynamic augmentation, the project used a generic 4-loop pressurized reactor (PWR) and 3-loop PWR as example plants. Station blackout (SBO) and large break loss of coolant accident SBLOCA) were selected as the example initiating events. Using the traditional event-tree (ET)/fault-tree (FT) methodology augmented by dynamic evet tree approach, the potential consequences of the initiating events were simulated with RELAP-3D and MELCOR/RASCAL codes to cover Level 1 through Level 3 of PRA. RAVEN and ADAPT software were used to generate Level 1 simulations with RELAP-3D and Level 2/3 simulations with MELCOR (Level 2)/RASCAL (Level 3), respectively. Example branching conditions (BCs) for SBO included AC power recovery time, valve repair failure time, reactor coolant pump leak time/break size and emergency power supply duration to a total of 9. Example BCs for LOCA included off-site power recovery time, diesel generator power recovery time, auxiliary feed water system operation time, safety relief valve failure to open upon demand, reactor coolant pump seal break time and size to a total of 21. Each RELAP-3D simulation (9,587 scenarios) was labelled OK or Core Damage based on the maximum allowed peak clad temperature (2,100oF). Each MELCOR simulation (4610 scenarios) was labeled as Bin over 10rem or Bin 0-10rem based on the dose at the site boundary. The scenarios were clustered based on the criteria above using the mean shift methodology. Classical PRA (CPRA) and DPRA results were compared to identify the ET sequences that need DPRA augmentation. Several approaches were proposed for the incorporation of these sequences into CPRA using clustering with the mean shift methodology, restructuring the CPRA ETs by adding new BCs/sequences, and using the concept of a limit surface. Procedures for decision making regarding the possible consequences of an initiating event (e.g. core damage or not, site evacuation or not) were developed using a convolutional neural network (CNN), a recurrent neural network (RNN) and a transformer neural network (TNN). The project has led to two PhD degrees, three archival journal papers and five refereed conference proceedings.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Scalable Techniques for Stochastic Power Flow Problems (Final Report)

The proposed research focuses on developing scalable algorithms for two-stage security-constrained OPF problems with AC power flow constraints, a class of problems complicated by (i) scale arising from a scenario representation; and (ii) the presence of nonlinearity, nonconvexity, and possibly second-stage discreteness or complementarity. Unfortunately, most existing solvers cannot contend with both challenges simultaneously; accordingly, the proposed research focuses on developing solution techniques that can both scale with the number of scenarios and contend with nonconvexity and second-stage complementarity. We consider three avenues for addressing such problems: (i) Variable sample-size SQP (VS-SQP) methods that combine sparse Quasi-Newton updates with a scalable variance-reduced stochastic gradient scheme for stochastic QP subproblems, allowing for contending with second-stage complementarity via regularization; (ii) Variable sample-size stochastic Interior-point (VS-sIP) schemes that propose a sampling-based regularized (to allow for contending with complementarity) interior-point schemes in which a Schur-complement technique is employed for decomposing the Newton direction computation step; (iii) Variable sample-size tractable ADMM (VS-tADMM) schemes combine variable sample-sizes with carefully designed techniques for resolving each of the nonconvex updates (by leveraging the QCQP structures). We intend to compare the three schemes using performance profiles in terms of solution quality, scalability, etc. and then select one scheme which will then be developed and further refined in Python for purposes of the GO competition.

42 ENGINEERING↗

Large Scale Caching and Streaming of Training Data for Online Deep Learning

The training of deep neural network models on large data remains a difficult problem, despite progress towards scalable techniques. In particular, there is a mismatch between the random but predetermined order in which AI flows select training samples and the streaming I/O patterns for which traditional HPC data storage (e.g., parallel file systems) are designed. In addition, as more data are obtained, it is feasible neither simply to train learning models incrementally, due to catastrophic forgetting (i.e., bias towards new samples), nor to train frequently from scratch, due to prohibitive time and/or resource constraints. In this paper, we study data management techniques that combine caching and streaming with rehearsal support in order to enable efficient access to training samples in both offline training and continual learning. We revisit state-of-art streaming approaches based on data pipelines that transparently handle prefetching, caching, shuffling, and data augmentation, and discuss the challenges and opportunities that arise when combining these methods with data-parallel training techniques. We also report on preliminary experiments that evaluate the I/O overheads involved in accessing the training samples from a parallel file system (PFS) under several concurrency scenarios, highlighting the impact of the PFS on the design of the data pipelines.

data pipelines↗