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.     

Stable drop shapes under disjoining pressure. II. Multiplicity and stability

The stability of the contact line region as affected by the disjoining pressure has been analyzed by solving the augmented Young-Laplace equation. Because of the results in Part I (Zhang, X., Neogi, P., and Ybarra, R. M., J. Colloid Interface Sci.), we have concentrated on obtaining multiple solutions for the same set of conditions. As many as five solutions were obtained: drops that end in a thin film with uniform thickness and where the film shape oscillates, drops that end with microscopic contact angles, as well as uniform thin films of two different thicknesses. The results of linear stability analysis were used to show that most cases were unstable to infinitesimal disturbances. Only two stable drop shapes for the particular disjoining pressure investigated are stable, a thin film of constant thickness and a thin drop that ends in a film of same thickness. Both multiplicity and stability have been discussed here for the first time and shed considerable light on the role of the attractive and repulsive forces.     

Polymeric stabilized emulsions: steric effects and deformation in soft systems

Polymeric stabilizers are used in a broad range of processes and products, from pharmaceuticals and engine lubricants to formulated foods and shampoos. In rigid particulate systems, the stabilization mechanism is attributed to the repulsive force that arises from the compression of the polymer coating or "steric brush" on the interacting particles. This mechanism has dictated polymer design and selection for more than thirty years. Here we show, through direct measurement of the repulsive interactions between immobilized drops with adsorbed polymers layers in aqueous electrolyte solutions, that the interaction is a result of both steric stabilization and drop deformation. Drops driven together at slow collision speeds, where hydrodynamic drainage effects are negligible, show a strong dependence on drop deformation instead of brush compression. When drops are driven together at higher collision speeds where hydrodynamic drainage affects the interaction force, simple continuum modeling suggests that the film drainage is sensitive to flow through the polymer brush. These data suggest, for drop sizes where drop deformation is appreciable, that the stability of emulsion drops is less sensitive to the molecular weight or size of the adsorbed polymer layer than for rigid particulate systems.     

Film drainage and interfacial instabilities in polymeric systems with diffuse interfaces

We report an experimental investigation on the effect of mutual diffusion in polymeric systems on film drainage between two captive drops. The main objective is to study the influence of diffuse interfaces on film drainage. This is done by using material combinations with different interfacial properties and interferometric visualization of the film between two interacting drops. For highly diffusive systems film drainage is observed to be, in contrast to immiscible systems, non-axisymmetric and unstable immediately after the film formation (at a few micrometers film thickness). Depending on whether the total thickness of the diffusion layers in the film is smaller or larger than the thickness of the film, Marangoni convection is found to enhance or delay film drainage. Enhanced film drainage is determined to be in order of 100 times faster than predicted by the current models, while delayed film drainage is observed after a drainage period where experimental and predicted results (assuming, a partially mobile interface) are in close agreement.     

Theory of non-equilibrium force measurements involving deformable drops and bubbles

Over the past decade, direct force measurements using the Atomic Force Microscope (AFM) have been extended to study non-equilibrium interactions. Perhaps the more scientifically interesting and technically challenging of such studies involved deformable drops and bubbles in relative motion. The scientific interest stems from the rich complexity that arises from the combination of separation dependent surface forces such as Van der Waals, electrical double layer and steric interactions with velocity dependent forces from hydrodynamic interactions. Moreover the effects of these forces also depend on the deformations of the surfaces of the drops and bubbles that alter local conditions on the nanometer scale, with deformations that can extend over micrometers. Because of incompressibility, effects of such deformations are strongly influenced by small changes of the sizes of the drops and bubbles that may be in the millimeter range. Our focus is on interactions between emulsion drops and bubbles at around 100 μm size range. At the typical velocities in dynamic force measurements with the AFM which span the range of Brownian velocities of such emulsions, the ratio of hydrodynamic force to surface tension force, as characterized by the capillary number, is ~ 10^{− 6} or smaller, which poses challenges to modeling using direct numerical simulations. However, the qualitative and quantitative features of the dynamic forces between interacting drops and bubbles are sensitive to the detailed space and time-dependent deformations. It is this dynamic coupling between forces and deformations that requires a detailed quantitative theoretical framework to help interpret experimental measurements. Theories that do not treat forces and deformations in a consistent way simply will not have much predictive power. The technical challenges of undertaking force measurements are substantial. These range from generating drop and bubble of the appropriate size range to controlling the physicochemical environment to finding the optimal and quantifiable way to place and secure the drops and bubbles in the AFM to make reproducible measurements. It is perhaps no surprise that it is only recently that direct measurements of non-equilibrium forces between two drops or two bubbles colliding in a controlled manner have been possible. This review covers the development of a consistent theory to describe non-equilibrium force measurements involving deformable drops and bubbles. Predictions of this model are also tested on dynamic film drainage experiments involving deformable drops and bubbles that use very different techniques to the AFM to demonstrate that it is capable of providing accurate quantitative predictions of both dynamic forces and dynamic deformations. In the low capillary number regime of interest, we observe that the dynamic behavior of all experimental results reviewed here are consistent with the tangentially immobile hydrodynamic boundary condition at liquid–liquid or liquid–gas interfaces. The most likely explanation for this observation is the presence of trace amounts of surface-active species that are responsible for arresting interfacial flow.     

Unconventional spinodal surface fluctuations on polymer films

We study the temporal growth pattern of surface fluctuations on a series of spinodally unstable polymer films where the instability can be adjusted with the film thickness, h0. For the most unstable film studied (whose /(h0 - h(sp))/h(sp)/ = 0.988; h(sp) is the thickness where the second derivative of the interfacial potential of the film equals zero), the growth rate function of the surface modes as a function of the wavevector fits well to the mean-field theory. When the film thickness is increased such that /(h0 - h(sp))/h(sp)/ < or = 0.977, the mean-field theory demonstrates marked disagreement with experiment, notwithstanding the provision of the known corrections from nonlinear effects and thermal noise. We show that the deviations arise from large-amplitude fluctuations induced by homogeneous nucleation, which are not considered in the conventional treatments.     

Hydrodynamic instability of a thin viscous film between two drops

Linear stability analysis is employed to derive analytical expressions for the growth rate of disturbances applied to a thin plane-parallel film trapped between two drops. From these expressions, the band of unstable wavenumbers and the "most dangerous" wavenumber are identified for systems in the absence and presence of insoluble surfactant. Marangoni effects are shown to exert a stabilizing influence and reasonably good agreement with experimental observations is found. Subsequent nonlinear analysis indicated amplification of the disturbance growth rate beyond that suggested by linear theory as the film proceeded toward rupture.     

Hydrodynamic interaction between spheres coated with deformable thin liquid films

In this article, we considered the hydrodynamic interaction between two unequal spheres coated with thin deformable liquids in the asymptotic lubrication regime. This problem is a prototype model for drop coalescence through the so-called "film drainage" mechanism, in which the hydrodynamic contribution comes dominantly from the lubrication region apart from the van der Waals interaction force. First, a general formulation was derived for two unequal coated spheres that experienced a head-to-head collision at a very close proximity. The resulting set of the evolution equations for the deforming film shapes and stress distributions was solved numerically. The film shapes and hydrodynamic interaction forces were determined as functions of the separation distance, film thickness, viscosity ratios, and capillary numbers. The results show that as the two spheres approach each other, the films begin to flatten and eventually to form negative curvature (or a broad dimple) at their forehead areas in which high lubrication pressure is formed. The dimple formation occurs earlier as the capillary number increases. For large capillary numbers, the film liquids are drained out from their forehead areas and the coated liquid films rupture before the two films "touch" each other. Meanwhile, for small capillary numbers, the gap liquid is drained out first and the two liquid films eventually coalesce.     

Detachment force of particles from air-liquid interfaces of films and bubbles

The detachment force required to pull a microparticle from an air-liquid interface is measured using atomic force microscopy (AFM) and the colloidal probe technique. Water, solutions of sodium dodecyl sulfate (SDS), and silicone oils are tested in order to study the effects of surface tension and viscosity. Two different liquid geometries are considered: the air-liquid interface of a bubble and a liquid film on a solid substrate. It was shown that detaching particles from liquid films is fundamentally different than from bubbles or drops due to the restricted flow of the liquid phase. Additional force is required to detach a particle from a film, and the maximum force during detachment is not necessarily at the position where the particle breaks away from the interface (as seen in bubble or drop systems). This is due to the dynamics of meniscus formation and viscous effects, which must be considered if the liquid is constrained in a film. The magnitude of these effects is related to the liquid viscosity, film thickness, and detachment speed.     

Effects of bubble size and approach velocity on bubble–particle interaction – a theoretical study

A theoretical analysis of the atomic force microscopy (AFM) approach–retract dynamic interaction between an air bubble and a hydrophilic silica plane was carried out based on the well-established Stokes–Reynolds–Young–Laplace model. An air bubble with different radii attached to the end of a cantilever approached the silica surface with different approach velocities in a 10^?3?M KCl solution. Results showed that with increasing approach velocity (0.1, 1, and 10?µm/s), the repulsive force, flattened area of the film, and hydrodynamic suction force between the 100-µm bubble and the silica plane increased. The film continued thinning at the initial stages of bubble retraction because of the attractive hydrodynamic pressure. When the bubble size decreased, the influence of hydrodynamic pressure was less evident. The final film thickness before bubble retraction was similar to the theoretical equilibrium thickness when the Laplace pressure was equal to the disjoining pressure.     

Spinodal instability and pattern formation in thin liquid films confined between two plates

The instability, morphology and pattern formation engendered by the van der Waals force in a thin liquid film of thickness h confined between two closely placed solid surfaces (at distance d > h) are investigated based on nonlinear 3D simulations. The initial and the final stages of dewetting and pattern formation are found to be crucially dependent on the volumetric (thickness) ratio of air and liquid and its deviation from the location of the maximum of the spinodal parameter versus volumetric ratio curve. On a low energy surface, relatively thinner films and wider air gaps favor initial dewetting of the lower plate by the formation of holes, whereas thicker films with thinner air gaps initially evolve by the formation of columns/bridges that join the upper plate. In the later stage of evolution, the initial holes in thinner films evolve into columns/drops, while a rapid coalescence of columns in the thicker films eventually causes formation of holes. Thus, a phase inversion, either from liquid-in-air to air-in-liquid dispersion or vice versa, occurs during the final stages of evolution. A thin film confined between two high-energy solid surfaces forms columns (bridges) only when its mean thickness, h0, is greater than a critical thickness (hc) or the air gap is smaller than a critical distance. The patterns can be aligned by using a topographically patterned confining surface. Conditions on pattern periodicity, amplitude, and the volumetric ratio of air and liquid in the gap are explored for the formation of various types of ordered patterns including annular rings of columns, concentric ripples, parallel channels and a rectangular array of complex features. The results are of significance in soft lithographies such as LISA, soft stamping and capillary force lithography.     

Stability of draining plane-parallel films containing surfactants

The stability of partially mobile draining thin liquid films with respect to axisymmetric fluctuations was studied. The material properties of the interfaces (Gibbs elasticity, surface and bulk diffusions) were taken into account. When studying the long wave stability of films, the coupling between the drainage and perturbation flows was considered and the lubrication approximation was applied. Two types of wave modes were examined: radially-bounded and unbounded waves. The difference between the thickness of loss of stability, h(st), the transitional thickness, h(tr), at which the critical wave causing rupture becomes unstable, and the critical thickness, h(cr), when the film ruptures, is demonstrated. Both the linear and the non-linear theories give h(st) > h(tr) > h(cr). The numerical results show that the interfacial mobility does not significantly influence the thickness of the draining film rupture. The interfacial tension and the disjoining pressure are the major factors controlling the critical thickness. The available experimental data for critical thicknesses of foam and emulsion films show excellent agreement with the theoretical predictions. The important role of the electromagnetic retardation term in the van der Waals interaction is demonstrated. Other published theories of the film stability are discussed.     

Instability and dewetting of ultrathin solid viscoelastic films on homogeneous and heterogeneous substrates

Instability and dewetting engendered by the van der Waals force in soft thin (<100 nm) linear viscoelastic solid (e.g., elastomeric gel) films on uniform and patterned surfaces are explored. Linear stability analysis shows that, although the elasticity of the film controls the onset of instability and the corresponding critical wavelength, the dominant length-scale remains invariant with the elastic modulus of the film. The unstable modes are found to be long-wave, for which a nonlinear long-wave analysis and simulations are performed to uncover the dynamics and morphology of dewetting. The stored elastic energy slows down the temporal growth of instability significantly. The simulations also show that a thermodynamically stable film with zero-frequency elasticity can be made unstable in the presence of physico-chemical defects on the substrate and can follow an entirely different pathway with far fewer holes as compared to the viscous films. Further, the elastic restoring force can retard the growth of a depression adjacent to the hole-rim and thus suppress the formation of satellite holes bordering the primary holes. These findings are in contrast to the dewetting of viscoelastic liquid films where nonzero frequency elasticity accelerates the film rupture and promotes the secondary instabilities. Thus, the zero-frequency elasticity can play a major role in imposing a better-defined long-range order to the dewetted structures by arresting the secondary instabilities.     

Contact instability of thin elastic films on patterned substrates

The free surface of a soft elastic film becomes unstable and forms an isotropic labyrinth pattern when a rigid flat plate is brought into adhesive contact with the film. These patterns have a characteristic wavelength, lambda approximately 3H, where H is the film thickness. We show that these random structures can be ordered, modulated, and aligned by depositing the elastic film (cross-linked polydimethylsiloxane) on a patterned substrate and by bringing the free surface of the film in increasing adhesive contact with a flat stamp. Interestingly, the influence of the substrate "bleeds" through the film to its free surface. It becomes possible to generate complex two-dimensional ordered structures such as an array of femtoliter beakers even by using a simple one-dimensional stripe patterned substrate when the instability wavelength, lambda approximately 3H, nearly matches the substrate pattern periodicity. The free surface morphology is modulated in situ by merely varying the stamp-surface separation distance. The free surface structures originating from the elastic contact instability can also be made permanent by the UV-ozone induced oxidation and stiffening.     

Nonflat equilibrium liquid shapes on flat surfaces

The hydrostatic pressure in thin liquid layers differs from the pressure in the ambient air. This difference is caused by the actions of surface forces and capillary pressure. The manifestation of the surface force action is the disjoining pressure, which has a very special S-shaped form in the case of partial wetting (aqueous thin films and thin films of aqueous electrolyte and surfactant solutions, both free films and films on solid substrates). In thin flat liquid films the disjoining pressure acts alone and determines their thickness. However, if the film surface is curved then both the disjoining and the capillary pressures act simultaneously. In the case of partial wetting their simultaneous action results in the existence of nonflat equilibrium liquid shapes. It is shown that in the case of S-shaped disjoining pressure isotherm microdrops, microdepressions, and equilibrium periodic films exist on flat solid substrates. Criteria are found for both the existence and the stability of these nonflat equilibrium liquid shapes. It is shown that a transition from thick films to thinner films can go via intermediate nonflat states, microdepressions and periodic films, which both can be more stable than flat films within some range of hydrostatic pressure. Experimental investigations of shapes of the predicted nonflat layers can open new possibilities of determination of disjoining pressure in the range of thickness in which flat films are unstable.     

Near-Contact Motion of Surfactant-Covered Spherical Drops: Ionic Surfactant

A lubrication analysis is presented for near-contact axisymmetric motion of spherical drops covered with an insoluble nondiffusing surfactant. The surfactant equation of state is arbitrary; detailed results are presented for ionic surfactants. The qualitative behavior of the system is determined by the dimensionless force parameter &Fcirc;, the external force normalized by the maximum resistance force generated by Marangoni stresses. For &Fcirc; > 1 drops coalesce on a time scale commensurate with the coalescence time tau0 for drops with clean interfaces. For &Fcirc; < 1, the system evolves on the time scale tau0 until Marangoni stresses approximately balance the external force; thereafter a slow evolution occurs on the Stokes time scale. In the long-time regime a self-similar surfactant concentration profile is attained that scales with the extent of the near-contact region. The gap width decreases exponentially with time but slower than for rigid particles because of surfactant backflow. For &Fcirc; < 1, drop coalescence does not occur without van der Waals attraction. Quantitative results depend only moderately on the surfactant equation of state. Copyright 1999 Academic Press.     

Faxen's Laws of a Composite Sphere under Creeping Flow Conditions

Under creeping flow conditions, Faxen's laws are derived for a composite sphere comprising a solid core covered by a permeable layer of arbitrary thickness. The derivations are carried out by applying reciprocal theorem in combination with fluid velocity and pressure distributions in certain simple flow as a comparison field. In this regard, the fluid velocity disturbances caused by a composite sphere subject to a simple shear flow and a rotational flow are solved individually. In the limiting case where the solid core vanishes, the resulting Faxen expressions for the drag force, torque, and stresslet compare very well with the existing Faxen's law for a porous sphere. It is found that when the porous layer is thick enough and its permeability is sufficiently low, the hydrodynamic behavior of a composite sphere can be approximated by that of a porous particle with equal permeability. This can be explained by the fact that the fluid cannot penetrate deeply into a porous layer of low permeability to flow through the pores near the core surface, and thereby the fluid can hardly feel the resistance from the core surface. Copyright 2000 Academic Press.     

Hydrodynamic forces acting on a microscopic emulsion drop growing at a capillary tip in relation to the process of membrane emulsification

Here, we calculate the hydrodynamic ejection force acting on a microscopic emulsion drop, which is continuously growing at a capillary tip. This force could cause drop detachment in the processes of membrane and microchannel emulsification, and affect the size of the released drops. The micrometer-sized drops are not deformed by gravity and their formation happens at small Reynolds numbers despite the fact that the typical period of drop generation is of the order of 0.1 s. Under such conditions, the flow of the disperse phase through the capillary, as it inflates the droplet, engenders a hydrodynamic force, which has a predominantly viscous (rather than inertial) origin. The hydrodynamic boundary problem is solved numerically, by using appropriate curvilinear coordinates. The spatial distributions of the stream function and the velocity components are computed. The hydrodynamic force acting on the drop is expressed in terms of three universal functions of the ratio of the pore and drop radii. These functions are computed numerically. Interpolation formulas are obtained for their easier calculation. It turns out that the increase in the viscosity of each of the two liquid phases increases the total ejection force. The results could find applications for the interpretation and prediction of the effect of hydrodynamic factors on the drop size in membrane emulsification.     

Wetting film stability and flotation kinetics

Single bubble experiments performed with different size fractions of quartz particles and different, but known, contact angles revealed two modes of flotation dynamics in superclean water. (1.) A monotonic increase of collection efficiency Ecoll with increasing particle size was observed at high particle hydrophobicity and, correspondingly, a low wetting film stability (WFS). (2.) At low particle hydrophobicity and, correspondingly, high WFS, an extreme dependence of Ecoll on particle size was observed. The use of superclean water in our experiments prevented the retardation of bubble surface movement caused by surfactants or other impurities that is usual for other investigations and where particle-bubble inertial hydrodynamic interactions are suppressed. In the present study the free movement of the bubble surface enhances particle-bubble inertial interaction, creating conditions for different flotation modes, dependent on WFS. At the instant of inertial impact, a particle deforms the bubble surface, which may cause its rebound. Where the stability of the thin water film, formed between opposing surfaces of a bubble and a particle, is low, its rupture is accompanied with three phase contact line extension and contact angle formation before rebound. This prevents rebound, i.e. the first collision is accompanied by attachment. A high WFS prevents rupture during an impact. As a result, a contact angle does not arise and rebound is not prevented. However, rebound is accompanied by a second collision, the kinetic energy of which is smaller and can cause attachment at repetitive collision. These qualitative considerations are confirmed by the model quantification and comparison with measured Ecoll. For the first time the Sutherland equation (SE) for Ecoll is confirmed by experiment for smaller particle sizes and, correspondingly, very small Stokes numbers. The larger the particle size, the larger is the measured deviation from the SE. The SE is generalized, accounting for the centrifugal force, pressing hydrodynamic force and drainage in the low WFS case and, correspondingly, attachment occurs at first collision or during sliding. The derived generalized Sutherland equation (GSE) describes experimental data at low WFS. However, its application without account for possible rebound does not explain the measured extreme dependence in the case of high WFS. The theory for drainage during particle impact and the beginning of rebound enables conditions for either attachment or rebound in terms of the normal component of the impact velocity and the critical film thickness to be derived. Combining this condition with the GSE allowed the equation for Ecoll to be derived, accounting for attachment area shrinkage and attachment during a repetitive collision. This equation predicts the extreme dependence. Thus the WFS determines the modes of flotation dynamics and, in turn, probably affects the mechanisms, which control the flotation domain. At low WFS its upper boundary is controlled by the stability of the particle-bubble aggregate. At high WFS the upper boundary can be controlled by rebound because the latter reduces the attachment efficiency by a factor of 30 or more even with repetitive collision.