Three‐dimensional simulation of liquid drop dynamics within unsaturated vertical Hele‐Shaw cells

[1] A three-dimensional, multiphase fluid flow model with volume of fluid-interface tracking was developed and applied to study the multiphase dynamics of moving liquid drops of different sizes within vertical Hele-Shaw cells. The simulated moving velocities are significantly different from those obtained from a first-order analytical approximation, based on simple force-balance concepts. The simulation results also indicate that the moving drops can exhibit a variety of shapes and that the transition among these different shapes is largely determined by the moving velocities. More important, there is a transition from a linear moving regime at small capillary numbers, in which the capillary number scales linearly with the Bond number, to a nonlinear moving regime at large capillary numbers, in which the moving drop releases a train of droplets from its trailing edge. The train of droplets forms a variety of patterns at different moving velocities.

[1]  D. B. Kothe,et al.  RIPPLE: A Computer Program for Incompressible Flows with Free Surfaces , 1991 .

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

[3]  J. Bikerman Sliding of drops from surfaces of different roughnesses , 1950 .

[4]  J. Brackbill,et al.  A continuum method for modeling surface tension , 1992 .

[5]  C. Furmidge,et al.  Studies at phase interfaces. I. The sliding of liquid drops on solid surfaces and a theory for spray retention , 1962 .

[6]  S. Wheatcraft,et al.  Gravity-driven infiltration instability in initially dry nonhorizontal fractures , 1994 .

[7]  R. G. Cox The dynamics of the spreading of liquids on a solid surface. Part 1. Viscous flow , 1986, Journal of Fluid Mechanics.

[8]  J. Wan,et al.  Water film flow along fracture surfaces of porous rock , 1997 .

[9]  C. Doughty,et al.  Liquid-release tests in unsaturated fractured welded tuffs: I. Field investigations , 2002 .

[10]  R. Hoffman A study of the advancing interface. I. Interface shape in liquid—gas systems , 1975 .

[11]  M. J. Nicholl,et al.  Infiltration into an Analog Fracture: Experimental Observations of Gravity‐Driven Fingering , 2005 .

[12]  L. Tanner,et al.  The spreading of silicone oil drops on horizontal surfaces , 1979 .

[13]  Reply to comment by Dani Or and Teamrat Ghezzehei on “Computer simulation of two‐phase immiscible fluid motion in unsaturated complex fractures using a volume of fluid method” by Hai Huang, Paul Meakin, and Moubin Liu , 2006 .

[14]  L. Limat,et al.  Shape and motion of drops sliding down an inclined plane , 2005, Journal of Fluid Mechanics.

[15]  Fracture intersections as integrators for unsaturated flow , 2002 .

[16]  K. Pruess,et al.  Solute transport along preferential flow paths in unsaturated fractures , 2001 .

[17]  J. Brock,et al.  Volume tracking of interfaces having surface tension in two and three dimensions , 1996 .

[18]  Maria Ines Dragila,et al.  Fluid motion through an unsaturated fracture junction , 2004 .

[19]  B. Andreotti,et al.  Boundary conditions in the vicinity of a dynamic contact line: experimental investigation of viscous drops sliding down an inclined plane. , 2005, Physical review letters.

[20]  R. G. Cox The dynamics of the spreading of liquids on a solid surface. Part 2. Surfactants , 1986, Journal of Fluid Mechanics.

[21]  R. Chow,et al.  On the ability of drops or bubbles to stick to non-horizontal surfaces of solids , 1983, Journal of Fluid Mechanics.

[22]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[23]  R. Glass,et al.  Liquid phase structure within an unsaturated fracture network beneath a surface infiltration event: Field experiment , 2002 .

[24]  P. Meakin,et al.  Computer simulation of two‐phase immiscible fluid motion in unsaturated complex fractures using a volume of fluid method , 2005 .

[25]  Teamrat A. Ghezzehei,et al.  Liquid fragmentation and intermittent flow regimes in unsaturated fractured media , 2005 .

[26]  Comment on “Computer simulation of two‐phase immiscible fluid motion in unsaturated complex fractures using a volume of fluid method” by Hai Huang, Paul Meakin, and Moubin Liu , 2006 .

[27]  Thomas Podgorski,et al.  Corners, cusps, and pearls in running drops. , 2001, Physical review letters.

[28]  O. Voinov Hydrodynamics of wetting , 1976 .

[29]  T. Inamuro,et al.  A lattice Boltzmann method for incompressible two-phase flows with large density differences , 2004 .

[30]  D. v. On the ability of drops or bubbles to stick to non-horizontal surfaces of solids. Part 2. Small drops or bubbles having contact angles of arbitrary size , 1985, Journal of Fluid Mechanics.

[31]  Karsten Pruess,et al.  Two‐Phase Flow Visualization and Relative Permeability Measurement in Natural Rough‐Walled Rock Fractures , 1995 .

[32]  N. Weisbrod,et al.  Parameters affecting maximum fluid transport in large aperture fractures , 2003 .

[33]  K. Pruess,et al.  Experimental studies of water seepage and intermittent flow in unsaturated, rough‐walled fractures , 1999 .

[34]  Robert Shaw,et al.  The Dripping Faucet As A Model Chaotic System , 1984 .