An Analysis of the Movement of Wetting and Nonwetting Fluids in Homogeneous Porous Media

The movement of wetting and nonwetting fluid flow in columns packed with glass beads is used to understand the more complicated flows in homogeneous porous media. The motion of two immiscible liquids (oil and water) is observed with different surfactants. Through dimensional analyses, fluid velocity is well correlated with interfacial tension and less dependent on particle size. In water–oil (W/O) experiments, finger pattern flows are observed if water is the displacing fluid that flows in an oil-filled porous media, whereas oil ganglia tend to form if oil is the displacing fluid in the water-wetted porous media. The results are well described by a simple model based on an earlier theory of flow in a tube.     

Wetting and nonwetting fluid displacements in porous media

The understanding of simple laminar flow in tubes has often been used to interpret the more complicated flow in porous media. A study of the motion of two immiscible liquids in closed tubes with relatively large diameter (> 0.3 cm i.d), was conducted in order to examine the influence of wetting and nonwetting liquids on the flow behavior. The results indicate that the wetting properties of the fluids with regard to the tube wall have a major efffect on the formation and motion of long bubbles. A physically based model was used to predict the velocity and the conditions for no motion of bubbles and drops in tubes. These results were used to interpret the nature of oil and water flow in porous media. Experiments in which the wetting liquid was displaced by the nonwetting, or vice versa, were conducted by injecting the displacing liquid at a constant flux at the center of a two-dimensional chamber saturated with the displaced liquid. The influence of wetting-nonwetting characteristics on the quantity of liquid displaced, the shape of the interface between the two liquids, and the interpretation of the no motion radius in a closed tube to the case of a porous medium are discussed. It would appear that the no motion radius gives a good indication of the minimum width of a nonwetting penetrating finger and the maximum width of nonwetting ganglia left by drainage.     

Wetting front instability in unsaturated porous media: A three-dimensional study in initially dry sand

Three-dimensional fingers caused by gravity driven, wetting front instability in unsaturated porous media were explored through laboratory experimentation. Two sets of experiments were conducted using initially dry sand in large 30 cm diameter columns to guide analytical development for finger velocity and diameter. The first set consisted of ponding water on a two-layer sand system with a fine sand overlying a coarse sand. Here, a complicated pattern of interaction among fingers was found to occur. In the second set, water was applied directly to the coarse layer as 2 cm diameter area sources, enabling systematic study of individual finger structure. Based on dimensional analysis and the experimental results, general relationships were found for finger velocity and cross-sectional area as a function of both the flux through the finger and hydraulic properties of the coarse layer represented by the sorptivity, saturated conductivity, and initial and saturated moisture contents. Unlike the two-dimensional case, air entrapment was an important factor in explaining the experimental results.     

Approximate Analytical Solution of the Nonlinear Diffusion Equation for Arbitrary Boundary Conditions

A general approximation for the solution of the one-dimensional nonlinear diffusion equation is presented. It applies to arbitrary soil properties and boundary conditions. The approximation becomes more accurate when the soil-water diffusivity approaches a delta function, yet the result is still very accurate for constant diffusivity suggesting that the present formulation is a reliable one. Three examples are given where the method is applied, for a constant water content at the surface, when a saturated zone exists and for a time-dependent surface flux.