Numerical simulations of capillary-driven flows in nonuniform cross-sectional capillaries.

In this study the wetting behavior of converging-diverging and diverging-converging capillaries is investigated numerically using an in-house written, finite-element code. An interface tracking procedure based on the predicted change in the total liquid volume, to update the interface location, and Cox's formulation, to determine the dynamic contact angle and the interface shape, is proposed and used. Flow simulations revealed that both converging-diverging and diverging-converging capillaries exhibit significantly slower wetting behavior than straight capillaries and that any deviation in the capillary diameter necessarily tends to slow the overall wetting speed. This behavior was attributed to local regions of very low capillary pressure and high viscous retardation force when the capillary diameter at the interface was significantly larger than the capillary diameter over the upstream fluid. Though the local wetting velocities were different, when equivalent capillaries were compared it was found that both converging-diverging and diverging-converging capillaries had the same total fill time independent of the number of irregular regions, suggesting that the simple model is sufficient for predicting the overall effect. The influence of surface tension and contact angle on the total wetting time was found to be similar for both straight and irregularly shaped capillaries.

[1]  J. N. Tilton.  The steady motion of an interface between two viscous liquids in a capillary tube , 1988 .

[2]  Raffi M. Turian,et al.  Capillary flow in a noncircular tube , 2000 .

[3]  G. Georgiou,et al.  The steady annular extrusion of a Newtonian liquid under gravity and surface tension , 2000 .

[4]  Krishnaswamy Nandakumar,et al.  A study of oil displacement on model surfaces , 1996 .

[5]  Norman R. Morrow,et al.  Capillary behavior of a perfectly wetting liquid in irregular triangular tubes , 1991 .

[6]  A. K. Chesters,et al.  Numerical investigation of the dynamic influence of the contact line region on the macroscopic meniscus shape , 1996, Journal of Fluid Mechanics.

[7]  Darrell W. Pepper,et al.  The Intermediate Finite Element Method: Fluid Flow and Heat Transfer Applications , 1999 .

[8]  Norman R. Morrow,et al.  Effect of contact angle on capillary displacement curvatures in pore throats formed by spheres , 1994 .

[9]  Nandakumar,et al.  A Model for Detachment of a Partially Wetting Drop from a Solid Surface by Shear Flow , 1997, Journal of colloid and interface science.

[10]  Dongqing Li,et al.  A model for a liquid drop spreading on a solid surface , 1998 .

[11]  Zhmud,et al.  Dynamics of Capillary Rise. , 2000, Journal of colloid and interface science.

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

[13]  James Lowndes,et al.  The numerical simulation of the steady movement of a fluid meniscus in a capillary tube , 1980, Journal of Fluid Mechanics.

[14]  John C. Slattery,et al.  Three‐dimensional, randomized, network model for two‐phase flow through porous media , 1982 .