Pinning-depinning transition of droplets on inclined substrates with a three-dimensional topographical defect

Droplets on inclined substrates can depin and slide freely above a critical substrate inclination angle. Pinning can be caused by topographical defects on the substrate, and understanding the influence of defect geometry on the pinning-depinning transition is important for diverse applications such as fog harvesting, droplet-based microfluidic devices, self-cleaning surfaces, and inkjet printing. Here, we develop a lubrication-theory-based model to investigate the motion of droplets on inclined substrates with a single three-dimensional Gaussian-shaped defect that can be in the form of a bump or a dent. A precursor-film/disjoining-pressure approach is used to capture contact-line motion, and a nonlinear evolution equation is derived which describes droplet thickness as a function of the position along the substrate and time. The evolution equation is solved numerically using an alternating direction implicit finite-difference scheme to study how the defect geometry influences the critical inclination angle and the shape of a pinned droplet. It is found that the critical substrate inclination angle increases as the defect becomes taller/deeper or wider along the direction lateral to the droplet-sliding direction. However, the critical inclination angle decreases as the defect becomes wider along the sliding direction. Below the critical inclination angle, the advancing contact line of the droplet at the droplet centerline is pinned to the defect at the point having maximum negative slope. Simple scaling relations that reflect the influence of defect geometry on the droplet retention force arising from surface tension are able to account for many of the trends observed in the numerical simulations.

Influence of Surface Roughness on Droplet Evaporation and Absorption: Insights into Experiments from Lubrication-Theory-Based Models

While solid substrates are often idealized as being perfectly smooth, all real surfaces possess some level of topographical and chemical heterogeneity. This heterogeneity can greatly influence droplet dynamics. Mathematical models based on lubrication theory that account for surface roughness reveal how topographical defects induce contact-line pinning and affect the deposition patterns of colloidal particles suspended in the droplet. Contact-line pinning profoundly changes the behavior of droplet evaporation on horizontal and inclined impermeable substrates and droplet absorption on horizontal permeable substrates. Models accounting for surface roughness yield predictions that are qualitatively consistent with experimental observations and also provide insight into the underlying physical mechanisms. These models are a foundation for the exploration of a rich array of problems concerning droplet dynamics which are of both fundamental and practical interest.

Acoustic, Phononic, Brillouin Light Scattering and Faraday Wave-Based Frequency Combs: Physical Foundations and Applications

Frequency combs (FCs)-spectra containing equidistant coherent peaks-have enabled researchers and engineers to measure the frequencies of complex signals with high precision, thereby revolutionising the areas of sensing, metrology and communications and also benefiting the fundamental science. Although mostly optical FCs have found widespread applications thus far, in general FCs can be generated using waves other than light. Here, we review and summarise recent achievements in the emergent field of acoustic frequency combs (AFCs), including phononic FCs and relevant acousto-optical, Brillouin light scattering and Faraday wave-based techniques that have enabled the development of phonon lasers, quantum computers and advanced vibration sensors. In particular, our discussion is centred around potential applications of AFCs in precision measurements in various physical, chemical and biological systems in conditions where using light, and hence optical FCs, faces technical and fundamental limitations, which is, for example, the case in underwater distance measurements and biomedical imaging applications. This review article will also be of interest to readers seeking a discussion of specific theoretical aspects of different classes of AFCs. To that end, we support the mainstream discussion by the results of our original analysis and numerical simulations that can be used to design the spectra of AFCs generated using oscillations of gas bubbles in liquids, vibrations of liquid drops and plasmonic enhancement of Brillouin light scattering in metal nanostructures. We also discuss the application of non-toxic room-temperature liquid-metal alloys in the field of AFC generation.

A mathematical model and mesh-free numerical method for contact-line motion in lubrication theory

ABSTRACT

We introduce a mathematical model with a mesh-free numerical method to describe contact-line motion in lubrication theory. We show how the model resolves the singularity at the contact line, and generates smooth profiles for an evolving, spreading droplet. The model describes well the physics of droplet spreading-including Tanner's Law for the evolution of the contact line. The model can be configured to describe complete wetting or partial wetting, and we explore both cases numerically. In the case of partial wetting, the model also admits analytical solutions for the droplet profile, which we present here.

ARTICLE HIGHLIGHTS

We formulate a mathematical model to regularize the contact-line singularity for droplet spreading.The model can be solved using a fast, accurate mesh-free numerical method.Numerical simulations confirm that the model describes the quantitative aspects of droplet spreading well.

Droplet evaporation on soft solid substrates

Droplet evaporation on soft solid substrates is relevant to applications such as fabrication of microlenses and controlled particle deposition. Here, we develop a lubrication-theory-based model to advance fundamental understanding of the important limiting case of a planar droplet evaporating on a linear viscoelastic solid. A set of partial differential equations describing the time evolution of the liquid-air and liquid-solid interfaces is derived and solved with a finite-difference method. A disjoining-pressure/precursor-film approach is used to describe contact-line motion, and the one sided model is used to describe solvent evaporation. Parametric studies are conducted to investigate the effect of solid properties (thickness, viscosity, shear modulus, wettability) and evaporation rate on droplet dynamics. Our results indicate that softer substrates speed up droplet evaporation due to prolonged pinning of the contact line. Results from our model are able to qualitatively reproduce some key trends observed in experiments. Due to its systematic formulation, our model can readily be extended to more complex situations of interest such as evaporation of particle-laden droplets on soft solid substrates.

The evaporation of surfactant-laden droplets: A comparison between contact line models

HYPOTHESIS

There are two different sharp-interface models for moving contact lines: slip models and precursor film models. While both models predict a mostly constant contact angle during the evaporation of pure droplets, it is expected that they behave differently when surfactants are present, because of the inherent dissimilarities in their respective interface definitions.

SIMULATIONS

Both contact line models are numerically implemented using lubrication theory to analyze evaporating droplets. A convection-diffusion equation is implemented for insoluble surfactants. For pure droplets the models are compared with experiments performed by Nguyen et al. (2012).

FINDINGS

The two contact line models show results comparable to the experiments with pure droplets. If insoluble surfactants are present, the slip model increasingly shows pinning-like behavior as the initial surfactant concentration is increased. This 'quasi-pinning' is found to be consistent with experimental results in literature. The precursor film model, in contrast, shows no significant change when surfactants are added. This lack of change is a result of surfactant flowing from the droplet into the precursor film and vice versa. While suggesting potential solutions to this unphysical behavior, it is concluded that in the context of surfactants, slip models are preferable over precursor film models given the current state of the art.

A thin-film model for droplet spreading on soft solid substrates

The spreading of droplets on soft solid substrates is relevant to applications such as tumor biophysics and controlled droplet condensation and evaporation. In this paper, we apply lubrication theory to advance fundamental understanding of the important limiting case of spreading of a planar droplet on a linear viscoelastic solid. The contact-line region is described by a disjoining-pressure/precursor-film approach, and nonlinear evolution equations describing how the liquid-air and liquid-solid interfaces evolve in space and time are derived and solved numerically. Parametric studies are conducted to investigate the effects of solid thickness, viscosity, shear modulus, and wettability on droplet spreading. Softer substrates are found to speed up spreading for perfectly wetting droplets but slow down spreading for partially wetting droplets. For perfectly wetting droplets, faster spreading is a result of more liquid being pumped toward the contact line due to a larger liquid-film thickness there arising from the repulsive component of the disjoining pressure. In contrast, slower spreading of partially wetting droplets is a result of less liquid being pumped toward the contact line due to a smaller liquid-film thickness there arising from the attractive component of the disjoining pressure. The model predictions for partially wetting droplets are qualitatively consistent with experimental observations, and allow us to disentangle the effects of substrate deformability and wettability on droplet spreading. Due to its systematic formulation, our model can readily be extended to more complex situations involving multiple droplets, substrate inclination, and droplet phase changes.

On Droplet Coalescence in Quasi-Two-Dimensional Fluids

Coalescence of droplets plays a crucial role in nature and modern technology. Various experimental and theoretical studies explored droplet dynamics in three-dimensional (3D) and on 2D solid or liquid substrates. In this paper, we demonstrate the complete coalescence of isotropic droplets in thin quasi-2D liquids-overheated smectic films. We observe the merging of micrometer-sized flat droplets using high-speed imaging and analyze the shape transformations of the droplets on the timescale of milliseconds. Our studies reveal the scaling laws of the coalescence time, which exhibits a different dependence on the droplet geometry from that in the case of droplets on a solid substrate. A theoretical model is proposed to explain the difference in behavior.

Non-axisymmetric elastohydrodynamic solid-liquid-solid dewetting: Experiments and numerical modelling

We have studied the dewetting dynamics of partially wetting liquid films confined between a soft elastic hemisphere and an elastomer layer by means of systematic experiments. We focused on the experimentally most relevant case of non-axisymmetric dewetting, which initiated at the locations of minimum film thickness near the perimeter of the contact area. We found the contact line speed to be highly anisotropic in this case. It is significantly faster in the azimuthal direction along the perimeter of the contact spot than in the radially inwards direction. We developed a three-dimensional, fully coupled numerical model that reproduces many features observed in the experiments.

Harmonic and subharmonic waves on the surface of a vibrated liquid drop

Liquid drops and vibrations are ubiquitous in both everyday life and technology, and their combination can often result in fascinating physical phenomena opening up intriguing opportunities for practical applications in biology, medicine, chemistry, and photonics. Here we study, theoretically and experimentally, the response of pancake-shaped liquid drops supported by a solid plate that vertically vibrates at a single, low acoustic range frequency. When the vibration amplitudes are small, the primary response of the drop is harmonic at the frequency of the vibration. However, as the amplitude increases, the half-frequency subharmonic Faraday waves are excited parametrically on the drop surface. We develop a simple hydrodynamic model of a one-dimensional liquid drop to analytically determine the amplitudes of the harmonic and the first superharmonic components of the linear response of the drop. In the nonlinear regime, our numerical analysis reveals an intriguing cascade of instabilities leading to the onset of subharmonic Faraday waves, their modulation instability, and chaotic regimes with broadband power spectra. We show that the nonlinear response is highly sensitive to the ratio of the drop size and Faraday wavelength. The primary bifurcation of the harmonic waves is shown to be dominated by a period-doubling bifurcation, when the drop height is comparable with the width of the viscous boundary layer. Experimental results conducted using low-viscosity ethanol and high-viscocity canola oil drops vibrated at 70Hz are in qualitative agreement with the predictions of our modeling.

Analysis of the oscillatory wetting-dewetting motion of a volatile drop during the deposition of polymer on a solid substrate

A recent experimental work revealed an oscillatory wetting-dewetting motion of the three phase contact line during the deposition of polymer from a volatile solution. Here we employ a theoretical model to explain the wetting-dewetting motion of the contact line by incorporating opposing evaporation and Marangoni induced flows in the deposition process. We take into account the contribution of polymer concentration to the surface tension of the volatile drop and show that by changing the different parameters of the system we are able to traverse the dynamics of the three phase contact line from a simple dewetting regime to the wetting-dewetting regime, observed in experiment. We further show that deposition patterns, which were previously attributed to stick and stick-slip modes of the contact line motion may be generated by the wetting-dewetting mode. We summarize our theoretical findings in phase diagrams, which show the expected regimes of contact line motion and the resulting types of patterned deposits which are to be obtained under different physical conditions.

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.

Biological applications of kinetics of wetting and spreading

Wetting and spreading kinetics of biological fluids has gained a substantial interest recently. The importance of these fluids in our lives has driven the pace of publications. Globally scientists have ever growing interest in understanding wetting phenomena due to its vast applications in biological fluids. It is impractical to review extremely large number of publications in the field of kinetics of complex biological fluids and cosmetic solutions on diverse surfaces. Therefore, biological and cosmetic applications of wetting and spreading dynamics are considered in the following areas: (i) Spreading of Newtonian liquids in the case of non-porous and porous substrates. It is shown that the spreading kinetics of a Newtonian droplet on non-porous and porous substrate can be defined through theoretical relations for droplet base radius on time, which agree well with the experimental results; (ii) Spreading of blood over porous substrates. It is shown that blood, which has a complex non-Newtonian rheology, can be successfully modelled with the help of simple power-law model for shear-thinning non-Newtonian liquids; (iii) Simultaneous spreading and evaporation kinetics of blood. This part enlightens different underlying mechanisms present in the wetting, spreading, evaporation and dried pattern formation of the blood droplets on solid substrates; (iv) Spreading over hair. In this part the wetting of hair tresses by aqueous solutions of two widely used by industry commercially available polymers, AculynTM 22 and AculynTM 33, are discussed. The influence of non-Newtonian rheology of these polymer solutions on the drainage of foams produced from these solutions is also briefly discussed.

Drying of Droplets of Colloidal Suspensions on Rough Substrates

In many technological applications, excess solvent must be removed from liquid droplets to deposit solutes onto substrates. Often, the substrates on which the droplets rest may possess some roughness, either intended or unintended. Motivated by these observations, we present a lubrication-theory-based model to study the drying of droplets of colloidal suspensions on a substrate containing a topographical defect. The model consists of a system of one-dimensional partial differential equations accounting for the shape of the droplet and depth-averaged concentration of colloidal particles. A precursor film and disjoining pressure are used to describe the contact-line region, and evaporation is included using the well-known one-sided model. Finite-difference solutions reveal that when colloidal particles are absent, the droplet contact line can pin to a defect for a significant portion of the drying time due to a balance between capillary-pressure gradients and disjoining-pressure gradients. The time-evolution of the droplet radius and contact angle exhibits the constant-radius and constant-contact-angle stages that have been observed in prior experiments. When colloidal particles are present and the defect is absent, the model predicts that particles will be deposited near the center of the droplet in a cone-like pattern. However, when a defect is present, pinning of the contact-line accelerates droplet solidification, leading to particle deposition near the droplet edge in a coffee-ring pattern. These predictions are consistent with prior experimental observations, and illustrate the critical role contact-line pinning plays in controlling the dynamics of drying droplets.

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.

Evaporation-triggered microdroplet nucleation and the four life phases of an evaporating Ouzo drop

Evaporating liquid droplets are omnipresent in nature and technology, such as in inkjet printing, coating, deposition of materials, medical diagnostics, agriculture, the food industry, cosmetics, or spills of liquids. Whereas the evaporation of pure liquids, liquids with dispersed particles, or even liquid mixtures has intensively been studied over the past two decades, the evaporation of ternary mixtures of liquids with different volatilities and mutual solubilities has not yet been explored. Here we show that the evaporation of such ternary mixtures can trigger a phase transition and the nucleation of microdroplets of one of the components of the mixture. As a model system, we pick a sessile Ouzo droplet (as known from daily life-a transparent mixture of water, ethanol, and anise oil) and reveal and theoretically explain its four life phases: In phase I, the spherical cap-shaped droplet remains transparent while the more volatile ethanol is evaporating, preferentially at the rim of the drop because of the singularity there. This leads to a local ethanol concentration reduction and correspondingly to oil droplet nucleation there. This is the beginning of phase II, in which oil microdroplets quickly nucleate in the whole drop, leading to its milky color that typifies the so-called "Ouzo effect." Once all ethanol has evaporated, the drop, which now has a characteristic nonspherical cap shape, has become clear again, with a water drop sitting on an oil ring (phase III), finalizing the phase inversion. Finally, in phase IV, all water has evaporated, leaving behind a tiny spherical cap-shaped oil drop.

Spreading of Thin Droplets of Perfect and Leaky Dielectric Liquids on Inclined Surfaces

The spreading of droplets may be influenced by electric fields, a situation that is relevant to applications such as coating, printing, and microfluidics. In this work we study the effects of an electric field on the gravity-driven spreading of two-dimensional droplets down an inclined plane. We consider both perfect and leaky dielectric liquids, as well as perfectly and partially wetting systems. In addition to the effects of electric fields, we examine the use of thermocapillary forces to suppress the growth of the capillary ridge near the droplet front. Lubrication theory is applied to generate a set of coupled partial differential equations for interfacial height and charge, which are then solved numerically with a finite-difference method. Electric fields increase the height of the capillary ridge in both perfect and leaky dielectric droplets due to electrostatic pressure gradients that drive liquid into the ridge. In leaky dielectrics, large interfacial charge gradients in the contact-line region create shear stresses that also enhance ridge growth and the formation of trailing minor ridges. The coalescence of these ridges can significantly affect the long-time thinning rate of leaky dielectric droplets. In partially wetting liquids, electric fields promote the splitting of smaller droplets from the primary droplet near the receding contact line due to the interplay between electrostatic forces and disjoining pressure. Cooling from below and heating from above generates thermocapillary forces that counteract the effects of electric fields and suppress the growth of the capillary ridge. The results of this work have important implications for manipulating the spreading of droplets down inclined surfaces.

A model for pattern deposition from an evaporating solution subject to contact angle hysteresis and finite solubility

We propose a model for the pattern deposition of the solute from an evaporating drop of a dilute solution on a horizontal substrate. In the model we take into account the three-phase contact angle hysteresis and the deposition of the solute whenever its concentration exceeds the solubility limit. The evaporating drop is governed by a film equation. We show that unless for a very small three-phase contact angle or a very rapid evaporation rate the film adopts a quasi-steady geometry, satisfying the Young-Laplace equation to leading order. The concentration profile is assumed to satisfy an advection diffusion equation subject to the standard Fick's law for the diffusive flux. We further use an integral boundary condition to describe the dynamics of the concentration in the vicinity of the three-phase contact line; we replace an exact geometric description of the vicinity of the contact line, which is usually assumed such that mathematical singularities are avoided, with general insights about the concentration and its flux. We use our model to explore the relationships between a variety of deposition patterns and the governing parameters, show that the model repeats previous findings, and suggest further insights.

Efficient modelling of droplet dynamics on complex surfaces

This work investigates the dynamics of droplet interaction with smooth or structured solid surfaces using a novel sharp-interface scheme which allows the efficient modelling of multiple dynamic contact lines. The liquid-gas and liquid-solid interfaces are treated in a unified context and the dynamic contact angle emerges simply due to the combined action of the disjoining and capillary pressure, and viscous stresses without the need of an explicit boundary condition or any requirement for the predefinition of the number and position of the contact lines. The latter, as it is shown, renders the model able to handle interfacial flows with topological changes, e.g. in the case of an impinging droplet on a structured surface. Then it is possible to predict, depending on the impact velocity, whether the droplet will fully or partially impregnate the structures of the solid, or will result in a 'fakir', i.e. suspended, state. In the case of a droplet sliding on an inclined substrate, we also demonstrate the built-in capability of our model to provide a prediction for either static or dynamic contact angle hysteresis. We focus our study on hydrophobic surfaces and examine the effect of the geometrical characteristics of the solid surface. It is shown that the presence of air inclusions trapped in the micro-structure of a hydrophobic substrate (Cassie-Baxter state) result in the decrease of contact angle hysteresis and in the increase of the droplet migration velocity in agreement with experimental observations for super-hydrophobic surfaces. Moreover, we perform 3D simulations which are in line with the 2D ones regarding the droplet mobility and also indicate that the contact angle hysteresis may be significantly affected by the directionality of the structures with respect to the droplet motion.

Stability and break-up of thin liquid films on patterned and structured surfaces

Solid surfaces with chemical patterning or topographical structure have attracted attention due to many potential applications such as manufacture of flexible electronics, microfluidic devices, microscale cooling systems, as well as development of self-cleaning, antifogging, and antimicrobial surfaces. In many configurations involving patterned or structured surfaces, liquid films are in contact with such solid surfaces and the issue of film stability becomes important. Studies of stability in this context have been largely focused on specific applications and often not connected to each other. The purpose of the present review is to provide a unified view of the topic of stability and rupture of liquid films on patterned and structured surfaces, with particular focus on common mathematical methods, such as lubrication approximation for the liquid flow, bifurcation analysis, and Floquet theory, which can be used for a wide variety of problems. The physical mechanisms of the instability discussed include disjoining pressure, thermocapillarity, and classical hydrodynamic instability of gravity-driven flows. Motion of a contact line formed after the film rupture is also discussed, with emphasis on how the receding contact angle is expected to depend on the small-scale properties of the substrate.

Nanodroplet Depinning from Nanoparticles

Nanoscale defects on a substrate affect the sliding motion of water droplets. Using in situ transmission electron microscopy imaging, we visualized the depinning dynamics of water nanodroplets from gold nanoparticles on a flat SiNx surface. Our observations showed that nanoscale pinning effects of the gold nanoparticle oppose the lateral forces, resulting in stretching, even breakup, of the water nanodroplet. Using continuum long wave theory, we modeled the dynamics of a nanodroplet depinning from a nanoparticle of comparable length scales, and the model results are consistent with experimental findings and show formation of a capillary bridge prior to nanodroplet depinning. Our findings have important implications on surface cleaning at the nanoscale.

Simulating droplet motion on virtual leaf surfaces

A curvilinear thin film model is used to simulate the motion of droplets on a virtual leaf surface, with a view to better understand the retention of agricultural sprays on plants. The governing model, adapted from Roy et al. (2002 J. Fluid Mech. 454, 235-261 (doi:10.1017/S0022112001007133)) with the addition of a disjoining pressure term, describes the gravity- and curvature-driven flow of a small droplet on a complex substrate: a cotton leaf reconstructed from digitized scan data. Coalescence is the key mechanism behind spray coating of foliage, and our simulations demonstrate that various experimentally observed coalescence behaviours can be reproduced qualitatively. By varying the contact angle over the domain, we also demonstrate that the presence of a chemical defect can act as an obstacle to the droplet's path, causing break-up. In simulations on the virtual leaf, it is found that the movement of a typical spray size droplet is driven almost exclusively by substrate curvature gradients. It is not until droplet mass is sufficiently increased via coalescence that gravity becomes the dominating force.

Sagging of evaporating droplets of colloidal suspensions on inclined substrates

A droplet of a colloidal suspension placed on an inclined substrate may sag under the action of gravity. Solvent evaporation raises the concentration of the colloidal particles, and the resulting viscosity changes may influence the sag of the droplet. To investigate this phenomenon, we have developed a mathematical model for perfectly wetting droplets based on lubrication theory and the rapid-vertical-diffusion approximation. Precursor films are assumed to be present, the colloidal particles are taken to be hard spheres, and particle and liquid dynamics are coupled through a concentration-dependent viscosity and diffusivity. Evaporation is assumed to be limited by how rapidly solvent molecules can transfer from the liquid to the vapor phase. The resulting one-dimensional system of nonlinear partial differential equations describing the evolution of the droplet height and particle concentration is solved numerically for a range of initial particle concentrations and substrate temperatures. The solutions reveal that the interaction between evaporation and non-Newtonian suspension rheology gives rise to several distinct regimes of droplet shapes and particle concentration distributions. The results provide insight into how evaporation and suspension rheology can be tuned to minimize sagging as well as the well-known coffee-ring effect, an outcome which is important for industrial coating processes.

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.

Instability of liquid Cu films on a SiO2 substrate

We study the instability of nanometric Cu thin films on SiO2 substrates. The metal is melted by means of laser pulses for some tens of nanoseconds, and during the liquid lifetime, the free surface destabilizes, leading to the formation of holes at first and then in later stages of the instability to metal drops on the substrate. By analyzing the Fourier transforms of the SEM (scanning electron microscope) images obtained at different stages of the metal film evolution, we determine the emerging length scales at relevant stages of the instability development. The results are then discussed within the framework of a long-wave model. We find that the results may differ whether early or final stages of the instability are considered. On the basis of the interpretation of the experimental results, we discuss the influence of the parameters describing the interaction of the liquid metal with the solid substrate. By considering both the dependence of dominant length scales on the film thickness and the measured contact angle, we isolate a model which predicts well the trends found in the experimental data.

On the moving contact line singularity: asymptotics of a diffuse-interface model

The behaviour of a solid-liquid-gas system near the three-phase contact line is considered using a diffuse-interface model with no-slip at the solid and where the fluid phase is specified by a continuous density field. Relaxation of the classical approach of a sharp liquid-gas interface and careful examination of the asymptotic behaviour as the contact line is approached is shown to resolve the stress and pressure singularities associated with the moving contact line problem. Various features of the model are scrutinised, alongside extensions to incorporate slip, finite-time relaxation of the chemical potential, or a precursor film at the wall.

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.

Thin-film flows with moving contact lines: An approach to reducing computing time

A numerical method to reduce the computing times of thin-film flows with moving contact lines is presented. The flows of a film and a droplet are calculated in a frame that moves with a nonconstant velocity U(t). The criterion employed to define this velocity is to reduce the maximum height change in the flow's most critical zone. The efficiency of the algorithm in reducing the CPU time is tested in gravity-driven flows, where the computing time is reduced by up to a factor of 13 depending on the parameters of the problem.

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.

Dynamics of a climbing surfactant-laden film II: stability

The linear and nonlinear stability of a spreading film of constant flux and a drop of constant volume, discussed in [1], are examined here. A linear stability analysis (LSA) is carried out to investigate the stability to spanwise perturbations, by linearisation of the two-dimensional (2-D) evolution equations derived in [1] for the film thickness and surfactant concentration fields. The latter correspond to convective-diffusion equations for the surfactant, existing in the form of monomers (present at the free surface and in the bulk) and micelles (present in the bulk). The results of the LSA indicate that the thinning region, present upstream of the leading front in the constant flux case, and the leading ridge in the constant volume case, are unstable to spanwise perturbations. Numerical simulations of the 2-D system of equations demonstrate that the above-mentioned regions exhibit finger formation; the effect of selected system parameters on the fingering patterns is discussed.

Droplet spreading on chemically heterogeneous substrates

Consider the spreading dynamics of a two-dimensional droplet over chemically heterogeneous substrates. Assuming small slopes and strong surface tension effects, a long-wave expansion of the Stokes equations yields a single evolution equation for the droplet thickness. The contact line singularity is removed by assuming slip at the liquid-solid interface. The chemical nature of the substrate is incorporated by local variations in the microscopic contact angle, which appear as boundary conditions in the governing equation. By asymptotically matching the flow in the bulk of the droplet with the flow in the vicinity of the contact lines, we obtain a set of coupled ordinary differential equations for the locations of the two droplet fronts. We verify the validity of our matching procedure by comparing the solutions of the ordinary differential equations with solutions of the full governing equation. The droplet dynamics is examined in detail via a phase-plane analysis. A number of interesting features that are not present in chemically homogeneous substrates are found, such as the existence of multiple equilibria, the pinning of the droplet fronts at localized chemical features, and the possibility for the droplet fronts to exhibit a stick-slip behavior.

Wetting on a spherical wall: influence of liquid-gas interfacial properties

We study the equilibrium of a liquid film on an attractive spherical substrate for an intermolecular interaction model exhibiting both fluid-fluid and fluid-wall long-range forces. We first reexamine the wetting properties of the model in the zero-curvature limit, i.e., for a planar wall, using an effective interfacial Hamiltonian approach in the framework of the well known sharp-kink approximation (SKA). We obtain very good agreement with a mean-field density functional theory (DFT), fully justifying the use of SKA in this limit. We then turn our attention to substrates of finite curvature and appropriately modify the so-called soft-interface approximation (SIA) originally formulated by Napiórkowski and Dietrich [Phys. Rev. B34, 6469 (1986)] for critical wetting on a planar wall. A detailed asymptotic analysis of SIA confirms the SKA functional form for the film growth. However, it turns out that the agreement between SKA and our DFT is only qualitative. We then show that the quantitative discrepancy between the two is due to the overestimation of the liquid-gas surface tension within SKA. On the other hand, by relaxing the assumption of a sharp interface, with, e.g., a simple "smoothing" of the density profile there, markedly improves the predictive capability of the theory, making it quantitative and showing that the liquid-gas surface tension plays a crucial role when describing wetting on a curved substrate. In addition, we show that in contrast to SKA, SIA predicts the expected mean-field critical exponent of the liquid-gas surface tension.

Identification of mobile species in cationic polymer lubricant layer on silicon oxide from AFM and XPS analyses

The nanoscale spreading of a cationic polymer lubricant (CPL) film consisting of polydimethylsiloxane with quaternary ammonium salt side chains on a SiO(2) surface was studied with the disjoining pressure measurements using atomic force microscopy. CPL shows a monotonic decrease in disjoining pressure as the film thickness increases from 1.3 to 4.5 nm, which suggests stable spreading in this thickness range. Comparing the spreading rates calculated from disjoining pressure and the viscosity of CLP to the self-healing time after tribo-contacts revealed that the ionic form may not be the main mobile species. The X-ray photoelectron spectroscopy analysis found that the CPL film on SiO(2) has about 30% of the quaternary ammonium salts (cationic groups) reduced to tertiary amines (neutral groups). The reduced CPL polymer has much lower viscosity than the original CPL polymer and yields a spreading rate consistent with that measured at the macroscale. Thus, the mobile component in the CPL/SiO(2) film responsible for self-healing is concluded to be the reduced tertiary amine components of CPL.

Dynamics of nanodroplets on wettability gradient surfaces

A lubrication model is used to study the dynamics of nanoscale droplets on wettability gradient surfaces. The effects of the gradient size, size of the nanodroplets and the slip on the dynamics have been studied. Our results indicate that the position of the center of mass of the droplets can be well described in terms of a third-order polynomial function of the time of the motion for all the cases considered. By increasing the size of the droplets the dynamics increases. It is also shown that the slip can considerably enhance the dynamics. The results have been compared with the results obtained using theoretical models and molecular dynamics simulations.

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.

Two-dimensional droplet spreading over random topographical substrates

We examine theoretically the effects of random topographical substrates on the motion of two-dimensional droplets via statistical approaches, by representing substrate families as stationary random functions. The droplet shift variance at both early times and in the long-time limit is deduced and the droplet footprint is found to be a normal random variable at all times. It is shown that substrate roughness inhibits wetting, illustrating also the tendency of the droplet to slide without spreading as equilibrium is approached. Our theoretical predictions are verified by numerical experiments.

Modeling the coalescence of sessile droplets

This paper proposes a simple scenario to describe the coalescence of sessile droplets. This scenario predicts a power-law growth of the bridge between the droplets. The exponent of this power law depends on the driving mechanism for the spreading of each droplet. To validate this simple idea, the coalescence is simulated numerically and a basic experiment is performed. The fluid dynamics problem is formulated in the lubrication approximation framework and the governing equations are solved in the commercial finite element software COMSOL. Although a direct comparison of the numerical results with experiment is difficult because of the sensitivity of the coalescence to the initial and operating conditions, the key features of the event are qualitatively captured by the simulation and the characteristic time scale of the dynamics recovered. The experiment consists of inducing coalescence by pumping a droplet through a substrate which grows and ultimately coalesces with another droplet resting on the substrate. The coalescence was recorded using high-speed imaging and also confirmed the power-law growth of the neck.