A general treatment for non-Darcy film condensation within a porous medium in the presence of gravity and forced flow

Non-Darcy film condensation over a vertical flat plate within a porous medium is considered. The Forchheimer extended Darcy model is adopted to account for the non-Darcy effects on film condensation in the presence of both gravity and externally forced flow. A general similarity transformation is proposed upon introducing a modified Peclet number based on the total velocity of condensate, resulting from both gravitational force and externally forced flow. This general treatment makes it possible to obtain all possible similarity solutions including the asymptotic results in the four different limiting regimes, namely, Darcy forced convection regime, Forchheimer forced convection regime, Darcy body force predominant regime and Forchheimer body force predominant regime. Appropriate dimensionless groups for distinguishing these asymptotic regimes are found to be the micro-scale Grashof and Reynolds numbers based on the square root of the permeability of the porous medium. Correspondingly, the non-Darcy effect on the heat transfer rate are investigated in terms of these micro-scale dimensionless numbers.     

Prediction of Non-Darcy Coefficients for Inertial Flows Through the Castlegate Sandstone Using Image-Based Modeling

Near wellbore flow in high rate gas wells shows the deviation from Darcy??s law that is typical for high Reynolds number flows, and prediction requires an accurate estimate of the non-Darcy coefficient (?? factor). This numerical investigation addresses the issues of predicting non-Darcy coefficients for a realistic porous media. A CT-image of real porous medium (Castlegate Sandstone) was obtained at a resolution of 7.57???m. The segmented image provides a voxel map of pore-grain space that is used as the computational domain for the lattice Boltzmann method (LBM) based flow simulations. Results are obtained for pressure-driven flow in the above-mentioned porous media in all directions at increasing Reynolds number to capture the transition from the Darcy regime as well as quantitatively predict the macroscopic parameters such as absolute permeability and ?? factor (Forchheimer coefficient). Comparison of numerical results against experimental data and other existing correlations is also presented. It is inferred that for a well-resolved realistic porous media images, LBM can be a useful computational tool for predicting macroscopic porous media properties such as permeability and ?? factor.     

Analytical quantification of coefficients in the Ergun equation for fluid friction in a packed bed

A new analytical derivation for momentum transport during laminar flow through granular porous media is discussed and some of its implied results described. In the very low Reynolds number regime fully developed laminar flow is assumed and in the higher laminar Reynolds number regime the Forchheimer (non-Darcy) effect is modelled through reference to form drag induced by the solid constituents of the porous medium. The results are compared to the Ergun equation, which is empirically based on experimental measurements, and the correspondence is shown to be remarkably close.     

A Criterion for Non-Darcy Flow in Porous Media

Non-Darcy behavior is important for describing fluid flow in porous media in situations where high velocity occurs. A criterion to identify the beginning of non-Darcy flow is needed. Two types of criteria, the Reynolds number and the Forchheimer number, have been used in the past for identifying the beginning of non-Darcy flow. Because each of these criteria has different versions of definitions, consistent results cannot be achieved. Based on a review of previous work, the Forchheimer number is revised and recommended here as a criterion for identifying non-Darcy flow in porous media. Physically, this revised Forchheimer number has the advantage of clear meaning and wide applicability. It equals the ratio of pressure drop caused by liquid–solid interactions to that by viscous resistance. It is directly related to the non-Darcy effect. Forchheimer numbers are experimentally determined for nitrogen flow in Dakota sandstone, Indiana limestone and Berea sandstone at flowrates varying four orders of magnitude. These results indicate that superficial velocity in the rocks increases non-linearly with the Forchheimer number. The critical Forchheimer number for non-Darcy flow is expressed in terms of the critical non-Darcy effect. Considering a 10% non-Darcy effect, the critical Forchheimer number would be 0.11.     

Inertia Effects in High-Rate Flow Through Heterogeneous Porous Media

The paper deals with the effects of large scale permeability–heterogeneity on flows at high velocities through porous media. The media is made of a large number of homogeneous blocks where the flow is assumed to be governed by the Forchheimer equation with a constant inertial coefficient. By assuming the validity of the Forchheimer equation at the large scale, an effective inertial coefficient is deduced from numerical simulations. Different media are investigated: serial-layers, parallel-layers and correlated media. The numerical results show that: (i) for the serial-layers, the effective inertial coefficient is independent of the Reynolds number and decreases when the variance and the mean permeability ratio increases; (ii) for the parallel-layers and the correlated media, the effective inertial coefficient is function of the Reynolds number and increases when the variance and the mean permeability ratio increases. Theoretical relationships are proposed for the inertial coefficient as function of the Reynolds number and the characteristics of the media.     

On Determining the Percolation Reynolds Number and the Characteristic Linear Dimension for Ideal and Fictitious Porous Media

Various versions of representations of the percolation Reynolds number for porous media with isotropic and anisotropic flow properties are considered. The formulas are derived and the variants are analyzed with reference to model porous media with a periodic microstructure formed by systems of capillaries and packings consisting of spheres of constant diameter (ideal and fictitious porous media, respectively). A generalization of the Kozeny formula is given for determining the capillary diameter in an ideal porous medium equivalent to a fictitious medium with respect to permeability and porosity and it is shown that the capillary diameter is nonuniquely determined. Relations for recalculating values of the Reynolds number determined by means of formulas proposed earlier are given and it is shown that taking the microstructure of porous media into account, as proposed in [1, 2], makes it possible to explain the large scatter of the numerical values of the Reynolds number in processing the experimental data.     

Optical measurements of pore geometry and fluid velocity in a bed of irregularly packed spheres

Imaging methods are proposed for the characterisation of liquid flows through transparent porous media of matched refractive index. The methods are based on the analysis of laser-illuminated slices, and specialized for the case in which the porous medium is composed of irregularly packed spheres. They include algorithms for the reconstruction of the three-dimensional (3D) sphere arrangement based on a laser scan of the packed bed, particle tracking velocimetry applied to the motions of micro-tracers in a laser-illuminated plane, and techniques for the co-registration of geometry and velocity measurements acquired from different slices. The methods are applied to a cylindrical flow cell filled with mono-sized spheres and operated at Reynolds number Re = 28. The data produced include the full 3D geometry of the packed spheres assembly, the 2D fluid velocity field in the axial centre-plane of the flow cell, and the corresponding porosity and velocity distributions.     

Predicting pressure drop in open-cell foams by adopting Forchheimer number

Interconnected struts arranged in 3-D foam structures pose a challenge in understanding fluid flow, which is significantly different from that in traditional porous media. Different flow regimes (Darcy, transition and weak inertia regimes) and thus, different flow laws in open-cell foams are used. The impact of characteristic lengths’ choices based on both, morphological and hydraulic parameters on flow law formulation has been studied. Ambiguities in definitions and measurements of several key parameters have been shown and limitations in the use of some parameters have been pointed out.An equivalent Reynolds number in the form of Forchheimer number (F_o) has been proposed to establish the friction factor relationship in order to avoid any morphological ambiguities. This number takes into account hydraulic characteristics of viscous and inertia regimes simultaneously. It has been observed that when F_o < 0.1, the flow through open-cell foams remains in the Darcy regime while the occurrence of weak inertia regime dominates when F_o > 1. Transition regime occurs in a narrow range of flow velocity when 0.1 < F_o < 1. The limits of transition for regime identification are found to be independent of foam morphologies. The form drag coefficient varies in relation with foam morphological parameters and is not a “universal” constant.Empirical correlations have been derived to predict hydraulic characteristics and friction factor data for different strut shapes and porosities. An excellent agreement has been obtained between predicted and numerical/experimental flow data.     

Flow-Induced Deformations Within Random Packed Beds of Spheres

Low Reynolds number flow-induced alterations of permeability of random packing of mono-sized spheres is studied. The number of spheres is several thousands and the porosities ranges between 0.4 and 0.6. The change of permeability is obtained for elastic deformations of the positions of the spheres using either of two methods. Each sphere is elastically attached to single points or the spheres that are connected via an elastic porous network. The system of spheres is divided into smaller volumes with Voronoi diagrams and the flow is derived by usage of a dual stream function. The local saturated flow fields are approximated as for close packed spheres and the overall flow pattern is obtained by minimising the dissipation rate of energy. The results show that the permeability for large random systems increases as a function of velocity and thus the deformation. The alteration is, however, much less than for two-dimensional cases like parallel cylinders. The relative increase in permeability becomes larger as the porosity increases from 0.4 to 0.6.     

Reduction in drag and vortex shedding frequency through porous sheath around a circular cylinder

A numerical study on the laminar vortex shedding and wake flow due to a porous‐wrapped solid circular cylinder has been made in this paper. The cylinder is horizontally placed, and is subjected to a uniform cross flow. The aim is to control the vortex shedding and drag force through a thin porous wrapper around a solid cylinder. The flow field is investigated for a wide range of Reynolds number in the laminar regime. The flow in the porous zone is governed by the Darcy–Brinkman–Forchheimer extended model and the Navier–Stokes equations in the fluid region. A control volume approach is adopted for computation of the governing equations along with a second‐order upwind scheme, which is used to discretize the convective terms inside the fluid region. The inclusion of a thin porous wrapper produces a significant reduction in drag and damps the oscillation compared with a solid cylinder. Dependence of Strouhal number and drag coefficient on porous layer thickness at different Reynolds number is analyzed. The dependence of Strouhal number and drag on the permeability of the medium is also examined. Copyright © 2009 John Wiley & Sons, Ltd.     

Polynomial Filtration Laws for Low Reynolds Number Flows Through Porous Media

In this study, we use the method of homogenization to develop a filtration law in porous media that includes the effects of inertia at finite Reynolds numbers. The result is much different than the empirically observed quadratic Forchheimer equation. First, the correction to Darcy’s law is initially cubic (not quadratic) for isotropic media. This is consistent with several other authors (Mei and Auriault, J Fluid Mech 222:647–663, 1991; Wodié and Levy, CR Acad Sci Paris t.312:157–161, 1991; Couland et al. J Fluid Mech 190:393–407, 1988; Rojas and Koplik, Phys Rev 58:4776–4782, 1988) who have solved the Navier–Stokes equations analytically and numerically. Second, the resulting filtration model is an infinite series polynomial in velocity, instead of a single corrective term to Darcy’s law. Although the model is only valid up to the local Reynolds number, at the most, of order 1, the findings are important from a fundamental perspective because it shows that the often-used quadratic Forchheimer equation is not a universal law for laminar flow, but rather an empirical one that is useful in a limited range of velocities. Moreover, as stated by Mei and Auriault (J Fluid Mech 222:647–663, 1991) and Barree and Conway (SPE Annual technical conference and exhibition, 2004), even if the quadratic model were valid at moderate Reynolds numbers in the laminar flow regime, then the permeability extrapolated on a Forchheimer plot would not be the intrinsic Darcy permeability. A major contribution of this study is that the coefficients of the polynomial law can be derived a priori, by solving sequential Stokes problems. In each case, the solution to the Stokes problem is used to calculate a coefficient in the polynomial, and the velocity field is an input of the forcing function, F, to subsequent problems. While numerical solutions must be utilized to compute each coefficient in the polynomial, these problems are much simpler and robust than solving the full Navier–Stokes equations.     

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.     

Intercomparison of boundary schemes in Lattice Boltzmann method for flow simulation in porous media

The Lattice Boltzmann method has been widely adopted to simulate flow in porous media. The choice of appropriate boundary schemes is essential to achieve simulation accuracy; however, the criteria for the most suitable boundary treatment in the simulation of flow in porous media flow remain unresolved. Here, three types of the most commonly used boundary conditions are tested: interpolation bounce back (IBB), partial saturated method (PSM), and immersed boundary method (IBM). The dimensionless drag of face-centered cubic (FCC) sphere array and the dimensionless permeability of a random closely packed (RCP) sphere array are calculated and compared at different viscosities and resolutions. In the FCC sphere array case where spheres are not contacted, the IBB and PSM exhibit the same accuracy and both are of the second-order convergence rate. The IBM is less accurate and is of the first-order convergence rate. In the RCP sphere array case where the spheres are contacted, the IBB shows finer results and a second-order convergence rate. PSM underestimates the dimensionless permeability and increases resolution only slightly improved the results of PSM. The IBM overestimates the dimensionless permeability. These results indicate that among the three methods, the IBB is the most accurate. The PSM has the same accuracy as the IBB when sediments are not contacted; however, it loses its accuracy in the simulation of flow in closely packed porous media. This work could serve as a benchmark for further research in choosing the most appropriate method in the simulation of flow in porous media.     

Effect of Confining Wall on Properties of Gas Flow Through Metal Foam: An Experimental Study

While the Darcy and Forchheimer relations are widely applied to determine the permeability and the form drag coefficient of open-cell metal foam, they both assume that the porous medium is infinite in all directions, i.e., large enough so that the effect of any confining walls is negligible. Many researchers, however, pay little or no attention to the size of metal foam samples in pressure-drop studies. The size of a foam sample perpendicular to the flow direction may be small enough such that wall effects are significant. This article experimentally investigates the wall effect on the permeability and form drag coefficient for two types of open-cell aluminum foam subjected to airflow entering the foam in the Forchheimer regime. The Forchheimer equation was recast in two different manners, which resulted in new non-dimensional numbers that correlated very well with the diameter of the foam samples measured in cells. The correlations are valid for a confining-tube-diameter-based Reynolds number ranging from approximately 13,000 to 105,000, and for diameters ranging from 12 to 36 cells and 24 to 60 cells for 10- and 20-pore per inch foam, respectively. The obtained correlations allow for determining pressure drop given only the velocity and the diameter of an aluminum foam sample.     

Analysis of Thermal Flow in a Rotating Porous U-Turn Duct Using Lattice Boltzmann Method

The lattice Boltzmann method is carried out to investigate the heat transfer enhancement in a U-turn duct which is partially filled with a porous media. The porous layer is inserted at the core of the duct and is modeled using the Brinkman–Forchheimer assumptions. In order to validate the results, first a channel flow problem without any porous layer is compared with available data. Second, the porous Couette flow and partially porous channel flow are successfully compared with the studies of other researchers. Then, fluid flow in a clear U-turn duct is studied looking carefully at the velocity, curvature and rotation effects. Finally, the effects of porous layer thickness on the rate of heat transfer and pressure drop are investigated. Parametric studies are conducted to evaluate the effects of various parameters (i.e., Reynolds number, Darcy number, rotation number), highlighting their influences on the thermo-hydrodynamics behavior of the flow. The optimum values of porous layer thickness are presented for specific flow parameters.     

Measurement and Correlation of Pressure Drop Characteristics for Air Flow Through Sintered Metal Porous Media

Sintered metal porous media are currently used to replace conventional orifices as restrictors in air-bearing systems. The flow properties in porous media are generally approximated by Darcy and Forchheimer regimes in different flow regions. In this study, an ISO expanded expression is proven defective when it is used to represent flow properties through porous media under slight pressure drops ( ${<}10$  kPa). A modified Forchheimer equation is therefore developed to correlate the pressure drop with flow rate, including compressibility and inertial effects. Experimental and theoretical investigations on pressure drop characteristics are conducted with a series of metal-sintered porous media. Permeability is first determined in a strict Darcy region with $Re<0.1$ , followed by the inertial coefficient with $Re>0.1$ , rather than determining these two simultaneously. The theoretical mass flow rate in terms of the modified Forchheimer equation provides close approximations to the experimental data.     

Tomography-Based Characterization and Optimization Of Fluid Flow Through Porous Media

Numerical simulations to characterize fluid flow through porous media have been carried out using tomography-derived real geometry data that has been manipulated using digital image processing techniques to obtain a wide range of porosities. Two kinds of porous media have been analyzed: (a) a reticulated porous ceramic (RPC) foam and (b) a packed bed of CaCO₃ particles. The porosity of the media is varied via morphological operations between 0.727 and 0.913 in case of the RPC and between 0.329 and 0.824 in case of the packed bed. A mesh generator based on the pore space indicator function is then used to generate unstructured tetrahedral grids from the processed tomography data. Fluid flow simulations are carried out for Reynolds numbers ranging from 0.1 to 200 and the results are used to determine the permeability and the Dupuit?CForchheimer coefficient in each case. The results are then compared with existing analytical models and the applicability of the models is examined. In the RPC case, the Happel?CBrenner (parallel-flow) model predicts the permeability with a normalized root mean square error (NRMSE) of 11.8 % across the porosity range and Modified Ergun (Macdonald et. al) model predicts the Dupuit?CForchheimer coefficient within a NRMSE of 13.5 %. In the packed-bed case, the Brinkman drag model predicts the permeability within a NRMSE of 8.26 % across the porosity range and the Modified Ergun model predicts the Dupuit?CForchheimer coefficient within an NRMSE of 5.94 %. For each material, an adjusted Kozeny constant is determined. For the RPC, the Kozeny constant is evaluated at 7.73 and for the CaCO₃ packed bed, it is found to be 6.10, leading to predictions of the permeability with an NRMSE of 4.16 and 3.37 %, respectively.     

Vortex structure of turbulence over permeable walls

In order to understand the effect of the wall permeability on the turbulent vortex structure near porous walls, based on PIV experimental data, a probability density analysis of fluctuating velocities, statistical quadrant and quadrant-hole analyses of the Reynolds shear stress are performed. The investigated flow fields are turbulent channel flows whose bottom walls are made of porous media. The porous media used are three kinds of foamed ceramics which have almost the same porosity (∼0.8) but different permeability. From the discussions on those analyses, a conceptual scenario of the development of the vortex structure near a permeable wall is proposed for a moderate permeability Reynolds number case. It explains the reason why the near-wall long streaky structure tends to vanish near a porous wall with increasing wall permeability.     

Predicting the Onset of Inertial Effects in Sandstone Rocks

This study presents a method to determine the onset of inertial effects at the microscopic level, to distinguish between Darcy and non-Darcy flow regions within porous media at the pore level, and to quantify the effects of retained polymer on gas mobility. Capillary pressure and polymer flood experiments were conducted using Elgin and Okesa sandstone samples. The pore-size distributions were used to study the high-velocity flow effects. A modified capillary-orifice model was used to determine the non-Darcy flow effects at the pore level, with and without residual polymer.The overall flow behavior at any flow rate may be described as the average of all contributions from the Darcy and the non-Darcy terms in all pores. Results of this study suggest that the conventional Reynolds number may lead to incorrect analysis of flow behavior when evaluating non-Darcy flow effects in porous media. The Forchheimer number, defined as the ratio of inertial forces to viscous forces, is found more adequate for analyzing microscopic flow behavior in porous media.