SEARCH · Engineering Papers

Results for “permeability models”

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 523 records · Page 29

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Coupled Hydrological-Thermal-Biogeochemical Modeling for Predicting Arctic Carbon Emissions (CH4PACE)

The Arctic, with nearly 1,700 billion metric tons of carbon stored in its soils, is experiencing significant warming trends that raise concerns about the potential release of greenhouse gases, particularly methane. This report investigates the conversion of carbon stored in Arctic soils and permafrost into methane and its subsequent atmospheric release over the coming decades. Recent data indicate a rise in global methane concentrations, highlighting the importance of understanding the mechanisms of methane generation under anaerobic conditions prevalent in thawing permafrost. Our study employs a multi-faceted approach that integrates fieldwork, benchtop experiments, and numerical simulations to explore the subsurface hydrological-thermal-biogeochemical systems influencing methane emissions. We focus on the relationship between microbial activity in decomposing organic matter under oxygen-limited conditions and the transport of methane out of these saturated wetland soils. Key hydrological parameters, including permeability, porosity, and thermal conductivity, were measured on soil samples from Fairbanks, Alaska, and used to parameterize PFLOTRAN, a coupled hydrological-thermal-reactive transport simulator. Our model enables prediction of the evolution of the active layer and its impact on methane production and release. Our findings underscore the critical need for a mechanistic understanding of methane dynamics in permafrost environments to improve predictions of methane emissions in Earth systems models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Accelerating CO2 Storage Site Characterization through a New Understanding of Favorable Formation Properties and the Impact of Core-Scale Heterogeneities

CO2 sequestration in deep geologic formations can permanently reduce atmospheric CO2 emissions and help to abate climate change. Target formations must undergo a time- and resource-intensive site evaluation process, assessing storage capacity, environmental safety, and suitability for CO2 trapping via reactive transport models based on data from a limited number of core samples. As such, simulations are often simplified and omit heterogeneities in formation properties that may be significant but are not well understood. To facilitate more rapid site assessment, this work first defines the aquifer properties of favorable storage formations through the analysis of promising and active storage sites. Data show quartz is the most prevalent formation mineral with carbonate minerals, highly reactive with injected CO2, present in over 75% of formations. Porosity and permeability data are highly clustered at 10–30% and 10–1000 mD. Field-scale reactive transport simulations are then constructed and used to analyze CO2 trapping efficiency. The models consider porosity and carbonate mineral heterogeneity as well as the impacts of typical temperature gradients. Simulated sequestration efficiencies are compared to results from a comparable homogenous model to understand the implications of aquifer non-uniformities. The results show a lower sequestration efficiency in the homogeneous model during the injection phase. During the post-injection phase, the homogenization of porosity and carbonate mineralogy results in a higher sequestration efficiency. Incorporating the temperature gradient also increases the sequestration efficiency. Importantly, the maximum deviation between the homogeneous and heterogeneous simulations at the end of the 50-year study period is only ~10%. Larger impacts may be incurred for properties outside the defined, promising ranges suggested here.

Environmental Sciences & Ecology↗

NASA NTRS1995

Fossilization Processes in Thermal Springs

To create a comparative framework for the study of ancient examples, we have been carrying out parallel studies of the microbial biosedimentology, taphonomy and geochemistry of modem and sub-Recent thermal spring deposits. One goal of the research is the development of integrated litho- and taphofacies models for siliceous and travertline sinters. Thermal springs are regarded as important environments for the origin and early evolution of life on Earth, and we seek to utilize information from the fossil record to reconstruct the evolution of high temperature ecosystems. Microbial contributions to the fabric of thermal spring sinters occur when population growth rates keep pace with, or exceed rates of inorganic precipitation, allowing for the development of continuous biofilms or mats. In siliceous thermal springs, microorganisms are typically entombed while viable. Modes of preservation reflect the balance between rates of organic matter degradation, silica precipitation and secondary infilling. Subaerial sinters are initially quite porous and permeable and at temperatures higher than about 20 C, organic materials are usually degraded prior to secondary infilling of sinter frameworks. Thus, organically-preserved microfossils are rare and fossil information consists of characteristic biofabrics formed by the encrustation and underplating of microbial mat surfaces. This probably accounts for the typically low total organic carbon values observed in thermal spring deposits. In mid-temperature, (approx. 35 - 59 C) ponds and outflows, the surface morphology of tufted Phormidium mats is preserved through mat underplating by thin siliceous: crusts. Microbial taxes lead to clumping of ceils and/or preferred filament orientations that together define higher order composite fabrics in thermal spring stromatolites (e.g. network, coniform, and palisade). At lower temperatures (less than 35 C), Calothrix mats cover shallow terracette pools forming flat carpets or pustular surfaces that produce palisade and "shrub" fabrics, respectively. At finer scales, composite fabrics are seen to consist distinctive associations of microstructures formed by the encrustation of individual cells and filaments. Composite fabrics survive the diagenetic transitions from primary opaline silica to quartz and are known from subaerial thermal spring deposits as old as Lower Carboniferous. However, fossil microorganisms tend to be rare in older deposits, and are usually preserved only where cells or sheaths have been stained by iron oxides. In subaqueous mineralizing springs at lower temperatures, early infilling leads to a more rapid and complete reduction in porosity and permeability. This process, along with the slower rates of microbial degradation at lower temperatures, creates a more favorable situation for organic matter preservation. Application of this taphonomic model to the Rhynie Chert, previously interpreted as subaerial, suggest it was probably deposited in a subaqueous spring setting at lower temperatures.

Farmer, Jack D.↗

DOE OSTI2024

Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets

This repository contains simulated outputs from superwell – a hydro-economic tool for long-term assessment of groundwater cost and supply – providing globally gridded groundwater extractable volumes and associated unit costs ($/km³) for accessible groundwater production, based on a variety of user-defined depletion and ponded depth scenarios. Key model documentation: Niazi, H., Ferencz, S. B., Graham, N. T., Yoon, J., Wild, T. B., Hejazi, M., Watson, D. J., & Vernon, C. R. (2025). Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1. Geoscientific Model Development, 18(5), 1737-1767. https://doi.org/10.5194/gmd-18-1737-2025 Find the source code of the superwell model on GitHub: https://github.com/JGCRI/superwell Repository Overview Main output: superwell_outputs.7z contains 6 files (4.5 GB) named as superwell_py_deep_all_0.*PD_0.*DL.csv. These files present superwell outputs of global groundwater extraction volumes and cost estimates on a 0.5° scale for six scenarios with different Ponded Depth (PD; 0.3 and 0.6 m) and Depletion Limit (DL; 5%, 25%, and 40% of available volume) targets over the entire pumping lifetime of a grid cell superwell_py_deep_all_0.3PD_0.25DL_sample_100.csv contains superwell outputs for 100 data points sampled to match the global inputs' distribution superwell_py_deep_all_0.3PD_0.25DL_Grid_72548.csv contains superwell output for a single grid cell concept_v5.png provides an overview of the superwell workflow Outputs Description year_number: year of pumping depletion_limit: set depletion limit (DL) as a volume fraction of total available groundwater Mappings: continent, country, gcam_basin_id, Basin_long_name, grid_id: geographic identifiers and basin information Inputs: grid_area (km²): area of the grid cell whyclass: hydrogeological classification of the aquifer permeability (m/day), porosity (%), total_thickness (m), depth_to_water (m): aquifer properties. The geo-processed input data has been published separately: https://doi.org/10.57931/2307831 Model outputs: orig_aqfr_sat_thickness (m), aqfr_sat_thickness (m): original and remaining/instantaneous saturated thickness of the aquifer hydraulic_conductivity (m/day), transmissivity (m²/day): hydraulic properties of the aquifer radius_of_influence (m), areal_extent (km²): well radius and area of influence from the center of the well number_of_wells (-): number of wells in a grid cell determined by a ratio of well area and grid area max_drawdown (m), drawdown (m), drawdown_interference (m): well and aquifer drawdown during extraction total_head (m): total lift for the groundwater (depth to water plus drawdown) total_well_length (m): total depth of wells drilled well_yield (m³/day): pumping rate or well yield power (kW), energy (kWh): power and energy required for pumping groundwater Volume Outputs: volume_produced_perwell (m³), cumulative_vol_produced_perwell (m³): production volume metrics per well volume_produced_allwells (m³), cumulative_vol_produced_allwells (m³): aggregate extraction volumes for all wells in a grid cell available_volume (m³): available groundwater in storage for the grid cell as determined by aquifer properties depleted_vol_fraction: fraction of total volume pumped over available volumes in a grid cell (same as depletion limit) Cost Outputs: well_installation_cost ($): well installation cost based on the hydrogeological complexity of the aquifer annual_capital_cost, maintenance_cost, nonenergy_cost ($): nonenergy costs energy_cost_rate ($/kWh): electricity rate energy_cost ($): energy cost of pumping groundwater total_cost_perwell ($), total_cost_allwells ($): total annual energy and non-energy cost for each and all wells in a grid cell a unit_cost ($/m³), unit_cost_per_km3 ($/km³), unit_cost_per_acreft ($/acre-ft): total cost of pumping a unit of groundwater, indicated for different spatial units Key Resources Model documentation: Niazi, H., Ferencz, S., Graham, N., Yoon, J., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024; In-prep). Long-term Hydro-economic Assessment Tool for Evaluating Global Groundwater Cost and Supply: Superwell v1. Geoscientific Model Development. Input data: Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MSD-LIVE Data repository. https://doi.org/10.57931/2307831 superwell source code: https://github.com/JGCRI/superwell Cite as Niazi, H., Ferencz, S., Yoon, J., Graham, N., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024). Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets. MSD-LIVE Data repository. https://doi.org/10.57931/2307832 Contact Reach out to Hassan Niazi or Stephen Ferencz or open an issue in the superwell repository for questions or suggestions.

Earth Systems↗

DOE OSTI2020

Probabilistic risk-based Area of Review (AoR) determination for a deep-saline carbon storage site

Regulatory oversight of a geologic carbon sequestration (GCS) project relies on iterative estimations, throughout the project lifetime, of the area where increased risks to underground sources of drinking water (USDWs) may occur due to injection of CO2. This area, referred to as Area of Review (AoR), is typically delineated by predicting the migration of fluid between the reservoir and the lowermost USDW via an open wellbore using predictions from physics-based reservoir simulators. The inherent uncertainty in input parameters used in reservoir modeling therefore affects the accuracy of determining the AoR for a project. Furthermore, the standard analytical approaches for calculating a critical pressure to delineate the risk area yield an infinite AoR for cases where the injection reservoir is overpressured relative to the USDW. A methodology is presented in this paper to better characterize the risk to USDWs while accounting for the uncertainty in reservoir modeling, with an application to a permitted GCS project with an overpressured injection formation, FutureGen 2.0. The methodology is demonstrated using the National Risk Assessment Partnership’s open-source integrated assessment model (NRAP-Open-IAM) to develop a probabilistic estimate of impact risk to USDW quality. CO2 and pressure predictions from the reservoir modeling conducted using the STOMP-CO2 simulator for the FutureGen 2.0 site are used in a NRAP-Open-IAM model with reservoir, wellbore, and aquifer components to: (1) assess the extent of potential leakage into the USDW for the predicted reservoir pressure conditions; (2) evaluate the extent of potential impact using “no-net-degradation” thresholds; and (3) account for uncertainty in reservoir permeabilities.

03 NATURAL GAS↗

NASA NTRS1973

Polarizability calculations on water, hydrogen, oxygen, and carbon dioxide

A semiclassical model of damped oscillators is used as a basis for the calculation of the dispersion of the refractive index, polarizability, and dielectric permeability in water, hydrogen, and oxygen in liquid and gaseous states, and in gaseous carbon dioxide. The absorption coefficient and the imaginary part of the refractive index are also calculated at corresponding wavelengths. A good agreement is obtained between the observed and calculated values of refractive indices, and between those of absorption coefficients in the region of absorption bands. The calculated values of oscillator strengths and damping factors are also discussed. The value of the polarizability of liquid water was about 2.8 times that of previous calculations.

Nir, S.↗

NASA NTRS2004

Aeroacoustic Analysis of a Simplified Landing Gear

A hybrid approach is used to investigate the noise generated by a simplified landing gear without small scale parts such as hydraulic lines and fasteners. The Ffowcs Williams and Hawkings equation is used to predict the noise at far-field observer locations from flow data provided by an unsteady computational fluid dynamics calculation. A simulation with 13 million grid points has been completed, and comparisons are made between calculations with different turbulence models. Results indicate that the turbulence model has a profound effect on the levels and character of the unsteadiness. Flow data on solid surfaces and a set of permeable surfaces surrounding the gear have been collected. Noise predictions using the porous surfaces appear to be contaminated by errors caused by large wake fluctuations passing through the surfaces. However, comparisons between predictions using the solid surfaces with the near-field CFD solution are in good agreement giving confidence in the far-field results.

Lockard, David P.↗

NASA NTRS1994

Permeation of membranes by the neutral form of amino acids and peptides: relevance to the origin of peptide translocation

The flux of amino acids and other nutrient solutes such as phosphate across lipid bilayers (liposomes) is 10(5) slower than facilitated inward transport across biological membranes. This suggest that primitive cells lacking highly evolved transport systems would have difficulty transporting sufficient nutrients for cell growth to occur. There are two possible ways by which early life may have overcome this difficulty: (1) The membranes of the earliest cellular life-forms may have been intrinsically more permeable to solutes; or (2) some transport mechanism may have been available to facilitate transbilayer movement of solutes essential for cell survival and growth prior to the evolution of membrane transport proteins. Translocation of neutral species represents one such mechanism. The neutral forms of amino acids modified by methylation (creating protonated weak bases) permeate membranes up to 10(10) times faster than charged forms. This increased permeability when coupled to a transmembrane pH gradient can result in significantly increased rates of net unidirectional transport. Such pH gradients can be generated in vesicles used to model protocells that preceded and were presumably ancestral to early forms of life. This transport mechanism may still play a role in some protein translocation processes (e.g. for certain signal sequences, toxins and thylakoid proteins) in vivo.

NASA Discipline Exobiology↗

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 OSTI2023

Generic Aquifer Component Training Dataset

Training dataset for the Generic Aquifer Component of NRAP-Open-IAM. Consists of the results of 62,500 STOMP-CO2 simulations of brine and CO2 leakage of varying rates and salinity into aquifers of varying thickness, depth, porosity, permeability, anisotropy and initial salinity. Described in report https://www.pnnl.gov/main/publications/external/technical_reports/PNNL-32590.pdf

Aquifer Impact Model↗

NASA NTRSDate not supplied

Computational Investigation of Oxidative Etch Pitting in FiberForm and Its Impact on Material Properties

Oxidation-driven carbon erosion does not occur uniformly but rather through the development of localized etch pits at active surface sites. These active sites form due to atomic defects on the carbon surface, making them significantly more reactive than the surrounding, non-defective areas. As a result, these sites are the first to react during ablation, leading to their removal. This process creates new defects in neighboring atoms, increasing their reactivity and causing localized carbon removal around these active sites. In this way, the highly reactive defective areas serve as nucleation points for the formation and growth of etch pits, which can have adverse effects on structural integrity of FiberForm. To better understand how these etch pits impact the material properties of carbon fiber microstructures, we have developed a new capability within the direct simulation Monte Carlo (DSMC) framework to capture the etch pit formation process. This capability, integrated into the DSMC code SPARTA (Stochastic Parallel Rarefied-gas Time-accurate Analyzer), models material removal in the presence of active sites, leading to the formation of etch pits. The current work focuses on studying the effects of these etch pits on the material properties of FiberForm, a widely used base material in thermal protection systems (TPS). The microstructure of virgin FiberForm, obtained via X-ray microtomography, is imported into SPARTA to generate the ablated geometries with etch pits. These modified microstructures are then analyzed using the Porous Microstructure Analysis (PuMA) software to compute various material properties, including elasticity, thermal conductivity, and permeability. We investigate the variation of these properties due to the complex surface topology changes caused by etch pit formation. Additionally, we compare the effects of pitting with the conventional model of shrinking fibers, traditionally used to simulate the ablation of carbon structures. Significant differences emerge between the two approaches. Consequently, this physically realistic model of material removal through etch pit formation offers improved accuracy in predicting the degradation of carbon-based TPS during oxidation. It also provides insights into other mechanisms, such as spallation, where chunks of material are removed into the flow due to etch pit growth. Ultimately, this model enhances our understanding of failure modes in these materials during ablation.

PuMA↗

DOE OSTI2022

Effects of Intrusive Confining Units on Groundwater Flow at Pahute Mesa, Nevada National Security Site

Geological, geophysical, hydrogeological, and hydrological evidence indicates that feeder dikes to lava flows play a role as confining units in the groundwater flow system underlying the Pahute Mesa area of underground nuclear testing at the Nevada National Security Site (NNSS). Of geological significance, the Pahute Mesa testing area is largely situated within the Silent Canyon caldera, which includes both tuffs and lavas that were fed by (and are connected to) intrusive rock units such as feeder dikes prior to and during volcanic eruptions (Orkild et al., 1968; Byers et al., 1976). Groundwater modeling results presented in this report support extension of “pool-dam” and “breach scenario” concepts presented in Jackson and Fenelon (2018) and Jackson et al. (2021) to low-permeability intrusive confining units (ICUs) such as feeder dikes that are expected to dissect and ring the caldera. As laterally extensive, steeply-dipping, and sub-parallel planar structures, ICUs can be expected to further compartmentalize groundwater flow in the volcanic aquifers systems at Pahute Mesa beyond the stratigraphically-defined tuff confining units (TCUs). ICUs provide an explanation for changes in water levels before and after the Cheshire test, detonated 2/14/76 (USDOE, 2015), and the unusually high horizontal gradient southwest of Cheshire. These ICU concepts are supported by groundwater modeling results presented in this report. The evidence and groundwater modeling results indicate that the downgradient extent of radionuclide contamination in some areas of Pahute Mesa could be overestimated by ignoring the confining effects of ICUs. On the other hand, contaminant extent could be underestimated by not considering transient flow effects and/or the re-equilibrated flow field caused by breaching of an ICU by an underground nuclear test.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Pore-scale CFD simulations of clay mobilization in natural porous media due to fresh water injection

The present work investigates mechanisms of permeability impairment as a result of low-salinity fluid injection into brine-saturated porous media containing dispersible clays. Here, we present a computational fluid dynamics model at the pore-scale to simulate detachment, migration and straining of fine particles in porous media. The model uses an immersed boundary method to simulate the motion of clay fines in a fluid. In addition to the hydrodynamic forces, we model the Derjaguin-Landau-Verwey-Overbeek forces (DLVO) between clay fines and grains. Our simulations show the impact of the injected fluid's salinity and velocity on the concentration of clay fines retained on the grain surface. We see clay particles dislodging from the grain surface in clusters of up to 12 particles. Our simulation results also demonstrate clogging of narrow pore spaces by the detached particles from the grain surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Impact of reservoir parameters and wellbore permeability uncertainties on CO 2 and brine leakage potential at the Shenhua CO 2 Storage Site, China

Here, this paper describes a system-level risk assessment for the Shenhua CO 2 storage site, China, using the National Risk Assessment Partnership Integrated Assessment Model for Carbon Storage (NRAP-IAM-CS). Here, we begin by determining the optimal number of Monte Carlo (MC) simulations to achieve CO 2 and brine leakage result convergence. Then, we calculate mass CO 2 and brine leakage to the atmosphere and a hypothetical shallow aquifer. Finally, we assess the geochemical impacts in the event of leakage as if there were a shallow freshwater aquifer at the Shenhua site. Simulation results show that leakage results tend to stabilize after 300 Monte Carlo simulations. When the three wells on site are assigned a permeability of 10 -11 m 2 (representing significantly leaking wells), moderate CO 2 and brine leakages occur, and the percentage of CO 2 leakage exceeds the threshold value we set based on the Intergovernmental Panel on Climate Change (IPCC). This is, however, unlikely to be the case for the Shenhua site. For all the other scenarios, the CO 2 leakage is trivial although there is still the possibility of CO 2 leakage into the groundwater aquifer exceeding the 1% threshold over 1000 years, assuming constant legacy wellbore permeability. For the significantly leaking scenario, there is a 10% probability to have a moderate (2.7 × 10 7 m 3 ) leakage-affected volume in the shallow aquifer with the pH below 6.5, and a large (1.4 × 10 8 m 3 ) volume with the total dissolved solids (TDS) above 500 ppm, hence pH and TDS may be considered for site monitoring plans. Based on the simulation results, there is a very low probability of significant CO 2 /brine leakage through the existing wells at the Shenhua CO 2 storage site.

58 GEOSCIENCES↗

DOE OSTI2026

Direction-specific enhanced diffusion of CO 2 in chiral hexagonal boron nitride nanotubes

To meet performance requirements, the next generation of gas separation membranes will need both high gas permeability and selectivity, attainable if we could coax adsorbates to minimize random Brownian motion and produce direction-specific diffusion along a desired axis. In this atomistic modeling study, we detail how direction-specific diffusion of CO 2 can be achieved in chiral hexagonal boron nitride nanotubes (hBNNTs) by means of a non-Knudsen diffusion mechanism. Our findings detail how this mechanism of diffusion is driven by interactions with the tube walls and enables the CO 2 molecules to diffuse along the nanotube’s z-axis with minimized collisions and directional changes. hBNNTs with chiral indices exhibit CO 2 diffusion rates faster than non-chiral tubes of comparable and larger diameters. Of the hBNNTs studied, a (7,3) tube appears to be ideally sized (3.7 Å radius) exhibiting CO 2 diffusion that is 3.4 times faster than diatomic N 2 . Applying this mechanism of diffusion to hypothetical sheet membranes prepared with aligned chiral (7,3) hBNNTs results in membranes with a calculated CO 2 /N 2 permselectivity of 170 and a CO 2 permeability limit of nearly 1.35 ×10 7 Barrer, readily surpassing the Robeson upper bound for CO 2 /N 2 separations.

CO2↗

DOE OSTI2020

Coupled Time-Lapse Full-Waveform Inversion for Subsurface Flow Problems Using Intrusive Automatic Differentiation

We describe a novel framework for estimating subsurface properties, such as rock permeability and porosity, from time-lapse observed seismic data by coupling full-waveform inversion (FWI), subsurface flow processes, and rock physics models. For the inverse modeling, we handle the back propagation of gradients by an intrusive automatic differentiation strategy that offers three levels of user control: (1) At the wave physics level, we adopted the discrete adjoint method in order to use our existing high-performance FWI code; (2) at the rock physics level, we used built-in automatic differentiation operators from the TensorFlow backend; (3) at the flow physics level, we implemented customized partial differential equation (PDE) operators for the multiphase flow equations. The three-level coupled inversion strategy strikes a good balance between computational efficiency and programming efforts, and when the gradients are chained together, it constitutes a coupled inverse system. Our numerical experiments demonstrate that the three-level coupled inverse problem is superior in terms of accuracy to a traditional decoupled inversion strategy. Additionally, our method is able to simultaneously invert for parameters in empirical relationships such as the rock physics models. Our proposed inverted model can be used for reservoir performance prediction and reservoir management/optimization purposes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Autonomous Fracture Conductivity Using Expandable Proppants in Enhanced Geothermal Systems

Summary Early thermal breakthrough in enhanced geothermal systems (EGS) due to the presence of preferential flow channels is a major challenge that endangers efficient and economic heat extraction in such systems. Previous studies mainly focused on adjusting circulation rates of the working fluid, which still leaves significant amounts of untapped heat behind. Currently, there is a lack of technologies for altering flow distribution within the fracture network to achieve uniform heat sweeping in the reservoir. This work presents a novel concept for making proppants to autonomously control fracture conductivity based on the surrounding temperature. Here, proppants with negative thermal expansion coefficients have demonstrated the capability for appropriate fracture conductivity adjustment as a function of temperature to achieve uniform flow across the fracture network. Particle-particle interactions governing such functions are explicitly modeled, and then the Lattice Boltzmann methods (LBM) is used to determine the potential impact of closure stress and temperature changes on the permeability of the proposed proppant packs. Microscale analyses are further used to determine the required material properties to achieve a certain improvement in the permeability of the proppant pack. Our analyses show an enhancement in permeability and the associated fracture conductivity by half of their initial values. Field-scale analysis further confirms the effectiveness of the proposed concept as 31.4% more heat can be extracted from EGS over 50 years of production when the proposed proppants are used. Such novel proppants may effectively delay thermal breakthrough, sweep heat from larger rock volumes, and elongate the life span of the EGS project.

Engineering↗

DOE OSTI2023

A Novel Approach to Map Permeability Using Passive Seismic Emission Tomography

Newly acquired magnetotelluric data and passive seismic data collected with tightly spaced geophone arrays are combined with historic drilling, active seismic, and potential fields data to generate 3-D permeability maps. A cooperative inversion methodology has been developed using active seismic, magnetotelluric, and gravity data in order to produce more robust velocity models for passive seismic data processing without requiring expensive 3-D active seismic surveys. The cooperative inversion estimates velocities from other geophysical data where no prior seismic velocity information is available at two geothermal sites in Nevada: San Emidio and Crescent Valley.

15 GEOTHERMAL ENERGY↗