Stokes flow near the contact line of an evaporating drop

The evaporation of sessile drops in quiescent air is usually governed by vapour diffusion. For contact angles below $90^\circ$, the evaporative flux from the droplet tends to diverge in the vicinity of the contact line. Therefore, the description of the flow inside an evaporating drop has remained a challenge. Here, we focus on the asymptotic behaviour near the pinned contact line, by analytically solving the Stokes equations in a wedge geometry of arbitrary contact angle. The flow field is described by similarity solutions, with exponents that match the singular boundary condition due to evaporation. We demonstrate that there are three contributions to the flow in a wedge: the evaporative flux, the downward motion of the liquid-air interface and the eigenmode solution which fulfils the homogeneous boundary conditions. Below a critical contact angle of $133.4^\circ$, the evaporative flux solution will dominate, while above this angle the eigenmode solution dominates. We demonstrate that for small contact angles, the velocity field is very accurately described by the lubrication approximation. For larger contact angles, the flow separates into regions where the flow is reversing towards the drop centre.Comment: Journal of Fluid Mechanics 709 (2012

On the effect of the atmosphere on the evaporation of sessile droplets of water

An experimental and theoretical study into the effect of the atmosphere on the evaporation of pinned sessile droplets of water is described. The experimental work investigated the evaporation rates of sessile droplets in atmospheres of three different ambient gases (namely, helium, nitrogen and carbon dioxide) at reduced pressure (from 40 to 1000 mbar) using four different substrates(namely, aluminium, titanium, Macor and PTFE) with a wide range of thermal conductivities.Reducing the atmospheric pressure increases the diffusion coefficient of water vapour in the atmosphere and hence increases the evaporation rate. Changing the ambient gas also alters the diffusion coefficient and hence also affects the evaporation rate. A mathematical model that takes into account the effect of the atmospheric pressure and the nature of the ambient gas on the diffusion of water vapour in the atmosphere and the thermal conductivity of the substrate is developed, and its predictions are found to be in encouraging agreement with the experimental results

A mathematical model of the evaporation of a thin sessile liquid droplet : comparison between experiment and theory

A mathematical model for the quasi-steady diffusion-limited evaporation of a thin axisymmetric sessile droplet of liquid with a pinned contact line is formulated and solved. The model generalises the theoretical model proposed by Deegan et al. [Contact line deposits in an evaporating drop, Phys. Rev. E, 62 (2000) 756-765] to include the effect of evaporative cooling on the saturation concentration of vapour at the free surface of the droplet, and the dependence of the coefficient of diffusion of vapour in the atmosphere on the atmospheric pressure. The predictions of the model are in good qualitative, and in some cases also quantitative, agreement with recent experimental results. In particular, they capture the experimentally observed dependence of the total evaporation rate on the thermal conductivities of the liquid and the substrate, and on the atmospheric pressure

Evaporation of a thin droplet on a thin substrate with a high thermal resistance

A mathematical model for the quasi-steady evaporation of a thin liquid droplet on a thin substrate that incorporates the dependence of the saturation concentration of vapour at the free surface of the droplet on temperature is used to examine an atypical situation in which the substrate has a high thermal resistance relative to the droplet (i.e. it is highly insulating and/or is thick relative to the droplet). In this situation diffusion of heat through the substrate is the rate-limiting evaporative process and at leading order the local mass flux is spatially uniform, the total evaporation rate is proportional to the surface area of the droplet, and the droplet is uniformly cooled. In particular, the qualitative differences between the predictions of the present model in this situation and those of the widely used 'basic' model in which the saturation concentration is independent of temperature are highlighted

The strong influence of substrate conductivity on droplet evaporation

We report the results of physical experiments that demonstrate the strong influence of the thermal conductivity of the substrate on the evaporation of a pinned droplet. We show that this behaviour can be captured by a mathematical model including the variation of the saturation concentration with temperature, and hence coupling the problems for the vapour concentration in the atmosphere and the temperature in the liquid and the substrate. Furthermore, we show that including two ad hoc improvements to the model, namely a Newton's law of cooling on the unwetted surface of the substrate and the buoyancy of water vapour in the atmosphere, give excellent quantitative agreement for all of the combinations of liquid and substrate considered

Analytical solution of Stokes flow inside an evaporating sessile drop: Spherical and cylindrical cap shapes

Exact analytical solutions are derived for the Stokes flows within evaporating sessile drops of spherical and cylindrical cap shapes. The results are valid for arbitrary contact angle. Solutions are obtained for arbitrary evaporative flux distributions along the free surface as long as the flux is bounded at the contact line. The field equations, E^4(Psi)=0 and Del^4(Phi)=0, are solved for the spherical and cylindrical cap cases, respectively. Specific results and computations are presented for evaporation corresponding to uniform flux and to purely diffusive gas phase transport into an infinite ambient. Wetting and non-wetting contact angles are considered with the flow patterns in each case being illustrated. For the spherical cap with evaporation controlled by vapor phase diffusion, when the contact angle lies in the range 0<theta_c<pi, the mass flux of vapor becomes singular at the contact line. This condition required modification when solving for the liquid phase transport. Droplets in all of the above categories are considered for the following two cases: the contact lines are either pinned or free to move during evaporation. The present viscous flow behavior is compared to the inviscid flow behavior previously reported. It is seen that the streamlines for viscous flow lie farther from the substrate than the corresponding inviscid ones.Comment: Revised version; in review in Physics of Fluid

Imaging internal flows in a drying sessile polymer dispersion drop using Spectral Radar Optical Coherence Tomography (SR-OCT)

In this work, we present the visualization of the internal flows in a drying sessile polymer dispersion drop on hydrophilic and hydrophobic surfaces with Spectral Radar Optical Coherence Tomography (SR-OCT).We have found that surface features such as the initial contact angle and pinning of the contact line, play a crucial role on the flow direction and final shape of the dried drop. Moreover, imaging through selection of vertical slices using optical coherence tomography offers a feasible alternative compared to imaging through selection of narrow horizontal slices using confocal microscopy for turbid, barely transparent fluids

Asymptotic analysis of the evaporation dynamics of partially wetting droplets

We consider the dynamics of an axisymmetric, partially wetting droplet of a one-component volatile liquid. The droplet is supported on a smooth superheated substrate and evaporates into a pure vapour atmosphere. In this process, we take the liquid properties to be constant and assume that the vapour phase has poor thermal conductivity and small dynamic viscosity so that we may decouple its dynamics from the dynamics of the liquid phase. This leads to a so-called ‘one-sided’ lubrication-type model for the evolution of the droplet thickness, which accounts for the effects of evaporation, capillarity, gravity, slip and kinetic resistance to evaporation. By asymptotically matching the flow near the contact line region and the bulk of the droplet in the limit of a small slip length and commensurably small evaporation and kinetic resistance effects, we obtain coupled evolution equations for the droplet radius and volume. The predictions of our asymptotic analysis, which also include an estimate of the evaporation time, are found to be in excellent agreement with numerical simulations of the governing lubrication model for a broad range of parameter regimes

The shielding effect extends the lifetimes of two-dimensional sessile droplets

We consider the diffusion-limited evaporation of thin two-dimensional sessile droplets either singly or in a pair. A conformal-mapping technique is used to calculate the vapour concentrations in the surrounding atmosphere, and thus to obtain closed-form solutions for the evolution and the lifetimes of the droplets in various modes of evaporation. These solutions demonstrate that, in contrast to in three dimensions, in large domains the lifetimes of the droplets depend logarithmically on the size of the domain, and more weakly on the mode of evaporation and the separation between the droplets. In particular, they allow us to quantify the shielding effect that the droplets have on each other, and how it extends the lifetimes of the droplets