On the lifetimes of two-dimensional droplets on smooth wetting patterns

We study the evolution and lifetime of droplets evaporating on a smooth chemical pattern, which is characterised by a spatially varying contact angle. We formulate a model that combines the evaporation rate of the droplet for a given volume with the static stability of the droplet as the volume changes in time quasi-statically. We derive an exact equation for the evaporation rate that is studied analytically under the limiting cases of nearly neutral wetting and highly hydrophilic conditions. We find that the evaporation rate of the droplet is highly dependent on the size and shape of the fictitious infinity where a far-field boundary conditions needs to be applied. We also study how the droplet’s lifetime depends on the averaged contact angle and strength ε\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varepsilon $$\end{document} of the chemical pattern, observing that the lifetime of the droplet is maximised for a droplet with average contact angle π/2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi /2$$\end{document} and with ε≲0.1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varepsilon \lesssim 0.1$$\end{document}.


Introduction
There has been a large amount of research devoted to understanding the motion of sessile droplets on solid surfaces.
A key aspect of this problem is to quantify the effects of patterned surfaces on the droplet's motion. For the case of non-evaporating droplets see, e.g. [1][2][3][4][5][6], and for droplet evaporation process see, for example, [7][8][9][10][11]. The breadth of this research has been driven by the wide variety of real-life situations in which a sessile droplet can occur, as well as a large amount of industrial applications, such as for example, in the drying of integrated circuits [12], the spreading of metallic inks that are used in inkjet printing [13,14], or pattern deposition from particulates in a droplet [15,16].
In the absence of an external flow and temperature gradients, droplet evaporation is driven purely by diffusion. Under these conditions, the model of droplet evaporation must also describe the vapour transport from the droplet to the surrounding atmosphere [7]. The vapour profile surrounding a static droplet is harmonic and has direct similarities to other fields, such as for example electric potential fields in electrostatic studies [17], fluid flow in porous materials [18], and dissolution of bubbles [19,20]. There are many variables which influence the evaporation rate of the droplet. The relationship between the contact angle and the evaporation rate has been studied in the threedimensional problem by making use of toroidal coordinates to study the lifetime of droplets under various conditions, for example in the stick-slip evaporation mode [21] or in the strongly hydrophobic limit [22]. These results have also been extended to include a surrounding shear flow [23].
In two dimensions, the difficulty occurs from the appropriate application of the far-field condition, i.e. that sufficiently far away from the droplet the presence of the droplet itself does not affect the vapour concentration. Indeed recently, an eloquent proof highlighted that in two dimensions, the far-field condition must be applied at a large but finite distance from the droplet to permit the existence of a solution [24]. To ensure that a solution exists, we shall impose a relaxed boundary condition following [25].
The classical understanding of droplet evaporation on a regular solid surface is that a droplet evaporates following a stick-slip model [26,27], also referred to as stick-slide [21,28], which has two phases -a stick phase and a slide phase that alternate from one another. During the stick phase, the droplet contact lines remain pinned to conserve the contact area and the contact angle decreases to accommodate the decreasing volume. In contrast, the slide phase conserves the contact angle while the contact line traverses. In addition, a second mode known as stick-jump has been recently reported, where intermittent contact line pinning causes the droplet to alternate between sticking and jumping [29][30][31]. It is generally well accepted that the emergence of stick-slip and stick-jump phenomena is a consequence of contact line pinning, which is itself a result of microscopic roughness or chemical defects of the solid surface, and hence generally difficult to control.
The study of droplet evaporation on smooth pinning-free surfaces has been sparked by the production of slippery surfaces, such as Slippery Liquid Infused Porous Surfaces [32][33][34] and Slippery Omniphobic Covalently Attached Liquid Like Surfaces [35]. These surfaces are characterised with a very low contact angle hysteresis and can be designed with topographical or wetting patterns [36]. The theory behind droplets atop smooth surface has been described in [37]. In addition, it has recently been shown that a droplet evaporating on a slippery surface that has a smooth macroscopic pattern follows a reproducible sequence of events that are characterised by bifurcations on the droplet's shape and position, and it has been referred to as snap evaporation, see [38] for topographical patterns and [39] for chemical patterns. In contrast to the stick-jump mode reported in [29], during snap evaporation, the droplet's contact line is never fully stationary but it moves very slowly except when a shape bifurcation occurs that forces the droplet to shift laterally and/or have a reduction of its base radius. In particular, the droplet may lose stability at certain critical values of its volume through snap-through and pitchfork bifurcations giving rise to dynamic snap events. A key assumption of this approach is a quasi-static evolution so that the state of the droplet is given by its equilibrium properties. Such approach has been used to quantify droplets evaporating on either smooth topographies [38] or chemical patterns [39], as well as for droplets with mass transfer moving on chemical patterns [40]. However, while all these works have focussed on the control of droplet motion on smooth patterns, little is known on the effect of patterning on the lifetime of evaporating droplets.
The primary focus of this paper is to investigate how the lifetime of the droplet is affected by the presence of a chemical smooth pattern, which is characterised by a spatially varying contact angle. Here, we propose a model that combines diffusion-limited evaporation with bifurcation theory. In particular, we formulate a framework which takes into account the evaporation rate of the droplet for a given shape while assuming that the shape and location of the droplet on the surface are dictated by the droplet's static stability as the volume changes in time. Importantly, we examine the stability from a quasi-static viewpoint and assume that the droplet moves from one stable solution to another significantly faster than any speed associated with the evaporation of the droplet. The full model yields an evolution equation for the droplet's volume which depends on the current configuration. Analysis of the evolution equation allows us to study both the evolution of the droplet and ultimately calculate the lifetime of the droplet for different wetting configurations. This paper is organised as follows. In Sect. 2, we formulate the problem mathematically and define the key parameters of the system. Section 3 goes on to investigate the evaporation rate of a static droplet using a previously documented conformal mapping technique [25]. The evaporation rate is then tested against an asymptotic approx-imation for a near neutral wetting droplet and the limiting case of a near vanishing contact angle [24]. Section 4 summarises the key results which are used to determine the stability of a droplet on a smooth patterned surface. The lifetime of a droplet is then examined in Sect. 5 by using the results from Sects. 3 and 4. We conclude in Sect. 6.

Problem formulation
Consider a two-dimensional droplet of liquid which has surface tension γ and density ρ at rest on a chemically patterned substrate. We define a Cartesian coordinate system, (x, y), with the x-axis running along the substrate and the y-axis normal to the surface. The droplet wets a region on the substrate which has half-width w and the centre of the droplet is located at x = d. In the absence of gravity, the equilibrium shape of the droplet is described by a circular segment with contact angle θ , and the volume of the droplet is given by We also assume that initially d = 0, i.e. the droplet is aligned with the centre of the coordinate system. Figure 1 shows a sketch of the droplet, where Θ(x) is the equilibrium contact angle imposed by the solid substrate. We consider periodic chemical patterns that lead to a spatial variation of the equilibrium contact angle, Θ(x), of the form where θ 0 is the average contact angle, ε is the ratio of the pattern's amplitude to its average contact angle, i.e. θ 0 ε is the amplitude of the pattern, and k = 2π/λ, where λ is the wavelength of the pattern. Here, we assume that typical length scale λ is comparable to the droplet's footprint. To ensure that the contact angle has physical relevance, we impose the constraint Θ(x) < π everywhere and hence we find bounds for ε, i.e.
Note that at equilibrium the contact angle of the droplet is given by θ = Θ(d ± w) with a volume given by Eq. (1). It is also important to remark, that while, in this paper, we focus on sinusoidal patterns, the methods presented here are applicable to any patterning that has an inherent length scale. We model the evaporation of the droplet as a quasi-static, diffusion-limited process and assume that the surrounding fluid is stationary. The quasi-static condition here means that the evaporation time scale is assumed to be much slower than the droplet's relaxation time scale, so that we can use a semi-circular approximation for the w x Θ(x) λ θ Fig. 1 A two-dimensional droplet resting on a flat solid surface with a smooth chemical pattern that is quantified in terms of the equilibrium contact angle Θ(x). The wavelength of the pattern is λ and the droplet's footprint is denoted by w. The circular shape of the droplet is quantified in terms of the contact angle θ, which is the same on both contact points droplet's shape. The evaporation rate is found by calculating the mass flux of liquid vapour across the boundary of the droplet. To this end, we will find first the liquid vapour profile in the surrounding fluid (although we note that in general this is not always necessary [41]). Under these assumptions, the liquid vapour concentration, c(x, y), is harmonic in the surrounding fluid and has a saturated value, c sat , on the boundary of the droplet. In the far-field limit, the liquid vapour concentration has a constant value that we denote c ∞ . To ensure that a solution exists, we apply the far-field condition on a finite curve [25]. Specifically, we enforce this condition on the perimeter of the upper-semi-circle with radius R w. In summary, the concentration of liquid vapour satisfies where Ω is the region of R 2 bounded by the curves y = 0, y = h(x) and y = √ Here, y = h(x) describes the droplet's interface, and (4d) is the standard no flux boundary condition applied on the solid surface. From the liquid vapour profile in the surrounding fluid, the evaporation rate, i.e. the volume change rate, is then calculated through where n is a coordinate normal to the surface and D represents the diffusion coefficient. It is convenient to nondimensionalise the system by defining the following variables: where λ is the wavelength of the patterning, i.e. λ = 2π/k, and the starred variables are all dimensionless. Henceforth, the stars will be dropped and all variables should be considered dimensionless, unless otherwise stated. With these new variables, the Laplace system (4) becomes and Eq. (5) for the evaporation rate becomes In the following, we find an exact equation which gives the evaporation rate of a droplet for a given (dimensionless) width and contact angle. While in Sect. 4, we shall study the stability of a droplet for a given (dimensionless) volume resting on a substrate with local contact angle given by (2). These results will allow us to form a first-order ordinary differential equation to describe the evolution of the droplet and hence determine its lifetime.

Evaporation rate
In the first part of this section, we summarise the work of [25], in which an exact equation for the droplet's evaporation rate is derived in terms of the droplets half-width, w, and contact angle, θ . We shall then investigate special cases where the evaporation rate can be written explicitly. We introduce the local inclination angle, φ ∈ [−π, π], which is defined as the angle formed against the horizontal, to describe the surface of the droplet as We note that this formulation is applicable for all contact angles, unlike a direct formulation using h(x). Since the shape of the droplet is a circular segment, we use the conformal mapping where z = x + iy is the complex variable in the physical plane, σ = 2(1 − θ/π), a = w/σ , and ζ is the complex variable in the transformed plane, see Fig. 2. We note that (10) maps the z-plane onto the ζ -plane, and in particular, a circular arc of width w and contact angle θ is mapped to the perimeter of the upper half circle with radius a. In the transformed space, we can solve the Laplace system (7) to find the complex potential, χ , (which satisfies Re(χ ) = c) given by Taking the derivative of χ with respect to z, we find the mass flux across the surface Fig. 2 Diagram showing how a circular segment with angle θ and width w is mapped to the upper semi-circle with radius a = w/σ using (10), where σ = 2(1 − θ/π). The grey shaded region in the left diagram is a depiction of Ω, and in the right diagram, it depicts the transformed domain where g(z) = (z − w)/(z + w). It follows that the normal derivative of the concentration is given by wheren x ,n y are the x, y components of the outward unit normal to the droplet's interface. Hence, we can use (13) to find the evaporation rate, (8), which, in these transformed coordinates, is given as with x(φ) and h(φ) given in (9). Equation (14) is the principle result of this section. It is an exact equation that predicts the evaporation rate of a droplet for a given width w and contact angle θ . In the following, we study two special cases for which it is possible to find asymptotic solutions of Eq. (14).

Limiting cases
We start by considering the case of nearly neutral wetting, i.e. the droplet interface is close to a semi-circle with the contact angle satisfying |θ − π/2| 1. Observe that the surface of the droplet is then described, in polar coordinates (r, ϕ), by where ε = θ − π/2. We consider a regular expansion for the concentration in terms of ε. By solving the resulting leading and first-order systems, we find that the concentration is given by (see Appendix A for a derivation) By noting that the outward unit normal vector satisfieŝ and applying Eq. (16) to (8) we find which gives the evaporation rate of a nearly neutral wetting setting. The second special case we consider is the extremely hydrophilic droplet, i.e. a droplet with asymptotically small contact angle. This case has been examined in [24], and here, we give a brief summary of that method.
In this method, the far-field condition is not applied on the boundary of a semi-circle; instead, an upper semiellipse is used. This choice of fictitious infinity is to allow a sufficiently tractable conformal mapping solution. Using the upper semi-ellipse with major-axis √ R 2 + 1 that is aligned with the x-axis and minor-axis R as the fictitious infinity, i.e.
we see that the domain is mapped to the rectangular domain (u, v) ∈ (−1, 1) × (0, S) by the conformal mapping where ζ = u + iv and S satisfies R = 2 sinh π S w . In the transformed space, we obtain the complex potential This result is then mapped back to the z-plane and applied to (8) to find

Numerical results for homogeneous surfaces
We present here numerical results of Eq. (14) for a given and constant contact angle. In particular, in this section, we examine the instantaneous evaporation rate for a defined droplet, while the examination of a droplet lifetime is discussed in Sect. 5. We note that to fully determine the solution of the droplet, we need to either fix the volume V of the droplet or its width w. We consider each of these two cases in turn. We start by fixing the width of the droplet to w = 1. As we vary the contact angle, the droplet's volume is then given by Eq. (1). Figure 3a shows numerical results of Eq. (14) for different values of R. The blue and red lines represent the results from the asymptotic limiting cases given by Eq. (22) (nearly complete wetting) and (18) (nearly neutral wetting), respectively. We see that in both cases, there is excellent agreement between numerical and asymptotic results. We note that for the case of θ = 0, the asymptotic result becomes more accurate as R increases (see inset panel of Fig. 3a). This is because in the analytical treatment of the zero contact angle case, the far-field conditions are applied on the elliptic curve (19) which tends to the circular arc of radius R as R → ∞. We also see that the asymptotic results derived here for nearly neutral wetting meet the exact results at a tangent to the point θ = π/2, as expected [24] .  Figure 3b shows the evaporation rate as a function of the contact angle when the volume of the droplet is fixed. Under this condition, the droplet's width is given by The red lines in Fig. 3b represent the asymptotic analysis given by Eq. (18), where we can see, again, that the two lines meet in tangent at θ = π/2. In order to conserve the volume as θ → 0, we observe that w → +∞, and therefore, the extremely hydrophilic case, given by (22), predicts that the volume change rate tends to negative infinity. This agrees with the numerical results shown in Fig. 3b.

Stability analysis: bifurcation diagrams
As is discussed in Sect. 1, the evolution of a droplet evaporating under quasi-static conditions can be described by its equilibrium properties. As the droplet slowly evaporates the shape and location of the droplet can be modelled by studying the stability of the equilibrium solutions. This approach has recently been used to describe droplets evaporating on chemical patterns [39] and droplet with mass transfer moving on chemical patterns [40]. This section reviews the work presented in these references to study the static stability of a droplet with known volume on a smooth patterned surface. This analysis will be used later in our evolution model, where we assume that the droplet will only follow stable points. Further, we shall assume that the droplet moves between stable solutions instantaneously, as discussed in the introduction. It is important to emphasise that the process of droplet evaporation is an intrinsically dynamic situation, which in general would require a dynamic stability analysis. However, here we are assuming that the evaporation process is sufficiently slow (quasi-static) such that we can use static stability analysis to effectively predict the evolution of the droplet. To ensure we have an equilibrium configuration the forces along the substrate must balance hence, and we define this contact angle to be θ . This relation shows that for sinusoidal patterns the droplets centre aligns with a maximum or minimum of the patterning, or that the width of the droplet (i.e. 2w) is a integer multiple of the pattern wavelength. It has been shown that only the droplets which align with a maximum or minimum of the patterning can be stable [39,40]. In addition, the stability of those solutions depend on the droplet's volume and solutions that are centred on a maximum or minimum of the patterning can be stable, partially stable (saddle node) or fully unstable. The term partially stable is used, in this context, to describe the stability to symmetric perturbations but unstable to asymmetric perturbations. Indeed, for small volumes or small ε all droplets aligned with a maximum or minimum of the patterning are stable to symmetric perturbations. As ε increases, solutions undergo a pitchfork bifurcation to lose their stability to symmetric perturbations [39]. The stability of a two-dimensional droplet can be determined from the interfacial energy, which in its dimensionless form is where the first term of the RHS describes the contribution from the liquid/gas interface and the second term quantifies the solid/fluid interfacial contribution to the energy. In the above equation, θ is the contact angle that satisfies (24) and, for a given volume V , is related to the half-width w through Eq. (1). The stability of a given (V, d, w) triple is then found by numerically evaluating the Hessian matrix of the energy, and determining the number of positive eigenvalues. If the Hessian matrix has only positive eigenvalues then the system is stable, while if there exists a negative eigenvalue then the system is unstable. By tracing the stable and unstable equilibrium positions as we change the volume, we are able to construct bifurcation diagrams, as shown in Fig. 4. Panel (a) shows two branches of solutions that correspond to droplets aligned with a minimum of the chemical pattern (blue curves) and a maximum (green curves). In both cases, the red regions denote solutions that are saddle nodes (i.e. they are unstable against asymmetric perturbations only). The key feature of these diagrams is the emergence of pitchfork bifurcation points at some critical volumes, which are denoted as black solid points in panel (a). Therefore, as the droplet evaporates and the volume decreases in time, the presence of asymmetric perturbations in the system can trigger lateral displacement whereby the droplet shifts from a maximum to a minimum of the chemical pattern, or vice versa (see black arrows in panel (a)). This has been referred to as snap evaporation [38]. In our model, we will assume that such transitions occur instantly.
When we increase the strength of the chemical pattern, ε, the stability diagrams show the existence of saddle-node bifurcations, which are denoted as empty circles in Fig. 4b. Under these conditions, if the system is not affected by any perturbation, it may be possible that, as the droplet evaporates along a stable branch, the droplet bypasses the pitchfork bifurcation and reaches the critical turning point of the saddle-node bifurcation. At this point there is no further solutions and the droplet radius rapidly decreases while the droplet remains at the same location. This is illustrated by the arrow in Fig. 4b.
In the following section, we use these ideas to construct the stability path of a droplet evaporating on a smooth pattern.

Droplet lifetimes
In this section, we investigate the lifetime of a droplet resting on a smooth patterned surface. As the droplet evaporates, it evolves according to the model described in Sect. 4, while the rate at which it evaporates is calculated from (14). We also assume that the shape of the droplet is always given by its equilibrium configuration with both contact angles equal. Hence, the droplet is assumed to evolve according to where at each time step, d and w correspond to an equilibrium solution for the corresponding volume. The evolution is then found by integrating (26) with the initial condition V (t = 0) = V 0 , and it is concluded at t = τ , when V < tol, for some small numerical tolerance, 0 < tol 1, which defines the lifetime τ of the droplet. Here, we have used tol = 10 −4 V 0 but other values were considered observing no effect in the results. We shall study different wetting cases, which lead to different evaporation modes.

No patterning (ε = 0)
We start by considering a substrate without a patterning, i.e. ε = 0. In this case we are simply in the classical slip evaporation mode. The lifetime of a droplet τ is depicted in Fig. 5 as a function of the contact angle. We see that the lifetime of the droplet, away from the extremely hydrophobic case, matches the general pattern of the evaporation rate, shown in Fig. 3b. In the extremely hydrophobic region, we see that the lifetime of the droplet increases dramatically with θ to a finite value as θ → π . We note that this is consistent with the results reported in [22] for highly hydrophobic axisymmetric droplets, although in that case the asymptotic limit θ → π is approached from above. The results presented in Fig. 5 show that the lifetime of a droplet is mostly independent of θ under partial wetting conditions, but it dramatically reduces or increases as the droplet becomes highly hydrophilic or highly hydrophobic, respectively.
We can compare these results with the lifetimes that have been predicted elsewhere under special conditions. First consider the case of the extremely hydrophilic droplet (θ 1) that was first detailed in [24]. In this case the volume of the droplet can be accurately approximated by V = 2w 2 θ/3. Hence, we obtain the differential equation (noting that θ is constant) This differential equation can be integrated implicitly over t ∈ (0, τ ), where τ gives the lifetime of the droplet, to yield which agrees with the result reported in [24]. The second case of interest is when the droplet has contact angle π/2. Under these conditions we obtain the differential equation where V = πw 2 /2. This, again, may be integrated implicitly to obtain which is consistent with the result obtained in [21,22,28] for axisymmetric 3D droplets, except for the dependency on the fictitious infinity R, which is not required in the 3D problem. Figure 5 shows the lifetime of the droplet as a function of the contact angle using the exact equation and the two special cases we have discussed here. We see that the exact equation matches the special cases in their respective regimes. Indeed, the small angle approximation agrees well with the exact form for contact angles less than 0.05π , and the hydroneutral solution gives a good approximation for contact angles between 0.4π and 0.8π .

Sinusoidal patterning (ε = 0)
We next consider the case with a smooth chemical pattern, ε = 0, and we will distinguish between two different modes. Firstly, consider a droplet which is unable to shift its position as it loses stability. This can be thought of as an analogue to a 3D droplet at the centre of a concentric rings patterning. In this case, the droplet will evaporate along the corresponding branch in the bifurcation diagram, but because the centre is fixed, the droplet will not snap to a new position when it loses stability (see Fig. 4b).
In the second case, we allow the droplet to shift at each pitchfork bifurcation leading to lateral snaps along the droplet's position (see Fig. 4a). Figure 6 compares the evolution of these two different cases alongside the homogeneous case (ε = 0). Panels (a) and (b) show the corresponding evolution of the droplet's width as a function of the volume, where we can see that when the droplet is fixed, the droplet's width follows a single bifurcation branch (see panel (a)), but when the droplet is free to traverse the surface it snaps at the pitchfork bifurcation point (see panel (b)). Panel (c) shows the time evolution of the volume for all cases. Interestingly, the volume of the droplet, given as a function of time, does not change significantly, and in particular, the droplet on the homogeneous surface has a marginally longer lifetime, while the freely moving droplet has a marginally shorter lifetime. This can be explained by considering Fig. 3a, in the region around θ = π/3. We see that the evaporation rate is an increasing function around θ = π/3. Hence, time-averaged contact angle for the droplets on the patterned surfaces must be greater than the spatially-averaged contact angle. Further, the droplet which is free to traverse the surface is able to maintain a higher contact angle for longer than the droplet which is fixed in position. We now consider how the droplet's evolution is affected by changing ε. In Fig. 7a we compare the lifetime of the droplet under the two regimes, as a function of ε, for various θ 0 . It is shown that increasing the patterning amplitude decreases the lifetime of the droplet. We also note that the difference between both regimes increases as ε 0.1, which can be attributed to the fact that above this value, there is the existence of saddle-node bifurcations that affect the dynamics of the droplet. It should be further noted that as the droplet becomes more hydrophobic (but not in the extreme region) the lifetimes drop, but follow the same trends as shown. Interestingly, there is a value of θ 0 between π/2 and 5π/12 where the droplet lifetime switches from being longer for a static droplet to being longer for a freely moving droplet. We also study how the time evolution of the droplet depends on the value of the outer boundary R, which has been shown to have an effect in similar type of two-dimensional problems [24,25] . In Fig.  7b we can see that increasing R increases the lifetime of the droplet, as expected. We also compare the evolution of the droplet under the two evaporation settings, where we note that the lifetime of the droplet varies marginally between the two regimes (difference between solid and dashed lines in Fig. 7b).

Discussion
We have investigated the evolution and lifetime of an evaporating two-dimensional sessile droplet resting atop a smooth patterned surface. The evaporation rate was calculated following the approach introduced in [25], which was then used to obtain an asymptotic expression for the limiting case of nearly neutral wetting. We also considered the limiting case of highly hydrophilic droplets, which was in agreement with the result reported in [24].
We found that the evaporation rate of the droplet is highly dependent on the size and shape of the fictitious infinity where the far-field boundary conditions are applied. Indeed, if the ratio of the radii of the fictitious infinity and the droplet is held constant, i.e. R = w(t) × const., then the evaporation rate remains almost constant. However, it makes physical sense to hold R constant as the droplet evaporates, rather than keeping the ratio of R and w constant. This dependence on the fictitious infinity was also reported in [24].
To estimate the lifetime of the droplet, we assumed that the evolution of the droplet's radius, w, as the volume V varies is dictated by the wetting equilibrium properties. This is quantified in terms of bifurcation diagrams w(R) which show the emergence of pitchfork and saddle-node bifurcations where solutions change stability [39,40]. Such hierarchy of bifurcations allows us to distinguish between two distinctive droplet's time evolution. In one case, the droplet is allowed to shift translationally at each pitchfork bifurcation, and on the other case, the location of the droplet is fixed but its width undergoes a transition at each saddle-node bifurcation. Hence, given a given chemical pattern, the time evolution of the volume is constructed by combining the information from the evaporation rate and the bifurcation diagrams.
We found that the lifetime of the droplet depends subtly on the contact angle, θ 0 , the parameter ε and the droplet's evolution path. Indeed, it was shown that the lifetime of the droplet is maximised for a droplet with average contact angle π/2 and with ε 0.1. This result is expected since the evaporation rate of the droplet is at its slowest for θ = π/2, as shown in Fig. 3b, and the contact angle is at this critical value for longest with ε = 0. Our results also showed that the droplet's lifetime always decreases with the strength ε of the pattern, indicating that droplets evaporate faster on wetting patterns than on homogeneous surfaces. It was further found that the two evaporation scenarios predict near equal lifetimes for ε 0.1 while when ε exceeds this value the two lifetimes begin to diverge. The case for which the droplet shifts laterally as it evaporates (snap evaporation) was found to have shorter lifetimes for droplets away from the nearly neutral case, whereas in the nearly neutral case, this snap mode increased the lifetime of the droplet.
It was shown that the lifetime does not change significantly as the contact angle changes (but remains away from the extreme cases), and the lifetime similarly only subtly changes with ε. The lifetime can be bound above (and approximated) by the case where the contact angle is always π/2, which yields or in dimensional terms, the lifetime of the droplet, T , may be bounded by where λ, D, c sat and c ∞ are the wavelength of the surface patterning, the diffusive coefficient, the saturated concentration of liquid vapour and natural concentration level in the surrounding fluid of the liquid vapour without the presence of a droplet, respectively. The ideas presented in this work could be used to explore three-dimensional configurations, although it should be noted that unlike the two-dimensional case presented here, a three-dimensional droplet on a patterned surface does not necessarily take a simply described shape, even in the absence of gravity. This presents a significant challenge to overcome in order to calculate the evaporation rate and then apply the result to find the lifetime of the droplet. In addition, alternative patterns could be considered, e.g. sharper variations on the contact angle that would allow to study transitions from slip to pinned configurations. We intend to address these and related issues in the future studies.
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/.
where k is a real constant. We see that the non-homogeneous boundary condition on r = w has been neglected from this system and shall be applied later. By solving these ordinary differential equations, we find where A 0 , A 1 , A 2 , . . . are constants, which are fixed by applying the condition on r = w. This condition becomes To then calculate A 0 , A 1 , A 2 , . . ., we multiply this condition by cos(mφ) for m = 0, 1, 2, . . . and integrate the resulting system over ϕ ∈ (0, π). This yields By combining these results, we obtain the desired result