Relative Permeability of scCO2 and Brine in Core from Potential CO2 Reservoirs
Description of efforts at NETL RIC measuring relative permeability and core properties associated with ongoing and potential CO2 storage locations.
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.
Description of efforts at NETL RIC measuring relative permeability and core properties associated with ongoing and potential CO2 storage locations.
Thin polymer films are being considered, as candidate materials to augment the permeation resistance of cryogenic hydrogen fuel tanks such as would be required for future reusable launch vehicles. To evaluate performance of candidate films after environmental exposure, an experimental study was performed to measure the thermal/mechanical and permeation performance of six, commercial-grade materials. Dynamic storage modulus, as measured by Dynamic Mechanical Analysis, was found over a range of temperatures. Permeability, as measured by helium gas diffusion, was found at room temperature. Test data was correlated with respect to film type and pre-test exposure to moisture, elevated temperature, and cryogenic temperature. Results indicated that the six films were comparable in performance and their resistance to environmental degradation.
Thermal Protection Materials Branch (TPMB) has been involved in various research programs to improve the properties and structural integrity of the existing aerospace high temperature materials. Specimens from various research programs were brought into the analytical laboratory for the purpose of obtaining and refining the material characterization. The analytical laboratory in TPMB has many different instruments which were utilized to determine the physical and chemical characteristics of materials. Some of the instruments that were utilized by the SJSU students are: Scanning Electron Microscopy (SEM), Energy Dispersive X-ray analysis (EDX), X-ray Diffraction Spectrometer (XRD), Fourier Transform-Infrared Spectroscopy (FTIR), Ultra Violet Spectroscopy/Visible Spectroscopy (UV/VIS), Particle Size Analyzer (PSA), and Inductively Coupled Plasma Atomic Emission Spectrometer (ICP-AES). The above mentioned analytical instruments were utilized in the material characterization process of the specimens from research programs such as: aerogel ceramics (I) and (II), X-33 Blankets, ARC-Jet specimens, QUICFIX specimens and gas permeability of lightweight ceramic ablators. In addition to analytical instruments in the analytical laboratory at TPMB, there are several on-going experiments. One particular experiment allows the measurement of permeability of ceramic ablators. From these measurements, physical characteristics of the ceramic ablators can be derived.
Abstract We present a synthetic multi‐scale, multi‐physics dataset constructed from the Kimberlina 1.2 CO 2 reservoir model based on a potential CO 2 storage site in the Southern San Joaquin Basin of California. Among 300 models, one selected reservoir‐simulation scenario produces hydrologic‐state models at the onset and after 20 years of CO 2 injection. Subsequently, these models were transformed into geophysical properties, including P‐ and S‐wave seismic velocities, saturated density where the saturating fluid can be a combination of brine and supercritical CO 2 , and electrical resistivity using established empirical petrophysical relationships. From these 3D distributions of geophysical properties, we have generated synthetic time‐lapse seismic, gravity and electromagnetic responses with acquisition geometries that mimic realistic monitoring surveys and are achievable in actual field situations. We have also created a series of synthetic well logs of CO 2 saturation, acoustic velocity, density and induction resistivity in the injection well and three monitoring wells. These were constructed by combining the low‐frequency trend of the geophysical models with the high‐frequency variations of actual well logs collected at the potential storage site. In addition, to better calibrate our datasets, measurements of permeability and pore connectivity have been made on cores of Vedder Sandstone, which forms the primary reservoir unit. These measurements provide the range of scales in the otherwise synthetic dataset to be as close to a real‐world situation as possible. This dataset consisting of the reservoir models, geophysical models, simulated time‐lapse geophysical responses and well logs forms a multi‐scale, multi‐physics testbed for designing and testing geophysical CO 2 monitoring systems as well as for imaging and characterization algorithms. The suite of numerical models and data have been made publicly available for downloading on the National Energy Technology Laboratory's (NETL) Energy Data Exchange (EDX) website.
Hydrogen-powered mining trucks are rapidly gaining traction as a key driver of decarbonization in the mining industry. This surge necessitates a decisive shift towards renewable energy sources, primarily for hydrogen production. However, the inherent intermittency of renewables poses a challenge for continuous operations, demanding efficient storage solutions for hydrogen fuel. Hydrogen, derived from renewables, requires large-scale storage, with underground hydrogen storage (UHS) in depleted oil and gas reservoirs projected to provide this capacity. An impermeable seal (caprock) overlying depleted reservoirs is crucial for preventing leakage and ensuring secure storage. Shale, a potential caprock, has been recognized for its effectiveness in containing natural gas for geological time-periods. This study investigates the impact of hydrogen exposure on the matrix permeability of Marcellus shale, a key factor influencing its sealing capacity. A novel time-resolved approach was employed, measuring matrix permeability after each exposure to understand the dynamic interactions between hydrogen and the shale matrix. Samples were exposed to hydrogen in pressure reactions at realistic reservoir conditions at varying time durations. Preliminary results reveal an initial increase in permeability, while subsequent exposures show a decrease, indicating potential microstructural changes.
The membrane-based electrolysis of water is a growing topic of interest due to the advantages of employing membranes in hydrogen production efficiency and system safety over the traditional alkaline water electrolysis. Ion-exchange membranes with low gas permeability are highly desirable for stable and safe operation of membrane-based water-splitting technologies, hence gas permeability through ion-exchange membranes needs to be properly assessed with standardized methods. We addressed three methods to measure gas permeability of ion-exchange membranes, a pressure permeation cell, chronoamperometry microelectrodes, and in situ testing of the membrane electrode assembly, and provide a guideline for choosing the appropriate method for the targeted operating conditions of the water electrolyzers.
Rigid fibrous refractory insulations (TPS tiles) are integral components of many spacecraft thermal protection systems. These materials are composed of refractory fibers With diameters on the order of 1 to 15 micrometers. They are lightweight and have an open, highly porous microstructure. Typical densities are less than 500 kilograms per cubic meters, and porosities generally exceed 0.8. Because of their open porosity, these materials are permeable to gas glow. There are numerous instances in which internal gas transport in a thermal protection system could be important; examples include the penetration of hot boundary-layer gases into the insulation, the flow of decomposition (pyrolysis) products from the interior, the use of convective flows to mitigate ice formation caused by cryopumping, and the design of refractory vents for pressure equilibration during atmospheric entry. Computational analysis of gas flow through porous media requires values of permeability which have not previously been available for the rigid fibrous insulations used in thermal protection systems. This paper will document measurements of permeability for a variety of insulations from NASA's LI, FRCI, and AETB families of lightweight ceramic ablators. The directional anisotropy of permeability and its dependence on gas pressure and material density will be presented. It will be shown that rarified-flow effects are significant in the flow through such materials. Connections will be drawn between the insulation microstructure and permeability. The paper will also include representative computations of flow through rigid fibrous insulations.
Two parachute fabrics, described by Parachute Industry Specifications PIA-C-7020D Type I and PIA-C-44378D Type I, were tested to obtain their permeabilities in air (i.e., flow-through volume of air per area per time) over the range of differential pressures from 0.146 psf (7 Pa) to 25 psf (1197 Pa). Both fabrics met their specification permeabilities at the standard differential pressure of 0.5 inch of water (2.60 psf, 124 Pa). The permeability results were transformed into an effective porosity for use in calculations related to parachutes. Models were created that related the effective porosity to the unit Reynolds number for each of the fabrics. As an application example, these models were used to calculate the total porosities for two geometrically-equivalent subscale Disk-Gap-Band (DGB) parachutes fabricated from each of the two fabrics, and tested at the same operating conditions in a wind tunnel. Using the calculated total porosities and the results of the wind tunnel tests, the drag coefficient of a geometrically-equivalent full-scale DGB operating on Mars was estimated.
Two parachute fabrics, described by Parachute Industry Specifications PIA-C-7020D Type I and PIA-C-44378D Type I, were tested to obtain their permeabilities in air (i.e., flow-through volume of air per area per time) over the range of differential pressures from 0.146 psf (7 Pa) to 25 psf (1197 Pa). Both fabrics met their specification permeabilities at the standard differential pressure of 0.5 inch of water (2.60 psf, 124 Pa). The permeability results were transformed into an effective porosity for use in calculations related to parachutes. Models were created that related the effective porosity to the unit Reynolds number for each of the fabrics. As an application example, these models were used to calculate the total porosities for two geometrically-equivalent subscale Disk-Gap-Band (DGB) parachutes fabricated from each of the two fabrics, and tested at the same operating conditions in a wind tunnel. Using the calculated total porosities and the results of the wind tunnel tests, the drag coefficient of a geometrically-equivalent full-scale DGB operating on Mars was estimated.
Mass spectrometry was used to measure oxygen and nitrogen permeabilities while polarography was used to measure oxygen permeabilities for several contact lens materials. Applicable sample holders were designed and fabricated to accommodate curved and flat specimens. A prepared standard was used to calibrate the mass spectrometric analyses. The oxygen permeability values determined by mass spectrometry were significantly greater than those determined by polarography. This was attributed to the phase boundary phenomena and the limiting oxygen permeance of water inherent in the polarographic technique. Polarographic values determined were in good agreement with proprietary values obtained by polarography, with the exception of one material.
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.
Abstract The Ghareb Formation is a shallowly buried porous chalk in southern Israel that is being considered as a host rock for a geologic nuclear waste repository. Setup and operation of a repository will induce significant mechanical, hydrological and chemical perturbations in the Ghareb. Developing a secure repository requires careful characterization of the rock behavior to different loads. To characterize hydromechanical behavior of the Ghareb, several short- and long-term deformation experiments were conducted. Hydrostatic loading tests were conducted both dry and water-saturated, using different setups to measure elastic properties, time-dependent behavior, and permeability. A set of triaxial tests were conducted to measure the elastic properties and rock strength under differential loading at dry and water-saturated conditions. The hydrostatic tests showed the Ghareb began to deform inelastically around 12–15 MPa, a relatively low effective pressure. Long-term permeability measurements demonstrated that permeability declined with increasing effective pressure and was permanently reduced by ~ 1 order of magnitude after unloading pressure. Triaxial tests showed that water saturation significantly degrades the rock properties of the Ghareb, indicating water-weakening is a significant risk during repository operation. Time-dependent deformation is observed during hold periods of both the hydrostatic and triaxial tests, with deformation being primarily visco-plastic. The rate of deformation and permeability loss is strongly controlled by the effective pressure as well. Additionally, during holds of both hydrostatic and triaxial tests, it is observed that when water-saturated, radial strain surpassed axial strain when above effective pressures of 13–20 MPa. Thus, deformation anisotropy may occur in situ during operations even if the stress conditions are hydrostatic when above this pressure range.
Two mechanisms have been proposed to account for solute permeation of lipid bilayers. Partitioning into the hydrophobic phase of the bilayer, followed by diffusion, is accepted by many for the permeation of water and other small neutral solutes, but transient pores have also been proposed to account for both water and ionic solute permeation. These two mechanisms make distinctively different predictions about the permeability coefficient as a function of bilayer thickness. Whereas the solubility-diffusion mechanism predicts only a modest variation related to bilayer thickness, the pore model predicts an exponential relationship. To test these models, we measured the permeability of phospholipid bilayers to protons, potassium ions, water, urea, and glycerol. Bilayers were prepared as liposomes, and thickness was varied systematically by using unsaturated lipids with chain lengths ranging from 14 to 24 carbon atoms. The permeability coefficient of water and neutral polar solutes displayed a modest dependence on bilayer thickness, with an approximately linear fivefold decrease as the carbon number varied from 14 to 24 atoms. In contrast, the permeability to protons and potassium ions decreased sharply by two orders of magnitude between 14 and 18 carbon atoms, and leveled off, when the chain length was further extended to 24 carbon atoms. The results for water and the neutral permeating solutes are best explained by the solubility-diffusion mechanism. The results for protons and potassium ions in shorter-chain lipids are consistent with the transient pore model, but better fit the theoretical line predicted by the solubility-diffusion model at longer chain lengths.
The large-scale deployment of geologic carbon storage will require minimizing the risk of CO 2 leakage. Here we report on a novel strategy to seal leakage pathways by injecting functionalized reactive calcium silicate (CaSiO 3 ) particles. The particles consist of a mineral silicate core within a temperature-sensitive Poly(N,N-dimethylaminoethyl methacrylate) coating that limits contact between the core and the surrounding CO 2 (aq) above a switching temperature making them relatively unreactive deep in the subsurface. If the CO 2 leaks to shallower depths and cooler temperatures, the coatings would extend and enable the silicate core to react with CO 2(aq) to form solid carbonate precipitates that reduce the permeability of the fracture or flow path. A combined laboratory- and simulation-based investigation was conducted to explore key early-stage research questions related to the feasibility of this leakage management strategy. Sand columns were injected with either coated or uncoated CaSiO 3 particles and the permeability was measured before and after exposure to CO 2 . Experiments were conducted at 35 °C and 70 °C, a range that spans the switching temperature for the coatings. Results show that at the cooler temperature the coating effectively enabled reaction and permeability reduction. At higher temperatures, the coated particles limited mineral carbonation and the permeability decrease was minimal compared to columns with uncoated particles. Here, the morphology and location of the precipitates provides insight about permeability evolution. Model simulations of 2D radial particle transport from an injection point provided insight into the effects of particle size, coating properties, and concentrations of particles in several representative rock formations.
Researchers at the NASA Glenn Research Center are evaluating high-temperature materials for Stirling heater heads for second- and third-generation Stirling radioisotope power systems that would help to increase the system efficiency to 30 to 35 percent and the system specific power to 8 to 10+ W/kg. Ceramic materials could make it possible for the convertor hot-end temperature to be increased to 1050 to 1200 C, in comparison to the current 650 C with an Inconel 718 heater head. A hermetically sealed Stirling heater head must retain a constant internal pressure of nearly 400-psi helium (He) throughout its useful life (120,000 hr) at the design operating temperature. Therefore, He permeability was measured for eight potential materials and compared with the permeability of the current heater head material, Inconel 718. The eight materials included silicon nitride (Si3N4), silicon dioxide (SiO2), both sintered and chemical vapor deposited (CVD) silicon carbide (SiC), alumina (Al2O3), two types of melt-infiltrated (MI) SiC/SiC composites, and a carbon/SiC composite (C/SiC). Glenn submitted samples of each material to Porous Materials, Inc., Ithaca, New York, for permeability analysis. At room temperature and 30-psi He, four materials--Si3N4, Al2O3, SiO2, and sintered SiC--demonstrated lower permeability than Inconel 718. The CVD SiC and all the composite materials were significantly more permeable to He than the baseline material.
We exploit nonlinear elastodynamic properties of fractured rock to probe the micro-scale mechanics of fractures and understand the relation between fluid transport and fracture aperture under dynamic stressing. Experiments were conducted on rough, tensile-fractured Westerly granite subject to triaxial stresses. We measure fracture permeability for steady-state fluid flow with deionized water. Pore pressure oscillations are applied at amplitudes ranging from 0.2 to 1 MPa at 1 Hz frequency. During dynamic stressing we transmit ultrasonic signals through the fracture using an array of piezoelectric transducers (PZTs) to monitor evolution of interface properties. We examine the influence of fracture aperture and contact area by conducting measurements at effective normal stresses of 10–20 MPa. Additionally, the evolution of contact area with stress is characterized using pressure sensitive film. These experiments are conducted separately with the same fracture and map contact area at stresses from 9 to 21 MPa. The measurements are a proxy for “true” contact area for the fracture surface and we relate them to elastic properties using the calculated PZT sensor footprints via numerical modeling of Fresnel zones. We compare the elastodynamic response of the fracture using the stress-induced changes in ultrasonic wave velocities for transmitter-receiver pairs to image spatial variations in contact properties. We show that nonlinear elasticity and permeability enhancement decrease with increasing normal stress. Additionally, post-oscillation wave velocity and permeability exhibit quick recoveries toward pre-oscillation values. Estimates of fracture contact area (global and local) demonstrate that the elastodynamic and permeability responses are dominated by fracture topology.
Pressure-permeation experiments were performed with the concentration-clamp cell. Streaming potentials and hydraulic permeabilities were measured for an AMF C-103 cation-exchange membrane bounded by 0.1 N NaCl solutions. The streaming potential calculated from the slope of the recorded potential differences versus the applied pressure, yields a value of 1.895 millivolt/dekabar. When comparison with other membranes of similar characteristics could be made, good agreement was found. The values of the hydraulic permeability varied somewhat with the applied pressure difference and are between 1.3 x 10 to the minus 8th power and 3.9 x 10 to the minus 8th power sq cm/dekabar-sec. The specific hydraulic permeabilities were also calculated and compared with data from the literature. Fair agreement was found. The diffusion coefficient of the chloride ion in the AMF C-103 membrane was calculated, using Fick's first law of diffusion based on ion concentrations calculated from the Donnan equilibrium concentration of Cl(-).
Reduced activation ferritic-martensitic (RAFM) steels have been developed for decades for use as fusion blanket structural materials, and have advantages in both mechanical properties and irradiation resistance following careful engineering of the microstructure. However, the hydrogen isotope behavior in these proposed fusion structural materials is not well understood, but is important to assess since it impacts the fusion reactor safety and self-sufficient tritium fuel cycle. Here, we investigated deuterium transport and retention in representative advanced RAFM steels, including castable nanostructured alloys (CNAs), and oxide-dispersion-strengthened (ODS) steels. A gas-driven permeation (GDP) system was used to measure the permeability, diffusivity and solubility of the studied materials, covering the temperature range from 623 K to 873 K, and the loading pressures from 1.8 x 10 4 to 1.0 x 10 5 Pa. The results indicated that the deuterium permeability has little material dependence. In contrast, the deuterium diffusivity of the studied materials showed significant variation. The deuterium diffusivity in ODS steels is one order of magnitude lower than that in RAFM steels and CNAs, and correspondingly, have an effective solubility that is 2–10 times larger than RAFM steels and CNAs. In addition, thermal desorption spectroscopy (TDS) measurements were performed to assess the deuterium retention and desorption of these materials following a static thermal deuterium charging at 723 Kfor 1 hour under the deuterium pressure of 1.0 x 10 5 Pa. It was found that ODS steels exhibit the highest deuterium retention and have broader desorption peaks. Microstructural features contributing to deuterium retention and impacting deuterium transport are discussed to rationalize the observed deuterium behavior in the studied RAFM steels.