Pore-scale imbibition experiments in dry and prewetted porous media

Two-phase imbibition behavior of immiscible fluids was studied in dry and prewetted porous media using a laser-induced fluorescence technique. Imbibition was first investigated in two-dimensional (2-D) systems under conditions comparable to those for a study of drainage [Ovdat H, Berkowitz B. Pore-scale study of drainage displacement under combined capillary and gravity effects in index-matched porous media. Water Resources Research 2006;42:W06411. doi:10.1029/2005WR004553] in the capillary-dominated regime. The effect of initial wetting saturation (IWS) was then explored in 2-D and 3-D porous media under the combined effect of gravity, capillary and viscous forces, within and outside the capillary-dominated regime. Parameters that describe maximum vertical advance, volumetric fraction, total surface area and specific surface area of the invading fluid were used to quantify the behavior. Comparison of 2-D drainage and imbibition patterns demonstrates significant qualitative differences under analogous viscosity ratio, buoyancy number, and capillary number values. However, quantitative analyses show strong pore-scale similarities between these patterns. Invasion structures in 3-D, prewetted (IWS ≈ 8% of the pore volume) porous media are ramified, with lateral branching and regions containing trapped residual fluid. These structures are qualitatively and quantitatively different from the compact, branchless structures that develop in dry (IWS = 0) porous media.     

Effect of structure and texture on infiltration flow pattern during flood irrigation

Understanding the flow pattern of water and solute in subsurface soils is critically important in the fields of agricultural and environmental sciences. Dye tracer tests using a flood irrigation of Brilliant Blue FCF solution (5 g l^-1) and excavation method was performed to investigate the effect of texture and structure on the infiltration pattern at three different field soils developed from granite (GR), gneiss (GN) and limestone (LS). The GR soil showed a homogeneous matrix flow in the surface soil with weak, medium granular structure and a macropore flow along pegmatitic vein and plant root in C horizon. The surface horizon (A1) of GN soil with moderate, medium granular structure and many fine roots had matrix flow. The fingering occurred at the interfaces of sandy loam A horizon and loamy sand C horizon in GR soil and loam A1 horizon and sandy loam A2 horizon in GN soil. The LS soil with strong, coarse prismatic structure and the finest texture showed a macropore flow along cracks and had the deepest penetration of the dye tracer. The macropore (crack and vein), layer interface and plant root induced the preferential flow in the studied soils.     

Fingered flow: The creator of sand columns in dune and beach sands

Sand columns, sand cones, sand mushrooms and other striking sand forms are frequently observed in the Dutch and German beach and dune sands. This paper aims to clarify the mechanism of sand column formation. Recently it has become evident that homogeneous beach and dune sands often become irregularly wetted by infiltrating rainwater. In otherwise dry sandy soils, wet preferential flow paths (‘fingers’) may develop. At two test sites the volumetric soil moisture content varied between 0·2 and 12·0 per cent. The wet fingers represent the premature state of sand columns. When the dry sand in between these fingers is blown away by the wind, the more resistant wet sand of the fingers will remain in its place and appear as sand columns at the surface. As a result of wind and erosive sand drifts, striking sand forms may be formed.     

Variable-density flow and transport in porous media: approaches and challenges

We review the state of the art in modeling of variable-density flow and transport in porous media, including conceptual models for convection systems, governing balance equations, phenomenological laws, constitutive relations for fluid density and viscosity, and numerical methods for solving the resulting nonlinear multifield problems. The discussion of numerical methods addresses strategies for solving the coupled spatio-temporal convection process, consistent velocity approximation, and error-based mesh adaptation techniques. As numerical models for those nonlinear systems must be carefully verified in appropriate tests, we discuss weaknesses and inconsistencies of current model-verification methods as well as benchmark solutions. We give examples of field-related applications to illustrate specific challenges of further research, where heterogeneities and large scales are important.