Abstract
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.
Similar content being viewed by others
1 Introduction
The evaporation of sessile droplets has been studied extensively in recent years [1]. Topics of particular interest include particle transport during drying [2,3,4] and the resulting instabilities [5], nano-scale phenomena [6], and droplet lifetimes [7,8,9].
A physically realistic model of droplet evaporation must describe vapour transport in the surrounding atmosphere. In the simplest setting, this requires us to solve Laplace’s equation for the vapour concentration \(\hat{c}\) in the atmosphere, subject to a saturation condition \(\hat{c}=\hat{c}_{\text {sat}}\) on the surface of the droplet and a no-flux condition on the substrate [10, 11]. The diffusive mass flux \(\hat{J}\) from the free surface of the droplet then controls the evolution, and hence the lifetime, of the droplet.
In the limit in which the droplet is thin [9, 12, 13], the problem simplifies further because the profile of the droplet may be neglected when imposing the boundary conditions on \(\hat{c}\). Thus, the mathematical problem typically becomes that of solving for \(\hat{c}\) in a half-space or other large domain, subject to appropriate mixed boundary conditions. Similar mixed boundary-value problems occur in physical contexts including elastostatics [14], electrostatics [15], thermostatics [16], and hydrodynamics [17].
A range of mathematical techniques can be deployed to solve such problems [14, 18]; contributions go back at least as far as the work of Weber [19], who presented what is effectively the vapour concentration field induced by a thin circular droplet. Subsequent work has employed methods including separation of variables [20], orthogonal polynomial expansions [21], Fourier or Hankel transforms [22, 23], and Green’s functions [24].
In two dimensions, additional techniques become available, notably conformal mapping [16, 25]. This makes two-dimensional analogues of droplet evaporation problems appealing from the modeller’s point of view: although two-dimensional problems may be somewhat artificial, their greater tractability allows more thorough analysis to be carried out. However, in two dimensions there is a fundamental difficulty concerning the specification of appropriate boundary conditions [14], which we will overcome, in the spirit of the work of Yarin et al. [26], by considering a suitably relaxed boundary condition.
In practice, droplets rarely occur in isolation, and so it is important to understand how droplets evaporate in the presence of other evaporating droplets. Previous studies of the evaporation of multiple sessile droplets have employed a variety of experimental, numerical and analytical approaches [27,28,29,30,31,32,33,34,35,36,37,38]. The critical difference between the evaporation of single and of multiple droplets is the occurrence of the shielding effect, namely that the presence of other evaporating droplets increases the local vapour concentration, and so each droplet evaporates more slowly than it would in isolation.
Again, analogous problems have been studied in other physical contexts. The first relevant study by Greenwood [39] examined the interaction of large numbers of microcontacts in electric contact theory, treating them as independent at leading order, and introducing an interaction term at higher order. Similar approaches have since been applied to elastic punches [40] and flow through pores [41], and have been put on a more rigorous asymptotic basis [42,43,44]. All these studies essentially considered the equivalent of thin circular droplets in three dimensions; recent work has used a variety of approaches to investigate the closely related problem of the dissolution of immersed nanobubbles and nanodroplets [45,46,47,48].
In this study we consider the evaporation of thin two-dimensional sessile droplets. In Sect. 2 we consider the one-droplet problem. We present the governing equations (Sect. 2.1), show that the most apparently natural problem does not have a solution (Sect. 2.2), and then show that by considering a suitably relaxed boundary condition we can obtain a physically acceptable solution via a conformal-mapping technique (Sect. 2.3). We validate this solution against numerical simulations (Sect. 2.4), and use it to obtain closed-form solutions for the evolution and lifetimes of the droplet in various modes of evaporation (Sect. 2.5). We then develop asymptotic expressions for these lifetimes in a large domain (Sect. 2.5.4). In Sect. 3 we consider the two-droplet problem. We obtain a solution to this problem (Sect. 3.1), which we again validate against numerical simulations (Sect. 3.2), before using it to obtain closed-form solutions for the evolution and lifetimes of the droplets (Sect. 3.3). We develop asymptotic expressions for these lifetimes (Sect. 3.3.3), and use these expressions to compare the lifetimes of a single droplet and a pair of droplets in dimensional terms (Sect. 4).
2 One-droplet problem
2.1 Model
Consider a thin two-dimensional sessile droplet with constant surface tension \(\hat{\sigma }\) and density \(\hat{\rho }\), evaporating in the diffusion-limited regime. (For simplicity, we shall refer to the fluid throughout as a droplet; viewed in three dimensions it is more accurately described as a ridge or line.) Let it have semi-width \(\hat{R}(\hat{t})\), contact angle \(\hat{\theta }(\hat{t})\) and cross-sectional area \(\hat{A}(\hat{t})\). Using Cartesian co-ordinates \((\hat{x},\hat{y})\) with origin at the centre of the base of the droplet, the droplet evaporates into a surrounding atmosphere with constant coefficient of vapour diffusion \(\hat{D}\), vapour saturation concentration \(\hat{c}=\hat{c}_\text {sat}\), and ambient vapour concentration \(\hat{c}=\hat{c}_\infty \, (<\hat{c}_\text {sat})\). The vapour concentration in the atmosphere is denoted by \(\hat{c}(\hat{x},\hat{y},\hat{t})\), and the diffusive mass flux from the surface of the droplet by \(\hat{J}(\hat{x},\hat{t})\).
Following the approach of [9, 12, 13], we nondimensionalise and scale according to
where \(\hat{R}_0=\hat{R}(0)\) and \(\hat{\theta }_0=\hat{\theta }(0)\).
The vapour concentration is assumed to be quasi-steady, and so c satisfies Laplace’s equation
throughout the atmosphere.
Assuming that the droplet is sufficiently small, the Eötvös–Bond number \(Eo = \hat{\rho }\hat{g}\hat{R}^2_0/\hat{\sigma }\) will be small; under these conditions the free surface of the droplet is approximately parabolic and its cross-sectional area is given by
The flux from the droplet is given by
which may be evaluated at \(y=0\) due to the thinness of the droplet. Similarly, the saturation condition, \(c=1\), on the surface of the droplet may also be imposed on \(y=0\).
The saturation condition on the droplet and the no-flux condition on the substrate thus become
respectively. To complete the problem we require a suitable boundary condition to be imposed in the “far field”; this turns out to be non-trivial to specify.
2.2 Absence of a solution in an infinite half-space
The simplest problem to specify is evaporation into an infinite half-space, so we aim to solve (2) subject to the far-field condition
as well as to a mixed boundary condition on \(y=0\) of the form
Applying a cosine transform to (2) and imposing the far-field condition (6) leads to a solution of the form
where the function A(u) is to be determined. Imposing the boundary condition (7) requires that
The work of Sneddon [14, Sect. 4.5] shows that requiring regularity of c at the contact line \(x=R\) imposes the condition
so specifying that the function f(x) is any positive constant is not an admissible boundary condition, and so, as could have been anticipated from the behaviour of the fundamental solution of Laplace’s equation in two dimensions, the problem specified by (2), (5) and (6) has no solution. We note that (11) precludes not only solutions to the simplest problem in which the saturation concentration is constant on the droplet, but also solutions to more general problems in which it varies along the droplet surface due, for example, to changes in temperature [11, 49].
2.3 Solution in a finite domain via conformal mapping
Since the most apparently natural problem does not have a solution, we instead look for a closely related analogue that does. We therefore consider a slightly modified problem in which the far-field condition (6) is replaced by a similar Dirichlet condition at a distant, but finite, boundary. We therefore aim to solve
subject to the standard boundary conditions on \(y=0\),
and the relaxed boundary condition
While it is difficult to find an analytical solution in a domain that is exactly semi-circular, we can obtain a solution in a semi-elliptical domain that approaches a semi-circular shape as it becomes large.
We proceed using conformal mapping. Let
Then the mapping
maps the semi-infinite strip \((u,v)\in (-1,1)\times (0,\infty )\) in the w-plane to the upper half of the z-plane. In particular, the rectangle \((-1,1)\times (0,S)\) shown in Fig. 1a is mapped to the semi-ellipse with semi-major axis length \(\sqrt{\varPsi ^2+R^2}\) and semi-minor axis length \(\varPsi \) shown in Fig. 1b and given by
where
An important point to note is that the shape of the semi-elliptical domain in the z-plane given by (17) depends on R as well as on \(\varPsi \). Thus, in general, for a droplet whose semi-width changes as it evolves, the shape of the domain also changes. However, in the regime of most interest, \(\varPsi \gg R\), in which the domain is large, Eq. (17) gives
and so the domain is semi-circular with radius \(\varPsi \) and independent of R up to \(O(\varPsi ^{-2}) \ll 1\).
In the rectangular domain in the w-plane we seek a harmonic function \(\varPhi (u,v)\) satisfying
Solving the problem for \(\varPhi \) in the rectangular domain immediately gives the corresponding solution for c in the semi-elliptical domain.
By inspection, the solution in the rectangular domain is given by
so
and the flux is given by
In particular, the flux satisfies
and so it has the same (integrable) square-root singularity at the contact line \(x=R\) as in the corresponding three-dimensional problem [10].
2.4 Numerical validation
In order to validate the solution obtained in Sect. 2.3 (i.e. in order to assess the accuracy of the solution and to quantify the effect of the non-circularity of the domain), we solved the problem in the semi-circular domain using COMSOL Multiphysics [50]. In Fig. 2 we compare these numerical solutions to the analytical solutions in the semi-elliptical domain given by (22) and (23).
Figure 2a–c shows solutions for the vapour concentration along the x-axis, c(x, 0); d–f shows solutions for the vapour concentration along the y-axis, c(0, y); g–i shows solutions for the flux, J(x). The first column [(a), (d), (g)] shows results for \(\varPsi =2\), the second column [(b), (e), (h)] shows results for \(\varPsi =10\), and the third column [(c), (f), (i)] shows results for \(\varPsi =100\). In all cases the solid lines denote the analytical solutions in the semi-elliptical domain, and the dashed lines denote the corresponding numerical solutions in the semi-circular domain. Figure 2a–f shows that the vapour concentration c takes its saturation value \(c=1\) on the surface of the droplet and decreases monotonically to its ambient value \(c=0\) at the edge of the domain, and Fig. 2g–i shows that the flux J increases monotonically with distance from a minimum value at the centre of the droplet and is singular at the contact line \(x=R\). Figure 2 also shows that the analytical solutions accurately capture the behaviour of the numerical solutions provided that \(\varPsi \) is sufficiently large, which is exactly as expected since it is for smaller domains that the semi-circular and semi-elliptical domains are most different.
2.5 Evolution and lifetime of the droplet
The rate of change of the cross-sectional area (3) is given by the flux (23) integrated over the surface of the droplet,
In order to integrate (25) to determine the evolution and lifetime of the droplet, we require additional information about the behaviour of R(t) and \(\theta (t)\), i.e. we need to specify the mode in which the droplet is evaporating. The works by Stauber et al. [7, 8] and Schofield et al. [9] describe various modes of evaporation for axisymmetric droplets. In the present work, we consider the two-dimensional analogues of three of these modes: the constant-radius (CR) mode, the constant-angle (CA) mode, and the stick–slide (SS) mode. (Throughout, for consistency with the three-dimensional problem, we will refer to modes in which R is fixed as “constant-radius” modes, although strictly R is not the radius but the semi-width of the two-dimensional droplet.)
2.5.1 Constant-radius (CR) mode
In the constant-radius (CR) mode, \(R(t)\equiv R_0=1\). Noting that \(\theta (0)=\theta _0=1\), we may immediately integrate (25) to obtain
Hence the lifetime of a single droplet evaporating in the CR mode is
Figure 3a, c, d shows the evolution and lifetime of a single droplet evaporating in the CR mode. The contact angle \(\theta \) and the area \(A\) both decrease linearly with time t (Fig. 3a, c). As the size of the domain \(\varPsi \) increases, the contact angle \(\theta \) and the area \(A\) decrease more slowly, and so the lifetime \(t_\text {CR}\) increases monotonically with \(\varPsi \) (Fig. 3d). This is because in two dimensions the distance to the outer boundary sets the distance over which the concentration decays to zero, and thus controls the concentration gradient close to the droplet, as seen in (23). This strong role of the outer boundary is a fundamental difference from the corresponding three-dimensional problem, in which the distance to the outer boundary becomes irrelevant for a sufficiently large domain, and so a far-field boundary condition can be safely imposed “at infinity”.
2.5.2 Constant-angle (CA) mode
In the constant-angle (CA) mode, \(\theta (t)\equiv \theta _0=1\). Noting that \(R(0)=R_0=1\), we may integrate (25) implicitly to obtain
Hence the lifetime of a single droplet evaporating in the CA mode is
which can be re-written as
Figure 3b–d shows the evolution and lifetime of a single droplet evaporating in the CA mode. In contrast to the CR mode, the semi-width R and the area \(A\) now both decrease nonlinearly with time t (Fig. 3b, c). However, as in the CR mode, the lifetime \(t_\text {CA}\) increases monotonically with \(\varPsi \) (Fig. 3d). Figure 3d also illustrates, as (30) also shows, that due to its pinned contact lines, a droplet evaporating in the CR mode always has a larger surface area, and hence a larger total flux and thus a shorter lifetime, than the same droplet evaporating in the CA mode, i.e. \(t_\text {CR}\le t_\text {CA}\) for all \(\varPsi \).
2.5.3 Stick–slide (SS) mode
In the stick–slide (SS) mode, the contact line is initially pinned, while the contact angle decreases until it reaches its critical de-pinning (receding) value \(\theta =\theta ^*\) (\(0\le \theta ^*\le 1\)) at the de-pinning time \(t=t^*\). After de-pinning, the contact angle remains at its critical value, while the semi-width decreases until it reaches zero. Thus we have
In the pinned (i.e. the CR) phase, \(0<t<t^*\), the droplet evolves according to (26), so that
while in the de-pinned (i.e. the CA) phase, \(t^*<t<t_\text {SS}\), the droplet evolves according to
Combining (32) and (33), the lifetime of a single droplet evaporating in the SS mode is
which can be re-written as
Figure 4 shows the evolution and lifetime of a single droplet evaporating in the SS mode. The de-pinning time \(t^*\) decreases linearly with \(\theta ^*\), while the lifetime \(t_\text {SS}\) increases linearly with \(\theta ^*\) (Fig. 4d). Comparing (27), (30) and (35) shows that, as might have been anticipated, the lifetime of a droplet evaporating in the SS mode always lies between those of the same droplet in the CR and CA modes, i.e. \(t_\text {CR}\le t_\text {SS}\le t_\text {CA}\) for all \(\varPsi \) and \(\theta ^*\). In the limit \(\theta ^*\rightarrow 1^-\) the SS mode approaches the CA mode and thus \(t_\text {SS}\rightarrow {t_\text {CA}}^-\), while in the limit \(\theta ^*\rightarrow 0^+\) the SS mode approaches the CR mode and thus \(t_\text {SS}\rightarrow t_\text {CR}{^+}\).
We note that in Fig. 4a all of the curves for which \(\theta ^*\ne 0\) intersect at \(t=t_\text {CR}\), and from (27), (32) and (33), the semi-width of the droplet at this time, \(R(t_\text {CR})\), satisfies
Note that \(R(t_\text {CR})\) is a monotonically decreasing function of \(\varPsi \) which takes its maximum value \(R(t_\text {CR}) = 1/2\) in the limit \(\varPsi \rightarrow 0^+\) and satisfies \(R(t_\text {CR}) \rightarrow 0^+\) as \(\varPsi \rightarrow \infty \).
2.5.4 Asymptotic behaviour of the lifetimes in a large domain, \(\varPsi \gg R\)
Consider the regime of most interest, \(\varPsi \gg R\), in which the domain is large and approximately semi-circular, and the condition at the outer boundary corresponds most closely to a far-field condition. From (27), (29) and (34) we obtain
respectively. Equations (37)–(39) show that the lifetimes of the droplets all depend logarithmically on the size of the domain, and differ by an amount of \(O(1)\) which depends on the mode of evaporation. The corrections at \(O(\varPsi ^{-2})\) are of the same order as the deviation of the domain from circularity.
3 Two-droplet problem
We now consider the analogous two-droplet problem in the \(\zeta \)-plane, where \(\zeta =\eta +\mathrm{i}\xi \). We assume that the droplets are identical, and use the initial semi-width of the droplets as the characteristic length scale in the non-dimensionalisation. The droplets are located so that they have inner contact lines at \(\eta =\pm \,I\) and outer contact lines at \(\eta =\pm \varOmega \), where \(\varOmega >I\), as shown in Fig. 5. The cross-sectional area of each droplet is then given by
3.1 Solution in a finite domain via conformal mapping
Consider the conformal map
from the z-plane to the \(\zeta \)-plane. This maps the real interval (0, R) in the z-plane to the real interval \((I,\varOmega )\) where \(\varOmega =\sqrt{I^2+R^2}\) in the \(\zeta \)-plane, preserving the saturation condition on the droplet. It maps the real interval \((R,\sqrt{\varPsi ^2+R^2})\) in the z-plane to the real interval \((\varOmega ,\sqrt{\varPsi ^2+\varOmega ^2})\) in the \(\zeta \)-plane, preserving the no-flux condition on the substrate. It maps the imaginary interval \((0,\mathrm{i}I)\) in the z-plane to the real interval \((0,I)\) in the \(\zeta \)-plane: the symmetry condition on the imaginary axis in the z-plane now becomes a no-flux condition on the real interval \((0,I)\) in the \(\zeta \)-plane. With a suitable choice of branch cut, the equivalent regions in the left half of the z-plane are mapped into corresponding regions in the left half of the \(\zeta \)-plane. The outer boundary of the domain, given by the rectangle in the z-plane and the semi-ellipse (17) in the z-plane, is mapped to the quasi-semi-elliptical curve shown in Fig. 5 and given by
However, as in the one-droplet problem, in the regime of most interest, \(\varPsi \gg I\), in which the domain is large, Eq. (42) gives
and so the domain is again semi-circular with radius \(\varPsi \) and independent of \(I\) and \(\varOmega \) up to \(O(\varPsi ^{-2}) \ll 1\).
The solution in the quasi-semi-elliptical domain is given by
Figure 6 shows the contours of the vapour concentration \(c(\eta ,\xi )\) for the two-droplet problem given by (44) and the corresponding contours of c(x, y) for the one-droplet problem given by (22). In both cases, far from the droplet(s) the contours approach the (semi-elliptical or quasi-semi-elliptical) shape of the outer boundary, and near to the droplet(s) the contours approach the flat shape(s) of the droplet(s). For the two-droplet problem the concentration in the region between the droplets is increased relative to that in the one-droplet problem, and near to the droplets the concentration falls away more gradually than it does in the one-droplet problem, resulting in the shielding effect described in Sect. 1.
The flux is given by
In particular, the flux satisfies
confirming that it again has square-root singularities at both contact lines.
3.2 Numerical validation
As we did in the one-droplet problem, in order to validate the solution obtained in Sect. 3.1, we solved the two-droplet problem in the semi-circular domain using COMSOL Multiphysics [50]. In Fig. 7 we compare these numerical solutions to the analytical solutions in the quasi-semi-elliptical domain given by (44) and (45).
Figure 7a–c shows solutions for the vapour concentration along the \(\eta \)-axis, \(c(\eta ,0)\); d–f shows solutions for the vapour concentration along the \(\xi \)-axis, \(c(0,\xi )\); g–i shows solutions for the flux, \(J(\eta )\). The first column [(a), (d), (g)] shows results for \(\varPsi =4\), the second column [(b), (e), (h)] shows the corresponding results for \(\varPsi =10\), and the final column [(c), (f), (i)] shows the corresponding results for \(\varPsi =100\). In all cases the solid lines denote the analytical solutions in the quasi-semi-elliptical domain, and the dashed lines denote the corresponding numerical solutions in the semi-circular domain. As in the one-droplet problem, Fig. 7 shows that c takes its saturation value on the surface of the droplets and decreases monotonically to its ambient value at the edge of the domain, that J is singular at the contact lines \(x=I\) and \(x=\varOmega \), and that the analytical solutions accurately capture the behaviour of the numerical solutions provided that \(\varPsi \) is sufficiently large.
However, Fig. 7a–f also shows that c decreases monotonically to an (unsaturated) minimum value between the droplets, and that this value is an increasing function of \(\varPsi \): this latter behaviour reflects the smaller concentration gradients, and thus the higher concentrations, which occur near to the droplets in larger domains. In addition, Fig. 7g–i clearly illustrate the shielding effect that the droplets have on each other. In particular, as (46) shows, the flux near to the outer contact line is suppressed less by the presence of the other droplet, and so remains larger than the flux near to the inner contact line, resulting in the skewed flux profiles shown in Fig. 7g–i. In particular, the minimum value of the flux no longer occurs at the centre of each droplet (as it does in the one-droplet problem).
3.3 Evolution and lifetime of the droplets
Using the solution for the flux given by (45), the evolution and lifetime of the droplets are determined by
In the one-droplet problem we investigated three modes of evaporation (namely the CR, CA and SS modes), but in the two-droplet problem there is a much richer variety of possible behaviours because any of the four contact lines may, in principle, move independently of the other three. In the present work, we consider four canonical behaviours, in each of which the droplets remain symmetric about the \(\xi \)-axis. Specifically, we consider the following modes of evaporation:
- 1.
The constant-inner-and-outer-contact-line (CIO) mode: the inner and outer contact lines of both droplets are pinned at \(\eta =\pm I(0)=\pm I_0\) and \(\eta =\pm \varOmega (0)=\pm \varOmega _0\) as the droplets evaporate.
- 2.
The constant-angle-centred (CAC) mode: both droplets evaporate with constant contact angle and remain centred at \(\eta =\pm (I+\varOmega )/2=\pm (I_0+\varOmega _0)/2\).
- 3.
The constant-angle and constant-inner-contact-lines (CAI) mode: both droplets again evaporate with constant contact angle, but now their inner contact lines are pinned at \(\eta =\pm I_0\).
- 4.
The constant-angle and constant-outer-contact-line (CAO) mode: both droplets again evaporate with constant contact angle, but now their outer contact lines are pinned at \(\eta =\pm \varOmega _0\).
3.3.1 Constant-inner-and-outer-contact-line (CIO) mode
In the CIO mode, the inner and outer contact lines of both droplets are pinned, \(I\equiv I_0\) and \(\varOmega \equiv \varOmega _0=I_0+2\). We may then immediately integrate (47) to obtain
Hence the lifetime of a pair of droplets evaporating in the CIO mode is
Figure 8 shows the evolution and the lifetime of a pair of droplets evaporating in the CIO mode. As for a single droplet in the CR mode, the contact angle \(\theta \) and the area \(A\) both decrease linearly with time t (Fig. 8a, b) and the lifetime \(t_\text {CIO}\) increases monotonically with \(\varPsi \) (Fig. 8c). In addition, since the shielding effect is weaker, and hence evaporation is faster, when the droplets are more widely separated, the lifetime \(t_\text {CIO}\) decreases monotonically with the separation between the droplets, \(2I_0\) (Fig. 8d).
3.3.2 Constant-angle (CAC, CAI, CAO) modes
In the CAC, CAI and CAO modes, the contact angle remains constant, \(\theta (t)\equiv \theta _0=1\). The three modes are distinguished by the different behaviours of the contact lines.
In the constant-angle-centred (CAC) mode, the droplets remain centred at \(\eta =\pm (I+\varOmega )/2=\pm (I_0+\varOmega _0)/2\). It is therefore convenient to write
where \(\varGamma =(I_0+\varOmega _0)/2\) is the position of the centre of the right-hand droplet and \(\varDelta (t)=(\varOmega -I)/2\) is its semi-width. We may then integrate (47) implicitly to obtain
Hence the lifetime of a pair of droplets evaporating in the CAC mode is
In the constant-angle-and-inner-contact-line (CAI) mode, the inner contact line is pinned, \(I\equiv I_0\). We may then integrate (47) implicitly to obtain
where
Hence the lifetime of a pair of droplets evaporating in the CAI mode is
In the constant-angle-and-outer-contact-line (CAO) mode, the outer contact line is pinned, \(\varOmega \equiv \varOmega _0\). We may then integrate (47) implicitly to obtain
whereFootnote 1
Hence the lifetime of a pair of droplets evaporating in the CAO mode is
Figure 9 shows the evolution and the lifetime of a pair of droplets evaporating in the three constant-angle modes. The difference between the modes is most clearly seen in Fig. 9a, which shows the inner and outer contact lines moving towards the centre of the droplet in the CAC mode, the outer contact line moving inward in the CAI mode, and the inner contact line moving outward in the CAO mode. Despite these differences, the evolution of the area \(A\), which decreases nonlinearly with t, is similar for all three modes (Fig. 9b). As in the CAI mode, the lifetimes \(t_\text {CAC}\), \(t_\text {CAI}\) and \(t_\text {CAO}\) increase monotonically with \(\varPsi \) (Fig. 9c) and decrease monotonically with the separation between the droplets, \(2I_{0}\) (Fig. 9d).
As Fig. 9c, d shows, the lifetimes of the three constant-angle modes are very similar, and it is only when the separation between the droplets is small that the difference between them becomes important. Specifically, Fig. 9d shows that the difference between \(t_\text {CAC}\), \(t_\text {CAI}\) and \(t_\text {CAO}\) becomes negligible when \(I_0 > rsim 5\) (i.e. when the droplet separation is several times the width of the droplets). As the droplets evaporate, the droplet separation is smallest in the CAI mode and largest in the CAO mode, resulting in the slowest evaporation, and hence the longest lifetime, in the CAI mode and the fastest evaporation, and hence the shortest lifetime, in the CAO mode. This point is further illustrated in Fig. 10, which shows the lifetimes \(t_\text {CAI}\), \(t_\text {CAC}\), \(t_\text {CAO}\) and \(t_\text {CIO}\) plotted as functions of \(\varPsi \). In particular, Fig. 10 shows that as the droplet separation increases the lifetimes of the three constant-angle modes (but not that of CIO mode) become virtually indistinguishable. We will discuss the latter behaviour in more detail in Sect. 3.3.3 below.
3.3.3 Asymptotic behaviour of the lifetimes in a large domain, \(\varPsi \gg I\)
The results shown in Fig. 10 motivate us to derive asymptotic expressions for the lifetimes of the droplets when \(\varPsi \gg I\). Noting the difference between closely-spaced and widely-separated droplets, we consider the regimes \(I_0 \ll 1\) and \(I_0 \gg 1\) separately.
In the regime \(I_0 \ll 1 \ll \varPsi \), the initial droplet separation is much less than the initial droplet semi-width. Equations (49), (52), (55) and (58) then yield
respectively.
On the other hand, in the regime \(1 \ll I_0 \ll \varPsi \), the initial droplet separation is much greater than the initial droplet semi-width. Equations (49), (52), (55) and (58) then yield
respectively.
Equations (59)–(66) show that, as in the one-droplet problem, in the regime \(\varPsi \gg I\), the lifetimes of all four modes depend logarithmically on the size of the domain.
When \(I_0 \ll 1\), Eqs. (59)–(62) show that the lifetimes depend on the mode of evaporation at \(O(1)\). The lifetime of the CIO mode is the shortest because, due to their pinned contact lines, the droplets in this mode have the largest surface area, and hence evaporate the fastest. Of the three constant-angle modes, the CAO mode has the shortest lifetime and the CAI mode the longest lifetime. This is because when the inner contact lines are pinned the droplets remain closer together and hence evaporate more slowly than in the CAC mode due to a stronger shielding effect, whereas when the outer contact lines are pinned the droplets move further apart and hence evaporate more quickly than in the CAC mode due to a weaker shielding effect.
On the other hand, when \(I_0 \gg 1\), Eqs. (63)–(66) show that the influence of the different behaviours of the contact lines on the three constant-angle modes is very weak and at \(O(1)\) the lifetimes of the CAC, CAI and CAO modes all coincide with each other, but differ from the lifetime of the CIO mode by \(1/(3\pi )\).
4 Comparison between the lifetimes of a single droplet and a pair of droplets
To provide further insight into the shielding effect the droplets have on each other, we compare the lifetimes of a single droplet and a pair of droplets in dimensional terms. For simplicity, we consider only the leading-order behaviour of the lifetimes in the regime of most interest in which the domain is large and approximately semi-circular with radius \(\hat{\varPsi }\).
Our reference point is the lifetime of a single droplet with initial semi-width \(\hat{R}_0\), which from (1), (37) and (38) is given by
A first natural comparison is with a pair of droplets, each with initial semi-width \(\hat{R}_0/2\), which together have the same total surface area as the single droplet (i.e. their initial separation is \(2\hat{I}_0=0\)). In this case the vapour concentration and flux are identical in the two problems. However, the lifetimes are not the same, because the cross-sectional area of the single droplet is twice the total cross-sectional area of the two droplets. Specifically, from (59)–(62) we obtain
so that, as expected, the lifetime of the pair of droplets is half that of the single droplet, i.e. \(\hat{t}_{\text {area}}\sim \hat{t}_{\text {single}}/2\).
Alternatively, we can consider the same total cross-sectional area of fluid, arranged either as two closely-spaced or two widely-separated droplets. In both cases the droplets have initial semi-width \(\hat{R}_0/\sqrt{2}\). If the droplets are closely spaced then from (59)–(62) we obtain
At leading order the lifetime of the pair of droplets is the same as that of the single droplet, but there is a negative \(O(1)\) correction because the two droplets have a larger total surface area than the single droplet. On the other hand, if the droplets are widely separated then from (63)–(66) we obtain
At leading order the lifetime of the pair of droplets is again the same as that of the single droplet, but now there is a larger negative \(O(\ln (\hat{I}_0/\hat{R}_0))\) correction due to a weaker shielding effect when the droplets are widely separated.
To illustrate these results we take the representative parameter values \(\hat{\rho }=998\) kg m\(^{-3}\), \(\hat{D} = 2.44\times 10^{-5}\) m\(^2\) s\(^{-1}\), \(\hat{c}_{\text {sat}} = 1.94\times 10^{-2}\) kg m\(^{-3}\) and \(\hat{c}_{\infty } = 7.76\times 10^{-3}\) kg m\(^{-3}\), corresponding to water at 295 K, evaporating into an environment with ambient vapour concentration \(\hat{c}_{\infty }=0.4\,\hat{c}_{\text {sat}}\) [11]. We further take \(\hat{\varPsi }=1\) m together with \(\hat{\theta }_0=0.1\) and \(\hat{R}_0 = 1\) mm, so that the droplet has a cross-sectional area of approximately \(6.7\times 10^{-8}\) m\(^2\).
With these parameter values, the timescale \(\hat{T} \approx 351\) s and the lifetime of a single droplet is \(\hat{t}_{\text {single}} \approx 567\) s. The lifetime of a pair of droplets with the same total surface area as the single droplet is \(\hat{t}_{\text {area}} \approx 283\) s. The lifetime of two closely-spaced droplets with the same total cross-sectional area as the single droplet is \(\hat{t}_{\text {close}} \approx 541\) s, whereas the lifetime of two widely-separated droplets with the same total cross-sectional area as the single droplet is \(\hat{t}_{\text {wide}} \approx 442\) s if the droplets are separated by \(2\hat{I}_0 = 2\) cm, and \(\hat{t}_{\text {wide}} \approx 356\) s if the droplets are separated by \(2\hat{I}_0 = 20 \hbox { cm}\).
5 Discussion and conclusion
In this contribution, we considered the diffusion-limited evaporation of thin two-dimensional sessile droplets either singly or in a pair. This two-dimensional problem is qualitatively different from the corresponding three-dimensional problem because, in contrast to in three dimensions, in two dimensions the size of the domain remains important even when it is much larger than the width of the droplets; it is therefore not possible to obtain a solution to the two-dimensional problem with a far-field boundary condition imposed “at infinity”. We therefore formulated a slightly modified problem in which the far-field condition was replaced by a relaxed condition at a distant, but finite, boundary. We then showed how a conformal-mapping technique may be used to calculate the vapour concentrations, and hence obtain closed-form solutions for the evolution and the lifetimes of the droplets in various modes of evaporation. These solutions demonstrate that 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 allowed us to quantify the shielding effect that the droplets have on each other, and how it extends the lifetimes of the droplets.
Although the present two-dimensional problem may be somewhat artificial, it has direct practical applications, for example to the inkjet printing of circuits [26], and may be realisable experimentally using a Hele-Shaw cell geometry. More fundamentally, it provides a rare opportunity to obtain a closed-form description of the behaviour of interacting droplets and to quantify the shielding effect. It therefore permits asymptotic and analytical insight into a class of problems of increasing scientific and industrial interest.
Notes
Note that, as a check on the correctness of (54) and (57), we may (up to an unimportant additive constant) recover \(F_\text {CAO}\) from \(F_\text {CAI}\) by replacing \(\varOmega \) with \(I\), \(I_0\) with \(\varOmega _0\) and \(\varPsi \) with \(\mathrm{i}\varPsi \) in (54), and \(F_\text {CAI}\) from \(F_\text {CAO}\) by replacing \(I\) with \(\varOmega \), \(\varOmega _0\) with \(I_0\) and \(\varPsi \) with \(\mathrm{i}\varPsi \) in (57).
References
Brutin D, Starov V (2018) Recent advances in droplet wetting and evaporation. Chem Soc Rev 47(2):558–585
Deegan RD, Bakajin O, Dupont TF, Huber G, Nagel SR, Witten TA (1997) Capillary flow as the cause of ring stains from dried liquid drops. Nature 389(6653):827–829
Routh AF (2013) Drying of thin colloidal films. Rep Prog Phys 76(4):046603
Larson RG (2014) Transport and deposition patterns in drying sessile droplets. AIChE J 60(5):1538–1571
Giorgiutti-Dauphiné F, Pauchard L (2018) Drying drops. Eur Phys J E 41:32
Lohse D, Zhang X (2015) Surface nanobubbles and nanodroplets. Rev Mod Phys 87(3):981–1035
Stauber JM, Wilson SK, Duffy BR, Sefiane K (2014) On the lifetimes of evaporating droplets. J Fluid Mech 744:R2
Stauber JM, Wilson SK, Duffy BR, Sefiane K (2015) On the lifetimes of evaporating droplets with related initial and receding contact angles. Phys Fluids 27(12):122101
Schofield FGH, Wilson SK, Pritchard D, Sefiane K (2018) The lifetimes of evaporating sessile droplets are significantly extended by strong thermal effects. J Fluid Mech 851:231–244
Popov YO (2005) Evaporative deposition patterns: spatial dimensions of the deposit. Phys Rev E 71(3):036313
Dunn GJ, Wilson SK, Duffy BR, David S, Sefiane K (2009) The strong influence of substrate conductivity on droplet evaporation. J Fluid Mech 623:329–351
Dunn GJ, Wilson SK, Duffy BR, David S, Sefiane K (2008) A mathematical model of the evaporation of a thin sessile liquid droplet: comparison between experiment and theory. Colloids Surf A 323(1–3):50–55
Dunn GJ, Wilson SK, Duffy BR, Sefiane K (2009) Evaporation of a thin droplet on a thin substrate with a high thermal resistance. Phys Fluids 21(5):052101
Sneddon IN (1966) Mixed boundary value problems in potential theory. Cambridge University Press, North-Holland
Jackson JD (1999) Classical electrodynamics, 3rd edn. Wiley, Hoboken
Mathews JH, Howell RW (2012) Complex analysis for mathematics and engineering. Jones & Bartlett Publishers, Burlington
Stoker JJ (1947) Surface waves in water of variable depth. Q Appl Math 5(1):1–54
Duffy DG (2008) Mixed boundary value problems. CRC Press, Boca Raton
Weber H (1873) Ueber die Besselschen Functionen und ihre Anwendung auf die Theorie der elektrischen Ströme. Journal für die reine und angewandete Mathematik 75:75–105
Whiteman JR (1968) Treatment of singularities in a harmonic mixed boundary value problem by dual series methods. Q J Mech Appl Mech 21(1):41–50
England AH (1979) Mixed boundary-value problems in potential theory. Proc Edinb Math Soc 22(2):91–98
Yang F, Li JCM (1993) Impression creep of a thin film by vacancy diffusion. I. Straight punch. J Appl Phys 74(7):4382–4389
Yang F, Li JCM (1993) Impression creep of a thin film by vacancy diffusion. II. Cylindrical punch. J Appl Phys 74(7):4390–4397
Yang F, Yao R (1996) The solution for mixed boundary value problems of two-dimensional potential theory. Indian J Pure Appl Math 27(3):313–322
Laporte O, Fowler RG (1966) Resistance of a plasma slab between juxtaposed disk electrodes. Phys Rev 148(1):170–175
Yarin AL, Szczech JB, Megaridis CM, Zhang J, Gamota DR (2006) Lines of dense nanoparticle colloidal suspensions evaporating on a flat surface: formation of non-uniform dried deposits. J Colloid Interface Sci 294(2):343–354
Lacasta AM, Sokolov IM, Sancho JM, Sagués F (1998) Competitive evaporation in arrays of droplets. Phys Rev E 57(5):6198–6201
Schäfle C, Bechinger C, Rinn B, David C, Leiderer P (1999) Cooperative evaporation in ordered arrays of volatile droplets. Phys Rev Lett 83(25):5302–5305
Kokalj T, Cho H, Jenko M, Lee LP (2010) Biologically inspired porous cooling membrane using arrayed-droplets evaporation. Appl Phys Lett 96(16):163703
Sokuler M, Auernhammer GK, Liu CJ, Bonaccurso E, Butt HJ (2010) Dynamics of condensation and evaporation: effect of inter-drop spacing. Europhys Lett 89(3):36004
Carrier O, Shahidzadeh-Bonn N, Zargar R, Aytouna M, Habibi M, Eggers J, Bonn D (2016) Evaporation of water: evaporation rate and collective effects. J Fluid Mech 798:774–786
Castanet G, Perrin L, Caballina O, Lemoine F (2016) Evaporation of closely-spaced interacting droplets arranged in a single row. Int J Heat Mass Transf 93:788–802
Shaikeea A, Basu S, Hatte S, Bansal L (2016) Insights into vapor-mediated interactions in a nanocolloidal droplet system: evaporation dynamics and affects on self-assembly topologies on macro-to microscales. Langmuir 32(40):10334–10343
Shaikeea AJD, Basu S (2016) Evaporating sessile droplet pair: insights into contact line motion, flow transitions and emergence of universal vaporisation pattern. Appl Phys Lett 108(24):244102
Shaikeea AJD, Basu S (2016) Insight into the evaporation dynamics of a pair of sessile droplets on a hydrophobic substrate. Langmuir 32(5):1309–1318
Hatte S, Pandey K, Pandey K, Chakraborty S, Basu S (2019) Universal evaporation dynamics of ordered arrays of sessile droplets. J Fluid Mech 866:61–81
Khilifi D, Foudhil W, Fahem K, Harmand S, Ben Jabrallah S (2019) Study of the phenomenon of the interaction between sessile drops during evaporation. Therm Sci 23(2B):1105–1114
Wray AW, Duffy BR, Wilson SK (2020) Competitive evaporation of multiple sessile droplets. J Fluid Mech 884:A45
Greenwood JA (1966) Constriction resistance and the real area of contact. Br J Appl Phys 17(12):1621–1632
Gladwell GML, Fabrikant VI (1982) The interaction between a system of circular punches on an elastic half space. J Appl Mech 49(2):341–344
Fabrikant VI (1985) On the potential flow through membranes. Zeitschrift für angewandte Mathematik und Physik 36(4):616–623
Argatov I, Sevostianov I (2009) On relations between geometries of microcontact clusters and their overall properties. Int J Eng Sci 47(10):959–973
Argatov II (2011) Electrical contact resistance, thermal contact conductance and elastic incremental stiffness for a cluster of microcontacts: asymptotic modelling. Q J Mech Appl Mech 64(1):1–24
Argatov II (2012) The contact problem for a periodic cluster of microcontacts. J Appl Math Mech 76:604–610
Dollet B, Lohse D (2016) Pinning stabilizes neighboring surface nanobubbles against Ostwald ripening. Langmuir 32(43):11335–11339
Laghezza G, Dietrich E, Yeomans JM, Ledesma-Aguilar R, Kooij ES, Zandvliet HJW, Lohse D (2016) Collective and convective effects compete in patterns of dissolving surface droplets. Soft Matter 12(26):5787–5796
Michelin S, Guérin E, Lauga E (2018) Collective dissolution of microbubbles. Phys Rev Fluids 3(4):043601
Xie Q, Harting J (2019) The effect of the liquid layer thickness on the dissolution of immersed surface droplets. Soft Matter 15(32):6461–2468
Sefiane K, Wilson SK, David S, Dunn GJ, Duffy BR (2009) On the effect of the atmosphere on the evaporation of sessile droplets of water. Phys Fluids 21(6):062101
COMSOL Multiphysics v. 5.4. http://www.comsol.com. COMSOL AB, Stockholm, Sweden
Acknowledgements
F.G.H.S. acknowledges the financial support of the UK Engineering and Physical Sciences Research Council (EPSRC) via EPSRC Doctoral Training Partnership Grant EP/N509760/1, the University of Strathclyde, and the University of Edinburgh. We are grateful to Dr Thomas Beck (MIT) for helpful discussions about the one-droplet problem, and to Dr Brian R. Duffy (University of Strathclyde) and Prof. Khellil Sefiane (University of Edinburgh) for insightful discussions about many aspects of droplet evaporation.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Schofield, F.G.H., Wray, A.W., Pritchard, D. et al. The shielding effect extends the lifetimes of two-dimensional sessile droplets. J Eng Math 120, 89–110 (2020). https://doi.org/10.1007/s10665-019-10033-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10665-019-10033-7