Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase object approximation”

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 127 records · Page 7

Initial Demonstration of New Griffin Technologies for Simulating the Running-In Phase of Pebble Bed Reactors

Griffin is a reactor multiphysics modeling application based on MOOSE (Multiphysics Object-Oriented Simulation Environment) and specifically targeting transient modeling of advanced reactors. Griffin has been used recently to model pebble-bed reactors for the Nuclear Regulatory Commission (NRC) Office of Nuclear Regulatory Research and the Advanced Reactor Technology program. This modeling work has focused thus far on the direct calculation of equilibrium cores. This report documents an initial demonstration of a new running-in simulation capability. The new running-in capability is verified using the existing direct equilibrium core calculation capability. A simplified pebble-bed reactor model is then used to demonstrate the running-in simulation capability. This demonstration shows that Griffin is able to simulate years of operation during the running-in phase efficiently with each depletion step taking only several seconds. Two new technologies are also presented in this report which have been developed in Griffin that will be essential for improved accuracy both of the direct equilibrium core computation and the new running-in simulation capability. The first technology is an online cross section generation capability specifically targeted for pebble-bed reactors. This will improve the accuracy of the depletion calculation as the cross sections are generated at the exact core status. This also avoids the difficult step of pre-generating a separate standalone multigroup cross section set. Secondly, a newly implemented discretization for discontinuous finite element method (DFEM) SN transport in cylindrical (RZ) coordinates, which can be solved efficiently using the existing SN sweep solver, is discussed and some results are shown demonstrating the usefulness of the additional accuracy transport provides over a diffusion approximation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Formulation and validation of a computational model for a dilute biomass slurry undergoing rotational mixing

In this paper we develop a computational model for the mixing and transport of a dilute biomass slurry. The objective was to create a sufficiently simple and efficient model for biomass transport that can be coupled with reaction models for the study of conversion of cellulosic biomass into fermentable sugars. Our target system is 5%-by-mass ..alpha..-cellulose, which is our proxy for more complex lignocellulosic biomass. In the authors' previous work, an experimental investigation with ..alpha..-cellulose under two vane-mixer configurations showed a bifurcation between a settling regime, for which settling effects dominate, and a suspended regime, for which solids are mostly suspended. Here, a mixed-fluid model was chosen, for which the model for the mixture-velocity field is the incompressible Navier-Stokes equations under the Boussinesq approximation for buoyancy. Solids transport includes solids motion due to diffusion, settling, advection, and shear. Comparison of simulated and experimental results show good agreement in the suspended regime, and in capturing the bifurcation rate. While the model captured well the distribution of solids in the settling regime, the model was incapable of capturing the high torque values seen in experiments with vanishing mixer rotation rate.

09 BIOMASS FUELS↗

Final Report to the National Energy Technology Laboratory on FY14- FY15 Cooperative Research with the Consortium for Electric Reliability Technology Solutions

The Power System Engineering Research Center (PSERC) engages in technological, market, and policy research for an efficient, secure, resilient, adaptable, and economic U.S. electric power system. PSERC, as a founding partner of the Consortium for Electric Reliability Technology Solutions (CERTS), conducted a multi-year program of research for U.S. Department of Energy (DOE) Office of Electricity Delivery and Energy Reliability (OE) to develop new methods, tools, and technologies to protect and enhance the reliability and efficiency of the U.S. electric power system as competitive electricity market structures evolve, and as the grid moves toward wide-scale use of decentralized generation (such as renewable energy sources) and demand-response programs. Phase I of OE’s funding for PSERC, under cooperative agreement DE-FC26-09NT43321, started in fiscal year (FY) 2009 and ended in FY2013. It was administered by DOE’s National Energy Technology Laboratory (NETL) through a cooperative agreement with Arizona State University (ASU). ASU provided sub-awards to the participating PSERC universities. This document is PSERC’s final report to NETL on the activities for OE, conducted through CERTS, from September 2015 through September 2017 utilizing FY 2014 to FY 2015 funding under cooperative agreement DE-OE0000670. PSERC is a thirteen-university consortium with over 30 industry members. Since 1996, PSERC has been engaged in research and education efforts with the mission of “empowering minds to engineer the future electric energy system.” Its work is focused on achieving: • An efficient, secure, resilient, adaptable, and economic electric power infrastructure serving society • A new generation of educated technical professionals in electric power • Knowledgeable decision-makers on critical energy policy issues • Sustained, quality university programs in electric power engineering. PSERC core research is funded by industry, with a budget supporting approximately 30 principal investigators and some 70 graduate students and other researchers. Its researchers are multi-disciplinary, conducting research in three principal areas: power systems, power markets and policy, and transmission and distribution technologies. The research is collaborative; each project involves researchers typically at two universities working with industry advisors who have expressed interest in the project. Examples of topics for recent PSERC research projects include grid integration of renewables and energy storage, new tools for taking advantage of increased penetration of real-time system measurements, advanced system protection methods to maintain grid reliability, and risk and reliability assessment of increasingly complex cyber-enabled power systems. A PSERC’s objective is to proactively address the technical and policy challenges of U.S. electric power systems. To achieve this objective, PSERC works with CERTS to conduct technical research on advanced applications and investigate the design of fair and transparent electricity markets; these research topics align with CERTS research areas 1 and 2: Real-time Grid Reliability Management (Area 1), and Reliability and Markets (Area 2). The CERTS research areas overlap with the PSERC research stems: Power Systems, Power Markets, and Transmission and Distribution Technologies, as described on the PSERC website (see http://www.pserc.org/research/research_program.aspx). The performers were with Arizona State University (ASU), Cornell University (CU), University of California at Berkeley (UCB), and University of Illinois at Urbana-Champaign (UIUC). PSERC research activities in the area of reliability and markets focused on electric market and power policy analyses. The resulting studies suggest ways to frame best practices using organized markets for managing U.S. grid assets reliably and to identify highest priority areas for improvement. PSERC research activities in the area of advanced applications focused on mid- to long-term software research and development, with anticipated outcomes that move innovative ideas toward real-world application. Under the CERTS research area of Real-time Grid Reliability Management, PSERC has been focused on Advanced Applications Research and Development (AARD), a subgroup of activities that works to develop advanced applications and tools to more effectively operate the electricity delivery system, by enabling advanced analysis, visualization, monitoring and alarming, and decision support capabilities for grid operators.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Planck intermediate results: XL. The Sunyaev-Zeldovich signal from the Virgo cluster

The Virgo cluster is the largest Sunyaev-Zeldovich (SZ) source in the sky, both in terms of angular size and total integrated flux. Planck’s wide angular scale and frequency coverage, together with its high sensitivity, enable a detailed study of this big object through the SZ effect. Virgo is well resolved by Planck, showing an elongated structure that correlates well with the morphology observed from X-rays, but extends beyond the observed X-ray signal. We find good agreement between the SZ signal (or Compton parameter, y c ) observed by Planck and the expected signal inferred from X-ray observations and simple analytical models. Owing to its proximity to us, the gas beyond the virial radius in Virgo can be studied with unprecedented sensitivity by integrating the SZ signal over tens of square degrees. In this paper, we study the signal in the outskirts of Virgo and compare it with analytical models and a constrained simulation of the environment of Virgo. Planck data suggest that significant amounts of low-density plasma surround Virgo, out to twice the virial radius. We find the SZ signal in the outskirts of Virgo to be consistent with a simple model that extrapolates the inferred pressure at lower radii, while assuming that the temperature stays in the keV range beyond the virial radius. The observed signal is also consistent with simulations and points to a shallow pressure profile in the outskirts of the cluster. This reservoir of gas at large radii can be linked with the hottest phase of the elusivewarm/hot intergalactic medium. Taking the lack of symmetry of Virgo into account, we find that a prolate model is favoured by the combination of SZ and X-ray data, in agreement with predictions. In conclusion, based on the combination of the same SZ and X-ray data, we constrain the total amount of gas in Virgo. Under the hypothesis that the abundance of baryons in Virgo is representative of the cosmic average, we also infer a distance for Virgo of approximately 18 Mpc, in good agreement with previous estimates.

79 ASTRONOMY AND ASTROPHYSICS↗

Progress Updates for the Hanford Waste Treatment and Immobilization Plant Project - 16281

The Waste Treatment and Immobilization Plant (WTP) Project is vital to the U.S. Department of Energy's (DOE) mission to clean up radioactive waste at the Hanford Site, located in southeastern Washington State. The overall WTP Project objective is to design, build, and commission the facilities and systems that will treat and immobilize approximately 56 million gallons of radioactive waste stored in Hanford's 177 underground storage tanks. The combination of technologies in the WTP makes it a first-of-a-kind facility at a scale larger than ever built for radioactive waste processing. The WTP includes several facilities that will contribute to waste processing: (1) the Pretreatment (PT) Facility, which is designed to receive waste retrieved from the Hanford tank farms and separate it into low-activity and high-level waste streams; (2) the Low-Activity Waste (LAW) Facility, which will vitrify the low-activity waste stream; (3) the High-Level Waste (HLW) Facility, which will vitrify the high-level waste stream; (4) the Analytical Laboratory (Lab), which will support both vitrification facilities' operations; and (5) the Balance of Facilities (BOF), which will provide industrial support capabilities and infrastructure, like electricity and steam, needed for WTP operations. In 2012, DOE suspended production work on the WTP's Pretreatment (PT) Facility and to a lesser degree the High-Level Waste (HLW) Facility because of unresolved technical issues associated with those facilities. DOE has established an approach that would facilitate immobilization of tank waste to begin as early as practicable without waiting for completion of work to resolve the technical issues associated with the PT and HLW Facilities. This approach enables LAW treatment to begin before the PT and HLW facilities are completed by transferring low activity waste directly from the Hanford tank farms to the LAW Facility - referred to as Direct Feed LAW (DFLAW). ORP has placed the highest priority on completing the Low Activity Waste Facility, Balance of Facilities, and Analytical Laboratory (collectively known as LBL) in support of the planned phased start to the Hanford Tank waste treatment mission. DOE has requested, and the WTP contractor has delivered contract modifications proposals for: (1) design of a permanent capability to accommodate DFLAW in WTP; (2) procurement and construction of modifications to enable DFLAW operations in WTP, and (3) completion of LBL work scope in the current contract through hot commissioning in a manner to accommodate a phased start toward completing the WTP mission. DOE is reviewing and approving contract modifications in a sequence and schedule that allows for an orderly completion of the LBL Facilities and to begin treating low-activity tank waste as soon as practicable. In August 2014, ORP authorized the WTP contractor to resume all engineering work necessary to finalize the design of the HLW Facility. Authorization for HLW equipment procurements and construction is forthcoming, based on conditions that must be satisfied for a full release of these activities, including completion of an update to the Preliminary Documented Safety Analysis. DOE has also established a comprehensive plan for resuming production activities for the PT Facility that is similar to the process established and implemented for the HLW Facility. Resolution of the technical issues and establishment of a revised safety basis are key prerequisites to resuming engineering, procurements, and construction. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Final Scientific/Technical Report: Characterizing the Response of the Cascadia Margin Gas Hydrate Reservoir to Bottom Water Warming Along the Upper Continental Slope

The objective of this project is to understand the response of the WA margin gas hydrate system to contemporary warming of bottom water along the upper continental slope. Through pre-cruise analysis and modeling of archive and recent geophysical and oceanographic data, we (1) inventoried bottom simulating reflectors along the WA margin and defined the upper limit of gas hydrate stability, (2) refined margin-wide estimates of heat flow and geothermal gradients, (3) characterized decadal scale temporal variations of bottom water temperatures at the upper continental slope of the Washington margin, and (4) used numerical simulations to provide quantitative estimates of how the shallow boundary of methane hydrate stability responds to modern environmental change. These pre-cruise results provided the context for a systematic geophysical and geochemical survey of methane seepage along the upper continental slope from 48° to 46°N during a 10-day field program on the R/V Thompson from October 10-19, 2014. This systematic inventory of methane emissions along this climate-sensitive margin corridor and comprehensive sediment and water column sampling program provided data and samples for Phase 3 of this project that focused on determining fluid and methane sources (deep-source vs. shallow; microbial, thermogenic, gas hydrate dissociation) within the sediment, and how they relate to contemporary intermediate water warming. During the 2014 research expedition, we sampled nine seep sites between ~470 and 520 m water depth, within the zone of predicted methane hydrate retreat over the past 40 years. We imaged 22 bubble plumes with heights commonly rising to ~300 meters below sea level with one reaching near the sea surface. We collected 22 gravity cores and 20 CTD/hydrocasts from the 9 seeps and at background locations (no acoustic evidence of seepage) within the depth interval of predicted downslope retreat of the methane hydrate stability zone. Approximately 300 pore water samples were extracted from the gravity cores, and the pore water was analyzed for a comprehensive suite of solutes, gases, and stable isotope ratios. This comprehensive geochemical dataset was used to characterize the fluid and gas source(s) at each of the seep sites surveyed. The primary results of this project are: 1) Bottom simulating reflector-derived heat flow values decrease from 95 mW/m2 10 km east of the deformation front to ~60 mW/m2 60 km landward of the deformation front, with anomalously low values of ~25 mW/m2 on a prominent mid-margin terrace off central Washington. 2) The temperature of the incoming sediment/ocean crust interface at the deformation front ranges between 164-179 oC off central Washington, and the 350 oC isotherm at the top of the subducting ocean crust occurs 95 km landward of the deformation front. Differences between BSR-derived heat flow and modeled conductive heat flow suggest mean upward fluid flow rates of 0.4 cm/yr across the margin, with local regions (e.g. fault zones) exhibiting fluid flow rates up to 3.5 cm/yr. 3) A compilation of 2122 high-resolution CTD, glider, and Argo float temperature profiles spanning the upper continental slope of the Washington margin from the years 1968 to 2013 show a long-term warming trend that ranges from 0.006-0.008 oC/yr. Based on this long-term bottom water warming, we developed a 2-D thermal model to simulate the change in sediment temperature distribution over this period, along with the downslope retreat of the methane hydrate stability field. Over the 43 years of the simulation, the thermal disturbance propagated 30 m into the sediment column, causing the base of the methane hydrate stability field to shoal ~13 m and to move ~1 km downslope. 4) A preliminary analysis of seafloor observations and mid-water column acoustic data to detect bubble plumes was used to characterize the depth distribution of seeps along the Cascadia margin. These results indicate high bubble plume densities along the continental shelf at water depths <180 m and at the upper limit of methane hydrate stability along the Washington margin. 5) The majority of the seeps cored during the 2014 research expedition on the R/V Thompson contained abundant authigenic carbonate indicating that they are locations of long-lived seepage rather than emergent seep systems related to methane hydrate dissociation. Despite the evidence for enhanced methane seepage at the upper limit of methane hydrate stability along the Washington margin, we found no unequivocal evidence for active methane hydrate dissociation as a source of fluid and gas at the seeps surveyed. The pore fluid and bottom water chemistry shows that the seeps are fed by a variety of fluid and methane sources, but that methane hydrate dissociation, if occurring, is not widespread and is only a minor source (below the detection limit of our methods). Collectively, these results provide a significant advance in our understanding of the thermal structure of the Cascadia subduction zone and the fluid and methane sources feeding seeps along the upper continental slope of the Washington-sector of the Cascadia margin. Though we did not find unequivocal evidence for methane hydrate dissociation as a source of water and methane at the upper pressure-temperature limit of methane hydrate stability at present, continued warming of North Pacific Intermediate Water in the future has the potential to impact the methane hydrate reservoir in sediments at greater depths along the slope. Thus, this study provides a strong foundation and the necessary characterization of the background state of seepage at the upper limit of methane hydrate stability for future investigations of this important process.

54 ENVIRONMENTAL SCIENCES↗

Investigating Spatial Variability of Aerosol, Cloud Condensation Nuclei, and Ice Nucleating Particles in Mountainous Terrain Field Campaign Report

The U.S. Department of Energy Atmospheric System Research (ASR)-supported Surface Atmosphere Integrated Field Laboratory (SAIL) campaign in the East River Watershed (ERW) of the Upper Colorado River Basin in southwestern Colorado ran from fall 2021 to spring 2023. Two monitoring sites were deployed in the ERW as part of SAIL. The two sites were the Aerosol Observation System (AOS) located on Crested Butte Ski Mountain, and the second ARM Mobile Facility (AMF2), located at the Rocky Mountain Biological Laboratory in Gothic, Colorado. To gain a more comprehensive understanding of aerosols in complex, mountainous terrain, Handix Scientific deployed SAIL-Net, a distributed network of six measurement nodes spanning the domain of the SAIL research area from October 2021 to July 2023. Each node measured aerosol particles between 140 nm and 3.4 μm in diameter using a small portable optical particle spectrometer (POPS; Gao et al. 2016), cloud condensation nuclei (CNN) using a miniature CCN counter (CloudPuck), and ice nucleating particles (INP) using the time-resolved aerosol filter sampler (TRAPS; Creamean et al. 2018). Our approach was similar to other studies that aimed to better characterize and understand aerosols and gas-phase pollutants using networks of lower-cost sensors (Caubel et al. 2019, Kelly et al. 2021, Asher et al. 2022). Such studies have identified neighborhood-level variations in pollutant concentrations (Schneider et al. 2017, Popoola et al. 2018, Caubel et al. 2019). Small-scale variations such as this are poorly represented in models and poorly measured by a single monitoring system (Caubel et al. 2019). Previous work has shown the representation error (the ability of measurements to represent a larger area) increases with complex orography, leading to decreases in model accuracy (Schutgens et al. 2017). The overall goal of SAIL-Net was to improve our understanding of the variability of aerosol in the ERW, thus increasing our knowledge of aerosol-cloud interactions in this region and informing the usefulness of distributed networks of measurements for future studies. We met this goal by answering the following science questions: 1. What is the aerosol temporal variability, and how does aerosol inhomogeneity vary seasonally? Is there significant seasonal variability in sources, or are short-term meteorological conditions the most important determining factor in sources for cloud nuclei? 2. What is the aerosol spatial variability? What are the aerosol characteristics at cloud base, presumably the particles most representative of those acting as cloud nuclei? 3. How should measurement networks be designed to capture aerosol-cloud interactions, and what do they need to measure? Can a single measurement site accurately represent aerosol properties in regions of complex terrain? SAIL-Net consisted of six measurement nodes spread across the ERW near Crested Butte, Colorado. The primary objective in site placement was to select locations that captured the vertical variation in aerosol properties while also spanning the domain of the SAIL campaign. The elevation of the sites ranged from roughly 2750 m along the valley floor of the ERW to approximately 3500 m near the top of Crested Butte Mountain, which is one of the taller peaks in the ERW. The farthest distance between sites was 14 km, while the closest two sites were approximately 1 km apart. Two of the sites were collocated with the ARM SAIL sites; our instruments sat on top of one of the trailers at AOS and another one of our sites was located in a meadow just above AMF2.

54 ENVIRONMENTAL SCIENCES↗

The linear instability of dilute ultrarelativistic e{sup ±} pair beams

The annihilation of TeV photons from extragalactic TeV sources and the extragalactic background light produces ultrarelativistic e{sup ±} beams, which are subject to powerful plasma instabilities that sap their kinetic energy. Here we study the linear phase of the plasma instabilities that these pair beams drive. To this end, we calculate the linear growth rate of the beam-plasma and oblique instability in the electrostatic approximation in both the reactive and kinetic regimes, assuming a Maxwell–Jüttner distribution for the pair beam. We reproduce the well-known reactive and kinetic growth rates for both the beam-plasma and oblique mode. We demonstrate for the oblique instability that there is a broad spectrum of unstable modes that grow at the maximum rate for a wide range of beam temperatures and wave-vector orientations relative to the beam. We also delineate the conditions for applicability for the reactive and kinetic regimes and find that the beam-plasma mode transitions to the reactive regime at a lower Lorentz factor than the oblique mode due to a combination of their different scalings and the anisotropy of the velocity dispersions. Applying these results to the ultrarelativistic e{sup ±} beams from TeV blazars, we confirm that these beams are unstable to both the kinetic oblique mode and the reactive beam-plasma mode. These results are important in understanding how powerful plasma instabilities may sap the energy of the ultrarelativistic e{sup ±} beams as they propagate through intergalactic space.

79 ASTRONOMY AND ASTROPHYSICS↗

RCCS Testing with the Water-based NSTF: Year-1 Single-Phase Results

Under support from the Department of Energy (DOE) and the Office of Advanced Reactor Technologies (ART), a large scale test facility has been constructed to generate validation data for passive decay heat removal systems of advanced reactors. Located at Argonne National Laboratory, the Natural convection Shutdown heat removal Test Facility (NSTF) reflects key features of a ½ scale, water-based, Reactor Cavity Cooling System (RCCS). This report serves as a summary of the single-phase testing results obtained during the programs first year of matrix testing. An introduction to the high-level program objectives is first presented, followed by a brief summary of facility design, instrumentation, and control systems. Details of the NQA-1 driven test method are included, which govern the programs quality assurance requirements and provide a background for the testing strategy and procedures. Checkout activities, including verification of gross system mass and energy balances, heat loss quantification, and isothermal facility characterization, were performed early in the first year and results presented in this report. These tests have identified relevant facility characteristics such as K-factors, frictional pressure losses, parasitic heat loss rates, and transient heating response. Results are then presented from four accepted matrix tests where the facility established a steady-state, single-phase mode of flow operation. All matrix tests were performed with an identical top-level procedure: first initiating a heater power ramp to target levels, and then holding power levels constant until the facilities structural members and network piping reached thermal equilibrium. Upon obtaining stable thermal and hydraulic conditions, active cooling was initiated off of the bulk storage tank to maintain constant riser inlet temperatures of 30°C. This steady-state period was then held for a minimum required duration, at which point the heaters were ramped down and test concluded. The two baseline tests completed reflect conditions prototypic to the full scale plant at both normal operation, and design basis Depressurized Loss of Forced Cooling (DLOFC) conditions. These parameters are defined by full RCCS decay heat loads of 1.4 MWt and 2.1 MWt, which translate to NSTF scale thermal power levels of 34.4 kWt and 51.6 kWt, respectively. Two additional tests were also conducted, which repeated baseline cases but included modification by single parameter variations in decay heat load, or heater power, at 700 kWt, and 2.8 MWt. Finally, a scoping test was performed to demonstrate the facilities ability to operate in a two-phase mode of operation. The facility successfully operated for a 3-hour period at this mode, reaching peak inventory temperatures of 105°C and condensing approximately 50 gallons of boil-off inventory. The testing parameters for this demonstration were based off the design basis DLOFC condition, and will be the reference for new baseline cases as the program moves forward with two-phase matrix testing.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Gas Stopper Developments for Improved Purity and Intensity of Low-Energy, Rare Isotope Ion Beams (Final Technical Report)

This final technical report summarizes the work of the Michigan State University (MSU) team supported by grant # DE-SC0021423 awarded by the Office of Nuclear Physics, Department of Energy. Objectives: The successful fulfillment of the FRIB science mission hinges on ensuring the availability of fast, stopped, and reaccelerated beams consisting of rare isotopes. This project's research and development focus was dedicated to supporting the advancement and creation of a cutting-edge linear gas stopper. The primary aim is to efficiently convert the high-intensity fast beams of rare isotopes provided by FRIB into high-quality, low-energy beams. These beams are essential for conducting stopped beam experiments or for subsequent reacceleration. The overarching goal is to advance technology, aiming to increase the beam rate capability of the linear gas stopper for medium-to-heavy-mass rare isotopes by more than tenfold compared to the currently most effective gas stopper in operation, and to improve the purity of the delivered beams. Project Description: The existing technology employed in gas stopping devices designed for low-energy, rare-isotope beams presents limitations in both the purity of extracted beams and the intensities of injected beams. These limitations are incompatible with the requirements of the recently commissioned rare isotope beam facility, FRIB. Our research and development efforts, aligned with the previously outlined objectives, focused on addressing the most critical aspects for enhancing beam-rate capability and purity. Specifically, advanced particle-in-cell simulations were developed and integrated into a simulation pipeline to explore the efficacy of multi-layer RF carpets on increasing ion transport efficiency with high incoming beam rates that generate space charge fields which can limit it. We also explored the possibility of using a collision-induced-dissociation (CID) gas cell to break up molecular contaminant ions that are generated during the stopping process. A prototype CID gas cell was constructed and tested with beams from an offline ion source, validating the concept with the successful demonstration of breaking of molecular ions. The outcome of this research enabled the formulation a conceptual design for a next-generation linear gas stopping device specifically tailored for FRIB. This device is envisioned to deliver rare-isotope-ion beams at a rate of 10 8 particles per second or higher, accompanied by advancements in purity. Methods employed: This project leverages advancements in technologies initially designed for the Advanced Cryogenic Gas Stopper (ACGS), the current state-of-the-art linear gas stopper, through the use of new simulations and beam purification via collision-induced-dissociation. The methods include: 1. Development of a prototype low-energy, low-pressure CID gas stopper. This prototype features a thin, approximately 20 nm, Si 3 N 4 entrance window designed for dissociating stable and rare-isotope molecular ions. The goal is to enhance beam purification and overall efficiency. 2. Creation of Particle-in-cell (PIC) simulations to assess the advantages of multi-layer RF carpets and multi-point extraction for ion transport efficiency. These simulations rely on the 3DCylPIC package, specifically designed for studying devices of this nature. The goal is to quantify and mitigate ion transport losses due to space charge generated in the stopping process of large numbers of ions. 3. Perform ion transport simulations across an RF carpet using an 8-phase travelling wave and evaluate its performance. Compared to the 4-phase RF carpets used in ACGS, the 8-phase carpets will double the wavelength of the generated traveling wave allowing for larger maximum RF amplitudes that could result in improved ion transport efficiency for high-intensity incoming beams when large space charge fields are present. Impact: Tackling the primary challenges associated with transforming high-energy projectile fragment beams into low-energy beams—specifically, addressing efficiency and purity—holds significant promise for advancing FRIB science. This advancement will particularly impact precision mass measurements, laser spectroscopy of short-lived nuclei, and studies in astrophysics and nuclear reactions using reaccelerated beams. These domains play a crucial role in addressing key questions outlined in the 2023 NSAC long-range plan, spanning nuclear structure, nuclear astrophysics, and fundamental symmetries. Additionally, they contribute to addressing 10 out of the 17 benchmarks identified by the NRC RISAC. The development of a next-generation gas stopping device capable of delivering low-energy, rare-isotope beams at a rate of 10 8 particles per second, or more, with high purity holds the potential to unlock experiments that would otherwise be unfeasible. Furthermore, it is expected to reduce the time required for experiments at FRIB, thereby maximizing scientific output. The research and development activities performed as part of this project bolstered essential competencies at FRIB in beam physics and ion source technologies, provided valuable training opportunities for junior scientists.

43 PARTICLE ACCELERATORS↗

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

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

01 COAL, LIGNITE, AND PEAT↗

A New Method for Low Cost Production of Titanium Alloys for Reducing Energy Consumption of Mechanical Systems

This project investigated an innovative manufacturing process intended to minimize the cost of production of titanium materials and components, and increase the adoption of Ti components for energy consuming applications, such as automobiles. A key innovation of the proposed manufacturing approach is a novel Ti powder sintering technology for making titanium materials with ultrafine grain microstructure in the as-sintered state with minimum, or an absence, of post-sintering processes. The new sintering technology is termed Hydrogen Sintering and Phase Transformations (HSPT), and constitutes a promising manufacturing technology that can be used to produce titanium (Ti) materials and components in a near-net-shape form, thus also minimizing machining costs. Our objective was to meet, or possibly surpass, the mechanical property levels for ASTM B348 Grade 5 for wrought Ti-6Al-4V. Although specific applications call for varying mechanical property requirements, ASTM B348 was created for the demanding applications of the aerospace industry, and is the established standard for Ti-6Al-4V. While the primary goal was to meet, or exceed this standard, the team also had the goal of demonstrating this could be done at a significantly lower cost of production. Interim goals of the project were to fully develop this novel sintering process, and provide sufficient baseline testing to make the method practical and attractive to industry. By optimizing the process parameters for the sintering of titanium hydride (TiH 2 ) powders in a hydrogen atmosphere and controlling the phase transformations during and after sintering, the HSPT process was expected to reduce the energy consumption, and thus cost, of making Ti alloys and fabricating Ti components. The process was designed such that no high temperature melting is required for producing Ti alloys; little or no post-sintering processing is needed for producing desired microstructures (and therefore enhanced mechanical properties), and finally, minimum machining is needed to fabricate finished Ti components. An energy analysis within this report provides more detail, but calculated values indicate that the HSPT process is less than half as energy intensive as conventional wrought processing, while producing mechanical properties that are comparable. In addition to the energy savings anticipated from the industrial production of Ti components, a second prong of energy savings resides in the use phase of components produced, primarily from use in the transportation sector. Titanium has a number of material qualities appropriate for the auto industry, particularly low mass and corrosion resistance. By reducing the weight of automobiles and other vehicles, energy costs and CO 2 production will be reduced over the lifetime of the vehicles, and components in corrosive environments on vehicles, such as exhaust systems and other under carriage parts, may not have to be replaced during a vehicle’s lifetime. Our analysis indicates that by replacing only 5.6 kg of steel parts in an auto with Ti components across the entire US fleet would save approximately 486 million gallons of gasoline per year. This correlates to a reduction of 3.6 million metric tons of CO 2 per year. The potential for replacing many more of the steel parts in automobiles with lighter weight titanium components is clear. The project was very successful overall, meeting all milestones and surpassing project goals in terms of mechanical properties and microstructures produced. In addition to tensile properties, fatigue properties were emphasized in the project work. Powder metallurgy processes often have porosity to some degree in their final microstructure, and porosity is a well-known cause of crack initiation and low fatigue performance. Although many automobile applications do not undergo fatigue stress regimes, many others do encounter cyclic stress, and design criteria in the latter case require good fatigue properties. Production and testing of HSPT parts showed excellent tensile properties and fracture toughness, and fatigue properties that exceeded all previously reported powder metallurgy Ti methods, overlapping with wrought processed values. Fatigue limits exceeded 500 MPa and tensile strength exceeded 1,000 MPa while maintaining good ductility. Microstructures produced during the project period easily surpassed pre-project expectations. In addition to producing very fine grains in the as-sintered state (without post sintered thermo-mechanical work), porosity was reduced and industrially relevant microstructures previously undemonstrated in any other powder metallurgy titanium method were produced using HSPT materials. These microstructures, both bi-modal and globularized, were produced with simple post-sinter heat treatments, but without the need for energy intensive mechanical work. The employed heat treatments expanded the available mechanical property range (tensile strength vs. ductility) of the HSPT system in Ti-6Al-4V. The project has resulted in the publication, thus far, of five refereed journal articles and five conference proceedings papers, as well as a patent application, two dissertations and a master’s thesis. Two additional journal articles are currently under review, and at least three others are currently in preparation, with several additional students anticipated to graduate within the coming year. Presentations and papers were a particular focus of the second half of the project, once significant experimentation had been performed and analyzed. As part of our efforts to disseminate information of our results, the Ti research teams within Prof. Fang’s and Prof. Chandran’s research groups had a strong presence at the 13th World Conference on Ti, August 16-20, 2015, in San Diego. Several research groups in the US and in Europe are now performing experiments using the HSPT process. Accompanying efforts to bring HSPT to the Ti community at large, and industry in particular, work has continued with our partners and with other interested industrial Ti users and producers, including Boeing and GKN (a major powder metallurgy parts manufacturer). Commercialization has been a central focus of the final phase of the project, and Reading Alloys signed a provisional licensing agreement in summer of 2015. They are currently seeking an appropriate customer with which to pursue initial parts manufacturing efforts. Other licensing options and partners are continuing to be pursued. The promise of lightweight, strong and corrosion resistant Ti alloys with long fatigue lifetimes for automobile or transportation applications has been the vision of the metal industry since titanium came to the attention of scientists and engineers. The sole limitation of realizing these goals has been cost, which is primarily a function of energy used in production. The HSPT process was shown through this work to be capable of realizing this goal, and facilitating the practical use of titanium in US automotive and other industries.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Testing- and Model- Based Optimization of Coal-fired Primary Heater Design for Indirect Supercritical CO 2 Power Cycles (Final Scientific and Technical Report)

The overall objective of this project was to perform the R&D necessary to mitigate the risk associated with the design of a primary heat exchanger for a solid-fired combustion system coupled with an indirect-fired closed-loop Brayton Cycle utilizing supercritical CO 2 . The key technological hurdle was the coupling of a solid-fuel firing system with the primary heater, which poses a singular challenge, which is the management of burner performance and operational conditions in a way to manage heat exchanger tube metal temperatures and temperature ramp rates in the absence of fluid phase change on the inside of the tubes. We designed and built the first ever pseudo power system employing a simple recuperated supercritical CO 2 closed-loop Brayton Cycle coupled to a solid-fuel fired system. Advanced coupled CFD and process modeling were used to design the primary heat exchanger (PHX), which consisted of both radiative and convective sections, to limit tube metal temperatures resulting from the heat release profile of the solid fuel flame near the radiative tubes. The heat exchanger was designed to produce finished CO 2 temperatures of 600 °C a pressure of 20.7 MPa and CO 2 flow of 5.5 kg/s. The constructed PHX was capable of 1.2 MWth heat uptake. During design of the PHX, the modeling showed that most variables influencing flame shape (burner stoichiometric ratio and register velocities and swirl) were not suitable to manage heat flux to the metal surfaces. This is because they substantially increased adiabatic flame temperature through the influence of localized stoichiometric ratio. Excess air and firing rate were the two most powerful variables that could be used to control tube surface temperatures. The coupled system was operated for a total of 407 hours, with the longest continuous run of 248 hours. For 62% of the operational time, the unit was unmanned and in automatic control. The fuels used for the testing included natural gas, two Utah Bituminous coals, woody biomass, and bagasse. During the testing we were able to verify the 1.2 MWth heat uptake and we operated at a finished CO 2 temperature of 607 °C and a pressure of 20.3 MPa simultaneously. The real-time corrosion rate of the Super 304H tube CO 2 surface in the region of the radiative section of the PHX were measured, at an approximate temperature of 550 °C. The two key variables related to corrosion rate are the pressure and flow rate of the CO 2 . A technoeconomic analysis was performed at a scale of 120 MWE. The updated analysis showed that the efficiency of an sCO 2 power producing plant will be related to the pressure drop of the PHX.

01 COAL, LIGNITE, AND PEAT↗

BWR Spent Nuclear Fuel Acquisition and Testing to Support DOE-NE High Burnup Spent Fuel Data Project

The Office of Spent Fuel and Waste Disposition (SFWD) within the US Department of Energy (DOE) Office of Nuclear Energy (NE) established the Spent Fuel and Waste Science and Technology (SFWST) campaign to conduct research and development (R&D) activities related to the storage, transportation, and disposal of spent nuclear fuel (SNF) and high-level radioactive waste. The SFWST program was created within SFWD to address issues of extended or long-term SNF storage and transportation. Some near-term objectives of SFWST are to use a science-based, engineering-driven approach to: Support the enhancement of the technical bases to support the continued safe and secure dry storage of SNF for extended periods; Support the enhancement of the technical bases for retrieving SNF after extended dry storage; Support the enhancement of the technical bases for transporting high burnup (HBU) fuel and transporting low burnup fuel and HBU fuel after dry storage DOE-NE, in partnership with the Electric Power Research Institute, developed the High Burnup Spent Fuel Data Project to perform a large-scale demonstration and laboratory-scale testing of HBU pressurized water reactor (PWR) fuels (exceeding 45 gigawatt-days per metric ton of uranium [GWd/MTU]). Under this project, 25 sister rods—which are rods that have the same design, power histories, and other characteristics—were removed from assemblies at the North Anna Nuclear Power Station and sent to Oak Ridge National Laboratory (ORNL) in January 2016. ORNL performed detailed nondestructive examination (NDE) on all 25 rods. The NDE consisted of visual examinations, gamma and neutron scanning, profilometry and rod length measurements, and eddy current examinations. After completing the NDE, 10 of the sister rods were delivered to Pacific Northwest National Laboratory (PNNL) in a NAC International, Inc. legal-weight truck cask in September 2018 for destructive examination (DE). To date, SFWD work has focused on the PWR fuel that is part of the Sister Rod Test program. No boiling water reactor (BWR) fuel has been tested in the program, and the data needs that were identified for the PWR fuel have not been collected for BWR fuel. The goal to obtain six to nine BWR rods and test them at ORNL will support closing this important data gap. BWR fuel comprises approximately 56% of the total fuel assemblies currently in storage at nuclear power plants in the United States. BWR nuclear fuel and cladding designs and manufacturing are significantly different from PWRs. Differences include the following: BWR fuel pellets are larger than PWR pellets; Variations of Zircaloy-2 (including liners) are used instead of the Zircaloy-4 cladding materials used in PWRs; Clad manufacturing and stress-relief processes are different between PWRs and BWRs; The fuel rod dimensions are different because larger rod diameters and thicker cladding are used in BWRs; BWR fuel typically has lower internal rod pressures and sees vastly different operating conditions than PWR fuel (i.e., two-phase flow); BWR assemblies are “canned,” meaning each assembly is surrounded by a metal fuel channel; BWR cladding is often composed of an inner pure Zr liner that has widely different mechanical properties than the Zircaloy-2 alloy and exhibits a stronger affinity for hydrogen; The BWR SNF generally has more total hydrogen in the cladding/liner than typical PWR fuel; The construction of the PWR and BWR assemblies is vastly different; BWR rods are solidly attached to the assembly nozzles and experience a much different vibration and shock load than PWR rods, which are “floating” within a grid system attached to guide tubes, and the rods sit loosely on the bottom end plates. These numerous differences will affect the way the BWR SNF responds under dry storage preparation processes (e.g., vacuum drying) and during transportation. The results collected in the PWR experimental program must be compared with a subset of similar data collected on BWR SNF to establish a technical basis for whether the larger PWR database is sufficient to bound the BWR SNF end-of-life conditions as is currently assumed for several fuel/clad properties. Changes that occur in both fuel types at HBU could exacerbate any mechanical property differences. As the fuel burnup increases, several changes occur that might affect the performance of the fuel, cladding, and assembly hardware in storage and transportation. These changes include increased cladding corrosion layer thickness, increased cladding hydrogen content, increased cladding creep strains, increased fission gas release, and the formation of the HBU structure at the surface of the fuel pellets. The Nuclear Regulatory Commission (NRC) limits the current maximum rod-averaged burnup to 62 GWd/MTU due to these changes and the lack of data at higher burnups.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Performance Assessment for Closure of the Calcined Solids Storage Facility at the INL Site - 19289

The Calcined Solids Storage Facility is composed of six discrete stainless steel 'bin sets' at the Idaho Nuclear Technology and Engineering Center at the INL Site. Approximately 4,440 m{sup 3} of the dry, granular calcine waste is stored in the facility. Closure of the facility will entail removal of the calcine waste followed by grouting. A 5-cm thick residual amount of calcine was assumed to remain at the bottom of the bins. Important radionuclides in the residual layer for the groundwater pathway were identified as the fission and activation products Tc-99, Se-79, I-129, Cs-135, H-3, and actinides Pu-239, and Pu-240. The stainless-steel bins prevent near-term release of grouted residual calcine waste to the environment. Future releases were conservatively based on the corrosion rate of the weld seams on the stainless-steel bins, which site-specific studies measured mean tunnel corrosion rates of 0.0132 mm/yr with a coefficient of variation of 0.333. Based on tunnel corrosion rates, mean bin hydrologic containment times were estimated to be between 238 years and 1,191 years, depending on the thickness of the bin. Any breach in the weld seam was assumed to result in water flow and radionuclide leaching through the bin. The structural integrity of the bin was assumed to be controlled by the general stainless-steel corrosion rate, which was much slower (3.5 x 10{sup -5} mm/yr), resulting in structural integrity times between 24,000 years to over 100,000 years. A mixing-cell model was used to compute water flow and transport through the hydrologically failed bins and the vadose zone. Credit was taken for the geochemical properties of the grout used to fill the bins at closure, which significantly slows the releases from the bins for key radionuclides like Tc-99. A semi-analytical model for three-dimensional dissolved-phase transport in an aquifer of constant unidirectional flow, infinite length and width, and finite thickness was used to compute concentrations in the aquifer across the model domain for a time varying source. Cumulative impacts from all bin sets were calculated by summation across the bin set sources that were superimposed on a uniform Cartesian grid. Effective doses were calculated for direct groundwater ingestion from wells located 100-m downgradient from each of the bin sets and ingestion of meat, milk, and vegetables irrigated with water derived from the well. Intruder doses were based on a well drilled through the grouted waste and cuttings spread around a home site after a 500-year institutional control period. Peak doses during the 1,000-year compliance period were less than 5% of the 250 μSv/yr (25 mrem/yr) performance objective and intruder doses were generally less 5% of the of the 5000 μSv/yr (500 mrem) and 1000 μSv/yr (100 mrem) performance objective for acute and chronic intruder exposure, respectively. Groundwater concentrations were all less than the maximum contaminant limits. In all cases, doses were highest for the oldest bin set because, compared to the other bin sets, it had the highest radionuclide concentrations in the calcine waste. The performance assessment results indicate that performance objectives can be met leaving considerably more than 5-cm of residual calcine waste in the bin sets. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modeling Mesoscale Processes of Scalable Synthesis (Final Report)

The main objective of the center is to develop tools that can help for the understanding of equilibrium and non-equilibrium properties of nanoparticles in aqueous electrolytes. To achieve this goal, we need to develop new computational methodologies that can integrate molecular level description such as the electronic structure analysis to mesoscale description such as generalized hydrodynamics. The challenge is to develop systematic coarse graining techniques that can adaptively choose the level of description needed for the particular task. For example, if bonding breaking is involved, then one should use appropriate models from quantum mechanics to accurately model the time scale and molecular mechanism for the breaking process. However, to model transport process or the process of self-assembly, much less detailed models can be used. There are two main components to this task. One is to develop accurate models that bridge different levels of description. The classical example is the QM/MM model that bridges quantum mechanics and molecular mechanics for studying chemical reactions of large molecules. This work was awarded the Nobel prize this year. For our purpose, we also need to develop methods that bridge molecular and continuum levels of description. Another main task is to develop systematic coarse-graining methods that can give rise to accurate mesoscale and macroscale models, starting from more detailed models at the atomic or molecular scale. None of these ideas are new. By now they have been fashionable topics for quite some time. But neither problem is satisfactorily resolved. The main problem is that uncontrolled approximations are used. As a result, we have very little idea about the validity and accuracy of these models. The main purpose of the center is to exploit the Mori-Zwanzig formalism to achieve the goals stated above. The Mori-Zwanzig formalism is a general tool for eliminating degrees of freedom from existing models. What is attractive about this formalism is its generality. The flip side is that it is also somewhat of a tautology. In order to make it a practical tool, one has to make approximations and this is where the trouble comes. Instead of focusing on Mori-Zwanzig formalism, the PI proposes to study the closely related renormalization group formalism (RG). Like Mori-Zwanzig, RG is a general strategy for eliminating small scale degrees of freedom in a model. Unlike Mori-Zwanzig, RG has the advantage that at least for some problems, one can write a PDE (partial differential equation) type of model for the exact RG flow, without any approximation. This gives us a solid starting point for discussing and analyzing approximations. The PI studied approximate RG flows for turbulent diffusion and related problems. However, extending it to more practical problems proved to be hard. The second line that the PI has pursued is to develop special numerical algorithms for studying particularly important events for dynamics processes in meso-scale systems. The PI succeeded in doing so for the problem of understanding phase transition in metallic systems.

97 MATHEMATICS AND COMPUTING↗

Integrated Commercial Carbon Capture and Storage (CCS) Prefeasibility Study at Rock Springs Uplift, Wyoming (Final Report)

To establish the potential of safe, long-term commercial-scale carbon capture & storage (CCS) in southwestern Wyoming, the University of Wyoming formed a CCS coordination team (CCT) to complete a Phase I CarbonSAFE pre-feasibility study adjacent to the region’s largest coal-fired power plant, the Jim Bridger Plant (JBP). The study’s primary objective was to evaluate the potential of storing 50+ million metric tons of anthropogenic carbon dioxide (CO 2 ) in stacked reservoirs over a period of 25 years. In addition, this integrated study evaluated related aspects of carbon capture utilization & storage (CCUS) at the study site, including: (1) the CO 2 character of post-combustion flue gas from JBP and the technical capacity for retrofitting proven commercial-scale CO 2 capture technology relative to CO 2 source character and technological constraints; (2) utilization of the existing CO 2 pipeline network in the immediate vicinity of JBP for enhanced oil recovery (CO 2 -EOR) and utilizing that network as part of a regional hub for CO 2 storage from other anthropogenic sources; (3) assessment of the challenges and benefits of meeting all of CarbonSAFE programmatic Phase I goals relative to Wyoming’s carbon management regulatory framework; and (4) assessment of the community and environmental factors that may impact CCS at JBP. With respect to storage, we analyzed the potential of stacked Mesozoic reservoirs at the study site, focusing on the Entrada and Nugget sandstones. The reservoir response to injection simulations varied greatly between formations, though they both were capable of CO 2 storage. The Entrada Sandstone is relative thin (~40’) and heterogeneous; storage using one injection and production well was limited to 1.2 Mt over 25 years at the study site. The Nugget Sandstone is much thicker (> 400’) and generally homogeneous; storage using one injection and production well approached 15.0 Mt over 25 years; this scenario suggests that the storage potential in the Nugget Sandstone is close to the estimated P90 calculation of 9.6 Mt/mi 2 as models indicate a storage potential of 8.3 Mt/mi 2 at saturation. Both reservoirs’ storage capacity benefited greatly from coupling an injection well with an offset production well; total storage capacity increased and reservoir pressure stayed below critical thresholds. This indicates that implementing a pressure management strategy at the study site would optimize storage and decrease risk. Seals associated with the targeted reservoirs were shown to be capable of retaining the potential storage volumes, and the site benefits from multiple seal redundancies (i.e. > 7,000’ of seal formations). Furthermore, assessment of different reservoir variables (e.g., fluid composition, pressure, etc.) indicates all reservoirs at the study site are isolated. Storage assessments from Mesozoic reservoirs were coupled with simulations from the site’s Paleozoic reservoirs (the Weber Sandstone and Madison Limestone), which were the focus of previous study. The fully integrated stacked reservoir pre-feasibility studies suggest that the Nugget and Madison formations alone could safely store 50 million metric tons of CO 2 within the study site. Furthermore, both formations are true saline reservoirs near JBP, making them ideal target reservoirs for CO 2 sequestration. Analysis of Wyoming’s laws and regulations confirm that the State has a favorable policy environment for CCS/CCUS. The project’s economic models estimated that the following revenues collectively are sufficient to finance the project over its lifespan if CO 2 utilization is coupled with CCS via the following strategies: (1) sales of CO 2 for EOR (approximately $69 million/year); (2) use of CO 2 tax credits such as amended §45Q and §48A (approximately $484 million total); and (3) sales of low-carbon electricity and marketable carbon offset/credits into carbon-constrained markets (approximately $11-$17 million/year) to West Coast consumers. These revenue estimates are broadly consistent with other studies that concluded arbitrage in JBP retail electricity sales between California and Wyoming could help to support the cost of deploying CCS at JBP. The models also indicate that maintaining a successful business case for long-term CCS -- wherein associated capture, transportation and operational costs are fully addressed at JBP -- is reliant on numerous external factors that complicate predictions. This study suggests that this site meets the CarbonSAFE program’s requirements of being able to feasibly store 50+ million metric tons of CO 2 over 25 years within the site’s stacked reservoirs, especially if coupled with pressure management strategies. This study suggests that the Nugget Sandstone and the Madison Limestone are two of the best reservoirs in the State with respect to overall storage capacity. In addition, the site’s proximity to existing CO 2 transportation networks, CO 2 -EOR opportunities, ability to sell to markets that value low-carbon electricity, Wyoming’s existing carbon regulatory framework, and a public that is well-educated with respect to energy markets and issues collectively increase the potential for implementing commercial-scale CCS adjacent to the JBP. With respect to the available technical and non-technical site data, and the reservoir quality at the study site, we suggest that this site is capable of moving directly into Phase III of the CarbonSAFE program.

01 COAL, LIGNITE, AND PEAT↗

Air-coupled tsunamis generated from impacts and airbursts: Our understanding before Hunga-Tonga Hunga-Ha'apai

The effort to prevent or mitigate the effects of an impact on Earth is known as planetary defense. A significant component of planetary defense research involves risk assessment. Much of our understanding of the risk from near-Earth objects comes from the geologic record in the form of impact craters, but not all asteroid impacts are crater-forming events. Small asteroids explode before reaching the surface, generating an airburst, and most impacts into the ocean do not penetrate the water to form a crater in the sea floor. The risk from these non-crater-forming ocean impacts and airbursts is difficult to quantify and represents a significant uncertainty in our assessment of the overall threat. We are currently working to better understand impact scenarios that can generate dangerous tsunamis. One of the suggested mechanisms for the production of asteroid–generated tsunamis is by direct coupling of the pressure wave to the water, analogous to the means by which a moving weather front can generate a meteotsunami. To test this hypothesis, we ran a series of airburst simulations and provided time-resolved pressure and wind profiles to use as source functions for tsunami models. We used the CTH hydrocode to model the various airburst scenarios to compare to the results of other simulations and provide time dependent boundary conditions as input to shallow-water wave propagation codes. The strongest and most destructive meteotsunamis are generated by atmospheric pressure oscillations with amplitudes of only a few hPa1 (mbar), corresponding to changes in sea level of a few cm. The resulting wave is strongest when there is a resonance between the ocean and the atmospheric forcing. A Proudman resonance takes place when the atmospheric disturbance’s translational speed (U) equals the longwave phase speed $\sqrt{gh}$ of shallow water wave. Coupling is strongest when the Froude number (Fr=U/c) is unity. A weather front propagates much slower than the speed of sound, so meteotsunamis are most common and dangerous in shallow bodies of water such as the Mediterranean Sea or Lake Michigan. By contrast, the blast wave from an airburst or crater-forming impact propagates at a speed faster than a tsunami in the deepest ocean, and a Proudman resonance cannot be achieved even though the overpressures are orders of magnitude greater. However, blast wave profiles are N-waves in which a sharp shock wave leading to overpressure is followed by a more gradual rarefaction to a much longer-duration underpressure phase. Even though the blast outruns the water wave it is forcing, the tsunami should continue to be driven by the out-of-resonance gradient associated with the suction phase, which may depend strongly on the details of the airburst or impact scenario. The open question is whether there are any conditions under which such an airburst-driven tsunami can be dangerous enough to contribute to the overall impact risk. We have also identified other potential mechanisms for airburst-generated tsunamis: 1) reaction force at the surface from the plume ejected into space, which carries significant momentum, 2) expanding toroidal vortices at the surface, which travel more slowly than the shock wave and can generate a Proudman resonance in relatively shallow ocean (such as continental shelf), and 3) steam explosion from seawater ablation by a “Type II” (Libyan Desert Glass-type) airburst in which the hot vapor jet descends to the surface. On January 15, 2022, the Hunga-Tonga Hunga-Ha’apai volcano, located approximately 60 km north of Tongatapu, the main island of Tonga, violently erupted with a powerful explosion, culminating the period of volcanic activity that started in December of 2021. This event and resulting tsunamis provided an existence proof for the air pressure wave coupling mechanism we proposed. It also suggests that it can be stronger and more significant over much greater distances than we contemplated, leading to global tsunamis associated with impact events on land as well as in the water. Large atmospheric explosions generate global Lamb waves with larger amplitudes, longer periods, and slower speeds than the local and regional blast waves we modeled prior to that event. This paper reviews our analysis and modeling of airburst-driven tsunamis prior to the 2022 Hunga-Tonga Hunga-Ha’apai tsunami, which was the subject of two presentations at the 2023 Planetary Defense Conference and is the subject of another paper currently in preparation.

54 ENVIRONMENTAL SCIENCES↗