Modeling of surface tension and contact angles with smoothed particle hydrodynamics.

A two-dimensional numerical model based on smoothed particle hydrodynamics (SPH) was used to simulate unsaturated (multiphase) flow through fracture junctions. A combination of standard SPH equations with pairwise fluid-fluid and fluid-solid particle-particle interactions allowed surface tension and three-phase contact dynamics to be simulated. The model was validated by calculating the surface tension in four different ways: (i) from small-amplitude oscillations of fluid drops, (ii) from the dependence of the capillary pressure on drop radius, (iii) from capillary rise simulations, and (iv) from the behavior of a fluid drop confined between parallel walls under the influence of gravity. All four simulations led to consistent values for the surface tension. The dependence of receding and advancing contact angles on droplet velocity was studied. Incorporation of surface tension and fluid-solid interactions allowed unsaturated flow through fracture junctions to be realistically simulated, and the simulation results compare well with the laboratory experiments of Dragila and Weisbrod.

[1]  H. Posch,et al.  Steady-state shear flows via nonequilibrium molecular dynamics and smooth-particle applied mechanics. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[2]  R. H. Dettre,et al.  Contact Angle Hysteresis. III. Study of an Idealized Heterogeneous Surface , 1964 .

[3]  J. Morris,et al.  Modeling Low Reynolds Number Incompressible Flows Using SPH , 1997 .

[4]  J. Monaghan,et al.  Smoothed particle hydrodynamics: Theory and application to non-spherical stars , 1977 .

[5]  J. Monaghan.  Simulating Free Surface Flows with SPH , 1994 .

[6]  Levine,et al.  Interface moving through a random background. , 1985, Physical review. B, Condensed matter.

[7]  M. Kardar,et al.  Dynamic relaxation of drifting polymers: A phenomenological approach. , 1993, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

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

[9]  P. Español,et al.  Statistical Mechanics of Dissipative Particle Dynamics. , 1995 .

[10]  J. Koelman,et al.  Simulating microscopic hydrodynamic phenomena with dissipative particle dynamics , 1992 .

[11]  J. Monaghan.  Smoothed particle hydrodynamics , 2005 .

[12]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[13]  Narayan,et al.  Threshold critical dynamics of driven interfaces in random media. , 1993, Physical review. B, Condensed matter.

[14]  Arturo A. Keller,et al.  Behavior of Nonaqueous Phase Liquids in Fractured Porous Media under Two‐Phase Flow Conditions , 2000 .

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

[16]  Shan,et al.  Lattice Boltzmann model for simulating flows with multiple phases and components. , 1993, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[17]  R. H. Dettre,et al.  Contact Angle Hysteresis. IV. Contact Angle Measurements on Heterogeneous Surfaces1 , 1965 .

[18]  J. Joanny,et al.  Motion of a contact line on a heterogeneous surface , 1990 .

[19]  James J. Feng,et al.  A diffuse-interface method for simulating two-phase flows of complex fluids , 2004, Journal of Fluid Mechanics.

[20]  L. Lucy.  A numerical approach to the testing of the fission hypothesis. , 1977 .

[21]  Wong,et al.  Contact line dynamics near the pinning threshold: a capillary rise and fall experiment , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[22]  Wm. G. Hoover.  Isomorphism linking smooth particles and embedded atoms , 1998 .

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

[24]  Lei-Han Tang,et al.  Dynamics of Interface Depinning in a Disordered Medium , 1992 .

[25]  James M. Keller,et al.  Lattice-Gas Cellular Automata: Inviscid two-dimensional lattice-gas hydrodynamics , 1997 .

[26]  H. Posch,et al.  Liquid drops and surface tension with smoothed particle applied mechanics , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[27]  Banavar,et al.  Lattice Boltzmann study of hydrodynamic spinodal decomposition. , 1995, Physical review letters.

[28]  J. Morris.  Simulating surface tension with smoothed particle hydrodynamics , 2000 .

[29]  G. Gonnella,et al.  Lattice Boltzmann simulations of lamellar and droplet phases , 1998 .

[30]  J. Vearncombe,et al.  Fault/fracture density and mineralization: a contouring method for targeting in gold exploration , 2004 .

[31]  P. Gennes,et al.  Dynamics of wetting with nonideal surfaces. The single defect problem , 1989 .

[32]  Patrick J. Fox,et al.  A pore‐scale numerical model for flow through porous media , 1999 .

[33]  Patrick J. Fox,et al.  Smoothed Particle Hydrodynamics Model for Diffusion through Porous Media , 2001 .

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

[35]  L Fan,et al.  Simulation of contact line dynamics in a two-dimensional capillary tube by the lattice Boltzmann model. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[36]  C. Radke,et al.  Kinetics of liquid/liquid capillary rise: I. Experimental observations , 1986 .

[37]  J. Timonen,et al.  Lattice-Boltzmann Simulation of Capillary Rise Dynamics , 2002 .