Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Exploring the Mass-loss Histories of the Red Supergiants

We report mid- to far-infrared imaging and photometry from 7 to 37 μm with SOFIA/FORCAST and 2 μm adaptive optics imaging with LBTI/LMIRCam of a large sample of red supergiants (RSGs) in four Galactic clusters: RSGC1, RSGC2 = Stephenson 2, RSGC3, and NGC 7419. The RSGs in these clusters cover their expected range in luminosity and initial mass from ≈9 to more than 25 M {sub ⊙}. The population includes examples of very late-type RSGs such as MY Cep, which may be near the end of the RSG stage, high-mass-losing maser sources, yellow hypergiants, and post-RSG candidates. Many of the stars and almost all of the most luminous have spectral energy distributions (SEDs) with extended infrared excess radiation at the longest wavelengths. To best model their SEDs, we use the DUSTY code with a variable radial density distribution function to estimate their mass-loss rates. Our M-dot –luminosity relation for 42 RSGs basically follows the classical de Jager curve, but at luminosities below 10{sup 5} L{sub ⊙}, we find a significant population of RSGs with M-dot below the de Jager relation. At luminosities above 10{sup 5} L {sub ⊙}, there is a rapid transition to higher mass-loss rates that approximates and overlaps the de Jager curve. We recommend that instead of using a linear relation or single curve, the empirical M-dot –luminosity relation is better represented by a broad band. Interestingly, the transition to much higher M-dot at about 10{sup 5} L {sub ⊙} corresponds approximately to an initial mass of 18–20 M {sub ⊙}, which is close to the upper limit for RSGs becoming Type II supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

GW190814: Gravitational Waves from the Coalescence of a 23 Solar Mass Black Hole with a 2.6 Solar Mass Compact Object

We report the observation of a compact binary coalescence involving a 22.2-24.3 M ⊙ black hole and a compact object with a mass of 2.50-2.67 M ⊙ (all measurements quoted at the 90% credible level). The gravitational-wave signal, GW190814, was observed during LIGO's and Virgo's third observing run on 2019 August 14 at 21:10:39 UTC and has a signal-to-noise ratio of 25 in the three-detector network. The source was localized to 18.5 deg 2 at a distance of 241 +41-41 Mpc; no electromagnetic counterpart has been confirmed to date. The source has the most unequal mass ratio yet measured with gravitational waves , 0.112+0.0090.008, , and its secondary component is either the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system. The dimensionless spin of the primary black hole is tightly constrained to ≤0.07. Tests of general relativity reveal no measurable deviations from the theory, and its prediction of higher-multipole emission is confirmed at high confidence. We estimate a merger rate density of 1-23 Gpc -3 yr -1 for the new class of binary coalescence sources that GW190814 represents. Astrophysical models predict that binaries with mass ratios similar to this event can form through several channels, but are unlikely to have formed in globular clusters. However, the combination of mass ratio, component masses, and the inferred merger rate for this event challenges all current models of the formation and mass distribution of compact-object binaries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Soil physical and thermal properties from soil samples at multiple locations at Teller Road Site, Seward Peninsula, Alaska, 2019

This dataset contains results of the laboratory analysis of soil samples collected across the watershed located at miles #27 along the Teller road, near Nome, AK. The dataset includes 124 soil samples retrieved from 52 locations during the first week of August 2019. The soil samples were collected in the top 1.2 m of soil with volume known samples, analyzed with a thermal analyser to infer thermal parameters (thermal conductivity, volumetric heat capacity and thermal diffusivity), and then processed in a lab for estimation of their bulk density, water content, carbon and nitrogen content, soil texture, and elemental fractions. This dataset has one csv file and one pdf file.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Unbinned extraction of $γ$ from $B\to DK$ with normalizing flows

We introduce an unbinned method for extracting the CKM angle $γ$ from the decay chain $B^\pm \to (D \to K_S π^+ π^-) K^\pm$ using normalizing flows (NFs). The NFs, trained on $D$ decay data, learn a faithful continuous representation of the amplitude and strong phase variation over the $D\to K_Sπ^+π^-$ Dalitz plot whose fidelity improves with increased data sample sizes. With this input, the $B$ decay data can be used to extract the parameters $r_B$, $δ_B$, and $γ$. We test the method on Monte Carlo generated data, where it successfully recovers the injected value of $γ$ within uncertainties. The present implementation propagates statistical uncertainties from finite training data via an ensemble of independently trained flows, and does not attempt to capture the effects of systematic experimental errors. We explore two versions of the method that differ in how the trigonometric constraint on phase variation is encoded, and comment on the possible extension to Bayesian NFs, which would provide direct uncertainty estimates on the learned densities without requiring ensemble training.

Grossman, Yuval [Cornell U., LEPP]↗
DOE OSTI2026

Computing the Critical Temperature of the Affine-Transformed $D=3$ Ising Model Using Masked Autoregressive Flow

The simple Ising model provides a rich environment to build and study lattice field theories. As part of an ongoing project to construct a conformal field theory (CFT) on an arbitrarily curved manifold, in this work we develop methods to measure the critical temperature $β_c$ of the affine-transformed Ising model on the face-centered cubic (FCC) lattice. The main challenge in this endeavor is finding a computationally efficient and accurate method of interpolating and extrapolating Monte Carlo observables with respect to coupling coefficients and temperature. Herein, we compare two such methods. A traditional statistical approach uses the multiple histogram (MH) method, while a newer machine learning approach uses a masked autoregressive flow (MAF) to estimate the underlying probability density function of a set of observables. While the MH method is specifically designed to interpolate and extrapolate Monte Carlo observables, we find that MAF is a viable alternative for measuring $β_c$ with a computational cost that scales more favorably. Furthermore, we comment on additional advantages of MAF relevant to our work, such as extrapolating in system volume.

Svenson, Kai [Texas U.]↗
DOE OSTI2021

Experimental 3-dimensional tracking of the dynamics of a single electron in the Fermilab Integrable Optics Test Accelerator (IOTA)

We present the results of experimental studies on the transverse and longitudinal dynamics of a single electron in the IOTA storage ring. IOTA is a flexible machine dedicated to beam physics experiments with electrons and protons. A method was developed to reliably inject and circulate a controlled number of electrons in the ring. A key beam diagnostic system is the set of sensitive high-resolution digital cameras for the detection of synchrotron light emitted by the electrons. With 60–130 electrons in the machine, we measured beam lifetime and derived an absolute calibration of the optical system. At exposure times of 0.5 s, the cameras were sensitive to individual electrons. Camera images were used to reconstruct the time evolution of oscillation amplitudes of a single electron in all 3 degrees of freedom. The evolution of amplitudes directly showed the interplay between synchrotron-radiation damping, quantum excitations, and scattering with the residual gas. From the distribution of measured single-electron oscillation amplitudes, we deduced transverse emittances, momentum spread, damping times, and beam energy. Estimates of residual-gas density and composition were calculated from the measured distributions of vertical scattering angles. Combining scattering and lifetime data, we also provide an estimate of the aperture of the ring. To our knowledge, this is the first time that the dynamics of a single electron are tracked in all three dimensions with digital cameras in a storage ring.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Aboveground biomass density models for NASA’s Global Ecosystem Dynamics Investigation (GEDI) lidar mission

NASA's Global Ecosystem Dynamics Investigation (GEDI) is collecting spaceborne full waveform lidar data with a primary science goal of producing accurate estimates of forest aboveground biomass density (AGBD). This paper presents the development of the models used to create GEDI's footprint-level (~25 m) AGBD (GEDI04_A) product, including a description of the datasets used and the procedure for final model selection. The data used to fit our models are from a compilation of globally distributed spatially and temporally coincident field and airborne lidar datasets, whereby we simulated GEDI-like waveforms from airborne lidar to build a calibration database. We used this database to expand the geographic extent of past waveform lidar studies, and divided the globe into four broad strata by Plant Functional Type (PFT) and six geographic regions. GEDI's waveform-to-biomass models take the form of parametric Ordinary Least Squares (OLS) models with simulated Relative Height (RH) metrics as predictor variables. From an exhaustive set of candidate models, we selected the best input predictor variables, and data transformations for each geographic stratum in the GEDI domain to produce a set of comprehensive predictive footprint-level models. We found that model selection frequently favored combinations of RH metrics at the 98th, 90th, 50th, and 10th height above ground-level percentiles (RH98, RH90, RH50, and RH10, respectively), but that inclusion of lower RH metrics (e.g. RH10) did not markedly improve model performance. Second, forced inclusion of RH98 in all models was important and did not degrade model performance, and the best performing models were parsimonious, typically having only 1-3 predictors. Third, stratification by geographic domain (PFT, geographic region) improved model performance in comparison to global models without stratification. Fourth, for the vast majority of strata, the best performing models were fit using square root transformation of field AGBD and/or height metrics. There was considerable variability in model performance across geographic strata, and areas with sparse training data and/or high AGBD values had the poorest performance. These models are used to produce global predictions of AGBD, but will be improved in the future as more and better training data become available.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Spatial and temporal overlap between hatchery- and natural-origin steelhead and Chinook Salmon during spawning in the Klickitat River, Washington, USA

Abstract Objective A goal of many segregated salmonid hatchery programs is to minimize potential interbreeding between hatchery- and natural-origin fish. Our objective was to assess this on the Klickitat River, Washington, USA. Methods We used radiotelemetry to evaluate spatiotemporal spawning overlap between hatchery- and natural-origin steelhead Oncorhynchus mykiss and spring Chinook Salmon O. tshawytscha. We estimated percentages of tagged fish that spawned naturally in the Klickitat River subbasin, emigrated from the Klickitat River, or died before spawning. A kernel density analysis was used to estimate probability of spatiotemporal overlap between hatchery- and natural-origin spawners. Result For steelhead, 12% of hatchery-origin and 50% of natural-origin fish spawned naturally. For spring Chinook Salmon, 18% of hatchery-origin and 44% of natural-origin fish spawned naturally. Tag loss may result in underestimates in these percentages. Most hatchery-origin steelhead (90%) spawned downstream of river kilometer (rkm) 32, and 75% spawned from November to mid-March. The majority of natural-origin steelhead (64%) spawned upstream of rkm 32, and 75% spawned from mid-March to late May. Spawn timing of hatchery-origin Chinook Salmon (early August to mid-September) overlapped with that of natural-origin Chinook Salmon (late July to late September), and fish of both origins spawned in the same 30-km reach of the river. We estimated the percentage of hatchery-origin spawners (pHOS) on the natural spawning grounds to be 12% for steelhead and 40% for spring Chinook Salmon across all study years. For steelhead, we estimated the overlap probability to be 25% (95% CI = 22.5–28%). For spring Chinook Salmon, tight spatial clustering of hatchery-origin fish resulted in a lower overlap estimate of 21% (13–31%). Conclusion We suggest adjusting pHOS estimates using these overlap estimates or similar spatiotemporal data on actual spawner proximity and possible interactions, and that these types of analyses be used in conjunction with gene flow analysis to accurately evaluate effects of individual hatchery programs.

Zendt, Joseph S.↗
DOE OSTI2025

Spatially and temporally resolved plasma parameter estimations of laser heated MagLIF relevant gas pipes at NIF

The ability to control laser pre-heat is an integral part of the inertial confinement fusion concept known as Magnetized Liner Inertial Fusion. This process is studied at the National Ignition Facility (NIF) where 4 of the 192 laser beams are propagated through a 1-cm long gas cell where they deposit >20 kJ of energy into the gaseous fuel via inverse bremsstrahlung absorption. This process ionizes the gas, producing a plasma that follows behind the laser front and expands over the radius of the cell. Emission from this plasma, as viewed by a gated x-ray detector, can be used to build spatially and temporally resolved estimations of the pre-heat plasma's density and temperature profiles. This can then be used to estimate the plasma pressure, internal energy, and radiation losses. Estimations show the evolution of the plasma in magnetized and unmagnetized gas cells filled with ambient temperature neopentane (C5H12) +1% Ar, as well as unmagnetized cryogenically cooled (32 K) deuterium +1% Ne filled targets. This analysis shows the effects of initial gas-fill density, composition, and axial magnetization on the time-dependent plasma parameters. Previously, these parameters at the NIF had not been experimentally characterized, and these estimations provided a potential new means of testing radiation magneto-hydrodynamic predictive capability models. Results in unmagnetized targets have strong agreement with simulations. However, in targets with a 19 T applied axial magnetic field, this method yields electron temperatures up to 100% hotter than those predicted by HYDRA codes.

Bremsstrahlung↗
DOE OSTI2023

Bulk viscosity from Urca processes: n p e μ matter in the neutrino-transparent regime

We study the bulk viscosity of moderately hot and dense, neutrino-transparent relativistic npeμ matter arising from weak-interaction direct Urca processes. This work parallels our recent study of the bulk viscosity of npeμ matter with a trapped neutrino component. The nuclear matter is modeled in a relativistic density functional approach with two different parametrizations—DDME2 (which does not allow for the low-temperature direct-Urca process at any density) and NL3 (which allows for low-temperature direct-Urca process above a low-density threshold). Here, we compute the equilibration rates of Urca processes of neutron decay and lepton capture, as well as the rate of the muon decay, and find that the muon decay process is subdominant to the Urca processes at temperatures T ≥ 3 MeV in the case of DDME2 model and T ≥ 1 MeV in the case of NL3 model. Thus, the Urca-process-driven bulk viscosity is computed with the assumption that pure leptonic reactions are frozen. As a result the electronic and muonic Urca channels contribute to the bulk viscosity independently and at certain densities the bulk viscosity of npeμ matter shows a double-peak structure as a function of temperature instead of the standard one-peak (resonant) form. In the final step, we estimate the damping timescales of density oscillations by the bulk viscosity. We find that, e.g., at a typical oscillation frequency f = 1~kHz, the damping of oscillation is most efficient at temperatures 3 ≤ T ≤ 5~MeV and densities n B ≤ 2n 0 where they can affect the evolution of the post-merger object.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Point Cloud Based Mapping of Understory Shrub Fuel Distribution, Estimation of Fuel Consumption and Relationship to Pyrolysis Gas Emissions on Experimental Prescribed Burns

Forest fires spread via production and combustion of pyrolysis gases in the understory. The goal of the present paper is to understand the spatial location, distribution, and fraction (relative to the overstory) of understory plants, in this case, sparkleberry shrub, namely its degree of understory consumption upon burn, and to search for correlations between the degree of shrub consumption to the composition of emitted pyrolysis gases. Data were collected in situ at seven small experimental prescribed burns at Ft. Jackson, an army base in South Carolina, USA. Using airborne laser scanning (ALS) to map overstory tree crowns and terrestrial laser scanning (TLS) to characterize understory shrub fuel density, both pre- and postburn estimates of sparkleberry coverage were obtained. Sparkleberry clump polygons were manually digitized from a UAV-derived orthoimage of the understory and intersected with the TLS point cloud-derived rasters of pre- and postburn shrub fuel bulk density; these were compared in relation to overstory crown cover as well as to ground truth. Shrub fuel consumption was estimated from the digitized images; sparkleberry clump distributions were generally found to not correlate well to the overstory tree crowns, suggesting it is shade-tolerant. Moreover, no relationship was found between the magnitude of the fuel consumption and the chemical composition of pyrolysis gases, even though mixing ratios of 25 individual gases were measured.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Construction of the damped Ly⁢𝛼 absorber catalog for DESI DR2 Ly⁢𝛼 BAO

We present the Damped Ly⁢𝛼 Toolkit for automated detection and characterization of damped Ly⁢𝛼 absorbers (DLAs) in quasar spectra. Our method uses quasar spectral templates with and without absorption from intervening DLAs to reconstruct observed quasar forest regions. The best-fitting model determines whether a DLA is present while estimating the redshift and HI column density. With an optimized quality cut on detection significance (Δ⁢𝜒$^{2}_{𝑟}$ >0.03), the technique achieves an estimated 80% purity and 79% completeness when evaluated on simulated spectra with S/N>2 that are free of broad absorption lines (BALs). We provide a catalog containing candidate DLAs from the DLA Toolkit detected in DESI DR1 quasar spectra, of which 21 719 were found in S/N>2 spectra with predicted log 10 ⁡(𝑁 𝙷𝙸 )>20.3 and detection significance Δ⁢𝜒$^{2}_{𝑟}$ >0.03. We compare the Damped Ly⁢𝛼 Toolkit to two alternative DLA finders based on a convolutional neural network and Gaussian process models. We present a strategy for combining these three techniques to produce a high-fidelity DLA catalog from DESI DR2 for the Ly⁢𝛼 forest baryon acoustic oscillation measurement. The combined catalog contains 41 152 candidate DLAs with log 10 ⁡(𝑁 𝙷𝙸 )>20.3 from quasar spectra with S/N>2. We estimate this sample to be approximately 85% pure and 79% complete when BAL quasars are excluded.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Assessment of Long-Term Degradation of Adsorbents for Direct Air Capture by Ozonolysis

Porous adsorbents are a promising class of materials for the direct air capture of CO 2 (DAC). Practical implementation of adsorption-based DAC requires adsorbents that can be used for thousands of adsorption–desorption cycles without significant degradation. We examined the potential degradation of adsorbents by a mechanism that appears to have not been considered previously, namely, ozonolysis by trace levels of ozone from ambient air. We focused on amine-appended metal–organic frameworks, specifically amine-functionalized Mg 2 (dobpdc), as a representative DAC adsorbent. Estimates based on the number of amine sites in these adsorbents and the ozone concentration in air suggest that degradation by ozone may be relevant over thousands of adsorption–desorption cycles if reactions with adsorbed ozone are fast. We used density functional theory calculations to estimate reaction rates for amine groups and carbon–carbon double bonds in amine-functionalized Mg 2 (dobpdc).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Bulk Viscosity of Relativistic npeμ Matter in Neutron-Star Mergers

We discuss the bulk viscosity of hot and dense npeμ matter arising from weak-interaction direct Urca processes. We consider two regimes of interest: (a) the neutrino-transparent regime with T ≤ T tr (T tr ≃ 5 ÷ 10 MeV is the neutrino-trapping temperature); and (b) the neutrino-trapped regime with T ≥ T tr . Nuclear matter is modeled in relativistic density functional approach with density-dependent parametrization DDME2. The maximum of the bulk viscosity is achieved at temperatures T ≃ 5 ÷ 6 MeV in the neutrino-transparent regime, then it drops rapidly at higher temperatures where neutrino-trapping occurs. As an astrophysical application, we estimate the damping timescales of density oscillations by the bulk viscosity in neutron star mergers and find that, e.g., at the oscillation frequency f = 10 kHz, the damping will be very efficient at temperatures 4 ≤ T ≤ 7 MeV where the bulk viscosity might affect the evolution of the post-merger object.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

High-compression implosions based on high density carbon ablator using modified drive and capsule dopant profiles

Laser-driven inertial fusion experiments have, for the first time, achieved a target gain greater than unity in a laboratory setting [Abu-Shawareb et al., Phys. Rev. Lett. 132, 065102 (2024)]. Despite this breakthrough, the burn-up fraction remains limited to about one-fourth of ideal estimates due to insufficient areal density, highlighting the potential for greater gains through enhanced compression. In our previous work, we demonstrated record-high compression of stagnated fuel in indirectly driven implosions using high-density carbon ablators. This was achieved by combining a continuous ramped pulse drive with a modified ablator dopant profile, which reduced mixing at the fuel–ablator interface and improved stability [Tommasini et al., Phys. Rev. Res. 5, L042034 (2023)]. Based on this foundation, the study presented here investigates the limits of compression achievable by combining the continuous ramped pulse drive with different dopant profiles to further minimize unstable interfaces and gradient discontinuities, thereby reducing fuel–ablator mixing. Our results demonstrate that the continuous ramped pulse consistently outperforms designs based on 3-shock drive pulses across all ablator profiles studied, with compression showing only a relatively modest dependence on dopant configurations that reduce the number of interfaces or eliminate discontinuities in the dopant gradient profile. Sub-scale experiments using the continuous ramped pulse achieved compression levels exceeding those of full-scale “HyE” implosions [Kritcher et al., Phys. Plasmas 28, 072706 (2021)] at similar adiabat, anticipating significant performance gains with increased scale, as supported by models and simulations. These findings underscore the critical role of the continuous ramped pulse in reducing mix and achieving improved compression. They also provide a foundation for future large-scale experiments to test the continuous ramped pulse design on deuterium–tritium fuel in the burn-wave propagation regime, leveraging the most effective combinations of continuous ramped pulse and dopant profiles identified in this study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Data-driven offshore CO 2 saline storage assessment methodology

The world produces approximately 50 billion tonnes of greenhouse gases annually. This is measured in CO 2 -equivalents, and geologic CO 2 storage has the potential to advance decarbonization and mitigate greenhouse gas emissions. New technologies to assess offshore carbon storage are needed to address resource, regulatory, and commercial needs. Although most efforts to assess storage resources focus on onshore criteria, offshore reservoirs offer significant storage potential and distinct development challenges. Potential advantages of offshore carbon storage include being further from human population centers and less potential to interact with groundwater. The U.S. Department of Energy's method for evaluating storage capacity in non-oil-bearing saline reservoirs has been enhanced to support assessments for offshore environments in the Offshore CO 2 Saline Storage methodology (OCSS). This methodology applies data-driven capabilities to estimate saline storage capacity while accounting for features specific to offshore reservoirs. Features include changes in CO 2 density and sedimentary differences that impact estimates of permanence and capacity. The Offshore CO 2 Saline Storage Calculator mechanizes OCSS to estimate storage capacity. This paper presents the methodology and estimates for 18 geologic domains in the Gulf of Mexico. Potential storage distributions, sensitivity analyses, and the incorporation of spatial data and tools to support safe site selection are also discussed.

58 GEOSCIENCES↗
DOE OSTI2024

Electrical Responses of Modified Mineral Surfaces as Observed With Spectral Induced Polarization and Atomic Force Microscopy

Abstract Atomic force microscopy (AFM) and spectral induced polarization (SIP) are widely used to investigate the electrical properties of mineral surfaces at vastly different scales of measurement. We compare AFM and SIP measurements made on two different materials (glass beads and silica gel) subjected to etching, deposition of iron oxide particles, and inclusion of calcite grains. We found that the treatments produced qualitatively consistent behaviors in the AFM and SIP data. Direct AFM measurements of surface charge density for silica and calcite surfaces were quantitatively compared to values estimated from the SIP results using a grain polarization model. No statistically significant difference (at a 95% confidence level) was found between the surface charge density of silica estimated by AFM (2.3 ± 6.6 mC/m 2 for glass beads and 1.6 ± 0.1 mC/m 2 for silica gel) versus SIP (5.4 ± 4.4 mC/m 2 for glass beads and 1.6 ± 0.5 mC/m 2 for silica gel). The surface charge density for calcite determined by AFM (43.5 ± 12.9 mC/m 2 ) was approximately 19 times higher than that found for silica. While the charge density of calcite surfaces determined by SIP was also generally higher than that found for silica, different treatments produced significantly different values between 4.7 and 258 mC/m 2 (with a maximum 95% CI of ±8.7 mC/m 2 ). Several possible explanations exist for the range of the observed SIP measurements, including aging of the calcite surfaces. Overall, this study suggests the potential for the complementary use of AFM and SIP measurements to constrain future investigations of polarization mechanisms in porous media.

Geochemistry & Geophysics↗