Shape analysis based critical Eotvos numbers for buoyancy induced partial detachment of oil drops from hydrophilic surfaces

Removal of oil drops from solid surfaces immersed in an aqueous medium is of interest in many applications. It has been shown that drop shape analysis can be used to predict conditions at which the stability limit of a lighter than water oil drop on a solid surface immersed in an aqueous bath is reached (Adv. Colloid Interface Sci. 98 (2002) 265). However the above analysis is restricted to cases where the contact angle made by the drop is below 90degrees and when the surface conditions result in a 'pinned' contact line. In this paper, it is shown that drop shape analysis can be used to predict the critical conditions at which drop stability limit is reached for drop contact angles of 90degrees and above, which is encountered with 'hydrophilic' surfaces. This critical condition can predict the occurrence of partial oil drop detachment, before complete removal due to 'roll-up', which occurs when the hydrophilic surface is adequately smooth which prevents 'pinning' of the contact line. The critical conditions at which partial drop detachment occurs can also be approximately predicted from simple force balances. It has been shown (Adv. Colloid Interface Sci. 98 (2002) 265) that for contact angles less than 90degrees, the critical limit based on shape analysis appears to resolve the differences that arise due to alternate expressions for capillary retention force. This paper shows that even for contact angles above 90degrees, the critical conditions predicted from the shape analysis resolves the differences in the predictions from the alternate force balances. Drop shape analysis used in this paper is based on the 'Arc-length' form of Young-Laplace or 'drop shape' equation, which is different from the 'Y vs X' form of the above equation that is used in Adv. Colloid Interface Sci. 98 (2002) 265. The above drop shape equation is solved by a fourth order Runge-Kutta technique and it is shown that for angles less than 90degrees, the two forms of the drop shape equation, predict almost identical values of the critical Eotvos number. This paper highlights the competing effects of interfacial tension lowering induced drop instability and 'roll-up', a term that is used to describe the retraction of the contact line of an oil drop on a surface, in being the primary c ause for drop detachment.     

Limiting conditions for applying the spherical section assumption in contact angle estimation

The shape of liquid drops on solid surfaces deviates from the spherical as tension decreases and gravity effects start affecting the drop shape. This paper attempts to define this deviation and estimates the dimensionless Eotvos number limits above which the deviation becomes "significant." The use of these limiting values can facilitate estimation of contact angle in the following manner. It is well known that the equilibrium contact angle made by a liquid drop on a solid surface can be estimated from measurements of two drop parameters. These parameters can be any two chosen from the drop volume, height, and wetted radius. In case the effect of gravity on the drop shape is negligible, simple algebraic relations derived from the spherical section assumption exist, from which the contact angle can be estimated. In systems where the "spherical section" assumption is invalid, the Laplace equation for the drop shape has been solved numerically with any two of the above parameters as the constraints, to obtain the contact angle. In this paper, Eotvos numbers at which the deviation of the drop profile from the spherical is significant enough to result in contact angle deviation of 1 degrees are estimated. The limiting values of Eotvos number, expressed as a function of the original contact angle made by the spherical profile, are obtained by solving the Laplace equation for the drop shape with the drop volume and wetted radius constraints for decreasing values of Interfacial tension. These limiting values are also estimated for different drop sizes and for cases where the drop phase is heavier (sessile) and lighter (buoyant) than the surrounding fluid. The independence of the Eotvos number estimates from the sign of the density difference as well as the drop size is shown. These Eotvos number limits can be used to check if the spherical section assumption, with the resulting simple algebraic relations, can be used for contact angle estimation and other shape-related analysis for a system.     

Wetting: Inverse Dynamic Problem and Equations for Microscopic Parameters

Movement of a liquid meniscus in a low-diameter capillary while it is being filled or emptied is considered. The liquid is nonvolatile. Assuming low Reynolds number and low capillary number, the liquid-gas interface shape is studied. Angles of inclination of this boundary to the solid near the contact line are small. Consideration is given to the inverse problem in wetting dynamics: to establish an analytic expression for the universal constant that influences the dynamics of a three-phase contact line. Inverse relations for microscopic parameters in terms of macroscopic measured values obtained in experiments with a meniscus moving through a capillary are derived. The inverse relations are substantiated independently. To do so, numerical experiments for a van der Waals liquid have been carried out, using the de Gennes model of partial wetting. General formulas for microparameters agree well with numerical experiments. The article provides the similarity criterion which influences the wetting in the case of a van der Waals liquid meniscus. The inverse dynamic problem for both an advancing and a receding meniscus is solved. A relation for the critical speed of meniscus recession is proposed. Two contact angles for a meniscus are discussed. Behavior of dynamic contact angles in the vicinity of the critical speed is studied. One of the angles is shown to vanish at less than the critical speed, and the other one, exactly at the critical speed. In the case of an advancing meniscus the equations for microparameters are valid for both partial and complete wetting. The proposed inverse expression for complete wetting allows determination of the maximum precursor film thickness and its dependence on the motion speed (also determination of the Hamaker constant in the case of a van der Waals liquid). Copyright 2000 Academic Press.     

Capillary forces between chemically different substrates

Motivated by experimental results, we present numerical and analytical calculations of the capillary force exerted by a capillary bridge spanning the gap between two parallel flat plates of asymmetric wettability. Depending on whether the sum of the two contact angles is smaller or larger than 180 degrees, the capillary force is either attractive or repulsive at small separations D between the plates. In either cases the magnitude of the force diverges as D approaches zero. The leading order of this divergence is captured by an analytical expression deduced from the geometry of the meniscus of a flat capillary bridge. The results for substrates with different wettability reveal an interesting behavior: with the sum of the contact angles fixed, the magnitude of the capillary force and the rupture separation decreases as the asymmetry in contact angles is increased. In addition, we present the rupture separation, i.e., the maximal extension of a capillary bridge, as a function of the contact angles. Our results provide an extensive picture of surface wettability effects on capillary adhesion.     

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.     

An analytical solution for determination of small contact angles from sessile drops of arbitrary size

An analytical solution to the capillary equation of Young and Laplace is derived that allows determination of the static contact angle based on the volume of a sessile drop and the wetted area of the substrate. This solution does not require numerical integration to determine the drop profile and accounts for surface deformation due to gravitational effects. Calculation of the static contact angle by this method is remarkably simple and accurate when the contact angle is less than 30 degrees. A natural scaling arises in the solution, which provides indication of when a drop is small enough so as to neglect gravitational influences on the surface shape which, for small contact angles, is generally less than 1 microl. The technique described has the simplicity of the spherical cap approximation but remains accurate for any size of sessile drop.     

Sliding behavior of liquid droplets on tilted Langmuir-Blodgett surfaces

The sliding behavior of liquid droplets on inclined Langmuir-Blodgett surfaces was investigated. The critical sliding angle defined as the tilt angle of the surface at which the drop slides down as well as the advancing and receding contact angles was measured for five different liquids on five surfaces. In addition, the contact line geometry was analyzed at critical sliding angle. The experimental relationship between the surface tension forces resulting from contact angle hysteresis and the weight of the drop was compared to theoretical predictions. Even though the shape of the drop bases was found as skewed ellipses, a model assuming parallel-sided elongated drops is shown to describe reasonably the experimental values. This result probably indicates the main influence of the capillary forces at the rear and front edges of the drop with respect to that exerted on the lateral sides.     

Assessing the accuracy of contact angle measurements for sessile drops on liquid-repellent surfaces

Gravity-induced sagging can amplify variations in goniometric measurements of the contact angles of sessile drops on super-liquid-repellent surfaces. The very large value of the effective contact angle leads to increased optical noise in the drop profile near the solid-liquid free surface and the progressive failure of simple geometric approximations. We demonstrate a systematic approach to determining the effective contact angle of drops on super-repellent surfaces. We use a perturbation solution of the Bashforth-Adams equation to estimate the contact angles of sessile drops of water, ethylene glycol, and diiodomethane on an omniphobic surface using direct measurements of the maximum drop width and height. The results and analysis can be represented in terms of a dimensionless Bond number that depends on the maximum drop width and the capillary length of the liquid to quantify the extent of gravity-induced sagging. Finally, we illustrate the inherent sensitivity of goniometric contact angle measurement techniques to drop dimensions as the apparent contact angle approaches 180°.     

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.     

The Neumann and Young equations for nematic contact lines

The Neumann and Young equations for three-phase nematic contact lines have been derived using the momentum balance equation and classical liquid crystal physics theories. The novel finding is the presence of bending forces, originating from the anchoring energy of nematic interfaces, and acting on the contact line. The classical Neumann triangle or tensile force balance becomes in the presence of a nematic phase the Neumann pentagon, involving the usual three tensile forces and two additional bending forces. The Young equation that describes the static contact angle of a fluid in contact with a rigid solid is again a tensile force balance along the solid, but for nematics it also involves an additional bending force. The effects of the bending forces on contact angles and wetting properties of nematic liquid crystals are thoroughly characterized. It is found that in terms of the spreading coefficient, bending forces enlarge the partial wetting window that exists between dewetting and spontaneous spreading. Bending forces also affect the behaviour of the contact angle, such that spreading occurs at contact angles greater than zero and dewetting at values greater than pi. Finally, the contact angle range in the partial wetting regime is always less than pi.     

Drop shape visualization and contact angle measurement on curved surfaces

The shape and contact angles of drops on curved surfaces is experimentally investigated. Image processing, spline fitting and numerical integration are used to extract the drop contour in a number of cross-sections. The three-dimensional surfaces which describe the surface-air and drop-air interfaces can be visualized and a simple procedure to determine the equilibrium contact angle starting from measurements on curved surfaces is proposed. Contact angles on flat surfaces serve as a reference term and a procedure to measure them is proposed. Such procedure is not as accurate as the axisymmetric drop shape analysis algorithms, but it has the advantage of requiring only a side view of the drop-surface couple and no further information. It can therefore be used also for fluids with unknown surface tension and there is no need to measure the drop volume. Examples of application of the proposed techniques for distilled water drops on gemstones confirm that they can be useful for drop shape analysis and contact angle measurement on three-dimensional sculptured surfaces.     

The Neumann and Young equations for nematic contact lines

The Neumann and Young equations for three-phase nematic contact lines have been derived using the momentum balance equation and classical liquid crystal physics theories. The novel finding is the presence of bending forces, originating from the anchoring energy of nematic interfaces, and acting on the contact line. The classical Neumann triangle or tensile force balance becomes in the presence of a nematic phase the Neumann pentagon, involving the usual three tensile forces and two additional bending forces. The Young equation that describes the static contact angle of a fluid in contact with a rigid solid is again a tensile force balance along the solid, but for nematics it also involves an additional bending force. The effects of the bending forces on contact angles and wetting properties of nematic liquid crystals are thoroughly characterized. It is found that in terms of the spreading coefficient, bending forces enlarge the partial wetting window that exists between dewetting and spontaneous spreading. Bending forces also affect the behaviour of the contact angle, such that spreading occurs at contact angles greater than zero and dewetting at values greater than pi. Finally, the contact angle range in the partial wetting regime is always less than pi.     

Computation of constant mean curvature surfaces: Application to the gas-liquid interface of a pressurized fluid on a superhydrophobic surface

The interface shape separating a gas layer within a superhydrophobic surface consisting of a square lattice of posts from a pressurized liquid above the surface is computed numerically. The interface shape is described by a constant mean curvature surface that satisfies the Young-Laplace equation with the three-phase gas-liquid-solid contact line assumed pinned at the post outer edge. The numerical method predicts the existence of constant mean curvature solutions from the planar, zero curvature solution up to a maximum curvature that is dependent on the post shape, size and pitch. An overall force balance between surface tension and pressure forces acting on the interface yields predictions for the maximum curvature that agree with the numerical simulations to within one percent for convex shapes such as circular and square posts, but significantly over predicts the maximum curvature for non-convex shapes such as a circular post with a sinusoidal surface perturbation. Changing the post shape to increase the contact line length, while maintaining constant post area, results in increases of 2 to 12% in the maximum computable curvature for contact line length increases of 11 to 77%. Comparisons are made to several experimental studies for interface shape and pressure stability.     

Effect of Uniform Magnetic Field on the Stability of the Rayleigh-Benard Convection in a Binary Electrolyte: A Theoretical Analysis

The Rayleigh-Benard stability in an electrochemical cell with a binary electrolyte in a uniform magnetic field is studied theoretically. The mathematical model is based on equations describing the motion of incompressible viscous liquid in the Boussinesq approximation, which account for the Archimedean buoyancy force and the Lorentz force, and equations of material balance with the electroneutrality condition. The set of equations for perturbation amplitudes is solved analytically by the Galerkin method and numerically, by the Runge-Kutta method. Dependences of the critical value of the Rayleigh number on the Hartmann number, which expresses the ratio between magnetic and viscous forces, are obtained for monotonic instability of limiting current at liquid and solid electrodes. The effect the variations in the electrolyte's electroconductance, which are caused by concentration variations, exert on the instability emergence is analyzed.     

Effect of Insoluble Surfactants on Drainage and Rupture of a Film between Drops Interacting under a Constant Force

The deformation, drainage, and rupture of an axisymmetrical film between colliding drops in the presence of insoluble surfactants under the influence of van der Waals forces is studied numerically at small capillary and Reynolds numbers and small surfactant concentrations. Constant-force collisions of Newtonian drops in another Newtonian fluid are considered. The mathematical model is based on the lubrication equations in the gap between drops and the creeping flow approximation of Navier–Stokes equations in the drops, coupled with velocity and stress boundary conditions at the interfaces. A nonuniform surfactant concentration on the interfaces, governed by a convection–diffusion equation, leads to a gradient of the interfacial tension which in turn leads to additional tangential stress on the interfaces (Marangoni effects). The mathematical problem is solved by a finite-difference method on a nonuniform mesh at the interfaces and a boundary-integral method in the drops. The whole range of the dispersed to continuous-phase viscosity ratios is investigated for a range of values of the dimensionless surfactant concentration, Peclét number, and dimensionless Hamaker constant (covering both “nose” and “rim” rupture). In the limit of the large Peclét number and the small dimensionless Hamaker constant (characteristic of drops in the millimeter size range) a fair approximation to the results is provided by a simple expression for the critical surfactant concentration, drainage being virtually uninfluenced by the surfactant for concentrations below the critical surfactant concentration and corresponding to that for immobile interfaces for concentrations above it.     

Contact angles of liquid drops on super hydrophobic surfaces: understanding the role of flattening of drops by gravity

Measurement of contact angles on super hydrophobic surfaces by conventional methods can produce ambiguous results. Experimental difficulties in constructing tangent lines, gravitational distortion or erroneous assumptions regarding the extent of spreading can lead to underestimation of contact angles. Three models were used to estimate drop shape and perceived contact angles on completely nonwetting super hydrophobic surfaces. One of the models employed the classic numerical solutions from Bashforth and Adams. Additionally, two approximate models were derived as part of this work. All three showed significant distortion of microliter-sized drops and similar trends in perceived contact angles. Liquid drops of several microliters are traditionally used in sessile contact angle measurements. Drops of this size are expected to and indeed undergo significant flattening on super hydrophobic surfaces, even if the wetting interactions are minimal. The distortion is more pronounced if the liquid has a lesser surface tension or greater density. For surfaces that are completely nonwetting, underestimation of contact angles can be tens of degrees. Our modeling efforts suggest that accurate contact angle measurements on super hydrophobic surfaces would require very small sessile drops, on the order of hundreds of picoliters.     

Indirect methods to measure wetting and contact angles on spherical convex and concave surfaces

In this work, a method was developed for indirectly estimating contact angles of sessile liquid drops on convex and concave surfaces. Assuming that drops were sufficiently small that no gravitational distortion occurred, equations were derived to compute intrinsic contact angles from the radius of curvature of the solid surface, the volume of the liquid drop, and its contact diameter. These expressions were tested against experimental data for various liquids on polytetrafluoroethylene (PTFE) and polycarbonate (PC) in the form of flat surfaces, spheres, and concave cavities. Intrinsic contact angles estimated indirectly using dimensions and volumes generally agreed with the values measured directly from flat surfaces using the traditional tangent method.     

Microscopic observation and characterization of the oil bridge between dehulled-roasted sesame seeds

The formation of an oil bridge between adhesive seeds was observed microscopically. The geometry of the oil bridge was affected by the shape of the adhesive seeds. The capillary force of the oil bridge was estimated from the image captured by the microscope. The average capillary force was 127 μN, which was five times higher than the average gravity of the seeds. It was observed that several oil bridges formed between two seeds. These results indicated an adequate ability of the seeds to adhere. The capillary force of the oil bridge increased with surface oil content. The probability of formation of an oil bridge increased with surface oil content when the surface oil content was above 0.63%. The probability of formation of an oil bridge markedly increased when sucrose was added to the seeds.     

Spreading of liquid drops over dry porous layers: complete wetting case

Spreading of small liquid drops over thin dry porous layers is investigated from both theoretical and experimental points of view. Drop motion over a porous layer is caused by an interplay of two processes: (a) the spreading of the drop over already saturated parts of the porous layer, which results in an expanding of the drop base; (b) the imbibition of the liquid from the drop into the porous substrate, which results in a shrinkage of the drop base and an expanding of the wetted region inside the porous layer. As a result of these two competing processes, the radius of the drop goes through a maximum value over time. A system of two differential equations is derived to describe the evolution with time of radii of both the drop base and the wetted region inside the porous layer. This system includes two parameters: one accounts for the effective lubrication coefficient of the liquid over the wetted porous substrate and the other is a combination of permeability and effective capillary pressure inside the porous layer. Two additional experiments are used for an independent determination of these two parameters. The system of differential equations does not include any fitting parameter after these two parameters are determined. Experiments were carried out on the spreading of silicone oil drops over various dry microfiltration membranes (permeable in both normal and tangential directions). The time evolution of the radii of both the drop base and the wetted region inside the porous layer are monitored. All experimental data fell on two universal curves if appropriate scales are used with a plot of the dimensionless radii of the drop base and of the wetted region inside the porous layer on dimensionless time. The predicted theoretical relationships are two universal curves accounting quite satisfactorily for the experimental data. According to our theory prediction, (i) the dynamic contact angle dependence on the same dimensionless time as before should be a universal function and (ii) the dynamic contact angle should change rapidly over an initial short stage of spreading and should remain a constant value over the duration of the rest of the spreading process. The constancy of the contact angle on this stage has nothing to do with hysteresis of the contact angle: there is no hysteresis in our system. These conclusions again are in good agreement with our experimental observations.     

Analysis of the equilibrium droplet shape based on an ellipsoidal droplet model

The extent of a droplet's spreading over a flat, smooth solid substrate and its equilibrium height in the presence of gravity are determined approximately, without a numerical solution of the governing nonlinear differential equation, by assuming that the droplet takes on the shape of an oblate spheroidal cap and by minimizing the corresponding free energy. The comparison with the full numerical evaluations confirms that the introduced approximation and the obtained results are accurate for contact angles below about 120° and for droplet sizes on the order of the capillary length of the liquid. The flattening effect of gravity is to increase the contact radius and decrease the height of the droplet, with these being more pronounced for higher values of the Bond number.