Numerical modelling of nonlinear effects in laminar flow through a porous medium

1988 ◽  
Vol 190 ◽  
pp. 393-407 ◽  
Author(s):  
O. Coulaud ◽  
P. Morel ◽  
J. P. Caltagirone

This paper deals with the introduction of a nonlinear term into Darcy's equation to describe inertial effects in a porous medium. The method chosen is the numerical resolution of flow equations at a pore scale. The medium is modelled by cylinders of either equal or unequal diameters arranged in a regular pattern with a square or triangular base. For a given flow through this medium the pressure drop is evaluated numerically.The Navier-Stokes equations are discretized by the mixed finite-element method. The numerical solution is based on operator-splitting methods whose purpose is to separate the difficulties due to the nonlinear operator in the equation of motion and the necessity of taking into account the continuity equation. The associated Stokes problems are solved by a mixed formulation proposed by Glowinski & Pironneau.For Reynolds numbers lower than 1, the relationship between the global pressure gradient and the filtration velocity is linear as predicted by Darcy's law. For higher values of the Reynolds number the pressure drop is influenced by inertial effects which can be interpreted by the addition of a quadratic term in Darcy's law.On the one hand this study confirms the presence of a nonlinear term in the motion equation as experimentally predicted by several authors, and on the other hand analyses the fluid behaviour in simple media. In addition to the detailed numerical solutions, an estimation of the hydrodynamical constants in the Forchheimer equation is given in terms of porosity and the geometrical characteristics of the models studied.

2005 ◽  
Author(s):  
C. A. Ortega Vivas ◽  
S. Barraga´n Gonza´lez ◽  
J. M. Garibay Cisneros

This study analyses the macroscopic flow through a two dimensional porous medium model by numerical and experimental methods. The objective of this research is to develop an empirical model by which the pressure drop can be obtained. In order to construct the model, a series of blocks are used as an idealized pressure drop device, so that the pressure drop can be calculated. The range of porosities studied is between 28 and 75 per cent. It is found that the pressure drop is a combination of viscosity and inertial effects, the later being more important as the Reynolds number is increased. The empirical equation obtained in this investigation is compared with the Ergun equation.


2014 ◽  
Vol 62 (2) ◽  
pp. 321-329
Author(s):  
R. Wojnar

Abstract The thermal effects of a stationary Stokesian flow through an elastic micro-porous medium are compared with the entropy produced by Darcy’s flow. A micro-cellular elastic medium is considered as an approximation of the elastic porous medium. It is shown that after asymptotic two-scale analysis these two approaches, one analytical, starting from Stoke’s equation and the second phenomenological, starting from Darcy’s law give the same result. The incompressible and linearly compressible fluids are considered, and it is shown that in micro-porous systems the seepage of both types of fluids is described by the same equations.


2017 ◽  
Vol 4 (20) ◽  
pp. 359-374
Author(s):  
Ryszard Wojnar

Brinkman’s law is describing the seepage of viscous fluid through a porous medium and is more acurate than the classical Darcy’s law. Namely, Brinkman’s law permits to conform the flow through a porous medium to the free Stokes’ flow. However, Brinkman’s law, similarly as Schro¨dinger’s equation was only devined. Fluid in its motion through a porous solid is interacting at every point with the walls of pores, but the interactions of the fluid particles inside pores are different than the interactions at the walls, and are described by Stokes’ equation. Here, we arrive at Brinkman’s law from Stokes’ flow equation making use of successive iterations, in type of Born’s approximation method, and using Darcy’s law as a zero-th approximation.


1996 ◽  
Vol 06 (08) ◽  
pp. 1051-1077 ◽  
Author(s):  
A. LYAGHFOURI

In this paper we study a fluid flow through a porous medium with Dirichlet boundary conditions and a general permeability. We establish the continuity of the free boundary and the uniqueness of the S3-connected solution.


Author(s):  
B Eitzinger ◽  
G Ederer

AbstractThis study investigates by nonlinear constitutive equations the influence of tipping paper, cigarette paper, filter, and tobacco rod on the degree of filter ventilation and draw resistance. Starting from the laws of conservation, the path to the theory of fluid dynamics in porous media and Darcy's law is reviewed and, as an extension to Darcy's law, two different nonlinear pressure drop-flow relations are proposed. It is proven that these relations are valid constitutive equations and the partial differential equations for the stationary flow in an unlit cigarette covering anisotropic, inhomogeneous and nonlinear behaviour are derived. From these equations a system of ordinary differential equations for the one-dimensional flow in the cigarette is derived by averaging pressure and velocity over the cross section of the cigarette. By further integration, the concept of an electrical analog is reached and discussed in the light of nonlinear pressure drop-flow relations. By numerical calculations based on the system of ordinary differential equations, it is shown that the influence of nonlinearities cannot be neglected because variations in the degree of filter ventilation can reach up to 20% of its nominal value.


Author(s):  
William G. Gray ◽  
Michael A. Celia

The mathematical study of flow in porous media is typically based on the 1856 empirical result of Henri Darcy. This result, known as Darcy’s law, states that the velocity of a single-phase flow through a porous medium is proportional to the hydraulic gradient. The publication of Darcy’s work has been referred to as “the birth of groundwater hydrology as a quantitative science” (Freeze and Cherry, 1979). Although Darcy’s original equation was found to be valid for slow, steady, one-dimensional, single-phase flow through a homogeneous and isotropic sand, it has been applied in the succeeding 140 years to complex transient flows that involve multiple phases in heterogeneous media. To attain this generality, a modification has been made to the original formula, such that the constant of proportionality between flow and hydraulic gradient is allowed to be a spatially varying function of the system properties. The extended version of Darcy’s law is expressed in the following form: qα=-Kα . Jα (2.1) where qα is the volumetric flow rate per unit area vector of the α-phase fluid, Kα is the hydraulic conductivity tensor of the α-phase and is a function of the viscosity and saturation of the α-phase and of the solid matrix, and Jα is the vector hydraulic gradient that drives the flow. The quantities Jα and Kα account for pressure and gravitational effects as well as the interactions that occur between adjacent phases. Although this generalization is occasionally criticized for its shortcomings, equation (2.1) is considered today to be a fundamental principle in analysis of porous media flows (e.g., McWhorter and Sunada, 1977). If, indeed, Darcy’s experimental result is the birth of quantitative hydrology, a need still remains to build quantitative analysis of porous media flow on a strong theoretical foundation. The problem of unsaturated flow of water has been attacked using experimental and theoretical tools since the early part of this century. Sposito (1986) attributes the beginnings of the study of soil water flow as a subdiscipline of physics to the fundamental work of Buckingham (1907), which uses a saturation-dependent hydraulic conductivity and a capillary potential for the hydraulic gradient.


Author(s):  
K. Yazdchi ◽  
S. Srivastava ◽  
S. Luding

Many important natural processes involving flow through porous media are characterized by large filtration velocity. Therefore, it is important to know when the transition from viscous to the inertial flow regime actually occurs in order to obtain accurate models for these processes. In this paper, a detailed computational study of laminar and inertial, incompressible, Newtonian fluid flow across an array of cylinders is presented. Due to the non-linear contribution of inertia to the transport of momentum at the pore scale, we observe a typical departure from Darcy’s law at sufficiently high Reynolds number (Re). Our numerical results show that the weak inertia correction to Darcy’s law is not a square or a cubic term in velocity, as it is in the Forchheimer equation. Best fitted functions for the macroscopic properties of porous media in terms of microstructure and porosity are derived and comparisons are made to the Ergun and Forchheimer relations to examine their relevance in the given porosity and Re range. The results from this study can be used for verification and validation of more advanced models for particle fluid interaction and for the coupling of the discrete element method (DEM) with finite element method (FEM).


1972 ◽  
Vol 54 (1) ◽  
pp. 153-161 ◽  
Author(s):  
Enok Palm ◽  
Jan Erik Weber ◽  
Oddmund Kvernvold

For convection in a porous medium the dependence of the Nusselt number on the Rayleigh number is examined to sixth order using an expansion for the Rayleigh number proposed by Kuo (1961). The results show very good agreement with experiment. Additionally, the abrupt change which is observed in the heat transport at a supercritical Rayleigh number may be explained by a breakdown of Darcy's law.


1998 ◽  
Vol 09 (08) ◽  
pp. 1491-1503 ◽  
Author(s):  
David M. Freed

An extension to the basic lattice-BGK algorithm is presented for modeling a simulation region as a porous medium. The method recovers flow through a resistance field with arbitrary values of the resistance tensor components. Corrections to a previous algorithm are identified. Simple validation tests are performed which verify the accuracy of the method, and demonstrate that inertial effects give a deviation from Darcy's law for nominal simulation velocities.


Sign in / Sign up

Export Citation Format

Share Document