A numerical study of the effect of insoluble surfactants on the stability of a viscous drop translating in a Hele-Shaw cell

A circular drop is a linearly stable solution for the buoyancy-driven motion of drops in a Hele-Shaw cell [Gupta et al. J. Colloid Interface Sci.218(1), 338 (1999)]. In the absence of surface-active agents, an initially prolate drop always goes to a steady circular shape while initially oblate drops exhibit complex dynamics [Gupta et al. J. Colloid Interface Sci.222, 107 (2000)]. In this study, the effect of insoluble surfactant impurities on the critical conditions for drop breakup is explored by using the Langmuir adsorption framework in conjunction with a physically based expression for the depth-averaged tangential stress exerted on a two-phase interface in a Hele-Shaw cell. It is shown that the presence of surfactants can have both a stabilizing and a destabilizing effect on the shape of the drop, depending on the Bond number, the magnitude of the initial perturbation, and the strength of surface convection. Similar to the clean drop dynamics, two marginally stable branches are found. Increasing the surface Peclet number results in the stabilization of the main branch while the secondary branch shifts to higher Bond numbers. The mode of breakup is also found to be strongly influenced by the strength of surface convection.     

Simulations of gravity-induced trapping of a deformable drop in a three-dimensional constriction

An efficient algorithm is developed to determine the three-dimensional shape of a deformable drop trapped under gravity in a constriction, employing an artificial evolution to a steady state. During the simulation, the drop surface is advanced using a rationally-devised normal "velocity", based on local deviation from the Young-Laplace equation and the adjacent solid shape, to approach the trapped drop shape. The artificial "time-dependent" evolution of the drop to the static, trapped shape requires that the free portions of the drop interface eventually satisfy the Young-Laplace equation, and the drop-solid contact portions of the drop interface conform to the solid surface. The significant advantage of this solution method is that a simple, numerically-efficient "velocity" is used to construct the evolution to the steady state; the coated areas where the drop is in near contact with solid boundaries of the constriction do not have to be specified a priori, but are found in the course of the solution. Alternative methods (e.g., boundary integral) based on realistic time-marching would be much more costly for determining the trapped state. Trapping conditions and drop shapes are studied for gravity-induced settling of a deformable drop into a three-dimensional constriction. For conditions near critical, where the trapped-drop steady state ceases to exist, severe surface-mesh distortions are treated by a combination of 'passive mesh stabilization', mesh relaxation and topological mesh transformations through node reconnections. For Bond numbers above a critical value, the drop is deformable enough to pass through the hole of the constriction, with no trapping. Critical Bond numbers are determined by linearly fitting minima of the root-mean-squared (rms) surface velocities versus corresponding Bond numbers greater than critical, and then extrapolating the Bond number to where the minimum rms velocity is zero (i.e., the drop becomes trapped). For ring and hyperbolic-tube constrictions, with axes parallel to the gravity vector, the results for trapped drops and critical Bond numbers are in close agreement with those obtained by the previous, highly-accurate axisymmetric method [1]. Also, the three-dimensional Young-Laplace and boundary-integral methods show good agreement for the static shape of a drop trapped in a tilted three-sphere constriction. For all constriction types studied, including circular rings, hyperbolic tubes and agglomerates of three and four spheres, the critical Bond number increases nearly linearly with an increase in the drop-to-hole size ratio. In contrast, the constriction type and tilt angle, which is the angle between the gravity vector and the normal to the plane of the constriction hole, have generally a weaker effect on the critical Bond number.     

Total Gaussian curvature,drop shapes and the range of applicability of drop shape techniques

Drop shape techniques are used extensively for surface tension measurement. It is well-documented that, as the drop/bubble shape becomes close to spherical, the performance of all drop shape techniques deteriorates. There have been efforts quantifying the range of applicability of drop techniques by studying the deviation of Laplacian drops from the spherical shape. A shape parameter was introduced in the literature and was modified several times to accommodate different drop constellations. However, new problems arise every time a new configuration is considered. Therefore, there is a need for a universal shape parameter applicable to pendant drops, sessile drops, liquid bridges as well as captive bubbles. In this work, the use of the total Gaussian curvature in a unified approach for the shape parameter is introduced for that purpose. The total Gaussian curvature is a dimensionless quantity that is commonly used in differential geometry and surface thermodynamics, and can be easily calculated for different Laplacian drop shapes. The new definition of the shape parameter using the total Gaussian curvature is applied here to both pendant and constrained sessile drops as an illustration. The analysis showed that the new definition is superior and reflects experimental results better than previous definitions, especially at extreme values of the Bond number.     

Critical Eotvos numbers for buoyancy-induced oil drop detachment based on shape analysis

Critical values of the Eotvos number, which is half the Bond number, above which buoyancy induced drop detachment occurs, are estimated based on force balance equations available in the literature [Colloids Surf. A: Physicochem. Eng. Aspects 178 (2001) 249]. Since there are two significantly different expressions of the capillary retention force responsible for holding oil drops on a solid substrate in an aqueous phase, the critical dimensionless number is estimated with these two distinct equations. The differential equation defining the drop shape, with the constraints of the drop volume and the 'pinned' or 'receding' contact line, is numerically solved. The equilibrium drop shapes predicted are shown to match the experimentally observed variations in drop shape. From the numerical solution, it is observed that for interfacial tension (IFT) values lower than a certain limit for a given drop size, no numerically estimated drop shape can fulfil the drop volume constraint. Similarly, for the dimensionless number above a critical value, no shape can meet all the constraints. These critical Eotvos numbers are estimated, based on the above numerical approach, for initial contact angles measured in oil varying from 20 degrees to 90 degrees. It is found that the critical Eotvos numbers estimated from the numerical shape analysis are between the critical values estimated from the two force-balance equations. Near 90 degrees, the critical values estimated from the drop shape analysis matches the values from one of the force balance estimates, but merges with the critical values of the dimensionless number, estimated from the other force balance model near 10 degrees. From this analysis, it appears that a combination of the two equations for the capillary retention force is required, with one dominating when the contact angles are high, while the other applies for low values of the contact angle.     

Liquid drops on vertical and inclined surfaces; I. An experimental study of drop geometry

Experiments have been conducted to investigate the geometric parameters necessary to describe the shapes of liquid drops on vertical and inclined plane surfaces. Two liquids and eight surfaces have been used to study contact angles, contact lines, profiles, and volumes of drops of different sizes for a range of surface conditions. The results show the contact-angle variation along the circumference of a drop to be best fit by a third-degree polynomial in the azimuthal angle. This contact-angle function is expressed in terms of the maximum and minimum contact angles of the drop, which are determined for various conditions. The maximum contact angle, thetamax, is approximately equal to the advancing contact angle, thetaA, of the liquid on the surface. As the Bond number, Bo, increases from 0 to a maximum, the minimum contact angle, thetamin, decreases almost linearly from the advancing to the receding angle. A general relation is found between thetamin/thetaA and Bo for different liquid-surface combinations. The drop contour can be described by an ellipse, with the aspect ratio increasing with Bo. These experimental results are valuable in modeling drop shape, as presented in Part II of this work.     

Retention forces and contact angles for critical liquid drops on non-horizontal surfaces

Retention forces and drop parameters are investigated for drops on the verge of sliding on vertical and inclined surfaces. Using earlier observations of drop geometry, the retentive-force factor relating surface-tension forces to contact-angle hysteresis is reliably determined. The retention force for a drop is found to be insignificantly affected by the aspect ratio of its contour. The maximum size of a drop is predicted with good accuracy, based on the two-circle method for approximating shapes of drops. The Bond number of a critical drop is found to be constant for a given surface and liquid. A general relation is proposed between the characteristic advancing and receding contact angles. The relation is supported by a large set of contact-angle data. In the absence of theta R data, the relation allows estimating the receding contact angle and the critical drop size, using only the advancing angle.     

Effect of emulsion drop-size distribution upon coalescence in simple shear flow: a population balance study

A population balance is used to examine the effect of the shape of the initial drop-size distribution of an emulsion upon its short and long-time evolution in simple shear flow. Initial distributions that are monodisperse, multidisperse, lognormal, bimodal, multimodal, and step functions are considered. At short times, it is shown that the rate of coalescence decreases by up to 25% for step distributions and up to 75% for lognormal distributions as the width of the distribution increases. Bimodal, multidisperse and multimodal distributions show intermediate decreases in the rate of coalescence, between these two values, with increases in the distribution width. Furthermore, it is found that the initial rate of coalescence is strongly dependent upon the presence of large drops. As the number fraction of large droplets within the distribution increases, the rate of coalescence also increases. At long times, all distributions move toward an asymptotic distribution shape in which the frequency of drops decreases algebraically with drop diameter at small drop diameters, and decreases exponentially with drop diameter at large drop diameters. Though portions of each distribution showed the expected asymptotic scaling behavior at long times, each asymptotic distribution nevertheless retains 'fingerprints' of the respective initial distribution. Overall, the rate of coalescence for a system is bounded by the initial rate, which is a function of the initial distribution shape, and the asymptotic rate, which is dependent upon the long-time scaling behavior. Finally, it is shown that the resolution with which the drop-size distribution of an emulsion is experimentally measured can have a significant effect upon predicted rates of coalescence.     

Development of a new methodology to study drop shape and surface tension in electric fields

Development of a new methodology for the study of both shape and surface tension of conducting drops in an electric field is presented. This methodology, called axisymmetric drop shape analysis-electric fields (ADSA-EF), generates numerical drop profiles in an electrostatic field, for a given surface tension. Then, it calculates the true value of the surface tension by matching theoretical profiles to the shape of experimental drops, using the surface tension as an adjustable parameter. ADSA-EF can be employed to simulate and study drop shapes in the electric field and to determine its effect on liquid surface tension. The method can also be used to measure surface tension in microgravity, where current drop-shape techniques are not applicable. The axisymmetric shape of the drop is the only assumption made in the development of ADSA-EF. The new scheme is applicable when both gravity and electrostatic forces are present. Preliminary measurements using ADSA-EF suggest that the surface tension of water increases by about 2% when an electric field with the magnitude of 10(6) V/m is applied.     

Surfactant-assisted spreading of a liquid drop on a smooth solid surface

Axisymmetric spreading of a liquid drop covered with an insoluble surfactant monolayer on a smooth solid substrate is numerically investigated. As the drop spreads, the adsorbed surfactant molecules are constantly redistributed along the air-liquid interface by convection and diffusion, leading to nonuniformities in surface tension along the interface. The resulting Marangoni stresses affect the spreading rate by altering the surface flow and the drop shape. In addition, surfactant accumulation in the vicinity of the moving contact line affects the spreading rate by altering the balance of line forces. Two different models for the constitutive relation at the moving contact line are used, in conjunction with a surface equation of state based on the Frumkin adsorption framework, to probe the surfactant influence. The coupled evolution equations for the drop shape and monolayer concentration profile are integrated using a pseudospectral method to determine the rate of surfactant-assisted spreading over a wide range of the dimensionless parameters governing the spreading process. The insoluble monolayer enhances spreading through two mechanisms; a reduction in the equilibrium contact angle, and an increase in the magnitude of the radial pressure gradient within the drop due to the formation of positive surface curvature near the moving contact line. Both mechanisms are driven by the accumulation of surfactant at the contact line due to surface convection. Although the Marangoni stresses induced at the air-liquid interface reduce the rate of spreading during the initial stages of spreading, their retarding effect is overwhelmed by the favorable effects of the aforementioned mechanisms to lead to an overall enhancement in the rate of spreading in most cases. The spreading rate is found to be higher for bulkier surfactants with stronger repulsive interactions. With the exception of monolayers with strong cohesive interactions which tend to retard the spreading process, the overall effect of an insoluble monolayer is to increase the rate of drop spreading. Simulation results for small Bond numbers indicate the existence of a power-law region for the time-dependence of the basal radius of the drop, consistent with experimental measurements.     

Desiccating colloidal sessile drop: dynamics of shape and concentration

Using lubrication theory, drying processes of sessile colloidal droplets on a solid substrate are studied. A simple model is proposed to describe temporal dynamics of both the shape of the drop and the volume fraction of the colloidal particles inside the drop. The concentration dependence of the viscosity is taken into account. It is shown that the final shapes of the drops depend on both the initial volume fraction of the colloidal particles and the capillary number. The results of our simulations are in a reasonable agreement with the published experimental data. Computations for the drops of aqueous solution of human serum albumin are presented.     

超疏水表面上冷凝液滴发生弹跳的机制与条件分析

使用液滴合并前后的体积和表面自由能守恒作为两个限制条件,确定了合并液滴的初始形状,即为偏离平衡态的亚稳态液滴,具有缩小其底半径而向平衡态液滴转变的推动力.进而分析了液滴变形过程中的推动力和三相线(TPCL)上的滞后阻力,建立了液滴变形的动态方程并进行了差分求解.如果液滴能够变形至底半径为0mm的状态,则根据该状态下液滴重心上移的速度确定液滴的弹跳高度.不同表面上冷凝液滴合并后的变形行为的计算结果表明,光滑表面上的液滴合并后,液滴只能发生有限的变形,一般都在达到平衡态之前就停止了变形,因此冷凝液滴不会发生弹跳;粗糙表面上的Wenzel态液滴的三相线上的滞后阻力更大,因而液滴更难以变形和弹跳;具有微纳二级结构表面上只润湿微米结构,但不润湿纳米结构的部分Wenzel态液滴能够变形至Cassie态,但没有明显的弹跳;只有在纳米或微纳二级结构表面上的较小Cassie态液滴合并后,液滴易于变形至底半径为0mm的状态并发生弹跳.因此,Cassie态合并液滴处于亚稳态,并且其三相线上的移动阻力很小,是导致冷凝液滴弹跳的关键因素.     

Microscopic treatment of a barrel drop on fibers and nanofibers

The microscopic approach of Berim and Ruckenstein (J. Phys. Chem. B 108 (2004) 19330, 19339) regarding the shape and stability of a liquid drop on a planar bare solid surface is extended to a liquid barrel drop on the bare surface of a solid cylinder (fiber) of arbitrary radius. Assuming the interaction potentials of the liquid molecules between themselves and with the molecules of the solid of the London-van der Waals form, the potential energy of a liquid molecule with an infinitely long fiber was calculated analytically. A differential equation for the drop profile was derived by the variational minimization of the total potential energy of the drop by taking into account the structuring of the liquid near the fiber. This equation was solved in quadrature and the shape and stability of the barrel drop were analyzed as functions of the radius of the fiber and the microscopic contact angle theta(0) which the drop profile makes with the surface of the fiber. The latter angle is dependent on the fiber radius and on the microscopic parameters of the model (strength of the intermolecular interactions, densities of the liquid and solid phases, hard core radii, etc.). Expressions for the evaluation of the microcontact angle from experimentally measurable characteristics of the drop profile (height, length, volume, location of inflection point) are obtained. All drop characteristics, such as stability, shape, are functions of theta(0) and a certain parameter a which depends on the model parameters. In particular, the range of drop stability consists of three domains in the plane theta(0)-a, separated by two critical curves a=a(c)(theta(0)) and a=a(c1)(theta(0)) [a(c)(theta(0))h(m1) cannot exist, whereas in the third domain (between those curves) the drop can have values of h(m) either smaller than h(m1) or larger than h(m2), where h(m2)>h(m1) is a second critical height. For sufficiently large fiber radii, R(f)1 >/= microm, the critical curves almost coincide and only two domains, the first and the second, remain. The smaller the radius, the larger is the difference between the critical curves and the larger is the second domain of drop stability. The shape of the drop depends on whether the point (theta(0),a) on the theta(0)-a plane is far from the critical curve or near it. In the first case the drop profile has generally a large circular part, while in the second case the shape is either almost planar or contains a long manchon that is similar to a film on the fiber.     

Evaluation of the surface tension measurement of axisymmetric drop shape analysis (ADSA) using a shape parameter

A quantitative criterion called “shape parameter” to evaluate the quality of surface tension measurement of Axisymmetric Drop Shape Analysis (ADSA) is presented. ADSA is a powerful technique for the measurement of interfacial tensions and contact angles of pendant drops, sessile drops, and bubbles. Despite the general success of ADSA, deficient results may be obtained for drops close to spherical shape. Therefore, the “shape parameter” was used to determine the range of drop shapes in which ADSA succeeds or fails. The “shape parameter” is a dimensionless parameter that expresses quantitatively the difference in shape between a given experimental profile and an inscribed circle. The surface tension measurements of ADSA were evaluated for both pendant drop and constrained sessile drop configurations using the shape parameter. Different shapes of the pendant drop were studied using different sizes and materials of holders. For each drop configuration, a “critical shape parameter” was defined based on the minimum value of the shape parameter that guarantees an error of less than ±0.1 mJ/m². Furthermore, the effects of the type of liquid and constellation on the “critical shape parameter” were studied.     

Polymer nanodroplets forming liquid bridges in chemically structured slit pores: a computer simulation

Using a coarse-grained bead-spring model of flexible polymer chains, the structure of a polymeric nanodroplet adsorbed on a chemically decorated flat wall is investigated by means of molecular dynamics simulation. We consider sessile drops on a lyophilic (attractive for the monomers) region of circular shape with radius R(D) while the remaining part of the substrate is lyophobic. The variation of the droplet shape, including its contact angle, with R(D) is studied, and the density profiles across these droplets also are obtained. In addition, the interaction of droplets adsorbed on two walls forming a slit pore with two lyophilic circular regions just opposite of one another is investigated, paying attention to the formation of a liquid bridge between both walls. A central result of our study is the measurement of the force between the two substrate walls at varying wall separation as well as the kinetics of droplet merging. Our results are compared to various phenomenological theories developed for liquid droplets of mesoscopic rather than nanoscopic size.     

Suspension of water droplets on individual pillars

We report results of extensive experimental and numerical studies on the suspension of water drops deposited on cylindrical pillars having circular and square cross sections and different wettabilities. In the case of circular pillars, the drop contact line is pinned to the whole edge contour until the drop collapses due to the action of gravity. In contrast, on square pillars, the drops are suspended on the four corners and spilling along the vertical walls is observed. We have also studied the ability of the two geometries to sustain drops and found that if we compare pillars with the same characteristic size, the square is more efficient in pinning large volumes, while if we normalize the volumes to pillar areas, the opposite is true.     

Modes of Nonaxisymmetry in the Stability of Fixed Contact Line Liquid Bridges and Drops

A method is presented for predicting the onset and stability character of nonaxisymmetric modes in liquid bridges and drops. The analysis applies to any fixed contact line axisymmetric interface in a steady force field. The onset and stability character of nonaxisymmetric equilibria in liquid bridges and drops is determined. Perturbation analysis is used to locate branches to nonaxisymmetry, and the configuration of the branches then gives the stability character. The number of unstable modes to both constant pressure and constant volume disturbances can be determined, so that changes in stability beyond the primary loss of stability may be examined. Although the first nonaxisymmetric mode tends to dominate higher order modes are significant for liquid bridges where length is less than radius and for drops at higher Bond numbers. At Bond numbers significantly greater than unity, the onset of the least unstable nonaxisymmetric modes tend to collapse between the fixed pressure and fixed volume axisymmetric modes of instability. For liquid bridges, two non-singular classes of nonaxisymmetric mode are distinguished: the predominant, classical shift mode; and a previously unreported tilt mode. The range over which the stability character of fixed contact line liquid bridges and drops is understood is significantly extended. Copyright 2000 Academic Press.     

Simulation analysis of contact angles and retention forces of liquid drops on inclined surfaces

A simulation study of liquid drops on inclined surfaces is performed in order to understand the evolution of drop shapes, contact angles, and retention forces with the tilt angle. The simulations are made by means of a method recently developed for dealing with contact angle hysteresis in the public-domain Surface Evolver software. The results of our simulations are highly dependent on the initial contact angle of the drop. For a drop with an initial contact angle equal to the advancing angle, we obtain results similar to those of experiments in which a drop is placed on a horizontal surface that is slowly tilted. For drops with an initial contact angle equal to the mean between the advancing and the receding contact angles, we recover previous results of finite element studies of drops on inclined surfaces. Comparison with experimental results for molten Sn-Ag-Cu on a tilted Cu substrate shows excellent agreement.     

Partial wetting of chemically patterned surfaces: the effect of drop size

Partial wetting of chemically heterogeneous substrates is simulated. Three-dimensional sessile drops in equilibrium with smooth surfaces supporting ordered chemical patterns are considered. Significant features are observed as a result of changing the drop volume. The number of equilibrated drops is found either to remain constant or to increase with growing drop volume. The shape of larger drops appears to approach that of a spherical cap and their three-phase contact line seems, on a larger scale, more circular in shape than that of smaller drops. In addition, as the volume is increased, the average contact angle of drops whose free energy is lowest among all equilibrium-shaped drops of the same volume appears to approach the angle predicted by Cassie. Finally, contrary to results obtained with two-dimensional drops, contact angle hysteresis observed in this system is shown to exhibit a degree of volume dependence in the advancing and receding angles. Qualitative differences in the wetting behavior associated with the two different chemical patterns considered here, as well as differences between results obtained with two-dimensional and three-dimensional drops, can possibly be attributed to variations in the level of constraint imposed on the drop by the different patterns and by the dimensionality of the system.     

Non-aqueous phase liquid drop formation within a water saturated fracture

This work focuses on the mechanisms of non-aqueous phase liquid (NAPL) drop formation within a single fracture fed from a NAPL reservoir by way of a circular orifice, such as a pore. The fracture is assumed to be fully saturated, the relative wettability of the system is assumed water-wet, and the water velocity profile within the fracture is described by a Poiseuille flow. The size of the NAPL drops is investigated for various water flow velocities and NAPL entrance diameters. A force balancing method was used to determine the radii of detached drops. The drop sizes calculated from the model developed here are shown to be in agreement with available experimental drop size data. It is shown that at low Reynolds numbers the buoyancy force is the dominant force acting on the drop during the formation process and at high Reynolds numbers the viscous forces dominate. A simplified expression relating the geometry of the fractured system to the drop radii is developed from the model equations, and it is shown to predict drop radii that match well with both the model simulations and the available experimental data.     

超疏水表面微纳二级结构对冷凝液滴最终状态的影响

从超疏水表面(SHS)上初始冷凝液核长大、合并、形成初始液斑开始,分析计算了冷凝液斑变形成为Wenzel或Cassie液滴过程中界面能量的变化,并以界面能曲线降低、是否取最小值为判据,确定冷凝液滴的最终稳定状态.计算结果表明:在只有微米尺度的粗糙结构表面上,冷凝液滴的界面能曲线一般都是先降低再升高,呈现Wenzel状态;而当表面具有微纳米二级粗糙结构,且纳米结构的表面空气面积分率较高时,冷凝液滴的能量曲线持续降低,直至界面能最小的Cassie状态,因此可以自发地形成Cassie液滴.还计算了文献中具有不同结构参数的SHS上冷凝液滴的状态和接触角,并与实验结果进行了比较,结果表明,计算的冷凝液滴状态与实验观察结果完全吻合.因此,微纳二级结构是保持冷凝液滴在SHS上呈现Cassie状态的重要因素.