SEARCH · Engineering Papers

Results for “SECONDARY INJECTION”

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.

108 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Unraveling the hysteretic behavior at double cations-double halides perovskite - electrode interfaces

Despite over a decade of research on metal halide perovskites (MHPs) in the context of photovoltaic applications, understanding the nature of electronic and ionic processes associated with current-voltage (I-V) hysteretic behavior has been limited. Here, we explore the hysteretic behavior in (FAPbI 3 ) 0.85 (MAPbBr 3 ) 0.15 perovskite devices with lateral Cr electrodes by applying first order reversal curve (FORC) bias waveform in I-V, Kelvin probe force microscopy (KPFM) measurements, and in-situ chemical imaging by time-resolved time-of-flight secondary ion mass spectrometry (tr-ToF-SIMS). In dark, we reveal pronounced hysteretic behaviors of charge dynamics in the off-field by probing time-dependent current and contact potential difference (CPD). Under illumination, transient and hysteretic behaviors are significantly reduced. The tr-ToF-SIMS results reveal that the hysteretic behaviors are strongly associated with accumulation of Br- ions at the interfaces. In addition, the low mobility MA + ions result in transient behavior and contribute to the hysteretic phenomena. It was shown that Pb 2+ ions can be reduced at the interfaces due to electrochemical reactions with the electrode in the presence of charge injection and photogenerated charges. Furthermore, these hysteretic behaviors associated with charge dynamics, ion migration, and interfacial electrochemical reaction are critical to further improve the performance and stability of MHPs photovoltaics and optoelectronics.

42 ENGINEERING↗

DOE OSTI2023

High-energy neutrino emission from magnetized jets of rapidly rotating protomagnetars

ABSTRACT Relativistic jets originating from protomagnetar central engines can lead to long duration gamma-ray bursts (GRBs) and are considered potential sources of ultra-high-energy cosmic rays and secondary neutrinos. We explore the propagation of such jets through a broad range of progenitors, from stars which have shed their envelopes to supergiants which have not. We use a semi-analytical spin-down model for the strongly magnetized and rapidly rotating protoneutron star (PNS) to investigate the role of central engine properties such as the surface dipole field strength, initial rotation period, and jet opening angle on the interactions and dynamical evolution of the jet-cocoon system. With this model, we determine the properties of the relativistic jet, the mildly relativistic cocoon, and the collimation shock in terms of system parameters such as the time-dependent jet luminosity, injection angle, and density profile of the stellar medium. We also analyse the criteria for a successful jet breakout, the maximum energy that can be deposited into the cocoon by the relativistic jet, and structural stability of the magnetized outflow relative to local instabilities. Lastly, we compute the high-energy neutrino emission as these magnetized outflows burrow through their progenitors. Precursor neutrinos from successful GRB jets are unlikely to be detected by IceCube, which is consistent with the results of previous works. On the other hand, we find that high-energy neutrinos may be produced for extended progenitors like blue and red supergiants, and we estimate the detectability of neutrinos with next generation detectors such as IceCube-Gen2.

Astronomy & Astrophysics↗

DOE OSTI2020

Structural Characterization of Potential Carbon Dioxide Reservoirs and Adjacent Strata within the Llandovery Silurian to Middle Devonian Strata of Ohio

The Midwest Regional Carbon Sequestration Partnership (MRCSP) has incorporated the work of geologic research teams (Geoteams) in its regional characterization, project planning and carbon dioxide (CO 2 ) injection implementation work since the partnership was established by the U.S. Department of Energy (DOE) in 2003. Over this 16-year period, the cohort of Geoteams has grown from five to ten states and has contributed to the characterization of geologic sequestration opportunities, refinement of reservoir and seal data, and supported injection efforts through both predictive and post-injection assessments. This report details the work led by the Ohio Department of Natural Resources, Division of Geological Survey to evaluate the structure of nine formations in the Appalachian Basin of Ohio to assess seal integrity and whether carbon dioxide (CO 2 ) migration pathways could pose a potential risk. As part of Phase II MRCSP work, the Pennsylvania Geological Survey led the characterization of the Llandovery Silurian (“middle Silurian”) to Middle Devonian regional geologic characterization CO 2 utilization and storage. This preliminary assessment included generalized, regional-scale mapping, as well as the development of relatively coarse structure and isochore maps of stratigraphic intervals of interest for the Appalachian Basin. Because this stratigraphic interval offers several enhanced oil recovery targets and secondary targets for CO 2 storage, the Ohio Geological Survey undertook efforts in Phase III to characterize the Ohio’s portion of the Appalachian Basin more precisely and to resolve inconsistencies within Ohio’s formation tops. This study refined the subsurface correlations of the Silurian-through-Devonian stratigraphic section in Ohio. Higher-resolution isochore mapping of the CO 2 storage reservoirs and seals provides greater detail of potential high-porosity zones and potential structural influence on the interval.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2024

An integrated experimental–modeling approach to identify key processes for carbon mineralization in fractured mafic and ultramafic rocks

Abstract Controlling atmospheric warming requires immediate reduction of carbon dioxide (CO2) emissions, as well as the active removal and sequestration of CO2 from current point sources. One promising proposed strategy to reduce atmospheric CO2 levels is geologic carbon sequestration (GCS), where CO2 is injected into the subsurface and reacts with the formation to precipitate carbonate minerals. Rapid mineralization has recently been reported for field tests in mafic and ultramafic rocks. However, unlike saline aquifers and depleted oil and gas reservoirs historically considered for GCS, these formations can have extremely low porosities and permeabilities, limiting storage volumes and reactive mineral surfaces to the preexisting fracture network. As a result, coupling between geochemical interactions and the fracture network evolution is a critical component of long-term, sustainable carbon storage. In this paper, we summarize recent advances in integrating experimental and modeling approaches to determine the first-order processes for carbon mineralization in a fractured mafic/ultramafic rock system. We observe the critical role of fracture aperture, flow, and surface characteristics in controlling the quantity, identity, and morphology of secondary precipitates and present where the influence of these factors can be reflected in newly developed thermo-hydro-mechanical–chemical models. Our findings provide a roadmap for future work on carbon mineralization, as we present the most important system components and key challenges that we are overcoming to enable GCS in mafic and ultramafic rocks.

58 GEOSCIENCES↗

DOE OSTI2024

Simulations of stand-off runaway electron beam termination by tungsten particulates for tokamak disruption mitigation

Abstract Stand-off runaway electron termination by injected tungsten particulates offers a plausible option in the toolbox of disruption mitigation. Tungsten is an attractive material choice for this application due to large electron stopping power and high melting point. To assess the feasibility of this scheme, we simulate runaway collisions with tungsten particulates using the MCNP program for incident runaway energies ranging from 1 to 10 MeV. We assess runaway termination from energetics and collisional kinematics perspectives. Energetically, the simulations show that 99% of runaway beam energy is removed by tungsten particulates on a timescale of 4–9 µ s. Kinematically, the simulations show that 99% of runaways are terminated by absorption or backscattering on a timescale of 3–4 µ s. By either metric, the runaway beam is effectively terminated before the onset of particulate melting. Furthermore, the simulations show that secondary radiation emission by tungsten particulates does not significantly impact the runaway termination efficacy of this scheme. Secondary radiation is emitted at lower particle energies than the incident runaways and with a broad angular distribution such that the majority of secondary electrons emitted will not experience efficient runaway re-acceleration. Overall, the stand-off runaway termination scheme is a promising concept as a last line of defense against runaway damage in ITER, SPARC, and other future burning-plasma tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Simulating Droplet-Resolved Haze and Cloud Chemistry Forming Secondary Organic Aerosols in Turbulent Conditions within Laboratory and Cloud Parcels

Most of our existing knowledge of cloud chemistry in regards to forming secondary organic aerosols (SOA) is based on measurements in bulk aqueous solutions. However, SOA reaction kinetics derived from bulk solution measurements might differ from the kinetics in actual cloud droplets, since turbulent mixing and ionic strengths, and ratio of surface area to volume in individual cloud droplets might vary substantially from bulk solutions in the real atmosphere. Three-dimensional models at various scales have been used to simulate aqueous chemistry. However, most of these models do not resolve turbulence down to the smallest length scales of 1 mm and do not simulate cloud chemistry in individual cloud droplets due to large computational costs. Here we incorporate the formation of isoprene epoxydiol SOA (IEPOX-SOA) in individual droplets within a one-dimensional explicit mixing parcel model (EMPM-Chem). We apply EMPM-Chem to simulate turbulence and droplet-resolved IEPOX-SOA formation using a configuration based on the Michigan Tech Pi chamber. We find that the dissolution of IEPOX gases is weighted more towards larger cloud droplets due to their large liquid water content (compared to smaller droplets), while the conversion of dissolved IEPOX to IEPOX-SOA is much greater within smaller deliquesced haze particles due to their higher acidity and ionic strengths compared to cloud droplets. We also find that as droplet residence times increase in the chamber, e.g., due to increasing ammonium bisulfate seed aerosol injection rates and/or increasing heights of the chamber, formation of IEPOX-SOA increases substantially. Thus, our EMPM-Chem model could be used to design future cloud chambers to maximize SOA production from cloud chemistry. We also apply the EMPM-Chem model to simulate how IEPOX-SOA formation evolves in individual cloud droplets within rising cloudy parcels in the atmosphere. We find that as subsaturated air is entrained into and turbulently mixed with the cloud parcel, evaporation causes a reduction in droplet sizes, which leads to corresponding increases in per droplet ionic strength and acidity. Increased droplet acidity in turn greatly accelerates the kinetics of IEPOX-SOA formation. In conclusion, our results provide key insights into single-cloud-droplet chemistry, suggesting that entrainment mixing may be an important process that increases SOA formation in the real atmosphere.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Dynamic Surface Incorporation of Pb 2+ Ions at the Actively Dissolving Calcite (104) Surface

The reaction of dissolved Pb 2+ with calcite surfaces at near-equilibrium conditions involves adsorption of Pb 2+ and precipitation of secondary heteroepitaxial Pb-carbonate minerals. A more complex behavior is observed under far-from-equilibrium conditions, including strong inhibition of calcite dissolution, development of microtopography, and near-surface incorporation of multiple monolayers (ML) of Pb 2+ without precipitation of secondary phases [where 1 ML ≡ 1 Ca/20.2 Å 2 , the crystallographic site density of the calcite (104) lattice plane]. However, the mechanistic controls governing far-from-equilibrium reactivity are not well understood. Here, in this study, we observe the interfacial incorporation of dissolved Pb 2+ during the dissolution of calcite (104) surfaces at pH ~3.7 in a flow-through reaction cell, revealing the formation of a ~1 nm thick Pb-rich calcite layer with a total Pb coverage of ~1.4 ML. These observations of the sorbed Pb distribution used resonant anomalous X-ray reflectivity, X-ray fluorescence, and nanoinfrared atomic force microscopy. We propose that this altered surface layer represents a novel sorption mode that is stabilized by conditions of sustained disequilibrium. This behavior may significantly impact the transport of dissolved metals during disequilibrium processes occurring in acid mine drainage and subsurface CO 2 injection and, if appropriately accounted for, could improve the predictive capability of geochemical reactive-transport models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Driving Zn-MnO2 grid-scale batteries: A roadmap to cost-effective energy storage

Highlights Zn-MnO 2 batteries promise safe, reliable energy storage, and this roadmap outlines a combination of manufacturing strategies and technical innovations that could make this goal achievable. Approaches such as improved efficiency of manufacturing and increasing active material utilization will be important to getting costs as low as $100/kWh, but key materials innovations that facilitate the full 2-electron capacity utilization of MnO 2 , the use of high energy density 3D electrodes, and the promise of a separator-free battery with greater than 2V potential offer a route to batteries at $50/kWh or less. Abstract Large-scale energy storage is certain to play a significant, enabling role in the evolution of the emerging electrical grid. Battery-based storage, while not a dominant form of storage today, has opportunity to expand its utility through safe, reliable, and cost-effective technologies. Here, secondary Zn–MnO 2 batteries are highlighted as a promising extension of ubiquitous primary alkaline batteries, offering a safe, environmentally friendly chemistry in a scalable and practical energy dense technology. Importantly, there is a very realistic pathway to also making such batteries cost-effective at price points of $50/kWh or lower. By examining manufacturing examples at the Zn–MnO 2 battery manufacturer Urban Electric Power, a roadmap has been created to realize such low-cost systems. By focusing on manufacturing optimization through reduced materials waste, scalable manufacturing, and effective materials selection, costs can be significantly reduced. Ultimately, though, coupling these approaches with emerging research and development advances to enable full capacity active materials utilization and battery voltages greater than 2V are likely needed to drive costs below a target of $50/kWh. Reaching this commercially important goal, especially with a chemistry that is safe, well-known, and reliably effective stands to inject Zn–MnO 2 batteries in the storage landscape at a critical time in energy storage development and deployment. Graphical abstract

Spoerke, Erik D. (ORCID:0000000322862560)↗

DOE OSTI2020

Characterizing Pore-Scale Geochemical Alterations in Eagle Ford and Barnett Shale from Exposure to Hydraulic Fracturing Fluid and CO 2 /H 2 O

As demand increases for an affordable energy source that is tied to an environmental obligation to reduce greenhouse gas emissions and water usage, there is a growing consideration in shale production utilizing processes such as 1) enhancing hydrocarbon recovery via carbon dioxide (CO 2 ) flooding, 2) using CO 2 as a fracturing agent to minimize water use, and 3) storing CO 2 in depleted shale formations to mitigate emissions to the atmosphere. Understanding the geochemical reactions and alterations that occur as shale is exposed to fluids and CO 2 is necessary to develop and optimize each of these processes for field applications. While the majority of shale formations are stimulated using traditional fracturing fluid, some may be fractured using CO 2 or other non-traditional means. We examine the effect fracturing fluid has on shale and how it behaves with secondary exposure to dry CO 2 or CO 2 -saturated water using in situ Fourier Transform infrared spectroscopy (FTIR), feature relocation scanning electron microscopy (SEM), and surface area and pore size analysis using volumetric gas sorption. These techniques were performed on Eagle Ford and Barnett shale samples that were exposed to fracturing fluid and unexposed (as received). Shales that have been exposed to traditional fracturing fluid experienced two reaction fronts. The first reaction front was formed during exposure to the fracturing fluid (pH of ~1.4). A secondary reaction front was formed as a result of CO 2 -saturated fluid exposure in the form of carbonic acid (pH ~5.6). These two different reaction mechanisms drove multiple dissolution and precipitation cycles which altered petrophysical properties of the shale and could lead to a significant impact on flow pathways. FTIR showed that equilibration of carbonate dissolution and precipitation cycles could take as long as 35 days. Samples exposed to fracturing fluid showed significantly less carbonate reactivity compared to those exposed to water. Pore size analysis results indicate exposure to fracturing fluid blocked small nanopores (10-15 nm) reducing BET surface area and total pore volume. SEM results show barite precipitated heavily during exposure to fracturing fluid. It appeared that carbonic acid was able to extract sulfur from organic matter to form gypsum evaporites. The mineralogical (barite precipitation and calcite dissolution/precipitation) and pore-scale alterations observed in these samples may lead to enhancement of flow pathways for injected CO 2 or produced hydrocarbons.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2022

Kimberlina 1.2 CCUS Geophysical Models and Synthetic Data Sets

This synthetic multi-scale and multi-physics data set was produced in collaboration with teams at the Lawrence Berkeley National Laboratory, National Energy Technology Laboratory, Los Alamos National Laboratory, and Colorado School of Mines through the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative. Data are associated with the following publication: Alumbaugh, D., Gasperikova, E., Crandall, D., Commer, M., Feng, S., Harbert, W., Li, Y., Lin, Y., and Samarasinghe, S., “The Kimberlina Synthetic Geophysical Model and Data Set for CO2 Monitoring Investigations”, The Geoscience Data Journal, 2023, DOI: 10.1002/gdj3.191. The dataset uses the Kimberlina 1.2 CO2 reservoir flow model simulations based on a hypothetical CO2 storage site in California (Birkholzer et al., 2011; Wainwright et al., 2013). Geophysical properties models (P- and S-wave seismic velocities, saturated density, and electrical resistivity) were produced with an approach similar to that of Yang et al. (2019) and Gasperikova et al. (2022) for 100 Kimberlina 1.2 reservoir models. Links to individual resources are provided below: [CO2 Saturation Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-co2-saturation-models); Resistivity Models – [part 1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-1), [part 2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-2), and [part 3](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-3); [Vp Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vp-velocity-models); [Vs Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vs-velocity-models); [Density Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-density-models). The 3D distributions of geophysical properties for the 33 time stamps of the SIM001 model were used to generate synthetic seismic, gravity, and electromagnetic (EM) responses for 33 times between zero and 200 years. Synthetic surface seismic data were generated using 2D and 3D finite-difference codes that simulate the acoustic wave equation (Moczo et al., 2007). 2D data were simulated for six point-pressure sources along a 2D line with 10 m receiver spacing and a time spacing of 0.0005 s. 3D simulations were completed for 25 surface pressure sources using a source separation of 1 km in both the x and y directions and a time spacing of 0.001 s. Links to individual resources are provided below: [2D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-velocity-models) and [2D surface seismic data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-surface-seismic-data). [3D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-velocity-models), and 3D seismic data [year0](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year0), [year1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year1), [year2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year2), [year5](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year5), [year10](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year10), [year15](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year15), [year20](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year20), [year25](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year25), [year30](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year30), [year35](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year35), [year40](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year40), [year45](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year45), [year49](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year49), [year50](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year50), [year51](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year51), [year52](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year52), [year55](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year55), [year60](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year60), [year65](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year65), [year70](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year70), [year75](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year75), [year80](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year80), [year85](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year85), [year90](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year90), [year95](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year95), [year100](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year100), [year110](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year110), [year120](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year120), [year130](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year130), [year140](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year140), [year150](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year150), [year175](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year175), [year200](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year200). The Python scripts to read these models and data are provided [here](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-python-scripts). EM simulations used a borehole-to-surface survey configuration, with the source located near the reservoir level and receivers on the surface using the code developed by Commer and Newman (2008). Pseudo-2D data for the source at [2500 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz2500m) and [3025 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz3025m), used a 2D inline receiver configuration to simulate a response over 3D resistivity models. The [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-csem-data) contain electric fields generated by borehole sources at monitoring well locations and measured over a surface receiver grid. Vector gravity data, both on the surface and in boreholes, were simulated using a modeling code developed by Rim and Li (2015). The simulation scenarios were parallel to those used for the EM: [pseudo-2D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were calculated along the same lines and within the same boreholes, and [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were simulated over 3D models on the surface and in three monitoring wells. A series of [synthetic well logs](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-well-logs) of CO2 saturation, acoustic velocity, density, and induction resistivity in the injection well and three monitoring wells are also provided at 0, 1, 2, 5, 10, 15, and 20 years after the initiation of injection. These were constructed by combining the low-frequency trend of the geophysical models with the high-frequency variations of actual well logs collected in the Kimberlina 1 well that was drilled at the proposed site. Measurements of permeability and pore connectivity were made on cores of Vedder Sandstone, which forms the primary reservoir unit: [CT micro scans](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-ct-micro-scans-of-vedder-formation) and [Industrial CT Images](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-industrial-ct-images-vedder-formation). These measurements provide the range of scales in the otherwise synthetic data set to be as close to a real-world situation as possible. References: Birkholzer, J.T., Zhou, Q., Cortis, A. and Finsterle, S., 2011. A sensitivity study on regional pressure buildup from large-scale CO2 storage projects. Energy Procedia, 4, 4371-4378. Commer, M., and Newman, G.A., 2008. New advances in three-dimensional controlled-source electromagnetic inversion, Geophysical Journal International, 172, 513-535. Gasperikova, E., Appriou, D., Bonneville, A., Feng, Z., Huang, L., Gao, K., Yang, X., Daley, T., 2022, Sensitivity of geophysical techniques for monitoring secondary CO2 storage plumes, Int. J. Greenh. Gas Control, Volume 114, 103585, ISSN 1750-5836, https://doi.org/10.1016/j.ijggc.2022.103585. Moczo, P., J.O. Robertsson and L. Eisner, 2007, The finite-difference time-domain method for modeling of seismic wave propagation: Advances in geophysics, 48, 421-516. Rim, H., and Y. Li, 2015, Advantages of borehole vector gravity in density imaging, Geophysics, 80, G1-G13. Wainwright, H. M.; Finsterle, S.; Zhou, Q.; Birkholzer, J. T., 2013. Modeling the Performance of Large-Scale CO2 Storage Systems: A Comparison of Different Sensitivity Analysis Methods. International Journal of Greenhouse Gas Control, 17, 189205. https://doi.org/10.1016/j.ijggc.2013.05.007, DOI: 10.18141/1603331. Yang, X., Buscheck, T.A., Mansoor, K., Wang, Z., Gao, K., Huang, L., Appriou, D., and Carroll, S.A., 2019. Assessment of geophysical monitoring methods for detection of brine and CO2 leakage in drinking water aquifers, International Journal of Greenhouse Gas Control, 90, 102803, https://doi.org/10.1016/j.ijggc.2019.102803.

CCUS↗

DOE OSTI2024

Dynamic Binary Complexes (DBC) as Super-Adjustable Viscosity Modifiers for Hydraulic Fracturing Fluids

In the preceding project year two, we refined three DBC formulations from a selection of over 50 different chemistries. The optimization study primarily encompassed testing for (i) reversibility extent, (ii) performance in the presence of chemical additives, (iii) adhesion and friction behavior during displacement in wellbores and pipelines, (iv) corrosion protection performance, and (v) injection performance with model fracture systems at the laboratory scale. Highly promising results obtained from all these tests signify the significant potential of DBCs in enhancing hydrocarbon recovery from unconventional reservoirs. The primary activities in the third project year included publishing experimental findings across multiple articles and conducting outreach initiatives. Throughout the year, we undertook tasks such as replicating experimental results, further optimizing various formulations and their associated experimental sets, and conducting additional tests to address missing components based on reviewer feedback and suggestions. We also explored the surfactant and friction-reduction aspects of selected formulations through drag reduction tests. In addition, we constructed an improved fracturing performance setup and performed flow injection tests. The specific DBC formulations focused on during this project period were A8/B1, A12/B5, and A10/B12. We also obtained results for additional DBC formulations and a select few commercial fracturing fluids for the purpose of comparison. Within the project's scope, we aim to enhance the experimental findings with the development of various models. The first two years focused on two key aspects: (i) the creation of a high-fidelity hydraulic fracturing model for non-Newtonian fluids to gain insights into the implementation of DBC fluids in fracking environments, and (ii) the development of a multiphase flow simulator for estimating total production, fluid saturation in the reservoir, and the creation of a fracture propagation model and kinetic Monte Carlo (kMC) models for diverse applications. In the third year, we delved into the fundamental nanostructural properties of DBCs, exploring aspects such as material chemistry, pH tunability, and control of DBC formation and stability. Subsequently, in the extension year, we conducted a systematic investigation of various building blocks containing primary, secondary, and tertiary amine functional groups to understand their impact on rheological and viscoelastic properties. Furthermore, we explored a Dissipative Particle Dynamics (DPD) model to simulate self-assembly processes with precision, creating a high-fidelity representation of relevant nanostructures. The tasks performed this year with the significant results obtained have been discussed in Section 2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Numerical study on spray collapse process of ECN spray G injector under flash boiling conditions

This work features a numerical study of the spray collapse process that occurs under flash-boiling, gasoline-direct-injection conditions. The Engine Combustion Network's Spray G injector was modeled, and iso-octane, n-hexane, and propane were used as fuels with the goal to investigate fuel property effects. The development of individual jets and jet-to-jet interaction were analyzed to reveal the physical mechanisms of spray collapse. The results agreed qualitatively well with the experimental data of propane available from the literature and the predictions revealed the model's ability to capture the near-nozzle shock structures witnessed in the experiments. Simulations showed that the flow development of flashing sprays was strongly tied with the fuel properties. Specifically, under-expanded jets were not observed with iso-octane due to the fuel's low vapor pressure. In the case of n-hexane, the jets were under-expanded, and a low-pressure core was observed inside the jets due to the flow expansion. With propane, the very strong expansion caused the flow to become transonic and generate a unique shock structure near the nozzle exit. The shock structure revealed the formation of primary shock cells within the jets and secondary shock cells induced by jet-to-jet interaction. Based on the results, the spray collapse under flash boiling conditions was mainly attributed to two effects, i.e., the decrease in pressure at the spray center induced by the low-pressure cores inside and between the jets, and the isolation effect due to the jet-to-jet interaction which contributed to maintain the pressure differential between the spray center and the ambient.

42 ENGINEERING↗

DOE OSTI2020

Surface Reactivity Analysis of the Crude Oil–Brine–Limestone Interface for a Comprehensive Understanding of the Low-Salinity Waterflooding Mechanism

Low-salinity waterflooding (LSWF) has proven to improve oil recovery in carbonate formations through rock wettability alteration, although the underlying mechanism remains elusive. Multivalent ionic exchange and calcite dissolution have usually been investigated using geochemical analysis in secondary coreflooding. In this work, coreflooding, in tertiary mode, coupled with a surface reactivity analysis approach was employed to investigate the interplay of wettability alteration mechanisms such as mineral dissolution, electrostatic bond attraction, and the effect of pH at in situ conditions. Improved oil recovery (IOR) in tertiary mode observed by coreflooding in Indiana limestone rocks showed an ionic strength dependence, that is, reducing brine ionic strength resulted in an increase in oil recovery. Coreflooding results showed that the seawater and low-salinity brines deprived of Mg 2+ ions resulted in the lowest IOR in tertiary mode, indicating the significance of Mg 2+ on IOR in limestone rocks. Similar results were observed through the contact angle measurement showing the limestone rock wettability state dependence on ionic strength and the effect of Mg2+ ions. Surface reactivity analysis showed an increase in solution pH, Ca 2+ and Mg 2+ ions concentration in the effluent solution from the coreflooding in tertiary mode using low salinity brines (about 40 and 20% increase in the effluent composition for Ca 2+ and Mg 2+ , respectively). These changes in solution composition were used to calculate the in situ oil–brine and rock–brine zeta potential using a validated surface complexation model, showing the changes of zeta potential as brine is injected into limestone rocks. The results show that using seawater-like brine in tertiary mode resulted in no mineral dissolution or ionic exchange. However, improved oil recovery (IOR) using such seawater-like brine was due to wettability alteration caused by reduced electrostatic bond attraction associated with Mg 2+ ions [from 2.6 × 10 –13 (mol/m 2 ) 2 for formation water salinity to 1.5 × 10 –13 (mol/m 2 ) 2 for seawater salinity]. Using low-salinity brines in tertiary mode improved oil recovery by mineral dissolution, resulting in oil desorption and an increase in solution pH. Finally, the increase in solution pH also resulted in reduced electrostatic bond attraction which lead to rock wettability alteration using low-salinity brines.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2024

Mechanisms of Waterflood Inefficiency: Analysis of Geological, Petrophysical and Reservoir History, a Field Case Study of FWU (East Section)

The petroleum reservoir represents a complex heterogeneous system that requires thorough characterization prior to the implementation of any incremental recovery technique. One of the most commonly utilized and successful secondary recovery techniques is waterflooding. However, a lack of sufficient investigation into the inherent behavior and characteristics of the reservoir formation in situ can result in failure or suboptimal performance of waterflood operations. Therefore, a comprehensive understanding of the geological history, static and dynamic reservoir characteristics, and petrophysical data is essential for analyzing the mechanisms and causes of waterflood inefficiency and failure. In this study, waterflood inefficiency was observed in the Morrow B reservoir located in the Farnsworth Unit, situated in the northwestern shelf of the Anadarko Basin, Texas. To assess the potential mechanisms behind the inefficiency of waterflooding in the east half, geological, petrophysical, and reservoir engineering data, along with historical information, were integrated, reviewed, and analyzed. The integration and analysis of these datasets revealed that several factors contributed to the waterflood inefficiency. Firstly, the presence of abundant dispersed authigenic clays within the reservoir, worsened by low reservoir quality and high heterogeneity, led to unfavorable conditions for waterflood operations. The use of freshwater for flooding exacerbated the adverse effects of sensitive and migratory clays, further hampering the effectiveness of the waterflood. In addition to these factors, several reservoir engineering issues played a significant role in the inefficiency of waterflooding. These issues included inadequate perforation strategies due to the absence of detailed hydraulic flow units (HFUs) and rock typing, random placement of injectors, and uncontrolled injected fresh water. These external controlling parameters further contributed to the overall inefficiencies observed during waterflood operations in the east half of the reservoir. A detailed understanding of the mechanistic factors of inefficient waterflood operation will provide adequate insights into the development of the improved recovery technique for the field.

Morgan, Anthony (ORCID:0000000211519153)↗

DOE OSTI2019

Single-fuel reactivity controlled compression ignition through catalytic partial oxidation reformation of diesel fuel

A single-fuel RCCI concept has been proposed to avoid the need for a secondary fuel system required for conventional RCCI by generating the secondary fuel from the primary fuel through catalytic partial oxidation (CPOX) reformation. In conventional RCCI, gasoline or natural gas can be used as the low-reactivity fuel, and diesel can be used as the high-reactivity fuel. In this study, two reformate mixtures generated by reforming diesel fuel at different operating conditions were used as the low-reactivity fuel, with the parent diesel as the high reactivity fuel. The combustion characteristics of reformate-diesel RCCI were compared with the conventional RCCI. A CFD model was also developed and validated against the experimental results. The model was then used to validate a necessary approximation of the reformate mixture's species concentrations. Compared to conventional RCCI fuel pairs, reformate-diesel RCCI shows marginally better thermal efficiency, approximately 10% better THC emissions, approximately 50% lower NOx emissions, and good controllability. Because the reformate mixture has a high concentration of diluents it displaces a large fraction of intake air and acts similarly to EGR. The combustion behavior of reformate-diesel RCCI is in between that of gasoline-diesel and natural gas-diesel conventional RCCI. From the results, it can be concluded that reformate-diesel RCCI is not overly sensitive to the reformation process itself and the exact species concentrations in the reformate mixture. Finally, a small change in the start of injection of diesel, blend ratio, and EGR fraction can be used to compensate for reformate mixture concentration differences.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2022

Wabash CarbonSAFE. Geologic Analysis of the Potosi Dolomite Reservoir Interval and Potential Confining Units, at the Wabash #1 Well, Vigo Co., Indiana, Subtask 7.2, Technical Report

The Wabash CarbonSAFE project drilled the Wabash #1 stratigraphic test well (ID# 168045) at the Wabash Valley Resources (WVR) IGCC facility in Vigo County, Indiana, to evaluate the feasibility of commercial-scale CO 2 storage near the site. The Mt. Simon Sandstone was the initial target for storage evaluation and was found to have generally poor reservoir qualities in the Wabash #1 well. Focus was then placed on the secondary target, the Potosi Dolomite. For storage and containment evaluation an extensive suite of geophysical logs was collected throughout the borehole, and whole cores were obtained from the Mt. Simon Sandstone, Eau Claire Formation, and Maquoketa Group. The Maquoketa Group is identified as a regionally extensive secondary sealing interval above the Mt. Simon Sandstone and as a primary seal for the Potosi Dolomite storage complex. Due to concerns over potential lost circulation zones in the Potosi Dolomite during drilling, no core samples or Formation Micro Imager (FMI) logs were acquired in the Potosi Dolomite. However, cased hole well testing was performed over a 20 feet (6 m) interval within the Potosi Dolomite during which a fluid swab sample was obtained. Approximately 35 miles of 2D seismic information was acquired in the project area near Wabash #1 to aid in reservoir and caprock characterization. The highly porous and permeable zones in the Potosi Dolomite and lower Oneota Dolomite observed within the Wabash #1 well are, for the purpose of this report, considered the same reservoir package and referred to as the Potosi reservoir interval. This report presents a regional depositional context for the Potosi reservoir interval, the overlying confining zones in the Knox Group, and Maquoketa Group primary seal. This context is developed through detailed information obtained from the Wabash #1 well such as lithologic data collected from cuttings and core, geophysical logging, geomechanical analysis of core samples, and well testing and fluid sampling within the Potosi reservoir interval. Information discussed in this report was compiled as a basis for the development of static and dynamic geologic models of the Potosi Dolomite storage complex. The static and dynamic modeling of CO 2 injection in the Potosi reservoir interval are discussed in a separate report (Dessenberger et al., 2022) under the Wabash CarbonSAFE project.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2023

Fuel Injection Dynamics and Composition Effects on RDE Performance

Rotating detonation engines (RDEs) provide a promising route to substantially increasing cycle efficiency in stationary gas turbines. Much of this increase relies on the ability to achieve consistent pressure gain within the combustor. In particular, the design of injectors that feed fuel and air into the detonation channel plays a crucial role. Such injectors have to ensure proper mixing of fuel and oxidizer, while minimizing backflow of detonation products into the feed plenums, and reduce susceptibility to the complex wave structures that exist within the combustor. From a practical perspective, such RDEs also need to operate with variable fuel composition. When fuel mixtures with components that possess vastly different oxidation pathways and time-scales are used, there could be additional losses through deflagrative burning instead of detonation-driven heat release. Such sensitivity to the complete flow path is akin to the physics of thermoacoustic instabilities in conventional gas turbines. In this sense, RDEs pose a unique research challenge: the performance of the device relies on the small-scale heat release process, which is highly dependent on the flow interactions within the full-scale system. As a result, canonical flow configurations, instrumented with detailed diagnostics or modeled using high-fidelity tools, but only focus on the small-scale processes will not contain the key system-level interactions. At the same time, macroscopic measurements and models that only capture system-level performance will not provide insight into the key sources of pressure losses. These couplings and sensitivities provide a formidable challenge to both experimental and simulation studies of the effects Thus, a joint experimental/computational program designed specifically to address these challenges was undertaken in this program. The focus of this program was on two key topics: a) the interaction between injector flow and the overall wave dynamics within the combustor, and b) the deflagration/detonation structure in multi-component fuels that are of practical interest. Both topics involve interaction of small-scale heat release processes with the geometry-dependent wave structure. Studies focused on the study of full-scale RDE systems, based on a 6-inch conventional annular geometry. Experimentally RDEs were studied using a combination of diagnostics. A combination of optical diagnostics and aero-thermo-acoustic analysis based on a combination of spectral and mode decomposition analysis was used to identify the dynamics of the detonation wave and other secondary waves that exist in the system. These studies have helped the identification and investigation of injector and detonation dynamics arising from coupling, and how they affect RDE mixing, detonation structure, operability and performance. Performance of RDEs was investigated through thrust stand measurements, which was used to evaluate the effective pressure gain generated by the system through the concept of equivalent available pressure. Optical diagnostics were developed and implemented to investigate the distribution of heat release, across the detonation wave. Novel optical diagnostics of NIR imaging was also developed and applied to investigate the high temperature / high pressure distribution across the detonation wave. In order to complement the experiments, the computational tools were geared to simulate the full experimental setup. GPU-based acceleration of the models and computations were developed to enable rapid simulation of the full system. In addition, the use of adaptive mesh refinement, and unstructured grid formulation, enabled the investigation of realistic geometries studied in the laboratory. The simulations produced a wealth of detail on the structure of the detonation wave under different operating conditions. Emphasis was placed on quantifying mixture pre-burning and the impact on wave propagation and structure.

03 NATURAL GAS↗

DOE OSTI2021

Spectral Induced Polarization-Biogeochemical Relationships for Remediation Amendment Monitoring

Geophysical tools such as electrical resistivity (ER) can indirectly monitor subsurface changes in response to remedial injections. These methods exhibit relatively low spatial resolution compared to sediment core characterization but are advantageous due to the ability to collect measurements non-intrusively over time across large volumes of the subsurface. Moreover, along with confirmatory groundwater or core sampling, geophysical methods can be used during active biogeochemical remedies to monitor short-term contaminant transformations and mobility, as well as part of an overall strategy for long-term monitoring of subsurface contaminated sites. The use of alternating current spectral induced polarization (SIP) provides significantly more information than conventional geophysical methods like direct current ER. SIP allows for monitoring changes in both solution and surface conductivity by separation of real and imaginary conductivity, respectively, as well as surface capacitance. In principle, SIP can measure indicators of remedy progression such as precipitation reactions that sequester contaminants, potentially providing a better indication of amendment delivery and reactivity as compared to conventional ER methods. However, multiple processes and material properties have overlapping (interacting) electrical responses across a range of frequencies. Hence, the purpose of this scoping study was to evaluate if SIP can be used to monitor (a) amendment delivery and (b) precipitation and reactivity of amendments under consideration at Hanford. SIP measurements were collected in fully saturated columns packed with sand and Hanford formation sediments containing (a) amendments that were highly conductive with significant capacitance (zero valent iron – ZVI, sulfur modified iron – SMI) and (b) amendments that exhibited low electrical conductivity with a small capacitance (calcite, apatite, bismuth). The sand was a quartz material with homogenous particle size that exhibited a relatively low surface conductivity. It was used as a control for comparison with the sediments from the Hanford Site, which have a greater surface conductivity because of their complex mineralogy and heterogeneous particle size distribution and may have complex interactions with amendments. The amendment mass fraction was varied to represent the different stages and subsurface locations associated with the delivery of a remedy. The primary objective was to identify the solution and solid surface changes associated with the delivery amendments and their secondary reactions within the subsurface. The figure below summarizes results for the amendments tested through this project. The SIP phase shift or imaginary conductivity change for the high conductivity amendments was more than 10 times that of the low conductivity amendments, highlighting the relative ease of detection of ZVI and SMI independent of the background signal from sand or Hanford sediments. The ZVI phase shift and imaginary conductivity changes occur primarily at high frequency (> 100 Hz) whereas SMI changes were at low frequency (0.01 to 10 Hz). SIP signals of SMI also increased over time and the maximum shifted to lower frequencies. The low conductivity amendments (calcite, abiotic and biotic apatite, bismuth subnitrate) exhibited relatively small phase and imaginary conductivity changes when added to sediments. The changes were above minimum detection limits (0.5 mrad for phase shift, 0.03 µS/cm for imaginary conductivity) for the highest concentration except for the commercial bismuth material. The lowest amendment concentration that can be detected is likely sediment specific, as minerals in sediments (clays, magnetite, Fe-oxides) have some capacitance and, therefore, exhibit a variable phase shift.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗