Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure”

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 829 records · Page 46

Monitoring Potential Transport of Radioactive Contaminants in Shallow Ephemeral Channels: FY2019

Desert Research Institute (DRI) conducted a field assessment of the potential for contaminated soil to be transported from the Smoky Site Contamination Area (CA) because of storm runoff. This activity supported U.S. Department of Energy (DOE) Environmental Management Nevada Program (EM NV) efforts to establish post-closure monitoring plans for the Smoky Site Soils Corrective Action Unit (CAU) 550. The work was intended to confirm the likely mechanism of transport and determine the meteorological conditions that might cause the movement of contaminated soils, as well as determine the particle size fraction most closely associated with transported radionuclide-contaminated soils. These data will facilitate the design of the appropriate post-closure monitoring program. In 2011, DRI installed a meteorological monitoring station on the west side of the Smoky Site CA and a hydrologic (runoff) monitoring station within the CA, near the east side. The meteorological station collected air temperature, wind speed, wind direction, relative humidity, precipitation, solar radiation, barometric pressure, soil temperature, and soil water content data. The maximum, minimum, and average or total values (as appropriate) for each of these parameters were recorded for each 10-minute interval. The maximum, minimum, and average water depth in the flume installed at the hydrologic station were also recorded for each 10-minute interval. This report presents the data collected from these stations during fiscal year (FY) 2019. During the FY2019 reporting period, the warmest months were June, July, and August and the coldest were December and February. The highest monthly solar radiation (i.e., sunshine) values were recorded in June, July, and August, whereas December through February had the lowest solar radiation values. Monthly mean wind speeds were highest in the spring (April) and early summer (June). Winds were predominantly from the north, except in February, July, August, and September when winds were mostly from the west. The monthly average relative humidity ranged from 17 percent to approximately 60 percent. Humidity was lowest in the summer and highest in the winter. Monthly total precipitation ranged from 0.0 during August and September to 2.22 inches (in) (56.39 millimeters [mm]) in February. Total precipitation for FY2019 was 7.83 in (198.88 mm). During the reporting period, a single, nonzero, flow event was recorded at the flume. This occurred over two days between March 6 and 7, 2019, in response to 1.16 in (29.46 mm) of precipitation falling within 9.3 hours on March 6, 2019. This event produced a peak water depth of 0.39 in (0.99 centimeters [cm]) in the flume, and runoff was measured for approximately 23.5 hours. However, the maximum recorded flow depth was less than the limit for flow discharge and velocity that can be calculated for the six-inch (15 cm) throat of a Parshall flume. The runoff event produced insufficient flow velocities through the natural channel to cause local erosion and transport of bedload materials. No bedload samples were collected for radiological analysis. Observing meteorological and environmental conditions that lead to storm runoff will help identify the parameters and threshold conditions that should be incorporated into post-closure monitoring plans for this and similar Soils CAU sites. Routine monitoring of meteorological and hydrologic parameters at the Smoky Site CA has been discontinued at the request of EM NV, effective at the conclusion of FY2019. Environmental data acquisition was completed and the meteorological and runoff monitoring stations were removed from the Smoky Site on October 15, 2019.

54 ENVIRONMENTAL SCIENCES↗

Features Events and Processes Relevant to DPC Disposal Criticality Analysis

The Department of Energy is evaluating the technical feasibility of disposal of spent nuclear fuel in dual-purpose canisters in various geologies. As part of ongoing research and development, the effect of potential post-closure criticality events on repository performance is being studied. Many different features, events, and processes (FEPs) could affect the potential for criticality or the extent of a criticality event. Additionally, a criticality event could affect other FEPs. This report uses existing lists of FEPs as a starting point to evaluate the FEPs that could affect or be affected by an in- package criticality event. The evaluation indicates that most of the FEPs associated with the waste form, the waste, or the EBS have some effect on post-closure criticality and/or are affected by the consequences of post-closure criticality. In addition, FEPs not previously considered are identified for further development.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Summary Report for the NEPA Impact Analysis. Revision 1

The Waste Isolation Pilot Plant (WIPP), located in southeastern New Mexico, has been developed by the U.S. Department of Energy (DOE) for the geologic (deep underground) disposal of defense-related transuranic (TRU) waste. Containment of TRU waste at the WIPP facility is derived from standards set forth in Title 40 of the Code of Federal Regulations (CFR), Part 191. The DOE assesses compliance with the containment standards according to the Certification Criteria in Title 40 CFR Part 194 by means of Performance Assessment (PA) calculations performed by Sandia National Laboratories (SNL). WIPP PA calculations estimate the probability of radionuclide releases from the repository to the accessible environment for a regulatory period of 10,000 years after facility closure. The DOE Carlsbad Field Office (CBFO) has initiated a National Environmental Policy Act (NEPA) action for a proposal to excavate and use additional transuranic (TRU) waste disposal panels at the WIPP facility. This report documents an analysis undertaken as part of an effort to evaluate the potential environmental consequences of the proposed action. Although not explicitly required for a NEPA analysis, evaluations of a dose indicator to hypothetical members of the public after final facility closure are presented in this report. The analysis is carried out in two stages: first, Performance Assessment (PA) calculations quantify the potential releases to the accessible environment over a 10,000-year post-closure period. Second, dose was evaluated for three hypothetical exposure pathways using the conservative radionuclide concentrations assumed to be released to the accessible environment.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

PORFLOW modeling of Vadose zone flow and transport for the E-area intermediate level vault

In support of the E-Area Performance Assessment, a two-dimensional model of water flow and radionuclide transport through the E-Area Intermediate Level Vault (ILV) and local vadose zone has been developed using the PORFLOWTM software. The purpose of the model is to calculate flux to the water table for radionuclides eluted from the ILV during its operational life, the period of institutional control, and times following site closure. Results of model calculations will be used by a three-dimensional PORFLOW model of transport through the aquifer to determine radionuclide concentrations at a hypothetical 100 meter well and at the site boundary where contaminated groundwater could be assessed by members of the public following site closure. While newly developed, the model structure closely follows that used in the 2008 PA while incorporating a refined computational mesh, updated material properties for the vadose zone soil and vault concrete, revised infiltration rates, and new solid-liquid distribution coefficients. The model also addresses degradation of concrete hydraulic properties by blending soil and concrete water retention curves over a 500 year period. That is, the hydraulic properties of ILV concrete start out as fresh concrete at the time of site closure and degrade to soil properties over a 500 year time period. This approach has not been used previously in E-Area PA’s. The model has been used to calculate water flow through the ILV vadose zone and, on a trial basis, transport of a limited number of radionuclides from the waste region inside the ILV to the water table. Figure ES-1 provides a view of the model and an example of water flow after the concrete has substantially degraded. Figure dimensions are in centimeters. The red area in the left hand side figure is the waste disposal region within the ILV. Colored bands in the figures indicate different soil and concrete regions. Figure ES-2 shows a calculation of flux to the water table (mol/yr) for U-234 and its one-year half-life decay chain (U-234 → Th-230 → Ra-226 → Pb-210). Pb-210 flux is too small to show on the same scale used for Th-230 and Ra-226 in Figure ES-2. The model has also been used to make a preliminary evaluation of the release of tritium from TPBAR disposal containers that will be place in the ILV. A separately calculated tritium release rate as a function of time was used as a source in the PORFLOW model. The TPBAR disposal schedule and tritium release rate differ from the values used in the 2008 PA. Results from these initial trial runs are compared to calculations made for the 2008 PA ILV analysis. In brief, a new PORFLOW model of flow and transport through the ILV vadose zone has been developed and tested. The model is ready to be used for PA calculations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Corrosion of Steel During Long-Term Exposure to Evolving Cementitious Environments

SRS is proceeding with closure of the H-Area Tank Farm (HTF) and F-Area Tank Farm (FTF). Closure consists of removing the bulk waste, heel removal, and filling the tank with tailored grout formulations. The long-term performance of the concrete and steel materials of construction and the closure grout is an important consideration for Performance Assessment (PA) of closed tanks in FTF and HTF. Over PA timeframes of hundreds, thousands, to tens of thousands of years, the chemical and physical properties of the concrete and steel materials will slowly degrade due to environmental exposure and material aging. In the interim these materials will provide a barrier to the leaching of radionuclides into the soil. Savannah River Remediation (SRR) is updating relevant portions of these analysis inputs to support an imminent revision to the HTF PA and a future update to the FTF PA. This analysis provides an update to the previous PA inputs for steel corrosion. The Central Scenario for the PA analysis includes three postulated modeling cases: 1) Realistic Case, 2) Compliance Case, and 3) Pessimistic Case. This analysis also included a Fast Flow Path Case. This latter case considered the circumstance where the initial condition of the concrete was a completely degraded state and that grout shrinkage exposed the steel to the soil environment immediately. In effect, the steel was unprotected by the concrete and grout. The cases have various degrees of conservatism considered. Chemical, physical and tank configuration parameters were investigated to understand their effects on the predicted time to release of the contaminants. The assessment reviewed the initial tank and steel configuration, service life degradation of the steel, and potential corrosion mechanisms associated with degradation of the concrete materials.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Performance Assessment for the E-Area Low Level Radioactive Waste Disposal Facility at the Savannah River Site: Chapter 1

This report documents the revised Performance Assessment (PA) analysis for the E-Area Low-Level Waste Facility (ELLWF) at the United States (U.S.) Department of Energy (DOE) Savannah River Site (SRS). A PA analysis is required for DOE-operated facilities that dispose of low-level radioactive waste. PA analyses simulate (1) the release of radionuclides from the disposal site after facility closure, (2) transport of those contaminants through the environment, and (3) exposure/impacts to potential receptors. The purpose of the PA analysis is to demonstrate that the facility is operated in a manner that ensures long-term environmental protection after facility closure, thereby providing for the protection of public health and safety in limiting doses to a hypothetical member of the public (MOP) or an inadvertent human intruder (IHI). DOE Manual (M) 435.1-1, Chg. 3, Radioactive Waste Management (U.S. DOE, 2021b) establishes quantitative post-closure environmental impact limits and requires a facility-specific PA analysis to demonstrate compliance with these limits for DOE low-level waste (LLW) disposed of after September 26, 1988. These limits are defined in terms of human health (e.g., dose limits) with respect to radioactive constituents in the waste.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

2020 Postclosure Groundwater Monitoring and Inspection Report, Central Nevada Test Area, Subsurface Corrective Action Unit 443

This report presents the groundwater monitoring data collected by the U.S. Department of Energy (DOE) Office of Legacy Management (LM) from the Central Nevada Test Area (CNTA), Nevada, Site, Subsurface Corrective Action Unit (CAU) 443 in Nye County, Nevada (Figure 1). The CNTA is the site of an underground nuclear test in 1968 that resulted in residual contamination near the detonation depth of 3200 feet (ft); the contamination requires long-term monitoring. Responsibility for the environmental restoration and long-term monitoring was transferred from DOE’s National Nuclear Security Administration, Nevada Field Office, to LM on October 1, 2006. The environmental restoration and site closure process were completed in 2015 in accordance with the amended 1996 Nevada Federal Facility Agreement and Consent Order (FFACO) (State of Nevada et al. 1996, as amended) and all applicable Nevada Division of Environmental Protection (NDEP) policies and regulations. The Closure Report, Central Nevada Test Area, Subsurface Corrective Action, Unit 443 (DOE 2018), also called the Closure Report, originally was completed in January 2016 and revised in October 2018; it describes LM’s plan for long-term postclosure monitoring. This includes monitoring of the radioisotopes of interest and water elevations, inspecting the site and maintaining the institutional controls (ICs), evaluating and reporting data, and documenting the site’s records and data management processes (DOE 2018).

54 ENVIRONMENTAL SCIENCES↗

UNF-ST&DARDS Enhancements for RCCA Data in As-loaded Dual Purpose Cask Models

This report summarizes the work performed to enable detailed modeling of rod cluster control assembly (RCCAs) in dual purpose canisters (DPCs) in as-loaded configurations using the Used Nuclear Fuel – Storage, Transportation & Disposal Analysis Resource and Data System (UNF-ST&DARDS). The goal of this project was to evaluate the reactivity impact that RCCAs have on k eff of DPCs with pressurized water reactor (PWR) fuel to potentially use the additional margin in future post-closure criticality safety analysis. This preliminary evaluation determined the number of DPCs that currently require compensatory actions prior to emplacement in a repository because of their high reactivity under post-closure criticality scenarios, which could be made acceptable by including the as-loaded RCCAs as specified in the Unified Database (UDB). This report briefly describes the modeling methods currently used within UNF-ST&DARDS and the modifications made to automate inclusion of the as-loaded RCCAs in the DPCs for post-closure criticality calculations for the loss of neutron absorber (NA) scenarios, or NA models. For the Zion site, the loss of basket scenarios, or degraded basket (DB) models, are also included. Discussion is also provided regarding the UDB data compared to available site-specific loading map data for the sites, which are specifically evaluated herein. This report compares the modifications to a similar evaluation for the Zion by site conducted by Walker as an initial validation. After the partial validation, all existing PWR sites with applicable NA models within the UDB were evaluated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Using National-Scale Data to Inform Coal Power Plant Redevelopment and Coal Community Revitalization Planning

In the United States, individual coal plants have high variation in their construction and operation. Coal plants commonly have multiple units, with different commissioning ages, operational position, and maintenance conditions – even multiple fuels. These units may be in changing states of utilization, with variation in productivity depending on economic contexts. Coal plants themselves may be temporarily offline, mothballed, retired, or decommissioned due to owner/operator portfolios or market conditions, and the current posture of the plant may be difficult to ascertain without local news sources (Tarekegne et al. 2021). Anticipated dates of plant closures are prospectively represented to regulators and to the Energy Information Administration (EIA), but may also experience uncertainty due to the engineering, environmental, economic and social complexity of closing a very large power plant. As noted in this report, closure of the first unit to the last within a single plant may span more than 20 years. Yet tracking and predicting the position of our nation’s coal fleets on a national scale is deeply important. In developing a case for federal datasets and an approach to building and curating these data, this report identified that data must have high fidelity at the plant-level to be meaningful in aggregate. Plant-scale detail is also a match for its application: national priorities in federal investment range from economic revitalization for coal plant communities to maintaining reliable electric grid operations. This document is organized into a statistical review of coal plant retirements, past and prospective, in the United States; trends in the relationship between closures and existing federal investment programs for community revitalization; and potential methods for accelerating site redevelopment through advanced geospatial analysis based on federal data.

01 COAL, LIGNITE, AND PEAT↗

Streamlined Approach for Environmental Restoration (SAFER) Plan for Corrective Action Unit 578: Miscellaneous Inactive Sites Nevada National Security Site, Nevada, (Rev. 0)

This Streamlined Approach for Environmental Restoration (SAFER) Plan addresses the actions needed to achieve closure for Corrective Action Unit (CAU) 578, Miscellaneous Inactive Sites, identified in the Federal Facility Agreement and Consent Order (FFACO). CAU 578 comprises the following 11 corrective action sites (CASs) located in Areas 2, 3, 4, 6, 7, 12, and 23 of the Nevada National Security Site: • 02-26-01, Lead (Concrete Box w/Lining) • 02-26-07, Gun Turret w/Lead in Barrel • 03-16-01, U-3auS Disposal Site • 04-26-02, Lead on Instrument Bunker • 06-99-07, Sump • 07-99-01, Miscellaneous Debris • 12-02-11, UST-P-2 • 12-06-01, Muckpile • 12-06-04, Muckpile • 23-12-01, Boiler (2) • 23-99-04, Salvage Yard This plan provides the methodology for field activities needed to gather the necessary information for closing each CAS. There is sufficient information and process knowledge from historical documentation and investigations of similar sites regarding the expected nature and extent of potential contaminants to recommend closure of CAU 578 using the SAFER process. Additional information will be obtained by conducting a field investigation before selecting the appropriate corrective action for each CAS. The results of the field investigation will support a defensible recommendation that no further corrective action is necessary. This will be presented in a closure report that will be prepared and submitted to the Nevada Division of Environmental Protection (NDEP) for review and approval.

54 ENVIRONMENTAL SCIENCES↗

Final Report for Center for Tokamak Transient Simulations at USU

Providing plasma fluid codes like NIMROD with continuum drift kinetic (CDK) physics that is quantitatively valid and computationally feasible throughout the spatial domain is difficult. Work at Utah State University (USU), in collaboration with the Center for Tokamak Transient Simulations (CTTS), focused on applying CDK closures in disruption-related calculations. Three examples where kinetic physics is paramount are (1) the electron stress tensor closure in Ohms law for accurately describing neoclassical tearing mode (NTM) evolution, (2) runaway electron (RE) density (nRE) and current (jRE) moments in NIMROD’s extended MHD model for self-consistent evolution of RE populations during disruptions and, (3) energetic ion effects on a myriad of MHD instabilities. While NTM simulations and continuum and PIC approaches to energetic ions in NIMROD have been a major goals of USU’s closure work for several years, the development of self-consistent CDK RE capability in NIMROD was started and extended considerably during the CTTS effort. Some goals of CDK RE in NIMROD are to explore the effects of the 2D relativistic phase space in 4D simulations and compare with NIMROD’s fluid RE model. Four publications and two PhD theses came out of the USU CTTS effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Beneficial Use of Harvested Ponded Fly Ash and Landfilled FGD Materials for High-Volume Surface Mine Reclamation

The overall motivation of this project was to demonstrate at laboratory, bench-scale, and full-scale demonstration levels that (a) coal ash surface impoundments can go through closure by removal as per USEPA and state regulations so that the material can be used as is (other than draining free water using CCRs piles) in high-volume beneficial applications, (b) FGD material from closed out FGD facilities can be excavated and recompacted for coal mine reclamation, and (c) harvested CCRs can be beneficially utilized (providing a net environmental gain) in large-volumes for reclamation at abandoned coal mine sites across the US, especially in the Eastern and Midwest coal mining regions. The objectives of this project were to: 1) promote the safe and cost-effective closure by removal of coal ash impoundments, 2) harvest landfilled FGD, and 3) promote the high-volume beneficial use of these harvested CCRs in the reclamation of abandoned surface coal mine sites across the eastern and midwestern coal mining regions of the United States. The major tasks carried out for this project are summarized below: 1) Conesville Full-Scale Demonstration Project: About 2 million tons of harvested CCR materials from the closure by removal of an inactive fly ash pond and an adjacent old FGD landfill were used for the full-scale demonstration project to fully reclaim a nearby partially completed abandoned surface coal mine. Site monitoring for the project duration was carried out and results are discussed. 2) Laboratory Testing: Geotechnical and environmental testing of harvested ponded fly ash and landfilled FGD material at the former Conesville power plant were carried out. Completing the laboratory testing allowed for QA/QC for the full-scale site construction and informed the formulation of the risk analysis. 3) Risk Analysis: We developed a reliable computational model for fate and transport. We used these models and the rich set of monitored data for the Conesville site to analyze risks to human health and ecological risks associated with high-volume surface mine reclamation using harvested CCRs. 4) GIS Siting Study: A Geographic Information System (GIS) study was carried out for three states in the Eastern coal mining region and two states in the Midwest coal region. This effort provided site specific GIS information for five states and allowed us to establish protocols that other states can follow in implementing their own state specific GIS study.

01 COAL, LIGNITE, AND PEAT↗

Effect of Anoxic Iron Corrosion on WIPP Brine Geochemistry FY23 Final Report (U)

A 280-day study was completed to evaluate the effect of zero-valent iron (Fe 0 ) on the Waste Isolation Pilot Plant (WIPP) brine geochemistry under anticipated reducing conditions. Hydrogen (H 2 ) gas is expected to be present in the repository after closure due to the anoxic corrosion of a vast quantity of iron contained in the waste forms disposed at WIPP; therefore, a background argon atmosphere containing H 2 was chosen for this study. WIPP groundwater brine pH and E h will impact the mobility and fate of plutonium within the repository. Modeling and laboratory results for Castile WIPP brine indicate that equilibrium fa values relative to the standard hydrogen electrode (SHE) are 40 mV more reducing (i.e., more negative) than those for Salado WIPP brine (-480 mV vs. -440 mV, respectively) because of the higher pH of the Castile brine (pH 9 .3 for Castile vs. pH 8.8 for Salado). The E h and pH data were corrected for the effects of high ionic strength. The experimental results for both brines are consistent with thermodynamic predictions using OLI Systems' Mixed Solvent Electrolyte chemical equilibrium model. The measured and corrected pH and E h data from this study are provided in Table ES-I and Table ES-2, respectively. The experimental study, with four test conditions in triplicate, was performed in a dual glovebox with a nominally 3 vol.% H 2 in argon atmosphere (target H 2 range: 3 ± I vol.%). Simulants containing MgO only ( experimental control) and MgO+Fe 0 (WIPP base case) were prepared for both the Salado and Castile brines. MgO was included in all simulants to account for the use of bulk magnesium oxide in the WIPP repository. Fe 0 was included in some simulants to incorporate the effects of the anoxic corrosion of iron and in-situ hydrogen generation in the study. The brine compositions were developed by Sandia National Laboratory (SNL; Xiong, 2008) and have been used in previous WIPP evaluations. The test method (agitation, etc.) is partially based on ASTM D3987-12. Twelve rounds of periodic measurements of pH and E h were performed over the course of the study. Chemical analysis results for liquids and solids (ICP-MS, ICP-ES, IC Anion, TIC, SEM-EDX) are consistent with the pH, E h , and thermodynamic modeling results. This study included the following conditions that deviate from anticipated post-closure conditions following brine intrusion, but were selected to facilitate bench-scale testing to validate modeling of pH and E h for the post-closure WIP P repository: an anoxic glove box atmosphere containing ≤ 4 vol. % H 2 vs. substantially higher H 2 gas concentrations assumed in the WIPP Performance Assessment (PA); a significantly higher liquid-to-solid test ratio compared to the much lower phase ratio anticipated in the WIP P repository; agitation of the simulant bottles to maximize mass transfer; and finally the use of Fe 0 reagents having a much greater surface area than expected in the WIP P repository. Non-representative conditions were chosen for various reasons such as: to provide bounding conservative results, to provide a margin of safety for testing, or to facilitate simulant sub-sampling and analysis. In a parallel effort, aqueous electrolyte thermodynamic models were developed for the synthetic Salado and Castile brines to inform the experimental design, facilitate laboratory data interpretation, and allow extension of evaluations beyond the parameters tested. Thermodynamic modeling simulations including the MgO and Fe 0 additives that are directly relevant to the experimental measurements (e.g., pH calibration curve, ORP corrections) are included in this report. The measured fa of the simulants was close to the OLI model predictions for both brines and was largely controlled by the background H 2 partial pressure in the vapor phase as well as H 2 generated in situ in the aqueous phase by the Fe 0 corrosion. The H 2 gas-phase concentration tested and thermodynamically evaluated was much lower than is assumed in the WIPP PA; however, H 2 (g) concentrations significantly below this level are still predicted to result in very reducing conditions. In conclusion: • The experimental results are consistent with thermodynamic model predictions for fa, pH, and the effects of high ionic strength. • Evidence to date suggests that the H2 concentration in the glovebox atmosphere ultimately determined the final E h values of the simulants and resulted in highly reducing conditions. As a result, little difference was observed between the control simulants containing only MgO and the WIPP base-case simulants that contained MgO and Fe 0 . • This test methodology is recommended for future studies evaluating WIPP repository conditions. The methodology includes: (1) background H 2 in argon with agitation ( or could alternatively include in-situ-generated H 2 in sealed bottles); (2) carefully measured and corrected ORP data ( with much effort focused on allowing the probes to fully stabilize); and (3) ionic-strength-corrected pH data. Other best practices, such as simulant sparging/handling, ORP probe replacement, etc., should also be considered. • The coupling of experimental studies and thermodynamic modeling is also highly recommended because these methods inform and direct one another leading to greater confidence in and understanding of the results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Collaborative Research: Enabling multi-scale studies of magnetic reconnection with interpretable data-driven models

The development of accurate reduced descriptions and improved closures for magnetic reconnection is an important and a long‐standing challenge in plasma physics. The four‐fluid approach, and associated closures, that were investigated have the potential to improve the accuracy of plasma fluid models, capturing physical effects which would otherwise require a kinetic description. If successful, this approach could have an important impact for the modeling of laboratory and space plasmas. The major goals of this project were to develop new machine learning (ML) tools based on sparse and symbolic regression techniques, and to extract interpretable and generalizable reduced models (e.g., in the form of partial differential equations - PDEs) from data generated by first principles plasma simulations. Preserving interpretability of such data‐driven models is key to addressing the long‐standing theoretical and numerical challenges. Prior proof‐of‐principle studies have demonstrated the enormous potential of this approach, by recovering the well‐established hierarchy of plasma equations (from Vlasov to MHD) from data produced by particle‐in‐cell (PIC) simulations. Our goal in this project was to extend and apply these new tools to construct better kinetic closures for magnetic reconnection; to derive better models of particle injection and acceleration by this fundamental plasma process; and to use this understanding to accelerate the development of multi‐scale plasma algorithms. While our immediate focus was on the problem of magnetic reconnection, the tools that were will developed are general and applicable to other areas of plasma physics, and more broadly to many‐body phenomena. We anticipate that the development of these multi‐scale models will have a significant impact across different areas of plasma science, from fusion to space and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multifidelity Approach to Sensitivity Estimation in Large-Eddy Simulation

A novel approach to compute affordable approximate sensitivities in a large-eddy simulation (LES) is proposed and assessed. The approach is based on solving a Reynolds-averaged Navier–Stokes (RANS) problem that has been linearized around the mean LES solution, with closure modeling required for the linearized changes in turbulent Reynolds stresses. In the present study, the closure modeling is based on a linearized algebraic turbulence model. The method is assessed for the flow over a NACA0012 airfoil at a fixed angle of attack, with the Reynolds number as the variable parameter. The results show that, provided an accurate linearized closure model, the method predicts the correct sensitivity of the skin friction coefficient and the mean velocity field at the cost of a linearized RANS, which provides an important proof-of-concept for this approach. Here, the linearized algebraic turbulence model with standard model coefficients produces a reasonably accurate sensitivity, but the results also suggest that accuracy could be gained from recalibrating the model coefficients for this new use of the model.

42 ENGINEERING↗

Do We Know How to Model Reionization?

I compare the power spectra of the radiation fields from two recent sets of fully-coupled simulations that model cosmic reionization: “Cosmic Reionization On Computers” (CROC) and “Thesan”. While both simulations have similar power spectra of the radiation sources, the power spectra of the photoionization rate are significantly different at the same values of cosmic time or the same values of the mean neutral hydrogen fraction. However, the power spectra of the photoionization rate can be matched at large scales for the two simulations when the matching snapshots are allowed to vary independently. I.e., on large scales, the clustering of the radiation field in two simulations evolves similarly, but the exact timing of this evolution is different in different simulations and is not parameterized by an easily interpretable physical quantity like the mean neutral fraction or the mean free path. On small scales, large differences are present and remain partially unexplained. Both CROC and Thesan use the Variable Eddington Tensor approximation for modeling radiative transfer, but adopt different closure relations (optically thin OTVET versus M1). The role of this key difference is tested by using smaller simulations with a new cosmological simulation code that implements both closure relations in a controlled environment (the same hydro, cooling, and gravity solvers and the star formation recipe). In these controlled tests, both the M1 closure and the OTVET ansatz follow the expected behavior from a simple analytical approximation, demonstrating that the differences in the 2-point function of the radiation field induced by the choice of the Eddington tensor are not dominant.

79 ASTRONOMY AND ASTROPHYSICS↗

In-Situ Stress Measurements at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) Site

A scientific injection campaign was conducted at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site in 2017 and 2019. The testing included pump-in/shut-in, pump-in/flowback, and step rate tests. Various methods have been employed to interpret the in-situ stress from the test dataset. This study focuses on methods to interpret the minimum in-situ stress from step rate, pump-in/extended shut-in tests data obtained during the stimulation of two zones in Well 58-32. This well was drilled in low-permeability granitoid. A temperature of 199 °C was recorded at the well’s total depth of 2297 m relative to the rotary Kelly bushing (RKB). The lower zone (Zone 1) consisted of 46 m of the openhole at the toe of the well. Fractures in the upper zone (Zone 2) were stimulated between 2123–2126 m measured depths (MD) behind the casing. The closure stress gradient variation depended on the depth and the injection chronology. The closure stress was found to increase with the pumping rate/volume. This stress variation could indicate that poroelastic effects (“back stress”) and the presence of adjacent natural fractures may play an important role in the interpretation of fracture closure stress. Further, progressively increasing local total stresses may, consequently, have practical applications when moderate volumes of fluid are injected in a naturally fractured or high-temperature reservoir. The alternative techniques that use pump-in/flowback tests and temperature signatures provide a valuable perspective view of the in-situ stress measurements.

15 GEOTHERMAL ENERGY↗

Development of explosively driven launcher for meteoroid studies

The results of a continuing program to develop an explosively driven 2-stage hypervelocity launcher capable of achieving velocities between 15 and 20 km/sec are described. Previous efforts had identified incomplete barrel collapse as a limiting factor in launcher performance. Correlation of experimental and computational results obtained in the present study indicate that boundary-layer gases within the barrel act to prevent complete closure. Simplified calculations suggest that in-contact explosives may have insufficient energy densities to collapse the barrel against a developed boundary layer. Higher energy densities, sufficient to produce complete closure, were obtained with the use of steel flyer plates accelerated by a phased explosive lens. However, when flat flyer plates were impacted on the barrel, the sides of the barrel were observed to rupture and leak gas prior to barrel closure. A promising solution to this problem (untested) is to produce a symmetrical collapse with a cylindrical tube around the barrel.

Baum, D. W.↗