Search NASA⌕ Search

SEARCH · Search NASA

Results for “Pore Network Model”

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

Striving to translate shale physics across ten orders of magnitude: What have we learned?

Shales will play an important role in the successful transition of energy from fossil-based resources to renewables in the coming decades. Aside from being a significant source of low-carbon intensity fuels, like natural gas, they also serve as geologic seals of subsurface formations that may be used to isolate nuclear waste, sequester CO 2 , or store intermittent energy (e.g., solar hydrogen). Despite their importance, shales pose significant engineering and environmental challenges due to their nanoporous structure and extreme heterogeneity that spans at least ~10 orders of magnitude in spatial scale. Two challenges inhibit a system-level understanding: (1) the physics of fluid flow and phase behavior in shales are poorly understood due to the dominant molecular interactions between minerals and fluids under confinement, and (2) the apparent lack of scale separation that prevents a reliable (closed) description of the physics at any single scale of observation. In this review, we focus on the latter issue and discuss scale translation, which in its broadest sense is transforming data or simulations from one spatiotemporal scale to another. While effective scale translation is not exclusive to shales, but all geologic porous media, the need for it is especially acute in shales given their high degree of heterogeneity. Classical theories like homogenization, while indispensable, fail when scales are not separated. Other methods, like numerical upscaling, scale-translate in only one direction: small to large, but not the reverse, called downscaling. However, the confluence of advances in three areas are bringing challenging problems such as shales within reach: increased computational power and scalable algorithms; high-resolution imaging and multi-modal data acquisition; and machine learning to process massive amounts of data. While these advances equip geoscientists with a wide array of experimental and computational tools, no individual tool can probe the entire gamut of heterogeneity in shales. Their effective use, therefore, requires an ability to bridge between various data types obtained at different scales. The aim of this review is to present a coherent account of computational and experimental methods that may be used to achieve just that, i.e., to perform scale translation. We provide a broader definition of scale translation, one that transcends classical homogenization and upscaling methods, but is consistent with them and accommodates notions like downscaling and data translation. After a brief introduction to homogenization, we review hybrid methods, numerical upscaling and its recent extensions, multiscale computing, high-resolution imaging, and machine learning. We place particular emphasis on multiscale computing and propose an algorithmic framework to bridge between the pore (micro) and Darcy (macro) scales. Throughout the paper, we draw comparisons between the various methods and highlight their (often hidden) similarities, differences, benefits, and pitfalls. We finally conclude with two case studies on shales that exemplify some of the methods presented.

58 GEOSCIENCES↗

Secondary Pore Network Design Optimization Analytical Diffusion Model for Lithium Ion Battery

Fast ion transport along the electrode thickness is required for Lithium ion battery (LIB) to sustain continuous high-rate discharge or charge currents, which is a desired feature for the fast charging of electric vehicle batteries. Advance in LIB electrode manufacturing allows to control, to a certain extent, the electrode microstructure architecture in order to improve the macroscopic performances of these devices. For instance, several groups have experimentally tested architectures tailored specifically to increase the through plane ionic diffusion in regards for fast charge application by adding channels or groove-like large pores oriented along the electrode thickness. These voids do not replace the existing pores of the electrode microstructure but consist in a second, coarser, pore network, called the secondary pore network (SPN). The model quantifies the impact of SPN on ionic diffusion with a composite electrode through a simple analytical approach, considering various pore channel geometries and comparing to standard electrodes with identical gravimetric and identical volumetric specific theoretical capacity. The goal is to find optimal SPN design parameters according to three optimization objectives aiming to balance the improved overall through-plane, thanks to the coarse aligned channels, and the degraded in-plane diffusion, due to the porous matrix densification required to maintain gravimetric and volumetric specific theoretical capacity. The analytical model must be considered as a pre-screening tool to quickly investigate large design space, that can be refined later with more complex electrochemical CPU-expensive numerical models.

Usseglio Viretta, Francois↗

Enabling fast charging of lithium-ion batteries through secondary- /dual- pore network: Part I - Analytical diffusion model

Battery performance is strongly correlated with electrode microstructural properties. Enabling fast charging of lithium-ion batteries requires improved through-plane ionic diffusion that can be achieved through, among other strategies, structured electrodes with a secondary- or dual-pore network (SPN). In this work, an analytical model investigates the impact of such an SPN on ionic diffusion with a composite electrode, considering various pore-channel geometries and comparing to standard electrodes with identical gravimetric- and volumetric-specific theoretical capacities. Relevant SPN design parameters and tortuosity coefficients are identified according to three optimization objectives that aim to balance the improved overall through-plane diffusion, thanks to the coarse aligned channels, and degraded in-plane diffusion because of the porous matrix densification required to maintain gravimetric- and volumetric-specific theoretical capacities. The model indicates that a relatively low amount of SPN is required and that electrodes with high through-plane tortuosity and low in-plane tortuosity benefit most from such architecture.

25 ENERGY STORAGE↗

Experimental Characterization of Hydrogen Diffusion in Shale Rocks for Geologic Storage Applications

As global energy systems undergo a transition to cleaner alternatives, geologic hydrogen storage has emerged as a promising solution for large-scale energy storage. A critical factor in determining the feasibility of this approach is the effectiveness of caprock formations, such as shale, in preventing hydrogen migration. This study investigates the diffusion behavior of hydrogen through shale to assess its suitability as a caprock for geologic hydrogen storage. Using a novel double-seal core holder design and a through-diffusion apparatus, hydrogen diffusion was measured through shale rock from the Eagle Ford and Wolfcamp Formations under dry conditions. These measurements were complemented by microstructural and mineralogical analyses using low-pressure nitrogen adsorption and X-ray diffraction. The effective diffusion coefficient of hydrogen in these shale caprocks ranged from 2.51 × 10 –8 to 9.85 × 10 –8 m 2 /s. Notably, we observed that the diffusion behavior was more related to the pore network structure and could not be attributed to differences in the total pore volume between shale types alone. Here, to further understand the role of pore network complexity, a fractal pore model was developed to correlate tortuosity with the fractal dimension of the pore structure (a measure of pore network complexity). The proposed model closely matched tortuosity values obtained from diffusion experiments, outperforming existing theoretical tortuosity–porosity correlations. These findings provide key quantitative parameters needed to assess the feasibility of geologic hydrogen storage as well as insights that can be applied to hydrogen storage in a range of geologic formations.

08 HYDROGEN↗

Soil porous microstructure control over soil organic matter mobility: A multimethod workflow for understanding chemistry-dependent organic matter binding in soil

Soil organic matter (SOM) has attracted a great deal of interest; particularly for its potential to mitigate human derived CO 2 emissions. Studies have demonstrated that SOM plays a critical role in carbon storage and CO 2 sequestration. However, the sorption properties of SOM, which influence its transport in pore water and stabilization within the soil, remain poorly understood. This study develops a workflow to: (1) examine compound-specific advective and diffusive transport and desorption behaviors, (2) quantify desorption rates through stop-flow and continuous-flow column experiments, and (3) evaluate the impact of soil microporosity on SOM mobility using high-resolution imaging and extractions. Intact core column experiments were conducted on Uncultivated (Natural) and Cultivated soil samples, both were arid soils, collected in Washington State. X-ray computed tomography was employed to measure porosity and pore connectivity, while Fourier-transform ion cyclotron resonance mass spectrometry was used to analyze SOM composition. The findings revealed that cultivation increased total carbon and nitrogen levels due to irrigation and fertilization, enhancing carbon capture potential in arid soils. In contrast, the Natural soil, characterized by higher porosity and connectivity, contained more oxidized carbon. Pore network analysis indicated that soil compaction in the Cultivated soil may lead to longer diffusion pathways, significantly influencing SOM transport and stability.

Hydraulic Properties↗

Pore-Scale Simulation of Spectral Induced Polarization

A new model was developed to simulate spectral induced polarization (SIP) in porous media. The model is based on a pore-network approximation of the soil or rock, in which bonds between pores are represented by equivalent circuits. In the circuit model, electrolytic and surface conduction occur in parallel, with the electrolyte behaving as a single resistor, and the mineral surface behaving as a resistor and capacitor in series. Based on the pore geometry, fluid conductivity, and mineral surface conductivity, an impedance network is generated. SIP spectra are calculated by solving the sinusoidal steady-state electrical problem for a range of frequencies. The new model provides a tool to bridge the gap between (1) our fundamental understanding of polarization processes at the scale of a single pore or grains and (2) existing simulation capabilities for macroscale systems that do not account for pore-scale properties. The model is implemented in Matlab and can solve two-dimensional and three-dimensional (2D and 3D) networks. Examples for regular cubic lattices are presented, but the method and code are flexible with respect to pore geometries and allow for future leveraging of existing open-source capabilities for generation of pore networks and extraction of network geometries from micro imaging of core samples.

58 GEOSCIENCES↗

Flow Reduction in Pore Networks of Packed Silica Nanoparticles: Insights from Mesoscopic Fluid Models

A modified many-body dissipative particle dynamics (mDPD) model is rigorously calibrated to achieve realistic fluid–fluid/solid interphase properties and applied for mesoscale flow simulations to elucidate the transport mechanisms of heptane liquid and water, respectively, through pore networks formed by packed silica nanoparticles with a uniform diameter of 30 nm. Two million CPU core hours were used to complete the simulation studies. Results show reduction of permeability by 54–64% in heptane flow and by 88–91% in water flow, respectively, compared to the Kozeny–Carman equation. In these nanopores, a large portion of the fluids are in the near-wall regions and thus not mobile due to the confinement effect, resulting in reduced hydraulic conductivity. Moreover, intense oscillations in the calculated flow velocities also indicate the confinement effect that contests the external driven force to flow. Here, the generic form of Darcy’s law is considered valid for flow through homogeneous nanopore networks, while permeability depends collectively on pore size and surface wettability. This fluid-permeability dependency is unique to flow in nanopores. In addition, potential dependence of permeability on pore connectivity is observed when the porosity remains the same in different core specimens.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessing pore network heterogeneity across multiple scales to inform CO2 injection models

Geologic heterogeneity is a key feature that must be considered when translations of scaled data are performed. This paper presents the assessment of geologic heterogeneity using a multiscale workflow that includes image analysis-based methods coupled with well log analysis to provide data in which fractals and machine learning methods estimate the carbon dioxide (CO 2 ) storage resource potential of a reservoir. The heterogeneity of rock properties of the complex Bell Creek reservoir in Montana, USA, was explored at the pore scale (~nm to mm), core scale (~mm to m), and well scale (~cm to m). The data used in this study included advanced image analysis of micro-CT (computed tomography) images (pore scale), thin sections (pore scale), plugs and core images (core scale) and well logs (well scale). The micro-CT images were segmented using a U-net segmentation approach into objects of pores and grains. Further, the segmented images were reconstructed into subvolumes of different sizes. Physical properties (porosity and permeability) and fractal dimensions were calculated for the various subvolumes, and Lorenz coefficient (Lc) values, a single parameter to describe the degree of heterogeneity within a pay zone section, were calculated from thin-section images and well logs. Porosity and fractal dimension values were used to estimate the 188-µm threshold of representative elementary volume (REV) in this study. Both the Lc and fractal dimension values were found to be negatively correlated. When these two parameters are combined, it is possible to discern differences in the complex porous networks of the samples analyzed in this study.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Solvent Transport in Disordered and Dynamic Membrane Pores: Implications for Reverse Osmosis and Nanofiltration Membranes

Pressure-driven separations with nanoporous membranes, such as reverse osmosis and nanofiltration, play a vital role in addressing water scarcity and enabling resource recovery. Understanding water or solvent transport in membrane pores is essential for advancing membrane separation technologies. A key question in transport modeling is to establish a relationship between solvent permeability and membrane porous structure properties, such as porosity or pore size. The nano- and subnanometer pores in polymeric membranes such as reverse osmosis and nanofiltration membranes are highly tortuous and dynamically connected, which challenges the conventional methods of transport modeling. This study addresses this challenge by developing a theoretical framework to describe solvent transport through membranes with dynamic and disordered porous structure. Specifically, we propose a lattice model to describe the pore network, while preserving the viscous nature of solvent permeation. We further establish a relationship between solvent permeability and membrane porosity or pore size, which is validated by molecular dynamics simulations and experimental data. By integrating this relationship into the solution-friction model, we define pore connectivity and local friction coefficient to quantify the impact of pore structure on solvent permeability. Our analysis highlights the dominant influence of pore connectivity on the permeability of reverse osmosis and nanofiltration membranes, particularly when the pore size approaches the dimensions of solvent molecules. Overall, this study provides critical insights into water and solvent transport mechanisms in nanoporous membranes, opening the door for strategies to substantially enhance membrane performance.

membranes↗

ThermoPore: Predicting part porosity based on thermal images using deep learning

Part qualification is often a critical and labor-intensive process in additive manufacturing, particularly in the detection of defects such as porosity, which stands to benefit significantly from advancements in machine learning. We present a deep learning approach for quantifying and localizing ex-situ porosity within Laser Powder Bed Fusion fabricated samples utilizing in-situ thermal image monitoring data. Our goal is to build the real time porosity map of parts based on thermal images acquired during the build. The quantification task builds upon the established Convolutional Neural Network model architecture to predict pore count and the localization task leverages the spatial and temporal attention mechanisms of the novel Video Vision Transformer model to indicate areas of expected porosity. Our model for porosity quantification achieved a R 2 score of 0.57 and our model for porosity localization produced an average Intersection over Union (IoU) score of 0.32 and a maximum of 1.0. This work is setting the foundations of part porosity “Digital Twins” based on additive manufacturing monitoring data and can be applied downstream to reduce time-intensive post-inspection and testing activities during part qualification and certification. In addition, we seek to accelerate the acquisition of crucial insights normally only available through ex-situ part evaluation by means of machine learning analysis of in-situ process monitoring data.

Deep learning↗

A Comparative Study of Deep Learning Models for Fracture and Pore Space Segmentation in Synthetic Fractured Digital Rocks

This study focuses on the comparative study of deep learning (DL) models for pore space and discrete fracture networks (DFNs) segmentation in synthetic fractured digital rocks, specifically targeting low-permeability rock formations, such as shale and tight sandstones. Accurate characterization of pore space and DFNs is critical for subsequent property analysis and fluid flow modeling. Four DL models, SegNet, U-Net, U-Net-wide, and nested U-Net (i.e., U-Net++), were trained, validated, and tested using synthetic datasets, including input and label image pairs with varying properties. The model performance was assessed regarding pixel-wise metrics, including the F1 score and pixel-wise difference maps. In addition, the physics-based metrics were considered for further analysis, including sample porosity and absolute permeability. Particularly, We first simulated the permeability of porous media containing only pore space and then simulated the permeability of porous media with DFNs added. The difference between these two values is used to quantify the connectivity of segmented DFNs, which is an important parameter for low-permeability rocks. The pixel-wise metrics showed that the nested U-Net model outperformed the rest of the DL models in pore space and DFNs segmentation, with the SegNet model exhibiting the second-best performance. Particularly, nested U-Net enhanced segmentation accuracy for challenging boundary pixels affected by partial volume effects. The U-Net-wide model achieved improved accuracy compared to the U-Net model, which indicated the influence of parameter numbers. Similarly, nested U-Net has the closest match to the ground truth of physics-based metrics, including the porosity of pore space and DFNs, and the permeability difference quantifying the connectivity of DFNs. The findings highlight the effectiveness of DL models, especially the U-Net++ model with nested architecture and redesigned skip connections, in accurately segmenting pore spaces and DFNs, which are crucial for pore-scale fluid flow and transport simulation in low-permeability rocks.

Wang, Hongsheng↗

Cavitation in a soft porous material

Abstract We study the collapse and expansion of a cavitation bubble in a deformable porous medium. We develop a continuum-scale model that couples compressible fluid flow in the pore network with the elastic response of a solid skeleton. Under the assumption of spherical symmetry, our model can be reduced to an ordinary differential equation that extends the Rayleigh–Plesset equation to bubbles in soft porous media. The extended Rayleigh–Plesset equation reveals that finite-size effects lead to the breakdown of the universal scaling relation between bubble radius and time that holds in the infinite-size limit. Our data indicate that the deformability of the porous medium slows down the collapse and expansion processes, a result with important consequences for wide-ranging phenomena, from drug delivery to spore dispersion.

Leng, Yu↗

Nanopore networks in colloidal silica assemblies characterized by XCT for confined fluid flow modeling

Fluid flow through nanopore structures has exhibited different behavior than that described by Darcy's law, derived from pore networks of micrometer scale. Here, pressure assemblies of colloidal silica spheres of specified diameters from 94 to 620 nanometers were prepared as model nano-porous media. The smaller silica spheres tended to form random packing, but some substantial ordering was found in the packing of the larger silica spheres. X-ray computed tomography (XCT), with a voxel resolution of 16 nm, was used to characterize the nanopore networks. The porosity of the Nano-XCT pore network (34%) was less than the total porosity of the assemblies as determined by X-ray attenuation (52%), possibly due to the limits of voxel resolution during segmentation. Fluid flow in the nanopore networks was simulated using the single-phase Lattice Boltzmann Method, and the simulated permeability was compared with the empirical and experimental values. Based on our characterization, although a well-ordered packing of silica spheres was not achieved, it was found that the nanopore networks in the colloidal silica assemblies had pore size distributions corresponding to the particle sizes. The simulated permeability was less than the experimental measurement for water flow, but the complex packing of silica spheres and surface chemistry issues need to be considered in future research.

02 PETROLEUM↗

Red Wine Fermentation Alters Grape Seed Morphology and Internal Porosity

During wine fermentation, grape berry components, including seeds, undergo extensive physical and chemical changes that result in the release of flavonoids, such as tannins, from seeds into wine. Understanding changes in seed morphology during fermentation is crucial for aiding the development of accurate prediction models for flavonoid extraction during winemaking, which enhances fermentation management and ensures consistency in wines from year to year. High-resolution x-ray microcomputed tomography (x-ray μCT) was used to investigate the effect of red wine fermentation on changes in grape seed morphology. Using a PyTorch-based implementation of a fully convolutional network with a Resnet-101 backbone for semantic segmentation of x-ray μCT images, we quantified extensive alteration to grape seed structure during fermentation. Image analyses revealed the development of a pore network breaking apart the seed endosperm by the end of fermentation, leading to an increase in surface area. Fermentation significantly altered grape seed morphology. Such alterations could enable transport of seed flavonoids from inside the endosperm and integument to outside the seed. Further research on the physical processes occurring in seeds during wine fermentation is necessary to build better physiochemical models.

59 BASIC BIOLOGICAL SCIENCES↗

Measuring Multicomponent Adsorption of Tracer Gases on Natural Zeolites

A natural clinoptilolite sample near the Nevada National Security Site was obtained to study adsorption and retardation on gas transport. Of interest is understanding the competition for adsorption sites that may reduce tracer gas adsorption relative to single-component measurements, which may be affected by the multi-scale pore structure of clinoptilolite. Clinoptilolite has three distinct domains of pore size distributions ranging from nanometers to micrometers: micropores with 0.4–0.7 nm diameters, measured on powders by CO2 adsorption at 273 K, representing the zeolite cages; mesopores with 4–200 nm diameters, observed using liquid nitrogen adsorption at 77 K; and macropores with 300–1000 nm diameters, measured by mercury injection on rock chips (~ 100 mesh), likely representing the microfractures. These pore size distributions are consistent with X-ray computed tomography (CT) and focused ion beam scanning electron microscope (FIB-SEM) images, which are used to construct the three-dimensional (3D) pore network to be used in future gas transport modeling. To quantify tracer gas adsorption in this multi-scale pore structure and multicomponent gas species environment, natural zeolite samples initially in equilibrium in air were exposed to a mixture of tracer gases. As the tracer gases diffuse and adsorb in the sample, the remaining tracer gases outside the sample fractionate. Using a quadrupole mass spectrometer to quantify this fractionation, the degree of adsorption of tracer gases in the multicomponent gas environment and multi-scale pore structure is assessed. The major finding is that Kr reaches equilibrium much faster than Xe in the presence of ambient air, which leads to more Kr uptake than Xe over limited exposure periods. When the clinoptilolite chips were exposed to humid air, the adsorption capability decreases significantly for both Xe and Kr with relative humidity (RH) as low as 3%. Both Xe and Kr reaches equilibrium faster at higher RH. The different, unexpected, adsorption behavior for Xe and Kr is due to their kinetic diameters similar to the micropores in clinoptilolite which makes it harder for Xe to access compared to Kr.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling the 4D discharge of lithium-ion batteries with a multiscale time-dependent deep learning framework

The lithium-ion battery (LIB) field is moving towards the direction of investigating spatially resolved physical phenomena in the 3D porous microstructure of electrodes. These pore-scale simulations give new insights into the local dynamics of lithiation/de-lithiation and charge transport, Nevertheless, the computational time of these simulations limits the integration of these models in optimization workflows of cycling conditions or electrode manufacturing processes. Machine learning models present a way of assessing in real-time the performance of materials. While several successful techniques for replicating simulations with machine learning have been proposed, this case study presents a more demanding problem, due to the necessity of understanding the behavior of heterogeneous 3D local data, as it evolves in time: this poses both a scientific and a technical challenge. To this end, we propose an autoregressive multiscale convolutional neural network model to predict relevant quantities at the pore-scale in the solid phase: the lithium concentration (in the active material) and potential (in the active material and carbon binder). Here, these are ultimately used to reconstruct the battery discharge curve. 3D images of the electrode microstructures are the input to the network, trained with a dataset of finite element method simulations to predict the discharge behavior of the cathode side in lithium ion batteries. We propose this machine learning model as a proof-of-concept of the applicability of multiscale networks for time-dependent physics problems. The trained model exhibits very high accuracy (with errors lower than 2 %) in forecasting the discharge behavior of new unseen cathodes.

25 ENERGY STORAGE↗

Data-scarce surrogate modeling of shock-induced pore collapse process

Understanding the mechanisms of shock-induced pore collapse is of great interest in various disciplines in sciences and engineering, including materials science, biological sciences, and geophysics. However, numerical modeling of the complex pore collapse processes can be costly. To this end, a strong need exists to develop surrogate models for generating economic predictions of pore collapse processes. Here, in this work, we study the use of a data-driven reduced-order model, namely dynamic mode decomposition, and a deep generative model, namely conditional generative adversarial networks, to resemble the numerical simulations of the pore collapse process at representative training shock pressures. Since the simulations are expensive, the training data are scarce, which makes training an accurate surrogate model challenging. To overcome the difficulties posed by the complex physics phenomena, we make several crucial treatments to the plain original form of the methods to increase the capability of approximating and predicting the dynamics. In particular, physics information is used as indicators or conditional inputs to guide the prediction. In realizing these methods, the training of each dynamic mode composition model takes only around 30 s on CPU. In contrast, training a generative adversarial network model takes 8 h on GPU. Moreover, using dynamic mode decomposition, the final-time relative error is around 0.3% in the reproductive cases. We also demonstrate the predictive power of the methods at unseen testing shock pressures, where the error ranges from 1.3 to 5% in the interpolatory cases and 8 to 9% in extrapolatory cases.

97 MATHEMATICS AND COMPUTING↗