An analytical homogenization method for heterogeneous porous materials

The objective of this work is to develop an analytical homogenization method to estimate the effective mechanical properties of fluid-filled porous media with periodic microstructure. The method is based on the equivalent inclusion concept of homogenization applied earlier for solid–solid mixture. It is assumed that porous media are described by the poroelastic constitutive law developed by Biot where porosity is a material parameter. By solving the governing equations of poroelasticity in Fourier transformed domain, the relation between periodic strain and eigenstrain in porous media is established. This relation is subsequently used in an average consistency condition involving both solid and fluid phase stresses and strains. The geometry of the porous microstructure is captured in the g-integral. This homogenization method can also be applied to estimate the equivalent properties of solid–fluid mixture where a pure solid and fluid can be modeled by assuming very low and high porosity, respectively. Several examples are considered to establish this new method by comparing with other existing analytical and numerical methods of homogenization. As an application, poroelastic properties of cortical bone fibril are estimated and compared with previously computed values.     

Determination of Nanoparticle Macrotransport Coefficients from Pore Scale Processes

Nanoparticle transport in porous media is modeled using a hierarchical set of differential equations corresponding to pore scale and macroscale. At the pore scale, movement and interaction of a single particle with a solid matrix is modeled using the advection–dispersion–sorption equation. A single nanoparticle entering the space encounters viscous, diffusion and surface forces. Surface forces (electrostatic and van der Waals forces) between nanoparticles and mineral grains appear as sorption propensity on solid matrix boundary condition. These local events are then transformed into a macroscale continuum by imposing periodic boundary conditions for contiguous unit cells representing porous media and using a scheme of moment analysis. At the macroscale, propagation and retention of particles are characterized by three position-independent coefficients: mean nanoparticle velocity vector \({\bar{\mathbf{U}}}^*\), macroscopic dispersion coefficient \({\bar{\mathbf{D}}}^*\), and mean nanoparticle retention rate constant \({\bar{K}}^*\). The modeling results are validated with a set of nanoparticle transport tests in porous microchips. We also present simulations of realistic porous media, where an actual image of sandstone samples is processed into binary tones. The representative unit cells are constructed from the resulting binary image by searching for areas within the sample with maximum similarities to the whole sample in terms of porosity and specific surface area, which are found to show strong correlations with the resulting \({\bar{\mathbf{U}}}^*\) and \({\bar{K}}^*\), respectively.     

Effective Behavior Near Clogging in Upscaled Equations for Non-isothermal Reactive Porous Media Flow

For a non-isothermal reactive flow process, effective properties such as permeability and heat conductivity change as the underlying pore structure evolves. We investigate changes of the effective properties for a two-dimensional periodic porous medium as the grain geometry changes. We consider specific grain shapes and study the evolution by solving the cell problems numerically for an upscaled model derived in Bringedal et al. (Transp Porous Media 114(2):371–393, 2016. doi: 10.1007/s11242-015-0530-9). In particular, we focus on the limit behavior near clogging. The effective heat conductivities are compared to common porosity-weighted volume averaging approximations, and we find that geometric averages perform better than arithmetic and harmonic for isotropic media, while the optimal choice for anisotropic media depends on the degree and direction of the anisotropy. An approximate analytical expression is found to perform well for the isotropic effective heat conductivity. The permeability is compared to some commonly used approaches focusing on the limiting behavior near clogging, where a fitted power law is found to behave reasonably well. The resulting macroscale equations are tested on a case where the geochemical reactions cause pore clogging and a corresponding change in the flow and transport behavior at Darcy scale. As pores clog the flow paths shift away, while heat conduction increases in regions with lower porosity.     

On the Permeability of Fractal Tube Bundles

The permeability of a porous medium is strongly affected by its local geometry and connectivity, the size distribution of the solid inclusions, and the pores available for flow. Since direct measurements of the permeability are time consuming and require experiments that are not always possible, the reliable theoretical assessment of the permeability based on the medium structural characteristics alone is of importance. When the porosity approaches unity, the permeability?Cporosity relationships represented by the Kozeny?CCarman equations and Archie??s law predict that permeability tends to infinity and thus they yield unrealistic results if specific area of the porous media does not tend to zero. The aim of this article is the evaluation of the relationships between porosity and permeability for a set of fractal models with porosity approaching unity and a finite permeability. It is shown that the tube bundles generated by finite iterations of the corresponding geometric fractals can be used to model porous media where the permeability?Cporosity relationships are derived analytically. Several examples of the tube bundles are constructed, and the relevance of the derived permeability?Cporosity relationships is discussed in connection with the permeability measurements of highly porous metal foams reported in the literature.     

Numerical Solution and Sensitivity Analysis of Pore Liquid Pressure and Cake Porosity to Model Parameters

Filter cake formation is important in groundwater and oil wells where drilling contains suspended mud particles. The accumulation of these mud particles on the borehole wall creates a pressure drop in the well. Furthermore, the migration of colloidal particles into adjacent porous rock could damage the formation and cause productivity decline. In this study, numerical solutions for pore liquid pressure variation across the cake with variable total stress and associated porosity variation are obtained. Mass equations for captured and suspended particles are averaged along the mud cake thickness, taking into account conditions on the cake surface and at the filter septum. The variability of total stress in soil consolidations problem is considered to determine the pore liquid pressure along the mud cake thickness. Then, the relation between porosity and pressure is studied to determine the mud cake porosity. Experimental data obtained by various researchers is used to compare and test the validity of numerical solutions to develop guidelines for model applications. Results show that the pore liquid pressure increases with the decrease of membrane impedance value (i.e. less pervious membrane). Also, the pressure profile has a cubic function of dimensionless cake thickness. The conclusions from the sensitivity analysis conducted in this study agree with earlier conclusions.     

Tomographic Study of Internal Erosion of Particle Flows in Porous Media

In particle-laden flows through porous media, porosity and permeability are significantly affected by the deposition and erosion of particles. Experiments show that the permeability evolution of a porous medium with respect to a particle suspension is not smooth, but rather exhibits significant jumps followed by longer periods of continuous permeability decrease. Their origin seems to be related to internal flow path reorganization by avalanches of deposited material due to erosion inside the porous medium. We apply neutron tomography to resolve the spatiotemporal evolution of the pore space during clogging and unclogging to prove the hypothesis of flow path reorganization behind the permeability jumps. This mechanistic understanding of clogging phenomena is relevant for a number of applications from oil production to filters or suffosion as the mechanisms behind sinkhole formation.     

Permeability of Pore Networks Under Compaction

The characteristic pore length fixes the scale of permeability of a porous medium. For pore networks undergoing strong random compaction, this length becomes singular at transition porosities, revealing a change in the microstructure of the porespace controlling the transport. Nodal balances and lattice Boltzmann simulations of flow in pore networks under compaction show that the scaling between permeability and porosity changes near the transition porosities. Simulation results are compared with experimental permeability data from transparent two-dimensional micromodels of networks decorated with the same pore size distribution. Permeability?Cporosity data of media undergoing smooth compaction is well described by a single power law. Under strong compaction, however, the scaling between permeability and porosity is possible by traits only, the scaling exponent changes notably at given transition porosities. These behaviors are common to a wealth of permeability?Cporosity data thus far unexplained.     

Boundary conditions at the interface between fluid layer and fibrous medium

The flow through a channel partially filled with fibrous porous medium was analyzed to investigate the interfacial boundary conditions. The fibrous medium was modeled as a periodic array of circular cylinders, in a hexagonal arrangement, using the boundary element method. The area and volume average methods were applied to relate the pore scale to the representative elementary volume scale. The permeability of the modeled fibrous medium was calculated from the Darcy's law with the volume‐averaged Darcy velocity. The slip coefficient, interfacial velocity, effective viscosity and shear jump coefficients at the interface were obtained with the averaged velocities at various permeabilities or Darcy numbers. Copyright © 2008 John Wiley & Sons, Ltd.     

Modified Particle Detachment Model for Colloidal Transport in Porous Media

Particle detachment from the rock during suspension transport in porous media was widely observed in laboratory corefloods and for flows in natural reservoirs. A new mathematical model for detachment of particles is based on mechanical equilibrium of a particle positioned on the internal cake or matrix surface in the pore space. The torque balance of drag, electrostatic, lifting and gravity forces, acting on the particle from the matrix and the moving fluid, is considered. The torque balance determines maximum retention concentration during the particle capture. The particle torque equilibrium is determined by the dimensionless ratio between the drag and normal forces acting on the particle. The maximum retention function of the dimensionless ratio (dislodging number) closes system of governing equations for colloid transport with particle release. One-dimensional problem of coreflooding by suspension accounting for limited particle retention, controlled by the torque sum, allows for exact solution under the assumptions of constant filtration coefficient and porosity. The explicit formulae permit the calculation of the model parameters (maximum retention concentration, filtration and formation damage coefficients) from the history of the pressure drop across the core during suspension injection. The values for maximum retention concentration, as obtained from two coreflood tests, have been matched with those calculated by the torque balance on the micro scale.     

Evaluation of the Darcy's law performance for two-fluid flow hydrodynamics in a particle debris bed using a lattice-Boltzmann model

In the present paper, multiphase flow dynamics in a porous medium are analyzed by employing the lattice-Boltzmann modeling approach. A two-dimensional formulation of a lattice-Boltzmann model, using a D2Q9 scheme, is used. Results of the FlowLab code simulation for single phase flow in porous media and for two-phase flow in a channel are compared with analytical solutions. Excellent agreement is obtained. Additionally, fluid-fluid interaction and fluid-solid interaction (wettability) are modeled and examined. Calculations are performed to simulate two-fluid dynamics in porous media, in a wide range of physical parameters of porous media and flowing fluids. It is shown that the model is capable of determining the minimum body force needed for the nonwetting fluid to percolate through the porous medium. Dependence of the force on the pore size, and geometry, as well as on the saturation of the nonwetting fluid is predicted by the model. In these simulations, the results obtained for the relative permeability coefficients indicate the validity of the reciprocity for the two coupling terms in the modified Darcy's law equations. Implication of the simulation results on two-fluid flow hydrodynamics in a decay-heated particle debris bed is discussed. Received on 1 December 1999     

Heat-Mass-Transport and Thermal Stresses in Porous Charring Materials

The aim of this paper is to develop a method of asymptotic averaging for processes occurring in porous charring materials under high temperatures. The advantage of the method is the ability to calculate not only averaged macrocharacteristics of the processes, namely internal gas generation, filtration and deforming processes, but also microcharacteristics, such as microstresses in phases of charring material, gas velocity in a pore, etc. To determine microcharacteristics, the method allows us to formulate special mathematical problems on a periodic cell. To calculate macrocharacteristics, such as pore pressure of filtrating gas, rate of charring and macrostresses, with the help of asymptotic averaging method, averaged global equations are formulated. Here effective characteristics of porous medium (gas permeability coefficient, rate of charring, elasticity modulus, thermal expansion coefficient) are determined not empirically, as in most works on porous materials, but on the basis of solving the local problems. Solution of these problems over the periodic cell allows us to derive analytically the law of the Darcy type for a gas phase flow in porous media, to obtain an expression for intensive mass transfer between solid and gas phases, to set the form of constitutive relations for charring porous media, and also to calculate microstresses in a vicinity of a growing pore. As an example of solving a global averaged problem, the problem on one-sided high-temperature heating of a plate made of epoxy binder has been solved numerically.     

Averaged Momentum Equation for Flow Through a Nonhomogenenous Porous Structure

This paper addresses the derivation of the macroscopic momentum equation for flow through a nonhomogeneous porous matrix, with reference to dendritic structures characterized by evolving heterogeneities. A weighted averaging procedure, applied to the local Stokes' equations, shows that the heterogeneous form of the Darcy's law explicitly involves the porosity gradients. These extra terms have to be considered under particular conditions, depending on the rate of geometry variations. In these cases, the local closure problem becomes extremely complex and the full solution is still out of reach. Using a simplified two-phase system with continuous porosity variations, we numerically analyze the limits where the usual closure problem can be retained to estimate the permeability of the structure.     

Impact of Non-uniform Properties on Governing Equations for Fluid Flows in Porous Media

The macroscopic governing equations of a compressible multicomponents flow with non-uniform viscosity and with mass withdrawal (due to heterogeneous reactions) in a porous medium are developed. The method of volume averaging was used to transform local (or microscopic) governing equations into averaged (or macroscopic) governing equations. The impacts of compressibility, non-uniform viscosity, and mass withdrawal on the form of the averaged equations and on the value of the macroscopic transport coefficients were investigated. The results showed that the averaged mass conservation equation is significantly affected by mass withdrawal when a specific criterion on the size of the domain is respected. The results also showed that the form of the averaged momentum equations is not affected by mass withdrawal, by compressibility effects or by non-uniform viscosity, provided that the Reynolds number at the pore level is small. Nonetheless, the velocity field is affected by the heterogeneous reaction via the averaged mass conservation equation, and also by viscosity variations due to the presence of the volume-averaged viscosity (which value changes with position) in the averaged momentum equations. A new closure variable definition was proposed to formulate the closure problem, which avoided the need to solve an integro-differential equation in the closure problem. This formulation was used to show that the permeability tensor only depends on the geometry of the porous medium. In other words, that tensor is independent on whether the fluid is compressible/incompressible, has uniform/non-uniform viscosities, and whether mass withdrawal due to heterogeneous reactions is present/absent.     

A linear dynamic model for a saturated porous medium

A linear isothermal dynamic model for a porous medium saturated by a Newtonian fluid is developed in the paper. In contrast to the mixture theory, the assumption of phase separation is avoided by introducing a single constitutive energy function for the porous medium. An important advantage of the proposed model is it can account for the couplings between the solid skeleton and the pore fluid. The mass and momentum balance equations are obtained according to the generalized mixture theory. Constitutive relations for the stress, the pore pressure are derived from the total free energy accounting for inter-phase interaction. In order to describe the momentum interaction between the fluid and the solid, a frequency independent Biot-type drag force model is introduced. A temporal variable porosity model with relaxation accounting for additional attenuation is introduced for the first time. The details of parameter estimation are discussed in the paper. It is demonstrated that all the material parameters in our model can be estimated from directly measurable phenomenological parameters. In terms of the equations of motion in the frequency domain, the wave velocities and the attenuations for the two P waves and one S wave are calculated. The influences of the porosity relaxation coefficient on the velocities and attenuation coefficients of the three waves of the porous medium are discussed in a numerical example.     

Surface porosity and permeability of porous media with a periodic microstructure

Various ways of determining the surface porosity, the relation between the porosity and the surface porosity and the representation of the permeability in terms of the characteristics of the microstructure of the porous medium are analyzed with reference to model porous media with a periodic microstructure. It is shown that it is necessary to distinguish between the geometric (scalar) and physical (tensor) suface porosities and that the geometric surface porosity, the physical surface porosity and the porosity are different characteristics of the porous medium.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 79–85, January–February, 1995.     

Motion equations of the dynamic theory of consolidation from the point of view of the theory of transient pipe flows

An elastic fluid-saturated porous medium is modeled as a bundle of parallel cylindrical tubes aligned in a direction parallel to the fluid movement. The pore space is filled with viscous compressible liquid. A cell model and the theory of transient pipe flow are used to derive one-dimensional governing equations in such media. All macroscopic constants in these equations are defined by the individual material constants of the fluid and solid. The interaction force includes an additional term not found in Biot's theory.     

Microdroplet absorption by two-layer porous media

Microdroplet absorption by two-layer porous media is studied both theoretically and experimentally. A two-dimensional model for liquid flow from a droplet into a porous medium is presented and veri.ed based on a simultaneous numerical solution of the Euler equations taking into account surface tension forces and the unsteady filtration equation. The effect of the structural parameters of the two-layer porous medium (pore size in the layers, and the thickness and porosity of the layers) on the droplet absorption is analyzed. It is shown that the presence of the second layer can have a significant effect on the droplet absorption rate and the liquid distribution in the medium. The pore size is found to be the main parameter that governs the effect of the second layer. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 121–130, January–February, 2007.     

Pore Network Analysis of Resistivity Index for Water-Wet Porous Media

Pore network analysis is used to investigate the effects of microscopic parameters of the pore structure such as pore geometry, pore-size distribution, pore space topology and fractal roughness porosity on resistivity index curves of strongly water-wet porous media. The pore structure is represented by a three-dimensional network of lamellar capillary tubes with fractal roughness features along their pore-walls. Oil-water drainage (conventional porous plate method) is simulated with a bond percolation-and-fractal roughness model without trapping of wetting fluid. The resistivity index, saturation exponent and capillary pressure are expressed as approximate functions of the pore network parameters by adopting some simplifying assumptions and using effective medium approximation, universal scaling laws of percolation theory and fractal geometry. Some new phenomenological models of resistivity index curves of porous media are derived. Finally, the eventual changes of resistivity index caused by the permanent entrapment of wetting fluid in the pore network are also studied.Resistivity index and saturation exponent are decreasing functions of the degree of correlation between pore volume and pore size as well as the width of the pore size distribution, whereas they are independent on the mean pore size. At low water saturations, the saturation exponent decreases or increases for pore systems of low or high fractal roughness porosity respectively, and obtains finite values only when the wetting fluid is not trapped in the pore network. The dependence of saturation exponent on water saturation weakens for strong correlation between pore volume and pore size, high network connectivity, medium pore-wall roughness porosity and medium width of the pore size distribution. The resistivity index can be described succesfully by generalized 3-parameter power functions of water saturation where the parameter values are related closely with the geometrical, topological and fractal properties of the pore structure.