Flow Regimes in Packed Beds of Spheres from Pre-Darcy to Turbulent

Characteristics of flow regimes in porous media, along the processes of energy dissipation in each regime, are critical for applications of such media. The current work presents new experimental data for water flow in packed steel spheres of 1- and 3-mm diameters. The porosity of the porous media was about 35 % for both cases. The extensive dataset covered a broad range of flow Reynolds number such that several important flow regimes were encountered, including the elusive pre-Darcy regime, which is rarely or never seen in porous-media literature and turbulent regime. When compared to previous information, the results of this study are seen to add to the divergence of available data on pressure drop in packed beds of spheres. The divergence was also present in the coefficients of Ergun equation and in the Kozeny–Carman constant. The porous media of the current work were seen to exhibit different values of permeability and Forchheimer coefficient in each flow regime. The current data correlated well using the friction factor based on the permeability (measured in the Darcy regime) and the Reynolds number based on the same length scale. An attempt was made to apply recent theoretical results regarding the applicability of the quadratic and cubic Forchheimer corrections in the strong and weak inertia regime.     

Letter to the Editor: Comments on the Paper “Inertia Effects in High-Rate Flow Through Heterogeneous Porous Media” by M. Fourar, R. Lenormand, M. Karimi-Fard, and R. Horne, Transport in Porous Media, DOI 10.1007/s11242-004-6800-6, 2005

In a recent article, Fourar et al. (Transp Porous Med, 2005, doi:10.1007/s11242-004-6800-6) analyzed the effect of heterogeneity in the permeability distribution on Forchheimer flow in porous media. They derived expressions to calculate the effective inertial coefficient in serial layers, parallel layers, and two-dimensional correlated media. Here, we highlight an inconsistency in their first-order expression for serial layers and extend their findings by providing closed-form expressions for the effective inertial coefficient in the case of a lognormal permeability distribution.     

Numerical Study of Heat Transfer in a Rectangular Channel with Porous Covering Obstacles

A numerical study is performed to analyze steady laminar forced convection in a channel in which discrete heat sources covered with porous material are placed on the bottom wall. Hydrodynamic and heat transfer results are reported. The flow in the porous medium is modeled using the Darcy–Brinkman–Forchheimer model. A computer program based on control volume method with appropriate averaging for diffusion coefficient is developed to solve the coupling between solid, fluid, and porous region. The effects of parameters such as Reynolds number, Prandtl number, inertia coefficient, and thermal conductivity ratio are considered. The results reveal that the porous cover with high thermal conductivity enhances the heat transfer from the solid blocks significantly and decreases the maximum temperature on the heated solid blocks. The mean Nusselt number increases with increase of Reynolds number and Prandtl number, and decrease of inertia coefficient. The pressure drop along the channel increases rapidly with the increase of Reynolds number.     

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.     

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.     

Numerical Simulation of Two-Phase Inertial Flow in Heterogeneous Porous Media

In this study, non-Darcy inertial two-phase incompressible and non-stationary flow in heterogeneous porous media is analyzed using numerical simulations. For the purpose, a 3D numerical tool was fully developed using a finite volume formulation, although for clarity, results are presented in 1D and 2D configurations only. Since a formalized theoretical model confirmed by experimental data is still lacking, our study is based on the widely used generalized Darcy–Forchheimer model. First, a validation is performed by comparing numerical results of the saturation front kinetics with a semi-analytical solution inspired from the Buckley–Leverett model extended to take into account inertia. Second, we highlight the importance of inertial terms on the evolution of saturation fronts as a function of a suitable Reynolds number. Saturation fields are shown to have a structure markedly different from the classical case without inertia, especially for heterogeneous media, thereby, emphasizing the necessity of a more complete model than the classical generalized Darcy’s one when inertial effects are not negligible.     

Inertial and Darcy’s Terms Ratio in Boundary Layer at Fluid–Porous Medium Interface

The study considers the forced boundary-layer flow overlying the Darcy–Brinkman porous medium and gives a quantitative analysis of the nonlinear inertial terms in the Brinkman filtration equation. The inertial terms are shown to be larger than the Darcy’s drag near the porous medium interface. The applicability range of boundary-layer approach is determined. It is suitable in high-permeable media with moderate velocities of an external flow. If it is slow enough, the inertial terms can be omitted in spite of interface effect. On the other hand, fast external flow produces the filtration with large pore-scale Reynolds number; therefore, the Forchheimer’s drag should be taken into account. It is shown the Brinkman term as well as inertial terms have a significant role in boundary-layer formation within the porous medium.     

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.     

Drag difference between steady and shocked gas flows passing through a porous body

Experimental and numerical investigations of gas flows through porous materials have been carried out. We have investigated steady and unsteady processes occurring when the gas flow interacts with porous materials. Densities and porosities of the four open-cell-type polyurethane foams which were investigated are kg/m and , with the foams having different structures. Experiments were conducted to determine the steady drag coefficient of the porous material at low Reynolds numbers, evaluated from the pressure drop. The Forchheimer equation was applied to determine the drag. Values of permeability coefficients () in the Forchheimer equation were estimated by comparing computed and experimental results. Results show that the drag coefficient is largely affected by the internal structure of the foam, and has almost no effect on the stress history, while the value of dominates the stress history variation. Differences of 1000 times exist between the steady flow and unsteady shock tube flow values. Received 15 May 1998/ Accepted 15 March 1999     

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.     

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.     

考虑气体压缩性的多孔材料渗透率和惯性系数的测定

多孔介质材料的渗透率和惯性系数是决定多孔介质中流体流动特性的重要参数，一般需要通过实验进行测定，在测定渗透律和惯性系数量，选用气体作为工作介质可以为实验带来极大的方便，然而通常的实验都将气体看作不可压缩流体，直接根据Darcy-Forchheimer定律得到这两个参数，这种近似对实验条件如样品厚度、工作压力等提出了很多要求，本文提出了在考虑气体压缩性的情况下测定渗透率和惯性系数的方法，该方法可以大大降低实验时对样品厚度、工作压力等条件的要求。本文还根据该方法对多孔材料PVF进行了渗透率和惯性系数的测定，并对测量结果进行了验证。     

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.     

Pulsatile magneto-biofluid flow and mass transfer in a non-Darcian porous medium channel

A numerical study of pulsatile flow and mass transfer of an electrically conducting Newtonian biofluid via a channel containing porous medium is considered. The conservation equations are transformed and solved under boundary conditions prescribed at both walls of the channel, using a finite element method with two-noded line elements. The influence of magnetic field on the flow is studied using the dimensionless hydromagnetic number, Nm, which defines the ratio of magnetic (Lorentz) retarding force to the viscous hydrodynamic force. A Darcian linear impedance for low Reynolds numbers is incorporated in the transformed momentum equation and a second order drag force term for inertial (Forchheimer) effects. Velocity and concentration profiles across the channel width are plotted for various values of the Reynolds number (Re), Darcy parameter (λ), Forchheimer parameter (Nf), hydro-magnetic number (Nm), Schmidt number (Sc) and also with dimensionless time (T). Profiles of velocity varying in space and time are also provided. The conduit considered is rigid with a pulsatile pressure applied via an appropriate pressure gradient term. Increasing the hydromagnetic number (Nm) from 1 to 15 considerably depresses biofluid velocity (U) indicating that a magnetic field can be used as a flow control mechanism in, for example, medical applications. A rise in Nf from 1 to 20 strongly retards the flow development and decreases the velocity, U, across the width of the channel. The effects of other parameters on the flowfield are also discussed at length. The flow model also has applications in the analysis of electrically conducting haemotological fluids flowing through filtration media, diffusion of drug species in pharmaceutical hydromechanics, and also in general fluid dynamics of pulsatile systems.     

The Forchheimer equation: A theoretical development

In this paper we illustrate how the method of volume averaging can be used to derive Darcy's law with the Forchheimer correction for homogeneous porous media. Beginning with the Navier-Stokes equations, we find the volume averaged momentum equation to be given by $$\langle v_\beta \rangle = - \frac{K}{{\mu _\beta }} \cdot (\nabla \langle p_\beta \rangle ^\beta - \rho _\beta g) - F\cdot \langle v_\beta \rangle .$$ The Darcy's law permeability tensor, K, and the Forchheimer correction tensor, F, are determined by closure problems that must be solved using a spatially periodic model of a porous medium. When the Reynolds number is small compared to one, the closure problem can be used to prove that F is a linear function of the velocity, and order of magnitude analysis suggests that this linear dependence may persist for a wide range of Reynolds numbers.     

The Use of Porous Fins for Heat Transfer Augmentation in Parallel-Plate Channels

In this study, a steady, fully developed laminar forced convection heat augmentation via porous fins in isothermal parallel-plate duct is numerically investigated. High-thermal conductivity porous fins are attached to the inner walls of two parallel-plate channels to enhance the heat transfer characteristics of the flow under consideration. The Darcy–Brinkman–Forchheimer model is used to model the flow inside the porous fins. This study reports the effect of several operating parameters on the flow hydrodynamics and thermal characteristics. This study demonstrates, mainly, the effects of porous fin thickness, Darcy number, thermal conductivity ratio, Reynolds number, and microscopic inertial coefficient on the thermal performance of the present flow. It is found that the highest Nusselt number is achieved at fully filled porous duct which requires the highest pumping pressure. The results show that using porous fins requires less pumping pressure with comparable high heat augmentation weight against fully filled porous duct. It is found that higher Nusselt numbers are achieved by increasing the microscopic inertial coefficient (A), the Reynolds number (Re), and the thermal conductivity of the porous substrate k ₂. The results show that heat transfer can be enhanced (1) with the use of high thermal conductivity fins, (2) by decreasing the Darcy number, and (3) by increasing microscopic inertial coefficient.     

Wavelet upscaling based on piecewise bilinear approximation of the permeability field

A method for upscaling of permeability in heterogeneous porous media is presented. The upscaled field takes the form K = e ^Y , where Y, in two dimensions, is a piecewise bilinear function. The method is tested on a number of random permeability fields, with different integral scale/correlation length and variance. The numerical results show that this method conserves much more of the heterogeneous fingering than classical block-based upscaling methods, e.g., geometric mean.     

岩体三维裂隙网络地质熵与渗透特性关系研究

裂隙网络是岩体地下水的主要流动通道,而工程岩体中裂隙网络错综复杂,裂隙网络的几何特征和连通性对其渗透性有着重要影响.为了综合量化裂隙迹长、间距、倾角、开度对裂隙网络连通性和渗透性的影响,基于信息熵原理,提出了三维裂隙网络地质熵理论和连通性指标-熵尺度,对比熵尺度与其他传统三维裂隙网络连通性指标,验证了熵尺度评价三维裂隙网络连通性和渗透性的合理性.结合锦屏一级水电站左岸边坡裂隙统计分布,建立三维裂隙网络渗流数值计算方法,分析不同裂隙迹长、倾角、间距、开度条件下三维裂隙面密度、无量纲逾渗密度、熵尺度和渗透系数的变化关系.结果表明:当体积率一定,考虑开度影响时,三维裂隙面密度和无量纲逾渗密度无法定量表征迹长和间距对裂隙网络连通性的影响;裂隙迹长与熵尺度和渗透系数呈负相关关系,裂隙间距和开度与熵尺度和渗透系数呈正相关关系,裂隙倾角变化对熵尺度和渗透系数影响较小;熵尺度与渗透系数的非线性关系近似满足二次多项式.