Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pore pressure”

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 415 records · Page 23

DOE Award No.: DE-FE0023919 Phase 4 Scientific/Technical Report

This is the Phase 4 Report for the ‘Deepwater Methane Hydrate Characterization and Scientific Assessment or Genesis of Methane Hydrates in the Gulf of Mexico (GOM2)’ research project (DOE Award No. DE-FE0023919). The report summarizes activities from October 1, 2019 to September 30, 2020. The project is led by the University of Texas at Austin (UT). The project objective is to gain insight into the nature, formation, occurrence and physical properties of methane hydrate-bearing sediments for the purpose of methane hydrate resource appraisal through the planning and execution of drilling, coring, logging, testing and analytical activities that assess the geologic occurrence, regional context, and characteristics of marine methane hydrate deposits in the Gulf of Mexico outer continental shelf (OCS). We published a dedicated American Association of Petroleum Geologists Bulletin volume describing initial results from the UT-GOM2-1 expedition in Sept. 2020. This is part 1 of a multi-volume commitment by AAPG to this project. We further confirmed that the natural gas in hydrate at GC-955 was formed by primary microbial processes (>76.1 %). The in-situ effective permeability hydrate-bearing sandy silts at the GC-955 reservoir ranges from 0.1 md (1.0×10-16 m2) to 2.4 md (2.4×10-15m2) in cores with 83% to 93% hydrate saturation. The intrinsic permeability (the single phase permeability) is estimated from reconstituted samples to be ~12 md (1.2×10-14 m2) to ~41 md (4.1×10-14 m2). We used observation and models to interpret that the core degradation that is found in pressure cores is due to dissociation of the methane hydrate in the outer circumference of the core and dissolution of that methane into the fresh pore water that the core is stored with. We are designing approaches to minimize this core loss in the future. We spent an enormous amount of effort to further improve the ability of the pressure coring tool (the PCTB) to pressure seal correctly. We completed upgrading the upper section of the PCTB to address poor pressure. We successful tested the modifications at Geotek’s test facility in Salt Lake City (Bench Test II). We completed a Land Test of the PCTB at the Schlumberger Cameron Test and Training Facility (CTTF). The tool did not seal in 6 out of 7 tests and we clearly demonstrated that cuttings were wedging in the ball valve assembly, keeping the ball valve from sealing. We reproduced the failure mechanism observed during the land test at Salt Lake City and confirmed the sensitivity of the ball valve assembly to grit. Geotek designed and tested 9 modifications to address this issue and the PCTB is now 100% successfully sealing in the presence of grit. Our science expedition is scheduled for spring 2022 and we are fully focused on preparing for this. UT and Ohio State completed a Shallow Hazard Assessment report for each proposed UT-GOM2-2 drilling location, pursuant to 30 CFR 250.214(f) and 250.244 (f). The Shallow Hazard Reports will accompany the UT-GOM2-2 Exploration Plan that is submitted to BOEM, and completes the geological and geophysical analysis for UT-GOM2-2 permitting efforts. We updated the UT-GOM2-2 Operations Plan (Version 1). We completed the UT-GOM2-2 Science and Sample Distribution Plan (Version 1). We evaluated the scope, budget, and schedule that would result from using a commercial vessel. We developed detailed drilling schedule, mud volume, and resource estimates. We developed a vessel specification document and a well plan, and sent these documents to prospective vessel contractors.

03 NATURAL GAS↗

Water transport in reverse osmosis membranes is governed by pore flow, not a solution-diffusion mechanism

We performed nonequilibrium molecular dynamics (NEMD) simulations and solvent permeation experiments to unravel the mechanism of water transport in reverse osmosis (RO) membranes. The NEMD simulations reveal that water transport is driven by a pressure gradient within the membranes, not by a water concentration gradient, in marked contrast to the classic solution-diffusion model. We further show that water molecules travel as clusters through a network of pores that are transiently connected. Permeation experiments with water and organic solvents using polyamide and cellulose triacetate RO membranes showed that solvent permeance depends on the membrane pore size, kinetic diameter of solvent molecules, and solvent viscosity. This observation is not consistent with the solution-diffusion model, where permeance depends on the solvent solubility. Motivated by these observations, we demonstrate that the solution-friction model, in which transport is driven by a pressure gradient, can describe water and solvent transport in RO membranes.

42 ENGINEERING↗

Pore size effect on selective gas transport in shale nanopores

In shale gas production, gas composition may vary over time. To understand this phenomenon, we use molecular dynamics simulations to study the permeation of CH 4 , C 2 H 6 and their mixture from a source container through a pyrophyllite nanopore driven by a pressure gradient. For a pure gas, the flow rate of CH 4 is always higher than that of C 2 H 6 , regardless of pore size. For a 1:1 C 2 H 6 : CH 4 mixture, however, C 2 H 6 :CH 4 flow rate ratio is higher than the compositional ratio in the container (i.e., 1:1) when the pore size is smaller than ~1.8 nm. The selective transport is caused by the competitive adsorption of C 2 H 6 over CH 4 in the nanopore. The selectivity is also determined by the interplay between the surface diffusion of the adsorbed molecules and the viscous flow in the center of the pore, and it diminishes as the viscous flow becomes to dominate the surface diffusion when the pore size becomes larger than 1.8 nm. Our work shows that compositional differentiation of shale gas in production is a consequence of nanopore confinement and therefore a key characteristic of an unconventional reservoir. Finally, the related compositional information can potentially be used for monitoring the status of a production well such as its recovery rate.

03 NATURAL GAS↗

Densification Pressure Optimization of MOF-808-Based Membranes for Lithium Metal Batteries

The use of metal–organic frameworks (MOFs) in hybrid electrolytes for lithium (Li) metal batteries has grown in prominence in recent years, primarily due to the chemical tunability of the MOF’s pore structures, which can directly influence Li–ion transport properties. The most attractive form factor for a MOF electrolyte is a thin, flexible membrane, which requires the application of pressure to increase the contact between the MOF particles. Herein, a systematic study of the influence of pressure on the properties of MOF-808-based membranes is presented. It is shown that when a dry, roll-pressed membrane is subjected to pressure ≥120 MPa, a total loss of crystallinity and a significant loss of porosity is observed. Alternatively, a slurry-cast membrane, compressed under controlled pressures, can maintain crystallinity and porosity while decreasing the interparticle void space. Interesting, the conductivity of the membranes infiltrated with liquid electrolyte is not greatly affected by the pressure applied, though ultimately it is shown that for cycling with Li metal, compressed membranes with compact particles are preferred. In conclusion, this study highlights the critical importance of controlling the pressure applied to MOF-based membranes during fabrication and during cell assembly and lays out the foundation for further investigation of how to optimize membrane fabrication for hybrid electrolytes that use MOFs as the dominate component.

25 ENERGY STORAGE↗

Calculation of Hydraulic Conductivity for a Porous Medium with Staggered Array of Solid Grains - 20087

For a porous medium with periodic array of solid grains. the flow of an incompressible viscous fluid in the pore space is considered with an aim of calculating the permeability of the medium. The solid grains are of ellipsoidal shape. For this purpose, a unit-cell boundary-value problem (Stokes problem) is solved numerically. With the assumption of periodicity of the medium structure and the fluid properties (in fact constant), the flow field and pressure distributions are determined numerically with inertial effects ignored. By applying the theory of homogenization to the viscous flow through pore space, a microscale boundary-value problem is defined. After solving the problem, the permeability is determined by taking the micro-cell average of the velocity components, Various porous structures are chosen to allow a wide range of porosity values. The permeability increases mildly for porosity values in the lower range. But it increases sharply for porosity values in the upper range of the porosity. The permeability of the present study agrees with the Kozeny-Carman relation excellently for smaller porosity. Certain discrepancy between the calculated permeability in this study and the Kozeny-Carman relation appears and increases with porosity for larger porosity. It is believed that the discrepancy is due to the replacement of the pores by a bundle of capillaries, which is unrealistic, in the Kozeny-Carman relation. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Injection data analysis using material balance time for CO 2 storage capacity estimation in deep closed saline aquifers

Estimating the ultimate storage capacity of deep saline aquifers is important to address the formation potential to store the envisioned large volumes of CO 2 . Injection data (i.e. injection rate, bottomhole pressure, and cumulative injected volume of CO 2 ) are routinely recorded during storage operations. These data contain valuable information on the subsurface (e.g. the reservoir pore volume and the formation storage capacity) that can be extracted. In this paper, we present a two-step graphical technique to infer the pore volume and the ultimate storage capacity of closed saline aquifers by analyzing the available injection data. First, the pore volume is inferred through adapting the concept of the material balance time. Material balance time is an approximate superposition time function developed to interpret production data from oil and gas wells operating at variable pressure/rate conditions during the boundary-dominated flow period. Using material balance techniques, the ultimate storage capacity is then estimated through linear extrapolation of the average pressure trend to the maximum allowable pressure the formation can withstand. The average pressure is not available in practice, but is can be obtained from the injection data. Two approaches are presented in this study to calculate the average pressure; namely the rigorous and the approximate approaches. Unlike the rigorous approach, the approximate approach does not require a prior knowledge of some reservoir properties (e.g. relative permeability, absolute permeability, formation porosity and thickness) to calculate the average pressure. To investigate its potential and reliability in analyzing CO 2 injection data, the proposed technique is applied to four synthetic cases representing different well operating conditions. Results indicate that the approximate approach consistently overestimates the actual (simulated) storage capacity as compared to the rigorous approach. The agreement - between the inferred and the simulated reservoir pore volume, and between the analytical and numerical estimates of storage capacity - validates the potential application of the technique to CO 2 storage in closed saline aquifers. The technique is further substantiated through application to a field data set utilized from a commercial-scale geological storage (CGS) project. Finally, field data interpretation shows that the proposed technique can be utilized to identify the degree of hydraulic continuity and reservoir compartmentalization within a target formation by interpreting the corresponding pressure and rate responses.

02 PETROLEUM↗

Numerical investigation and experimental validation of water condensation in the gas diffusion layer with different properties

Liquid water in porous media hinders the transport of reactant to catalyst layer, where the electrochemical reactions occur, which in turn, reduces performance and lifetime of a proton exchange membrane fuel cell (PEMFC). In this study, two types of gas diffusion layer (GDL) materials are used to study the two-phase saturation distribution in a PEMFC. The correlations between the effective and anisotropic transport properties and porosity of the two GDLs are solved by the pore scale model. In addition, the anisotropic liquid water permeability and the relationship between saturation pressure and capillary pressure is determined by using the Lattice-Boltzmann method. Fuel cell performance as well as liquid water distribution in the GDL using neutron radiography are obtained for validation. Under wet conditions, the cell with Freudenberg GDL performs better than that using Toray GDL, especially at high current densities. The results of the two-phase model simulation show that the peak water saturation for Toray GDL occurs in the catalyst layer and can reach 40%, while the peak water saturation for Freudenberg GDL occurs inside the GDL directly under the ribs and only reaches 25% saturation. The combined results provide key insights to enable high power density operation of a PEMFC through GDL material optimization.

42 ENGINEERING↗

Reaction process involving capillary condensation within a microporous catalyst

Described herein is a catalytic reaction process including introducing one or more gas-phase reactants into a reactor comprising a microporous catalyst having a pore size less than or equal to 2 nm and adjusting the temperature and/or the pressure of the reactor such that one or more of the gas-phase reactants condense within the micropores of the catalyst thereby causing the catalytic reaction to take place in a liquid phase. Additionally, a process for engineering defects on a carboxylate-based metal organic framework (MOF) catalyst is described. The process includes providing a carboxylate-based MOF catalyst; and heating the carboxylate-based MOF catalyst in an inert gas atmosphere at temperatures between about 150° C. and about 900° C.

Aguirrezabal, Iker↗

Melting of Magnesium Borohydride under High Hydrogen Pressure: Thermodynamic Stability and Effects of Nanoconfinement

The thermodynamic stability and melting point of magnesium borohydride were probed under hydrogen pressures up to 1000 bar (100 MPa) and temperatures up to 400 °C. At 400 °C, Mg(BH4)2 was found to be chemically stable between 700 and 1000 bar H2, whereas under 350 bar H2 or lower pressures, the bulk material partially decomposed into MgH2 and MgB12H12. The melting point of solvent-free Mg(BH4)2 was estimated to be 367-375 °C, which was above previously reported values by 40-90 °C. Our results indicated that a high hydrogen backpressure is needed to prevent the decomposition of Mg(BH4)2 before measuring the melting point and that molten Mg(BH4)2 can exist as a stable liquid phase between 367 and 400 °C under hydrogen overpressures of 700 bar or above. The occurrence of a pure molten Mg(BH4)2 phase enabled efficient melt-infiltration of Mg(BH4)2 into the pores of porous templated carbons (CMK-3 and CMK-8) and graphene aerogels. Both transmission electron microscopy and small-angle X-ray scattering confirmed efficient incorporation of the borohydride into the carbon pores. The Mg(BH4)2@carbon samples exhibited comparable hydrogen capacities to bulk Mg(BH4)2 upon desorption up to 390 °C based on the mass of the active component; the onset of hydrogen release was reduced by 15-25 °C compared to the bulk. Importantly, melt-infiltration under hydrogen pressure was shown to be an efficient way to introduce metal borohydrides into the pores of carbon-based materials, helping to prevent particle agglomeration and formation of stable closo-polyborate byproducts.

White, James L.↗

High Capacity Step-Shaped Hydrogen Adsorption in Robust, Pore-Gating Zeolitic Imidazolate Frameworks (Final Technical Report)

The use of porous adsorbents for the densification of H 2 under relatively mild pressures and temperatures has been long pursued through various avenues of synthetic organic frameworks, namely metal–organic frameworks. However, the overwhelming majority of these materials are macroscopically rigid, exhibiting adsorption–desorption profiles that require appreciable energetic input for their full payload capacity to be stored and subsequently delivered. This reduced “deliverable” capacity is a core inefficiency in the use of porous adsorbents for H 2 densification. Our core mission is to determine a material-based strategy that fundamentally alters the adsorption¬–desorption profile in a manner that obviates these deleterious energetic requirements. Towards this, we have focused on metal–organic frameworks that exhibit “cooperative flexibility”, wherein the material undergoes reversible macroscopic phase changes in response to changes to external temperature and/or adsorbate pressure (i.e., concentration). Where these reversible switches occur between states of disparate enough accessible porosity, cooperatively flexible frameworks can exhibit “step-shaped” adsorption–desorption profiles that enable use of the delivery gaseous payload without the aforementioned additional energy requirements. However, before our work, such the observation of this phenomenon with very weakly adsorbing H 2 had not been observed with the necessary profile and P/T conditions desired for H 2 storage and delivery. During this project we have leveraged a deep understanding of the phase change in family of cooperatively flexible metal–organic frameworks to systematically alter their adsorption–desorption profiles to achieve step-shaped adsorption and desorption of H 2 in a precise P/T regime that enables a high deliverable capacity. Specifically, the thermodynamics of the adsorption/desorption induced phase changes in the baseline framework CdIF-13 (sod–Cd(benzimidazolate) 2 ) were tuned using the mixed-linker, or multivariate, approach to framework derivatization. Follow this, we have pursued a similar approach to tuning the adsorption–desorption profile of a second framework family, MIL-53(Al), which contains only Al(III) and theoretically exhibits a much higher deliverable capacity. In total, our work demonstrates that by first understanding the structural consequences of phase changes in these materials, their behavior can be systematically altered through ligand substitutions, tuning the pressure and temperature conditions at which H 2 is adsorbed and desorbed. Thus, enabling high deliverable capacities of H 2 with minimized energetic requirements.

08 HYDROGEN↗

Multiaxis sensing system for detection of airborne molecules and radionuclides

Various technologies presented herein relate to a method and equipment for detecting both airborne radioisotope and molecular effluent gases. Multi-Axis Sensing can be conducted by utilizing a pressurized MOF sorbent, such as a scintillating Metal-Organic Frameworks (S-MOFs). These MOFs are crystalline nanoporous materials that have synthetic versatility that allow adjustment of pore size, chemical environment, and luminescence properties. A method for detecting an analyte in a fluid sample is provided that comprises: loading a sorbent with a sample fluid, wherein the sorbent comprises a MOF material; pressurizing the sample fluid to increase the fluid in the sorbent thereby making a pressurized sorbent; and detecting ionizing radiation or a chemical property of the analyte in the pressurized sorbent.

Allendorf, Mark D.↗

Investigation of the Effect of Framework Flexibility on CO 2 Adsorption in SIFSIX-3-Cu Using a Machine-Learned Force Field

Metal–organic frameworks (MOFs) offer promise as selective CO 2 sorbents, but successful MOF sorbent materials need high CO 2 binding affinity and selectivity for CO 2 over water. This work focuses on the use of machine-learned force fields (MLFFs) to model CO 2 adsorption in flexible MOFs, with a focus on SIFSIX-3-Cu, an anion-pillared MOF known for its high CO 2 affinity. A preliminary high-throughput screening of over 900 anion-pillared MOFs was performed using rigid UFF+DDEC6 force fields to predict zero-loading heats of adsorption for CO 2 and H 2 O. SIFSIX-3-Cu was selected for further computational study due to its predicted CO 2 heat of adsorption and experimental relevance. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies and forces, with an iterative sampling scheme combining molecular dynamics, geometry optimization, random geometric insertion, and NVT Monte Carlo-based configuration generation to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included, contrasting with previous models that approximated the MOF as rigid. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations with the MLFF produced CO 2 adsorption isotherms in good agreement with experimental data at direct air capture (DAC) pressures (e.g., 40 Pa), in contrast to previous overestimations of CO 2 sorption by models with rigid structures. Bond and angle histogram analysis showed that MOF flexibility increased the variance of fluorine–fluorine diagonal distances at adsorption sites, resulting in a lower predicted sorption for flexible, asymmetric SIFSIX-3-Cu pore geometries compared to the rigid, symmetric DFT-optimized SIFSIX-3-Cu pore geometry. A detailed description of flexibility afforded by the MLFF resulted in an accurately predicted CO 2 uptake (0.88 mmol/g) at low pressure (40 Pa) compared to the experimentally measured value (1.24 mmol/g). In conclusion, these results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low-pressure applications.

adsorption↗

Analysis of gas storage and transport in Eagle Ford shale using pressure pulse decay measurements with He, Kr and CO 2

Detailed characterization of gas transport and storage in low-permeability and organic-rich shales is associated with an array of challenges due to the complex morphology of the pore space, representing a broad range of pore sizes, combined with the heterogeneous fabric of the shale matrix. These factors, and their interplay during gas transport and sorption, complicate a) the analysis of shale samples at the laboratory scale (~ up to a ft. In length), and b) the prediction of natural gas production and carbon sequestration potential at larger scales. Here, in this work, tandem experiments with inert (helium – He) and adsorbing (krypton - Kr and carbon dioxide - CO 2 ) gases were performed and analyzed to develop an efficient workflow for characterizing and modeling transport and sorption at the laboratory scale. In particular, pressure pulse-decay (PPD) measurements were conducted on an Eagle Ford shale core sample at room temperature using He, Kr, and CO 2 . PPD measurements with He (a non-sorbing gas) were employed to probe the overall porosity, including natural fractures, microcracks, mesopores, and micropores. A triple-porosity model (TPM) was adopted to interpret the gas transport in the shale sample: The pore space is represented by three interacting continua, including macropores (larger fractures), mesopores (including microcracks), and micropores. To facilitate the application of the TPM, a modified analytical approach is introduced to extract effective transport parameters (in terms of the characteristic time for transport at relevant porosity levels) directly from the PPD measurements with He. To validate the modeling approach, model parameters extracted from two He PPD measurements are demonstrated to provide excellent agreement with a 3rd He PPD measurement performed at a higher pressure. The effective transport parameters, extracted from the He experiments, were subsequently converted for application to the Kr and CO 2 PPD experiments, by accounting for relevant transport modes and differences in fluid properties. Excess adsorption isotherms were extracted from the equilibrium pressures of the Kr and CO 2 PPD experiments using the He pore volume as a baseline. These adsorption isotherms were then integrated into the TPM to predict combined gas transport and sorption for Kr and CO 2 , with transport coefficients translated from the He measurements. The predictions for Kr and CO 2 are demonstrated to be in excellent agreement with the experimental observations. This, in turn, demonstrates that the proposed analytical approach provides for an effective characterization of mass transfer rates in shales, that can be applied directly in a TPM representation of mass transfer and sorption.

58 GEOSCIENCES↗

Small-Angle Neutron Scattering Investigation of Oil Recovery in Mineralogically Distinct Wolfcamp Shale Strata

Understanding and improving hydrocarbon yields during enhanced oil recovery (EOR) in unconventional reservoirs is complicated by the intrinsic mineralogical and geochemical heterogeneity of shale formations. Here, in this study, we utilized small-angle neutron scattering (SANS) and ultra-small-angle neutron scattering (USANS) to investigate the degree of oil retention and its location in the nanoporous shale matrix for two mineralogically distinct shale samples. The two samples, dubbed “dark” and “light” based on their color, were taken from adjacent strata in a Wolfcamp shale core. While both samples contained kerogen, the dark sample contained more kerogen and clay (43.7 wt %) while the light sample contained more calcite (54.9 wt %). Samples were presaturated with decane, a model hydrocarbon, prior to pressure cycling with methane. Results showed significantly more retention of decane in 1.5–10 nm radius pores of both, likely indicating that oil is retained within kerogen nanopores. Although the dark sample had a higher porosity of 8.7%, versus 3.3% for the light sample, more pores were accessible to decane and a higher percentage of the imbibed decane was removable from the light sample compared to the dark sample. The majority of decane was not recoverable for the dark sample, indicating that EOR with methane can be challenging. These new findings can help to model expected recoveries of in-place oil from heterogeneous shale formations, as well as inform improved EOR strategies.

04 OIL SHALES AND TAR SANDS↗

Localized growth and remodelling drives spongy mesophyll morphogenesis

The spongy mesophyll is a complex, porous tissue found in plant leaves that enables carbon capture and provides mechanical stability. Unlike many other biological tissues, which remain confluent throughout development, the spongy mesophyll must develop from an initially confluent tissue into a tortuous network of cells with a large proportion of intercellular airspace. How the airspace in the spongy mesophyll develops while the tissue remains mechanically stable is unknown. Here, we use computer simulations of deformable polygons to develop a purely mechanical model for the development of the spongy mesophyll tissue. By stipulating that cell wall growth and remodelling occurs only near void space, our computational model is able to recapitulate spongy mesophyll development observed in Arabidopsis thaliana leaves. We find that robust generation of pore space in the spongy mesophyll requires a balance of cell growth, adhesion, stiffness and tissue pressure to ensure cell networks become porous yet maintain mechanical stability. The success of this mechanical model of morphogenesis suggests that simple physical principles can coordinate and drive the development of complex plant tissues like the spongy mesophyll.

Science & Technology - Other Topics↗

The physical basis for solvent flow in organic solvent nanofiltration

Organic solvent nanofiltration (OSN) is an emerging membrane technology that could revolutionize chemical separations in numerous vital industries. Despite its significance, there remains a lack of fundamental understanding of solvent transport mechanisms in OSN membranes. Here, we use an extended Flory-Rehner theory, nonequilibrium molecular dynamic simulations, and organic solvent transport experiments to demonstrate that solvent flow in OSN membranes is driven by a pressure gradient. We show that solvent molecules migrate as clusters through interconnected pathways within the membrane pore structure, challenging the widely accepted diffusion-based view of solvent transport in OSN. We further reveal that solvent permeance is dependent on solvent affinity to the OSN membrane, which, in turn, controls the membrane pore structure. Our fundamental insights lay the scientific groundwork for the development of next-generation OSN membranes.

42 ENGINEERING↗

The Role of Surface Hydrophobicity on the Structure and Dynamics of CO2 and CH4 Confined in Silica Nanopores

Advancing a portfolio of technologies that range from the storage of excess renewable natural gas for distributed use to the capture and storage of CO 2 in geological formation are essential for meeting our energy needs while responding to challenges associated with climate change. Delineating the surface interactions and the organization of these gases in nanoporous environments is one of the less explored approaches to ground advances in novel materials for gas storage or predict the fate of stored gases in subsurface environments. To this end, the molecular scale interactions underlying the organization and transport behavior of CO 2 and CH 4 molecules in silica nanopores need to be investigated. To probe the influence of hydrophobic surfaces, a series of classical molecular dynamics (MD) simulations are performed to investigate the structure and dynamics of CO 2 and CH 4 confined in OH-terminated and CH 3 -terminated silica pores with diameters of 2, 4, 6, 8, and 10 nm at 298 K and 10 MPa. Higher adsorption extents of CO 2 compared to CH 4 are noted on OH-terminated and CH 3 -terminated pores. The adsorbed extents increase with the pore diameter. Further, the interfacial CO 2 and CH 4 molecules reside closer to the surface of OH-terminated pores compared to CH 3 -terminated pores. The lower adsorption extents of CH 4 on OH-terminated and CH 3 -terminated pores result in higher diffusion coefficients compared to CO 2 molecules. The diffusivities of both gases in OH-terminated and CH 3 -terminated pores increase systematically with the pore diameter. The higher adsorption extents of CO 2 on OH-terminated and CH 3 -terminated pores are driven by higher van der Waals and electrostatic interactions with the pore surfaces, while CH 4 adsorption is mainly due to van der Waals interactions with the pore walls. These findings provide the interfacial chemical basis underlying the organization and transport behavior of pressurized CO 2 and CH 4 gases in confinement.

Mohammed, Sohaib↗