On the Validity of Darcy's Law for Stable High-Concentration Displacements in Granular Porous Media

Recently, it has been suggested that Darcy's Law might not be applicable for modelling miscible, density-dependent flow in porous media. To investigate this, three sets of careful laboratory column experiments were performed on coarse and medium sands, consisting of upward displacement of water by sodium chloride solutions with concentrations ranging from 5 to 200g/l. Data on salt concentrations and water pressures were collected in horizontal transects along the flow direction. Salt concentration data were also collected in the influent and exit lines. The experimental data were analysed using a simplified approach based on Darcy's Law alone, applied with the assumption of a sharp interface. Darcy's Law was used to estimate porous medium permeability by fitting predictions to experimental data. Consistent estimates of permeability were obtained for each set of experiments. The results indicate that Darcy's Law adequately describes high concentration displacements through saturated coarse- and medium-grained porous media.     

On the transient non-Fickian dispersion theory

The Fickian dispersion equation is the basic relationship used to describe the nonconvective mass flux of a solute in a porous medium. This equation prescribes a linear relationship between the dispersive mass flux and the concentration gradient. An important characteristic of the Fickian relationship is that it is independent of the history of dispersion (e.g. the time rate of change of the dispersion flux). Also, the dispersivities are supposed to be medium constants and invariant with temporal and spatial scales of observation. It is believed that in general these restrictions do not hold. A number of authors have proposed various alternative relationships. For example, differential equations have been employed that prescribe a relationship between the dispersion flux and its time and space derivatives. Also, stochastic theories result in integro-differential equations in which dispersion tensor grow asymptotically with time or distance. In this work, three different approaches, which lead to three different non-Fickian equations with a transient character, are discussed and their primary features and differences are highlighted. It is shown that an effective dispersion tensor defined in the framework of the transient non-Fickian theory, grows asymptotically with time and distance; a result which also follows from stochastic theories. Next, principles of continuum mechanics are employed to provide a solid theoretical basis for the non-Fickian transient dispersion theory. The equation of motion of a solute in a porous medium is used to provide a rigorous derivation of various dispersion relationships valid under different conditions. Under various simplifying assumptions, the generalized theory is found to agree with the conventional Fickian theory as well as several other non-Fickian relationships found in the literature. Moreover, it is shown that for nonconservative solutes, the traditional dispersion tensor is affected by the rate of mass exchange of the solute.Also with National Institute of Public Health and Environmental Protection (RIVM), PO Box 1; 3720BA Bilthoven, The Netherlands     

Macro-Scale Dynamic Effects in Homogeneous and Heterogeneous Porous Media

It is known that the classical capillary pressure-saturation relationship may be deficient under non-equilibrium conditions when large saturation changes may occur. An extended relationship has been proposed in the literature which correlates the rate of change of saturation to the difference between the phase pressures and the equilibrium capillary pressure. This linear relationship contains a damping coefficient, \tau, that may be a function of saturation. The extended relationship is examined at the macro-scale through simulations using the two-phase simulator MUFTE-UG. In these simulations, it is assumed that the traditional equilibrium relationship between the water saturation and the difference in fluid pressures holds locally. Steady-state and dynamic numerical experiments are performed where a non-wetting phase displaces a wetting phase in homogeneous and heterogeneous domains with varying boundary conditions, domain size, and soil parameters. From these simulations the damping coefficient can be identified as a (non-linear) function of the water saturation. It is shown that the value of increases with an increased domain size and/or with decreased intrinsic permeability. Also, the value of for a domain with a spatially correlated random distribution of intrinsic permeability is compared to a homogeneous domain with equivalent permeability; they are shown to be almost equal.     

Effective Stress in Unsaturated Soils: A Thermodynamic Approach Based on the Interfacial Energy and Hydromechanical Coupling

In recent years, the effective stress approach has received much attention in the constitutive modeling of unsaturated soils. In this approach, the effective stress parameter is very important. This parameter needs a correct definition and has to be determined properly. In this paper, a thermodynamic approach is used to develop a physically-based formula for the effective stress tensor in unsaturated soils. This approach accounts for the hydro-mechanical coupling, which is quite important when dealing with hydraulic hysteresis in unsaturated soils. The resulting formula takes into account the role of interfacial energy and the contribution of air?Cwater specific interfacial area to the effective stress tensor. Moreover, a bi-quadratic surface is proposed to represent the contribution of the so-called suction stress in the effective stress tensor. It is shown that the proposed relationship for suction stress is in agreement with available experimental data in the full hydraulic cycle (drying, scanning, and wetting).     

Horizontal Redistribution of Two Fluid Phases in a Porous Medium: Experimental Investigations

Classical models for flow and transport processes in porous media employ the so-called extended Darcy’s Law. Originally, it was proposed empirically for one-dimensional isothermal flow of an incompressible fluid in a rigid, homogeneous, and isotropic porous medium. Nowadays, the extended Darcy’s Law is used for highly complex situations like non-isothermal, multi-phase and multi-component flow and transport, without introducing any additional driving forces. In this work, an alternative approach by Hassanizadeh and Gray identifying additional driving forces were tested in an experimental setup for horizontal redistribution of two fluid phases with an initial saturation discontinuity. Analytical and numerical solutions based on traditional models predict that the saturation discontinuity will persist, but a uniform saturation distribution will be established in each subdomain after an infinite amount of time. The pressure field, however, is predicted to be continuous throughout the domain at all times and is expected to become uniform when there is no flow. In our experiments, we also find that the saturation discontinuity persists. But, gradients in both saturation and pressure remain in both subdomains even when the flow of fluids stops. This indicates that the identified additional driving forces present in the truly extended Darcy’s Law are potentially significant.     

A novel deep reactive ion etched (DRIE) glass micro-model for two-phase flow experiments

In the last few decades, micro-models have become popular experimental tools for two-phase flow studies. In this work, the design and fabrication of an innovative, elongated, glass-etched micro-model with dimensions of 5 × 35 mm(2) and constant depth of 43 microns is described. This is the first time that a micro-model with such depth and dimensions has been etched in glass by using a dry etching technique. The micro-model was visualized by a novel setup that allowed us to monitor and record the distribution of fluids throughout the length of the micro-model continuously. Quasi-static drainage experiments were conducted in order to obtain equilibrium data points that relate capillary pressure to phase saturation. By measuring the flow rate of water through the flow network for known pressure gradients, the intrinsic permeability of the micro-model's flow network was also calculated. The experimental results were used to calibrate a pore-network model and test its validity. Finally, we show that glass-etched micro-models can be valuable tools in single and/or multi-phase flow studies and their applications.     

Insights into the Relationships Among Capillary Pressure,Saturation, Interfacial Area and Relative Permeability Using Pore-Network Modeling

To gain insight in relationships among capillary pressure, interfacial area, saturation, and relative permeability in two-phase flow in porous media, we have developed two types of pore-network models. The first one, called tube model, has only one element type, namely pore throats. The second one is a sphere-and-tube model with both pore bodies and pore throats. We have shown that the two models produce distinctly different curves for capillary pressure and relative permeability. In particular, we find that the tube model cannot reproduce hysteresis. We have investigated some basic issues such as effect of network size, network dimension, and different trapping assumptions in the two networks. We have also obtained curves of fluid–fluid interfacial area versus saturation. We show that the trend of relationship between interfacial area and saturation is largely influenced by trapping assumptions. Through simulating primary and scanning drainage and imbibition cycles, we have generated two surfaces fitted to capillary pressure, saturation, and interfacial area (P _c –S _w –a _nw ) points as well as to relative permeability, saturation, and interfacial area (k _r –S _w –a _nw ) points. The two fitted three-dimensional surfaces show very good correlation with the data points. We have fitted two different surfaces to P _c –S _w –a _nw points for drainage and imbibition separately. The two surfaces do not completely coincide. But, their mean absolute difference decreases with increasing overlap in the statistical distributions of pore bodies and pore throats. We have shown that interfacial area can be considered as an essential variable for diminishing or eliminating the hysteresis observed in capillary pressure–saturation (P _c –S _w ) and the relative permeability–saturation (k _r –S _w ) curves.     

Modeling Kinetic Interphase Mass Transfer for Two-Phase Flow in Porous Media Including Fluid–Fluid Interfacial Area

Interphase mass transfer in porous media takes place across fluid–fluid interfaces. At the field scale, this is almost always a kinetic process and its rate is highly dependent on the amount of fluid–fluid interfacial area. Having no means to determine the interfacial area, modelers usually either neglect kinetics of mass transfer and assume local equilibrium between phases or they estimate interfacial area using lumped parameter approaches (in DNAPL pool dissolution) or a dual domain approach (for air sparging). However, none of these approaches include a physical determination of interfacial area or accounts for its role for interphase mass transfer. In this work, we propose a new formulation of two-phase flow with interphase mass transfer, which is based on thermodynamic principles. This approach comprises a mass balance for each component in each phase and a mass balance for specific interfacial area. The system is closed by a relationship among capillary pressure, interfacial area, and saturation. We compare our approach to an equilibrium model by showing simulation results for an air–water system. We show that the new approach is capable of modeling kinetic interphase mass exchange for a two-phase system and that mass transfer correlates with the specific interfacial area. By non-dimensionalization of the equations and variation of Peclet and Damköhler number, we make statements about when kinetic interphase mass transfer has to be taken into account by using the new physically based kinetic approach and when the equilibrium model is a reasonable simplification.