Fluid flow through porous media

Last updated

In fluid mechanics, fluid flow through porous media is the manner in which fluids behave when flowing through a porous medium, [1] for example sponge or wood, or when filtering water using sand or another porous material. As commonly observed, some fluid flows through the media while some mass of the fluid is stored in the pores present in the media.

Contents

Classical flow mechanics in porous media assumes that the medium is homogenous, isotropic, and has an intergranular pore structure. It also assumes that the fluid is a Newtonian fluid, that the reservoir is isothermal, that the well is vertical, etc. Traditional flow issues in porous media often involve single-phase steady state flow, multi-well interference, oil-water two-phase flow, natural gas flow, elastic energy driven flow, oil-gas two-phase flow, and gas-water two-phase flow. [2]

The physicochemical flow process will involve various physical property changes and chemical reactions in contrast to the basic Newtonian fluid in the classical flow theory of porous system. Viscosity, surface tension, phase state, concentration, temperature, and other physical characteristics are examples of these properties. Non-Newtonian fluid flow, mass transfer through diffusion, and multiphase and multicomponent fluid flow are the primary flow issues. [3]

Governing laws

The movement of a fluid through porous media is described by the combination of Darcy's law with the principle of conservation of mass in order to express the capillary force or fluid velocity as a function of various other parameters including the effective pore radius, liquid viscosity or permeability. [4] However, the use of Darcy's law alone does not produce accurate results for heterogeneous media like shale, and tight sandstones, where there is a huge proportion of nanopores. This necessitates the use of a flow model that considers the weighted proportion of various flow regimes like Darcy flow, transition flow, slip flow, and free molecular flow. [5]

Darcy's law

SymbolDescription
Volumetric flow rate [m3/s]
Permeability of porous medium [m2]. The permeability is a function of material type, and also varies with stress, temperature, etc.
Fluid viscosity [Pa.s]
Cross-sectional area of Porous medium [m2]
Pressure drop across medium [Pa]
Length of sample [m]

The basic law governing the flow of fluids through porous media is Darcy's Law, which was formulated by the French civil engineer Henry Darcy in 1856 on the basis of his experiments on vertical water filtration through sand beds. [6]

According to Darcy's law, the fluid's viscosity, effective fluid permeability, and fluid pressure gradient determine the flow rate at any given location in the reservoir. [7]

For transient processes in which the flux varies from point to-point, the following differential form of Darcy’s law is used.

Darcy's law is valid for situation where the porous material is already saturated with the fluid. For the calculation of capillary imbibition speed of a liquid to an initially dry medium, Washburn's or Bosanquet's equations are used.

Mass conservation

Mass conservation of fluid across the porous medium involves the basic principle that mass flux in minus mass flux out equals the increase in amount stored by a medium. [8] This means that total mass of the fluid is always conserved. In mathematical form, considering a time period from to , length of porous medium from to and being the mass stored by the medium, we have

Furthermore, we have that , where is the pore volume of the medium between and and is the density. So where is the porosity. Dividing both sides by , while , we have for 1 dimensional linear flow in a porous medium the relation

In three dimensions, the equation can be written as

The mathematical operation on the left-hand side of this equation is known as the divergence of , and represents the rate at which fluid diverges from a given region, per unit volume.

Diffusion Equation

Material TypeCompressibility (m2N−1 or Pa−1) [9]
Clay10−6 - 10−8
Sand10−7 - 10−9
Gravel10−8 - 10−10
Jointed rock10−8 - 10−10
Sound Rock10−9 - 10−11
Water (beta)4.4 x 10−10

Using product rule(and chain rule) on right hand side of the above mass conservation equation (i),

Here, = compressibility of the fluid and = compressibility of porous medium. [10] Now considering the left hand side of the mass conservation equation, which is given by Darcy's Law as

Equating the results obtained in & , we get

The second term on the left side is usually negligible, and we obtain the diffusion equation in 1 dimension as

where . [11]

Related Research Articles

In fluid mechanics, the Rayleigh number (Ra, after Lord Rayleigh) for a fluid is a dimensionless number associated with buoyancy-driven flow, also known as free (or natural) convection. It characterises the fluid's flow regime: a value in a certain lower range denotes laminar flow; a value in a higher range, turbulent flow. Below a certain critical value, there is no fluid motion and heat transfer is by conduction rather than convection. For most engineering purposes, the Rayleigh number is large, somewhere around 106 to 108.

<span class="mw-page-title-main">Stress–energy tensor</span> Tensor describing energy momentum density in spacetime

The stress–energy tensor, sometimes called the stress–energy–momentum tensor or the energy–momentum tensor, is a tensor physical quantity that describes the density and flux of energy and momentum in spacetime, generalizing the stress tensor of Newtonian physics. It is an attribute of matter, radiation, and non-gravitational force fields. This density and flux of energy and momentum are the sources of the gravitational field in the Einstein field equations of general relativity, just as mass density is the source of such a field in Newtonian gravity.

In the calculus of variations, a field of mathematical analysis, the functional derivative relates a change in a functional to a change in a function on which the functional depends.

The primitive equations are a set of nonlinear partial differential equations that are used to approximate global atmospheric flow and are used in most atmospheric models. They consist of three main sets of balance equations:

  1. A continuity equation: Representing the conservation of mass.
  2. Conservation of momentum: Consisting of a form of the Navier–Stokes equations that describe hydrodynamical flow on the surface of a sphere under the assumption that vertical motion is much smaller than horizontal motion (hydrostasis) and that the fluid layer depth is small compared to the radius of the sphere
  3. A thermal energy equation: Relating the overall temperature of the system to heat sources and sinks
<span class="mw-page-title-main">Einstein field equations</span> Field-equations in general relativity

In the general theory of relativity, the Einstein field equations relate the geometry of spacetime to the distribution of matter within it.

Permeability in fluid mechanics and the Earth sciences is a measure of the ability of a porous material to allow fluids to pass through it.

Darcy's law is an equation that describes the flow of a fluid through a porous medium. The law was formulated by Henry Darcy based on results of experiments on the flow of water through beds of sand, forming the basis of hydrogeology, a branch of earth sciences. It is analogous to Ohm's law in electrostatics, linearly relating the volume flow rate of the fluid to the hydraulic head difference via the hydraulic conductivity.

<span class="mw-page-title-main">Electromagnetic tensor</span> Mathematical object that describes the electromagnetic field in spacetime

In electromagnetism, the electromagnetic tensor or electromagnetic field tensor is a mathematical object that describes the electromagnetic field in spacetime. The field tensor was first used after the four-dimensional tensor formulation of special relativity was introduced by Hermann Minkowski. The tensor allows related physical laws to be written very concisely, and allows for the quantization of the electromagnetic field by Lagrangian formulation described below.

In physics and fluid mechanics, a Blasius boundary layer describes the steady two-dimensional laminar boundary layer that forms on a semi-infinite plate which is held parallel to a constant unidirectional flow. Falkner and Skan later generalized Blasius' solution to wedge flow, i.e. flows in which the plate is not parallel to the flow.

<span class="mw-page-title-main">Newman–Penrose formalism</span> Notation in general relativity

The Newman–Penrose (NP) formalism is a set of notation developed by Ezra T. Newman and Roger Penrose for general relativity (GR). Their notation is an effort to treat general relativity in terms of spinor notation, which introduces complex forms of the usual variables used in GR. The NP formalism is itself a special case of the tetrad formalism, where the tensors of the theory are projected onto a complete vector basis at each point in spacetime. Usually this vector basis is chosen to reflect some symmetry of the spacetime, leading to simplified expressions for physical observables. In the case of the NP formalism, the vector basis chosen is a null tetrad: a set of four null vectors—two real, and a complex-conjugate pair. The two real members often asymptotically point radially inward and radially outward, and the formalism is well adapted to treatment of the propagation of radiation in curved spacetime. The Weyl scalars, derived from the Weyl tensor, are often used. In particular, it can be shown that one of these scalars— in the appropriate frame—encodes the outgoing gravitational radiation of an asymptotically flat system.

In fluid dynamics, the Buckley–Leverett equation is a conservation equation used to model two-phase flow in porous media. The Buckley–Leverett equation or the Buckley–Leverett displacement describes an immiscible displacement process, such as the displacement of oil by water, in a one-dimensional or quasi-one-dimensional reservoir. This equation can be derived from the mass conservation equations of two-phase flow, under the assumptions listed below.

The intent of this article is to highlight the important points of the derivation of the Navier–Stokes equations as well as its application and formulation for different families of fluids.

The Kozeny–Carman equation is a relation used in the field of fluid dynamics to calculate the pressure drop of a fluid flowing through a packed bed of solids. It is named after Josef Kozeny and Philip C. Carman. The equation is only valid for creeping flow, i.e. in the slowest limit of laminar flow. The equation was derived by Kozeny (1927) and Carman from a starting point of (a) modelling fluid flow in a packed bed as laminar fluid flow in a collection of curving passages/tubes crossing the packed bed and (b) Poiseuille's law describing laminar fluid flow in straight, circular section pipes.

In nonideal fluid dynamics, the Hagen–Poiseuille equation, also known as the Hagen–Poiseuille law, Poiseuille law or Poiseuille equation, is a physical law that gives the pressure drop in an incompressible and Newtonian fluid in laminar flow flowing through a long cylindrical pipe of constant cross section. It can be successfully applied to air flow in lung alveoli, or the flow through a drinking straw or through a hypodermic needle. It was experimentally derived independently by Jean Léonard Marie Poiseuille in 1838 and Gotthilf Heinrich Ludwig Hagen, and published by Poiseuille in 1840–41 and 1846. The theoretical justification of the Poiseuille law was given by George Stokes in 1845.

<span class="mw-page-title-main">Diffusion</span> Transport of dissolved species from the highest to the lowest concentration region

Diffusion is the net movement of anything generally from a region of higher concentration to a region of lower concentration. Diffusion is driven by a gradient in Gibbs free energy or chemical potential. In the context of Quantum Physics, diffusion refers to spreading of wave packets. In simplest example, a Gaussian wave packet will spread along the spatial dimensions, as time progresses, resulting in diffusion of the wave packet energy. It is possible to diffuse "uphill" from a region of lower concentration to a region of higher concentration, like in spinodal decomposition. Diffusion is a stochastic process due to the inherent randomness of the diffusing entity and can be used to model many real-life stochastic scenarios. Therefore, diffusion and the corresponding mathematical models are used in several fields, beyond physics, such as statistics, probability theory, information theory, neural networks, finance and marketing etc.

Nuclear magnetic resonance (NMR) in porous materials covers the application of using NMR as a tool to study the structure of porous media and various processes occurring in them. This technique allows the determination of characteristics such as the porosity and pore size distribution, the permeability, the water saturation, the wettability, etc.

Lagrangian field theory is a formalism in classical field theory. It is the field-theoretic analogue of Lagrangian mechanics. Lagrangian mechanics is used to analyze the motion of a system of discrete particles each with a finite number of degrees of freedom. Lagrangian field theory applies to continua and fields, which have an infinite number of degrees of freedom.

<span class="mw-page-title-main">Finite volume method for one-dimensional steady state diffusion</span> Scientific Technique

The Finite volume method in computational fluid dynamics is a discretization technique for partial differential equations that arise from physical conservation laws. These equations can be different in nature, e.g. elliptic, parabolic, or hyperbolic. The first well-documented use of this method was by Evans and Harlow (1957) at Los Alamos. The general equation for steady diffusion can be easily be derived from the general transport equation for property Φ by deleting transient and convective terms.

Morris Muskat et al. developed the governing equations for multiphase flow in porous media as a generalisation of Darcy's equation for water flow in porous media. The porous media are usually sedimentary rocks such as clastic rocks or carbonate rocks.

The porous medium equation, also called the nonlinear heat equation, is a nonlinear partial differential equation taking the form:

References

  1. Mohammadizadeh, SeyedMehdi; Moghaddam, Mehdi Azhdary; Talebbeydokhti, Naser (2021). "Analysis of Flow in Porous Media using Combined Pressurized-Free surface Network". Journal of Porous Media. Begel House Inc. 24 (10): 1–15. doi:10.1615/JPorMedia.2021025407. S2CID   235877042.
  2. "Classical flow theory in porous media", Fluid Flow in Porous Media, WORLD SCIENTIFIC, pp. 69–197, 2020-02-28, doi:10.1142/9789811219535_0003, ISBN   978-981-12-1952-8, S2CID   242538024 , retrieved 2022-07-01
  3. "Physicochemical flow in porous media", Fluid Flow in Porous Media, WORLD SCIENTIFIC, pp. 249–285, 2020-02-28, doi:10.1142/9789811219535_0005, ISBN   978-981-12-1952-8, S2CID   241660780 , retrieved 2022-07-01
  4. Tailoring porous media for controllable capillary flow Journal of colloid and interface science (2019), 539, 379-387
  5. Saurabh, Suman (2022-02-15). "Flow regime evolution and stress-dependent permeability in nanoporous rocks". Fuel. Elsevier. 310: 122413. doi:10.1016/j.fuel.2021.122413. S2CID   243474407.
  6. Whitaker, Stephen (1986). "Flow in porous media I: A theoretical derivation of Darcys law". Transport in Porous Media. 1: 3–25. doi:10.1007/BF01036523. S2CID   121904058.
  7. Johns, Russell T. (1 January 2004). "Oil Recovery". Encyclopedia of Energy. Elsevier: 701–713. doi:10.1016/b0-12-176480-x/00255-2. ISBN   9780121764807.
  8. Bear, Jacob (2013-02-26). Dynamics of fluids in Porous Media. ISBN   9780486131801.
  9. https://eng.ucmerced.edu/people/jfisher/.../EnveEss110_20081110.pdf%5B%5D
  10. Ohirhian, Peter. "Steady State Compressible Fluid Flow in Porous Media".{{cite journal}}: Cite journal requires |journal= (help)
  11. Zimmerman, Dr. R.W. "Flow in porous media".

Further reading