SEARCH · Engineering Papers

Results for “Injection data analysis”

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 721 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Free-spray characteristics and spray-wall interactions of methanol on a gasoline direct injector under flash-boiling and non-flash-boiling conditions

Methanol is considered a promising alternative fuel for internal combustion engines (ICEs) due to its high-octane number, fast laminar flame speed, and elevated latent heat of vaporization, all of which support higher compression ratios and improved thermal efficiency. However, its substantial latent heat of vaporization also poses cold-start challenges, such as misfire and fuel film deposition. This study aims to investigate methanol spray morphology and spray-wall interaction using the Spray M injector from the Engine Combustion Network within a constant-pressure flow vessel. A recently developed unified numerical framework capable of modeling both flash and non-flash boiling sprays is validated against experimental liquid volume fraction data acquired via 3-D computed tomography. Here, the results reveal that flash boiling significantly alters the spray morphology, leading to smaller droplets and spray collapse due to enhanced air-entrainment-induced turbulence. Quantitative agreement between experiments and simulations confirms this behavior. Coupled 0-D equilibrium and 3-D computational fluid dynamics analyses show that flash boiling accelerates evaporation and reduces fuel residence time, while non-flash conditions maintain a persistent liquid core more susceptible to wall wetting. Wall temperature diagnostics reveal that spray collapse alters heat transfer patterns by shifting cooling effects. Mixture fraction analysis indicates that evaporation is primarily governed by shear-layer turbulence, though deviations from adiabatic equilibrium mixing emerge under low-turbulence conditions. Finally, increasing fuel, ambient, and wall temperatures reduces wall wetting and film thickness, mitigating cold-start risks. These findings enhance the understanding of methanol sprays’ behavior and support its adoption as a viable, alternative fuel for ICEs.

Engine Combustion Network↗

DOE OSTI2023

Total metals, carbon, nitrogen & anion concentration data; Slate River & East River floodplains, Crested Butte, CO; May 2022-October 2022

This data package includes processed and undiluted measurements for metal, total carbon, total nitrogen, and anion concentrations from pore water (groundwater) and surface water samples from the Slate River and East River floodplains of Crested Butte, CO, focus field sites for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? Samples were collected between May and October of 2022. These measurements were all recorded at the Arizona Laboratory for Emerging Contaminants (ALEC) at the University of Arizona located in Tucson, AZ. Groundwater samples were extracted from a network of installed rhizon (Rhizosphere Research Products, part no. 19.60.21F, 0.6 micrometer mesh size) and piezometer wells within the river floodplain. All water samples were shaded from sun exposure during extraction from the subsurface and preserved at 4C until measured at ALEC.Analysis by ICP-MS (metals):Measurements for total metals were made on the Agilent 7700x ICP-MS (for total metals) – Agilent Technologies, Santa Clara, CA. The analytical QA/QC protocol was adapted from US EPA Method 200.8 for analysis by ICP-MS. Calibration standards were prepared from multi-element stock solutions (SPEX Certiprep, Metuchen, NJ). Calibration curves include at least 7 points with correlation coefficients > 0.995. The QC protocol includes a continuing calibration blank (CCB), a continuing calibration verification (CCV) solution and at least one quality control sample (QCS) to be analyzed just after calibration and again after every 12 samples and at the completion of the run. The QCS solutions are from an independent source, such as NIST SRM 1643e - Trace elements in water, or QCS solutions from High Purity Standards (Charleston, SC). Acceptable QC responses must be between 90 and 110% of the certified value. Lastly, a suitable internal standard (usually Rh, In, Ga or Ge) is added using on-line addition into the sample line and mixing tee.Analysis by Shimadzu TOC-L (TOC/TN):The TOC-L system is a combustion technique where liquid samples are injected and combusted into CO2 for carbon detection by non-dispersive infrared (NDIR) and NO for detection by chemiluminescence. A calibration curve using five standard solutions between 0.1 and 7 ppm for carbon and 0.05 and 3.5 ppm for nitrogen is made for each type of measurement with a linearity >0.99. All samples, standards, and QC’s are prepared in 24mL scintillation vials that have been baked for 4hrs at 475 Cº and made using RO water (18.2mΩ). QC’s include a calibration blank check (CCB), continuing calibration check (CCC), and a certified reference material check (CRM). All QC’s are within ±10% error and are run before and after each batch of samples. Samples are diluted and rerun if any measurement concentrations are above the highest standard.Analysis by Ion Chromatography (Anions):The instrument used is the Thermo Scientific Dionex ICS-6000 using AS+AG22 column set for anion analysis with sodium carbonate eluent. A calibration curve using five standard solutions between 5 and 250 umol/L is made with a linearity >0.99. Standards and QC’s are prepared in 15mL polypropylene conical tubes, pipetted along with the samples into 1.5mL polypropylene vials. Dilutions are made using RO water (18.2mΩ). QC’s include a calibration blank check (CCB), continuing calibration check (CCC), and a certified reference material check (CRM). All QC’s are within ±10% error and are run before and after each batch of samples. Samples are diluted and rerun if any measurement concentrations are above the highest standard.All files are in csv format.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1993

CFD simulation of coaxial injectors

The development of improved performance models for the Space Shuttle Main Engine (SSME) is an important, ongoing program at NASA MSFC. These models allow prediction of overall system performance, as well as analysis of run-time anomalies which might adversely affect engine performance or safety. Due to the complexity of the flow fields associated with the SSME, NASA has increasingly turned to Computational Fluid Dynamics (CFD) techniques as modeling tools. An important component of the SSME system is the fuel preburner, which consists of a cylindrical chamber with a plate containing 264 coaxial injector elements at one end. A fuel rich mixture of gaseous hydrogen and liquid oxygen is injected and combusted in the chamber. This process preheats the hydrogen fuel before it enters the main combustion chamber, powers the hydrogen turbo-pump, and provides a heat dump for nozzle cooling. Issues of interest include the temperature and pressure fields at the turbine inlet and the thermal compatibility between the preburner chamber and injector plate. Performance anomalies can occur due to incomplete combustion, blocked injector ports, etc. The performance model should include the capability to simulate the effects of these anomalies. The current approach to the numerical simulation of the SSME fuel preburner flow field is to use a global model based on the MSFC sponsored FNDS code. This code does not have the capabilities of modeling several aspects of the problem such as detailed modeling of the coaxial injectors. Therefore, an effort has been initiated to develop a detailed simulation of the preburner coaxial injectors and provide gas phase boundary conditions just downstream of the injector face as input to the FDNS code. This simulation should include three-dimensional geometric effects such as proximity of injectors to baffles and chamber walls and interaction between injectors. This report describes an investigation into the numerical simulation of GH2/LOX coaxial injectors. The following sections will discuss the physical aspects of injectors, the CFD code employed, and preliminary results of a simulation of a single coaxial injector for which experimental data is available. It is hoped that this work will lay the foundation for the development of a unique and useful tool to support the SSME program.

Landrum, D. Brian↗

DOE OSTI2020

Mercury Speciation During Vitrification of LAW

n this work, laboratory and engineering-scale tests were conducted that mimic the reactive environment of the Hanford Tank Waste Treatment and Immobilization Plant (WTP) low activity waste (LAW) melter and off-gas system in order to assess the speciation of mercury at select points in the processing system. The experimental protocols, test equipment, feed materials used, and rationale for testing are detailed in the Test Plan for this work. Particular attention was paid to the amount and speciation of mercury in submerged bed scrubber (SBS) solutions, which are intended to be recycled back to the melter feed at the WTP. Results are presented from detailed mercury analysis of these solutions using Environmental Protection Agency (EPA) Method 1630 and modifications to this method for dimethyl mercury and methyl mercury as well as EPA Method 1631 and modifications to this method for total, dissolved, elemental, suspended, and ionic mercury. Testing was conducted in two separate sets of experiments: (1) One set of tests (crucible-scale furnace tests) that involved heating small batches of mercury-spiked melter feed in crucibles to determine the amount of mercury retained in the glass and the species of mercury in the exhaust gases; (2) A second set of tests (DM10 melter tests) was planned that involved creating LAW melter plenum gas compositions using the DM10 melter system, injecting mercury into the off-gas stream, and passing that stream though a reactor that simulates various plenum gas conditions. Operational issues led to the need to use one of the VSL DM100 melters in place of the planned DM10 melter for these tests. The speciation of mercury after exposure to those conditions was monitored. In both sets of tests, the exhaust gases were run through a scrubber that was intended to mimic the LAW SBS in order to determine how much of the mercury exiting the melter would be retained in the primary off-gas system fluids and in what form. After passing through the SBS, the exhaust stream was analyzed to determine particulate, ionic, elemental, and total mercury passing downstream of the SBS. In tests employing gases derived from the DM100 melter that were spiked with elemental mercury, the processing system provided gas temperatures and residence times that are representative of the WTP LAW vitrification system in order to assess the effect of those conditions on mercury speciation. The Decontamination Factor (DF) across the system in the tests with mercury-spiked DM100 melter exhaust was determined using the analytical data from EPA Method 30B (Fluegas Adsorption Mercury Speciation (FAMS TM )) exhaust samples and the amount of mercury detected in the SBS solutions. Mercury species used in melter feed crucible scale tests were divalent (chloride and iodide), monovalent (chloride and fluoride), and elemental mercury. Individual tests included only a single form of mercury in the feed. The results for total mercury mass balances in the crucible tests are also presented.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS2017

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

NASA NTRS2017

Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

DOE OSTI2022

Effect of injected flux and current temporal phasing on self-organization in the HIT-SI3 experiment

The HIT-SI3 device at the University of Washington uses three oscillating inductive helicity injectors to form and sustain spheromak plasma equilibria. By adjusting the temporal phase of the injector waveforms with respect to each other, the toroidal spectrum of the imposed perturbations can be controlled. Using a recently implemented GPU-based control system, the available mode spectra were explored experimentally by scanning the space of relative injector phasing. In this space, significant variation in the toroidal mode spectrum ($n$ = 1, 2, 3) of the perturbations was observed. Additionally, variation in characteristics of driven equilibria was also observed, including a ≈ $30$% range in toroidal current gain ($I$ $Φ$ / $I$ $Inj$ ). Experimental results are compared with both a composite-equilibria and nonlinear dynamic model, including extended MHD simulations using the NIMROD code and composite Taylor state equilibria computed using the PSI-Tet code. In conclusion, qualitative agreement is seen with the nonlinear models, but not with composite-equilibria models, suggesting the use of nonlinear models to better capture observed plasma dynamics and provide predictive use for future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

COMBLE-ISLAS water isotopes in precipitation (COMBLEISLASISO)

During the measurement period, precipitation was sampled on daily to sub-daily basis with a sampling kit consiting of a clear plastic box (60x40x40 cm) mounted to a ground structure. Daily sampling lasted from 22 Feb to 24 Mar 2020. At a fixed time of day (10 UTC), samples were collected and the sampling box cleaned for the next sample. Additionally, during intense observations periods (IOPs) announced by ISLAS, ARMS scientists increased sampling frequency to sub-daily sampling. During the campaign period, 2 IOPs were conducted, with in total 12 samples taken. IOP1 lasted from 27 Feb 10 to 02 Mar 2020, and IOP2 lasted from 12 Mar to 14 Mar 2020. Solid precipitation in the box was then melted, transferred to an 8 ml glass vial, sealed and shipped to FARLAB (University of Bergen, Norway) for stable water isotope analysis. Measured stable isotope composition of precipitation samples will be related to conditions at other sampling locations, and at the evaporation site to inform about conservation of water isotope quantities. In the period from 22 Feb to 24 Mar, 23 daily samples were collected at 10 UTC for stable water isotope analysis. During the entire sampling period, a high variability of precipitation phase was encountered. Most samples were recovered in the form of graupel/snow (12 samples), as rain (7 samples), and the remaining samples fell as graupel, snow, rain or frozen rain. Precipitation amounts varied between <1 mm (26 Feb) and 140 mm (25 Feb) within the sampling box. Furthermore, higher-resolution sampling was conducted during 2 intense observation periods (IOPs). During IOP1, lasting from 27 Feb 10 UTC to 02 Mar 10 UTC, 7 samples were collected, all of them as solid phase precipitation. IOP2 lasted from 12 Mar 10 UTC to 14 Mar 09 UTC, and 5 samples (all as graupel/snow) were collected. Samples were processed according to FARLAB standard measurement procedures. In short, samples were transferred to 1.5 ml glass vials with rubber/PTFE septa (part #548-0907, VWR, USA). An autosampler (A0325, Picarro Inc) transferred ca. 2µl per injection into a high-precision vapourizer (A0211, Picarro Inc, USA) heated to 110°C. After blending with dry N2 (< 5 ppm H2O) the gas mixture was directed into the measurement cavity of a Cavity-Ring Down Spectrometer (L2140-i, Picarro Inc) for about 7 min with a typical water concentration of 20 000 ppm. Memory effects were reduced by two times measuring a vapour mixture at a mixing ratio of 50 000 ppm, obtained from 2 injections of 2 µl for 5 min at the beginning of each new sample vial. Thereafter, another 6 injections of 2 µl per sample were measured individually as described above, and averages of the last 5 injections were used for further processing. Three standards were measured at the beginning and end of each batch, including a drift standard DI2 (δD: -50.72±0.73 permil, δ18O: -7.63±0.10 permil), and for calibration the laboratory standards GLW (δD: -307.79±0.75 permil, δ18O: -40.02±0.07 permil) and EVAP2 (δD: 9.52±0.65 permil, δ18O: 1.81±0.13 permil). A detailed calibration report is included with the final uploaded data set.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

GeoThermalCloud: Cloud Fusion of Big Data and Multi-Physics Models using Machine Learning for Discovery, Exploration, and Development of Hidden Geothermal Resources

The primary goals of this project are exploring hidden geothermal resources in the U.S.A. and designing profitable enhanced geothermal systems (EGS). Many processes and parameters control geothermal exploration and energy production from geothermal fields. Diverse datasets (e.g., geology, geochemistry, geophysics, satellite, airborne geophysics) are available to help characterize subsurface geothermal conditions. Sparse and multi-scale characteristics of these datasets prohibit properly leveraging these datasets for geothermal exploration and profitable EGS design. Recent advancements in machine learning (ML) promise to resolve these issues. The tremendous challenges and risks of geothermal exploration and production bring the demand for novel ML methods and tools that can (1) analyze large field datasets, (2) assimilate model simulations (large inputs and outputs), (3) process sparse datasets, (4) perform transfer learning (between sites with different exploratory levels), (5) extract hidden geothermal signatures in the field and simulation data, (6) label geothermal resources and processes, (7) identify high-value data acquisition targets, and (8) guide geothermal exploration and production by selecting optimal exploration, production, and drilling strategies. To address these necessities, ML-based geothermal resources exploration and enhanced geothermal systems (EGS) design tools have been developed. The exploration tool is called GeoThermalCloud and EGS design tool is called GeoDT-ML. GeoThermalCloud (https://github.com/SmartTensors/GeoThermalCloud.jl) utilizes a LANL unsupervised ML platform called SmartTensors (https://tensors.lanl.gov/) to automate data analyses and interpretations by extracting hidden signatures to identify geothermal prospects. Also, it enables the identification of critical measurements needed to identify geothermal resource signatures. Alternatively, GeoDT-ML (https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS) is an ML-based alternative to GeoDT (https://github.com/GeoDesignTool/GeoDT.git), a fast, simplified multi-physics solver to evaluate EGS project designs in uncertain geologic systems. GeoDT-ML leverages recent advances in deep learning and high-performance computing. It is a faster and simpler version of GeoDT. To make this project a success, we used capabilities of LANL, PNNL, Google, Stanford, and Julia Computing. We analyzed eight datasets of the U.S.A. using GeothermalCloud and demonstrated potential highly prospective geothermal resources and identified key factors defining highly prospective sites. The first data set includes 44 locations in southwest New Mexico and 18 geological, hydrogeological, geophysical, geothermal, geochemical attributes. We defined low- and medium-temperature hydrothermal systems and discovered a new highly prospective site. The second data set analyzed 18 shallow water chemistry attributes at 14,342 locations in the Great Basin. It demarcated modestly, moderately, and highly prospective sites including key attributes for each type of prospectivity. The third data set analyzed Utah FORGE data including satellite (InSAR), geophysical (gravity, seismic), geochemical, and geothermal attributes. Here, we performed prospectivity analysis to identify future drilling locations using geological, geochemical, and geophysical attributes. Maps of temperature at depth and heat flow are constructed based on the available data. Prospectivity maps were generated, and drilling locations were proposed for future geothermal field exploration. The fourth data set analyzed 21 attributes at 120 locations in Tularosa Basin, New Mexico; data comes from past play fairway analyses in this region. ML analyses identified geothermal signatures associated with modestly, moderately, and highly hydrothermal systems. We also defined dominant attributes and spatial distribution of the geothermal signatures. The fifth, sixth, seventh, and eighth datasets include Tohatchi Springs, New Mexico, Hawaii, Brady site, Nevada, and EGS Collab, respectively. Moreover, we coupled GeothermalCloud and magnetotellurics data to pinpoint drilling locations for developing geothermal projects in the Tularosa Basin, New Mexico. GeothermalCloud found potential prospective locations for geothermal resources near White Sands Missile Range and McGregor Range at Fort Bliss. Magnetotellurics data determined the potential depth (~1800m) of geothermal prospects at McGregor Range based on apparent resistivity structures/layers in the subsurface. The McGregor Range consists of three resistivity layers and two resistivity structures. Magnetotellurics data also helps identify that the western portion of the McGregor Range has thick and low-resistivity earth materials. The low resistivity to the west is most likely for a fault system. Assuming temperature is consistent with a geothermal reservoir, the west-central part of the McGregor Range has the highest geothermal potential because of the increase in porosity and associated permeability attributed to the interpreted fault system. Also, we devised a coupling strategy between a process model and GeothermalCloud to characterize hydrogeological conditions and geothermal conditions, respectively. The process model characterizes hydrogeological and geothermal conditions on highly prospective geothermal sites provided by GeothermalCloud. We developed a physics-informed neural network (PINN) version of the Burns equation that can be easily coupled with GeothermalCloud. Furthermore, we performed an optimal design decision maximizing the economic value of an EGS power plant. This study optimized the range of well spacing between injection and production wells maximizing net present value in dollars (NPV). For this task, we used the GeoDT to simulate the Utah FORGE EGS development cycle from the initial well design to the end of production. Next, we accomplished another crucial task, which is predicting permeability of geothermal reservoirs. Predicting permeability of geothermal reservoirs is a non-trivial task because of huge computational runtime of simulation and lack of measurements. To avoid these limitations, we used easy-to-measure chemical concentrations in the subsurface as measurement data and convolutional neural network based ML model of a high-fidelity model. Next, we predicted permeability using Markov chain Monte Carlo simulation. We found that Markov chain Monte Carlo simulation predicts permeability with a high certainty if the prediction zone in the simulation area has chemical concentration data. Finally, we analyzed the DOE funded INGENIOUS and GeoDAWN projects data. For discovering hidden geothermal systems in the Great Basin, the INGENIOUS project accumulated old data, collected new data, and released them in 2022. The dataset includes a total of 24 geological, geophysical, and geochemical attributes. Data resolution and scale significantly vary prohibiting an appropriate usage. To avoid such limitations, we brought all data in the same resolution and scale by applying the inverse distance weighting interpolation technique for predicting data in unsampled locations. Subsequently, we analyzed LiDAR data of the GeoDAWN project. We received data in tiles format. The DOE’s overarching goal is to use ML on LiDAR data for finding favorable geological structures (e.g., step up faults in Brady, Nevada). To serve the purpose, we need to label favorable geologic structures that correspond to LiDAR data. We wrote an algorithm to label the LiDAR data with the favorable geologic structures.

15 GEOTHERMAL ENERGY↗

NASA NTRS2013

The XRD Amorphous Component in John Klein Drill Fines at Yellowknife Bay, Gale Crater, Mars

Drill fines of mudstone (targets John Klein and Cumberland) from the Sheepbed unit at Yel-lowknife Bay were analyzed by MSL payload elements including the Chemistry and Mineralogy (CheMin), APXS (Alpha Particle X-Ray Spectrometer), and Sample Analysis at Mars (SAM) instruments. CheMin XRD results show a variety of crystalline phases including feldspar, pyroxene, olivine, oxides, oxyhydroxides, sulfates, sulfides, a tri-octahedral smectite, and XRD amorphous material. The drill fines are distinctly different from corresponding analyses of the global soil (target Rocknest) in that the mudstone samples contained detectable phyllosilicate. Here we focus on John Klein and combine CheMin and APXS data to calculate the chemical composition and concentration of the amorphous component. The chemical composition of the amorphous plus smectite component for John Klein was calculated by subtracting the abundance-weighted chemical composition of the individual XRD crystalline components from the bulk composition of John Kline as measured by APXS. The chemical composition of individual crystalline components was determined either by stoichiometry (e.g., hematite and magnetite) or from their unit cell parameters (e.g., feldspar, olivine, and pyroxene). The chemical composition of the amorphous + smectite component (approx 71 wt.% of bulk sample) and bulk chemical composition are similar. In order to calculate the chemical composition of the amorphous component, a chemical composition for the tri-octahedral smectite must be assumed. We selected two tri-octahedral smectites with very different MgO/(FeO + Fe2O3) ratios (34 and 1.3 for SapCa1 and Griffithite, respectively). Relative to bulk sample, the concentration of amorphous and smectite components are 40 and 29 wt.% for SapCa1 and 33 and 36 wt.% for Griffithite. The amount of smectite was calculated by requiring the MgO concentration to be approx 0 wt.% in the amorphous component. Griffithite is the preferred smectite because the position of its 021 diffraction peak is similar to that reported for John Klein. In both cases, the amorphous component has low SiO2 and MgO and high FeO + Fe2O3, P2O5, and SO3 concentrations relative to bulk sample. The chemical composition of the bulk drill fines and XRD crystalline, smectite, and amorphous components implies alteration of an initially basaltic material under near neutral conditions (not acid sulfate), with the sulfate incorporated later as veins of CaSO4 injected into the mudstone.

Morris, Richard V.↗

NASA NTRS2001

Delta L: An Apparatus for Measuring Macromolecule Crystal Growth Rates in Microgravity

Strongly diffracting high quality macromolecule crystals of suitable volume are keenly sought for X-ray diffraction analysis so that high-resolution molecular structure data can be obtained. Such data is of tremendous value to medical research, agriculture and commercial biotechnology. In previous studies by many investigators microgravity has been reported in some instances to improve biological macromolecule X-ray crystal quality while little or no improvement was observed in other cases. A better understanding of processes effecting crystal quality improvement in microgravity will therefore be of great benefit in optimizing crystallization success in microgravity. In ground based research with the protein lysozyme we have previously shown that a population of crystals grown under the same solution conditions, exhibit a variation in X-ray diffraction properties (Judge et al., 1999). We have also observed that under the same solution conditions, individual crystals will grow at slightly different growth rates. This phenomenon is called growth rate dispersion. For small molecule materials growth rate dispersion has been directly related to crystal quality (Cunningham et al., 1991; Ristic et al., 1991). We therefore postulate that microgravity may act to improve crystal quality by reducing growth rate dispersion. If this is the case then as different, Materials exhibit different degrees of growth rate dispersion on the ground then growth rate dispersion could be used to screen which materials may benefit the most from microgravity crystallization. In order to assess this theory the Delta L hardware is being developed so that macromolecule crystal growth rates can be measured in microgravity. Crystal growth rate is defined as the change or delta in crystal size (defined as a characteristic length, L) over time; hence the name of the hardware. Delta L will consist of an optics, a fluids, and a data acquisition sub-assemblies. The optics assembly will consist of a video microscope camera mounted on three axis computer controlled translation stages. The fluids assembly consists of macromolecule and precipitant reservoirs, a temperature controlled growth cell and waste container, The data acquisition is achieved by using a frame-gabber, with images being stored on a hard drive. In operation, macromolecule and precipitant solution will be injected into the temperature controlled growth cell. As macromolecule crystals grow, the video microscope camera controlled by the translation stages, will be used to locate and record images of individual crystals, returning to the same crystals at specific time intervals. The images will be stored on the hard drive and used to calculate the crystal growth rate. To prevent vibrations interfering in the crystal growth rate measurements (Snell et al., 1997) Delta L will be used in connection with the Glovebox Integrated Microgravity Isolation Technology (g-LIMIT) inside the Microgravity Science Glovebox (MSG), onboard the International Space Station (ISS).

Judge, Russell A.↗

NASA NTRS2017

Board Level Proton Testing Book of Knowledge for NASA Electronic Parts and Packaging Program

This book of knowledge (BoK) provides a critical review of the benefits and difficulties associated with using proton irradiation as a means of exploring the radiation hardness of commercial-off-the-shelf (COTS) systems. This work was developed for the NASA Electronic Parts and Packaging (NEPP) Board Level Testing for the COTS task. The fundamental findings of this BoK are the following. The board-level test method can reduce the worst case estimate for a board's single-event effect (SEE) sensitivity compared to the case of no test data, but only by a factor of ten. The estimated worst case rate of failure for untested boards is about 0.1 SEE/board-day. By employing the use of protons with energies near or above 200 MeV, this rate can be safely reduced to 0.01 SEE/board-day, with only those SEEs with deep charge collection mechanisms rising this high. For general SEEs, such as static random-access memory (SRAM) upsets, single-event transients (SETs), single-event gate ruptures (SEGRs), and similar cases where the relevant charge collection depth is less than 10 μm, the worst case rate for SEE is below 0.001 SEE/board-day. Note that these bounds assume that no SEEs are observed during testing. When SEEs are observed during testing, the board-level test method can establish a reliable event rate in some orbits, though all established rates will be at or above 0.001 SEE/board-day. The board-level test approach we explore has picked up support as a radiation hardness assurance technique over the last twenty years. The approach originally was used to provide a very limited verification of the suitability of low cost assemblies to be used in the very benign environment of the International Space Station (ISS), in limited reliability applications. Recently the method has been gaining popularity as a way to establish a minimum level of SEE performance of systems that require somewhat higher reliability performance than previous applications. This sort of application of the method suggests a critical analysis of the method is in order. This is also of current consideration because the primary facility used for this type of work, the Indiana University Cyclotron Facility (IUCF) (also known as the Integrated Science and Technology (ISAT) hall), has closed permanently, and the future selection of alternate test facilities is critically important. This document reviews the main theoretical work on proton testing of assemblies over the last twenty years. It augments this with review of reported data generated from the method and other data that applies to the limitations of the proton board-level test approach. When protons are incident on a system for test they can produce spallation reactions. From these reactions, secondary particles with linear energy transfers (LETs) significantly higher than the incident protons can be produced. These secondary particles, together with the protons, can simulate a subset of the space environment for particles capable of inducing single event effects (SEEs). The proton board-level test approach has been used to bound SEE rates, establishing a maximum possible SEE rate that a test article may exhibit in space. This bound is not particularly useful in many cases because the bound is quite loose. We discuss the established limit that the proton board-level test approach leaves us with. The remaining possible SEE rates may be as high as one per ten years for most devices. The situation is actually more problematic for many SEE types with deep charge collection. In cases with these SEEs, the limits set by the proton board-level test can be on the order of one per 100 days. Because of the limited nature of the bounds established by proton testing alone, it is possible that tested devices will have actual SEE sensitivity that is very low (e.g., fewer than one event in 1 × 10(exp 4) years), but the test method will only be able to establish the limits indicated above. This BoK further examines other benefits of proton board-level testing besides hardness assurance. The primary alternate use is the injection of errors. Error injection, or fault injection, is something that is often done in a simulation environment. But the proton beam has the benefit of injecting the majority of actual SEEs without risk of something being missed, and without the risk of simulation artifacts misleading the SEE investigation.

Guertin, Steven M.↗

DOE OSTI2023

HFTS-1 (Final Report)

Despite the long history of hydraulic fracturing, the optimal number of fracturing stages during multi-stage fracture stimulation in horizontal wells can be a complex subject. In addition to the increased expense of multistage fracturing in horizontal wells, increasing the number of fracturing stages does not always correlate with a commensurate increase in production. The problem is the application of a uniform fracture stimulation design to all stages with no consideration for geological variations along the wellbore. The result is an inefficient use and costly waste of energy and water. Optimization of the fracturing process requires an understanding of the cause-and-effect relationship between fracturing parameters and local geological properties at a given location along the wellbore. Realizing that the generalized rock mechanics theories and hypotheses are not truly applicable to fractured and laminated shales, quantifiable impacts of a shale’s geomechanical and depositional features are a prerequisite for design and implementation of optimized hydraulic fractures. The overarching goal of this project is to understand and define the relationships of shale geology and fracture dynamics using detailed field data that includes coring of the fracture domain. Finally, resource recovery from shale formations is estimated to be less than fifteen percent. Research proposed in this project will establish the foundations for investigating enhanced recovery techniques for increased resource recovery in existing fracture treated wells. Natural gas as an EOR fluid is in the initial stages of being used broadly in the oil and gas industry and as such, many aspects of the process need to be researched and addressed prior to widespread acceptance. In addition, re-fracturing techniques, particularly using cemented liners also need to be further investigated. This project explores several of these EOR techniques in greater detail, not just at the theoretical level, but through actual field test sites. This report delves into the research program, identified tasks, several studies, observations, results, and conclusions under the Hydraulic Fracturing Test Site (HFTS) project. This test project is subdivided into three distinct phases and each phase is discussed in detail under separate sections. The original test site for Phase 1 as well as the Phase 2 site are in Midland Basin in Texas and operated by Laredo Petroleum. In the first two phases, the program targeted the Wolfcamp formation in Midland Basin. Phase 3 of the study targeted the EagleFord play in South Texas with Devon Energy as the operator. Significant conclusions from this study include: 1) lateral and vertical fracture growth mapping as well as proppant transport behavior, 2) Use of field gas and its applicability for huff n puff EOR injection in the Midland Basin, 3) Development and validation of proppant log, both at core scale as well as in-fill well implementation by analyzing drilling mud return samples, 4) Applicability of liner refracturing including best practices for enhanced unit level productivity. This report will detail all these studies as well as technical results and conclusions from said studies. These include general subsurface characterization work, completion designs as implemented, studies looking at fracture geometry, high resolution microseismic study, advanced crosswell strain mapping studies including drainage characterization, core characterization results, proppant analysis (core and observation wells), geochemistry, etc. to name a few.

02 PETROLEUM↗

NASA NTRS1999

Effect of color coding and subtraction on the accuracy of contrast echocardiography

BACKGROUND: Contrast echocardiography may be used to assess myocardial perfusion. However, gray scale assessment of myocardial contrast echocardiography (MCE) is difficult because of variations in regional backscatter intensity, difficulties in distinguishing varying shades of gray, and artifacts or attenuation. We sought to determine whether the assessment of rest myocardial perfusion by MCE could be improved with subtraction and color coding. METHODS AND RESULTS: MCE was performed in 31 patients with previous myocardial infarction with a 2nd generation agent (NC100100, Nycomed AS), using harmonic triggered or continuous imaging and gain settings were kept constant throughout the study. Digitized images were post processed by subtraction of baseline from contrast data and colorized to reflect the intensity of myocardial contrast. Gray scale MCE alone, MCE images combined with baseline and subtracted colorized images were scored independently using a 16 segment model. The presence and severity of myocardial contrast abnormalities were compared with perfusion defined by rest MIBI-SPECT. Segments that were not visualized by continuous (17%) or triggered imaging (14%) after color processing were excluded from further analysis. The specificity of gray scale MCE alone (56%) or MCE combined with baseline 2D (47%) was significantly enhanced by subtraction and color coding (76%, p<0.001) of triggered images. The accuracy of the gray scale approaches (respectively 52% and 47%) was increased to 70% (p<0.001). Similarly, for continuous images, the specificity of gray scale MCE with and without baseline comparison was 23% and 42% respectively, compared with 60% after post processing (p<0.001). The accuracy of colorized images (59%) was also significantly greater than gray scale MCE (43% and 29%, p<0.001). The sensitivity of MCE for both acquisitions was not altered by subtraction. CONCLUSION: Post-processing with subtraction and color coding significantly improves the accuracy and specificity of MCE for detection of perfusion defects.

Multicenter Study↗

DOE OSTI2021

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs Results

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells - increasing or decreasing the fluid flow rates across the wells - and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. Data and supporting literature from a study describing a new approach combining reservoir modeling and machine learning to produce models that enable strategies for the mitigation of decreased heat and power production rates over time for geothermal power plants. The computational approach used enables translation of sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy and discovery of optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an "open-source" reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 hours, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 seconds. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs. Includes a synthetic, yet realistic, model of a geothermal reservoir, referred to as open-source reservoir (OSR). OSR is a 10-well (4 injection wells and 6 production wells) system that resembles Brady Hot Springs (a commercially operational geothermal field in Nevada, USA) at a high level but has a number of sufficiently modified characteristics (which renders any possible similarity between specific characteristics like temperatures and pressures as purely random). We study OSR through CMG simulations with a wide range of flow allocation scenarios. Includes a dataset with 101 simulated scenarios that cover the period of time between 2020 and 2040 and a link to the published paper about this project, where we focus on the Machine Learning work for predicting OSR's energy production based on the simulation data, as well as a link to the GitHub repository where we have published the code we have developed (please refer to the repository's readme file to see instructions on how to run the code). Additional links are included to associated work led by the USGS to identify geologic factors associated with well productivity in geothermal fields. Below are the high-level steps for applying the same modeling + ML process to other geothermal reservoirs: 1. Develop a geologic model of the geothermal field. The location of faults, upflow zones, aquifers, etc. need to be accounted for as accurately as possible 2. The geologic model needs to be converted to a reservoir model that can be used in a reservoir simulator, such as, for instance, CMG STARS, TETRAD, or FALCON 3. Using native state modeling, the initial temperature and pressure distributions are evaluated, and they become the initial conditions for dynamic reservoir simulations 4....

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

Relocating microseismicity from downhole monitoring of the Decatur CCS site using a modified double-difference algorithm

SUMMARY The injection of CO2 at the Decatur carbon capture and storage site has generated significant microseismic activity, which occurs in distinct spatial clusters up to approximately 2.2 km from the primary injection well. Accurate and precise event locations are vital for the characterization of the microseismicity to help understand the reservoir response to the CO2 injection, whilst enabling the identification of minor faults and fractures below the resolution of conventional active seismic imaging. However, microseismic monitoring of fluid injection sites, such as Decatur, is often performed using a network of borehole sensors often from a single well. While these downhole sensors have excellent detection capabilities, their poor azimuthal coverage limits the ability to precisely determine event locations. We have developed a modified double-difference relocation algorithm suitable for both 1-D and 3-D velocity models, and which incorporates differential back azimuth observations to allow the benefits of the original double-difference algorithm to be applicable to a downhole microseismic monitoring setting. Applying the modified double-difference algorithm to the microseismicity at Decatur, we have successfully relocated 4293 events. The relocation included over 59 million observations for 757 285 event pairs, split across seven geographic regions. Despite the majority of observations being recorded in only two boreholes, with an almost identical azimuthal coverage, the results have shown to be reliable with significantly reduced residuals and low uncertainties associated with the final locations. We have analysed the residuals in terms of their association with each geographic region, data type, station and individual events, to fully appreciate their influence in the inversion and the fit of the data to the final set of event locations. For each region, the relocated seismicity has become less diffuse with improved clustering, and with newly visible linear features often orientated in a NE–SW direction. These results show the potential improvements that can be made to microseismic event locations recorded by a borehole network with a limited and variable azimuthal distribution.

58 GEOSCIENCES↗

DOE OSTI2020

Erosion-Corrosion Detection in Carbon Steel Pipe Loops using SRNLs Thickness and Mass Loss Measurement Coupons - 20464

Structural health of waste transfer infrastructure is of utmost importance to DOE and communities surrounding their active sites. A Fitness-for-Service (FFS) program for the Waste Transfer System has been implemented by the US Department of Energy (DOE) for the purposes of evaluating and quantifying the degradation of critical components. Currently, the method of approach is to section and evaluate pipelines that have been removed at end service. In order to augment data collection methods applied to making lifespan estimates of waste transfer system components, Florida International University (FIU) has been investigating new-technology based strategies for quantifying the rate of degradation. Present research considers the option of a coupon system installed on a section of pipe at points of interest that will provide easy access to high resolution measurements of erosion/corrosion rates during operation. Researchers at Savannah River National Laboratory (SRNL) have developed this novel solution for in situ erosion and corrosion measurements in many applicable locations such as tanks, pipes, ducts, and injection points. The coupons are approximately 12.7 mm (1/2 inch) in diameter and 25.4 mm (1 inch) in length and can be placed in the pipe sections using Swagelok compression fittings. The coupon design allows for monitoring of the total degradation resulting from corrosion and erosion. High-resolution Ultrasonic Thickness (UT) testing allows in-situ monitoring of both degradation mechanisms. When removed from the system mass loss resulting from corrosion and erosion may also be evaluated. Experimental testing at FIU implemented the SRNL coupon system in a 6.10 meter (20 foot) pipe loop system comprised of 50.8 mm (2 inch) and 76.2 mm (3 inch) pipe sections in order to replicate expected use cases at the sites. Operational conditions were emulated by pumping abrasive materials through the pipeline at full capacity and collecting measurements from the SRNL coupons. Initial testing was conducted by circulating water followed by a mixture of sand and water for uniform wear rates. Thickness changes and mass loss were measured, and the surface exposed to the process fluid was observed for visual evidence of characteristic wear patterns resulting from an erosive environment. The test matrix consisted of several hours of operation and using fine, medium and coarse grain sands. Test results were analyzed to check for all three parameters (mass loss/gain, thickness change and surface patters) which could provide estimates of the remaining useful life of the pipe sections. Thus, the research presents the findings of this evaluation of the SRNL coupons as a means for conducting high accuracy, cost effective in-situ analysis of pipeline and waste transfer systems for erosion and corrosion susceptibility. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2026

Hot-Carrier Injection and Millisecond Charge Separation from a Robust Heteroleptic Iron(II) Chromophore Immobilized on TiO 2

The synthesis, spectroscopic characterization, computational analysis, and photoelectrochemical behavior of a new iron-based chromophore, [(Cpy) 2 Fe- (deeb)](PF 6 ) 2 (Fe(Cpy) 2 (deeb)), where Cpy is 1-methyl-3-(2-pyridyl)imidazole and deeb is 4,4’-(CO 2 CH 2 CH 3 ) 2 -2,2’-bipyridine, is reported. Electrochemically reversible waves assigned to a metal-centered E o (Fe III/II ) = +0.48 and a ligand-centered E o (Fe 2+/+ ) = −1.47 V vs Fc +/0 reduction were evident in cyclic voltammetry measurements. The combination of a strong σ-donor and a π-acceptor lowered the energy of the metal-toligand charge-transfer (MLCT) excited state relative to the metal-centered state. Two MLCT transitions appear in the visible region at 424 and 580 nm. TDDFT calculations revealed that the lower-energy band was well formulated as Fe(II)→deeb, and the higher-energy transition was charge transfer to both the deeb and Cpy ligands. Resonance Raman spectroscopy supports these findings showing enhanced deeb vibrational modes with 532 nm excitation, both deeb and Cpy modes with 473 nm excitation, and exclusively Cpy with 405 nm excitation. Ultrafast spectroscopy reveals a short-lived (∼2 ps) MLCT excited state and a longer-lived (∼20 ps) metal-centered state. Efficient methods to deprotect the ester groups and anchor the complex to mesoporous TiO 2 (anatase) thin films in high surface coverages, Fe(Cpy) 2 (dcb)|TiO 2 σ = 3 × 10 −8 mol/cm 2 , were established. Pulsed light excitation of Fe(Cpy) 2 (dcb)|TiO 2 resulted in rapid excited state injection (k inj > 10 8 s −1 ) and formation of a charge-separated state, Fe III (Cpy) 2 (dcb)|TiO 2 (e), which persists on the millisecond time scale before returning cleanly to the ground state with secondorder kinetics. Injection yields measured 50 ns after light excitation were found to double from Φ = 0.15 with green (532 nm) light to 0.30 with blue (457 nm) light excitation. Incident photon-to-current efficiency (% IPCE) measurements as a function of excitation wavelength in a 0.5 M LiI/I 2 /CH 3 CN electrolyte provide clear evidence for band-selective “hot carrier” injection from the remote Cpy-localized excited state. Collectively, the spectroscopic and photoelectrochemical data indicate that a semiconductor can intercept hot electrons from iron chromophores even when the excited-state dipole is oriented away from the surface-anchoring ligand.

Charge transfer↗