SEARCH · Engineering Papers

Results for “UNDERGROUND”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Simulating High-Frequency Seismograms in Realistic Earth Models to Better Understand Source Discrimination Based on Differential Magnitudes ( M L− M c)

Discriminating low-yield underground nuclear explosions from small earthquakes is a key task in monitoring nuclear test ban treaties. P/S amplitude ratios have been an effective discriminant for moderate-sized events recorded at regional distances, but it is unclear if they are as effective in discriminating small seismic events recorded at local distances (<150 km). The difference between local magnitude (M L ) and coda duration magnitude (M c ) has been proposed as a new discriminant that may complement P/S amplitude ratios at local distances. Here, in this work, we calculate high-frequency (up to ∼4 Hz) synthetic seismograms at epicentral distances of 0–30 km in realistic models of the Salt Lake basin (Utah, United States) to better understand how variations in source type and depth affect M L −M c values. The Earth models incorporate simplified 1D and deterministic 3D structures, small-wavelength stochastic velocity perturbations, and surface topography. Coda waves are enhanced for the more complicated models compared to the base 1D model, but still underpredict observed durations by about a factor of two, which results in overprediction of amplitude to duration ratios (i.e., M L −M c values) for a near-surface explosion and a 7 km deep earthquake. For both source types, the predicted M L and M c values decrease as source depth increases, and M L −M c shows only minor variation with depth; however, M L −M c is on average ∼0.5 units smaller for explosions than earthquakes. This finding may imply that M L −M c has sensitivity to source type, in addition to being a depth discriminant, but more modeling is needed given the limitations of the current study. Future modeling should incorporate higher-frequency (≳5 Hz) simulations over a larger distance range (0–150 km), where M L and M c are commonly measured, while honoring low shear velocities (<300 m/s) near the surface and sampling a wider range of earthquake and explosion source mechanisms.

Hutchings, Sean J. [Univ. of Utah, Salt Lake City,↗

DOE OSTI2026

Seismic Event Characterization Using Full Moment Tensors on the Hypersphere

Moment tensor solutions provide insights into the deformation that has occurred in the source region of a seismic event and are therefore of great value in identifying different types of seismic sources, such as when monitoring for underground nuclear tests. Despite this utility, inversion of waveforms recorded by seismometers for their full seismic moment tensor is not yet routine, and development of robust methods to classify events based on this information is in its infancy. Here, we assemble an inventory of 1405 full moment tensor solutions that include explosive, earthquake, and collapse events, and investigate the use of anisotropic probability distribution functions on the 5D hypersphere to discriminate between these sources. Using a Bayesian classifier, we obtain optimal success rates of 98.4% across all events and demonstrate that modification of the prior probabilities provides a natural way to alter the balance between not missing desirable events (such as explosions) versus misclassifying large numbers of undesired events (such as earthquakes). The approach is specifically designed to progress from traditional, bipolar event screening metrics to more generalized event identification across multiple types of seismic sources. Despite current databases containing insufficient numbers of events to definitively demonstrate at present, we also find intriguing evidence of subgroupings within individual source populations on the hypersphere, for example, between chemical and nuclear explosions, raising the potential possibility of discriminating between these event types in the future.

Geosciences↗

DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

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

aquifer↗

DOE OSTI2023

WELLS Database

The Wellbore Exploration and Location Logistic System (WELLS) is a living national wellbore database - created and maintained by the National Energy Technology Laboratory (NETL). This resource contains more than seven million public wellbore records from state, federal, and tribal resources. Sourced from over 65 authoritative, yet disparate resources, the WELLS Database combines and synthesizes well data from oil, gas, underground injection, research, geothermal, geotechnical, groundwater and other types of wells in a single, unified system. This resource can be explored and visualized through the WELLS Interactive Application, also on EDX: https://edx.netl.doe.gov/dataset/wells-interactive-application The WELLS Database (formerly titled CO2-Locate) is an integrated national well dataset, representing open-source wellbore data from disparate state, tribal, and federal entities. The database provides publicly available well header data with key attributes such as well age, depth, and status. The database contains a fully integrated CSV file with all values in numerical columns, such as depth, converted into numbers. This version has a NETL derived API (American Petroleum Institute) number column and has been handled for redundancies, resulting in one record for every unique API number. The database also contains a fully integrated CSV file, where all original data are kept as text values. Additionally, the database includes a shapefile containing key attributes and coordinates from the integrated dataset, reformatted public wells CSV files, and a proprietary well density grid shapefile. Notes for consideration: The WELLS Database will be updated periodically with new datasets and information. A field dictionary with field (i.e., attribute) coverage across acquired public well resources, and the resulting integrated public well datasets are available in the spreadsheet, WELLS_Field_Dictionary.xlsx. Summary layers provided in this database are derived from proprietary layers and do not always contain key features (status, type, true vertical depth, or spud year) and therefore might not be shown when data are queried for those features.

AS↗

DOE OSTI2025

WELLS Interactive Application

The Wellbore Exploration and Location Logistic System (WELLS) Interactive Application is an interactive tool to enable easy exploration and visualization of the living national wellbore database (WELLS Database (https://edx.netl.doe.gov/dataset/wells_database)). The tool and underlying database were created and are maintained by the National Energy Technology Laboratory (NETL), providing visualization of the more than six million public wellbore records from more than 65 authoritative state, federal, and tribal resources. The WELLS Interactive Application serves up wellbore data from oil, gas, underground injection, research, geothermal, geotechnical, groundwater, and other types of wells in a single, standardized, unified system. In addition to the surface location of these wells, the underlying database combines select key attributes for features such as well age, depth, and operating status. The system also provides users with references back to the original sources used in this unified platform. The underlying data can be accessed through the WELLS Database: https://edx.netl.doe.gov/dataset/wells_database Additional Information: The WELLS Interactive Application (formerly titled CO2-Locate) enables visualization and access to the public wellbore records through an intuitive web-based mapping tool. The WELLS Interactive Application was designed to help users visualize, query, analyze, and download wellbore records. Public wellbore points are included as a layer in the Map page, called Public Wells. Additionally, a multivariate hexagon grid summarizing well density from proprietary well data, called Well Density, is included to identify data gaps between the public and proprietary well data. Filtering functionalities in the tool allow these two layers to be spatially filtered by state, county, or basin as well as by status, type, true vertical depth, and spud year. The WELLS Interactive Application also contains a Near Me tool can be used to search and explore wellbore data within a user-defined distance of a specified location on the map, which can also be downloaded. The Query tool allows users to query the selected or filtered wells in the Public Wells layer and export the data. For additional information on these tool functionalities, see the help documentation on the About page of the tool. Notes for Consideration: The Well Density layer provided in this application is derived from proprietary wellbore data, the records of which do not always contain values for key features (status, type, true vertical depth, or spud year). Therefore, data might not be available when layers are queried for all filter combinations. Additionally, visualizing layers and applying filters may take additional time to load (i.e., draw on the map) due to the large size of the data.

ccs↗

DOE OSTI2024

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

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

3-D seismic↗

DOE OSTI2021

Argonaut - A Robotic System for Cryogenic Environments

Fermilab and the HEP community invest significant resources into liquid argon detectors. The largest and most expensive of these detectors will be located in the Deep Underground Neutrino Experiment (DUNE). However, recent experiences have shown that there are limited avenues of monitoring, intervention, and interaction in the internal liquid environment. This proposal shows a technological path that could provide a valuable tool to ensure or at least improve the management of these HEP detectors. The development of a robotic system named Argonaut will demonstrate several technologies including 1) demonstration of suitable mobility of a small robotic device at liquid argon temperatures, 2) demonstration of wireless communication, 3) demonstration of improved diagnostics capabilities - such as tunable optics with motion control, 4) demonstration of interconnectivity of a robotic system with hardware residing within the detector. This initial research will be a seed for extended development in cold robotics and associated technologies. This work will allow FNAL to contribute a significant technology capability to recent efforts to cryogenic detector operations.

43 PARTICLE ACCELERATORS↗

DOE OSTI2021

Physics Studies for the LBNF Graphite Target Design

We present the simulated physics performance of the Long-Baseline Neutrino Facility (LBNF) graphite target that is being designed by the RAL High Power Targets Group for the Deep Underground Neutrino Experiment (DUNE). We first compare three conceptual cylindrical target design options as a function of target length (up to 2.2 m): downstream supported, two individual targets and an upstream-supported cantilever. Choosing the cantilever design as the baseline, we show the effect of widening the upstream inner conductor of the first focusing horn to provide extra space for supporting the target. We also give estimates of the expected performance of the 1.5 m prototype and 1.8 m production cantilevered targets. Furthermore, we show the effects of the main engineering updates made to the other two focusing horns since the DUNE TDR.

43 PARTICLE ACCELERATORS↗

DOE OSTI2021

Simulation Studies on the Interactions of Electron Beam with Wastewater

The manufactured chemical pollutants, like 1,4 dioxane and PFAS (per- and polyfluroralkyl substances), found in the underground water and/or drinking water are challenging to be removed or biodegraded. Energetic electrons are capable of mediating and removing them. This paper utilizes FLUKA code to evaluate the beam-wastewater interaction effects with different energy, space and divergence distributions of the electron beam. With 8 MeV average energy, the electron beam exits from a 0.0127 cm thick titanium window, travels through a 4.3 cm distance air and a second 0.0127 cm thick stainless water container window with 2.43 cm radius, and finally is injected into the water area, where the volume of water is around 75 cubic cm. The distribution parameters of the electron beam are from the GPT (General Particle Tracer) simulations for UITF (Upgraded Injector Test Facility) in Jefferson lab. By varying the distributions, several measurements including the dose (or energy deposition) distribution, electron fluence, photon fluence are scored and compared. Taking the comparisons into consideration, this paper is aiming to find better electron beams for the wastewater irradiation.

Li, X.↗

DOE OSTI2026

AI-driven neutrino diagnostics and radiation-hard beam instrumentation for next-generation neutrino experiments

The Long Baseline Neutrino Facility (LBNF) at Fermilab will deliver a high-intensity, multi-megawatt neutrino beam to the Deep Underground Neutrino Experiment (DUNE), enabling precision tests of the three-neutrino paradigm, CP violation searches, neutrino mass ordering determination, and supernova neutrino studies. To accelerate DUNE’s physics reach and ensure robust beam operations, we propose an integrated AI-driven framework with real-time diagnostics and radiation-hardened instrumentation. At its core is a Real-Time Beam Integrity Monitor using a physics-informed Digital Twin. By reconstructing pion phase space from muon profiles and exploiting magnetic horn optic linearity, it enables spill-by-spill beam correction and flux stabilization. By using this approach, flux-related systematics could be reduced from 5% to 1%, potentially accelerating the discovery of CP violations by four to six years. Complementing this, a US–Japan R&D effort will deploy a LAPPD-based muon monitor in the NuMI beamline. ToF measurements can be acquired with picosecond precision using this radiation-hard system, enhancing sensitivity to horn chromatic effects. Simulations confirm strong response to these effects. ML models predict beam quality and horn current to sub-percent accuracy from muon data, enhancing anomaly detection and stability. This scalable, AI-enabled strategy improves beam fidelity, reduces systematics, and sets a new standard for high-power accelerator operations.

Ganguly, Sudeshna [Fermilab] (ORCID:00000003163482↗

DOE OSTI2020

Fast and Accurate Core Analysis by the Full-Immersion Pressure-Pulse Decay: Part 1—Theory

Summary Core-scale measurements are considered the ground truth that oil and gas or electric utility companies use to predict the migration of fluids such as oil, natural gas, carbon dioxide, or brine deep underground during their extraction or injection operations. To provide a greater understanding of petrophysical properties of low-permeability geologic formations such as shales and tight gas sandstones, this study introduces a novel core-analysis procedure. The technique follows conventional pressure-pulse-decay permeametry, where the pressure in an inlet chamber adjacent to a cylindrical core plug undergoes a rapid pressurization, the system is shut in, and the pressure reaches a new equilibrium. However, unlike a standard unidirectional pressure-pulse decay, the full-immersion pressure-pulse decay (Hannon 2019) applies a pressure disturbance to the entire outer surface area of the sample. This article covers the numerical simulator designed to model flow through an anisotropic porous sample in this scenario. The model assumes distinct but uniform permeabilities along the radial and axial directions of the cylindrical plug sample. When extracting a plug vertically (or perpendicular to bedding), the permeability along the radial direction associates with the horizontal permeability (i.e., parallel to bedding), whereas in the axial direction, flow occurs perpendicular to bedding (a vertical permeability). Investigations of these model outputs demonstrate an approximately 20-fold decrease in time to complete a full-immersion experiment compared with conventional pressure-pulse decay. Furthermore, the pressure-decay curves resulting from the full-immersion method have slightly different shapes than those resulting from other unidimensional transient methods. These differences begin to demonstrate that, under achievable experimental conditions, the analysis of pressure data from one full-immersion test could enable the simultaneous estimation of the apparent permeabilities parallel and perpendicular to bedding of a cylindrical sample in addition to its porosity. A follow-up article finalizes the proof of this capability with a parameter-estimation procedure and presents experimental verification through a proof-of-concept study. Described in greater detail in that article, the parameter estimator requires multiple forward simulations to analyze the data, which behooves minimizing the compute time for each simulation. By using an alternating direction implicit time-marching scheme and a structured but variably spaced grid, the numerical simulator built for this purpose provides a forward model output with suitable accuracy in approximately 0.5 seconds.

Energy & Fuels↗

DOE OSTI2020

The new K, Pb-bearing uranyl-oxide mineral kroupaite: Crystal-chemical implications for the structures of uranyl-oxide hydroxy-hydrates

Abstract Kroupaite (IMA 2017-031), ideally KPb0.5[(UO2)8O4(OH)10]·10H2O, is a new uranyl-oxide hydroxylhydrate mineral found underground in the Svornost mine, Jáchymov, Czechia. Electron-probe micro-analysis (WDS) provided the empirical formula (K1.28Na0.07)Σ1.35(Pb0.23Cu0.14Ca0.05Bi0.03Co0.02Al0.01)Σ0.48 [(UO2)7.90(SO4)0.04O4.04(OH)10.00]·10H2O, on the basis of 40 O atoms apfu. Sheets in the crystal structure of kroupaite adopt the fourmarierite anion topology, and therefore kroupaite belongs to the schoepite-family of minerals with related structures differing in the interlayer composition and arrangement, and charge of the sheets. Uptake of dangerous radionuclides (90Sr or 135Cs) into the structure of kroupaite and other uranyl-oxide hydroxy-hydrate is evaluated based on crystal-chemical considerations and Voronoi-Dirichlet polyhedra measures. These calculations show the importance of these phases for the safe disposal of nuclear waste.

Geochemistry & Geophysics↗

DOE OSTI2021

Hydroxylpyromorphite, a mineral important to lead remediation: Modern description and characterization

Abstract Hydroxylpyromorphite, Pb5(PO4)3(OH), has been documented in the literature as a synthetic and naturally occurring phase for some time but has not previously been formally described as a mineral. It is fully described here for the first time using crystals collected underground in the Copps mine, Gogebic County, Michigan. Hydroxylpyromorphite occurs as aggregates of randomly oriented hexagonal prisms, primarily between about 20–35 μm in length and 6–10 μm in diameter. The mineral is colorless and translucent with vitreous luster and white streak. The Mohs hardness is ~3½–4; the tenacity is brittle, the fracture is irregular, and indistinct cleavage was observed on {001}. Electron microprobe analyses provided the empirical formula Pb4.97(PO4)3(OH0.69F0.33Cl0.06)Σ1.08. The calculated density using the measured composition is 7.32 g/cm3. Powder X-ray diffraction data for the type material is compared to data previously reported for hydroxylpyromorphite from the talc mine at Rabenwald, Austria, and from Whytes Cleuch, Wanlockhead, Scotland. Hydroxylpyromorphite is hexagonal, P63/m, at 100 K, a = 9.7872(14), c = 7.3070(10) Å, V = 606.16(19) Å3, and Z = 2. The structure [R1 = 0.0181 for 494 F>4σ(F) reflections] reveals that hydroxylpyromorphite adopts a column anion arrangement distinct from other members of the apatite supergroup due to the presence of fluorine and steric constraints imposed by stereoactive lone-pair electrons of Pb2+ cations. The F– anion sites are displaced slightly from hydroxyl oxygen anions, which allows for stronger hydrogen-bonding interactions that may in turn stabilize the observed column-anion arrangement and overall structure. Our modern characterization of hydroxylpyromorphite provides deeper understanding to a mineral useful for remediation of lead-contaminated water.

Geochemistry & Geophysics↗

DOE OSTI2023

Finchite, Sr(UO 2 ) 2 (V 2 O 8 )·5H 2 O, a new uranyl sorovanadate with the francevillite anion topology

Finchite (IMA2017-052), Sr(UO 2 ) 2 (V 2 O 8 )·5H 2 O, is the first uranium mineral known to contain essential Sr. The new mineral occurs as yellow-green blades up to ~10 µm in length in surface outcrops of the calcrete-type uranium deposit at Sulfur Springs Draw, Martin County, Texas, U.S.A. Crystals of finchite were subsequently discovered underground in the Pandora mine, La Sal, San Juan County, Utah, U.S.A., as diamond-shaped golden-yellow crystals reaching up to 1 mm. In this work, the crystal structure of finchite from both localities was determined using single-crystal X-ray diffraction and is orthorhombic, Pcan, with a = 10.363(6) Å, b = 8.498(5) Å, c = 16.250(9) Å, V = 1431.0(13) Å 3 , Z = 4 (R 1 = 0.0555) from Sulfur Springs Draw; and a = 10.3898(16), b = 8.5326(14), c = 16.3765(3) Å, V = 1451.8(4) Å 3 , Z = 4 (R 1 = 0.0600) from the Pandora mine. Electron-probe microanalysis provided the empirical formula (Sr 0.88 K 0.17 Ca 0.10 Mg 0.07 Al 0.03 Fe 0.02 ) Σ1.20 (UO 2 ) 2 (V 2.08 O 8 )·5H 2 O for crystals from Sulfur Springs Draw, and (Sr 0.50 Ca 0.28 Ba 0.22 K 0.05 ) Σ0.94 (U 0.99 O 2 ) 2 (V 2.01 O 8 )·5H 2 O for crystals from the Pandora mine, based on 17 O atoms per formula unit. The structure of finchite contains uranyl vanadate sheets based upon the francevillite topology. Finchite is a possible immobilization species for both uranium and the dangerous radionuclide 90 Sr because of the relative insolubility of uranyl vanadate minerals in water.

58 GEOSCIENCES↗

DOE OSTI2020

XXVIIth International Conference on Supersymmetry and Unification of Fundamental Interactions (SUSY 2019) (Final Report)

Supersymmetry (SUSY) is one of most elegant extensions of the Standard Model (SM) and explains the puzzles of the SM by providing a candidate to explain the dark matter content of the universe, allowing scientists to understand the origin of the electroweak scale requiring the top mass to be around 170 GeV and leading to the unification of forces at a grand unified scale. Further the minimal supersymmetric standard model (MSSM) predicts the Higgs boson mass to be less than 135 GeV. The discovery of the Higgs Boson with mass around 125 GeV at the LHC has provided a major support to SUSY ideas. Searches for SUSY are ongoing at the Large Hadron Collider (LHC). Direct and indirect dark matter experiments are searching for a particle dark matter candidate which arises most naturally in SUSY models. Proton decay predicted by SUSY grand unified theories is being searched for at deep underground experiments. In addition, recent advances in neutrino and dark matter physics, observational astrophysics, precision cosmology and the promising new window into the cosmos opened by the direct detection of gravitational waves, have brought new ideas on the potential connections between new fundamental particles and our understanding of their impact on the early universe and its evolution. At present, the major questions include: Is SUSY still the best candidate for models beyond the SM? Do we have any well motivated alternative to SUSY? Have we exhausted all possibilities to search for new physics at high and low energy scales? XXVIIth International Conference on Supersymmetry and Unification of Fundamental Interactions (SUSY 2019), hosted by Texas A&M University – Corpus Christi during May 20-24, 2019, provided a unique venue to discus and understand the status of SUSY, connection between particle physics and cosmology, supersymmetry and its alternative, Higgs sector, neutrino sector, flavor sector, dark matter, electroweak phase transition, astroparticle physics, gravitational waves and string theory. Discussion of results from the LHC, recent neutrino experiments and observations, direct and indirect dark matter detection experiments, detection of gravitational waves, data from particle colliders, as well as measurements of the CMB and Large Scale Structure were an integral part of SUSY 2019. To ensure the younger participants will benefit from the conference the most, the conference was preceded by the 4 day long pre-SUSY summer school for graduate students and postdocs. The invited speakers were leading scientists in the fields of SUSY interest. The school took place on Texas A&M University – Corpus Christi campus during the week prior the SUSY 2019 conference (May 15 – 18, 2019). Since its inception in 1993, SUSY has become one of the most important and widely attended international meetings in high energy physics, devoted to new ideas in fundamental particle physics. SUSY 2019 brought together approximately 250 scientists, theorists, phenomenologists, experimentalists and cosmologists, (including over 60 graduate students and 70 postdocs) representing 22 nations: Australia, Belgium, Canada, Chile, China, Colombia, France, Germany, India, Italy, Japan, Mexico, Peru, Portugal, Romania, South Korea, Spain, Sweden, Switzerland, Taiwan, United Kingdom and United States. SUSY 2019 provided a stimulating venue for the exchange of scientific ideas among experts in dark matter, neutrino physics, particle physics, astrophysics and cosmology. The following scientific topics were delivered during SUSY 2019 in form of 44 plenary talks and over 200 parallel talks: Unification of Forces; Electroweak, Top and Higgs Physics; Precision Calculations and MC tools; BSM in Flavor Physics; Neutrino Masses: Models and Phenomenology; Cosmology and Gravitational Waves; Dark Matter, Astroparticle Physics; Formal Field Theory and Strings; Alternatives to Supersymmetry; Quantum Information: Machine Learning/Big Data. 28 talks were given during the pre-SUSY program related to the following topics: Neutrino Physics; Big Data; Collider Physics & SUSY; String Phenomenology; Cosmology; Dark Matter; SUSY Models and Phenomenology

43 PARTICLE ACCELERATORS↗

DOE OSTI2020

Xyce Case Study

With the elimination of underground nuclear testing and declining defense budgets, science-based stockpile stewardship requires increased reliance on high performance modeling and simulation of weapon systems. Today's weapon systems are comprised of various electrical components and systems. As a result, there is a need for tools that will allow the use of massively parallel modeling and simulation techniques on high performance computers in existing and future weapons' electrical systems models. The Xyce Parallel Electronic Simulator is a SPICE (Simulation Program with Integrated Circuit Emphasis)- compatible circuit simulator designed to run on large-scale parallel computing platforms, though it can also execute efficiently on a variety of architectures including single processor workstations. As a mature platform for large-scale parallel circuit simulation, Xyce supports standard capabilities available in commercial simulators, in addition to various devices and models specific to Sandia's needs. Specifically, Xyce aids in the design and verification of electrical and electronic circuits and systems prior to weapons' manufacturing and deployment.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2020

North Dakota Integrated Carbon Storage Complex Feasibility Study. Final report

In spring 2017, the Energy & Environmental Research Center (EERC) initiated an effort to determine the feasibility of developing a commercial-scale CO 2 geologic storage complex able to store 50+ million tonnes (Mt) of CO 2 in central North Dakota safely, permanently, and economically. The objective was to fulfill the goals of the U.S. Department of Energy (DOE) Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative and address technical and nontechnical challenges specific to commercial-scale deployment of a CO2 storage project. The findings clearly show that the concept of capturing CO 2 from a lignite-fired electrical generation facility in central North Dakota and safely and permanently storing the CO 2 in the deep subsurface is indeed technically, economically, and socially feasible. This project evaluated two study areas and their respective geologic storage complexes located adjacent to separate coal-fired facilities in North Dakota: the Basin Electric Power Cooperative (BEPC)-owned Great Plains Synfuels Plant (GPSP) and the Minnkota Power Cooperative (Minnkota)-owned Milton R. Young Station (MRYS). These locations, one with CO 2 capture in place and an existing CO 2 pipeline, are bolstered by progressive North Dakota pore space ownership and long-term liability laws. These elements and a motivated team created an ideal synergistic scenario for ensuring success of the CarbonSAFE Initiative and promoting North Dakota’s statewide vision for carbon management. The project included drilling two new geologic characterization wells, integrating an existing 3-D seismic survey, creating a geologic model subsequently used for injection simulation, a risk assessment, public outreach, and generating a site development plan based on results. In addition, the performance of select National Risk Assessment Partnership tools was evaluated. The geologic characterization wells were drilled ~5600 feet deep to the Broom Creek Formation; ~350 feet of core was retrieved from each well. The core included the Broom Creek (targeted injection zone) and a portion of the overlying Opeche Shale (seal). The Flemmer-1 well, west of Beulah, North Dakota, yielded 169 feet of sandstone. The BNI-1 well located south of Center, North Dakota, yielded 124 feet of sandstone. In each case, laboratory analysis of the sandstone showed permeability in the 300–1000-mD range, with porosity of 20%–30%. The Flemmer-1 well was sited within the boundaries of an existing 3-D seismic survey. Colocating the well with the seismic survey maximized the relationship between new and legacy data and developed a first-of-its-kind interpretation of the geologic fabric of the Broom Creek. Geologic characterization data were integrated into a 5544-mi 2 geocellular model that encompassed both new wells and stratigraphy from the surface to the Amsden Formation (underlying the Broom Creek). The model was later expanded vertically to include the deeper Black Island–Deadwood interval to examine its potential viability as a storage target. The geocellular model provided the foundation for dynamic simulation of CO 2 into the Broom Creek. Results of the simulation suggest that the Broom Creek could accept the DOE target rate of 2 Mt/yr of CO 2 into as few as two wells. To bracket the expected capture from MRYS, simulations were also investigated for a 4-Mt/yr rate near MRYS. Although more wells are needed (two additional), the Broom Creek still has the storage resource to accept the CO 2 at the increased rate. A risk assessment exercise was conducted to identify and assess technical and nontechnical risks that could prevent potential candidate storage complexes within the study area from serving as commercial storage sites. The assessment identified and evaluated six technical risk categories: 1) CO 2 injectivity, 2) storage capacity, 3) lateral migration of CO 2 , 4) lateral pressure propagation, 5) vertical migration of CO 2 or formation brine, and 6) induced seismicity. Following two rounds of analysis and scoring, no risks were determined to preclude continued efforts to develop carbon capture, utilization, and storage (CCUS) in central North Dakota. The risk assessment results will be used to guide future site characterization, modeling and simulation, and monitoring activities. A specific nontechnical strategic risk based on challenges that may be realized in amalgamation of pore space resulted in vertically expanding the geologic model to incorporate the potential for stacked storage in multiple saline reservoirs. By including the Black Island–Deadwood interval (the basal sedimentary reservoir in this region), the amalgamated areal extent could be reduced by as much as 45%. Working with a smaller geographic area reduces risks and costs associated with monitoring and pore space leasing. An economic evaluation incorporating capture; transport (<5 mi); Class VI wells; permitting; and monitoring, verification, and accounting suggests implementing commercial-scale CCUS is economically attractive if the federal tax benefits of 45Q are included. This is validated by Minnkota’s continued pursuit of CO 2 capture and geologic storage at MRYS through its Project Tundra initiative, indicating that there is a business case for CCUS in central North Dakota. Currently, North Dakota is the only state with underground injection control (UIC) Class VI primacy. Built into the North Dakota Century Code is a series of regulatory requirements that guide the process to obtain a Class VI CO 2 storage facility permit. As part of this project, a site development plan was compiled to assure compliance with North Dakota’s requirements to permit a commercial-scale CO 2 storage operation and includes a prospective time line encompassing a general breakdown of activities. In total, an estimated 30 months would be needed to execute the necessary steps to attain a North Dakota CO 2 storage facility permit. Outreach was an integral part of the project and encompassed any project-related activity that had contact or exposure beyond the project team. The goals of outreach were to foster an environment from which stakeholders could make informed decisions about the project and gauge community receptiveness to a CCUS project. A consistent set of messages and outreach products were developed in conjunction with an outreach advisory board that integrated project partners and team members. To gauge public acceptability of geologically storing CO 2 , 5611 households in the project area were invited to participate in an online survey. The survey results indicate that the public attitude regarding CCUS is neutral to positive, with strong sentiment that CO 2 capture and storage may be an approach to maintain the economic vitality of the region. To achieve project objectives, critical support in the form of financial backing, engineering evaluations, site access, outreach collaboration, operations data, risk assessment/evaluation, and software access was provided by BEPC, the North Dakota Industrial Commission Lignite Research Council, ALLETE Clean Energy, BNI Energy, North American Coal Corporation, Minnkota, Prairie Public Broadcasting, Computer Modelling Group Ltd., and Schlumberger.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2020

Xyce Case Study

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗