scholarly journals On the apparent permeability of porous media in rarefied gas flows

2017 ◽  
Vol 822 ◽  
pp. 398-417 ◽  
Author(s):  
Lei Wu ◽  
Minh Tuan Ho ◽  
Lefki Germanou ◽  
Xiao-Jun Gu ◽  
Chang Liu ◽  
...  

The apparent gas permeability of a porous medium is an important parameter in the prediction of unconventional gas production, which was first investigated systematically by Klinkenberg in 1941 and found to increase with the reciprocal mean gas pressure (or equivalently, the Knudsen number). Although the underlying rarefaction effects are well known, the reason that the correction factor in Klinkenberg’s famous equation decreases when the Knudsen number increases has not been fully understood. Most of the studies idealize the porous medium as a bundle of straight cylindrical tubes; however, according to the gas kinetic theory, this only results in an increase of the correction factor with the Knudsen number, which clearly contradicts Klinkenberg’s experimental observations. Here, by solving the Bhatnagar–Gross–Krook equation in simplified (but not simple) porous media, we identify, for the first time, two key factors that can explain Klinkenberg’s experimental results: the tortuous flow path and the non-unitary tangential momentum accommodation coefficient for the gas–surface interaction. Moreover, we find that Klinkenberg’s results can only be observed when the ratio between the apparent and intrinsic permeabilities is ${\lesssim}30$; at large ratios (or Knudsen numbers) the correction factor increases with the Knudsen number. Our numerical results could also serve as benchmarking cases to assess the accuracy of macroscopic models and/or numerical schemes for the modelling/simulation of rarefied gas flows in complex geometries over a wide range of gas rarefaction. Specifically, we point out that the Navier–Stokes equations with the first-order velocity-slip boundary condition are often misused to predict the apparent gas permeability of the porous medium; that is, any nonlinear dependence of the apparent gas permeability with the Knudsen number, predicted from the Navier–Stokes equations, is not reliable. Worse still, for some types of gas–surface interactions, even the ‘filtered’ linear dependence of the apparent gas permeability with the Knudsen number is of no practical use since, compared to the numerical solution of the Bhatnagar–Gross–Krook equation, it is only accurate when the ratio between the apparent and intrinsic permeabilities is ${\lesssim}1.5$.

Author(s):  
Angelos G. Klothakis ◽  
Georgios N. Lygidakis ◽  
Ioannis K. Nikolos

During the past decade considerable efforts have been exerted for the simulation of rarefied gas flows in a wide range of applications, like the flow over suborbital vehicles, in microelectromechanical systems, etc. Such flows appear to be significantly different from those at the continuum regime, making the Navier-Stokes equations to fail without further amendment. In this study an in-house academic CFD solver, named Galatea, is modified appropriately to account for rarefied gases. The no-slip condition on solid walls is no longer valid, hence, velocity slip and temperature jump boundary conditions are applied instead. Additionally, a second-order accurate slip model has been incorporated, namely, this of Beskok and Karniadakis, increasing the accuracy in the same area but avoiding simultaneously the numerical difficulties, entailed by the computation of the second derivative of slip velocity when complex geometries and unstructured grids are coupled. The proposed solver is validated against rarefied laminar flow over a suborbital shuttle, designed by the Azim’UTBM team. The obtained results are compared with those extracted with the parallel open-source kernel SPARTA, which is based on the DSMC method. A satisfactory agreement is reported between the two methodologies, demonstrating the potential of the modified solver to simulate effectively such flows.


2013 ◽  
Vol 135 (10) ◽  
Author(s):  
Y. S. Muzychka ◽  
R. Enright

Analytical solutions for slip flows in the hydrodynamic entrance region of tubes and channels are examined. These solutions employ a linearized axial momentum equation using Targ's method. The momentum equation is subjected to a first order Navier slip boundary condition. The accuracy of these solutions is examined using computational fluid dynamics (CFD) simulations. CFD simulations utilized the full Navier–Stokes equations, so that the implications of the approximate linearized axial momentum equation could be fully assessed. Results are presented in terms of the dimensionless mean wall shear stress, τ⋆, as a function of local dimensionless axial coordinate, ξ, and relative slip parameter, β. These solutions can be applied to either rarefied gas flows when compressibility effects are small or apparent liquid slip over hydrophobic and superhydrophobic surfaces. It has been found that, under slip conditions, the minimum Reynolds number should be ReDh>100 in order for the approximate linearized solution to remain valid.


2016 ◽  
Vol 792 ◽  
pp. 5-35 ◽  
Author(s):  
Giuseppe A. Zampogna ◽  
Alessandro Bottaro

The interaction between a fluid flow and a transversely isotropic porous medium is described. A homogenized model is used to treat the flow field in the porous region, and different interface conditions, needed to match solutions at the boundary between the pure fluid and the porous regions, are evaluated. Two problems in different flow regimes (laminar and turbulent) are considered to validate the system, which includes inertia in the leading-order equations for the permeability tensor through a Oseen approximation. The components of the permeability, which characterize microscopically the porous medium and determine the flow field at the macroscopic scale, are reasonably well estimated by the theory, both in the laminar and the turbulent case. This is demonstrated by comparing the model’s results to both experimental measurements and direct numerical simulations of the Navier–Stokes equations which resolve the flow also through the pores of the medium.


2020 ◽  
Vol 8 (2) ◽  
pp. 87 ◽  
Author(s):  
Paran Pourteimouri ◽  
Kourosh Hejazi

An integrated two-dimensional vertical (2DV) model was developed to investigate wave interactions with permeable submerged breakwaters. The integrated model is capable of predicting the flow field in both surface water and porous media on the basis of the extended volume-averaged Reynolds-averaged Navier–Stokes equations (VARANS). The impact of porous medium was considered by the inclusion of the additional terms of drag and inertia forces into conventional Navier–Stokes equations. Finite volume method (FVM) in an arbitrary Lagrangian–Eulerian (ALE) formulation was adopted for discretization of the governing equations. Projection method was utilized to solve the unsteady incompressible extended Navier–Stokes equations. The time-dependent volume and surface porosities were calculated at each time step using the fraction of a grid open to water and the total porosity of porous medium. The numerical model was first verified against analytical solutions of small amplitude progressive Stokes wave and solitary wave propagation in the absence of a bottom-mounted barrier. Comparisons showed pleasing agreements between the numerical predictions and analytical solutions. The model was then further validated by comparing the numerical model results with the experimental measurements of wave propagation over a permeable submerged breakwater reported in the literature. Good agreements were obtained for the free surface elevations at various spatial and temporal scales, velocity fields around and inside the obstacle, as well as the velocity profiles.


2011 ◽  
Vol 2011 ◽  
pp. 1-10 ◽  
Author(s):  
S. Islam ◽  
Hamid Khan ◽  
Inayat Ali Shah ◽  
Gul Zaman

The flow between two large parallel plates approaching each other symmetrically in a porous medium is studied. The Navier-Stokes equations have been transformed into an ordinary nonlinear differential equation using a transformationψ(r,z)=r2F(z). Solution to the problem is obtained by using differential transform method (DTM) by varying different Newtonian fluid parameters and permeability of the porous medium. Result for the stream function is presented. Validity of the solutions is confirmed by evaluating the residual in each case, and the proposed scheme gives excellent and reliable results. The influence of different parameters on the flow has been discussed and presented through graphs.


2018 ◽  
Vol 17 (01) ◽  
pp. 85-116 ◽  
Author(s):  
Marc Briant ◽  
Sara Merino-Aceituno ◽  
Clément Mouhot

We study the Boltzmann equation on the [Formula: see text]-dimensional torus in a perturbative setting around a global equilibrium under the Navier–Stokes linearization. We use a recent functional analysis breakthrough to prove that the linear part of the equation generates a [Formula: see text]-semigroup with exponential decay in Lebesgue and Sobolev spaces with polynomial weight, independently of the Knudsen number. Finally, we prove well-posedness of the Cauchy problem for the nonlinear Boltzmann equation in perturbative setting and an exponential decay for the perturbed Boltzmann equation, uniformly in the Knudsen number, in Sobolev spaces with polynomial weight. The polynomial weight is almost optimal. Furthermore, this result only requires derivatives in the space variable and allows to connect solutions to the incompressible Navier–Stokes equations in these spaces.


2012 ◽  
Vol 698 ◽  
pp. 406-422 ◽  
Author(s):  
Thomas Veltzke ◽  
Jorg Thöming

AbstractIn microducts deviation from continuum flow behaviour of a gas increases with rarefaction. When using Navier–Stokes equations to calculate a flow under slightly and moderately rarefied conditions, slip boundary conditions are used which in turn refer to the tangential momentum accommodation coefficient (TMAC). Here we demonstrate that, in the so-called slip and transition regime, the flow in microducts can be reliably described by a consistently non-empirical model without considering the TMAC. We obtain this equation by superposition of convective transport and Fickian diffusion using two-dimensional solutions of Navier–Stokes equations and a description for the Knudsen diffusion coefficient as derived from kinetic theory respectively. For a wide variety of measurement series found in the literature the calculation predicts the data accurately. Surprisingly only size of the duct, temperature, gas properties and inlet and outlet pressure are necessary to calculate the resulting mass flow by means of a single algebraic equation. From this, and taking the discrepancies of the TMAC concerning surface roughness and nature of the gases into account, we could conclude that neither the diffusive proportions nor the total mass flow rates are influenced by surface topology and chemistry at Knudsen numbers below unity. Compared to the tube geometry, the model slightly underestimates the flow rate in rectangular channels when rarefaction increases. Likewise, the dimensionless mass flow rate and the diffusive proportion of the total flow are distinctly higher in a tube. Thus the cross-sectional geometry has a significant influence on the transport mechanisms under rarefied conditions.


Sign in / Sign up

Export Citation Format

Share Document