Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geologic”

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 1,369 records · Page 76

Transformational faulting in Mn 2 GeO 4 from olivine to wadsleyite structure: Implications for physical mechanism of deep-focus earthquakes

High-pressure and temperature deformation experiments interfaced with acoustic emission (AE) monitoring have been conducted to study transformational faulting in Mn 2 GeO 4 olivine, which transforms to the β phase, isostructural to wadsleyite. Metastable Mn 2 GeO 4 olivine exhibits a marked embrittlement behavior at temperatures between 800 and 1100 K, emitting numerous AEs. At each temperature, brittle deformation is characterized by a two-stage process: (1) a “preparation” stage with numerous diffusedly located low-magnitude AEs and large b values (>2), and (2) a failure stage where larger-magnitude AEs form a planar distribution with b values about 1. Microstructure analysis reveals extensive kink band development in olivine grains in the recovered samples. Kink band boundaries (KBBs), with a typical thickness of ∼100 nm, are filled with a nanometric β-Mn 2 GeO 4 “gouge”. A dense array of secondary shear localizations is often present within the kink bands, suggesting significant shear deformation therein. The combined observations suggest that faulting in metastable Mn 2 GeO 4 olivine is a self-similar process, from grain-scale to the sample-scale. Both observed embrittlement behavior and the microstructure of metastable Mn 2 GeO 4 olivine are essentially identical to those in Mg 2 GeO 4 olivine we have reported previously, indicating that the physical mechanism of faulting in metastable olivine is insensitive to the specific crystallographic structure of the high-pressure phase. The low b values (about 1) observed in the faulting process in our experiments are similar to those of deep focus earthquakes in cold subduction zones. Our observed mechanism explains deep focus seismicity in cold metastable mantle wedges, provided that the self-similarity assumption holds to geological scales.

58 GEOSCIENCES↗

Reductive biomining of pyrite by methanogens

Pyrite (FeS 2 ) is the most abundant iron sulfide mineral in Earth’s crust. Until recently, FeS 2 has been considered a sink for iron (Fe) and sulfur (S) at low temperature in the absence of oxygen or oxidative weathering, making these elements unavailable to biology. However, anaerobic methanogens can transfer electrons extracellularly to reduce FeS 2 via direct contact with the mineral. Reduction of FeS 2 occurs through a multistep process that generates aqueous sulfide (HS – ) and FeS 2 -associated pyrrhotite (Fe 1-x S). Subsequent dissolution of Fe 1-x S provides Fe(II) (aq) , but not HS – , that rapidly complexes with HS – (aq) generated from FeS 2 reduction to form soluble iron sulfur clusters [nFeS (aq) ]. Cells assimilate nFeS (aq) to meet Fe/S nutritional demands by mobilizing and hyperaccumulating Fe and S from FeS 2 . As such, reductive dissolution of FeS 2 by methanogens has important implications for element cycling in anoxic habitats, both today and in the geologic past.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coupled hydro-thermal flow and radionuclide transport driven by spatial variation of heat-generating radioactive wastes in shale formations

Deep geologic disposal of multiple nuclear waste packages with various heat sources can induce nonuniform hydro-thermal behaviors in the near-field of the repository, consequently influencing the long-term radionuclide transport in the far-field once waste form breach initiates. Here, this study looks into three cases with variation in the spatial order of six groups of heat sources (10th, 50th, 75th, 90th, 95th, and 99th percentiles of heat outputs generated from 1,981 as-loaded dual-purpose canisters in the field site) in a shale-hosted repository with respect to the uni-directional groundwater flow (from west to east): (1) cooler waste packages from west to east, (2) hotter waste packages from west to east, and (3) hottest waste packages in the middle of the repository. Our field-scale PFLOTRAN simulation represents heat-driven multiphysics coupled mechanisms, including multiphase flow, heat transfer, and chemical/radioactive transport, and also, calculates the onset of waste form breach based on temperature-dependent canister vitality. The results from this sensitivity study will quantify the short- (less than 1 × 10 3 years) and long-term (up to 1 × 10 6 years) impacts of sporadic heat pulses from waste package on the spatio-temporal perturbation in hydro-thermal flow quantities and the rate of radionuclide transport in both near- and far-field of the repository system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Transmission electron microscopy with in-situ ion irradiation: Facilities and community

Whilst there is a clear scientific and technological need for the technical capabilities of transmission electron microscopes with in-situ ion irradiation, it also requires a collaborative community of international researchers to support such facilities in successfully meeting this demand. Instruments of this type serve to provide fundamental understanding of the mechanisms which drive changes in materials important to nuclear fission and fusion energy, the semiconductor industry, quantum information systems, space travel, astronomy, geology and many more applications. As these areas continue to evolve and the instrumentation possibilities expand, the capacity of in-situ ion irradiation facilities must also develop hand-in-hand with the user community to deliver an ever-greater diversity of high-fidelity extreme-environment experimentation. Future directions for the field, such as miniaturization from MEMS/microfluidic devices and advanced controls with ML-based analysis, continuously emerge to advance both the hardware and software which support the coupling of TEMs with ion beams. This review sets out to provide up-to-date insights into the community and advancement of current, and development of future, facilities which have the potential to further unlock access to the nanoscale exploration of coupled extreme environments crucial to many of the important science and engineering challenges we face today.

In-situ irradiation↗

Plutonium reactive transport in fractured granite: Multi-species experiments and simulations

Plutonium (Pu) in the subsurface environment can transport in different oxidation states as an aqueous solute or as colloidal particles. The transport behavior of Pu is affected by the relative abundances of these species and can be difficult to predict when they simultaneously exist. This study investigates the concurrent transport of Pu intrinsic colloids, Pu(IV) (aq) and Pu(V-VI) (aq) through a combination of controlled experiments and semi-analytical dual-porosity transport modeling. Pu transport experiments were conducted in a fractured granite at high and low flow rates to elucidate sorption processes and their scaling behavior. In the experiments, Pu(IV) (aq) was the least mobile of the Pu species, Pu(V-VI) (aq) had intermediate mobility, and the colloidal Pu, which consisted mainly of precipitated and/or hydrolyzed Pu(IV), was the most mobile. The semi-analytical modeling revealed that the sorption of each Pu species was rate-limited, as the sorption could not be described by assuming local equilibrium in the experiments. The model was able to describe the sorption of the different Pu species that occurring either on fracture surfaces, in the pores of the rock matrix, or simultaneously in both locations. While equally good fits to the data could be achieved using any of these assumptions, a fracture-dominated process was considered to be the most plausible because it provided the most reasonable estimates of sorption rate constants. Importantly, a key result of this work is that the sorption rate constant of all Pu species tends to decrease with increasing time scales, which implies that Pu will tend to be more mobile at longer time scales than observations at shorter time scales suggest. This result has important implications for predicting the environmental impacts of Pu in the safety assessments of geologic repositories for radioactive waste disposal, and we explore potential mechanistic bases for upscaling the sorption rate constants to time and distance scales that cannot be practically evaluated in experiments.

58 GEOSCIENCES↗

Parts-Per-Million Carbonate Mineral Quantification with Thermogravimetric Analysis–Mass Spectrometry

Mitigating the deleterious effects of climate change requires the development and implementation of carbon capture and storage technologies. To expand the monitoring, verification, and reporting (MRV) capabilities of geologic carbon mineralization projects, we developed a thermogravimetric analysis–mass spectrometry (TGA–MS) methodology to enable quantification of <100 ppm calcite (CaCO 3 ) in complex samples. We extended TGA–MS calcite calibration curves to enable a higher measurement resolution and lower limits of quantification for evolved CO 2 from a calcite–corundum mixture. We demonstrated <100 ppm carbonate mineral quantification with TGA–MS for the first time, an outcome applicable across earth, environmental, and materials science fields. We applied this carbonate quantification method to a suite of Columbia River Basalt Group (CRBG) well cuttings recovered in 2009 from Pacific Northwest National Laboratory’s Wallula #1 Well. Our execution of this new combined calcite and calcite–corundum calibration curve TGA–MS method on our CRBG sample suite indicated average carbonate contents of 0.050 wt % in flow interiors (caprocks) and 0.400 wt % in interflow zones (reservoirs) in the upper 1250 m of the Wallula #1 Well. Finally, by advancing our knowledge of continental flood basalt-hosted carbonates in the mafic subsurface and reaching new TGA–MS quantification limits for carbonate minerals, we expand MRV capabilities and support the commercial-scale deployment of carbon mineralization projects in the Pacific Northwest United States and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Single-Particle Isotopic Analysis Using the Liquid Sampling-Atmospheric Pressure Glow Discharge Microplasma Coupled to an Orbitrap Mass Spectrometer

Isotope ratio (IR) determinations on individual particles can provide valuable information relative to the sample, including processing and dating, which could be valuable to geological and nuclear material analysis communities. In comparison to “bulk” measurements, in which all particle information is lost and homogenized within a sample, “particle” measurements allow for fingerprinting and can provide unique insights related to the sample’s nano- and micro- compositions. Furthermore, the complex nature of particles, with potential isobaric interferences from either the matrix or other elements within the same particle, poses a significant analytical challenge for conventional mass spectrometry platforms employed for elemental analysis due to their limited mass resolution. Presented here is the first demonstration of an ultrahigh resolution method for single particle (SP) isotopic analysis using the liquid sampling-atmospheric pressure glow discharge (LS-APGD) microplasma ionization source coupled to an Orbitrap mass spectrometer and FTMS Booster X2T acquisition/processing system. For proof of concept, well-characterized CeO 2 microparticles with a nominal diameter of 1.0 μm were analyzed at a mass resolution of ∼330,000. Important instrumental parameters were optimized to enable the detection of single particles. The 142 Ce/ 140 Ce ratio of individual particles and the population average were determined and compared to an SP inductively coupled plasma time-of-flight mass spectrometry (ICP-TOF-MS) analysis. The isotopic accuracy and precision were determined by two methods and found to be in good agreement with the SP-ICP-TOF-MS method, with the linear regression slope method providing a more accurate ratio, showing a ∼1.4% relative difference from the ratio determined from a particle digest. The encouraging performance of the method was further supported by a determined detection limit of 2.9 fg ( 142 Ce). The developed method is anticipated to overcome many of the challenges posed by isobaric interferences encountered in conventional mass spectrometry techniques when analyzing real-world samples for particle isotopic composition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Opposing Effects of Impurity Ion Sr 2+ on the Heterogeneous Nucleation and Growth of Barite (BaSO 4 )

Nucleation of barite (BaSO 4 ) has broad implications in geological, environmental, and materials sciences. While impurity metals are common, our understanding of how they impact nucleation remains dim. In this work, we used classical optical microscopy compared to fast X-ray nanotomography (XnT) to investigate heterogeneous nucleation of barite on silica in situ with Sr 2+ as an impurity ion. The observed barite nucleation rates were consistent with classical nucleation theory (CNT), where barite crystals displayed a nonuniform size distribution, exhibiting distinct morphologies and incubation periods in Sr-free solutions. While undetectable with optical microscopy, nanotomography revealed that addition of Sr2+ enhances nucleation rates driven by the pre-factor in CNT, likely because both adsorbed Ba 2+ and Sr 2+ act as precursor sites on which nucleation occurs. Sr 2+ simultaneously inhibits growth, however, leading to a homogeneous distribution of smaller crystals. This finding will enable an improved predictive understanding of nucleation in natural and synthetic environments.

36 MATERIALS SCIENCE↗

Recent Advances in Nonclassical Crystallization: Fundamentals, Applications, and Challenges

We report crystals are ubiquitous and play an important role in our daily life. Recent studies on properties and formation of crystals have proven the attachment of particle ranging from ion pairs to well-crystallized nanocrystal is one of the most important ways for crystallization. Compared with the monomer-by-monomer (i.e., atoms, ions, or molecules) addition described in longstanding classical theories, so-called nonclassical crystallization pathways are more complex as the non-specific interaction of variation in system free energy and reaction kinetics. In view of the detailed insights into nonclassical crystallization pathways has broadened the scope of research across various disciplines such as geological events, biomineralization mechanisms, environmental remediation, and the development of functional materials, this paper reviews the current state of the art in nonclassical nucleation (the prenucleation clusters pathway, aggregation pathway, etc.) and growth pathways (oriented attachment and random attachment) based on several spotlights. In addition, the observed specific phenomena on crystallization in the last 10 years as well as recent advancements and achievements made on the synthesis of mesocrystals are also elucidated. With implications and perspectives of nonclassical crystallization researches proposed in this paper, this review aims to summarize research advancement regarding crystallizing systems. Emphasis will be placed on novel insights into nonclassical crystallization of materials, which are expected to provide completely new perspectives on the underlying crystallization mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hybrid Double Perovskite Containing Helium: [He 2 ][CaZr]F 6

Perovskites are of great technological and geological importance, in large part, due to their considerable compositional and structural flexibility. However, the formation of perovskites with neutral species on their A-sites is very unusual. The formation, phase transitions and properties of [He 2 ][CaZr]F 6 , which is the first noble gas containing perovskite to be made, are reported. It is likely that a large family of related materials can also be prepared. On compression in neon, the negative thermal expansion (NTE) material CaZrF 6 amorphizes at ~ 0. 5 GPa. However, on compression in helium at room temperature, gas is inserted into the structure to form a perovskite with helium on the A-site. This suppresses the amorphization until > 3 GPa. Volume versus pressure and Raman measurements suggest that filling of the A-site, to give [He 2 ][CaZr]F 6 , is complete at > 1 GPa. The presence of helium on the A-site in this perovskite leads to a reduction in the magnitude of negative thermal expansion when compared to the parent phase CaZrF 6 , likely due to steric impediment of the transverse vibrational motion of fluoride. The helium also leads to considerable stiffening of the structure. At room temperature and ~2.5 GPa, the helium containing hybrid perovskite has a bulk modulus of ~ 47 GPa, whereas CaZrF 6 has a bulk modulus of ~ 40 GPa at ambient. At 15 K, the cubic perovskite [He 2 ][CaZr]F 6 undergoes a structural phase transition, that likely involves a cooperative tilting of framework octahedra to give a lower-symmetry phase, which is tentatively assigned as tetragonal.

36 MATERIALS SCIENCE↗

Nonenergy Biomass Carbon Removal and Storage (BiCRS): Assessing Durability of Nongaseous Carbon Products Across Terrestrial Storage Fates

Biomass Carbon Removal and Storage, or BiCRS, pathways use plants or algae that remove carbon dioxide from the atmosphere through photosynthesis and store it underground or in long-lived products. While some BiCRS approaches generate an energy product, all BiCRS approaches generate a carbon product. A new subset of BiCRS approaches focus on the storage of these raw or converted carbon products for generation of carbon credits. However, the durability of these approaches is highly variable as carbon products vary widely in their “form” and the conditions of their “fate.” We organize our thinking about carbon products and their durability around these two primary axes. The durability of carbon product “forms” is mediated by chemical recalcitrance and ranges substantially across agricultural residues, municipal solid waste, woody biomass, and nongaseous products of thermochemical conversion (e.g., biochars and bio-oils). Meanwhile, terrestrial storage “fates” vary in the mechanism employed to stall decay, including surface storage, dry storage, shallow anoxic storage, and deep or geologic anoxic storage (or injection). Each mechanism has different implications for suitability with different feedstock forms as well as long-term risks. We present a framework for assessing durability of solid or liquid raw and conversion carbon products under terrestrial storage fates, highlighting knowns, unknowns, and research priorities moving forward.

09 BIOMASS FUELS↗

Probing the Core–Shell Organization of Nanoconfined Methane in Cylindrical Silica Pores Using In Situ Small-Angle Neutron Scattering and Molecular Dynamics Simulations

Determining the structure of nanoconfined fluids is essential for predicting the fate of these fluids in subsurface geologic formations with nanoporous features and for engineering novel nanoporous materials for storing compressed gases. In this study, we probe the structure of nanoconfined methane at pressures in the range of 15–100 bar using in situ small-angle neutron scattering (SANS) measurements and molecular dynamics (MD) simulations. The structure of methane is probed in MCM-41 and SBA-15 with cylindrical pores and diameters of 3.3 and 6.8 nm, respectively. In situ SANS measurements and MD simulations showed that the confined methane molecules are organized in a core–shell structure, with the shell arising from the adsorption of methane molecules on the silica surface. The shell thicknesses of the adsorbed deuterated methane (CD4) molecules in MCM-41 obtained by SANS measurements are 1.4 ± 0.5, 2.1 ± 0.1, 3.2 ± 0.8, 4.0 ± 0.3, and 6.0 ± 0.7 Å at equilibrated pressures of 15.6, 35.6, 55.5, 73.3, and 95.7 bar, respectively. The shell thicknesses of the adsorbed CD4 layer in SBA-15 pores are 2.7 ± 0.5, 4.3 ± 0.7, 6.6 ± 0.7, 10.2 ± 0.3, and 14.6 ± 0.8 Å at equilibrated pressures of 15.5, 32.7, 52.4, 70, and 99.5 bar, respectively. These experimental results are in close agreement with the results predicted from MD simulations. Adsorption of methane molecules on the silica surfaces is primarily driven by van der Waals interactions between the methane molecules and the hydroxyl groups on the silica surface, while electrostatic interactions play a minor role. In conclusion, the experimental and simulation approaches described in this study provide fundamental insights into the organization of confined gases using methane as a specific example in the context of compressed fluid storage in natural and engineered materials for adaptive energy use.

03 NATURAL GAS↗

Structure and Shape of Surface-Mediated Assembly of Surfactants

Achieving controls on the self-assembly of cationic and block copolymer micelles has significant implications for advancing novel material systems for energy, environmental, and biological applications and for tuning fluid flow behavior in subsurface geologic environments. While it is known that the aggregation of micelles is influenced by the composition of the surrounding fluid, the influence of a solid interface has been less explored. In this study, we probe the organization of cetyltrimethylammonium bromide (CTAB) micelles in the absence and presence of the Pluronic P123 block copolymer and quartz substrate using transmission and grazing-incidence small-angle X-ray scattering measurements and classical molecular dynamics (MD) simulations. In the absence of the quartz interface, CTAB with and without P123 molecules assemble as ellipsoid core–shell micelles. Densification of the core on adding P123 is noted from the decrease in the core radius from 15.9 to 14.6 Å. Here, the presence of a quartz substrate causes the micelles to elongate, which is noted by the emergence of a power-law slope in the low q region (<0.02 Å –1 ). Moreover, in the presence of both P123 and the quartz substrate, the micelle shape changes from ellipsoid core–shell to cylindrical core–shell, and a significantly higher number of aggregates are formed. The higher number of aggregates, and faster aggregation kinetics are linked to the organization of the solvent structures as noted from the MD simulations. These findings demonstrate the effects of adding a block copolymer and including a quartz substrate on the self-assembly of cationic (CTAB) micelles, which can be used to inform the design of functional material systems that harness surfactant chemistries.

42 ENGINEERING↗

Confinement-Driven Heterogeneous Benzene Crystallization in Silica Nanopores

Nanoconfinement alters the thermodynamics, dynamics, and kinetics of fluids hosted in nanoscale solid nanopores to an extent that depends on the characteristics of the confining space and the chemistry of the confined fluids. Confinement-induced alterations in the phase behavior of confined energetic fluids under high pressure or low temperature are highly relevant to subsurface and subsea phenomena such as fluid flow in porous media, hydrate formation and dissociation, and gas storage capacity. Although extensive efforts have been directed toward understanding the phase behavior of confined fluids, the role of solid-liquid interfaces in the phase transitions of organic liquids has not been resolved yet. Here, we explore the onset and growth of benzene crystallization confined in 6 nm sized SBA-15 silica nanopores in the temperature range 300-200 K using in situ extended range small-angle and wide-angle X-ray scattering (SAXS/WAXS) measurements and atomistic classical molecular dynamics (MD) simulations. The crystallization onset of confined benzene depresses to 265 K compared to the freezing point of bulk benzene ( ~278 K), followed by the continuous growth of the emerged crystals in the pore space with complete crystallization at 200 K. The orientation of the emerged benzene crystals is dominated by parallel (π-π stacking) and perpendicular (T-shape stacking) orientations along the cylindrical pore radius and pore length, respectively. The onset of benzene crystals occurs heterogeneously on the pore surface and grows continuously toward the pore center. Further, confined benzene undergoes a dynamical crossover from fragile to strong dynamics behavior, inferred from the rotational and translational diffusion. The insights provided by this study have significant implications for the phase transitions of confined organic liquids that are relevant to a wide range of applications in the biological, geological, environmental, and chemical fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of Supercritical CO 2 Injection on the Shale Pore Structures and Mass Transport Rates

Characterizing the pore structures and transport properties of low-permeability shales is critical for evaluating these formations as potential seals or storage sites for geological CO 2 sequestration. Here, we use low-pressure gas adsorption in conjunction with nuclear magnetic resonance (NMR) to characterize the pore-size distribution of shales before and after injection of supercritical CO 2 . Nitrogen gas was used as the detecting phase for the adsorption experiments and pentane liquid was used for the NMR experiments. We also performed time-resolved NMR and gravimetric microbalance measurements to observe mass transport during desaturation. We use these data to estimate the self-diffusion coefficient of pentane and changes in the saturation state of the pore network. We analyzed samples with a range of compositions from the Wolfcamp shale before and after exposure to supercritical CO 2 for 3 days. Integrating the gas adsorption and NMR data shows how supercritical CO 2 injection alters the pore-size distribution for pore sizes <1 nm to 1 mm. Finally, our results provide insights on how the pore structure and mass transport properties of different shale lithologies may evolve during storage of supercritical CO 2 .

58 GEOSCIENCES↗

CO 2 -Soluble Nonionic Surfactants for Enhanced CO 2 Storage via In Situ Foam Generation

Geologic carbon storage (GCS) is a rapidly evolving technology with the potential to reduce the environmental impact of fos-sil fuels. Saline aquifers, which comprise a sandstone matrix with brine in the pores, make up much of the pore space availa-ble for CO 2 storage in the United States. When CO 2 is injected in saline aquifers, however, capillary fingering occurs and only a small percentage of the pore space is filled with CO 2 . This is due to the low viscosity of CO 2 , which is roughly ten times less viscous than brine. To address this problem, we tested the ability of inexpensive, commercially available nonionic surfactants to be dissolved in the injected CO 2 and increase the apparent viscosity of CO 2 by generating CO 2 -in-water foams in situ.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantification of the Impact of Acidified Brine on Fracture-Matrix Transport in a Naturally Fractured Shale Using in Situ Imaging and Modeling

Understanding flow, transport, chemical reactions, and hydromechanical processes in fractured geologic materials is key for optimizing a range of subsurface processes including carbon dioxide and hydrogen storage, unconventional energy resource extraction, and geothermal energy recovery. Flow and transport processes in naturally fractured shale rocks have been challenging to characterize due to experimental complexity and the multiscale nature of quantifying continuum scale descriptions of mass exchange between micrometer-scale fractures and nanometer-scale pores. In this study, we use positron emission tomography (PET) to image the transport of a conservative tracer in a naturally fractured Wolfcamp shale core before and after the core was exposed to low pH brine conditions. Image-based experimental observations are interpreted by fitting an analytical transport model to fracture-containing voxels in the core. Results of this analysis indicate subtle increases in matrix diffusivity and a slightly more uniform fracture velocity distribution following exposure to low pH conditions. These observations are compared with a multicomponent one-dimensional reactive transport model that indicates the capacity for a 10% increase in porosity at the fracture-matrix interface as a result of the low pH brine exposure. This porosity change is the result of the dissolution of carbonate minerals in the shale matrix to low pH conditions. Further, this image-based workflow represents a new approach for quantifying spatially resolved fracture-matrix transport processes and provides a foundation for future work to better understand the role of coupled transport, reaction, and mechanical processes in naturally fractured rocks.

58 GEOSCIENCES↗

Hydrogen Diffusion in Slit Pores: Role of Temperature, Pressure, Confinement, and Roughness

Diffusion of hydrogen (H 2 ) is important to understand the leakage risk and transport behavior for H 2 geologic storage. We applied molecular dynamics simulations to investigate the influencing factors of H2 diffusion in the slit pores of calcite, hematite, and quartz, owing to their abundance. It is revealed that the H2 self-diffusion coefficient increases with the temperature, regardless of the type of pore minerals. The diffusion of H 2 in the 20 nm slit pores falls into the bulk diffusion regime when the pressure is 10 MPa. The self-diffusion of H 2 decreases with pressure in all three types of slit pores, following a power law model with the exponents ranging from -0.825 to -0.964. Furthermore, the impact of confinement on H 2 diffusion is more pronounced for the slit pores with stronger interactions with H 2 -like calcite. The role of surface roughness in H 2 diffusion depends on the slit aperture. The rough surface enhances H 2 diffusion in the larger slit pores due to the enlarged effective pore space, whereas it weakens H 2 diffusion in the small slit pores due to stronger adsorption. These findings will fill the knowledge gap on the coupling effect of different factors influencing H 2 diffusion.

08 HYDROGEN↗