Computer simulation of the effect of deformation on the morphology and flow properties of porous media

An Analytical Model for Hysteretic Pressure-Sensitive Permeability of Nanoporous Media

Hysteretic pressure-sensitive permeability of nanohybrids composed of substantial nanopores is critical to characterizing fluid flow through nanoporous media. Due to the nanoscale effect (gas slippage), complex and heterogeneous pore structures of nanoporous media, the essential controls on permeability hysteresis of nanohybrids are not determined. In this study, a hysteretic pressure sensitive permeability model for nitrogen flow through dry nanoporous media is proposed. The derived model takes into account the nanoscale effect and pore deformation due to effective stress. The model is validated by comparing it with the experimental data. The results show that the calculated permeability and porosity are consistent with the measured results with the maximum relative error of 6.08% and 0.5%, respectively. Moreover, the hysteretic pressure-sensitive permeability of nanohybrids is related to effective stress, gas slippage, pore microstructure parameters, grain quadrilateral angle, and the loss rate of grain quadrilateral angle. The nanoscale effect is crucial to the permeability of nanoporous media. In addition, as impacted by the comprehensive impact of multiple relevant influential parameters, permeability during the pressure unloading process is not a monotonous function but presents complicated shapes. The proposed model can explain, quantify, and predict the permeability hysteresis effect of nanoporous media reasonably well.

Influence of pore-size distributions and pore-wall mechanics on the mechanical behavior of cellular solids like aerogels

The pore-size distributions play a critical role in the determination of the properties of nanoporous cellular materials like aerogels. In this paper, we propose a micromechanical model, and by further designing artificial normal pore-size distributions, we inspect their effect on the macroscopic stress-strain curves. We show that the location of the mean pore size as well as the broadness of the distribution strongly affects the overall macroscopic behavior. Moreover, we also show that by using different damage criteria within the proposed model, the elastic, inelastic, and brittle nature of the macroscopic material can be captured. The damage criteria are based on the different modes of deformation in the pore walls, namely, elastic buckling, irreversible bending and brittle collapse under compression, and combined bending and stretching under tension. The proposed model approach serves as a reverse engineering tool to develop cellular solids with desired mechanical properties.

Numerical Simulation of the Aerosol Particle Motion in Granular Filters with Solid and Porous Granules

In this work, a study was carried out to compare the filtering and hydrodynamic properties of granular filters with solid spherical granules and spherical granules with modifications in the form of micropores. We used the discrete element method (DEM) to construct the geometry of the filters. Models of granular filters with spherical granules with diameters of 3, 4, and 5 mm, and with porosity values of 0.439, 0.466, and 0.477, respectively, were created. The results of the numerical simulation are in good agreement with the experimental data of other authors. We created models of granular filters containing micropores with different porosity values (0.158–0.366) in order to study the micropores’ effect on the aerosol motion. The study showed that micropores contribute to a decrease in hydrodynamic resistance and an increase in particle deposition efficiency. There is also a maximum limiting value of the granule microporosity for a given aerosol particle diameter when a further increase in microporosity leads to a decrease in the deposition efficiency.

A Coarse Grained Model for Viscoelastic Solids in Discrete Multiphysics Simulations

Viscoelastic bonds intended for Discrete Multiphysics (DMP) models are developed to allow the study of viscoelastic particles with arbitrary shape and mechanical inhomogeneity that are relevant to the pharmaceutical sector and that have not been addressed by the Discrete Element Method (DEM). The model is applied to encapsulate particles with a soft outer shell due, for example, to the partial ingress of moisture. This was validated by the simulation of spherical homogeneous linear elastic and viscoelastic particles. The method is based on forming a particle from an assembly of beads connected by springs or springs and dashpots that allow the sub-surface stress fields to be computed, and hence an accurate description of the gross deformation. It is computationally more expensive than DEM, but could be used to define more effective interaction laws.

Lattice Boltzmann-Discrete Element Modeling Simulation of SCC Flowing Process for Rock-Filled Concrete

Since invented in 2003, rock-filled concrete (RFC) has gained much attention and has been successfully applied in more and more civil and hydraulic projects in China. This study developed a numerical framework to simulate self-compacting concrete (SCC) flows in the voids among rocks of RFC, which couples the lattice Boltzmann method and discrete element method (DEM). The multiple-relaxation-time scheme is used to simulate self-compacting mortar (SCM) for better stability while the motion of coarse aggregates in SCC is simulated with DEM. The immersed moving boundary method is incorporated to deal with the interactions between coarse aggregates and SCM. After validation, the coupled framework is applied to study SCC flows in a single channel and in porous media with multi-channels. A passing factor PF was proposed and calculated to describe quantitatively the passing ability of SCC through a single channel. The study found that jamming of SCC occurs when the ratio Ar of the gap width to particle diameter is smaller than 2.0 and the jamming risk increases with solid particles fraction while the passing ability has a weak relation with the pressure gradient. Further, SCC flow is self-tuning in porous media with multi-channels and it is prone to go through larger or wider gaps. Exceeded existence of narrow gaps will significantly increase the jamming risk.

Pore-scale characteristics of multiphase flow in heterogeneous porous media using the lattice Boltzmann method

This study provides a pore-scale investigation of two-phase flow dynamics during primary drainage in a realistic heterogeneous rock sample. Using the lattice Boltzmann (LB) method, a series of three-dimensional (3D) immiscible displacement simulations are conducted and three typical flow patterns are identified and mapped on the capillary number (Ca)-viscosity ratio(M) phase diagram. We then investigate the effect of the viscosity ratio and capillary number on fluid saturation patterns and displacement stability in Tuscaloosa sandstone, which is taken from the Cranfield site. The dependence of the evolution of saturation, location of the displacement front, 3D displacement patterns and length of the center of mass of the invading fluid on the viscosity ratio and capillary number have been delineated. To gain a quantitative insight into the characteristics of the invasion morphology in 3D porous media, the fractal dimension D_f of the non-wetting phase displacement patterns during drainage has been computed for various viscosity ratios and capillary numbers. The logarithmic dependence of D_f on invading phase saturation appears to be the same for various capillary numbers and viscosity ratios and follows a universal relation.

Fluid flow through packings of elastic shells

Fluid transport in porous materials is commonly studied in geological samples (soil, sediments, etc.) or idealized systems, but the fluid flow through compacted granular materials, consisting of substantially strained granules, remains relatively unexplored. As a step toward filling this gap, we study a model of liquid transport in packings of deformable elastic shells using finite-element and lattice-Boltzmann methods. We find that the fluid flow abruptly vanishes as the porosity of the material falls below a critical value, and the flow obstruction exhibits features of a percolation transition. We further show that the fluid flow can be captured by a simplified permeability model in which the complex porous material is replaced by a collection of disordered capillaries, which are distributed and shaped by the percolation transition. To that end, we numerically explore the divergence of hydraulic tortuosity τ_{H} and the decrease of a hydraulic radius R_{h} as the percolation threshold is approached. We interpret our results in terms of scaling predictions derived from the percolation theory applied to random packings of spheres.

Image-based modeling of gas adsorption and deformation in porous media

Understanding adsorption of CO₂ in porous formations is crucial to its sequestration in geological formations. We describe a model for adsorption of CO₂ and the deformation that it induces in a sandstone formation over wide ranges of temperature and pressure. The model couples the thermodynamics of sorption with elastic deformation of the solid. Finite-element computations are then used in order to compute CO₂ adsorption isotherms along with the induced strain in the formation. We also compute the Darcy permeability of the porous medium using the lattice-Boltzmann method. All the computations are carried out with a three-dimensional image of a core sample from Mt. Simon sandstone, the target porous formation for a pilot CO₂ sequestration project that is currently being carried out by Illinois State Geological Survey. Thus, no assumptions are made regarding the shape and sizes of the pore throats and pore bodies. The computed CO₂ sorption isotherm at 195 K is in excellent agreement with our experimental data. The computed permeability is also in good agreement with the measurement. As a further test we also compute the sorption isotherm of N₂ in the same formation at 77.3 K, and show that it is also in good agreement with our experimental data. The model is capable of predicting adsorption of CO₂ (or any other gas for that matter) in porous formations at high pressures and temperatures. Thus, it is used to study the effect of hydrostatic pressure on adsorption and deformation of the porous formation under various conditions. We find that the effect of the confining pressure is more prominent at higher temperatures. Also computed is the depth-dependence of the capacity of the formation for CO₂ adsorption, along with the induced volumetric strain.