Flow in porous media III: Deformable media

Stokes flow in a deformable medium is considered in terms of an isotropic, linearly elastic solid matrix. The analysis is restricted to steady forms of the momentum equations and small deformation of the solid phase. Darcy's law can be used to determine the motion of the fluid phase; however, the determination of the Darcy's law permeability tensor represents part of the closure problem in which the position of the fluid-solid interface must be determined.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A ^* interfacial area of the- interface contained within a unit cell, m² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - E Young's modulus for the-phase, N/m² - e i unit base vectors (i = 1, 2, 3) - g gravity vector, m²/s - H height of elastic, porous bed, m - k unit base vector (=e ₃) - characteristic length scale for the-phase, m - L characteristic length scale for volume-averaged quantities, m - n unit normal vector pointing from the-phase toward the-phase (n = -n ) - p pressure in the-phase, N/m² - P p – g·r, N/m² - r ₀ radius of the averaging volume, m - r position vector, m - t time, s - T total stress tensor in the-phase, N/m² - T ⁰ hydrostatic stress tensor for the-phase, N/m² - u displacement vector for the-phase, m - V averaging volume, m₃ - V volume of the-phase contained within the averaging volume, m³ - v velocity vector for the-phase, m/s Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - shear coefficient of viscosity for the-phase, Nt/m² - first Lamé coefficient for the-phase, N/m² - second Lamé coefficient for the-phase, N/m² - bulk coefficient of viscosity for the-phase, Nt/m² - T –T ⁰ , a deviatoric stress tensor for the-phase, N/m²     

Flow of Maxwell fluids in porous media

In this paper we analyze the flow of a Maxwell fluid in a rigid porous medium using the method of volume averaging. We first present the local volume averaged momentum equation which contains Darcy-scale elastic effects and undetermined integrals of the spatial deviations of the pressure and velocity. A closure problem is developed in order to determine the spatial deviations and thus obtain a closed form of the momentum equation that contains a time-dependent permeability tensor. To gain some insight into the effects of elasticity on the dynamics of flow in porous media, the entire problem is transformed to the frequency domain through a temporal Fourier transform. This leads to a dynamic generalization of Darcy's law. Analytical results are provided for the case in which the porous medium is modeled as a bundle of capillary tubes, and a scheme is presented to solve the transformed closure problem for a general microstructure.     

Transport in ordered and disordered porous media IV: Computer generated porous media for three-dimensional systems

In the method of volume averaging, the difference between ordered and disordered porous media appears at two distinct points in the analysis, i.e. in the process of spatial smoothing and in the closure problem. In theclosure problem, the use of spatially periodic boundary conditions isconsistent with ordered porous media and the fields under consideration when the length-scale constraint,r ₀L is satisfied. For disordered porous media, spatially periodic boundary conditions are an approximation in need of further study.In theprocess of spatial smoothing, average quantities must be removed from area and volume integrals in order to extractlocal transport equations fromnonlocal equations. This leads to a series of geometrical integrals that need to be evaluated. In Part II we indicated that these integrals were constants for ordered porous media provided that the weighting function used in the averaging process contained thecellular average. We also indicated that these integrals were constrained by certain order of magnitude estimates for disordered porous media. In this paper we verify these characteristics of the geometrical integrals, and we examine their values for pseudo-periodic and uniformly random systems through the use of computer generated porous media.

Nomenclature

Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - a _i i=1, 2, 3 gaussian probability distribution used to locate the position of particles - I unit tensor - L general characteristic length for volume averaged quantities, m - L characteristic length for , m - L characteristic length for , m - characteristic length for the -phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1, 2, 3 lattice vectors, m - m convolution product weighting function - m _v special convolution product weighting function associated with the traditional volume average - n _i i=1, 2, 3 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - n _e outwardly directed unit normal vector at the entrances and exits of the-phase - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - r position vector, m - r _m support of the weighting functionm, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume,, m³ - x positional vector locating the centroid of an averaging volume, m - x ₀ reference position vector associated with the centroid of an averaging volume, m - y position vector locating points relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - /L, small parameter in the method of spatial homogenization - standard deviation ofa _i - _r standard deviation ofr - _r intrinsic phase average of      

Transport in ordered and disordered porous media II: Generalized volume averaging

In this paper we develop the averaged form of the Stokes equations in terms of weighting functions. The analysis clearly indicates at what point one must choose a media-specific weighting function in order to achieve spatially smoothed transport equations. The form of the weighting function that produces the cellular average is derived, and some important geometrical theorems are presented.Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - A _p surface area of a particle, m² - d _p 6V _p/A_p, effective particle diameter, m - g gravity vector, m/s² - I unit tensor - K _m permeability tensor for the weighted average form of Darcy's law, m² - L general characteristic length for volume averaged quantities, m - L _p general characteristic length for volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - l characteristic length (pore scale) for the-phase - l _i i=1, 2, 3 lattice vectors, m - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - m _v special convolution product weighting function associated with the traditional averaging volume - m _g general convolution product weighting function - m _V unit cell convolution product weighting function - m _C special convolution product weighting function for ordered media which produces the cellular average - m _D special convolution product weighting function for disordered media - m _M master convolution product weighting function for ordered and disordered media - n unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p_m superficial weighted average pressure, N/m² - p _m intrinsic weighted average pressure, N/m² - p traditional intrinsic volume averaged pressure, N/m² - p p p _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function, m - r position vector, m - r position vector locating points in the-phase, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - V velocity vector in the-phase, m/s - vm superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - V volume of the-phase contained in the averaging volume, m³ - V _p volume of a particle, m³ - v traditional superficial volume averaged velocity, m/s - v v p _m spatial deviation velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - _m m * . weighted average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m² - V /V, volume fraction of the-phase     

Transport in ordered and disordered porous media V: Geometrical results for two-dimensional systems

In this paper we continue the geometrical studies of computer generated two-phase systems that were presented in Part IV. In order to reduce the computational time associated with the previous three-dimensional studies, the calculations presented in this work are restricted to two dimensions. This allows us to explore more thoroughly the influence of the size of the averaging volume and to learn something about the use of anon-representative region in the determination of averaged quantities.

Nomenclature

Roman Letters A interfacial area of the– interface associated with the local closure problem, m² - a _i i=1, 2, gaussian probability distribution used to locate the position of particles - l unit tensor - characteristic length for the-phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1,2,3 lattice vectors, m - m convolution product weighting function - m _V special convolution product weighting function associated with a unit cell - n _i i=1, 2 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume,V, m³ - x position of the centroid of an averaging area, m - x ₀ reference position of the centroid of an averaging area, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters V /V, volume average porosity - _{a i} standard deviation ofa _i - _r standard deviation ofr - intrinsic phase average of      

Origin and quantification of coupling between relative permeabilities for two-phase flows in porous media

An extended formulation of Darcy's two-phase law is developed on the basis of Stokes' equations. It leads, through results borrowed from the thermodynamics of irreversible processes, to a matrix of relative permeabilities. Nondiagonal coefficients of this matrix are due to the viscous coupling exerted between fluid phases, while diagonal coefficients represent the contribution of both fluid phases to the total flow, as if they were alone. The coefficients of this matrix, contrary to standard relative permeabilities, do not depend on the boundary conditions imposed on two-phase flow in porous media, such as the flow rate. This formalism is validated by comparison with experimental results from tests of two-phase flow in a square cross-section capillary tube and in porous media. Coupling terms of the matrix are found to be nonnegligible compared to diagonal terms. Relationships between standard relative permeabilities and matrix coefficients are studied and lead to an experimental way to determine the new terms for two-phase flow in porous media.     

Transport in ordered and disordered porous media I: The cellular average and the use of weighting functions

In this work we consider transport in ordered and disordered porous media using singlephase flow in rigid porous mediaas an example. We defineorder anddisorder in terms of geometrical integrals that arise naturally in the method of volume averaging, and we show that dependent variables for ordered media must generally be defined in terms of thecellular average. The cellular average can be constructed by means of a weighting function, thus transport processes in both ordered and disordered media can be treated with a single theory based on weighted averages. Part I provides some basic ideas associated with ordered and disordered media, weighted averages, and the theory of distributions. In Part II a generalized averaging procedure is presented and in Part III the closure problem is developed and the theory is compared with experiment. Parts IV and V provide some geometrical results for computer generated porous media.Roman Letters A interfacial area of the- interface contained within the macroscopic region, m² - A_e area of entrances and exits for the-phase contained within the macroscopic system, m² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor, m² - L general characteristic length for volume averaged quantities, m - characteristic length (pore scale) for the-phase - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - N unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p0 reference pressure in the-phase, N/m² - p traditional intrinsic volume averaged pressure, N/m² - r0 radius of a spherical averaging volume, m - r position vector, m - r position vector locating points in the-phase, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - v velocity vector in the-phase, m/s - v traditional superficial volume averaged velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V/V, volume average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m²     

Viscous coupling in two-phase flow in porous media and its effect on relative permeabilities

An idealized model of a porous rock consisting of a bundle of capillary tubes whose cross-sections are regular polygons is used to assess the importance of viscous coupling or lubrication during simultaneous oil-water flow. Fluids are nonuniformly distributed over tubes of different characteristic dimension because of the requirements of capillary equilibrium and the effect of interfacial viscosity at oil-water interfaces is considered. With these assumptions, we find that the importance of viscous coupling depends on the rheology of the oil-water interface. Where the interfacial shear viscosity is zero, viscous coupling leading to a dependence of oil relative permeability on oil-water viscosity ratio for viscosity ratios greater than one is important for a range of pore cross-section shapes and pore size distributions. For nonzero interfacial shear viscosity, viscous coupling is reduced. Using values reported in the literature for crude oil-brine systems, we find no viscous coupling.     

Modeling species transport by concentrated brine in aggregated porous media

Basic equations governing the transport of species by concentrated brine flowing through an aggregated porous medium are developed. Some simple examples are solved numerically. The medium is considered to be composed of porous rock aggregates separated by ‘macropores’ through which the brine flows and transport of salt and low-concentration species takes place. The aggregates contain dead-end pores, cracks, and stationary pockets collectively called ‘micropores’. The micropore space does not contribute to the flow, but it serves as a storage for salt and species. Adsorption of fluid species takes place at internal surface of aggregates where it is assumed that a linear equilibrium isotherm describes the process. The effects of high salt concentrations are accounted for in the brine density relation, the viscosity relation, Darcy's and Fick's laws, and the rate of mass transfer between macropores and micropores. Mass balance equations, supplemented by extended forms of Darcy's and Fick's laws, are employed to arrive at two sets of equations. One set consists of seven coupled equations for the salt mass fraction and fluid density in macropores, salt mass fraction in micropores, fluid velocity vector, and the fluid pressure. The other set consists of two coupled equations to be solved for the mass fractions of low-concentration species in micropores and macropores. Based on these equations, a mathematical model called TORISM is developed. Using this model, the potential significance of modifications to Darcy's Law are demonstrated.     

Fractal porous media III: Transversal Stokes flow through random and Sierpinski carpets

The transversal Stokes flow of a Newtonian fluid through random and Sierpinski carpets is numerically calculated and the transversal permeability derived. In random carpets derived from site percolation, the average macroscopic permeability varies as (- _c)^3/2, close to the critical porosity _c. This exponent is found to be slightly different from the conductivity exponent. Results for Sierpinski carpets are presented up to the fourth generation. The Carman equation is not verified in these two model porous media.     

Motion through spherical droplet with non-homogenous porous layer in spherical container

The problem of the creeping flow through a spherical droplet with a nonhomogenous porous layer in a spherical container has been studied analytically. Darcy’s model for the flow inside the porous annular region and the Stokes equation for the flow inside the spherical cavity and container are used to analyze the flow. The drag force is exerted on the porous spherical particles enclosing a cavity, and the hydrodynamic permeability of the spherical droplet with a non-homogeneous porous layer is ca...     

Percolation theory approach to transport phenomena in porous media

The percolation theory approach to static and dynamic properties of the single- and two-phase fluid flow in porous media is described. Using percolation cluster scaling laws, one can obtain functional relations between the saturation fraction of a given phase and the capillary pressure, the relative permeability, and the dispersion coefficient, in drainage and imbibition processes. In addition, the scale dependency of the transport coefficient is shown to be an outcome of the fractal nature of pore space and of the random flow pattern of the fluids or contaminant.     

Flow in porous media I: A theoretical derivation of Darcy's law

Stokes flow through a rigid porous medium is analyzed in terms of the method of volume averaging. The traditional averaging procedure leads to an equation of motion and a continuity equation expressed in terms of the volume-averaged pressure and velocity. The equation of motion contains integrals involving spatial deviations of the pressure and velocity, the Brinkman correction, and other lower-order terms. The analysis clearly indicates why the Brinkman correction should not be used to accommodate ano slip condition at an interface between a porous medium and a bounding solid surface.The presence of spatial deviations of the pressure and velocity in the volume-averaged equations of motion gives rise to aclosure problem, and representations for the spatial deviations are derived that lead to Darcy's law. The theoretical development is not restricted to either homogeneous or spatially periodic porous media; however, the problem ofabrupt changes in the structure of a porous medium is not considered.Roman Letters A interfacial area of the - interface contained within the macroscopic system, m² - A _e area of entrances and exits for the -phase contained within the macroscopic system, m² - A interfacial area of the - interface contained within the averaging volume, m² - A ^* interfacial area of the - interface contained within a unit cell, m² - Ae area of entrances and exits for the -phase contained within a unit cell, m² - B second order tensor used to represent the velocity deviation (see Equation (3.30)) - b vector used to represent the pressure deviation (see Equation (3.31)), m^–1 - d distance between two points at which the pressure is measured, m - g gravity vector, m/s² - K Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the -phase (see Figure 2), m - characteristic length scale for the -phase (see Figure 2), m - n unit normal vector pointing from the -phase toward the -phase (n =–n ) - n _e unit normal vector for the entrances and exits of the -phase contained within a unit cell - p pressure in the -phase, N/m² - p intrinsic phase average pressure for the -phase, N/m² - p –p , spatial deviation of the pressure in the -phase, N/m² - r ₀ radius of the averaging volume and radius of a capillary tube, m - v velocity vector for the -phase, m/s - v phase average velocity vector for the -phase, m/s - v intrinsic phase average velocity vector for the -phase, m/s - v –v , spatial deviation of the velocity vector for the -phase, m/s - V averaging volume, m³ - V volume of the -phase contained within the averaging volume, m³ Greek Letters V/V, volume fraction of the -phase - mass density of the -phase, kg/m³ - viscosity of the -phase, Nt/m² - arbitrary function used in the representation of the velocity deviation (see Equations (3.11) and (B1)), m/s - arbitrary function used in the representation of the pressure deviation (see Equations (3.12) and (B2)), s^–1     

Fractal porous media IV: Three-dimensional stokes flow through random media and regular fractals

The three-dimensional Stokes flow of a Newtonian fluid through random and/or fractal media is numerically determined. The permeability of these media is derived. Results relative to these structures are presented and discussed. The validity of the Carman equation and of a simple scaling argument is questioned.     

Fractal porous media I: Longitudinal stokes flow in random carpets

Two-dimensional porous media whose random cross-sections are derived from site percolation are constructed. The longitudinal flow of a Newtonian fluid in the Stokes approximation is then computed and the longitudinal permeability is obtained. Two methods are used and yield the same result when porosity is low. The Carman equation is shown to apply within ±7% when porosity is within the range from 0 to 0.75. Finally, random structures derived from stick percolation are investigated; results are qualitatively the same, but the Carman equation yields a poorer approximation.     

Flow in porous media II: The governing equations for immiscible,two-phase flow

The Stokes flow of two immiscible fluids through a rigid porous medium is analyzed using the method of volume averaging. The volume-averaged momentum equations, in terms of averaged quantities and spatial deviations, are identical in form to that obtained for single phase flow; however, the solution of the closure problem gives rise to additional terms not found in the traditional treatment of two-phase flow. Qualitative arguments suggest that the nontraditional terms may be important when / is of order one, and order of magnitude analysis indicates that they may be significant in terms of the motion of a fluid at very low volume fractions. The theory contains features that could give rise to hysteresis effects, but in the present form it is restricted to static contact line phenomena.Roman Letters (, = , , and ) A interfacial area of the- interface contained within the macroscopic system, m² - A _e area of entrances and exits for the -phase contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A ^* interfacial area of the- interface contained within a unit cell, m² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - g gravity vector, m²/s - H mean curvature of the- interface, m^–1 - H area average of the mean curvature, m^–1 - H – H , deviation of the mean curvature, m^–1 - I unit tensor - K Darcy's law permeability tensor, m² - K permeability tensor for the-phase, m² - K viscous drag tensor for the-phase equation of motion - K viscous drag tensor for the-phase equation of motion - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the-phase, m - n unit normal vector pointing from the-phase toward the-phase (n = –n ) - p _c p – P , capillary pressure, N/m² - p pressure in the-phase, N/m² - p intrinsic phase average pressure for the-phase, N/m² - p – p , spatial deviation of the pressure in the-phase, N/m² - r ₀ radius of the averaging volume, m - t time, s - v velocity vector for the-phase, m/s - v phase average velocity vector for the-phase, m/s - v intrinsic phase average velocity vector for the-phase, m/s - v – v , spatial deviation of the velocity vector for the-phase, m/s - V averaging volume, m³ - V volume of the-phase contained within the averaging volume, m³ Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - viscosity of the-phase, Nt/m² - surface tension of the- interface, N/m - viscous stress tensor for the-phase, N/m² - / kinematic viscosity, m²/s     

Transient oscillatory and constantly accelerated non-Newtonian flow in a porous medium

This paper looks at the magnetohydrodynamic (MHD) analysis for transient flow of an Oldroyd-B fluid in a porous medium. The presented analysis takes into account the modified Darcy's law. The flow is induced due to constantly accelerated and oscillating plate. Expressions for the corresponding velocity field and the adequate tangential stress are determined by means of the Fourier sine transform. The influence of various parameters of interest on the velocity and tangential stress has been shown and discussed. A comparison for different kinds of fluids is also provided.     

Stokes’ first problem for a second grade fluid in a porous half-space with heated boundary

Based on a modified Darcy's law, Stokes’ first problem was investigated for a second grade fluid in a porous half-space with a heated flat plate. Exact solutions of the velocity and temperature fields were obtained using Fourier sine transforms. In contrast to the classical Stokes’ first problem, there is a steady-state solution for the second grade fluid in the porous half-space, which is a damping exponential function with respect to the distance from the flat plate. The well-known solutions for Newtonian fluids in non-porous or porous half-space appear in limiting cases of our solutions.     

The thermogravitational effect in porous media: A modelling approach

The thermogravitational effect may induce large concentration contrasts, particularly in porous media. This phenomenon arises from a coupling of the Soret effect and convection currents in a temperature field. The present study of this phenomenon is motivated by the safety assessment of nuclear waste repositories, which are sources of thermal energy. Here, we present a modelling approach of laboratory experiments carried out at the University of Toulouse. The results of this model, though more adequate than the analytical solution to account for the influence of permeability, remain far from the experimental ones. In conclusion, it appears that the research must now focus on both a comprehensive phenomenology of the transport processes and experiments with new dimensional constraints.Notation C concentration of the solute - c mass fraction of the solute - D molecular diffusion coefficient - D coefficient of thermodiffusion - F ₁ external forces on the solute - F ₂ external forces on the solvent - g gravity - J _q heat flux - J _l solute mass flux - J _x horizontal mass flux of the solute - k intrinsic permeability - k _L longitudinal permeability - k _T transverse permeability - T ₀ reference temperature - U Darcy velocity - thermal expansion coefficient of the mixture - thermal conductivity - porosity of the medium - dynamic viscosity of the mixture - _T ⁰ dynamic viscosity of pure water atT - _1T dynamic viscosity of 1 atT - _2T dynamic viscosity of 2 atT - ₁ chemical potential of the solute - ₂ chemical potential of the solvent - specific mass of the fluid - ₀ specific mass of pure water atT ₀