Evaporation of Small Sessile Drop Deposited on Horizontal Solid Surface: New Exact Solutions and Approximations

1. Introduction

Evaporating liquid sessile drop is an important object both for theoretical investigations (evaporation dynamics, microfluidics, self-organization of solute, etc.) and wide variety of applications (printing technologies for functional coatings, medical diagnostics, food science, geophysics, etc.) [1,2,3,4,5,6,7,8,9,10].

The model of diffusion-limited quasi-stationary evaporation of the spherical droplet in quiescent air environment was originally proposed by J.C. Maxwell [11]. He calculated the diffusive drift of vapor from the surface of the evaporating droplet into the air, assuming that the vapor concentration at the surface of the droplet determines by saturated vapor density. This is valid for a drop radius much larger than the mean free path length of vapor molecules in the air. For example, this is not true for droplets smaller than 100 nm in the normal conditions.

Within the framework of the Maxwellian model (quasi-stationary evaporation), the diffusion equation practically turns into the Laplace equation for the vapor concentration,${\nabla }^{2}n=0$, with the following initial boundary conditions: on the drop surface $n(S)=n_s$. Outside the drop the vapor concentration is determined by the asymptotic value of the vapor concentration in the atmosphere (for aqueous solution, it is the relative humidity of the air, $n(\infty )=n_{\infty }$).

It is considered that the liquid-gas transition layer is infinitely thin compared to the droplet size. Moreover, for correctness of the Maxwellian model it is necessary to assume that density of air is much greater than vapor density, so that diffusion of vapor is determined by the vapor diffusion coefficient in the air. In particular, this requirement means that we must consider liquid at a temperature much smaller than its boiling point at a given air pressure.

The smallness of the convective Stefan flux also has to be satisfied. Stephan showed for the first time [12] that near the surface of an evaporating drop there is an air current directed away from the surface, because in order to maintain the constancy of the total pressure in the medium under conditions of vapor production by the liquid surface, along with the gradient of the vapor density there must exist an equal and opposite in direction gradient of the partial pressure of the other components of the air.

Fuchs showed [13] that the relative contribution of the Stefan flux to the evaporation process is given by factor $\frac{p_s+p_{\infty }}{2p}$, where p is the pressure of the air, $p_s$and $p_{\infty }$ are the saturated vapor pressure and the vapor pressure in the air far from the droplet, respectively. That factor for a drop of water under normal conditions does not exceed 1-2 percents (%). This is one of essential limitations of the Maxwellian model accuracy that indicates acceptable correctness of approximate solutions with respect to exact solutions of the diffusion problem. There are also other error factors: inaccurate determination of such parameters as drop surface temperature, vapor diffusion coefficient in air at a given temperature, saturated vapor pressure, etc. The capillary surface oscillations, air movement near the drop surface also are the factors that introduce uncertainties. Therefore, the approximate solutions whose accuracy is of the order of 5 percents can be considered acceptable.

In the following, we consider the small drop deposited on a flat solid horizontal substrate (so-called sessile drop). It is easy to derive that the equilibrium shape of a sessile drop of a slowly evaporating liquid, the size of which is much smaller than the capillary constant (Bond number, Bo<<1), approximately corresponds to a spherical segment with the given contact (wetting) angle. Sessile water droplets with a height of less than 1 mm satisfy this criterion quite well.

The Maxwellian model was apparently first applied to such sessile drop in the paper of Picknett and Bexon [14]. All solutions described below in this paper are given within the framework of this approach.

Exact analytical solutions have been obtained to describe the total evaporation rate (mass loss per unit time) and evaporation flux density (mass loss per unit surface area per unit time) of a small sessile liquid droplet having the shape of an axisymmetric spherical segment deposited on a horizontal substrate [15,16,17,18,19,20,21,22,23,24].

There are currently two alternative expressions for the droplet evaporation flux density that are mathematically equivalent. First, the following solution was proposed [19,20]

$$\begin{array}{l}J(\alpha )=D\frac{n_s-n_{\infty }}{R}\times \\ \times \left[\frac{\sin \theta }{2}+\sqrt{2}{(\cosh \alpha +\cos \theta )}^{3/2}{\displaystyle \underset{0}{\overset{\infty }{\int }}\frac{\cosh \theta \tau }{\cosh \pi \tau }\tanh \left[(\pi -\theta )\tau \right]P_{-1/2+i\tau }(\cosh \alpha )\tau d\tau }\right]\hspace{0.17em},\end{array}$$

(1)

where is the Legendre function of the first kind.  Here, toroidal coordinate $\alpha$ ranges in the interval from 0 (top of the drop) to ∞ (contact line). So, this coordinate is related to the cylindrical coordinate r by

$$r=\frac{R\sinh \alpha \hspace{0.17em}}{\hspace{0.17em}\cosh \alpha +\cos \theta }$$

(2)

where $\theta$ is the contact angle of the droplet (Figure 1). Equation (1) is a double integral, being an implicit function of cylindrical coordinate r, which makes formula (1) extremely difficult to use in calculations. A simpler from a computational point of view, but mathematically equivalent expression in polar coordinates was also obtained [21,22]. It allows to calculate the flux density at the surface of the droplet as a function of the polar angle φ explicitly (Figure 1):

$$J(\phi )=\frac{D(n_S-n_{\infty })}{R}f(\phi )$$

(3)

where

$$f(\phi )=\frac{\pi }{{(\pi -\theta )}^2}\frac{{\sin }^3\theta }{\cos \phi -\cos \theta }{\displaystyle \underset{\phi }{\overset{\theta }{\int }}\frac{{(1-\cos (\theta +\beta ))}^{\frac{\pi }{2(\pi -\theta )}}-{(1-\cos (\theta -\beta ))}^{\frac{\pi }{2(\pi -\theta )}}}{{\left({(1-\cos (\theta +\beta ))}^{\frac{\pi }{2(\pi -\theta )}}+{(1-\cos (\theta -\beta ))}^{\frac{\pi }{2(\pi -\theta )}}\right)}^2}\frac{{(\cos \beta -\cos \theta )}^{\frac{\pi }{2(\pi -\theta )}-\frac{1}{2}}d\beta }{\sqrt{\cos \phi -\cos \beta }}}.$$

(4)

Unfortunately, unlike expression (1), equations (3)-(4) are not yet known to the wider community and are completely ignored in the latest topical review [23].

For the total evaporation rate, the following expression was obtained [20]:

$$W(\theta )=\pi RD(n_s-n_{\infty })\left[\frac{\sin \theta }{1+\cos \theta }+4{\displaystyle \underset{0}{\overset{\infty }{\int }}\frac{1+\cosh 2\theta \tau }{\sinh 2\pi \tau }\tanh \left[(\pi -\theta )\tau \right]d\tau }\right]\hspace{0.17em}$$

(5)

where D is a diffusion coefficient of the vapor in the air, n is a vapor volume concentration outside the drop with the boundary conditions n = n_s at the drop air-liquid surface and far from the drop, R is the radius of contact line. However, as shown in [22], although the formula (5), correctly describes almost the entire dependence $W(\theta )$, where $\theta \in [0,\hspace{0.17em}\pi )$, it gives wrong result in the limit at $\theta =\pi$, which can be determined by the direct calculation.

From the point of view of application in computer simulations, this universal analytical expression that describes the evaporation flux density over the entire range of contact angles 0– π (0–180 degrees), is still quite complex. It requires significant computational resources. To accelerate the calculations of the evaporation flux density, it is reasonable to use simplified approximate expressions.

For example, there is a very good approximation for the integral evaporation flux proposed by Picknett and Bexon [14]:

$$W(\theta )=2\pi \rho D(n_s-n_{\infty })g(\theta )$$

(6)

where

$$\begin{array}{l}g(\theta )=0.6366\hspace{0.17em}\theta +0.09591\hspace{0.17em}{\theta }^2-0.06144\hspace{0.17em}{\theta }^3,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}0\le \theta \le 0.175;\\ g(\theta )=0.00008957+0.6333\theta +0.1160\hspace{0.17em}{\theta }^2-0.08878\hspace{0.17em}{\theta }^3+0.01033\hspace{0.17em}{\theta }^4,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}0.175\le \theta \le \pi .\end{array}$$

(7)

This expression has a maximum error of about 0.2% and looks much more preferable for simulation than the exact analytical solution (5). In Section 3 of this paper, we propose a much simpler approximate solution in place of equations (6)-(7) and apply it to calculate the droplet evaporation time under different contact line motion scenarios.

Similarly, for the evaporation flux density, instead of the exact formulas (1) or (3), approximate expressions can be proposed for selected narrow ranges of droplet contact angles. So, earlier an expression for the evaporation flux density was proposed, applicable in the case of small contact angles [24]:

$$J(r)=J_0(\theta ){\left(1-\frac{r^2}{R^2}\right)}^{-\lambda (\theta )}, \mathsf{\lambda }(\theta )=\frac{1}{2}-\frac{\theta }{\mathsf{\pi }},$$

(8)

where $J_0(\theta )=\frac{D(n_s-n_{\infty })}{R}(0.27{\theta }^2+1.30)(0.6381-0.2239{(\theta -0.25\pi )}^2)$.

Equation (8), being represented in the form of equations (3)-(4), can be rewritten as

$$\hspace{0.17em}J(\phi )=\frac{D(n_s-n_{\infty })}{R}\stackrel{⌣}{f}(\phi )$$

(9)

where

$$\stackrel{⌣}{f}(\phi )={\left(1-\frac{{\sin }^2\phi }{{\sin }^2\theta }\right)}^{-\lambda (\theta )}(0.27{\theta }^2+1.30)(0.6381-0.2239{(\theta -0.25\pi )}^2)$$

(10)

This expression gives a quite good description for contact angles smaller than 30 degrees (π/6). The graph that is represented in Figure 2 shows a comparison of the approximate expression (10) with the exact solution (3) for the contact angles of 1 degree and 30 degrees.

It can be seen from the comparison that the error of the approximate expression (10) increases as angle φ approaches θ. Except for the region of φ angles near θ, the deviation of the approximate formula (10) from the exact solution (4) practically does not exceed 2%, which makes formula (10) suitable for describing droplets with sharp contact angles in the range approximately 0-30 degrees.

However, the question of which simplified formulas would be appropriate to apply in other ranges of contact angles, for example, in the case of water drop deposited on hydrophobic substrates, still remains open. In Section 3 of this paper, we propose an approximate variant of exact expressions (3)- (4) in the whole range of contact angles $\theta \in [0,\hspace{0.17em}\hspace{0.17em}\pi )$) that answers this question.

2. New exact solutions for some values of contact angles

Previously [21,22], an expression was obtained that describes the vapor concentration near an evaporating drop:

$$\begin{array}{l}n(\omega ,\xi )=n_S-\frac{2(n_S-n_{\infty })}{\pi -\theta }{(\cosh \omega -\cos \xi )}^{1/2}\cos \frac{\pi \xi }{2(\pi -\theta )}\times \\ \times {\displaystyle \underset{\varsigma =\omega }{\overset{\infty }{\int }}{\left(\cosh \frac{\pi \varsigma }{(\pi -\theta )}-\cos \frac{\pi \xi }{(\pi -\theta )}\right)}^{-1}\frac{\sinh \frac{\pi \varsigma }{2(\pi -\theta )}d\varsigma }{\sqrt{\cosh \varsigma -\cosh \omega }}}\end{array}$$

(11)

This expression can also be represented as:

$$\begin{array}{l}n(\omega ,\xi )=n_s-\frac{(n_s-n_{\infty })}{2\gamma }{(\cosh \omega -\cos \xi )}^{1/2}\times \\ \times {\displaystyle \underset{\omega }{\overset{\infty }{\int }}\left[\frac{\sinh \frac{\pi \varsigma }{2\gamma }}{\cosh \frac{\pi \varsigma }{2\gamma }-\cos \frac{\pi \xi }{2\gamma }}-\frac{\sinh \frac{\pi \varsigma }{2\gamma }}{\cosh \frac{\pi \varsigma }{2\gamma }+\cos \frac{\pi \xi }{2\gamma }}\right]\frac{d\varsigma }{\sqrt{\cosh \varsigma -\cosh \omega }}}.\end{array}$$

(11a)

Here $\gamma =\pi -\theta$. Integral (11a) can be represented as the sum of a finite number of terms for some specific contact angles [25,26]. It was shown that expression (11a) under the condition

$$\xi =\frac{\pi }{2j},\mathrm{where}j=1,2,3,\dots$$

(12)

can be rewritten as

$$\begin{array}{l}n(\omega ,\xi )=n_s-(n_s-n_{\infty }){(\cosh \omega -\cos \xi )}^{1/2}\times \\ \times {\displaystyle \sum _{k=0}^{k=j-1}\left\{\frac{1}{\sqrt{\cosh \omega -\cos \left[\xi +2k\pi j^{-1}\right]}}-\frac{1}{\sqrt{\cosh \omega -\cos \left[\xi +(2k-1)\pi j^{-1}\right]}}\right\}}\end{array}$$

(13)

Then, evaporation flux density is given by [21,22]

$$J=-D\frac{\cosh \omega -\cos \xi }{R}\frac{\partial n(\omega ,\xi )}{\partial \xi }=\frac{(n_s-n_{\infty })D}{R}(\cosh \omega -\cos \xi )\frac{\partial \psi (\omega ,\xi )}{\partial \xi }$$

(14)

where

$$\begin{array}{l}\psi (\omega ,\xi )={(\cosh \omega -\cos \xi )}^{1/2}\times \\ \times {\displaystyle \sum _{k=0}^{k=j-1}\left\{\frac{1}{\sqrt{\cosh \omega -\cos \left[\xi +2k\pi j^{-1}\right]}}-\frac{1}{\sqrt{\cosh \omega -\cos \left[\xi +(2k-1)\pi j^{-1}\right]}}\right\}}\end{array}$$

By differentiating in the formula (13), we get

$$J=\frac{(n_s-n_{\infty })D}{R}\tilde{J}$$

(15)

where

$$\begin{array}{l}\tilde{J}={(\cosh \omega -\cos \xi )}^{1/2}\frac{\sin (0.5\pi j^{-1})}{2}\times \\ \times {\displaystyle \sum _{k=0}^{k=j-1}\left\{\frac{1}{\sqrt{\cosh \omega -\cos \left[\xi +2k\pi j^{-1}\right]}}-\frac{1}{\sqrt{\cosh \omega -\cos \left[\xi +(2k-1)\pi j^{-1}\right]}}\right\}}-\\ -\frac{{(\cosh \omega -\cos \xi )}^{3/2}}{2}\times \\ \times {\displaystyle \sum _{k=0}^{k=j-1}\left\{\frac{\sin \left[\xi +2k\pi j^{-1}\right]}{{\left(\cosh \omega -\cos \left[\xi +2k\pi j^{-1}\right]\right)}^{\frac{3}{2}}}-\frac{\sin \left[\xi +(2k-1)\pi j^{-1}\right]}{{\left(\cosh \omega -\cos \left[\xi +(2k-1)\pi j^{-1}\right]\right)}^{\frac{3}{2}}}\right\}}\end{array}$$

(16)

It can be shown that the first sum in formula (16) gives identically zero (can be verified by direct calculation), and the terms in brackets in the second sum are equal in absolute value. With this in mind, the expression is greatly simplified

$$\tilde{J}={(\cosh \omega -\cos \xi )}^{3/2}{\displaystyle \sum _{k=0}^{k=j-1}\frac{\sin \left[\xi +2k\pi j^{-1}\right]}{{\left(\cosh \omega -\cos \left[\xi +2k\pi j^{-1}\right]\right)}^{\frac{3}{2}}}}$$

(17)

Taking into account (11), we get

$$\tilde{J}(j)={(\cosh \omega -\cos 0.5\pi j^{-1})}^{3/2}{\displaystyle \sum _{k=0}^{k=j-1}\frac{\sin \left[0.5\pi j^{-1}+2k\pi j^{-1}\right]}{{\left(\cosh \omega -\cos \left[0.5\pi j^{-1}+2k\pi j^{-1}\right]\right)}^{\frac{3}{2}}}}$$

(18)

Thus, for any j, the evaporation flux density is given by [25]

$$J(j)=\frac{(n_s-n_{\infty })D}{R}{(\cosh \omega -\cos 0.5\pi j^{-1})}^{3/2}{\displaystyle \sum _{k=0}^{k=j-1}\frac{\sin \left[0.5\pi j^{-1}+2k\pi j^{-1}\right]}{{\left(\cosh \omega -\cos \left[0.5\pi j^{-1}+2k\pi j^{-1}\right]\right)}^{\frac{3}{2}}}}$$

(19)

If k=j, the term under summation in (19) has the form

$$\frac{\sin \left[0.5\pi j^{-1}+2j\pi j^{-1}\right]}{{\left(\cosh \omega -\cos \left[0.5\pi j^{-1}+2j\pi j^{-1}\right]\right)}^{\frac{3}{2}}}=\frac{\sin \left[0.5\pi j^{-1}\right]}{{\left(\cosh \omega -\cos \left[0.5\pi j^{-1}\right]\right)}^{\frac{3}{2}}}$$

(20)

If k=0, the term under summation in (19) has the same form

$$\frac{\sin \left[0.5\pi j^{-1}\right]}{{\left(\cosh \omega -\cos \left[0.5\pi j^{-1}\right]\right)}^{\frac{3}{2}}}=\frac{\sin \left[0.5\pi j^{-1}\right]}{{\left(\cosh \omega -\cos \left[0.5\pi j^{-1}\right]\right)}^{\frac{3}{2}}}$$

(21)

Taking into account (20) and (21), equation (19) can be transformed as

$$J(j)=\frac{(n_s-n_{\infty })D}{R}{(\cosh \omega -\cos 0.5\pi j^{-1})}^{3/2}{\displaystyle \sum _{k=1}^{k=j}\frac{\sin \left[0.5\pi j^{-1}+2k\pi j^{-1}\right]}{{\left(\cosh \omega -\cos \left[0.5\pi j^{-1}+2k\pi j^{-1}\right]\right)}^{\frac{3}{2}}}}$$

(22)

Expression (22) is the equivalent to (19).

To apply expressions (19) or (22) for calculations, it is necessary to take into account the formula (11) and following geometric relationship [21,22]:

$$\cosh \omega =\frac{{\sin }^2\xi }{\cos \phi +\cos \xi }+\cos \xi$$

(23)

To establish a relationship between the parameter j and the corresponding contact angle $\theta$, one has to use the geometric relation

$$\xi =\pi -\theta$$

(24)

It means that

$$\theta =\pi \left(1-\frac{1}{2j}\right)=\pi \frac{2j-1}{2j},\mathrm{where}j=1,2,3,\dots$$

(25)

Table 1 shows the first two solutions corresponding to j=1, 2, 3. It's obvious that

$$\underset{j\to \infty }{\lim }\theta =\pi$$

(26)

First three solution of the expression (22) with j=1,2,3 are placed into Table 1. These mean that the general solution$f(\theta )$given by equation (4) can be represented as:

$$f_{\theta ={90}^{\mathrm{o}}}(\phi )=1$$

(27)

$$f_{\theta ={135}^{\mathrm{o}}}(\phi )=\frac{1}{\sqrt{2}}-\frac{0.25}{{\left(\frac{3}{2}+\sqrt{2}\cos \phi \right)}^{\frac{3}{2}}}$$

(28)

$$f_{\theta ={150}^{\mathrm{o}}}(\phi )=\frac{1}{2}\left(1+{\left(7+4\sqrt{3}\cos \phi \right)}^{-\frac{3}{2}}-2{\left(4+2\sqrt{3}\cos \phi \right)}^{-\frac{3}{2}}\right)$$

(29)

Figure 3 represents the dependences of the evaporation flux density (dimensionless) on the polar angle, which is given by formulas (28) and (29).

It is easy to verify by direct calculation that formula (4) gives the same curves, which confirms the correctness of the mathematical transformation that led to the formula (19) or the expression (22).

3. New approximate solutions

3.1. Total evaporation rate

Let us consider compact approximate expression for the total evaporation rate of the form

$$\stackrel{⌣}{W}(\theta )=4(\ln 2)\pi \rho D(n_s-n_{\infty })\sin \frac{\theta }{2}$$

(30)

and compare formula (30) with the very accurate approximation of Picknett and Bexon (6). For this purpose, one can represent equation (30) in the form (6):

$$\stackrel{⌣}{W}(\theta )=2\pi \rho D(n_s-n_{\infty })\stackrel{⌣}{g}(\theta )$$

(31)

where

$$\stackrel{⌣}{g}(\theta )=2\ln 2\sin \frac{\theta }{2}$$

(32)

Figure 4 shows a graph of the function $\stackrel{⌣}{g}(\theta )$ together with the exact solution $g(\theta )$. We see that visually the approximate solution (30) practically does not differ from the exact one. The maximum relative error reaches 9 percent, and the average error over the entire contact angle change interval, calculated by the formula

$$\frac{1}{\pi }{\displaystyle \underset{0}{\overset{\pi }{\int }}\frac{\stackrel{⌣}{g}(x)-g(x)}{g(x)}dx\cdot 100\raisebox{1ex}{$0$}\!\left/ \!\raisebox{-1ex}{$0$}\right.}$$

(33)

does not exceed 1.5 percent (%).

Formula (30) can be recommended for estimating drop evaporation time if its contact angle varies over a wide range during evaporation. For example, this scenario occurs when a drop of aqueous solution evaporates on a hydrophobic substrate under contact line pinning conditions. Theoretically, equation (30) allows us to calculate the evaporation time of a sessile droplet at any character of dependence of the droplet radius on the contact angle during evaporation. Two extreme scenarios of contact line behavior are distinguished: the constant contact line radius mode (CCR) and the constant contact angle mode (CCA) [14].

The change in volume of a drop for time dt is defined by

$$dV=-W(\theta )n_L{}^{-1}dt$$

(34)

where $n_L$ is volume concentration of evaporating molecules in the liquid. The volume of the droplet with contact angle θ and spherical segment radius ρ (Figure 1) is determined by

$$V=\pi {\rho }^3\left(\frac{2}{3}-\cos \theta +\frac{{\cos }^3\theta }{3}\right)$$

(35)

Substituting equations (30) and (35) into expression (34), one can obtain the differential equation

$$\pi \frac{d}{dt}{\rho }^3\left(\frac{2}{3}-\cos \theta +\frac{{\cos }^3\theta }{3}\right)=-4\pi \rho D(n_s-n_{\infty })n_L{}^{-1}\ln 2\sin \frac{\theta }{2}$$

(36)

which allows to calculate the rate of drop size change under any scenario of contact line motion $\theta (t)$. In particular, the condition that is satisfied in the CCR mode can be written as follows

$$R=\rho \sin \theta =\mathrm{const}$$

(37)

In that case

$$dV=\frac{\pi R^{3}d\theta }{{(1+cos\theta )}^2}$$

By expressing $\rho$ from equation (37), substituting it into (36), and performing trigonometric transformations, we get the evaporating time of the droplet of radius R and contact angle θ for CCR mode

$$t=\frac{n_L{R}^2}{4\ln 2(n_s-n_{\infty })D}{\displaystyle \underset{0}{\overset{\raisebox{1ex}{$\theta $}\!\left/ \!\raisebox{-1ex}{$2$}\right.}{\int }}\frac{d\xi }{{\cos }^3\xi }}$$

(38)

or

$$t=\frac{n_L{R}^2}{8\ln 2(n_s-n_{\infty })D}\left[\frac{\sin \frac{\theta }{2}}{{\cos }^2\frac{\theta }{2}}+\ln (\tan (\frac{\pi +\theta }{4}))\right], \theta \in [0,\hspace{0.17em}\hspace{0.17em}\pi )$$

(39)

If we consider a drop of aqueous solution on a hydrophobic surface with an initial contact angle close to π, then its dimensions are conveniently characterized by the radius of the spherical segment and the contact angle at the initial instant of time. Then formula (39) can be conveniently rewritten as

$$t_{CCR}=\frac{n_L{\rho }^2{\sin }^2\theta }{8\ln 2(n_s-n_{\infty })D}\left[\frac{\sin \frac{\theta }{2}}{{\cos }^2\frac{\theta }{2}}+\ln (\tan (\frac{\pi +\theta }{4}))\right], \theta \in [0,\hspace{0.17em}\hspace{0.17em}\pi ) , \rho \sin \theta =\mathrm{const}.$$

(40)

Since, as follows from equation (35), the segment radius is described by the formula

$$\rho =V^{\frac{1}{3}}{\pi }^{-\frac{1}{3}}{\left(\frac{2}{3}-\cos \theta +\frac{{\cos }^3\theta }{3}\right)}^{-\frac{1}{3}}$$

then expression (40) can be rewritten as

$$t_{CCR}=\frac{n_L{V}^{\frac{2}{3}}{\sin }^2\theta }{8{\pi }^{\frac{2}{3}}\ln 2(n_s-n_{\infty })D}\left[\frac{\sin \frac{\theta }{2}}{{\cos }^2\frac{\theta }{2}}+\ln (\tan (\frac{\pi +\theta }{4}))\right]{\left(\frac{2}{3}-\cos \theta +\frac{{\cos }^3\theta }{3}\right)}^{-\frac{2}{3}}$$

(40a)

Similarly, for the evaporation time of a drop in the CCA mode one can obtain

$$t_{CCA}=\frac{3V^{\frac{2}{3}}{n}_L}{8{\pi }^{\frac{2}{3}}\ln 2D(n_s-n_0)\sin \frac{\theta }{2}}{\left(\frac{2}{3}-\cos \theta +\frac{{\cos }^3\theta }{3}\right)}^{\frac{1}{3}}, \theta =\mathrm{const}.$$

(41)

Let's look at a practical example [27]: we have the water droplet with initial volume 1 μL , and initial θ = 110 deg, vapor diffusion coefficient, D = 25.41 × 10⁻⁶ m²/s, saturated vapor concentration n_s = 0.0175 kg/m3, relative humidity (H)= 0.29. Calculations according to formulas (40a) and (40) give the following values: $t_{CCR}=860\hspace{0.17em}s$ and $t_{CCA}=974\hspace{0.17em}s$. Figure 5 shows the graphs of the evaporation time of a drop with the indicated parameters as a function of the initial volume parameter in both modes. We see that at equal initial volumes and initial contact angles, evaporation in the CCR mode is faster than in the CCA mode, as it should be, based on the general ideas about evaporation.

The authors of the paper [27] calculated the evaporation time in this case on the basis of the Popov model [20], i.e., using the equation (5). They obtained the following values: $t_{CCR}=870\hspace{0.17em}s$ and $t_{CCA}=954\hspace{0.17em}s$. We see that these values are insignificantly different from the values obtained by the formulas (40a) and (41). At the same time, as indicated in the paper [26], the experiment that was made for such a drop drying in the CCA mode gives the evaporation time of 1009 s. We can conclude that, with respect to the real experiment, the simplified formula (30) gives approximately the same error as the complex expression (5).

3.1. Evaporation flux density for acute and obtuse contact angles

The previously derived new expression for the evaporation flux density gives new possibilities for finding approximate relations. Let us consider the expression for the evaporation flux density in bipolar coordinates [21,22]:

$$J(\omega )=\frac{\pi D(n_S-n_{\infty })}{2{(\pi -\theta )}^{2}R}{(\cosh \omega -\cos (\pi -\theta ))}^{3/2}{\displaystyle \underset{\omega }{\overset{\infty }{\int }}\frac{\sinh \frac{\pi \varsigma }{2(\pi -\theta )}}{{\cosh }^2\frac{\pi \varsigma }{2(\pi -\theta )}}\frac{d\varsigma }{\sqrt{\cosh \varsigma -\cosh \omega }}}$$

(42)

where $\cosh \omega$is determined by the equation (23).

Obviously, for obtuse contact angles, the value of $(\pi -\theta )$ is small. In this case, the argument of hyperbolic functions containing $(\pi -\theta )$ in the denominator under the integral will be large, so that $\tanh \frac{\pi \varsigma }{2(\pi -\theta )}\to 1$. Consequently, the asymptotics of expression (42) is determined by ${\cosh }^{-1}\frac{\pi \omega }{2(\pi -\theta )}$. Therefore, taking into account equation (23), the asymptotic expression can be written as follows

$$J(\phi )\propto \frac{{\cosh }^{-1}\frac{\pi \omega (\phi )}{2(\pi -\theta )}}{\cos \phi -\cos \theta }$$

(43)

where

$$\omega (\phi )=\mathrm{arc}\cosh \left[\frac{1-\cos \phi \cos \theta }{\cos \phi -\cos \theta }\right]$$

(44)

Taking into account (23), one can obtain

$$\sinh \omega =\frac{\sin \phi \sin \theta }{\cos \phi -\cos \theta }$$

(45)

Performing the transformations, we have

$$\cosh \frac{\pi \omega (\phi )}{2(\pi -\theta )}=\frac{1}{2}{\left[\frac{1-\cos (\phi +\theta )}{\cos \phi -\cos \theta }\right]}^{\frac{\pi }{2(\pi -\theta )}}+\frac{1}{2}{\left[\frac{1-\cos (\phi -\theta )}{\cos \phi -\cos \theta }\right]}^{\frac{\pi }{2(\pi -\theta )}}$$

(46)

Substituting equation (46) into (43), one can obtain

$$J(\phi )\propto \frac{2{\left[\cos \phi -\cos \theta \right]}^{\frac{\pi }{2(\pi -\theta )}-1}}{{\left[1-\cos (\phi +\theta )\right]}^{\frac{\pi }{2(\pi -\theta )}}+{\left[1-\cos (\phi -\theta )\right]}^{\frac{\pi }{2(\pi -\theta )}}}$$

(47)

Comparing formula (43) with expressions (3)-(4), taking into account (47), let us write an approximate expression for function (4) in the following form

$$\stackrel{⌣}{f}(\phi )=\frac{\pi {\sin }^3\theta }{{(\pi -\theta )}^2\sin \frac{\theta }{2}}\frac{{\left[\cos \phi -\cos \theta \right]}^{\frac{\pi }{2(\pi -\theta )}-1}\times P(\theta )}{{\left[1-\cos (\phi +\theta )\right]}^{\frac{\pi }{2(\pi -\theta )}}+{\left[1-\cos (\phi -\theta )\right]}^{\frac{\pi }{2(\pi -\theta )}}}$$

(48)

where $P(\theta )$is correction factor that, in general, depends on the contact angle.

Our studies have shown that, in the first approximation, to obtain a relatively good description of the evaporation flux density in the whole range of contact angles $\theta \in [0,\hspace{0.17em}\hspace{0.17em}\pi )$using formula (48), it is convenient to choose a correction factor in the form of a constant

$$P\approx 1.166$$

(49)

Taking into account (49), expression (48) can be rewritten as

$$\stackrel{⌣}{f}(\phi )=\frac{1.166\hspace{0.17em}\hspace{0.17em}\pi {\sin }^3\theta }{{(\pi -\theta )}^2\sin \frac{\theta }{2}}\frac{{\left[\cos \phi -\cos \theta \right]}^{\frac{\pi }{2(\pi -\theta )}-1}}{{\left[1-\cos (\phi +\theta )\right]}^{\frac{\pi }{2(\pi -\theta )}}+{\left[1-\cos (\phi -\theta )\right]}^{\frac{\pi }{2(\pi -\theta )}}}$$

(50)

Evaporation flux density can be calculated by a formula similar to equation (3):

$$J(\phi )=\frac{D(n_S-n_{\infty })}{R}\stackrel{⌣}{f}(\phi )$$

(51)

Figure 6 represents the calculation results of the dimensionless evaporation flux density, determined by the approximate formula (50), together with the exact solution, determined by expression (4), for some acute and obtuse contact angles of a sessile droplet.

In general, Figure 6 shows that the approximate formula (50) can adequately describe the behavior of the evaporation flux density over the entire range of contact angles$\theta \in [0,\hspace{0.17em}\hspace{0.17em}\pi )$ , more accurately for acute angles than for obtuse angles. A more detailed approach requires breaking the entire range of contact angles into smaller sub-intervals and finding correction functions $P(\theta )$for each sub-interval, possibly in a polynomial representation. This procedure is beyond the scope of this publication.

4. Discussion

The arsenal of formulas for calculating the slow evaporation of an axisymmetric drop of capillary dimensions deposited on a flat solid surface is reviewed. Such characteristics as vapor density, evaporation flux density, total evaporation rate are considered. Exact solutions obtained in the framework of the Maxwellian model, in which the evaporation process of the drop is limited by vapor diffusion from the drop surface to the surrounding air, are presented.

Along with the long-known solutions published by Popov et al. [19,20] during the last two decades, the existence of alternative expression to describe evaporation flux density is pointed out. This alternative equation depends explicitly on the polar angle and is a one-dimensional integral (3)-(4), while the corresponding mathematical equivalent expression of Popov et al. (1) is a double integral with implicit dependence on the cylindrical coordinate. We draw the attention of researchers to the paper [22], which, apparently, remained unknown to the authors of the newest review of Wilson and D’Ambrosio [23] on drop evaporation.

New complex solution (22) was derived for the evaporation flux density of a small liquid droplet having the shape of an axisymmetric spherical segment deposited on a horizontal substrate for the set of discrete contact angles $\theta =\pi \left(1-\frac{1}{2j}\right)$, where j=1,2,3… As an example, very simple exact expressions (28) and (29) were obtained explicitly for the evaporation flux density for droplets with contact angles $\frac{3\pi }{4}=135$deg and $\frac{5\pi }{6}=150$ deg that do not contain integral function. They can also be used as approximate expressions for a narrow range of contact angles around the specified values.

Also, new approximate solutions are presented for the first time: equation (30) - total evaporation rate and expression (50) - mass loss per unit surface area per unit time in the whole range of contact angles $\theta \in [0,\hspace{0.17em}\hspace{0.17em}\pi )$). These expressions are described through elementary functions and do not contain integrals. Thus, they can be used in modeling without requiring significant computational resources.

Expression (50), taking into account (48), contains significant potential for successive improvements in accuracy through the breakdown of the contact angle determination domain into intervals and the introduction of individual correction factors. That may be a further task to advance work in this direction.