Search NASA⌕ Search

SEARCH · Search NASA

Results for “Permeability”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins

This repository of global hydrogeologic datasets contains aquifer properties on 0.5° scale, including depth to groundwater (Fan et al., 2013), aquifer thickness (de Graaf et al., 2015), WHYMap aquifer classes (Richts et al., 2011), recharge (Döll and Fiedler, 2008; Gleeson et al., 2016), lakes (Messager et al., 2016), porosity and permeability (Gleeson et al., 2014), digitized and geo-processed from their respective sources. Globally gridded aquifer properties could be used independently to estimate global groundwater availability or used as critical inputs to the superwell model to simulate groundwater extraction and provide estimates of pumped volumes and unit costs under user-specific scenarios. Key resources related to this data are: 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 superwell model repository which uses this data to simulate groundwater extraction and provides estimates of the global extractable volumes and unit-costs ($/km3) of accessible groundwater production under user-specified extraction scenarios. Repository Overview Main output: aquifer_properties_rec.csv contains all processed outputs, including aquifer properties like porosity, permeability, recharge, lake areas, aquifer thickness, and depth to groundwater. shapefiles.zip: contains all digitized GIS databases and shapefile for all aquifer properties prep_inputs.R and prep_inputs_recharge_lakes.R: R scripts that process the shapefiles to produce the aquifer_properties_rec.csv file plot_inputs.R: R script for plotting the maps and conducting preliminary analysis on the available groundwater volume basin_to_country_mapping.csv, basin_country_region_mapping.csv and continent_county_mapping.csv provide the mapping between continents, 32 energy-economic macro regions, countries, and water basins for post-processing aquifer_properties_rec.csv Maps: Each map visualizes the spatial distribution of one of the aquifer properties across the globe map_in_Porosity.png map_in_Permeability.png map_in_Aquifer_thickness.png map_in_Depth_to_water.png map_in_Recharge.png map_in_Grid_area_km.png map_in_Lake_area_km.png map_in_WHYClass.png Sample inputs sample_inputs.py: this script samples inputs from the aquifer_properties_rec dataset, ensuring the sampled and original inputs maintain the same distributions sampled_data_100.csv contains 100 sampled data points and sampled_data_100.png compares their distributions Dataset Overview The main outputs are consolidated in a comprehensive aquifer_properties_rec.csv file and include the following fields: GridCellID: Unique identifier for each (roughly 0.5°) grid cell Continent: Continent name Country: Country name GCAM_basin_ID: Identifier for GCAM hydrologic basin Basin_long_name: Full name of the basin WHYClass: Hydrogeologic classification based on WHYMap aquifer classes (Richts et al., 2011) Porosity: Soil porosity (%) (Gleeson et al., 2014) Permeability: Soil permeability (in square meters; Gleeson et al., 2014) Aquifer_thickness: Thickness of the aquifer (in meters; de Graaf et al., 2015) Depth_to_water: Depth to groundwater (in meters; Fan et al., 2013) Recharge: long-term annual averaged recharge rates (in m/yr; Döll and Fiedler, 2008; Gleeson et al., 2016) Grid_area: Area of the grid cell (in square meters) Lakes_area: Area of inland lakes (in square meters; Messager et al., 2016) Key References The datasets are digitized versions of global hydrogeologic properties from the following key literature sources: Depth to Groundwater: Fan, Y., Li, H., & Miguez-Macho, G. (2013). Global Patterns of Groundwater Table Depth. Science, 339(6122), 940-943. https://doi.org/10.1126/science.1229881 Aquifer Thickness: de Graaf, I. E. M., Sutanudjaja, E. H., van Beek, L. P. H., & Bierkens, M. F. P. (2015). A high-resolution global-scale groundwater model. Hydrol. Earth Syst. Sci., 19(2), 823-837. https://doi.org/10.5194/hess-19-823-2015 Porosity and Permeability: Gleeson, T., Moosdorf, N., Hartmann, J., & van Beek, L. P. H. (2014). A glimpse beneath earth's surface: GLobal HYdrogeology MaPS (GLHYMPS) of permeability and porosity. Geophysical Research Letters, 41(11), 3891-3898. https://doi.org/10.1002/2014GL059856 Aquifer classes: Richts, A., Struckmeier, W. F., & Zaepke, M. (2011). WHYMAP and the Groundwater Resources Map of the World 1:25,000,000. In J. A. A. Jones (Ed.), Sustaining Groundwater Resources: A Critical Element in the Global Water Crisis (pp. 159-173). Springer Netherlands. https://doi.org/10.1007/978-90-481-3426-7_10 Recharge: Döll, P., & Fiedler, K. (2008). Global-scale modeling of groundwater recharge. Hydrol. Earth Syst. Sci., 12(3), 863-885. https://doi.org/10.5194/hess-12-863-2008; Gleeson, T., Befus, K. M., Jasechko, S., Luijendijk, E., & Cardenas, M. B. (2016). The global volume and distribution of modern groundwater. Nature Geoscience, 9(2), 161-167. https://doi.org/10.1038/ngeo2590 Inland Lakes: Messager, M. L., Lehner, B., Grill, G., Nedeva, I., & Schmitt, O. (2016). Estimating the volume and age of water stored in global lakes using a geo-statistical approach. Nature Communications, 7(1), 13603. https://doi.org/10.1038/ncomms13603 Cite as 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. MultiSector Dynamics-Living, Intuitive, Value-adding, Environment. https://doi.org/10.57931/2484226 Contact Reach out to Hassan Niazi or open an issue in the superwell repository for questions or suggestions.

aquifer thickness↗
DOE OSTI2024

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1981

Porous tooling process for manufacture of graphite/polyimide composites

A porous tooling system was selected for the processing of Graphite/PMR-15 Polyimide laminates in thickness up to 3.2 mm. (0.125 inch). This tool system must have a reasonable strength, permeability dimensional stability, and thermal conductivity to accomplish curing at 600 F and 200 psi and 200 psi autoclave temperature and pressure. A permeability measuring apparatus was constructed and permeability vs. casting water level determined to produce tools at three different permeability levels. On these tools, laminates of 5, 11, and 22 plies (.027, .060, and 0.121 inch) were produced and evaluated by ultrasonic, mechanical, and thermal tests to determine the effect of the tool permeability on the cured laminates. All tools produced acceptable laminates at 5 and 11 plies but only the highest permeability produced acceptable clear ultrasonic C-Scans. Recommendations are made for future investigations of design geometry, and strengthening techniques for porous ceramic tooling.

Smiser, L. W.↗
NASA NTRS2003

Intracellular localization of a group II chaperonin indicates a membrane-related function

Chaperonins are protein complexes that are believed to function as part of a protein folding system in the cytoplasm of the cell. We observed, however, that the group II chaperonins known as rosettasomes in the hyperthermophilic archaeon Sulfolobus shibatae, are not cytoplasmic but membrane associated. This association was observed in cultures grown at 60 degrees C and 76 degrees C or heat-shocked at 85 degrees C by using immunofluorescence microscopy and in thick sections of rapidly frozen cells grown at 76 degrees C by using immunogold electron microscopy. We observed that increased abundance of rosettasomes after heat shock correlated with decreased membrane permeability at lethal temperature (92 degrees C). This change in permeability was not seen in cells heat-shocked in the presence of the amino acid analogue azetidine 2-carboxylic acid, indicating functional protein synthesis influences permeability. Azetidine experiments also indicated that observed heat-induced changes in lipid composition in S. shibatae could not account for changes in membrane permeability. Rosettasomes purified from cultures grown at 60 degrees C and 76 degrees C or heat-shocked at 85 degrees C bind to liposomes made from either the bipolar tetraether lipids of Sulfolobus or a variety of artificial lipid mixtures. The presence of rosettasomes did not significantly change the transition temperature of liposomes, as indicated by differential scanning calorimetry, or the proton permeability of liposomes, as indicated by pyranine fluorescence. We propose that these group II chaperonins function as a structural element in the natural membrane based on their intracellular location, the correlation between their functional abundance and membrane permeability, and their potential distribution on the membrane surface.

Cell Membrane/physiology/ultrastructure↗
NASA NTRS2009

Amorphous Silk Fibroin Membranes for Separation of CO2

Amorphous silk fibroin has shown promise as a polymeric material derivable from natural sources for making membranes for use in removing CO2 from mixed-gas streams. For most applications of silk fibroin, for purposes other than gas separation, this material is used in its highly crystalline, nearly natural form because this form has uncommonly high tensile strength. However, the crystalline phase of silk fibroin is impermeable, making it necessary to convert the material to amorphous form to obtain the high permeability needed for gas separation. Accordingly, one aspect of the present development is a process for generating amorphous silk fibroin by treating native silk fibroin in an aqueous methanol/salt solution. The resulting material remains self-standing and can be prepared as thin film suitable for permeation testing. The permeability of this material by pure CO2 has been found to be highly improved, and its mixed-gas permeability has been found to exceed the mixed-gas permeabilities of several ultrahigh-CO2-permeable synthetic polymers. Only one of the synthetic polymers poly(trimethylsilylpropyne) [PTMSP] may be more highly permeable by CO2. PTMSP becomes unstable with time, whereas amorphous silk should not, although at the time of this reporting this has not been conclusively proven.

Aberg, Christopher M.↗
DOE OSTI2026

Evaluating the Transport Mechanism for Lignocellulosic Inhibitors in Zymomonas mobilis

The conversion of lignocellulosic biomass to biofuels and bioproducts is limited by small inhibitory molecules generated during the pretreatment process. The diverse compounds within the cell wall hydrolysate including aldehydes, carboxylic acids, phenolics, and alcohols inhibit microbial growth by disrupting microbial membranes. Although Zymomonas mobilis combines high ethanol productivity with hopanoid-rich membranes that confer solvent tolerance, hydrolysate compounds still inhibit microbial growth. One proposed mitigation strategy is to pump these molecules out through active transport processes, reducing direct toxicity. However, if passive permeation rates are high, active transport would create a futile cycle where the exported molecules diffuse back in, creating a net drag on fitness. To investigate this hypothesis, we measure the passive permeation for thirty-three lignocellulose derived compounds using atomistic molecular dynamics simulations with enhanced sampling. Permeability coefficients determined spanned more than eight orders of magnitude driven largely by the balance between hydrophobicity and polarity. Hydrophobic inhibitors such as aromatic acids, aldehydes, and short chain alcohols readily partitioned into the membrane and showed high permeability consistent with rapid passive influx and strong toxicity. Most polar molecules have lower but still appreciable permeation rates, with only the most polar sugars having low permeabilities. By integrating these predictions together with hydrolysate composition information, we find that most of these molecules are so permeable that futile cycles where passive permeation overwhelms efflux is highly likely. These findings provide a predictive framework for strain engineering, highlighting that, for highly permeable inhibitors, reducing extracellular inhibitor concentrations through hydrolysate detoxification, pretreatment optimization, or dilution may be necessary to mitigate futile cycling and improve microbial productivity in industrial biorefineries.

Zymomonas mobilis↗
DOE OSTI2024

On Caprock Seal Integrity of Tuscaloosa Mudstone at Cranfield, MS (USA), CO 2 Injection Site

Saline formations deep underground are among the most promising targets for large-scale geologic carbon storage. The caprock seal integrity evaluation is an important component of commercial-scale CO 2 sequestration projects. Measurements of the porosity and permeability of mudstone samples from the NETL-supported Cranfield Project were performed using a helium porosimeter and a core flow apparatus, before and after exposure in a CO 2 -saturated brine environment. The permeability of the core samples rapidly decreased with the increase in confining pressure and did not fully recover after decompression. On the other hand, exposure to CO 2 led to an increase in the permeability by at least an order of magnitude. The porosity changes after the exposure were not substantial. The post-exposure increase in permeability was subsequently offset by its rapid decrease during the higher-pressure confinement. Extrapolation of the observed permeability trends to in situ reservoir conditions suggests that Tuscaloosa mudstone can effectively serve as a natural seal.

58 GEOSCIENCES↗
DOE OSTI2025

Ultra-selective CANAL Polymers for Hydrogen-based Membrane Separations after Long-Term Aging

Hydrogen is a critically important molecule in the chemicals and energy industries, driving the need for efficient purification technologies such as membrane-based separations. However, polymer membranes often suffer from physical aging, leading to decreased gas permeability over time. This study leverages the unique aging behavior of contorted microporous polymers, synthesized via catalytic arene-norbornene annulation (CANAL) polymerization, to enhance molecular sieving for H 2 separations. While prior studies demonstrated a 1000% increase in H 2 /CH 4 selectivity with minimal permeability loss over 150 days, we extend this investigation to hyperaged (∼1 year) CANAL polymers and evaluate their performance under various conditions of industrial relevance. We report a remarkable 6100% increase in H 2 /CH 4 selectivity and a 2400% improvement in H 2 /N 2 selectivity, with only a 55% decline in H 2 permeability. These membranes exhibit excellent long-term stability in binary and ternary gas mixtures and at elevated temperatures. Additionally, extended aging enhances H 2 /CO 2 selectivity beyond the 2008 upper bound due to improved size-sieving. High-temperature permeation experiments and grand canonical Monte Carlo simulations reveal that H 2 permeability increases from 570 to 3500 barrer as the temperature rises from 35 °C to 190 °C, while H 2 /CO 2 selectivity declines from 14 to 3 due to the higher diffusion activation energy of CO 2 . Despite this tradeoff, hyperaged CANAL membranes maintain performance beyond the 200 °C upper bound, having selectivity–permeability performance comparable to many state-of-the-art membranes for H 2 /CO 2 separations. As a result, this study highlights the robustness and industrial viability of ultra-selective CANAL polymers for hydrogen purification, contributing to advancements in clean energy and sustainable separation technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Relationship between microporous structure and light gas transport through glassy polymeric membranes revealed by molecular simulations

Microporous glassy polymers are attractive materials for gas separation membranes, due to their high permeability and tailorable selectivity, provided physical aging can be delayed. The archetypal microporous glassy polymer, PTMSP, can be blended with a hyper-crosslinked isatin–triptycene porous polymer network (PPN) to delay physical aging. However, while PPN is effective at reducing physical aging, it also affects the permeability of light gases through PTMSP. Molecular dynamics simulations were used here to shed fundamental light on the mechanisms responsible for these effects. Atomistic models are developed that satisfactorily reproduce experimental observations such as matrix density and cavity size distributions for neat PTMSP as well as for PTMSP–PPN blends. Analysis of the simulation results suggests that physical aging is delayed because the PPN inclusions slow down PTMSP relaxation while reducing the connectivity between free volume pockets. To understand how PPN inclusions affect light gas permeability, the atomistic models developed are used to probe CO 2 and CH 4 diffusion and sorption. These simulations are conducted for PTMSP matrices exhibiting varying density, towards reproducing experimental permeability data. Interrogating the simulation trajectories, it is discovered that while CH 4 travels preferentially through the free volume cavities, CO 2 preferentially interacts with the available surfaces, especially in the presence of PPN. These differences help interpret experimental observations. The transport of C 2 H 6 and H 2 S through PTMSP matrices was also investigated. These gases, and in particular H 2 S, were found to absorb within PTMSP, yielding very low diffusion coefficients. The ability to predict differences in diffusion pathways for various natural gas components unveils the possibility of engineering membranes to control permeability and selectivity towards large-scale applications.

Bao Le, Tran Thi [Univ. of Oklahoma, Norman, OK (U↗
DOE OSTI2024

Impact of Pressure-Dependent Interfacial Tension and Contact Angle on Capillary Heterogeneity Trapping of CO2 in Storage Aquifers

Summary Carbon dioxide (CO2) capillary trapping increases the total amount of CO2 that can be effectively immobilized in storage aquifers. This trapping, manifesting itself as accumulated CO2 columns at a continuum scale, is because of capillary threshold effects that occur below low-permeability barriers. Considering that capillary pressure is dictated by heterogeneous pore throat size, the trapped CO2 column height and associated CO2 saturation will vary spatially within a storage aquifer. This variation will be influenced by two pressure-dependent interfacial parameters—CO2/brine interfacial tension (IFT) and CO2/brine/rock contact angle. Our objective is to understand how the pressure dependence of these two parameters affects the heterogeneity of capillary trapped CO2 at a continuum scale. Our conceptual model is a 1D two-zone system with the upper zone being a flow barrier (low permeability) and the lower zone being a flow path (high permeability). The inputs to this model include microfacies-dependent capillary pressure vs. saturation curves and permeability values. The input capillary pressure curves were collected in the literature that represents carbonate microfacies (e.g., dolograinstone) in a prevalent formation in the Permian Basin. We then used the Leverett j-function to scale the capillary pressure curve for the two zones that are assigned with the same or different microfacies. During scaling, we considered the influence of pressure on both the IFT and contact angle of CO2/brine/dolomite systems. We varied the zone permeability contrast ratio from 2 to 50. We then assumed capillary gravity equilibriums and calculated the CO2 saturation buildup corresponding to various trapped CO2 column heights. The CO2 saturation buildup is defined as the CO2 saturation in the lower layer minus that in the upper one. We found that the saturation buildup can be doubled when varying pressure in a storage aquifer, after considering pressure-dependent IFT and contact angles. Thus, assuming these two parameters to be constant across such aquifers would cause large errors in the quantification of capillary trapping of CO2. The whole study demonstrates the importance of considering pressure-dependent interfacial properties in predicting the vertical distribution of capillary trapped CO2. It has important implications in developing a better understanding of leakage risks and consequent storage safety.

Engineering↗
DOE OSTI2025

Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties

This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

University of Hawai‘i, Shallow Geothermal Resources: Energy Technology Innovation Partnership Project (Final Report)

Scientists at Lawrence Berkeley National Laboratory (Berkeley Lab) have teamed up with the University of Hawai‘i at Manoa (UH Manoa) through the U.S. Department of Energy’s Energy Technology Innovation Partnership Project to evaluate the technological and market feasibility of shallow geothermal heat exchanger (GHE) technology. UH requested this analysis to evaluate opportunities in building cooling, energy efficiency, and emissions reduction applications in Hawai‘i. UH has an abundance of geologic and geothermal data and is looking to the national labs’ expertise to execute this analysis. UH is also interested in investigating policy, regulatory, and business conditions advantageous for implementation of a pilot project and more broad deployment of this technology in Hawai‘i. In many locations around the world, the demands for heating and cooling are roughly balanced over the course of the year, so GHEs do not cause significant long-term changes in subsurface temperature. This is not the case in Hawai’i, where the demand for heating is very small, meaning that, over time, GHEs will add heat to the subsurface. If temperatures increase significantly, GHE systems will not work as designed. Regional groundwater flow has the potential to sweep heated water away from boreholes, thereby maintaining the functionality of the GHE system. Significant regional groundwater flow requires two things: a sufficiently large driving hydraulic head gradient (usually closely related to surface topography), and sufficient porosity and permeability to enable groundwater to flow in large enough quantities to enable near-borehole temperatures to be maintained at ambient values. Hawai‘i’s volcanic terrain offers ample surface topographic variation. The lava itself shows an extremely large range of porosity and permeability, so sites with large enough values of these properties must be selected. Numerical modeling of coupled groundwater and heat flow can be used to determine how large is large enough. Primarily, closed-loop systems have been investigated. Other options considered are open-loop systems and using cool seawater as the chilling source. Project work investigated the feasibility of GHE technology at two scales. At the island scale, GIS layers of various attributes relevant for GHE were combined to develop an overall favorability map for employing GHE in Hawai‘i. At the local scale, a hydrogeologic model for the subsurface component of a closed-loop system was developed for the Stan Sheriff Center at the UH Manoa campus. This site is considered promising because the rock below and immediately downgradient of the borefield is highly permeable, consisting of a subsurface karst system (limestone containing high-permeability open channels), which is underlain by a thick, high-permeability fractured basalt. Moreover, the site is near the base of the Ko‘olau Range, providing a large hydraulic head gradient. Thus, groundwater flow through the site is expected to be large, enabling efficient removal of heated groundwater. A full-GHE-system model of the site was also developed, with a simplified representation of the subsurface, in which groundwater flow is not considered and heat transfer is purely by conduction. Using the building cooling load data provided by UH, simulation results show that with groundwater flow present, a GHE can operate successfully for at least 10 years, but with no groundwater flow, the subsurface begins to heat up after only one year of operation, making the GHE unviable within 2-6 years. The team also developed a techno-economic model for this site to compare the cost of cooling using a GHE system with the costs of operating the current air-conditioning system. The GHE system is advantageous economically if favorable tax incentives and interest rates can be obtained.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

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

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

15 GEOTHERMAL ENERGY↗
NASA NTRS1976

Hydrogen transport in iron and steel

The permeabilities of protium, deuterium, and tritium in foil specimens of Marz grade iron, 4130 steel, Armco iron, HP-9-4-20, and T-1 steels were studied at hydrogen pressures between 0.02 and 0.5 MPa over the temperature range 260-700 K. The permeability was measured by a pressure-rise method, deuterium counting with a detector, and radioactivity counting. Good agreement is found between the measurement techniques used. It is shown that the permeabilities of protium, deuterium, and tritium in iron and T-1 steel at temperatures as low as 260 K are in good agreement with the equation proposed by Gonzalez (1967). However, the permeabilities of HP-9-4-20 and 4130 steel to hydrogen are typically lower than predicted. The isotope effect on hydrogen permeability of HP-9-4-20, 4130 and T-1 steels, and high-purity iron can be estimated by an inverse square root of mass correction.

Louthan, M. R., Jr.↗
NASA NTRS2002

Controls on Martian Hydrothermal Systems: Application to Valley Network and Magnetic Anomaly Formation

Models of hydrothermal groundwater circulation can quantify limits to the role of hydrothermal activity in Martian crustal processes. We present here the results of numerical simulations of convection in a porous medium due to the presence of a hot intruded magma chamber. The parameter space includes magma chamber depth, volume, aspect ratio, and host rock permeability and porosity. A primary goal of the models is the computation of surface discharge. Discharge increases approximately linearly with chamber volume, decreases weakly with depth (at low geothermal gradients), and is maximized for equant-shaped chambers. Discharge increases linearly with permeability until limited by the energy available from the intrusion. Changes in the average porosity are balanced by changes in flow velocity and therefore have little effect. Water/rock ratios of approximately 0.1, obtained by other workers from models based on the mineralogy of the Shergotty meteorite, imply minimum permeabilities of 10(exp -16) sq m2 during hydrothermal alteration. If substantial vapor volumes are required for soil alteration, the permeability must exceed 10(exp -15) sq m. The principal application of our model is to test the viability of hydrothermal circulation as the primary process responsible for the broad spatial correlation of Martian valley networks with magnetic anomalies. For host rock permeabilities as low as 10(exp -17) sq m and intrusion volumes as low as 50 cu km, the total discharge due to intrusions building that part of the southern highlands crust associated with magnetic anomalies spans a comparable range as the inferred discharge from the overlying valley networks.

Harrison, Keith P.↗
NASA NTRS2010

Microfluidic device, and related methods

A method of making a microfluidic device is provided. The method features patterning a permeable wall on a substrate, and surrounding the permeable wall with a solid, non-permeable boundary structure to establish a microfluidic channel having a cross-sectional dimension less than 5,000 microns and a cross-sectional area at least partially filled with the permeable wall so that fluid flowing through the microfluidic channel at least partially passes through the permeable wall.

Wong, Eric W.↗
NASA NTRS2016

Wind Tunnel Test of Subscale Ringsail and Disk-Gap-Band Parachutes

A subsonic wind tunnel test was conducted to determine the drag and static aerodynamic coefficients, as well as to capture the dynamic motions of a new Supersonic Ringsail parachute developed by the Low Density Supersonic Decelerator Project. To provide a comparison against current Mars parachute technology, the Mars Science Laboratory's Disk-Gap-Band parachute was also included in the test. To account for the effect of fabric permeability, two fabrics ("low" and "standard" permeability) were used to fabricate each parachute canopy type, creating four combinations of canopy type and fabric material. A wide range of test conditions were covered during the test, spanning Mach numbers from 0.09 to 0.5, and static pressures from 103 to 2116 pounds per square inch (psf) (nominal values). The fabric permeability is shown to have a first-order effect on the aerodynamic coefficients and dynamic motions of the parachutes. For example, for a given parachute type and test condition, models fabricated from "low" permeability fabric always have a larger drag coefficient than models fabricated from "standard" permeability material. This paper describes the test setup and conditions, how the results were analyzed, and presents and discusses a sample of the results. The data collected during this test is being used to create and improve parachute aerodynamic databases for use in flight dynamics simulations for missions to Mars.

Zumwalt, Carlie H.↗
NASA NTRSDate not supplied

Effects of Porous Gap Fillers on 30P30N Leading-Edge SlatNoise. Part I: Surface Pressure and Acoustics

The leading edge slat of a high-lift system is one of the main contributors to airframe noise during approach. In a previous experimental study, we assessed the performance of an impermeable slat gap filler as a passive flow control device to reduce the slat noise associated with a two-dimensional, three-element high-lift airfoil. The present paper, the first of two parts,represents a follow-on investigation to assess the relative efficacy of permeable gap fillers thatallow successively higher amounts of flow to pass through the gap. To evaluate the influence of this passive flow control device on the acoustics generated by the unsteady flow near the slat, experiments are conducted in an anechoic wind tunnel by mounting a gap filler to the slat element of the two-dimensional 30P30N high-lift configuration. Measurements are performed at a single geometric angle of attack (𝛼𝑘=8◦) and three different flow speeds that correspond to Reynolds numbers of𝑅𝑒𝑐=1.2𝑒6,1.5𝑒6, and1.71𝑒6, respectively. Steady surface pressure measurements are used to gauge the influence of the permeable gap filler treatments on the overall lift. The effect of each treatment on the radiated noise is analyzed via acoustic array measurements, followed by delay-and-sum beamforming to locate the slat noise sources and to provide the integrated acoustic spectra. The porous gap fillers are found to eliminate the narrowband peaks in the acoustic spectra and, also, to yield a 10 dB reduction in the broadband noise in comparison with the baseline case with no gap filler. The porous gap filler with thelowest permeability acts similar to the impermeable gap filler examined previously. However, the aerodynamics and noise reduction both degrade with increasing permeability. An accompanying abstract describes the particle image velocimetry measurements of the flowfield within the slat cove and on either side of the permeable gap f

slat noise↗