Pore-scale Model Research Articles

Heat transfer analysis of subcooled flow boiling in copper foam helical coiled heat exchanger – A pore-scale numerical study

High performance thermal energy systems are crucial in fulfilling the growing demand of thermal energy. Consequently, numerous heat transfer enhancements have been proposed. Some of the promising methods are boiling flow heat transfer, helical coiled tube and metal foam insert. Although the performance of each strategy has been extensively studied, based on our literature review, the effectiveness of combining these strategies has not been investigated. This missing information may impede further development of efficient heat transfer systems. This study is therefore performed to investigate subcooled flow boiling heat transfer in helical tube with copper foam insert. A high-resolution three-dimensional pore scale model that considers the shape, porosity and pore density of the metal foam is developed and validated against the published experimental results. Using the model, the effect of inlet temperature, mass flux, wall heat flux, porosity, pore density and half-filled metal foam are studied. The results suggest that operating parameters have mild effects on the heat transfer performance. In contrast, porosity and pore density show significant influence in heat transfer coefficient and pressure drop. It has also been found that total heat transfer surface area and base area are better in characterizing heat transfer in pore-scale level than porosity and pore density at pore-scale level. Moreover, implementing metal foam at outer half of the tube is observed to offer slightly better heat transfer than implementing it at the inner side by around 4.8%. Overall, this study indicates that helical metal foam tubes offer the best performance at mass flux of 200 kg/m2∙s with performance index of 0.72, while straight tubes with 0.90 porosity has the highest performance index of 0.66. Empirical correlations to relate dimensionless parameters to the Nusselt number and friction factor is generated. Finally, the model developed in this study can serve to provide guidelines for the upcoming studies in boiling in metal foam structure and assist in designing metal foam heat exchangers.

Upscaling reactive transport models from pore-scale to continuum-scale using deep learning method

Reactive transport modeling of subsurface environments plays an important role in addressing critical problems of geochemical processes, such as dissolution and precipitation of minerals. Current transport models for porous media span various scales, ranging from pore-scale to continuum-scale. In this study, we established an upscaling method connecting pore-scale and continuum-scale models by employing a deep learning methodology of Convolutional Neural Networks (CNNs). We applied Darcy-Brinkmann-Stokes (DBS) method to simulate the fluid flow and reactive transport in pore-scale models, which would act as constituents of a continuum-scale model. The datasets of spatial pore distribution of subvolume samples were used as the input for the deep learning model, and the continuum (Darcy)-scale parameters such as permeability, effective surface area, and effective diffusion coefficient were figured out as outputs (i.e., labels). By applying the trained models of the subvolumes in the entire sample volume, we generated the initial field of porosity, permeability, effective diffusion coefficient, and effective surface area for continuum-scale simulation of a mineral dissolution problem. We took an acid dissolution case as an example to utilize the outcomes of trained deep learning models as input data in the continuum-scale simulation. This work presents a comprehensive upscaling workflow, as bridging the findings of microscale simulations to the continuum-scale simulations of a reactive transport problem.

Development of stochastically reconstructed 3D porous media micromodels using additive manufacturing: numerical and experimental validation

We propose an integrated methodology for the design and fabrication of 3D micromodels that are suitable for the pore-scale study of transport processes in macroporous materials. The micromodels, that bear the pore-scale characteristics of sandstone, such as porosity, mean pore size, etc, are designed following a stochastic reconstruction algorithm that allows for fine-tuning the porosity and the correlation length of the spatial distribution of the solid material. We then construct a series of 3D micromodels at very fine resolution (i.e. 16μ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$16\\,\\upmu $$\\end{document}m) using a state-of-the-art 3D printing infrastructure, specifically a ProJet MJP3600 3D printer, that utilizes the Material Jetting technology. Within the technical constraints of the 3D printer resolution, the fabricated micromodels represent scaled-up replicas of natural sandstones, that are suitable for the study of the scaling between the permeability, the porosity and the mean pore size. The REV- and pore-scale characteristics of the resulting physical micromodels are recovered using a combination of X-ray micro-CT and microfluidic studies. The experimental results are then compared with single-phase flow simulations at pore-scale and geostatistic models in order to determine the effects of the design parameters on the intrinsic permeability and the spatial correlation of the velocity profile. Our numerical and experimental measurements reveal an excellent match between the properties of the designed and fabricated 3D domains, thus demonstrating the robustness of the proposed methodology for the construction of 3D micromodels with fine-tuned and well-controlled pore-scale characteristics. Furthermore, a pore-scale numerical study over a wider range of 3D digital domain realizations reveals a very good match of the measured permeabilities with the predictions of the Kozeny–Carman formulation based on a single control parameter, k0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$k_0$$\\end{document}, that is found to have a practically constant value for porosities ϕ≥0.2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\phi \\ge 0.2$$\\end{document}. This, in turn, enables us to customize the sample size to meet REV constraints, including enlarging pore morphology while considering the Reynolds number. It is also found that at lower porosities there is a significant increase in the fraction of the non-percolating pores, thus leading to different k0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$k_0$$\\end{document}, as the porosity approaches a numerically determined critical porosity value, ϕc\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\phi _c$$\\end{document}, where the domain is no longer percolating.

State of the Art on Fe Precipitation in Porous Media: Hydrogeochemical Processes and Evolving Parameters

The mixing of terrestrial groundwater and seawater creates dynamic reaction zones in intertidal areas, where land-derived Fe(II) is oxidized to Fe(III) and then precipitates as Fe hydroxides at the groundwater–seawater interface. These hydrogeochemical processes contribute to the formation of iron bands at the saltwater wedge (SW) and beneath the upper saline plume (USP). This study provides a comprehensive review of physical and geochemical processes at field scale in coastal areas, explores the impact of mineral precipitation on pore structure at pore scale, and synthesizes reactive transport modeling (RTM) approaches for illustrating continuum-scale soil physio-chemical parameters during the evolution of porous media. Upon this review, knowledge gaps and research needs are identified. Additionally, challenges and opportunities are presented. Therefore, we reach the conclusion that the incorporation of observational data into a comprehensive physico-mathematical model becomes imperative for capturing the pore-scale processes in porous media and their influence on groundwater flow and solute transport at large scales. Additionally, a synergistic approach, integrating pore-scale modeling and non-invasive imaging, is equally essential for providing detailed insights into intricate fluid–pore–solid interactions for future studies, as well as facilitating the development of regional engineering-scale models and physio-chemical coupled models with diverse applications in marine science and engineering.

Effects of Nafion content in the catalyst layer of PEMFC on the transport phenomenon among nanoscale particles

Nafion, as a key component in the catalyst layer, plays a crucial role in the mass transport and electrochemical reactions within the catalyst layer. In this study, a numerical method is employed to reconstruct the three-dimensional nanoscale microstructure of the catalyst layer, including pores, Nafion, carbon particles, and Pt phases. A pore-scale model in which the Knudsen effect is taken into account is used to simulate the impact of Nafion contents on the interparticle transport and electrochemical reactions within the catalyst layer. The results show that the pore-scale simulation based on the reconstructed model can clearly reveal the mass transport process and highly uneven distribution of each physical quantity in the nanoscale particles and pores; When the Nafion content is relatively higher, the proton current distribution at the Pt/Nafion interface is relatively more uniform and the proton supply is sufficient. However, there is an inverse correlation between the oxygen mole concentration on the Pt surface and the Nafion content, which is the dominant factor affecting the electrochemical reactions on the Pt particle surface. When the Nafion content is lower, the proton supply becomes the dominant factor in the electrochemical reaction on the surface of Pt particles due to the weak connectivity of Nafion between the agglomerates. The combined effect of oxygen and proton transport result in the optimal Nafion volume fractions of 0.296 and 0.262 at voltages of 0.8 V and 0.7 V, respectively. The optimal Nafion volume fraction ranges between 0.224 and 0.262 at 0.6 V. As the voltage decreases, the effect of oxygen transport becomes more significant, and the optimal Nafion content decreases accordingly.

Pore-scale modeling of effects of multiphase reactive transport on solid dissolution in porous media with structural heterogeneity

Fundamental understanding of the multiphase reactive transport with solid dissolution is crucial for a wide range of scientific and engineering problems. In this study, a pore-scale model is improved to study coupled processes of multiphase flow, mass transport, chemical reaction and solid dissolution in porous media with increasing structural heterogeneity. The results show that with the existence of multiphase flow, under the reactive transport conditions studied (Pe = 5, Da = 10; and Pe = 5, Da = 1), more compact dissolution pattern is generated due to the negative feedback between reactant transport and solid dissolution, while the feedback is positive in single-phase condition leading to wormhole dissolution. Besides, while it is found that the structural heterogeneity amplifies the wormhole dissolution in single-phase reactive transport processes or promotes the formation of fingering for multiphase flow, such heterogeneity effect is homogenized in the multiphase case. Finally, based on the pore-scale results, a new correlation between porosity, saturation and interfacial length is proposed which can be upscaled into the continuum-scale models.

Macroscopic transport characteristics in packed bed with a multi-physics inversion strategy for solid breeding blanket of HCCB TBM

The solid breeding blanket with typical packed bed structure is a key component for tritium breeding and thermal conversion in future fusion reactors. The pore scale model for multi-physical field solution of solid breeding blanket is fine but very expensive. The simplified packed-bed scale model is more favored in engineering especially for the system-level calculation of the blanket, but it needs to be inputted appropriate macroscopic transport parameters. In this work, the multi-physics calculations of blanket in Chinese HCCB TBM (Helium-Cooled Ceramic-Breeder Test Blanket Module) at pore scale and packed-bed scale were performed. It was found that the packed-bed scale model with traditional macroscopic transport parameters is not accurate enough, especially the prediction deviation of tritium concentration distribution could reach about 25 %. Hence, a novel multi-physics inversion strategy was proposed to predict the macroscopic transport parameters of the packed bed under the multi-field coupling conditions. The inversion results show that the convection of helium significantly enhances the heat and tritium diffusion in the packed bed, and the effective thermal conductivity and effective diffusivity could be increased by 18 % and 46 %, respectively. Then based on the results by multi-physics inversion strategy, the modified effective thermal conductivity and effective diffusivity correlation models for solid breeding blanket packed bed are constructed. With the modified correlation models, the relative prediction deviations of pressure drop, temperature distribution and concentration distribution are within 2.5 %, 0.6 % and 10 % respectively at different inlet velocities (0.05 m/s ∼ 0.25 m/s) and different inlet temperatures (300 K ∼ 900 K). The study in this paper would assist the thermal-hydraulics analysis of the solid breeding blanket in the fusion reactor.

Solid–liquid phase transition simulated by the lattice Boltzmann model: from pore scale to representative elementary volume scale

Purpose This paper aims to establish a lattice Boltzmann (LB) method for solid-liquid phase transition (SLPT) from the pore scale to the representative elementary volume (REV) scale. By applying this method, detailed information about heat transfer and phase change processes within the pores can be obtained, while also enabling the calculation of larger-scale SLPT problems, such as shell-and-tube phase change heat storage systems. Design/methodology/approach Three-dimensional (3D) pore-scale enthalpy-based LB model is developed. The computational input parameters at the REV scale are derived from calculations at the pore scale, ensuring consistency between the two scales. The approaches to reconstruct the 3D porous structure and determine the REV of metal foam were discussed. The implementation of conjugate heat transfer between the solid matrix and the solid−liquid phase change material (SLPCM) for the proposed model is developed. A simple REV-scale LB model under the local thermal nonequilibrium condition is presented. The method of bridging the gap between the pore-scale and REV-scale enthalpy-based LB models by the REV is given. Findings This coupled method facilitates detailed simulations of flow, heat transfer and phase change within pores. The approach holds promise for multiscale calculations in latent heat storage devices with porous structures. The SLPT of the heat sinks for electronic device thermal control was simulated as a case, demonstrating the efficiency of the present models in designing and optimizing SLPT devices. Originality/value A coupled pore-scale and REV-scale LB method as a numerical tool for investigating phase change in porous materials was developed. This innovative approach allows for the capture of details within pores while addressing computations over a large domain. The LB method for simulating SLPT from the pore scale to the REV scale was given. The proposed method addresses the conjugate heat transfer between the SLPCM and the solid matrix in the enthalpy-based LB model.

Thermal management of PV based on latent energy storage of composite phase change material: A system-level analysis with pore-scale model

Perovskite Solar Cell (PSC) have recently emerged as exciting new candidates of Photovatics (PVs) for solar-to-electrical energy conversion. Nevertheless, one huge obstacle to its commercialization is how to improve its tolerance to operation at elevated temperatures. To address this issue, Composite Phase Change Material (CPCM) incorporated with a porous skeleton structure is proposed to integrate with PVs. In this structure, an important theoretical problem is that the temperature dependence of PV and the temperature control behavior of CPCMs interact with each other (bidirectional relationship). In this work, a pore-scale lattice Boltzmann phase change model is established to describe the dynamic thermal behavior of CPCMs and explore the roles of the bidirectional relationship between PVs and CPCMs, the change of power conversion efficiency influenced by temperature variation and glass transmissivity is discussed by comparing the overall liquid fraction of CPCMs, and the average temperature and power conversion efficiency of PVs. Furthermore, the situation of time-varying solar irradiation is discussed to mimic real application conditions. Results indicate that the neglecting temperature dependency of PVs leads to underestimation of PV surface temperatures and overestimation of power conversion efficiencies, while glass transmissivity provides reverse effects. These phenomena not only happen in constant solar irradiation but also in varying solar irradiation, calling for increasing attention to designing PV/CPCMs for varying heat flux conditions. The findings of this work establish a comprehensive manner toward integrated design of CPCMs for PV thermal control.

A Mini Review on Electroosmotic Phenomena in Porous Media

Electroosmosis phenomenon in porous media finds widespread applications in various fields such as microfluidic systems, polymer electrolyte membrane fuel cells, oil and gas engineering, wastewater sludge dewatering and groundwater dynamics etc. Therefore, the electroosmotic flow mechanism in porous media have attracted broad interests from multiple disciplines. This paper provides an overview on the physical mechanisms and mathematical models for electroosmosis in porous media. The background of electric double layer theory and state-of-art research progress on pore-scale models for electroosmotic flow through porous media are reviewed. Two typical and significant research topics, electroosmosis under pressure coupling effect and nanoscale electroosmotic phenomena, are then focused on. The advances in theoretical analysis, numerical simulation and experimental measurements are summarized. Finally, the potential research directions for the electroosmotic flow in porous media are addressed out.

Pore-Scale Modeling of Pressure-Driven Flow and Spontaneous Imbibition in Fracturing-Shut-In-Flowback Process of Tight Oil Reservoirs

Tight oil reservoirs are characterized by multiple pore spaces where nano-micropores and multiscale fractures coexist, and each type of medium varies in scale, implying a tight coupling of multiscale fractures with matrix and giving rise to extremely complicated flow patterns. To further investigate its flow mechanism, we first construct three two-dimensional (2-D) fracture-pore geometry models based on microfocus computed tomography (CT) imaging of a typical tight rock. A pore-scale modeling workflow is thereafter developed using the Shan-Chen lattice Boltzmann model (SC-LBM) to simulate the pressure-driven flow and spontaneous imbibition. The influence of fracture-pore geometry on the pore-scale fluid exchange dynamics in the fracturing-shut-in-flowback process is clearly clarified. Results show that, for the porous medium model without fracture, the fracturing fluid can displace and replace some crude oil by spontaneous imbibition while a large amount of crude oil droplets remains unexploited away from the oil/water contact line, resulting in a low oil imbibition recovery. The injected fracturing fluid migrates into the deep position along the fractures, only a few entering the matrix pore space near the fractures. Small pores are the main channel for fracturing fluid to imbibe into the matrix pores, and the replaced crude oil droplets flow into fractures through large pores as intermittent or continuous pipe flow. The complexity of fracture network typically exerts a significant impact on the fluid exchange dynamics during pressure-driven flow and spontaneous imbibition. The more complex the fracture network, the larger the volume of fracturing fluid injected, the easier for oil droplets replaced from the matrix pores, and the more difficult for the flowback of fracturing fluid in the fracture-pore geometry model. Our understanding will provide a basis for explaining the underlying mechanisms of oil replacement by pressure-driven flow and spontaneous imbibition in the fracturing-shut-in-flowback process.

Pore-Scale Modeling of Heterogeneous Carbonate Rock Subjected to Modified Salinity Waterflooding

Pore-Scale Modeling of Heterogeneous Carbonate Rock Subjected to Modified Salinity Waterflooding

A parallel programming application of the A* algorithm in digital rock physics

Pore-scale analysis and characterization of reservoir rocks provide valuable information for the definition and management of underground hydrogen storage and CO2 storage or sequestration.This article presents an optimized implementation of the A* algorithm, the most popular pathfinding method in the presence of obstacles. In this context, the algorithm is applied to recognize the minimum length connected paths through each flow direction of 3D images of a porous medium representative of a reservoir rock. The identification of the main paths allows the characterization of the pore space and the calculation of fundamental petrophysical properties such as tortuosity and effective porosity, which can be used for permeability estimation. Compared to other algorithms available for pore-scale characterization, such as the pore centroid, A* provides a better approximation of the pore space available for the flow and, therefore, a reliable characterization of the petrophysical properties. On the other hand, path identification is significantly consuming in terms of time and memory. In this paper, an efficient and optimized implementation based on C++/OpenMP programming language is presented.The proposed implementation aims to the analysis of large-scale models profiting from parallelization, memory optimization, and enhanced managing of dead paths. Three test cases of increasing sizes are presented, to analyze the advantages and the disadvantages of the algorithm as the number of explored points increases. The 3D binary images analyzed are related to a synthetic domain (1 million voxels) and two actual sandstone samples (about 4 and 64 million voxels respectively). The code is validated against a Matlab serial implementation, showing a significant improvement in efficiency. Remarkable test cases of several millions of voxels were afforded, overcoming the memory and execution slowness issues. Moreover, the proposed implementation is suitable for large pore-scale models run in HPC environments.

Multiscale Model for Hydrogen Chloride Uptake on Aluminum for Spacecraft Fire Safety

An electrical fire within spacecraft environments often produces poisonous gases, such as hydrogen chloride (HCl). These gases adsorb on the interior materials of the spacecraft, primarily treated aluminum surfaces. A safe design of the spacecraft can be facilitated by modeling the spacecraft compartment. In this study, a multiscale model is developed to predict HCl uptake by anodized aluminum, which has a thick oxide layer covering its surface compared to other treatments. X-ray photoelectron spectroscopy revealed that aluminum chloride is found deep in the oxide layer of anodized aluminum, implying that the pore-scale transport and reaction of HCl within the oxide layer covering the aluminum surface must be considered. Hence, a multiscale approach is warranted. A pore-scale model is first developed, wherein the HCl reacts with the aluminum oxide pore walls to create an aluminum chloride layer. The HCl then diffuses through that layer for further reaction, until the layer (diffusion resistance) grows too large and the uptake of HCl stops. This pore-scale model is coupled with a reactor-scale model to predict overall HCl uptake. Results show that this new multiscale approach predicts the uptake of HCl far more accurately when compared to experimental data than previously used single (reactor)-scale models.

Pore-scale mushy layer modelling

Equations describing mushy systems, in which solid and liquid are described as a single continuum, have been extensively studied. Most studies of mushy layers have assumed them to be ‘ideal’, such that the liquid and solid were in perfect thermodynamic equilibrium. It has become possible to simulate flows of passive porous media at the pore scale, where liquid and solid are treated as separate continua. In this contribution, we study the simplest possible mushy layers at the pore scale, modelling a single straight cylindrical pore surrounded by a cylindrical annulus representing the solid matrix. Heat and solute can be exchanged at the solid–liquid boundary. We consider harmonic temperature and concentration perturbations and examine their transport rates due to advection and diffusion and the melting and solidification driven by this transport. We compare the results of this numerical model with those of a one-dimensional ideal mushy layer and with analytical solutions valid for ideal mushy layers for small temperature variations. We demonstrate that for small values of an appropriately defined Péclet number, the results of the pore-scale model agree well with ideal mushy layer theory for both transport rates and melting. As this Péclet number increases, the temperature and concentration profiles with radius within the pore differ significantly from constant, and the behaviour of the pore-scale model differs significantly from that of an ideal mushy layer. Some effects of mechanical dispersion arise naturally in our pore-scale model and are shown to be important at high Péclet number.

Multiscale modeling of smectite illitization in bentonite buffer of engineered barrier system

With the increasing usage of nuclear energy, how to properly dispose nuclear waste becomes a critical issue. In this study, a multiscale modeling approach combining the experimental findings is presented to address the illitization process, its impact on transport properties, and system behavior of bentonite buffer in engineered barrier systems (EBS). Through the pore-scale modeling, reactive transport properties such as illite generation rate and effective diffusion coefficient of potassium ion as a function of porosity and temperature are quantified by employing the findings of hydrothermal reaction experiments of Bentonil-WRK. The capability of pore-scale modeling has been developed based on the Darcy-Brinkmann-Stokes equation, involving the processes of smectite illitization and clay swelling. Obtained reactive transport properties are utilized as input parameters for the macroscale modeling to predict the long-term behavior of bentonite buffer in EBS. As such, this study involves the whole workflow of quantifying the reaction parameters of smectite illitization through the hydrothermal reaction experiments, and numerically modeling the reactive transport process of smectite illitization in bentonite buffer of EBS from pore-scale to macroscale. The presented multiscale modeling findings are expected to provide reliable solution for safe nuclear waste disposal with EBS.

Pore-scale modeling of preformed particle gel (PPG) extrusion in porous media

Preformed particle gel (PPG) is a promising profile control agent for enhanced oil recovery in heterogeneous reservoirs. For an improved sweep, the remigration behavior of PPG once blocked by pore throats is of crucial importance. In this paper, the process of a hyper-elastic particle passing through a pore throat is simulated, allowing investigation of the extrusion mechanism of PPG with deformation under a certain hydraulic pressure. The results reveal that PPG is subjected to a varying resistance during the process of entering, penetrating and leaving the throat. The maximal radial resistance occurs when the particle center arrives at the narrowest throat, whereas the maximal axial resistance occurs before that. In the meantime, it is shown that the particle gel bears over 80% body strain during the extrusion process through frictional pore throat. The hydraulic pressure needed for the PPG extrusion through the pore throat depends on the mechanical behaviors (hyper-elasticity of PPG, and friction coefficient of the PPG-throat interface) and geometrical properties (the particle size, and the shape of throat wall) of the PPG-throat system. Based on the simulation results, an empirical model is finally established to predict the required driving pressure. Eventually, the pore-scale model is upscaled and employed to attain effective pressure windows for profile-control in a reservoir scale formation with 10 times permeability heterogeneity.

Numerical investigation of effect of mechanical compression on the transport properties of fuel cell microporous layer using a pore-scale model

The microporous layer (MPL) plays an important role in water and thermal management of proton exchange membrane fuel cells (PEMFCs). An in-depth investigation of the mechanical compression effect on transport properties in the MPL can help optimize cell performance. In this work, the microstructure of the MPL is numerically reconstructed and the finite element method is applied to simulate mechanical behavior. Besides, the distribution of stress-strain, porosity, and pore size in the MPL under ten different levels of mechanical compression strains are studied. Lastly, the pore-scale model is employed to investigate the effective transport properties of the MPL as a function of compression strain. The analysis reveals that as the MPL strain increases from 0% to 40%, there is a 29% decrease in porosity, a 50% reduction in average pore diameter, a 60% decrease in effective gas diffusivity, a 100% increase in tortuosity, and an 80% increase in electrical and thermal conductivity. With the escalation of mechanical compression, both the magnitude and uniformity of stress-strain-displacement concurrently rise. Mechanical compression strains below 20% exhibit a lesser impact on transport properties. Beyond this threshold, exceeding the 20% compression strain point, mechanical stress assumes a critical role in influencing MPL transport properties.

Using X-ray computed tomography and pore-scale numerical modeling to study the role of heterogeneous rock surface wettability on hydrogen-brine two-phase flow in underground hydrogen storage

Underground hydrogen storage (UHS) is receiving increasing attention to address the challenges in hydrogen storage. A crucial aspect of UHS is understanding the transport of hydrogen in subsurface porous media. In this work, hydrogen core flooding experiments were conducted in a sandstone sample and then the pore-scale distribution of hydrogen and brine was visualized using high-resolution X-ray micro-computed tomography (micro-CT). After CT image processing, we measured the surface contact angles (CAs) on rock surfaces and found that the measured CAs followed a log-normal distribution. To investigate the influence of rock surface wettability on the transport properties in the hydrogen-brine-sandstone system, the lattice Boltzmann (LB) method was used to simulate two-phase flows with different surface CA distributions; hydrogen-brine relative permeability and capillary pressure curves through the primary drainage, imbibition, and secondary drainage processes were obtained by LB simulations. X-ray CT scanning showed that hydrogen resided in both large pores and small pores and pore throats after the primary drainage process (i.e., hydrogen displacing brine), whereas after the imbibition process (i.e., brine displacing hydrogen) hydrogen stayed primarily in large pores where the capillary pressure barriers were low. The LB two-phase flow modeling illustrated that water’s relative permeability increased whereas hydrogen’s relative permeability decreased when the core flooding moved from primary drainage to imbibition; in contrast, when the core flooding moved from imbibition to secondary drainage, water’s relative permeability decreased whereas hydrogen’s relative permeability increased. Decreasing the surface CA increased the capillary pressure and reduced the hydrogen residual rate, which indicates high hydrogen retrievability and thus is favorable for UHS practice. The change in CA’s standard deviation did not cause noticeable changes in water’s relative permeability curves, whereas it resulted in noticeable changes in hydrogen’s relative permeability curves. This work is the first study that utilizes X-ray micro-CT scanning and pore-scale multiphase flow modeling to quantitatively and comprehensively investigate hydrogen-brine two-phase flow properties under different rock surface wettability distributions. The research findings from this study will advance the understanding of hydrogen transport in an underground storage system and provide essential data for large-scale field studies in UHS.

The Effect of Porous Media Grain Size on non-Darcy Flow Behavior using Pore Scale Simulation

The impact of grain size of the porous media of the 3D images of a beadpack on the permeability values and non-Darcy flow behaviour is analyzed using pore-scale flow modelling utilizing the finite volume method. In this study, the sample used was a 3D image from random tight packing of spheres (beadpack) with a size of 250×250×250 voxels and a porosity of 0.36. Variation in grain size is carried out by upscaling the 3D image sample, so that this treatment affects the specific surface area value but does not change the porosity and tortuosity of the sample. Several variations of grain diameter include 100 μm, 125 μm, 150 μm, 175 μm, and 200 μm. In this study, we adopted PIMPLE algorithm, which combines SIMPLE and PISO algorithms to accomplish the Navier-Stokes equations of fluid dynamics in such porous media. The simulation results show that the permeability value exhibits an inverse relationship with the specific surface area. This correlation is in accordance with the Kozeny-Carman equation. In addition, the initiation of non-Darcy flow occurs at the reference length-based Reynolds number of 2.06 and the permeability-based Reynolds number of 0.049.