An Air Permeability Study of Anisotropic Glass Wool Fibrous Products

This article reports an experimental study of the permeability of fibrous porous media based on air flow. Two glass wools have been analyzed in three different flow directions to evaluate their anisotropic structure. The usual flow theories such as Darcy’s law, valid for granular materials, have been shown to apply also to these fibrous media which is a more complex structure (fiber arrangement, high porous media,…). The intrinsic permeabilities and the anisotropic factors have been determined from an air permeability measurement based on a standard test for industrial fabrics. The experimental results have been compared to different models, showing the limits of some simplistic models when applied to our highly porous structures. One model, Kyan’s model based on the flow around submerged objects, was, however, found to agree with our experimental data. This study also demonstrates the interesting use of an experimental technique, which is air permeability measurement (water permeability is not adapted), in the investigation of flow phenomena in a fibrous medium.     

Determination of permeability in fibrous porous media using the lattice Boltzmann method with application to PEM fuel cells

The lattice Boltzmann method (LBM) is used to simulate the flow through an idealized proton exchange membrane fuel cell (PEMFC) porous transport layer (PTL) geometry generated using a Monte Carlo method. Using the calculated flow field, Darcy's law is applied and the permeability is calculated. This process is applied in both through‐ and in‐plane directions of the paper as both of these permeability values are important in computational fluid dynamics models of PEMFCs. It is shown that the LBM can be used to determine permeability in a random porous media by solving the flow in the microstructure of the material. The permeability in the through‐ and in‐plane directions is shown to be different and the anisotropic nature of the geometry creates anisotropic permeability. It is also found that fiber arrangement plays a large role in the permeability of the PTL. New correlations are presented for in‐ and though‐plane permeabilities of fibrous porous media with (0.6<ε<0.8). Copyright © 2008 John Wiley & Sons, Ltd.     

Characterization of bubble mobility in channel flow with fibrous porous media walls

In composites processing, resin is introduced into a fibrous domain to cover all the empty spaces between the fibers. It is important to extract air bubbles from the domain before the resin solidifies. Failure to do so will entrap these voids in the final part, which is detrimental to its performance. Hence, there is a need to understand bubble motion in a fibrous porous domain in which the bubbles move with the resin in channels surrounded by fibrous walls. A rising bubble model is presented that consists of a single spherical void in a cylindrical axisymmetric two-phase domain of resin and air surrounded by porous media boundaries. The motion of a bubble in a channel flow with porous boundaries is modeled by replacing the walls with a slip velocity. Focus is on how the porous media permeability influences the bubble motion. A parameter called bubble mobility is defined as the ratio of bubble rise velocity to the resin free surface velocity. Results suggest that fabric permeability and fluid properties can be optimized to increase bubble mobility and ultimately lead to reduction in void content during composites processing.     

Smoothed Particle Hydrodynamics Modeling of Transverse Flow in Randomly Aligned Fibrous Porous Media

The Lagrangian smoothed particle hydrodynamics (SPH) method is employed to obtain a meso-/micro-scopic pore-scale insight into the transverse flow across the randomly aligned fibrous porous media in a 2D domain. Fluid is driven by an external body force, and a square domain with periodic boundary conditions imposed at both the streamwise and transverse flow direction is assumed. The porous matrix is established by randomly embedding a certain number of fibers in the square domain. Fibers are represented by position-fixed SPH particles, which exert viscous forces upon, and contribute to the density variations of, the nearby fluid particles. An additional repulsive force, similar in form to the 12-6 Lennard-Jones potential between atoms, is introduced to consider the no-penetrating restraint prescribed by the solid pore structure. This force is initiated from the fixed solid material particle and may act on its neighboring moving fluid particles. Fluid flow is visualized by plotting the local velocity vector field; the meandering fluid flow around the porous microstructures always follow the paths of least resistance. The simulated steady-state flow field is further used to calculate the macroscopic permeability. The dimensionless permeability (normalized by the squared characteristic dimension of the fiber cross section) exhibits an exponential dependence on the porosity within the intermediate porosity range, and the derived dimensionless permeability—porosity relation is found to have only minor dependence on either the relative arrangement condition among fibers or the fiber cross section (shape or area).     

Numerical Approach of the Permeability of a Macroporous Bioceramic with Interconnected Spherical Pores

Manufacturing a hybrid bone substitute requires a dynamic culture of the cells preliminarily seeded in a scaffold through a flow of physiological fluid. The velocity, pressure, and the distribution of fluid flow in this kind of macroporous medium are the important keys. Because of the difficulties in determining these parameters by experiment, a numerical approach has been chosen. One of the primary step of this study consists in the determination of permeability K. In this article, two types of structure of macroporous bioceramics are concerned. One is the interconnected pore spheres arranged either simple cubic, body-centered cubic or face-centered cubic systems. The other is the interconnected pore spheres randomly arranged. Based on Darcy??s law, the permeability K was calculated for many cases (type, porosity) by simulating the fluid flow through a small representative volume. These results are compared with some previous models such as Ergun, Carman?CKozeny, Rumpf?CGupte, and Du Plessis. The limits of Darcy??s law and the above-mentioned models have been determined using numerical simulation. The result showed that the porous media with spherical interconnected pores of BCC systems can be used to replace a complex random system in a range of porosity from 0.71 to 0.76 (i.e., porosity of our scaffolds). This assumption is validated for a pressure gradient lower the 1,000?Pa m^?C1 and a simple polynomial relation linking permeability and porosity (0.71?C0.76) has been established.     

A Simple Model for the Variation of Permeability due to Partial Saturation in Dual Scale Porous Media

The main focus of this work is to model macroscopically the effects of partial saturation upon the permeability of dual scale fibrous media made of fiber bundles when a Newtonian viscous fluid impregnates it. A new phenomenological model is proposed to explain the discrepancies between experimental pressure results and analytical predictions based on Darcy's law. This model incorporates the essential features of relative permeability but without the necessity of measuring saturation of the liquid for its prediction. The model is very relevant for the small scale industrial systems where a liquid is forced to flow through a fibrous porous medium. It requires four parameters. Two of them are the two permeability values based on the two length scales. One length scale is of the order of magnitude of the individual fiber radius and corresponds to the permeability of the completely staurated medium, the other is of the order of magnitude of the distance between the fiber bundles and corresponds to the permeability of the partially saturated medium. The other two parameters are the lengths of the two partially saturated regions of the flow domain. The two lengths of the partially saturated region and the permeability of the fully saturated flow domain can be directly measured from the experiments. The excellent agreement between the model and the experimental results of inlet pressure profile with respect to time suggests that this model may be used to describe the variation of the permeability behind a moving front in such porous media for correct pressure prediction. It may also be used to characterize the fibrous medium by determining the two different permeabilities and the relative importance of the unsaturated portion of the flow domain for a given architecture.     

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.     

Advanced Study About the Permeability for Micro-Porous Structures Using the Lattice Boltzmann Method

The micro flows through two-dimensional (2D) and three-dimensional (3D) granular porous media at various Knudsen numbers are studied by using the lattice Boltzmann method. For 2D cases, the correlation between the permeability, the porosity, and the Knudsen number is derived. For 2D cases, the correlation can estimate the permeability well except for the staggered square cylinder. The permeability of the porous media, which have the inclusions of different sizes, is calculated. For 3D cases, simulations for the uniform overlapping and non-uniform non-overlapping granular media are carried out. The results are compared with the correlation of previous study. The effect of rarefaction on the permeability is also discussed.     

A Potential Formulation of Non-Linear Models of Flow through Anisotropic Porous Media

A theoretical analysis, based on the search for a normal dissipation potential, is performed in order to generalize the empirical non-Darcy one-dimensional flow models to 3-D flows through anisotropic porous media. In an abstract framework, it is proven that a large number of heuristic non-linear equations governing the multidimensional flow through isotropic porous media can be derived starting from a potential strictly related to the mechanical power dissipated by the fluid. Such a formulation allows to define, for the tensor permeability case, a wide class of filtration models according to the Onsager's generalized theory of dissipative mechanical systems. A consistent generalization to anisotropic permeability case of polynomial flow models is proposed. Both primal and dual mixed variational formulations associated to the proposed quadratic and incomplete cubic flow models are introduced and discussed.     

Microstructural effects on the permeability of periodic fibrous porous media

An analytical–numerical approach is presented for computing the macroscopic permeability of fibrous porous media taking into account their microstructure. A finite element (FE) based model for viscous, incompressible flow through a regular array of cylinders/fibers is employed for predicting the permeability associated with this type of media. High resolution data, obtained from our simulations, are utilized for validating the commonly used semi-analytical models of drag relations from which the permeability is often derived. The effect of porosity, or volume fraction, on the macroscopic permeability is studied. Also microstructure parameters like particle shape, orientation and unit cell staggered angle are varied. The results are compared with the Carman–Kozeny (CK) equation and the Kozeny factor (often assumed to be constant) dependence on the microstructural parameters is reported and used as an attempt to predict a closed form relation for the permeability in a variety of structures, shapes and wide range of porosities.     

The transverse permeability of disordered fiber arrays: a statistical correlation in terms of the mean nearest interfiber spacing

In the porous media literature, unidirectional fibrous systems are broadly categorized as ordered or disordered. The former class, easily tractable for analysis purposes but limited in its relation to reality, involves square, hexagonal and various staggered arrays. The latter class involves everything else. While the dimensionless hydraulic permeability of ordered fibrous media is known to be a deterministic function of their porosity ϕ, the parameters affecting the permeability of disordered fiber arrays are not very well understood. The objective of this study is to computationally investigate flow across many unidirectional arrays of randomly placed fibers and derive a correlation between K and some measure of their microstructure. In the process, we explain the wide scatter in permeability values observed computationally as well as experimentally. This task is achieved using a parallel implementation of the Boundary Element Method (BEM). Over 600 simulations are carried out in two-dimensional geometries consisting of 576 fiber cross-sections placed within a square unit cell by a Monte Carlo procedure. The porosity varies from 0.45 to 0.90. The computed permeabilities are compared with earlier theoretical results and experimental data. Analysis of the computational results reveals that the permeability of disordered arrays with ϕ < 0.7 is reduced as the non-uniformity of the fiber distribution increases. This reduction can be substantial at low porosities. The key finding of this study is a direct correlation between K and the mean nearest inter-fiber spacing , the latter depending on the microstructure of the fibrous medium.     

Investigation of the characteristics of particulate flows through fibrous filters using the lattice Boltzmann method

A fibrous filter is one of the most common systems used to separate suspended particles from air.Two important factors(i.e.,the pressure drop and capture efficiency) are usually used to evaluate the performance of this type of filter.This study considers three two-dimensional arrangements of fibers(parallel,staggered,and random) to geometrically model fibrous media.The lattice Boltzmann method is employed to numerically simulate fluid flow through the filter.The Lagrangian form of the equation of motion of a particle is numerically solved to track the path of each particle in the flow field,where a one-way interaction between the fluid and particles is considered.The effects of pertinent parameters such as the fiber arrangement,solid volume fraction,particle-to-fiber diameter ratio,particle-to-fluid density ratio,Reynolds number,Stokes number,and size of the fibrous medium on the pressure drop and capture efficiency are studied.The obtained results are compared with existing empirical and theoretical findings and discussed.     

Practical Considerations for Selection of Representative Elementary Volumes for Fluid Permeability in Fibrous Porous Media

We investigate the lower bound of the area of a square-shaped representative elementary volume (REV) for the permeability tensor for transverse Stokes flow through randomly packed, parallel, and monodisperse cylinders. The investigation is significant to flow models using small calculation regions for fibrous porous media, such as modeling defect formation during directional solidification in the mushy zone of dendritic alloys. Using 90 ensembles of 1,000 domains, where each ensemble comprises domains with the same number and size of cylinders, we develop correlations between the permeability tensor invariants and macroscopic features of the domain. We find that for ensembles of domains with fewer than 200 cylinders, the eigenvectors of the permeability tensors exhibit preferential alignment with the domain axes, demonstrating that the estimated permeability is significantly affected by the periodic boundary conditions for these cases. Our results also suggest that the anisotropy of the permeability tensor may not be insignificant even for large sampling volumes. These results provide a practical lower bound for the calculation volumes used in permeability simulations in fibrous porous media, and also suggest that modelers should consider using an anisotropic tensor for small calculation volumes if phenomena such as channeling are important.     

Low Pressure Gas Percolation Characteristic in Ultra-low Permeability Porous Media

Low pressure gas percolation characteristic in ultra-low permeability porous media is investigated in this article through core flow experiments. The results show that the wall-slip layer covers more than 10% of the average porous channel radius on account of minimum pore size when the permeability is below 0.1 × 10^?3μ m ² order, and seepage behavior is contrasted to that in mid-high permeability pore media. When the gas pressure is not high enough, the flow regime turns into transitional flow instead of slip flow, and nonlinear relationship between the measured gas permeability and the reciprocal of average pressure exists. The gas measuring permeability experiment would be influenced by the non-linear relationship. If Klinkenberg-corrected method is applied to speculate the equivalent liquid permeability, the extrapolated value will become less or minus. Simultaneously, actual gas flow velocity at the outlet is beyond the calculated value with Klinkenberg formula. A new gas seepage model based on the general slip boundary condition is derived from the homogenization technique in this article. At last the flow model is examined to be suitable for representing the gas flow behavior in ultra-low permeability media and estimating the absolute permeability from single-point, steady-states measurements.     

Flow near the permeable boundary of a porous medium: An experimental investigation using LDA

 Fluid flow at the interface of a porous medium and an open channel is the governing phenomenon in a number of processes of industrial importance. Traditionally, this has been modeled by applying the Brinkman’s modification of Darcy’s law to obtain the velocity profile in terms of an additional parameter known as the “apparent viscosity” or the “slip coefficient”. To test this ad hoc approach, a detailed experimental investigation of the flow was conducted using Laser Doppler Anemometry (LDA) in the close vicinity of the permeable boundary of a porous medium. The porous medium used in the experiments consisted of a network of continuous glass strands woven together in a random fashion. A Hele–Shaw cell was partially filled with a fibrous preform such that an open channel flow is coupled with the Darcy flow inside the preform through the permeable interface of the preform. The open channel portion of the Hele–Shaw cell also acts as an ideal porous medium of known in-plane permeability which is much higher than the permeability of the fibrous porous medium. A viscous fluid is injected at a constant flow rate through the above arrangement and a saturated and steady flow is established through the cell. Using LDA, steady state velocity profiles are accurately measured by traversing across the cell in the direction perpendicular to the flow. A series of experiments were conducted in which fluid viscosity, flow rate, solid volume fraction of the porous medium and depth of the Hele–Shaw cell were varied. For each and every case in which the conditions for Hele–Shaw approximation were valid, the depth of the boundary layer zone or the screening length inside the fibrous preform was found to be of the order of the channel depth. This is much larger as compared to the Brinkman’s prediction of the screening length which is of the order of √K, where K is the permeability of the fibrous porous medium. Based on this finding, we modified the boundary condition in the Brinkman’s solution and found that the velocity profile results compared well with the experimental data for the planar geometry and the fibrous preforms for volume fractions of 7%, 14% and 21% for Hele–Shaw cell depths of 1.6 and 3.175 mm. For a cell depth of 4.8 cm, in which the Hele–Shaw approximation was not valid, the boundary layer thickness or the screening length was found to be less than the mold or channel depth but was still much larger than the Brinkman’s prediction. Received: 10 May 1996 / Accepted: 26 August 1996     

Finite-Difference Approximation for Fluid-Flow Simulation and Calculation of Permeability in Porous Media

We introduce a finite-difference method to simulate pore scale steady-state creeping fluid flow in porous media. First, a geometrical approximation is invoked to describe the interstitial space of grid-based images of porous media. Subsequently, a generalized Laplace equation is derived and solved to calculate fluid pressure and velocity distributions in the interstitial space domain. We use a previously validated lattice-Boltzmann method (LBM) as ground truth for modeling comparison purposes. Our method requires on average 17 % of the CPU time used by LBM to calculate permeability in the same pore-scale distributions. After grid refinement, calculations of permeability performed from velocity distributions converge with both methods, and our modeling results differ within 6 % from those yielded by LBM. However, without grid refinement, permeability calculations differ within 20 % from those yielded by LBM for the case of high-porosity rocks and by as much as 100 % in low-porosity and highly tortuous porous media. We confirm that grid refinement is essential to secure reliable results when modeling fluid flow in porous media. Without grid refinement, permeability results obtained with our modeling method are closer to converged results than those yielded by LBM in low-porosity and highly tortuous media. However, the accuracy of the presented model decreases in pores with elongated cross sections.     

A Fractal Model for Oil Transport in Tight Porous Media

Tight porous media are mainly composed of micro/nano-pores and throats, which leads to obvious microscale effect and nonlinear seepage characteristics. Based on the capillary bundle model and the fractal theory, a new nonlinear seepage equation was deduced, and a further fractal permeability model was obtained for oil transport in tight porous media by considering the effect of the boundary layer. The predictions of the model were then compared with experimental data to demonstrate that the model is valid. This model clarifies the oil transport mechanisms in tight porous media: the effective permeability is no longer a constant value and is governed by properties of tight porous media and oil. Furthermore, parameters influencing effective permeability were analyzed. The model can accurately present the seepage characteristics of the oil in tight porous media and provide a reliable basis for the development of unconventional reservoirs.     

Knudsen’s Permeability Correction for Tight Porous Media

Various flow regimes including Knudsen, transition, slip and viscous flows (Darcy’s law), as applied to flow of natural gas through porous conventional rocks, tight formations and shale systems, are investigated. Data from the Mesaverde formation in the United States are used to demonstrate that the permeability correction factors range generally between 1 and 10. However, there are instances where the corrections can be between 10 and 100 for gas flow with high Knudsen number in the transition flow regime, and especially in the Knudsen’s flow regime. The results are of practical interest as gas permeability in porous media can be more complex than that of liquid. The gas permeability is influenced by slippage of gas, which is a pressure-dependent parameter, commonly referred to as Klinkenberg’s effect. This phenomenon plays a substantial role in gas flow through porous media, especially in unconventional reservoirs with low permeability, such as tight sands, coal seams, and shale formations. A higher-order permeability correlation for gas flow called Knudsen’s permeability is studied. As opposed to Klinkenberg’s correlation, which is a first-order equation, Knudsen’s correlation is a second-order approximation. Even higher-order equations can be derived based on the concept used in developing this model. A plot of permeability correction factor versus Knudsen number gives a typecurve. This typecurve can be used to generalize the permeability correction in tight porous media. We conclude that Knudsen’s permeability correlation is more accurate than Klinkenberg’s model especially for extremely tight porous media with transition and free molecular flow regimes. The results from this study indicate that Klinkenberg’s model and various extensions developed throughout the past years underestimate the permeability correction especially for the case of fluid flow with the high Knudsen number.     

Porous Medium Flow and an Overlying Parallel Flow: PIV Interrogation Area and Overlaps, Interfacial Location, and Depth Ratio Effects

This paper presents an experimental work aimed at studying the effects of particle image velocimetry (PIV) interrogation area and overlaps, location of the interface, as well as depth ratios on the flow at the interface between a model porous medium and an overlying free flow. The porous media were modeled using square arrays of circular rods of diameter and porosity 0.88, filling fraction ranging from 0.47 to 0.75, and depth-to-porous medium pore ratio ranging from 5.75 to 13.69. Using a pressure-driven refractive-index matched viscous fluid, the bulk Reynolds number was kept approximately constant at a regime in which inertia was not a factor. PIV measurements were made across various streamwise-transverse planes of the test section. For the present tests, it was observed that PIV interrogation area (IA) and overlap effects on the interfacial velocities are negligible when the IA sizes in dimensionless units ranged from 0.017 to 0.145 in flow parameters and 0.036 to 0.300 in porous media parameters. Other dimensionless slip parameters are however significantly affected. Interfacial slip parameters of porous media models were found to change by as much as 120 % with change in the interfacial location. The interfacial location sensitivity was also found to be dependent on the direction of deviation, the type of porous medium, and depth ratios. Volume averaged results showed that for flows over models of porous media, the depth-to-porous medium pore ratio effects are more prominent compared with the filling fraction effects, for both two- and three-dimensional porous media.     

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.