Analysis of the microfluid flow in an evaporating sessile droplet

The axisymmetric time-dependent flow field in an evaporating sessile droplet whose contact line is pinned is studied numerically and using an analytical lubrication theory with a zero-shear-stress boundary condition on the free surface of the droplet at low capillary and Reynolds numbers. A finite element algorithm is developed to solve simultaneously the vapor concentration and flow field in the droplet under conditions of slow evaporation. The finite element solution confirms the accuracy of the lubrication solution, especially when terms of higher order in the droplet flatness ratio (the ratio of droplet height to radius, h/R) are included in the lubrication theory to account more accurately for the singular flow near the contact line.