Film rupture in the diffuse interface model coupled to hydrodynamics

Analytical solutions for the profile of two-dimensional droplets with finite-length precursor films

By means of the lubrication approximation we obtain the full family of static bidimensional profiles of a liquid resting on a substrate under partial-wetting conditions imposed by a disjoining-conjoining pressure. We show that for a set of quite general disjoining-conjoining pressure potentials, the free surface can adopt only five nontrivial static patterns; in particular, we find solutions when the height goes to zero which describe satisfactorily the complete free surface for a finite amount of fluid deposited on a substrate. To test the extension of the applicability of our solutions, we compare them with those obtained when the lubrication approximations are not employed and under conditions where the lubrication hypothesis are not strictly valid, and also with axisymmetric solutions. For a given disjoining-conjoining potential, we report a new analytical solution that accounts for all the five possible solutions.

Films, layers, and droplets: The effect of near-wall fluid structure on spreading dynamics

We present a study of the spreading of liquid droplets on a solid substrate at very small scales. We focus on the regime where effective wetting energy (binding potential) and surface tension effects significantly influence steady and spreading droplets. In particular, we focus on strong packing and layering effects in the liquid near the substrate due to underlying density oscillations in the fluid caused by attractive substrate-liquid interactions. We show that such phenomena can be described by a thin-film (or long-wave or lubrication) model including an oscillatory Derjaguin (or disjoining or conjoining) pressure and explore the effects it has on steady droplet shapes and the spreading dynamics of droplets on both an adsorption (or precursor) layer and completely dry substrates. At the molecular scale, commonly used two-term binding potentials with a single preferred minimum controlling the adsorption layer height are inadequate to capture the rich behavior caused by the near-wall layered molecular packing. The adsorption layer is often submonolayer in thickness, i.e., the dynamics along the layer consists of single-particle hopping, leading to a diffusive dynamics, rather than the collective hydrodynamic motion implicit in standard thin-film models. We therefore modify the model in such a way that for thicker films the standard hydrodynamic theory is realized, but for very thin layers a diffusion equation is recovered.

Stability of ultrathin nanocomposite polymer films controlled by the embedding of gold nanoparticles

Thin and ultrathin polymer films combined with nanoparticles (NPs) are of significant interest as they are used in a host of industrial applications. In this paper we describe the stability of such films (hpoly ≤ 30 nm) to dewetting, specifically, how the development of a spinodal instability in a composite NP-polymer layer is controlled by the embedding of Au NPs. At working temperatures (T = 170 °C) above the polymer glass transition temperature (Tg ≈ 100 °C) the absence of Au NPs leads to film rupture by nucleation dewetting, while their presence over a large surface area enhances the development of a spinodal instability without destroying the film continuity. When the NPs embed, the surface undulations are suppressed. The dynamics change from an unstable to a stable state, and the thin composite NP-polymer layer returns to a flat configuration, while the wavelength of the pattern remains constant. Moreover, we demonstrate from a thermodynamic perspective that NPs will remain on the surface or embed in the polymer film depending on their free energy, which is determined by the NP interactions with the underlying polymer, the native SiOx layer, and the Si substrate.

Interaction of 3D dewetting nanodroplets on homogeneous and chemically heterogeneous substrates

Long-time interaction of dewetting nanodroplets is investigated using a long-wave approximation method. Although three-dimensional (3D) droplets evolution dynamics exhibits qualitative behavior analogous to two-dimensional (2D) dynamics, there is an extensive quantitative difference between them. 3D dynamics is substantially faster than 2D dynamics. This can be related to the larger curvature and, as a consequence, the larger Laplace pressure difference between the droplets in 3D systems. The influence of various chemical heterogeneities on the behavior of droplets has also been studied. In the case of gradient surfaces, it is shown how the gradient direction could change the dynamics. For a chemical step located between the droplets, the dynamics is enhanced or weakened depending on the initial configuration of the system.

Patterned deposition at moving contact lines

When a simple or complex liquid recedes from a smooth solid substrate it often leaves a homogeneous or structured deposit behind. In the case of a receding non-volatile pure liquid the deposit might be a liquid film or an arrangement of droplets depending on the receding speed of the meniscus and the wetting properties of the system. For complex liquids with volatile components as, e.g., polymer solutions and particle or surfactant suspensions, the deposit might be a homogeneous or structured layer of solute--with structures ranging from line patterns that can be orthogonal or parallel to the receding contact line via hexagonal or square arrangements of drops to complicated hierarchical structures. We review a number of recent experiments and modelling approaches with a particular focus on mesoscopic hydrodynamic long-wave models. The conclusion highlights open question and speculates about future developments.

Final state of a perturbed liquid film inside a container under the effect of solid-liquid molecular forces and gravity

We investigate theoretically the possible final stationary configurations that can be reached by a laterally confined uniform liquid film inside a container. The liquid is under the action of gravity, surface tension, and the molecular interaction with the solid substrate. We study the case when the container is in an upright position as well as when it is turned upside down. The governing parameters of the problem are the initial thickness of the film, the size of the recipient that contains the liquid, and a dimensionless number that quantifies the relative strength of gravity with respect to the molecular interaction. The uniform film is always a possible final state and depending on the value of the parameters, up to three different additional final states may exist, each one consisting in a droplet surrounded by a thin film. We derive analytical expressions for the energy of these possible final configurations and from these we analyze which state is indeed reached. A uniform thin film may show three different behaviors after a perturbation: The system recovers its initial shape after any perturbation, the fluid evolves towards a drop (if more than one is possible, it tends toward that with the thinnest precursor film) for any perturbation, or the system ends as a uniform film or a drop depending on the details of the perturbation.

Stability of thin liquid films and sessile droplets under confinement

The stability of nonvolatile thin liquid films and of sessile droplets is strongly affected by finite size effects. We analyze their stability within the framework of density functional theory using the sharp kink approximation, i.e., on the basis of an effective interface Hamiltonian. We show that finite size effects suppress spinodal dewetting of films because it is driven by a long-wavelength instability. Therefore nonvolatile films are stable if the substrate area is too small. Similarly, nonvolatile droplets connected to a wetting film become unstable if the substrate area is too large. This instability of a nonvolatile sessile droplet turns out to be equivalent to the instability of a volatile drop which can attain chemical equilibrium with its vapor.

Deformation and dewetting of thin liquid films induced by moving gas jets

We study the deformation of thin liquid films subjected to impinging air-jets that are moving with respect to the substrate. The height profile and shape of the deformed liquid film is evaluated experimentally and numerically for different jet Reynolds numbers and translation speeds, for different liquids and substrate materials. Experiments and numerical results are in good agreement. On partially wetting substrates film rupture occurs. We imaged the appearance of dry spots and emergence of droplet patterns by high-speed, dual-wavelength interference microscopy. We systematically evaluated the resulting average droplet size and droplet density as a function of the experimental conditions. We show that within experimental accuracy the distribution of dry spots is dependent only on the residual film thickness and is not directly influenced by the shear stress and pressure gradients of the air-jet, nor by the speed of the substrate.

Closed-form expression for the profile of partially wetting two-dimensional droplets under gravity

Analytical solutions for the shape of both hanging and sitting droplets under the effects of gravity and surface tension are presented. The modeling also includes the action of molecular forces arising between the liquid and the substrate, which are responsible for the formation of a stable nanometric film in the region close to the droplet contact line. The shape of the droplet is completely described by an analytical solution that also accounts for the pancake-shaped droplets as a limiting case. We find expressions that relate microscopic and nanoscopic aspects, such as the strengths of the molecular forces and the thickness of the nanometric film, to macroscopic quantities, such as the cross-sectional area and the width of the droplet. We study the effect of gravity on the contact angle and find that for small droplets the contact angle follows a power law with the droplet's size. For sitting droplets we find that the there is an upper limit for the value of the gravity.

Rupture of thin liquid films induced by impinging air-jets

Thin liquid films on partially wetting substrates are subjected to laminar axisymmetric air-jets impinging at normal incidence. We measured the time at which film rupture occurs and dewetting commences as a function of diameter and Reynolds number of the air-jet. We developed numerical models for the air flow as well as the height evolution of the thin liquid film. The experimental results were compared with numerical simulations based on the lubrication approximation and a phenomenological expression for the disjoining pressure. We achieved quantitative agreement for the rupture times. We found that the film thickness profiles were highly sensitive to the presence of minute quantities of surface-active contaminants.

Dynamical model for the formation of patterned deposits at receding contact lines

We describe the formation of deposition patterns that are observed in many different experiments where a three-phase contact line of a volatile nanoparticle suspension or polymer solution recedes. A dynamical model based on a long-wave approximation predicts the deposition of irregular and regular line patterns due to self-organized pinning-depinning cycles corresponding to a stick-slip motion of the contact line. We analyze how the line pattern properties depend on the evaporation rate and solute concentration.

Rayleigh and depinning instabilities of forced liquid ridges on heterogeneous substrates

Depinning of two-dimensional liquid ridges and three-dimensional drops on an inclined substrate is studied within the lubrication approximation. The structures are pinned to wetting heterogeneities arising from variations of the strength of the short-range contribution to the disjoining pressure. The case of a periodic array of hydrophobic stripes transverse to the slope is studied in detail using a combination of direct numerical simulation and branch-following techniques. Under appropriate conditions the ridges may either depin and slide downslope as the slope is increased, or first break up into drops via a transverse instability, prior to depinning. The different transition scenarios are examined together with the stability properties of the different possible states of the system.

Thin film evolution equations from (evaporating) dewetting liquid layers to epitaxial growth

In the present contribution we review basic mathematical results for three physical systems involving self-organizing solid or liquid films at solid surfaces. The films may undergo a structuring process by dewetting, evaporation/condensation or epitaxial growth, respectively. We highlight similarities and differences of the three systems based on the observation that in certain limits all of them may be described using models of similar form, i.e. time evolution equations for the film thickness profile. Those equations represent gradient dynamics characterized by mobility functions and an underlying energy functional. Two basic steps of mathematical analysis are used to compare the different systems. First, we discuss the linear stability of homogeneous steady states, i.e. flat films, and second the systematics of non-trivial steady states, i.e. drop/hole states for dewetting films and quantum-dot states in epitaxial growth, respectively. Our aim is to illustrate that the underlying solution structure might be very complex as in the case of epitaxial growth but can be better understood when comparing the much simpler results for the dewetting liquid film. We furthermore show that the numerical continuation techniques employed can shed some light on this structure in a more convenient way than time-stepping methods. Finally we discuss that the usage of the employed general formulation does not only relate seemingly unrelated physical systems mathematically, but does allow as well for discussing model extensions in a more unified way.

Modelling approaches to the dewetting of evaporating thin films of nanoparticle suspensions

We review recent experiments on dewetting thin films of evaporating colloidal nanoparticle suspensions (nanofluids) and discuss several theoretical approaches to describe the ongoing processes including coupled transport and phase changes. These approaches range from microscopic discrete stochastic theories to mesoscopic continuous deterministic descriptions. In particular, we describe (i) a microscopic kinetic Monte Carlo model, (ii) a dynamical density functional theory and (iii) a hydrodynamic thin film model. Models (i) and (ii) are employed to discuss the formation of polygonal networks, spinodal and branched structures resulting from the dewetting of an ultrathin 'postcursor film' that remains behind a mesoscopic dewetting front. We highlight, in particular, the presence of a transverse instability in the evaporative dewetting front, which results in highly branched fingering structures. The subtle interplay of decomposition in the film and contact line motion is discussed. Finally, we discuss a simple thin film model (iii) of the hydrodynamics on the mesoscale. We employ coupled evolution equations for the film thickness profile and mean particle concentration. The model is used to discuss the self-pinning and depinning of a contact line related to the 'coffee-stain' effect. In the course of the review we discuss the advantages and limitations of the different theories, as well as possible future developments and extensions.

Analytical solutions for partially wetting two-dimensional droplets

We present a new analytical solution for the static shape of a two-dimensional droplet in equilibrium with a surrounding thin film on a solid substrate. The modeling includes the effects of capillarity and disjoining-conjoining pressure accounting for intermolecular forces between the solid and the liquid. We derive new analytical solutions for the shape of the droplet, the cross-sectional area, the half-width, and the maximum curvature and inflection points. We study the effects of the size of the droplet on the apparent contact angle. The shape of the droplet in the contact line region is compared with profiles obtained by employing approximations suggested in the literature, and the observed differences are discussed. Finally, we present the time evolution to the steady state to show how the whole profile, including the thin film, evolves to the corresponding stationary configuration.

Dewetting dynamics of stressed viscoelastic thin polymer films

Ultrathin polymer films that are produced, e.g., by spin coating are believed to be stressed since polymers are "frozen in" into out-of-equilibrium configurations during this process. In the framework of a viscoelastic thin-film model, we study the effects of lateral residual stresses on the dewetting dynamics of the film. The temporal evolution of the height profiles and the velocity profiles inside the film as well as the dissipation mechanisms are investigated in detail. Both the shape of the profiles and the importance of frictional dissipation vs viscous dissipation inside the film are found to change in the course of dewetting. The interplay of the nonstationary profiles, the relaxing initial stress, and the changes in the dominance of the two dissipation mechanisms caused by nonlinear friction with the substrate is responsible for the rich behavior of the system. In particular, our analysis sheds a different light on the occurrence of the unexpected maximum in the rim width obtained recently in experiments on polystyrene-polydimethylsiloxane systems.

Front instabilities in evaporatively dewetting nanofluids

Various experimental settings that involve drying solutions or suspensions of nanoparticles-often called nanofluids-have recently been used to produce structured nanoparticle layers. In addition to the formation of polygonal networks and spinodal-like patterns, the occurrence of branched structures has been reported. After reviewing the experimental results we use a modified version of the Monte Carlo model first introduced by Rabani [Nature 426, 271 (2003)] to study structure formation in evaporating films of nanoparticle solutions for the case that all structuring is driven by the interplay of evaporating solvent and diffusing nanoparticles. After introducing the model and its general behavior we focus on receding dewetting fronts which are initially straight but develop a transverse fingering instability. We analyze the dependence of the characteristics of the resulting branching patterns on the driving effective chemical potential, the mobility and concentration of the nanoparticles, and the interaction strength between liquid and nanoparticles. This allows us to understand the underlying instability mechanism.

Emulsions stability, from dilute to dense emulsions -- role of drops deformation

The present paper starts with a review of fundamental descriptions based on physico-chemical laws derived for emulsions with a special interest for eventual evidences of drops deformation. A critical analysis of theories and experiments is given that leads the authors to propose new static and dynamic models for the approach to flocculation and coalescence of two deformable drops in dense and dilute environments of other neighboring drops. The model developed is based on an old paper by Albers and Overbeek for W/O dense emulsions with non-deformable particles, that has been improved recently first by Sengupta and Papadopoulos and then by Mishchuk et al. to account for all the interaction forces (electrostatic, van der Waals and steric). The basic idea here rests in the assumption that the flat surface area of the two coalescing drops, interacting in the field of other particles, increases when the distance between the particles decreases according to an exponential law with a characteristic length related to the disjoining force in the inter-particle film and to the capillary pressure that opposes flattening. The difficulty lies, indeed, in manifold interpretations on experimental observations so that no clear conclusion can be derived on mechanisms responsible for the deformation of droplets. This is why, from a pure theoretical and physical point of view, according to rather complicated models, we propose a much more simple approach that permits to define a capillary length as part of virtual operations. In a static approach, this length is based on analogy with electricity, namely repulsion leads to flatness while attraction to hump. Therefore this brings us to a definition of a length depending on the maximum value of the disjoining pressure in competition with the capillary pressure. Gravity also promotes flocculation, therefore we compare the maximum values of the surface forces acting between the surfaces of two floculating particles to gravity. Finally, considering that in most publications on emulsions foams and colloidal systems, much attention is paid on the role of the drainage in the stability process, we devote the last section to the drainage between flattened drops. We first describe briefly Taylor's approach and extend Reynolds revisited formulae taking into account the viscous friction, the disjoining pressure, the film elasticity and the wetting angle weighting the capillary pressure through the characteristic length. Our calculated values are compared to some experimental data. In conclusion to make this long paper as useful as possible for research purposes, we have the hope that our understanding of emulsion stability is not only based on knowledge of numerous theoretical and experimental works sometimes controversial given in a critical way but that it gives a new approach based on an interpretation of the drop deformation in terms of a characteristic length linked to a deformation number analogous to a Bond number.

Ratchet-driven fluid transport in bounded two-layer films of immiscible liquids

We present a two-layer thin film model allowing us to study the behavior of a general class of 'wettability ratchets' that can be employed to transport a continuous phase. Brownian ratchets, in contrast, are normally used to transport particles or molecules within a continuous carrier fluid without transporting the fluid itself. The wettability ratchet is based on a switchable, spatially asymmetric, periodic interaction of the free surface of the film and the walls. To illustrate the general concept, we focus on an electrical dewetting mechanism based on the effective force exercised by a static electric field on the liquid-liquid interface between two dielectric liquids. In particular, we analyse (i) an on-off ratchet with a constant lateral force resulting in a dewetting-spreading cycle, (ii) a ratchet switching between two shifted potentials that shows a transition between oscillating and sliding drops, and (iii) a flashing external force ratchet. For the three cases, the macroscopic transport is studied in its dependence on spatial and temporal characteristics of the ratchet, and physical properties and volume of the liquids.

Coarsening of unstable thin films subject to gravity

Thin films of viscous fluids coating hydrophobic substrates are unstable to dewetting instabilities, and long-time evolution leads to the formation of an array of near-equilibrium droplets connected by ultrathin fluid layers. In the absence of gravity, previous use of lubrication theory has shown that coarsening dynamics will ensue-the system will evolve by successively eliminating small drops to yield fewer larger drops. While gravity has only a weak influence on the initial thin film, we show that it has a significant influence on the later stages of the coarsening dynamics, dramatically slowing the rate of coarsening for large drops. Small drops are relatively unaffected, but as coarsening progresses, these aggregate into larger drops whose shape and dynamics are dominated by gravity. The change in the mean drop shape causes a corresponding gradual transition from power-law coarsening to a logarithmic behavior.

Formation and mobility of droplets on composite layered substrates

A mesoscale fluid film placed on a solid support may break up and form droplets. In addition, droplets may exhibit spontaneous translation by modifying the wetting properties of the substrate, resulting in asymmetry in the contact angles. We examine mechanisms for droplet formation and motion on uniform and terraced landscapes, i.e., composite substrates. The fluid film stability, droplet formation and velocity are studied theoretically in the isothermal case using a lubrication approach in one spatial dimension. The droplet properties are found to involve contributions from both the terraced layer thickness and molecular interactions via the disjoining potential.

Morphology changes in the evolution of liquid two-layer films

We consider a thin film consisting of two layers of immiscible liquids on a solid horizontal (heated) substrate. Both the free liquid-liquid and the liquid-gas interface of such a bilayer liquid film may be unstable due to effective molecular interactions relevant for ultrathin layers below 100-nm thickness, or due to temperature-gradient-caused Marangoni flows in the heated case. Using a long-wave approximation, we derive coupled evolution equations for the interface profiles for the general nonisothermal situation allowing for slip at the substrate. Linear and nonlinear analyses of the short- and long-time film evolution are performed for isothermal ultrathin layers, taking into account destabilizing long-range and stabilizing short-range molecular interactions. It is shown that the initial instability can be of a varicose, zigzag, or mixed type. However, in the nonlinear stage of the evolution the mode type, and therefore the pattern morphology, can change via switching between two different branches of stationary solutions or via coarsening along a single branch.

Modelling thin-film dewetting on structured substrates and templates: bifurcation analysis and numerical simulations

We study the dewetting process of a thin liquid film on a chemically patterned solid substrate (template) by means of a thin-film evolution equation incorporating a space-dependent disjoining pressure. Dewetting of a thin film on a homogeneous substrate leads to fluid patterns with a typical length scale, that increases monotonously in time (coarsening). Conditions are identified for the amplitude and periodicity of the heterogeneity that allow to transfer the template pattern onto the liquid structure ("pinning") emerging from the dewetting process. A bifurcation and stability analysis of the possible liquid ridge solutions on a periodically striped substrate reveal parameter ranges where pinning or coarsening ultimately prevail. We obtain an extended parameter range of multistability of the pinning and coarsening morphologies. In this regime, the selected pattern depends sensitively on the initial conditions and potential finite perturbations (noise) in the system as we illustrate with numerical integrations in time. Finally, we discuss the instability to transversal modes leading to a decay of the ridges into rows of drops and show that it may diminish the size of the parameter range where the pinning of the thin film to the template is successful.

Coarsening dynamics of dewetting films

Lubrication theory for unstable thin liquid films on solid substrates is used to model the coarsening dynamics in the long-time behavior of dewetting films. The dominant physical effects that drive the fluid dynamics in dewetting films are surface tension and intermolecular interactions with the solid substrate. Instabilities in these films lead to rupture and other morphological changes that promote nonuniformity in the films. Following the initial instabilities, the films break up into near-equilibrium droplets connected by an ultrathin film. For longer times, the fluid will undergo a coarsening process in which droplets both move and exchange mass on slow time scales. The dynamics of this coarsening process will be obtained through the asymptotic reduction of the long-wave PDE governing the thin film to a set of ODEs for the evolution of the droplets. From this, a scaling law that governs the coarsening rate is derived.

Dewetting of thin films on heterogeneous substrates: pinning versus coarsening

We study a model for a thin liquid film dewetting from a periodic heterogeneous substrate (template). The amplitude and periodicity of a striped template heterogeneity necessary to obtain a stable periodic stripe pattern, i.e., pinning, are computed. This requires a stabilization of the longitudinal and transversal modes driving the typical coarsening dynamics during dewetting of a thin film on a homogeneous substrate. If the heterogeneity has a larger spatial period than that of the critical dewetting mode, weak heterogeneities are sufficient for pinning. Our results imply a large region of coexistence between coarsening dynamics and pinning.

Sliding drops in the diffuse interface model coupled to hydrodynamics

Using a film thickness evolution equation derived recently combining long-wave approximation and diffuse interface theory [L. M. Pismen and Y. Pomeau, Phys. Rev. E 62, 2480 (2000)] we study one-dimensional surface profiles for a thin film on an inclined plane. We discuss stationary flat film and periodic solutions including their linear stability. Flat sliding drops are identified as universal profiles, whose main properties do not depend on mean film thickness. The flat drops are analyzed in detail, especially how their velocity, advancing and receding dynamic contact angles and plateau thicknesses depend on the inclination of the plane. A study of nonuniversal drops shows the existence of a dynamical wetting transition with hysteresis between droplike solutions and a flat film with small amplitude nonlinear waves.