Nonlocal interface dynamics and pattern formation in gravity-driven unsaturated flow through porous media.

Existing continuum models of multiphase flow in porous media are unable to explain why preferential flow (fingering) occurs during infiltration into homogeneous, dry soil. Following a phase-field methodology, we propose a continuum model that accounts for an apparent surface tension at the wetting front and does not introduce new independent parameters. The model reproduces the observed features of fingered flows, in particular, the higher saturation of water at the tip of the fingers, which is believed to be essential for the formation of fingers. From a linear stability analysis, we predict that finger velocity and finger width both increase with infiltration rate, and the predictions are in quantitative agreement with experiments.

[1]  Stokes,et al.  Dynamic capillary pressure in porous media: Origin of the viscous-fingering length scale. , 1987, Physical review letters.

[2]  J. Selker,et al.  Evaluation of hydrodynamic scaling in porous media using finger dimensions , 1998 .

[3]  Eric A. Davidson,et al.  Control of cation concentrations in stream waters by surface soil processes in an Amazonian watershed , 2001, Nature.

[4]  J. Parlange,et al.  Theoretical analysis of wetting front instability in soils , 1976 .

[5]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[6]  D. R. Nielsen,et al.  Water movement in glass bead porous media: 2. Experiments of infiltration and finger flow , 1994 .

[7]  William G. Gray,et al.  Thermodynamic basis of capillary pressure in porous media , 1993 .

[8]  Ruben Juanes,et al.  Nonmonotonic traveling wave solutions of infiltration into porous media , 2008 .

[9]  Heike Emmerich,et al.  Advances of and by phase-field modelling in condensed-matter physics , 2008 .

[10]  J. Langer,et al.  Pattern formation in nonequilibrium physics , 1999 .

[11]  Andrea L. Bertozzi,et al.  Linear stability and transient growth in driven contact lines , 1997 .

[12]  Tammo S. Steenhuis,et al.  Wetting front instability: 2. Experimental determination of relationships between system parameters and two‐dimensional unstable flow field behavior in initially dry porous media , 1989 .

[13]  M. Blunt,et al.  Determination of finger shape using the dynamic capillary pressure , 2000 .

[14]  Tammo S. Steenhuis,et al.  Fingered flow in two dimensions: 2. Predicting finger moisture profile , 1992 .

[15]  R. Glass,et al.  On the porous‐continuum modeling of gravity‐driven fingers in unsaturated materials: Extension of standard theory with a hold‐back‐pile‐up effect , 2002 .

[16]  Y. Yortsos,et al.  INVASION PERCOLATION WITH VISCOUS FORCES , 1998 .

[17]  M. Sahimi.  Flow phenomena in rocks : from continuum models to fractals, percolation, cellular automata, and simulated annealing , 1993 .

[18]  J. Nieber,et al.  Dynamic Capillary Pressure Mechanism for Instability in Gravity-Driven Flows; Review and Extension to Very Dry Conditions , 2005 .

[19]  R. Lenormand.  Pattern growth and fluid displacements through porous media , 1986 .

[20]  Herbert E. Huppert,et al.  Flow and instability of a viscous current down a slope , 1982, Nature.

[21]  G. Taylor,et al.  The penetration of a fluid into a porous medium or Hele-Shaw cell containing a more viscous liquid , 1958, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[22]  P. Gennes.  Wetting: statics and dynamics , 1985 .

[23]  Liquid Conservation and Nonlocal Interface Dynamics in Imbibition , 1999, cond-mat/9907394.

[24]  Chao-Yang Wang,et al.  Fundamental models for fuel cell engineering. , 2004, Chemical reviews.

[25]  Anne E. Trefethen,et al.  Hydrodynamic Stability Without Eigenvalues , 1993, Science.