SEARCH · Search NASA
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.
Permeability and porosity of Precambrian Onverwacht cherts and other low permeability rocks, discussing geological origin of rock-contained organic compounds
OBJECTIVES: To investigate the physicochemical effect of silver diamine fluoride (SDF) by correlating permeability with mineral density and elemental composition of hypomineralized enamel and carious dentin. METHODS: Enamel and dentin from human carious primary teeth with and without SDF treatment in-vivo, and hypomineralized enamel from permanent molars with and without SDF treatment in-vitro were scanned using micro X-ray computed tomography. Spatial maps of biometals (calcium, zinc), phosphorus, and silver were generated using X-ray fluorescence microprobe. Permeabilities were computed using Porous Microstructure Analysis software. RESULTS: The intrinsic permeability of SDF-treated carious dentin was 14.3 % lower than untreated sound dentin (6.39e-15 ± 3.01e-15 m² vs 7.46e-15 ± 1.82e-15 m²; P < 0.0001), while untreated carious dentin was 98.4 % higher (1.48e-14 ± 7.11e-15 m²; P < 0.0001). SDF-treated and untreated transparent dentin showed similar reduced permeabilities (75.6 % and 78.4 % lower than untreated sound dentin, respectively; P = 0.93). Severely hypomineralized enamel showed permeability reaching 108.1 % of adjacent sound dentin (5.71e-15 ± 2.04e-15 m² vs 5.28e-15 ± 1.30e-15 m²; P = 0.1409) and was significantly higher than mildly hypomineralized enamel (1.39e-15 ± 1.04e-15 m²; P < 0.0001). SDF treatment did not significantly impact the permeability of severely hypomineralized enamel (12.4 % reduction; P = 0.07). Principal component regression identified Zn level as a significant effector of tissue permeabilities in carious primary teeth (P < 0.0001). SIGNIFICANCE: This study introduces a computational method to measure dental tissue permeability, and demonstrates that SDF significantly reduces permeability in carious dentin but not intact hypomineralized enamel. The study reveals biometal Zn localization can alter dentin and enamel permeabilities, providing new insights into pathobiological mechanisms underlying caries and hypomineralization.
Directional permeability variations, which govern directional fluid flow in porous media with anisotropy, are important to accurately predict flow behavior, reactive transport, and fluid–solid interactions for various processes such as enhanced geothermal systems, energy storage devices, and biological systems. However, the intricate architecture of porous media makes it difficult to predict directional permeabilities. In this work, we present a novel machine learning (ML) framework, AnisONet, built upon an integration of a convolutional neural network, Swin transformer, and the deep operator network architecture, designed to predict anisotropic permeability and upscale predictions to larger spatial domains. First, AnisONet was evaluated with three classes of two-dimensional (2D) porous media, including synthetic circular and elliptical grains and natural sandstone grains from micro-computed tomography images. A lattice Boltzmann model (LBM) was used to calculate directional permeabilities at every 10° angle, producing 19 data points per image of porous media. AnisONet is then trained to predict permeability as a function of rotation angle. AnisONet showed strong predictive capability of directional permeability. Second, we tested our model for five upscaling cases with a large image size in the finite-element method (FEM) for 2D Darcy flow with various permeability tensor construction methods. Overall, upscaled permeability tensors in FEM simulations produce a reasonably good match with LBM results, highlighting the importance of selecting appropriate tensor formation strategies for accurate permeability upscaling. AnisONet, as a directional permeability estimator, could be further developed for more complex geometries, with the potential to develop a foundational ML model for various applications in porous media.
Gas hydrate, a solid combination of gas (mostly methane in nature) and water molecules stable at low temperatures and elevated pressures, occurs naturally in marine and permafrost-associated environments. Gas hydrate reservoirs, such as those in the Alaska North Slope, have been considered potential energy resources for gas production. To understand the petrophysical and geo-mechanical characteristics of the reservoir, core samples retrieved from the site of the JOGMEC-DOE-USGS collaborative gas hydrate R&D project have been analyzed in the laboratory for their hydraulic and mechanical properties. This paper focuses on both seal and reservoir samples associated with the B1 and D1 sands, which are evaluated for index properties (including porosity, grain size distribution, liquid and plastic limits, specific surface area, and specific gravity), consolidation, permeability, and water retention. Furthermore, the reservoir core samples were tested with pore-filling, laboratory-grown tetrahydrofuran hydrate, in order to assess reservoir behavior during gas production from hydrates. Under simulated in situ stress conditions, the seal and hydrate-free reservoir cores had a permeability anisotropy ratio of k h /k v = 3.0−5.0, and k h /k v = 2.4−3.0 for the reservoir tetrahydrofuran hydrate-bearing cores. The data suggest that depressurizing the reservoir to induce hydrate dissociation alters the reservoir effective permeability in three ways: permeabilities decrease due to porosity lost (e.g., the initial reservoir thickness can decrease by up to 5% upon 7 MPa depressurization), permeability increases due to the loss of solid hydrate in the pore space, and permeability anisotropy k h /k v decreases in response to the evolving pore-space geometry. We show that given the simulated in situ gas hydrate saturations (i.e., S h = 32% in core 7P-2E and S h = 21% in core 20P-4), gas production from the dissociation of tetrahydrofuran hydrate in the two tested cores results in a net increase in effective permeability and a decrease in k h /k v . This study highlights the importance of investigating seal and reservoir sediments and the impacts of depressurization on the porosity and permeability responses during production.
Carboxysomes are protein microcompartments found in cyanobacteria, whose shell encapsulates rubisco at the heart of carbon fixation in the Calvin cycle. Carboxysomes are thought to locally concentrate CO 2 in the shell interior to improve rubisco efficiency through selective metabolite permeability, creating a concentrated catalytic center. However, permeability coefficients have not previously been determined for these gases, or for Calvin-cycle intermediates such as bicarbonate (HCO$^–_3$), 3-phosphoglycerate, or ribulose-1,5-bisphosphate. Starting from a high-resolution cryogenic electron microscopy structure of a synthetic β-carboxysome shell, we perform unbiased all-atom molecular dynamics to track metabolite permeability across the shell. The synthetic carboxysome shell structure, lacking the bacterial microcompartment trimer proteins and encapsulation peptides, is found to have similar permeability coefficients for multiple metabolites, and is not selectively permeable to HCO$^–_3$ relative to CO 2 . To resolve how these comparable permeabilities can be reconciled with the clear role of the carboxysome in the CO 2 -concentrating mechanism in cyanobacteria, complementary atomic-resolution Brownian Dynamics simulations estimate the mean first passage time for CO 2 assimilation in a crowded model carboxysome. Despite a relatively high CO 2 permeability of approximately 10 -2 cm/s across the carboxysome shell, the shell proteins reflect enough CO 2 back toward rubisco that 2,650 CO 2 molecules can be fixed by rubisco for every 1 CO 2 molecule that escapes under typical conditions. The permeabilities determined from all-atom molecular simulation are key inputs into flux modeling, and the insight gained into carbon fixation can facilitate the engineering of carboxysomes and other bacterial microcompartments for multiple applications.
An experimental investigation was made to determine the permeability and strength characteristics of a number of sintered, porous, hollow turbine rotor blades and to determine the effectiveness of the blade fabrication method on permeability control. The test blades were fabricated by the American Electro Metal Corporation under a contract with the Office of Naval Research, Department of the Navy, and were submitted to the NACA for testing. Of the 22 test blades submitted, ten were sintered but not coined, five were sintered and coined, and seven were sintered and not coined but contained perforated reinforcements integral with the blade shells. Representative samples of each group of blades were tested. Large variations in permeability in both chordwise and spanwise directions were found. Local deviations as large as 155 to -85 percent from prescribed values were found in chordwise permeability. Only one blade, an uncoined one, had a chordwise permeability variations which reasonably approached that specified. Even for this blade, local deviations exceeded 10 percent. Spanwise permeability, specified to be held constant, varied as much as 50 percent from root to tip for both an uncoined and a coined blade. Previous NACA analyses have shown that in order to maintain proper control of blade wall temperatures, permeability variations must not exceed plus or minus 10 percent. Satisfactory control of permeability in either the chordwise or the spanwise direction was not achieved in the blades tested. Spin tests made at room temperature for six blades revealed the highest material rupture strength to be 8926 pounds per square inch. This value is about one third the strength required for rotor blades in present-day turbojet engines. The lowest value of blade strength was 1436 pounds per square inch.
Artificial water channels (AWCs) have emerged as a promising framework for stable water permeation, with water transport rates comparable to aquaporins (3.4–40.3 × 10 8 H 2 O/channel/s). In this study, we probe the influence of ring-size and side-chain length on the water permeability observed within a class of AWCs termed ligand-appended pillar[n]arenes (LAPs) that have an adjustable ring-size (m) and side-chain length (n). Through all-atom molecular dynamics simulations, we calculate the permeability of these channels using the collective diffusion model and find their permeabilities. We characterize the mechanistic influence of pillar[n]arene ring-size and side-chain length on the channel water permeability by analyzing the characteristics of the internal permeating water-wire and the surrounding channel structure. We observe that water permeability decreases as a function of increasing ring-size due to increases in hydrophilic contacts between the permeating water-wire and the oxygen groups on the channel wall. Further, we observe an increase in water permeability as a function of side-chain length due to increased partitioning of the channel terminal groups into the hydrophilic blocks of the surrounding bilayer. For the LAP6 channel, with increase in side-chain length, the distance between terminal groups increases and leads to an increase in pore size, thereby enhancing water permeability. In the case of LAP5, as side-chain length increases, the channel displays a compensatory effect between tilt and bend angle due to the flexible side-chains. Such flexibility leads to higher terminal group partitioning in the hydrophilic blocks of the bilayer and extends the permeating water-wire. Furthermore, this increase in water-wire length and hydrophilic block access overcomes the nonmonotonic pore size trend in pillar[5]arene channels.
Abstract Fractures and faults represent planes of weakness and compliance in rock masses that serve as focal points for both microearthquakes and fluid transport, with seismicity and permeability evolution closely linked. Contact stiffness is highly stress‐sensitive and directly influences permeability. We explore the co‐evolution of specific stiffness and permeability of rough fractures under normal stress and shear offset using numerical simulations. Individual rough fractures are represented by variable amplitude (Root mean square) and wavelength ( λ ) using a granular mechanics model. Contacting rough surfaces are mated, offset in shear, and then compacted in displacement mode. The compacting fractures generate stress‐dependent changes in contact porosity, which govern both permeability and stiffness evolution. We establish a universal dimensionless relationship linking specific stiffness and permeability that inherently incorporates the effects of surface roughness, shear offset, and microcracking. The observed cracking effect—where local stress redistribution and pressure‐driven microcrack propagation dynamically alter the aperture field—introduces a nonlinear permeability response at high stress. Increased roughness amplitude and larger shear offsets reduce stiffness while dampening permeability sensitivity to stress, demonstrating a strong interplay between surface texture and hydro‐mechanical behavior. While the model captures this behavior effectively, deviations emerge at very low porosities due to extreme aperture sensitivity in this limit.
Optical fiber in a borehole can be interrogated with distributed acoustic sensors (DAS) to capture fracture displacements with the potential to map surrounding fracture networks. We designed a laboratory experiment to test the capability of DAS to determine borehole flow characteristics, and we show that for the first time DAS can be used to remotely estimate permeability. Optical fiber was wrapped around a bead filled pipe and the pressure drop and flow velocity were measured to directly calculate permeability. A machine learning model using statistical features from continuous DAS estimated the bulk permeability. Fluid interactions with the permeable material demonstrate insufficient resolution using DAS amplitude-based measurements for estimating pressure drop to infer permeability. Variations in the spectral domain relate DAS measurements to the pressure drop and provide consistent permeability estimates. Resolution with DAS is sufficient to estimate permeability and provides a reliable method to monitor at depth in borehole conditions.
Supported ionic liquid membranes (SILMs), containing phosphonium ionic liquids with aprotic N-heterocyclic anions (AHA ILs) in an inert inorganic support, were tested under both dry and humidified (40 % RH) mixed-gas conditions down to 420 ppm CO 2 in N 2 at 35 °C. In the dry case, the best performing IL, triethyl(octyl)phosphonium 4-bromopyrazolide ([P 2228 ][4-BrPyra]) exhibited mixed-gas CO 2 permeabilities and CO 2 /N 2 permeability selectivities as high as 26,800 barrer and 7,000, respectively. In the presence of humidity, the CO 2 permeability and CO 2 /N 2 selectivity increased to 49,100 barrer and 13,200, respectively, and these are the highest reported combination in the literature. Humidity amplifies CO 2 permeabilities and CO 2 /N 2 permeability selectivities through increases in CO 2 capacity due to bicarbonate formation and through faster mobility of the mobile carrier from decreased viscosity. Furthermore, N 2 permeability stayed roughly invariant in the presence of humidity, likely from competing effects of viscosity reduction and lower N 2 solubility.
Caprock formations are a crucial part of subsurface-engineered systems. Composed largely of shale, caprocks act as natural barriers that prevent the upward migration of fluids, thereby ensuring the containment of stored substances in subsurface formations. Fractures in these formations are potential leakage pathways for stored fluids. Mineral precipitation reactions in these fractures, particularly calcite, can significantly restrict the fluid permeability, reducing leakage potential. However, predictive capabilities of mineral precipitation in fractures and associated permeability evolution are limited due to a lack of fundamental understanding of such reactions in natural samples, complicated by mineral heterogeneity and the complexity of the fracture structure. In this study, 3D-printed fracture samples are used to understand the impact of carbonate mineralization on fracture aperture and permeability evolution. Samples were printed using a digital light processing (DLP) 3D printer and commercial liquid resin. Calcite precipitation was first tested on printed 2D films before conducting plug flow column experiments aimed to understand fracture permeability changes due to mineral precipitation. Contact angle measurement and Fourier transform infrared (FTIR) spectroscopy on printed 2D films show evidence of a substantial amount of surface energy for calcite nucleation and precipitation. Surface topography analysis of printed fractured surfaces reveals comparable values, highlighting the high replicability of the printed samples. During the column experiments, the permeability reduces exponentially due to a decrease in fracture aperture. Reductions in fracture aperture estimated from effluent concentration and 3D X-ray computed tomography (CT) show comparable results. Moreover, 3D X-ray CT images suggest the impact of local flow velocities on precipitation. The insights gained from this research contribute to a deeper understanding of the permeability evolution due to carbonate mineralization in caprock formations.
A recent chemical explosive test in P-Tunnel at the Nevada National Security Site, Nevada, USA, was conducted to better understand how signals propagate from explosions in the subsurface. A primary signal of interest is the migration of gases that can be used to differentiate chemical from nuclear explosions. Gas migration is highly dependent on the rock permeability which is notoriously difficult to determine experimentally in the field due to a potentially large dependence on the scale over which measurements are made. Here, we present pre-explosion permeability estimates to characterize the geologic units surrounding the recent test. Permeability measurements were made at three scales of increasing size: core samples (≈2 cm), borehole packer system tests (≈1 m), and a pre-shot cavity pressurization test (> 10 m) across ten tuff units. Permeability estimates based on core measurements showed little difference from borehole packer tests. However, permeability in most rock units calibrated from cavity pressurization tests resulted in higher permeability estimates by up to two orders of magnitude. Here, we demonstrate that the scale of the measurement significantly impacts the characterization efforts of hydraulic properties in volcanic tuff, and that local-scale measurements (< 10 m scale) do not incorporate enough heterogeneity to accurately predict field-scale flow and mass transport.
The widely used Darcy’s law specifies a linear relation between the Darcy velocity of fluid flow and the pressure gradient that drives the flow. However, studies have shown that Darcy velocity can exhibit a nonlinear dependence on the pressure gradient in low-permeability porous media such as clay and shale when the pressure gradient is adequately low. This phenomenon is referred to as low-velocity non-Darcian flow or pre-Darcy flow. This paper provides a comprehensive review of the theories, experimental data, and modeling methods for pre-Darcy flow in low-permeability porous media. The review begins by outlining the fundamental mechanisms underlying pre-Darcy flow that regulate the unique characteristics such as nonlinear dependence of the Darcy velocity on the pressure gradient and its relevance to fluid–rock interactions. The review then proceeds to present a thorough compilation of experimental investigations performed in various low-permeability geomaterials including tight sandstones, shales, and clays. Next, empirical and theoretical models and simulation methods that have been developed to fit and interpret experimental data are reviewed. Finally, the review underscores the challenges encountered in conducting and interpreting pre-Darcy flow experiments and suggests future research directions. By analyzing previous experimental investigations, this review aims to offer a valuable resource for researchers and practitioners seeking to enhance their understanding of fluid dynamics in low-permeability geomaterials. This provides insights into the application of pre-Darcy flow in numerous natural and engineered processes such as shale oil and gas recovery, contaminant transport in low-permeability aquifers, and geological disposal of nuclear waste.
Accurate permeability estimation is essential across subsurface engineering applications but remains challenging due to the complex pore structures of natural geomaterials. Traditional empirical methods and simplified theoretical models often inadequately capture the role of pore size distribution and connectivity. Here, this study develops a statistical and simulation-informed permeability model that collapses pore-scale complexity into a compact scaling of the form k = αϕμ d 2 , where ϕ is porosity, μ d is mean pore size, and α is a weakly varying coefficient. By combining pore network simulations with statistical analysis of unimodal and bimodal pore size distributions, we identify three key findings: (i) permeability is much more sensitive to mean pore size than to porosity; (ii) across extensive datasets, the ratio σ d /μ d (standard deviation to mean) clusters around a characteristic value ∼0.4, allowing the effects of the full pore size distribution to be represented by μ d and a narrowly varying α ≈ 0.05; and (iii) for bimodal systems, there exists a critical fraction of small pores ∼0.78 above which flow becomes small-pore dominated, enabling the definition of an effective flow-controlling pore population and facilitating simplified permeability estimation for such systems. The resulting model, which requires only porosity and a representative mean pore size as inputs, is validated against comprehensive experimental datasets (>1700 samples) spanning diverse soils and rocks and achieves good predictive accuracy. Overall, this work provides a physically grounded yet practically simple permeability estimator suitable for subsurface engineering, environmental protection, and resource management applications.
Relative permeability is a crucial two-phase property in porous media that can be significantly impacted by wettability conditions. While traditional research has predominantly examined homogeneous wettability, this work explores the less studied pore-size dependent (PSD) wettability, featured by a pore-size dependent wettability distribution. Leveraging high-fidelity Lattice Boltzmann simulations on CT-scanned porous samples, we demonstrate how PSD wettability would impact relative permeability at the pore scale. Our findings reveal that the deviation of relative permeability from the homogeneous wettabilitity induced by PSD wettability can be 5%–20%. The deviation of relative permeability curves increases as the spanning range of the contact angle increases. We also find that this impact is less pronounced as the capillary number increases. By adopting a pore-size-dependent contact angle relationship, our approach provides a more accurate and nuanced understanding of how PSD wettability would impact two-phase flow. These findings discovered at the pore scale may also provide valuable insights on relative permeability at the core to reservoir scales.
The effect of oxide films developed in situ from CO/CO2 doped hydrogen on high pressure hydrogen permeability at 820 C was studied on N-155, A-286, IN 800, 19-9DL, Nitronic 40, HS-188, and IN 718 tubing in a Stirling materials simulator. The hydrogen permeability decreased with increasing dopant levels of CO or CO2 and corresponding decreases in oxide porosity. Minor reactive alloying elements strongly influenced permeability. At high levels of CO or CO2, a liquid oxide formed on alloys with greater than 50 percent Fe. This caused increased permeability. The oxides formed on the inside tube walls were analyzed and their effective permeabilities were calculated.
In this study, we link the multiphase flow measurements with controls of sedimentary structures (e.g., heterogeneity and anisotropy) on relative permeability to variations in ultrasonic velocities for two deep saline aquifer formations (Entrada and Bluff Sandstones) in the San Juan Basin of the Southwestern USA. The rock specimens were extracted from outcrop sites near Durango, CO, USA. They have distinct differences in grain size, cementation composition, and individual chemical amounts, despite both formations being eolian sandstones. We performed a series of unsteady- and steady-state CO 2 -brine relative permeability experiments under capillary-controlled displacement rates. Unsteady-state experiments were conducted at 71 °C and 9.65 MPa; steady-state experiments were conducted at 85 °C and 22.8 MPa and 89 °C and 24.1 MPa for the Bluff and Entrada Sandstones, respectively. During the unsteady experiments, X-ray computed tomography was used to visualize multiphase flow in porous media and quantify saturations during brine drainage under various flow rates. Scan images and saturation profiles indicate that the CO 2 distribution in the pore volume was strongly impacted by the presence of high-angle cross-laminations, heterogeneous rock structure, and direction of bedding orientation. Those factors contribute to dramatic and quick initial breakthroughs and affect the overall saturation dynamics. Steady-state relative permeability tests were conducted at net flow rates of 1 mL/min for both brine drainage and imbibition. During the experimental steps, the CO 2 fractional flow was increased and decreased for both drainage and imbibition scenarios to mimic the front when CO 2 contacts brine and behind the front when brine enters space previously occupied by CO 2 . It was found that compressional velocity decreased, while shear waves slightly increased as brine saturation decreased. The hysteresis effects for the relative permeability and acoustic velocities were distinct. Furthermore, the CO 2 /brine front stability is quantified by applying a mobility ratio approach to spot saturations at which the boundary line between fluids becomes uneven. The results presented in this work can potentially boost the quality and precision of forecasts for the CO 2 storage projects in which the vertical and horizontal core-scale heterogeneity and anisotropy impact the plume migration within host reservoirs.
Developing a fundamental understanding of the effects of varying ligand chemistries on mass transport rates is key to designing membranes with solute-specific selectivity. While permeation cells offer a robust method to characterize membrane performance, they are limited to assessing a single membrane chemistry or salt solution per test. As a result, investigating the effects of varying ligand chemistries on membrane performance can be a tedious process, involving both the preparation of multiple samples and numerous, time-consuming permeation tests. This study uses digital light processing (DLP) 3D printing to fabricate a millifluidic flow-based permeation device made from a hydrogel active ester network that can be easily functionalized with ion-selective ligands. Without the need for bonding or assembly steps, ligands can be introduced and tested in the permeation device by simply injecting a small volume of a ligand solution. Various salt concentrations and molecular species can be cycled through a single device by switching the solution feeding into the salt reservoir, thereby reducing the number of samples needed for permeability and selectivity screenings. This research sets the groundwork for formulation development and postprocessing methods to 3D-print functional millifluidic devices capable of assessing solute selectivity in membranes and polymer adsorbents for aqueous separations. In this work, comparable salt permeability trends were observed with both 3D-printed devices and traditional assays. Devices were functionalized with an imidazole ligand to investigate salt permeability and selectivity of monovalent and divalent salts. Measurements showed increasing permeability for monovalent salts (NaCl) relative to divalent salts (MgCl 2 , CuCl 2 ) in functionalized membranes, with higher monovalent/divalent selectivity at increasing imidazole grafting densities. Here, the methods and findings described here represent a step toward developing higher-throughput methods with 3D-printed devices for screening the effects of ligand chemistry on mass transport rates in membrane materials.