Numerical Analysis of the Dispersion and Deposition of Particles in Evaporating Sessile Droplets.

Evaporating sessile droplets containing dispersed particles are used in different technological applications, such as 3D printing, biomedicine, and micromanufacturing, where an accurate prediction of both the dispersion and deposition of the particles is important. Furthermore, the interaction between the droplet and the substrate must be taken into account: the motion of the contact line, in particular, must be modeled carefully. To this end, studies have typically been limited to either pinned or moving contact lines to simplify the underlying mathematical models and numerical methods, neglecting the fact that both scenarios are observed during the evaporation process. Here, a numerical algorithm considering both contact line regimes is proposed whereby the regimes are distinguished by predefined threshold contact angles. After a detailed validation, this new algorithm is applied to study the influence of both regimes on the dispersion and deposition of particles in an evaporating sessile droplet. In particular, the presented analysis focuses on the influence of (i) the contact line motion characteristics by varying the limiting contact angle and spreading speed, (ii) the Marangoni number, characterizing the importance of thermocapillarity, (iii) the evaporation number, which quantifies the importance of evaporation, (iv) the Damköhler number, a measure of the particle deposition rate, and (v) the Peclet number, which compares the convection and diffusion of the particle concentration. When thermocapillarity becomes dominant or the limiting contact angle is larger, the particle accumulation near the contact line decreases, which, in turn, means that more particles are deposited near the center of the droplet. In contrast, increasing the evaporation number supports particle accumulation near the contact line, while a larger Damköhler number and/or smaller Peclet number yield more uniform final deposition patterns. Finally, a larger characteristic speed of spreading results in fewer particles being deposited at the center of the droplet.