Calculating the Anisotropic Permeability of Porous Media Using the Lattice Boltzmann Method and X-ray Computed Tomography

A lattice Boltzmann (LB) method is developed in this article in a combination with X-ray computed tomography to simulate fluid flow at pore scale in order to calculate the anisotropic permeability of porous media. The binary 3D structures of porous materials were acquired by X-ray computed tomography at a resolution of a few microns, and the reconstructed 3D porous structures were then combined with the LB model to calculate their permeability tensor based on the simulated velocity field at pore scale. The flow is driven by pressure gradients imposed in different directions. Two porous media, one gas diffusion porous layer used in fuel cells industry and glass beads, were simulated. For both media, we investigated the relationship between their anisotropic permeability and porosity. The results indicate that the LB model is efficient to simulate pore-scale flow in porous media, and capable of giving a good estimate of the anisotropic permeability for both media. The calculated permeability is in good agreement with the measured date; the relationship between the permeability and porosity for the two media is well described by the Kozeny–Carman equation. For the gas diffusion layer, the simulated results showed that its permeability in one direction could be one order of magnitude higher than those in other two directions. The simulation was based on the single-relaxation time LB model, and we showed that by properly choosing the relaxation time, it could give similar results to those obtained using the multiple-relaxation time (MRT) LB method, but with only one third of the computational costs of MRTLB model.     

Investigation of thermal dispersion and intra-pore turbulent heat flux in porous media

In the present study, the importance of the thermal dispersion and the turbulent heat flux in porous media and their effects on the macroscopic distribution of thermal energy are investigated. To this end, turbulent flow and heat transfer within five unit-cells mimicking porous media are solved using large eddy simulation. It is shown that the thermal dispersion and the turbulent heat flux are negligible as compared to the convection term in the macroscopic energy equation. When further scrutinizing this equation, it is revealed that except for the longitudinal components of the thermal dispersion, the other components of thermal dispersion and turbulent heat flux may be neglected away from the boundaries as compared to the interfacial heat transfer. Visualizations of vortices show that the size of the turbulence structures within the cells is of the same order as the size of the pores; therefore, the turbulent heat flux is limited to the intra-pore level. Finally, a discussion is provided on the accuracy of the gradient type diffusion model commonly used for turbulent heat flux in porous media in the absence of macroscopic turbulence. It is shown that the intra-pore turbulence does not affect the macroscopic transport of thermal energy within the porous media studied.     

Analysis of the moisture diffusion transfer through fibrous porous membrane used for waterproof breathable fabrics

In this paper, we propose a new fractal model to determine the moisture effective diffusivity of porous membrane such as expanded polytetrafluorethylene membrane, by taking account of both parallel and perpendicular channels to diffusion flow direction. With the consideration of both the Knudsen and bulk diffusion effect, a relationship between micro-structural parameters and effective moisture diffusivity is deduced. The effective moisture diffusivities predicted by the present fractal model are compared with moisture diffusion experiment data and calculated values obtained from other theoretical models.     

Upscaling Forchheimer law

We investigate the high velocity flow in heterogeneous porous media. The model is obtained by upscaling the flow at the heterogeneity scale where the Forchheimer law is assumed to be valid. We use the method of multiple scale expansions, which gives rigorously the macroscopic behaviour without any prerequisite on the form of the macroscopic equations. We show that Forchheimer law does not generally survive upscaling. The macroscopic flow law is strongly non-linear and anisotropic. A 2-point Padé approximation of the flow law in the form of a Forchheimer law is given. However, this approximation is generally poor. These results are illustrated in two particular cases: a layered composite porous media and a composite constituted by a square array of circular porous inclusions embedded in a porous matrix. We show that non-linearities are sources of anisotropy.     

Upscaling of Nonwetting Phase Residual Transport in Porous Media: A Network Approach

Based on the volume averaging method, a macroscopic model is developed for the upscaling of NAPL transport in a porous medium idealised by a network model. Under the assumption of local mass non-equilibrium, a macroscopic equation involving a dispersion tensor, additional convective terms and a linear form for the interfacial mass flux is obtained. The resolution of the two local closure problems obtained allow the determination of the local properties without adjustable parmeters. These problems are solved in a semi-analytical, semi-numerical manner on the network. The originality of this work is the association of the upscaling by volume averaging method with the network approach. The local properties, including the dispersion tensor and the mass exchange coefficient, can therefore be calculated over a large number of pore-bodies and pore-throats in a computationaly tractable manner, thus leading to more significant results. Results are presented for 3D, spatially periodic models of porous media.     

Effect of Mean Network Coordination Number on Dispersivity Characteristics

In this study, we investigate the role of topology on the macroscopic (centimeter scale) dispersion characteristics derived from pore-network models. We consider 3D random porous networks extracted from a regular cubic lattice with coordination number distributed in accordance with real porous structures. We use physically consistent rules including ideal mixing in pore bodies, molecular diffusion, and Taylor dispersion in pore throats to simulate transport at the pore-scale level. Fundamental properties of porous networks are based on statistical distributions of basic parameters. Theoretical calculations demonstrate strong correspondence with data obtained from numerical experiments. For low coordination numbers, we observe normal transport in porous networks. Anomalous effects expressed by tailing in concentration evolution are seen for higher coordination numbers. We find that the mean network coordination number has significant influence on averaged characteristics of porous networks such as geometric and hydraulic dispersivity, while other topological properties are of less significance. We give an explicit formula that describes the decrease of geometric dispersivity with growing mean coordination number. The results demonstrate the importance of network topology for models for flow and transport in porous media.     

Modelling turbulent and dispersion heat fluxes in turbulent porous medium flows using the resolved LES data

Resolved large eddy simulations (LESs) of turbulent conjugate heat transfer in porous media are performed by the lattice Boltzmann method (LBM) for modelling turbulent and dispersion heat flux terms of the double-averaged energy equation. The considered porous structures are square rod arrays, staggered cube arrays and body centred cubic foam. In the LBM, the double-distribution function method which solves the distribution functions for the velocity and the internal energy is used. For the velocity and thermal fields, the D3Q27 multiple-relaxation-time method and the regularized D3Q19 single-relaxation-time method are applied, respectively. A priori tests using the LES data suggest that the trends of the sum of the dispersion and volume-averaged turbulent heat fluxes can be well captured by the second order gradient diffusion model.     

Analysis of hydrodynamic and thermal dispersion in porous media by means of a local approach

A pore scale analysis is implemented in this numerical study to investigate the behavior of microscopic inertia and thermal dispersion in a porous medium with a periodic structure. The macroscopic characteristics of the transport phenomena are evaluated with an averaging technique of the controlling variables at a pore scale level in an elementary cell of the porous structure. The Darcy–Forchheimer model describes the fluid motion through the porous medium while the continuity and Navier–Stokes equations are applied within the unit cell. An average energy equation is employed for the thermal part of the porous medium. The macroscopic pressure loss is computed in order to evaluate the dominant microscopic inertial effects. Local fluctuations of velocity and temperature at the pore scale are instrumental in the quantification of the thermal dispersion through the total effective thermal diffusivity. The numerical results demonstrate that microscopic inertia contributes significantly to the magnitude of the macroscopic pressure loss, in some instances with as much as 70%. Depending on the nature of the porous medium, the thermal dispersion may have a marked bearing on the heat transfer, particularly in the streamwise direction for a highly conducting fluid and certain values of the Peclet number.     

Magnetic Resonance Imaging of Hydrodynamic Dispersion in a Saturated Porous Medium

By using nuclear magnetic resonance imaging (NMRI) we have been able to analyse dispersion at the microscopic scale during steady-state flow through water-saturated glass beads. The flow rate through the porous medium was chosen high enough in order to neglect the influence of molecular diffusion on dispersion. Velocity statistics were measured, by NMRI, within slices of increasing thickness perpendicular to the direction of flow. It took more than two bead diameters before a representative elementary volume (REV) for the mean velocity was reached. This was in a region in the middle of the column that was not influenced by the boundary conditions. There the velocity variance decreased exponentially as a function of the slice thickness, due we consider to the formation of an interconnecting streamline network. The exponential decrease in the velocity variance reflects the transition from a local pattern of stochastic–convective flow to a convective–dispersion regime at the scale of the REV. We found that the point-like preferential influx and efflux boundary condition increased velocity variances and thus enhanced longitudinal hydrodynamic dispersion. Using the transverse correlation length of longitudinal velocity variance, we derived a mean transverse dispersivity that agreed well with Saffmans (1959) model. So we have been able to provide for the first time a direct observation verification of a part of Saffmans (1959) conjectures. By NMRI we observed this value to be independent of the observation scale of the slice thickness.     

Measurement of 3D pore-scale flow in index-matched porous media

We present the experimental analysis of fluid flow at the pore-scale of a transparent porous medium with matched refractive indices of the solid and liquid phases. The planar laser-induced fluorescence (PLIF) technique described is the first to simultaneously visualize the 3D pore-scale flow of two immiscible liquid phases in porous media. Through the application of a highly precise index matching method and the employment of up-to-date CCD imaging hardware, the system features a high spatial resolution and sampling rate. The method was used to study the dispersion of a tracer dye in single-phase flow and the displacement of oil by water in an imbibition process.     

Numerical analysis of fluid flow and mass transfer in a channel with a porous bottom wall

The flow of a solution between parallel plates is considered. The bottom plate is porous, while the top one is an impermeable solid. A computer program based on the control volume approach was developed to analyse the flow and concentration fields. The effects of the slip at the porous wall on the velocity and particle concentration distributions were investigated. It was observed that as the slip increases, the concentration on the porous wall decreases and the maximum velocity moves towards the porous wall. The concentration on the porous wall increases in the flow direction. This increase in the particle concentration along the porous wall may cause a reduction of the porosity and hence a variation in the suction rate along the porous wall. In order to take this effect into account, a linearly varying transverse velocity along the porous wall was considered. The results were compared with the data available in the literature.     

Permeability mapping in porous media by magnetization prepared centric-scan SPRITE

The ability of porous media to transmit fluids is commonly referred to as permeability. The concept of permeability is central for hydrocarbon recovery from petroleum reservoirs and for studies of groundwater flow in aquifers. Spatially resolved measurements of permeability are of great significance for fluid dynamics studies. A convenient concept of local Darcy’s law is suggested for parallel flow systems. The product of porosity and mean velocity images in the plane across the average flow direction is directly proportional to permeability. Single Point Ramped Imaging with T ₁ Enhancement (SPRITE) permits reliable quantification of local fluid content and flow in porous media. It is particularly advantageous for reservoir rocks characterized by fast magnetic relaxation of a saturating fluid. Velocity encoding using the Cotts pulsed field gradient scheme improves the accuracy of measured flow parameters. The method is illustrated through measurements of 2D permeability maps in a capillary bundle, glass bead packs and composite sandstone samples.     

Upscaling LBM-TPM simulation approach of Darcy and non-Darcy fluid flow in deformable,heterogeneous porous media

This paper presents a numerical approach for the simulation of fluid flow through porous media by proposing a theoretical and numerical meso-to-macro multiscale framework, which combines the advantages of the lattice Boltzmann method (LBM) with the continuum Theory of Porous Media (TPM) to efficiently and accurately model fluid transport in heterogeneous porous media. In particular, LBM presents an alternative to experiments by studying the flow from a mesoscopic perspective, which in turn, allows the derivation of the material parameters needed for simulating the flow in the macroscopic TPM model. In this work, a meso-macro hierarchic upscaling scheme is applied to investigate the deformation-dependent intrinsic permeability properties and the Darcy/non-Darcy fluid flow regime. Concerning the mesoscale, the intrinsic permeability of the porous domain is computed by means of the LBM model at the first stage. Subsequently, deformation of the medium takes place in furtherance of determining the relation of the aforementioned deformation dependency. Thereupon, these findings are input into the TPM model in order to compute the primary unknown variables, where special focus is laid on the stability challenges in the compaction and near compaction states. With respect to the criteria of non-Darcy fluid flow, the conditions of its onset, i.e. the induced pressure gradient and mean fluid flow velocity, are computed as well using the LBM solver and conveyed afterwards to the macroscopic TPM model. Herein, the non-Darcy intrinsic permeability has been investigated in the TPM approach based on the Forchheimer equation. Simulations done on a synthetic porous micro-structure show that the combined framework proved to stand well between the two approaches.     

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.     

Free Flow at the Interface of Porous Surfaces: A Generalization of the Taylor Brush Configuration

A solution to the problem of shallow laminar water flow above a porous surface is essential when modeling phenomena such as erosion, resuspension, and mass transfer between the porous media and the flow above it. Previous studies proposed theoretical, experimental, and numerical insight with no single general solution to the problem. Many studies have used the Brinkman equation, while others showed that it does not represent the actual interface flow conditions. In this paper we show that the interface macroscopic velocity can be accurately modeled by introducing a modification to the Brinkman equation. A moving average approach was proved to be successful when choosing the correct representative elementary volume and comparing the macroscopic solution with the average microscopic flow. As the size of the representative elementary volume was found to be equal to the product of the square root of the permeability and an exponential function of the porosity, a general solution is now available for any brush configuration. Given the properties of the porous media (porosity and permeability), the flow height and its driving force, a complete macroscopic solution of the interface flow is obtained.     

Dependency of Tortuosity and Permeability of Porous Media on Directional Distribution of Pore Voids

Single-phase fluid flow in porous media is usually direction dependent owing to the tortuosity associated with the internal structures of materials that exhibit inherent anisotropy. This article presents an approach to determine the tortuosity and permeability of porous materials using a structural measure quantifying the anisotropic distribution of pore voids. The approach uses a volume averaging method through which the macroscopic tortuosity tensor is related to both the average porosity and the directional distribution of pore spaces. The permeability tensor is derived from the macroscopic momentum balance equation of fluid in a porous medium and expressed as a function of the tortuosity tensor and the internal structure of the material. The analytical results generally agree with experimental data in the literature.     

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.     

Acoustical characterization of fluid-saturated porous media with local heterogeneities: Theory and application

A linear dynamic model of fully saturated porous media with local (either microscopic or mesoscopic) heterogeneities is developed within the context of Biot’s theory of poroelasticity. Viscoporoelastic behavior associated with local fluid flow is characterized by the notion of the dynamic compatibility condition on the interface between the solid and the fluid. Complex, frequency-dependent material parameters characterizing the viscoporoelasticity are derived. The complex properties can be obtained through determining the quasi-static poroelastic parameters, the properties of individual constituents, and the relaxation time of the dynamic compatibility condition on the interface. Relationships among various quasi-static poroelastic parameters are developed. It is shown that local fluid flow mechanism is significant only in the porous media with local heterogeneities. The relaxation time of the compatibility condition on the interface depends upon the details of local structure of porous media that control local fluid pressure diffusion. The new model is used to describe the velocity dispersion and attenuation in fully saturated porous media. The proposed model provides a theoretical framework to simulate the acoustical behavior of fully saturated porous media over a wide range of frequencies without making any explicit assumption about the structure of local heterogeneities.     

A Lattice Gas Automaton Simulation of the Nonlinear Diffusion Equation: A Model for Moisture Flow in Unsaturated Porous Media

We investigate a two-dimensional lattice gas automaton (LGA) for simulating the nonlinear diffusion equation in a random heterogeneous structure. The utilility of the LGA for computation of nonlinear diffusion arises from the fact that, the diffusion coefficient in the LGA depends on the local density of fluid particles which statistically determines the collision rate and thus, the mean free path of the particles at the microscopic scale. The LGA may therefore be used as a physical analogue to simulate moisture flow in unsaturated porous media. The capability of the LGA to account for unsaturated flow is tested through a set of numerical experiments simulating one-dimensional infiltration in a simplified semi-infinite homogenous isotropic porous material. Different mechanisms of interactions are used between the fluid and the solid phase to simulate various fluid–solid interfaces. The heterogeneous medium, initially at low density is submitted to a steep density gradient by continuously injecting fluid particles at high concentration and zero velocity along one face of the model. The propagation of the infiltration front is visualized at different time steps through concentration profiles parallel to the applied concentration gradient and the infiltration rate is measured continuously until steady-state flow is reached. The numerical results show close agreement with the classical theory of flow in unsaturated porous media. The cumulative absorption exhibits the expected t ^1/2 dependence. The evolution of the effective diffusion coefficient with the particle concentration is estimated from the measured density profiles for the various porous materials. Depending on the applied fluid–solid interactions, the macroscopic effective diffusivity may vary by more than two orders of magnitude with density.