Transport in chemically and mechanically heterogeneous porous media—III. Large-scale mechanical equilibrium and the regional form of Darcy's law

1998 ◽  
Vol 21 (7) ◽  
pp. 617-629 ◽  
Author(s):  
Michel Quintard ◽  
Stephen Whitaker
2019 ◽  
Vol 22 (9) ◽  
pp. 1131-1139
Author(s):  
G. Z. Liu ◽  
F. L. Liu ◽  
M. Li ◽  
X. Jin ◽  
W. F. Lv ◽  
...  

Hydrology ◽  
2019 ◽  
Vol 7 (1) ◽  
pp. 4
Author(s):  
R. William Nelson ◽  
Gustavious P. Williams

We present a rigorous mathematical treatment of water flow in saturated heterogeneous porous media based on the classical Navier-Stokes formulation that includes vorticity in a heterogeneous porous media. We used the mathematical approach proposed in 1855 by James Clark Maxwell. We show that flow in heterogeneous media results in a flow field described by a heterogeneous complex lamellar vector field with rotational flows, compared to the homogeneous lamellar flow field that results from Darcy’s law. This analysis shows that Darcy’s Law does not accurately describe flow in a heterogeneous porous medium and we encourage precise laboratory experiments to determine under what conditions these issues are important. We publish this work to encourage others to perform numerical and laboratory experiments to determine the circumstances in which this derivation is applicable, and in which the complications can be disregarded.


2021 ◽  
Author(s):  
Marco Dentz ◽  
Alexandre Puyguiraud ◽  
Philippe Gouze

<p>Transport of dissolved substances through porous media is determined by the complexity of the pore space and diffusive mass transfer within and between pores. The interplay of diffusive pore-scale mixing and spatial flow variability are key for the understanding of transport and reaction phenomena in porous media. We study the interplay of pore-scale mixing and network-scale advection through heterogeneous porous media, and its role for the evolution and asymptotic behavior of hydrodynamic dispersion. In a Lagrangian framework, we identify three fundamental mechanisms of pore-scale mixing that determine large scale particle motion: (i) The smoothing of intra-pore velocity contrasts, (ii) the increase of the tortuosity of particle paths, and (iii) the setting of a maximum time for particle transitions. Based on these mechanisms, we derive an upscaled approach that predicts anomalous and normal hydrodynamic dispersion based on the characteristic pore length, Eulerian velocity distribution and Péclet number. The theoretical developments are supported and validated by direct numerical flow and transport simulations in a three-dimensional digitized Berea sandstone sample obtained using X-Ray microtomography. Solute breakthrough curves, are characterized by an intermediate power-law behavior and exponential cut-off, which reflect pore-scale velocity variability and intra-pore solute mixing. Similarly, dispersion evolves from molecular diffusion at early times to asymptotic hydrodynamics dispersion via an intermediate superdiffusive regime. The theory captures the full evolution form anomalous to normal transport behavior at different Péclet numbers as well as the Péclet-dependence of asymptotic dispersion. It sheds light on hydrodynamic dispersion behaviors as a consequence of the interaction between pore-scale mixing and Eulerian flow variability. </p>


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):  
Michel Quintard ◽  
Stephen Whitaker

Most porous media of practical importance are hierarchical in nature; that is, they are characterized by more than one length-scale. When these length-scales are disparate, the hierarchical structure can be analyzed by the method of volume averaging (Anderson and Jackson, 1967; Marie, 1967; Slattery, 1967; Whitaker, 1967). In this approach, macroscopic quantities at a given length-scale are defined in terms of a boundary value problem that describes the phenomena at a smaller length-scale, and information is filtered from one scale to another by a series of volume and area integrals. Other methods, such as ensemble averaging (Matheron, 1965; Dagan, 1989) or homogenization theory (Bensoussan et al, 1978; Sanchez-Palencia, 1980), have been used to study hierarchical systems, and developments specific to the problems under consideration in this chapter can be found in Bourgeat (1984), Auriault (1987), Amaziane and Bourgeat (1988), and Sáez et al. (1989). The transformation from the Darcy scale to the large scale is a recurrent problem in reservoir and aquifer engineering. A detailed description of reservoir properties is available through geostatistical analysis (Journel, 1996) on a fine grid with a length-scale much smaller than the scale of the blocks in the reservoir simulator. “Effective” or “pseudo” properties are assigned to the coarse grid blocks, while the forms of the large-scale equations are required to be the same as those used at the Darcy scale (Coats et al., 1967; Hearn, 1971; Jacks et al., 1972; Kyte and Berry, 1975; Dake, 1978; Killough and Foster, 1979; Yokoyama and Lake, 1981; Kortekaas, 1983; Thomas, 1983; Kossack et al., 1990). A detailed discussion of the comparison between the several approaches is beyond the scope of this chapter; however, one can read Bourgeat et al. (1988) for an introductory comparison between the method of volume averaging and the homogenization theory, and Ahmadi et al. (1993) for a discussion of the various classes of pseudofunction theories.


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).


Author(s):  
Tian-Chyi Yeh ◽  
Raziuddin Khaleel ◽  
Kenneth C. Carroll

Author(s):  
Stefan Doser ◽  
Sang-Joon John Lee

This work investigates the special case of in-plane fluid flow of a Newtonian incompressible fluid at low Reynolds numbers across a paper-thin porous medium in a confined conduit. Fluid transport in sheets with these characteristics are used in emerging devices such as microscale paper-based analytical devices (μPADs) and “e-paper” displays. Darcy’s law is applied and tested to determine if experimentally measured pressures at two flow rates of 5 μL/min and 10 μL/min agree with predicted values. A test device was designed using kinematic design principles to ensure a deterministic 318 μm gap that directs prescribed flow, unidirectionally across porous filter paper. The paper used was Grade 50 Whatman filter paper with an average pore size of 2.7 μm. Pressure was measured along the direction of flow over a 125 mm distance by six pressure ports placed at uniform increments of 25 mm to determine a profile of pressure along the flow path. Measurements were recorded at discrete time intervals over a period up to 48 hours with at least four replicates. Experimental measurements of the pressure profile show a linear relationship as predicted by Darcy’s law, allowing material permeability to be calculated. Among replicates measured under the same set of controllable conditions, experimental data also show a nonlinear relationship. The nonlinearity suggests evidence of transition into an inertia region, providing insight into the factors and behavior of the Darcy-Forchheimer transition for this special case of porous media flow.


Sign in / Sign up

Export Citation Format

Share Document