Search NASA⌕ Search

SEARCH · Search NASA

Results for “Pore restriction”

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 37 records · Page 2

Multiphysics Simulations of Nanofibrous High Power Targets for the Inference of Failure Modes

Fermilab s High Power Targetry Research and Development (HPT R&D) group have been developing and studying an electrospun nanofiber target concept to support the need for robust targets in future fixed target facilities. These nanofiber mats have demonstrated resistance to radiation damage, and the free motion of the individual fibers is expected to mitigate the cyclic stresses induced by a pulsed, high-power beam. To evaluate the efficacy of this concept, nanofiber mat samples have been sent by the HPT R&D group to the HiRadMat facility at CERN for prototypic thermal shock testing; the outcomes of these experiments show that the survivability of a nanofiber target depends on its construction parameters, in particular the packing density of the fibers. Samples with higher packing densities have consistently been destroyed by exposure to the HiRadMat beam, with a hole at the beam center visible, and layers of nanofibers peeled away from the center hole, whereas samples with lower densities have survived with limited damage. The exact reason for the failure of the higher density targets was unclear at the time of the original experiments, but the results of our recent multiphysics simulations which recreate the experiments support the hypothesis that the expansion and pressurization of the air inside the target after being heated by the pulsed beam is the cause; in a high density nanofiber mat, the motion of air through the pores of the mat is much more restricted, and induces a larger pressure on the fibers, blowing the mat apart. In this talk, we ll share the results of these simulations and discuss how they support this hypothesis.

Asztalos, Will↗

Iron Oxide Reduction Rate Affects Iron Product Morphology

Magnetite concentrates from different parts of the Mesabi Range, exhibiting varying tendencies to encapsulate, were reduced using thermogravimetric analysis (TGA) and a high-temperature confocal scanning laser microscope (CSLM). Gas-phase mass-transfer conditions were intentionally varied by employing two reactors with different geometries: the TGA produced significantly slower reduction, whereas the CSLM setup yielded reduction rates up to 20 times faster under mass-transfer-controlled conditions. A preliminary mechanism is proposed to explain the transition from a porous iron product to a dense encapsulating layer. Encapsulation is interpreted as the outcome of competition between pore creation, driven by the reduction rate, and pore elimination, driven by surface diffusion. When pore formation is insufficient to counteract pore coarsening, a dense iron layer develops and restricts further reduction.

encapsulation↗

Wettability of ultra-small pores of carbon electrodes by size-asymmetric ionic fluids

Recently, we studied the phase behavior of ionic fluids under confinement using the classical density functional theory within the framework of the restricted primitive model. Additionally, the theoretical results indicate that narrowing the pore size may lead to a drastic reduction in the electric double layer capacitance, while increasing the surface electrical potential would improve the ionic accessibility of micropores. In this work, we extend the theoretical investigation to systems containing size-asymmetric electrolytes that may exhibit a vapor-liquid like phase transition in the bulk phase. The effects of pore size and surface electric potential on the phase diagram and microscopic structures of the confined electrolytes were studied over a broad range of parameters. We found that decreasing the pore size or increasing the surface potential could destabilize the liquid phase in micropores, and capillary evaporation could occur regardless of the size asymmetry between cations and anions. Compared to that in a symmetric ionic system, the vapor-liquid phase separation is more likely to take place as the size asymmetry becomes more pronounced. The phase transition would alter the “accessibility” of ions to micropores and lead to coexisting micropores with different surface charge densities as identified by Monte Carlo simulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enhanced Selectivity and Stability of Finned Ferrierite Catalysts in Butene Isomerization

Abstract Designing zeolite catalysts with improved mass transport properties is crucial for restrictive networks of either one‐ or two‐dimensional pore topologies. Here, we demonstrate the synthesis of finned ferrierite (FER), a commercial zeolite with two‐dimensional pores, where protrusions on crystal surfaces behave as pseudo nanoparticles. Catalytic tests of 1‐butene isomerization reveal a 3‐fold enhancement of catalyst lifetime and an increase of 12 % selectivity to isobutene for finned samples compared to corresponding seeds. Electron tomography was used to confirm the identical crystallographic registry of fins and seeds. Time‐resolved titration of Brønsted acid sites confirmed the improved mass transport properties of finned ferrierite compared to conventional analogues. These findings highlight the advantages of introducing fins through facile and tunable post‐synthesis modification to impart material properties that are otherwise unattainable by conventional synthesis methods.

Dai, Heng↗

A refined design concept for sulfur-tolerant Pd catalyst supported on zeolite by shape-selective exclusion and hydrogen spillover for hydrogenation of aromatics

Sulfur poisoning of noble metal catalysts has been a major challenge for decades. This work demonstrates the superior sulfur tolerance of a hybrid zeolite-supported Pd catalyst for hydrogenation of tetralin containing benzothiophene. The hybrid catalyst consists of Pd supported on a large-pore acidic zeolite Y (Pd/HY) and a small-pore zeolite A without (Pd/HA) or with surface metal passivation by chemical vapor deposition of SiO 2 (SiO 2 -Pd/HA) and with potassium ion exchange (SiO 2 -Pd/KA). Pd/HY is very active for tetralin hydrogenation but quickly deactivates after exposure to benzothiophene at high concentration of 100 ppm sulfur. The SiO 2 -coated Pd/HA shows no activity for tetralin hydrogenation, but continues to activate H 2 and serve as a source of hydrogen spillover from Pd/HA since these metal sites are protected from thiophenic sulfur due to size-selective exclusion. Adding K ion-exchanged and SiO 2 -coated catalyst SiO 2 -Pd/KA to Pd/HY is even more effective for enhancing sulfur tolerance, both for tetralin hydrogenation and for isomerization of cis-decalin to trans-decalin in the presence of benzothiophene. The turnover frequencies (TOF) of all the catalysts are similar at the exposed sulfur/Pd atomic ratio of around 0.5. With further increase in sulfur/Pd ratio, rapid decline in TOF was observed on Pd/HY, but SiO 2 -Pd/KA + Pd/HY hybrid catalyst shows a significantly higher TOF than that over Pd/HY, even though SiO 2 -Pd/KA or SiO 2 -Pd/HA alone shows no TOF for tetralin hydrogenation. Even after high-dose poisoning at 400 ppm sulfur, the SiO 2 -Pd/KA + Pd/HY hybrid recovered activity more quickly than SiO 2 -Pd/HA + Pd/HY, while the latter in turn is much better than Pd/HY alone. These results point to the higher activity for hydrogen spillover from SiO 2 -Pd/KA where Pd metal sites are protected from not only thiophenic sulfur but also inorganic sulfur H 2 S by shape-selective exclusion. The present work further establishes the validity of shape-selective exclusion and hydrogen spillover in the design concept proposed for a sulfur-tolerant bimodal acidic zeolite-supported metal catalyst for hydrogenation of naphthalene [C.S. Song, Chemtech, 29 (1999) 26–30], and refines the design concept further for limiting the size of small pore in hybrid catalysts to~3 Å for restricting access of H 2 S but allowing entrance of H 2 . Further, this refined design concept is applicable to other supported monometallic and bimetallic catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enhanced Selectivity and Stability of Finned Ferrierite Catalysts in Butene Isomerization

Designing zeolite catalysts with improved mass transport properties is crucial for restrictive networks of either one- or two-dimensional pore topologies. In this paper we demonstrate the synthesis of finned ferrierite (FER), a commercial zeolite with two-dimensional pores, where protrusions on crystal surfaces behave as pseudo nanoparticles. Catalytic tests of 1-butene isomerization reveal a 3-fold enhancement of catalyst lifetime and an increase of 12 % selectivity to isobutene for finned samples compared to corresponding seeds. Electron tomography was used to confirm the identical crystallographic registry of fins and seeds. Time-resolved titration of Brønsted acid sites confirmed the improved mass transport properties of finned ferrierite compared to conventional analogues. These findings highlight the advantages of introducing fins through facile and tunable post-synthesis modification to impart material properties that are otherwise unattainable by conventional synthesis methods.

2D zeolite↗

Development of Self-Assembly Supports Enabling Transformational Membrane Performance for Cost-Effective Carbon Capture

This final technical report describes work conducted by Membrane Technology and Research, Inc. (MTR) for the U.S. Department of Energy (DOE), National Energy Technology Lab (NETL) on the development of membranes with transformational performance for carbon capture under award number DE-FE0031596. The work was performed from June 1, 2018 through May 31, 2024. For more than a decade, MTR has worked in partnership with DOE to develop an innovative membrane-based CO 2 capture process. This effort has included the first test of membrane modules with coal-fired flue gas at the Arizona Public Services (APS) Cholla plant in 2010; the accumulation of >11,000 hours of flue gas operation for Polaris modules on a bench-scale 1 tonne/day (TPD) system at the National Carbon Capture Center (NCCC); scale-up of first-generation (Gen-1) Polaris to a 20 TPD small pilot system, and successful operation of this system on a flue gas slipstream at NCCC and in integrated boiler testing at Babcock & Wilcox (B&W). Through continued development efforts, a second-generation (Gen-2) version of the Polaris membrane has been scaled-up to pilot production. This membrane offers 70% higher CO 2 permeance with similar selectivity to the base case Polaris. MTR also developed planar modules designed specifically for the low-pressure, high-volumetric flow rate process conditions of flue gas operation. These new modules have significantly lower pressure-drop values compared to the type originally used (spiral-wound modules), which results in significant energy savings. The goal of the work described in this report was to improve on the Polaris Gen-2 membrane with the ultimate aim to reduce the cost of carbon capture. The majority of the effort was to develop improved support membranes for the multi-layer composite structure of MTR’s Polaris membrane. Earlier work at MTR had identified the support structure as limiting membrane permeances, not because the support itself represents a permeation resistance, but because the distribution of pores at the surface of the support imposes a geometric restriction to diffusion in the layers above it. Support membranes were prepared from a range of polymers, including commercially available block copolymers and a custom synthesized block copolymer alternative. The best support membranes developed in this project reduced the geometric restriction by a factor of two to three. These supports then were used to produce Polaris composite membranes with improved permeances. The second topic was to create a high-selectivity version of the Polaris membrane. The high-selectivity version uses a novel selective polymeric material and high selectivities were confirmed in experiments at MTR. The material is not easily made into very thin films. Consequently, the permeances are significantly lower than the Polaris Gen-2 membrane. The utility of this membrane is therefore limited to the carbon dioxide purification step that produces liquid CO 2 . A Technical and Economic Analysis (TEA) was performed for a carbon capture system that uses both advanced membrane types. The TEA shows the novel advanced membranes reduce the cost of capture by 10%, from $63.32/tonne CO 2 to $56.90/tonne CO 2 (2022 USD). Most of the development work was carried out with laboratory-scale casting and coating equipment. A number, but not all, of the improvements identified have been implemented on commercial-scale manufacturing equipment. The focus of future work at MTR is to incorporate the advancements made into the Polaris membrane manufacturing process.

01 COAL, LIGNITE, AND PEAT↗

Assessment of oil and gas fields in California as potential CO 2 storage sites

California's total annual greenhouse gas (GHG) emissions (425.3 MtCO 2 e) in 2018 were about 6.4% of the US total (6,677 MtCO 2 e) and around 1% of global emissions. About 39% of 2018 GHG emissions in California were from the industrial and electrical sectors. Many of these emissions were from large stationary point sources and were suitable for carbon capture retrofit with subsequent storage of the captured carbon dioxide (CO 2 ) in geological formations. Previous studies of California found suitable geology and CO 2 storage resource. This study refines and furthers prior work using a three-stage screening process of oil fields, gas fields, and underground natural gas storage (UGS) sites by combining criteria from previous studies while excluding sites that pose technical risk or are located in regions with surface restrictions including sensitive habitats and dense populations. In the first stage, 129 CO 2 storage sites in California were identified using qualification criteria based upon formation properties including geological conditions and pore pressure. The second stage identified sensitive sites by applying conservative screens including seismic activity, faulting, population density, restricted lands, and sensitive habitats. During the third stage, 61 CO 2 potential storage sites were identified by subtraction of stage 2 areas from stage 1. The potential storage volume in the third stage ranged from 1.0 to 2.0 GtCO 2 . Finally, we applied a scoring system with seven parameters to rank the 61 potential sites based on subsurface technical criteria. The scored sites are classified as high priority, medium priority, and sites for future study. Prospective CO 2 storage sites with high and moderate priority were selected and linked to CO 2 sources. There are 14 prospective sites (above 20 MtCO 2 storage resource per site) with a total storage resource of 1024 MtCO 2 distributed in Northern and Southern California. Of these sites, there are 9 potential CO 2 -EOR sites and 1 depleted oil field with a total estimated CO 2 storage volume of ~800 MtCO 2 in the Southern San Joaquin and Ventura Basin. These 10 prospective sites with a storage resource greater than 20 MtCO 2 could potentially deliver more than 20 years of storage with an average injection rate of 40 MtCO 2 /year. The remaining 4 highly prospective sites are in Northern California. Additionally, study results suggest that saline formations should be re-evaluated in concert with storage in oil, gas, and natural gas storage reservoirs.

58 GEOSCIENCES↗

PILBCP-IL Composite Ionomers for High Current Density Performance

Wide-spread commercialization of fuel cell electric vehicles using proton exchange membrane fuel cell (PEMFC) power sources requires that several existing limitations be addressed. These include: (1) a reduction in platinum (Pt) loading in the catalytic electrodes, (2) improvements in reactant and electronic mobility throughout the catalytic electrodes, (3) reduction in the reliance on materials derived from polluting “forever chemicals”, and (4) a significant improvement in the operational longevity of catalytic electrode components. In this project, a team of two universities, Drexel University and Texas A&M University, one national lab, National Renewable Energy Laboratory, and one company, General Motors, collaborated to develop a new cathode ionomer chemistry that would address these limitations and result in an improvement in performance over existing ionomer materials. The key technology developed through this collaborative project was a composite cathode ionomer encompassing an ionic liquid interlayer between Pt catalysts and a sulfonated polymerized ionic liquid block co-polymer (S-PILBCP) that possess the orthogonal properties of protonic conductivity and ionic liquid enhanced kinetics and durability (see schematic in Figure 1). The composite S-PILBCP ionomer eliminates many of the existing issues with perfluorosulfonic acid-based ionomers including active site blocking by sulfonate specific adsorption, restricted O 2 transport through ionomer films, limited humidity tolerance and active area loss for carbon pore confined catalyst particles, and use of polluting “forever chemicals”. Following successful integration of the developed composite ionic liquid into a PEMFC cathode catalyst layer, we demonstrate enhanced performance over Nafion containing cathodes with Pt/C and PtCo/C at both low and high current density. The performance with our composite S-PILBCP ionomer meets the Department of Energy (DOE) targets for light duty vehicle applications

08 HYDROGEN↗

The weathered bedrock vadose zone: A hidden control on water availability in the western United States (Final Technical Report)

This report summarizes the research completed during the project entitled "The weathered bedrock vadose zone: A hidden control on water availability in the western United States" (Award: DE-SC0018039). The report covers work completed between 9/15/17 and 9/14/20. The objective of the exploratory project work was to evaluate the overarching hypothesis that the subsurface underlying actively eroding landscapes is systematically structured as a result of a balance between weathering and erosion and that this structure leads to dynamic moisture storage within weathered rock. To accomplish project objectives, we characterized hillslope weathering profile structure and conducted subsurface moisture monitoring across four study sites with comparable bedrock types but differing erosion rate and hydroclimate: The Angelo Coast Range Reserve (ACRR), Sagehorn-Russell Ranch (SRR), Antelope Valley Ranch (AVR), and the LBNL Watershed Function Scientific Focus Area (SFA). New characterization and monitoring infrastructure at the AVR and SFA sites was established through this project and enables long-term monitoring of hydrologic and carbon fluxes at the sites for comparison with long-term monitoring in other sites. Deep drilling, downhole logging, and characterization of chemistry and porosity of rock recovered from drilling reveals two weathering fronts across the sites. A shallower front is marked by the extent of pervasive fracturing and oxidation of pyrite and organic carbon. A deeper weathering front is marked by the appearance of open fractures and discoloration. At the AVR site, multiple ridge tops of different hillslope length were investigated to establish the relationship between the weathering front and topography. Drilling reveals that the fraction of the hillslope relief that is weathered scales with hillslope length, consistent with the theoretical framework that guided the project relating groundwater drainage with landscape evolution. With the establishment of new boreholes, we have documented significant water storage in weathered bedrock using low-field borehole nuclear magnetic resonance (NMR) across all four sites. In the upper 4 m of bedrock, the average water contents at ACRR, AVR, SRR, and SFA are approximately 25%, 16%, 25%, 23% respectively. At the ACRR site, where long-term monitoring via neutron probe surveys and time-domain transmission sensors are available, we have demonstrated that NMR reliably records moisture content and importantly, moisture content dynamics in weathered bedrock. Across all sites, our NMR monitoring has also revealed that a considerable fraction of the water storage occurs in fractures or pores that are significantly larger than the fine-grained bedrock matrix of the parent rock. To identify the proportion of water storage occurring in fractures and large pores, we analyzed the NMR signal by evaluating both the sum of echoes and an inversion of the NMR signal for a distribution of T2 relaxation times. Both methods reveal that at least 20% of the water storage occurs in fractures and larger pores across all sites. To constrain the interpretation of the in-situ field-scale NMR measurements, we compare the NMR response in the field to independent pore-structure information obtained via helium pycnometry. Results support the inference that dynamic seasonal water storage is dominantly restricted to the fracture network, which can reach 7% of the total volume. This dynamic range in water content is comparable to that of some soils, underscoring the need to understand how this region of the Earth System is structured, and functions as a hydrologic reservoir.

58 GEOSCIENCES↗

Impact of downhole pressure and fluid-access on the effectiveness of wellbore cement expansion additives

Autogenous shrinkage of wellbore cement widely impairs zonal isolation. MgO-based cement expansion additives (CEAs) can mitigate this shrinkage, or even impart net expansion, by creating porosity through displacive crystal growth-processes. However, both MgO hydration and autogenous shrinkage behaviour depend strongly on stress state. Evaluation of CEA performance in wellbore cements therefore requires testing under elevated pressures representative for subsurface environments. We report experiments addressing the chemical and bulk volume changes that occur in cement hydrating at 10 MPa confining pressure and 90 °C. Volumetric response was investigated as function of MgO concentration, external water supply, and pore pressure decrease through water consumption during reaction. Results show the bulk expansion achieved using MgO-based CEAs diminishes markedly with increasing effective confining pressure or, equivalently, upon restricting fluid supply. This reduced expansion-potential under pressure has profound implications for slurry design, notably regarding CEA-concentrations required to counteract micro-annulus formation while maintaining low cement permeability.

36 MATERIALS SCIENCE↗

Electromechanics of lipid-modulated gating of potassium channels

Experimental studies reveal that the anionic lipid phosphatidic acid (POPA), non-phospholipid cholesterol, and cationic lipid DOTAP inhibit the gating of voltage-sensitive potassium (Kv) channels. In this work, we develop a continuum electromechanical model to investigate the interaction of these lipids with the ion channel. Our model suggests that: (i) POPA lipids may restrict the vertical motion of the voltage-sensor domain through direct electrostatic interactions; (ii) cholesterol may oppose the radial motion of the pore domain of the channel by increasing the mechanical rigidity of the membrane; and (iii) DOTAP can reduce the effect of electrostatic forces by regulating the dielectric constant at the channel–lipid interface. The electromechanical model predictions for the three lipid types match well with the experimental observations and provide mechanistic insights into lipid-dependent gating of Kv channels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Insights into scale translation of methane transport in nanopores

Accurate prediction of flow behavior in shale matrix is critical for efficient development of shale gas reservoirs. In these systems, the majority of pores are in the nano-size range. As a result, continuum-based approaches may not be appropriate to simulate flow in such systems. Molecular dynamics (MD) simulations are capable of capturing the relevant microscale physics. Their relatively high computational expense, however, restricts MD simulations to rather small systems and domains. This limitation creates a gap between computational need of macroscale systems and capabilities of MD simulations. The lattice Boltzmann method (LBM) is a suitable candidate to bridge this gap. In this work, the multiple-relaxation-time (MRT)-LBM is used to study methane transport in nano-size pores. Adsorption effects near solid boundaries, as well as non-ideal behavior of fluids, are accounted for via incorporating appropriate force terms in LBM. In this work, parameters associated with the force terms in the equation of state are studied in detail, and a workflow is proposed to determine optimal values of these parameters for gas flow in slit pores. Specifically, we establish these parameters such that the range of density values that the model is able to simulate is maximized. We demonstrate this workflow by simulating gas flow where velocity and density profiles from MD simulations are used as reference data. Results from LBM simulations are in good agreement with MD reference data for pores that are 4 nm in width or larger. Moreover, we propose a preconditioning scheme to improve the stability of LBM in dealing with complex geometries. The robustness of this scheme is demonstrated by simulating several roughness geometries. This work motivates the use of LBM in scale translation of the physics of mass transport in more complex permeable media.

03 NATURAL GAS↗

RockFlow: Fast Generation of Synthetic Source Rock Images Using Generative Flow Models

Image-based evaluation methods are a valuable tool for source rock characterization. The time and resources needed to obtain images has spurred development of machine-learning generative models to create synthetic images of pore structure and rock fabric from limited image data. While generative models have shown success, existing methods for generating 3D volumes from 2D training images are restricted to binary images and grayscale volume generation requires 3D training data. Shale characterization relies on 2D imaging techniques such as scanning electron microscopy (SEM), and grayscale values carry important information about porosity, kerogen content, and mineral composition of the shale. Here, we introduce RockFlow, a method based on generative flow models that creates grayscale volumes from 2D training data. We apply RockFlow to baseline binary micro-CT image volumes and compare performance to a previously proposed model. We also show the extension of our model to 2D grayscale data by generating grayscale image volumes from 2D SEM and dual modality nanoscale shale images. The results show that our method underestimates the porosity and surface area on the binary baseline datasets but is able to generate realistic grayscale image volumes for shales. With improved binary data preprocessing, we believe that our model is capable of generating synthetic porous media volumes for a very broad class of rocks from shale to carbonates to sandstone.

58 GEOSCIENCES↗

Deep Learning Super-Resolution X-Ray Computed Tomography Algorithms for Additive Manufacturing

Industrial X-ray computed tomography (XCT) is a nondestructive method for inspection and characterization of additively manufactured (AM) materials and parts. In practice, the resolution of XCT can be limited by factors such as detector binning, restricted field of view for large-scale objects, system blur, motion during scanning, and acquisition settings. These limitations can reduce the detectability of critical flaws such as pores, cracks, and lack of fusion. Super-resolution (SR) techniques offer a promising solution for improving the effective resolution and image quality of XCT reconstructions without the need for expensive hardware upgrades or laborious, time-consuming scans. In particular, deep learning-based SR methods have garnered attention in recent years as powerful tools for reconstructing high-resolution volumes from low-resolution inputs. In this work, a novel deep learning-based SR method is proposed for XCT scans of AM parts, and compared against several existing state-of-the-art (SOTA) methods. The proposed method, Simurgh-SR, is built on the pre-existing Simurgh framework and consists of a 2.5D U-Net trained to map low-quality inputs containing noise and artifacts to high-quality reconstructions characterized by higher flaw contrast, better noise texture, and reduced artifacts. The experimental results demonstrate superior performance of Simurgh-SR in performing 4× SR on real industrial XCT scans of thick 316L components, enhancing the structural similarity score and peak signal-to-noise ratio (>7dB) compared to the LR counterpart while improving the F1-score for flaw detection by more than 2.3× when compared to alternative SOTA SR methods. This improvement enables more accurate and significantly faster characterization of metal AM components. Additionally, Simurgh-SR was trained for both 2X and 4X SR and performs effectively at both levels, enabling the use of a single model for various SR factors.

Rahman, Obaid [ORNL] (ORCID:0000000277810840)↗

Influence of cementation on the yield surface of rocks numerically determined from digital microstructures

Digital Rock Physics has reached a level of maturity on the characterisation of primary properties that depend on the microstructure – such as porosity, permeability or elastic moduli – by numerically solving field equations on μCT scan images of rock. After small deformations or at depth though, most rocks eventually reach their limit of elasticity and the complementary plastic properties are needed to describe the full mechanical behaviour. Currently, determination of a rock’s yield surface from its microstructure is often restricted to semi-analytical criteria derived by limit analysis or numerical simulations performed on idealised geometries. Such simplification lacks representativeness, particularly for processes that affect directly the pore-grain interface such as the cementation phenomenon, happening during diagenesis. Eventually, only direct numerical simulation of elasto-plasticity performed on digitalised microstructures can be used to assess the strength of different cemented materials and its evolution with the alteration of the microstructure. In this study, we provide a comprehensive parametric study on the impact of cementation on rock strength for real microstructures of cemented granular materials. Compared to most previous studies, the whole yield surface is determined numerically (using Finite Element Method) in order to assess the influence of cementation for different stress-paths. The previously known tendency of rock to strengthen with increasing cementation volume is verified. New results on the influence of cement property namely Young’s modulus, friction and cohesion on the rock’s yield surface are explored. The envelopes obtained are compared to the ones obtained by experimental data and existing models. Finally, the framework presented in this study showcases the wider possibility of determining any rock’s or porous material’s yield surface from its microstructure.

36 MATERIALS SCIENCE↗

Role of heterogeneous surface wettability on dynamic immiscible displacement, capillary pressure, and relative permeability in a CO 2 -water-rock system

Numerical examination of CO 2 displacement in brine filled rock pores with varying wettability. Results show that as the CO 2 contact angles on the rock surface increased (intermediate water wetting conditions) negative entry pressures for CO 2 displacement were observed, indicating that when heterogeneous surface wettability is accounted for, greater mobility of CO 2 into rock pores is expected. In addition, simulated relative permeability CO 2 curves with heterogeneous intermediate water wetting conditions exhibited a greater amount of variation, indicative of a complex interrelation of geometric capillary restrictions and surface wettability.

58 GEOSCIENCES↗

Efficient reconstruction and validation of heterogeneous microstructures for energy applications

The digital reconstruction of microstructures is necessary for simulations in fields ranging from geology to electrochemistry, but the state-of-the-art digital reconstruction techniques often compromise between resolution and field of view. It is challenging to retain detailed microstructure information in large-scale reconstructions. Here, this study investigates different aspects of the Yeong-Torquato algorithm based on correlation functions to make it more efficient. We achieve this goal by reducing the computational complexity of the chord-length distribution function and the two-point correlation function, applying the random sphere-packing method as the initial condition, and restricting potential voxel swaps to interfaces. In addition, a novel superposition parallel scheme is introduced to aid in searching for potential voxel swaps. The algorithm proposed is validated by comparing the pore-size distributions of reconstructed 3D custom battery electrodes from a sample dataset obtained from transmission X-ray microscopy. From a sample image with 200 x 200 pixels, the code can reconstruct a 300 x 300 x 300 structure in under 22 h and reconstruct a 400 x 400 x 400 structure in 43 h with eight cores.

42 ENGINEERING↗