Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Pore Space”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Factors controlling pore network development of thermally mature Early Palaeozoic mudstones from the Baltic Basin (N Poland)

Understanding the formation of pore space, especially in low porosity shales (as source rocks and as unconventional resources), is critical to the oil and gas industry, since pores control the space available for hydrocarbon and participate in hydrocarbon transport. We examined 87 Ordovician and Silurian mudstone samples collected from four wells located in the Pomeranian part of the Baltic Basin (northern Poland), one of the primary Polish targets for hydrocarbon exploration. These samples represent the Pelplin, the Paslek, the Jantar, the Prabuty, and the Sasino Formations, which still requires more detailed porosity studies. Our study aimed to identify factors controlling porosity development, by applying bulk techniques (organic petrology and TOC analyses, quantitative mineralogy, and porosimetry) as well as nano- to microscale techniques (thin section petrography, electron microscopy). The studied samples are mainly argillaceous mudstones. The results of porosimetry measurements, combined with image analysis, indicate that the pores of all studied rocks are dominated by micropores (pores < 2 nm in diameter), mesopores (2-50 nm in diameter) and small macropores. The SEM images showed three main pore types: a) voids related to clay mineral aggregates, b) pores inside organic matter particles, and c) pores between other mineral grains. In the Jantar and Sasino mudstones, the organic matter content and its thermal maturity control porosity. The occurrence of solid bitumen in the rocks from these formations reduces samples’ mesoporosity because of the pore-clogging effect. In contrast, in the Paslek and Prabuty Formations, there is low organic matter content and specific surface area and the volume of mesopores increase with clay minerals content. In the Pelplin mudstones, there are no prevailing factors controlling porosity. As a result, we suggest that a combination of SEM image analysis and dual liquid porosity (DLP) measurements is a powerful method to assess porosity available for petroleum flow in mudstones.

03 NATURAL GAS↗
DOE OSTI2021

Unveiling stimulation fluid-driven alterations in shale pore architecture through combined interpretation of TD-NMR and multi-component gas adsorption

The hydraulic fracturing stimulation process can alter petrophysical and flow properties of the shale matrix in unconventional reservoirs. These alterations can result from interactions with native and exogenous fluids, and can be detrimental to reservoir productivity. Different experimental procedures have been proposed to study reactive processes and thereby understand the aforementioned interactions. Most published works do not consider the rock initial saturation with formation fluid, or the effect of organic matter (OM). To gain insights into these effects, two reactive experimental sets, differing in experimental conditions, were conducted. In the first set, at 125°C and 45 MPa, a stimulation fluid surrogate was put in contact with Baxter Shale samples, previously saturated and equilibrated with synthetic formation fluid. A baseline consisted of exposure to formation fluid exclusively. In the second set, at 0.08 MPa and 25°C, contact with stimulation fluid took place as in the first set. Two additional assays in this second set either skipped the initial formation fluid saturation or removed the OM, respectively. Petrophysical properties analyses, focused on the pore architecture, relied on the combined interpretation of gas adsorption and Time-Domain Nuclear Magnetic Resonance (TD-NMR). Results show that the presence of the formation fluid significantly restricts pore space accessibility to the stimulation fluid and mitigates the pore-architecture alteration to the macro-pore region. Mineral dissolution and precipitation dominate over surface or wettability alteration. Next, results suggest that the OM restricts accessibility to some pores. Ultimately, the equilibration process with synthetic formation fluid drives mainly surface alteration, unhindered by OM.

58 GEOSCIENCES↗
DOE OSTI2024

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↗
DOE OSTI2020

Micro- and nano-porosity of the active Alpine Fault zone, New Zealand

Porosity reduction in rocks from a fault core can cause elevated pore fluid pressures and consequently influence the recurrence time of earthquakes. We investigated the porosity distribution in the New Zealand's Alpine Fault core in samples recovered during the first phase of the Deep Fault Drilling Project (DFDP-1B) by using two-dimensional nanoscale and three-dimensional microscale imaging. Synchrotron X-ray microtomography-derived analyses of open pore spaces show total microscale porosities in the range of 0.1 %–0.24 %. These pores have mainly non-spherical, elongated, flat shapes and show subtle bipolar orientation. Scanning and transmission electron microscopy reveal the samples' microstructural organization, where nanoscale pores ornament grain boundaries of the gouge material, especially clay minerals. Our data imply that (i) the porosity of the fault core is very small and not connected; (ii) the distribution of clay minerals controls the shape and orientation of the associated pores; (iii) porosity was reduced due to pressure solution processes; and (iv) mineral precipitation in fluid-filled pores can affect the mechanical behavior of the Alpine Fault by decreasing the already critically low total porosity of the fault core, causing elevated pore fluid pressures and/or introducing weak mineral phases, and thus lowering the overall fault frictional strength. We conclude that the current state of very low porosity in the Alpine Fault core is likely to play a key role in the initiation of the next fault rupture.

58 GEOSCIENCES↗
DOE OSTI2021

Effects of Pore Connectivity on the Sorption of Fluids in Nanoporous Material: Ethane and CO2 Sorption in Silicalite

Adsorption of fluids in nanoporous materials is important for several applications including gas storage and catalysis. The pore network in natural, as well as engineered, materials can exhibit different degrees of connectivity between pores. While this might have important implications for the sorption of fluids, the effects of pore connectivity are seldom addressed in the studies of fluid sorption. We have carried out Monte Carlo simulations of the sorption of ethane and CO2 in silicalite, a nanoporous material characterized by sub-nanometer pores of different geometries (straight and zigzag channel like pores), with varied degrees of pore connectivity. The variation in pore connectivity is achieved by selectively blocking some pores by loading them with methane molecules that are treated as a part of the rigid nanoporous matrix in the simulations. Normalized to the pore space available for adsorption, the magnitude of sorption increases with a decrease in pore connectivity. The increased adsorption in the systems where pore connections are removed by blocking them is because of additional, albeit weaker, adsorption sites provided by the blocker molecules. By selectively blocking all straight or zigzag channels, we find differences in the absorption behavior of guest molecules in these channels.

Gautam, Siddharth (ORCID:0000000314435382)↗
DOE OSTI2023

Quantifying Pore Morphology in Spray-Formed Tantalum Using X-ray Micro-computed Tomography

In order to establish quantitative process–structure–property relationships in thermal spray coatings, a robust framework for defining (micro)structural characteristics is needed. Here, we present a quantitative characterization of the three-dimensional morphology of porosity in spray-formed tantalum samples based on high-resolution X-ray micro-computed tomography. Using synchrotron facilities, we acquired dozens of high-resolution scans, enabling a statistically meaningful comparison across multiple samples, different regions within samples, and spray processes. We quantify the spatial distribution, size, and topology of porous inclusions, with a significant focus on variability across samples and different spray processes (plasma and cold sprayed), as well as sensitivity to image segmentation and resolution. Based on a typical segmentation, we report porosities ranging from 0.9 to 1.7 pct for all samples tested, with significant sensitivity due to image segmentation resulting in estimates as low as 0.6 pct and as high as 4.8 pct. For the complex pore space morphology observed in these materials, we argue that a conventional analysis based on identifying individual pores is not well suited, and propose an alternative approach based on morphological metrics with a rich history in porous media literature, such as spatial correlations, local pore thickness, and scale-dependent sub-sampling. Spatial correlations indicate anisotropic splat structures, but only mildly anisotropic pores. Various measures of pore size show a wide distribution of sizes, ranging from sub-micron to 10-micron length scales. Scale-dependent variations in porosity suggest that representative volumes of several hundred microns are required for convergence of morphological metrics, with larger volumes for cold-spray materials. This work provides a robust quantitative basis for describing three-dimensional pore structure in thermal spray coatings.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A wide-range thermodynamically consistent constitutive model for dynamic loading of partially saturated porous media

SUMMARY This paper presents a novel thermodynamically consistent constitutive model for partially saturated porous rocks across a wide range of conditions. The material states generated behind the shock wave from an explosive source can vary significantly, ranging from crushed and melted rock near the source to a poroelastic medium in the far field. In the model, rock strength is determined by the effective pressure, which is calculated using two independent equations of state: one for the solid rock and another for the pore fluid. The model accounts for shock-induced liquefaction resulting from fluid pressure buildup in the pore spaces near the explosive source. Simultaneously, it describes the increase in wave propagation speed due to elastic pore contraction in both dry and partially saturated rocks. This model is applied to investigate how fluid saturation affects the amplitude and shape of the generated waves, as well as the residual stress surrounding the cavity formed by spherical explosions.

Vorobiev, Oleg [Lawrence Livermore National Labora↗
DOE OSTI2024

Approaches for the Simulation of Coupled Processes in Evolving Fractured Porous Media Enabled by Exascale Computing

Models have historically represented fractured porous media with continuum descriptions that characterize the media using bulk parameters. The impact of small-scale features is not captured in these models, although they may be controlling the performance of subsurface applications. Pore-scale models can simulate processes in small-scale features by representing the pore space geometry explicitly but are computationally expensive for large domains. The alternative multiscale approach entails the combination of pore-scale and continuum-scale descriptions in a single framework. We use Chombo-Crunch, a computational capability that discretizes complex geometries with an adaptive, embedded boundary method to contrast these two approaches. Chombo-Crunch takes advantage of recent computational performance and memory bandwidth improvements resulting from the emergence of exascale computing resources. These combined improvements enable the efficient simulation of reactive transport in fractured media with a high degree of fidelity and the ability to capture the control small-scale processes exert on the overall medium evolution.

42 ENGINEERING↗
DOE OSTI2022

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↗
DOE OSTI2021

Dissolution of periodic arrays of grains: Upscaling of pore-scale simulations with fast reactions

Upscaling methods are frequently used to derive transport equations at the macroscopic scale from more fundamental equations formulated at the pore scale. These methods typically give a suitable structure for the macroscopic equations and can also provide explicit expressions for the constitutive parameters, such as permeability and dispersion coefficients. Introducing chemical reactions complicates upscaling in at least two important ways. First, the interplay between chemical reactions and transport processes introduces a new length scale, which can be much smaller than the convective or dispersive length scales. A small reactive length scale breaks one of the key assumptions in upscaling; that there is a significant separation in length between the pore-scale and macro-scale processes. The second complication is that if reactions take place at mineral surfaces (dissolution or precipitation) then the pore space itself is evolving in time. In this paper we suggest ways in which these difficulties can be approached, based on analysis of pore-scale simulation data. First, we noticed that the concentration field in successive unit cells has an almost identical spatial variation, with a single scaling factor for each unit cell that is proportional to the incoming reactant flux. Using pore-scale simulations to determine the mass transfer coefficient in a few unit cells, we can calculate the concentration field in the whole domain, even when dissolution is entirely transport limited. Second, we have noticed a time-dependent mapping of the grain shapes from different unit cells. From these observations, we can deduce constitutive relations where the only time-varying parameter is the porosity. Here, we show that a model based on these ideas can quantitatively account for the pore-scale simulation data.

58 GEOSCIENCES↗
DOE OSTI2025

Interactions of Polar and Nonpolar Groups of Alcohols in Zeolite Pores

Understanding the quantitative interactions among zeolite pore walls, Bro̷nsted acid sites, and molecules with both polar and nonpolar regions is essential for scoping out the potential of zeolites as sorbents and catalysts. Purely siliceous zeolites (MFI and Beta in the present study) are hydrophobic, whereas those containing aluminum are considered hydrophilic, preferentially adsorbing organic molecules even in aqueous environments. To characterize these interactions, we use primary alcohols of increasing molecular weight, quantifying their specific interactions in the confined pore space of the alkyl (CH x ) and OH groups. Three types of interactions were identified: (i) alkyl CH x groups interacting with the zeolite pore walls (approximately 10 kJ mol −1 per carbon), (ii) alcohol OH groups interacting with the pore walls (30−35 kJ mol −1 ), and (iii) alcohol OH groups interacting with Bro̷nsted acid sites (37 kJ mol −1 ). All three interactions were well mirrored by computational simulations. The contribution of the alkyl CH x groups was inferred from the incremental increase in sorption enthalpy with increasing molecular weight; the interaction strength of the OH groups was determined by extrapolating the global adsorption enthalpy of the alcohols to a hypothetical OH group without an alkyl group. This value was identical to the adsorption enthalpy of water. The experiments demonstrated that only water has an adsorption enthalpy on zeolite pore walls lower than its condensation enthalpy (30−35 kJ mol −1 vs 45 kJ mol −1 ), limiting the concentration of water that can be adsorbed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A Monte Carlo Approach to Approximating the Effects of Pore Geometry on the Phase Behavior of Soil Freezing

Abstract Freezing in porous media is associated with a host of dynamic phenomena that stem from the presence and mobility of premelted liquid at subzero temperatures. Accurate assessments of the progressive liquid‐ice phase transition is required for predictive models of frost damage, glacier‐till coupling, and many other cold regions processes, as well as for evaluating the capacity for water storage in near‐surface extraterrestrial environments. We use a Monte Carlo approach to sample the pore space in a synthetic 3D packing of poly‐dispersed spherical particles and evaluate local geometrical constraints that allow us to assess changes in the relative proportions of pore fluid and ice. By approximating the phase boundary geometry in fine‐grained pores while considering both the curvature of the liquid‐ice interface and wetting interactions with matrix particles, our model predicts changes in phase equilibrium in granular media over a broad temperature range, where present accounting for the colligative effects of chloride and perchlorate solutes. In addition to formulating the constitutive behavior needed to better understand properties and processes in frozen soils, our results also provide insight into other aspects of phase equilibria in porous media, including the formation of methane hydrates in permafrost and marine sediments, and the partitioning between liquid water and vapor in the vadose zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Pore‐Scale Modeling of Reactive Transport with Coupled Mineral Dissolution and Precipitation

Abstract We present a new pore‐scale model for multicomponent advective‐diffusive transport with coupled mineral dissolution and precipitation. Both dissolution and precipitation are captured simultaneously by introducing a phase transformation vector field representing the direction and magnitude of the overall phase change. An effective viscosity model is adopted in simulating fluid flow during mineral dissolution‐precipitation that can accurately capture the velocity field without introducing any empirical parameters. The proposed approach is validated against analytical solutions and interface tracking simulations in simplified structures. After validation, the proposed approach is employed in modeling realistic rocks where mineral dissolution and precipitation are dominant at different locations. We have identified three regimes for mineral dissolution‐precipitation coupling: (a) compact dissolution‐precipitation where dissolution is dominant near the inlet and precipitation is dominant near the outlet, (b) wormhole dissolution with clustered precipitation where dissolution generates wormholes in the main flow paths and precipitation clogs the secondary flow paths, and (c) dissolution dominant where all solid grains are gradually dissolved. In the three regimes, the proposed approach provides reliable porosity‐permeability relationships that cannot be described well by traditional macroscale models. We find that the permeability can increase while the overall porosity decreases when the main flow paths are expanded by dissolution and adjacent pore spaces are clogged by precipitation.

58 GEOSCIENCES↗
DOE OSTI2023

Pore-scale simulation of multiphase flow and reactive transport processes involved in geologic carbon sequestration

Multiphase flow and reactive transport are two essential physicochemical processes that govern the effectiveness of geological carbon sequestration (GCS). The interaction and feedback among different phases and components during intricate physicochemical processes hold great significance in understanding CO 2 sequestration. Pore-scale simulations can account for multiphase flow and reactive transport processes in porous media and obtain spatial distributions of parameters (density, velocity, concentration, etc.) in the pore space as well as their temporal evolutions. This proves especially valuable considering that experiments can be hindered by constraints in spatial and temporal resolution. The comprehensive insights garnered from pore-scale research can be leveraged for continuum modeling using the representative elementary volume (REV) concept. In this contribution, four sequential mechanisms of CO 2 -brine-rock interaction in three zones delineated by CO 2 saturation are elaborated to elucidate complicated physicochemical processes involved in GCS, which are followed by general descriptions of mathematical equations and pore-scale numerical methods. In addition, as interested and commonly encountered processes, leakage risks during GCS and CO 2 -enhanced oil recovery (CO 2 -EOR) processes are presented. The existing challenges and future directions are discussed for both the performance of the pore-scale models and the current gaps in the field of GCS. Importantly, we expect that this review will prove beneficial for researchers interested in pore-scale simulations, GCS, and related disciplines.

58 GEOSCIENCES↗
DOE OSTI2024

A mineral precipitation model based on the volume of fluid method

A novel volume of fluid method is presented for mineral precipitation coupled with fluid flow and reactive transport. The approach describes the fluid-solid interface as a smooth transitional region, which is designed to provide the same precipitation rate and viscous drag force as a sharp interface. Specifically, the governing equation of mineral precipitation is discretized by an upwind scheme, and a rigorous effective viscosity model is derived around the interface. The model is validated against analytical solutions for mineral precipitation in channel and ring-shaped structures. It also compares well with interface tracking simulations of advection-diffusion-reaction problems. Here, the methodology is finally employed to model mineral precipitation in fracture networks, which is challenging due to the low porosity and complex geometry. Compared to other approaches, the proposed model has a concise algorithm and contains no free parameters. In the modeling, only the pore space requires meshing, which improves the computational effciency especially for low-porosity media.

58 GEOSCIENCES↗
DOE OSTI2020

Understanding Stabilization in Nanoporous Intermetallic Alloy Anodes for Li-Ion Batteries Using Operando Transmission X-ray Microscopy

Tin-based alloying anodes are exciting due to their high energy density. Unfortunately, these materials pulverize after repetitive cycling due to the large volume expansion during lithiation and delithiation; both nanostructuring and intermetallic formation can help alleviate this structural damage. Here, these ideas are combined in nanoporous antimony–tin (NP-SbSn) powders, synthesized by a simple and scalable selective-etching method. The NP-SbSn exhibits bimodal porosity that facilitates electrolyte diffusion; those void spaces, combined with the presence of two metals that alloy with lithium at different potentials, further provide a buffer against volume change. This stabilizes the structure to give NP-SbSn good cycle life (595 mAh/g after 100 cycles with 93% capacity retention). Operando transmission X-ray microscopy (TXM) showed that during cycling NP-SbSn expands by only 60% in area and then contracts back nearly to its original size with no physical disintegration. The pores shrink during lithiation as the pore walls expand into the pore space and then relax back to their initial size during delithiation with almost no degradation. Importantly, the pores remained open even in the fully lithiated state, and structures are in good physical condition after the 36th cycle. The results of this work should thus be useful for designing nanoscale structures in alloying anodes.

25 ENERGY STORAGE↗
DOE OSTI2021

Pore accessibility and trapping of methane in Marcellus shale

Accessibility of pore space in unconventional reservoirs is an important factor influencing both methane storage capacity and the kinetics of methane desorption. The determination of open (accessible) versus closed (inaccessible) porosity is therefore vital for the prediction of gas production potential. This study investigates accessibility of pores to methane in overmature middle Devonian Marcellus shale samples (cut parallel and perpendicular to bedding) using small and ultra-small neutron scattering (SANS and USANS) with contrast matching (CM), supplemented by other complementary techniques, such as mercury injection capillary pressure (MICP) and low pressure gas (N2 and CO2) adsorption. Our results demonstrate that for the samples studied, only about 6 % of pores with diameter 25-500 nm are accessible to methane. The accessibility fraction for pores larger than 500 nm is 35%. For nanopores smaller than 25 nm, pore accessibility could not be quantitatively determined due to increased methane density and condensation effects in confinement. Our observations indicate that methane penetrates the accessible small mesopores and micropores down to at least 1 nm in diameter, and the density of confined deuterated methane (CD4) is 0.68 g/cm3 for pores of diameter 25 nm and gradually increases with the decreased pore size. Moreover, elevated gas pressure causes formation of additional high-density methane nano-clusters. These clusters have a form of slightly anisotropic polydisperse discs oriented along bedding plane, about 1-12 nm in diameter and with average thickness of 3.6 nm. Utilizing samples cut parallel and perpendicular to the bedding, this study also briefly addresses anisotropy of pores. Based on the iso-size intensity ratio R, our SANS and USANS results demonstrate anisotropy in the out-of-bedding direction and suggest that the degree of anisotropy depends on the pore size. Specifically, for pore diameters ~2.5 to 250 nm, the extent of anisotropy is smaller than for pores ~500 nm to 6 µm in diameter. Finally, comparison of pore size distribution results calculated from SANS/USANS to those obtained using MICP shows good agreement at low pressures, but large difference at pressures above 1000 bar. Here, this discrepancy requires further testing; it is possible that the high mercury pressure used in MICP alters the mesopore and micropore structure of shales.

03 NATURAL GAS↗
DOE OSTI2025

Organic molecules are deterministically assembled in variably inundated river sediments, but drivers remain unclear

Dissolved organic matter (DOM) is central to ecosystem function. A fundamental challenge is understanding the processes leading to variation in the chemistry of organic molecules that comprise DOM. Here we study these processes in variably inundated riverbed sediments, as an understudied, yet ubiquitous component of rivers. Using null-model approaches adopted from community ecology, we found that within-site variation in environmental conditions caused non-random (i.e., deterministic) shifts in DOM chemistry. Deterministic shifts were observed across diverse biomes, though the strength of determinism varied substantially. We found that the strength of determinism decreased with increasing sediment moisture, but in the form of a constraint space. Many systems fell below the upper constraint boundary, however. We propose a conceptual model based on our results and other publications in which DOM assemblages are hypothesized to be increasingly deterministic across the continuum from the river water column to saturated sediment pore spaces to unsaturated and dry soils/sediments.

59 BASIC BIOLOGICAL SCIENCES↗