화학공학소재연구정보센터
Langmuir, Vol.36, No.1, 340-353, 2020
A Numerical Study of Micro-Droplet Spreading Behaviors on Wettability-Confined Tracks Using a Three-Dimensional Phase-Field Lattice Boltzmann Model
Wettability-confined tracks have been extensively used in open-surface microfluidic devices for their high capacity of transporting droplet pumplessly. Inspired by the experimental work of Sen et al. [Langmuir 2018, 34, 1899-1907], in the present study, a three-dimensional phase-field lattice Boltzmann model is developed and used to investigate the spreading behaviors of microdroplet on a series of wettability-confined tracks. The experimental findings are successfully reproduced through our simulation, where three distinct stages of droplet spreading on the horizontal wettability-confined diverging track are fairly exhibited, that is, the initial stage with droplet front spreading quickly, the intermediate stage with both droplet front and bulge moving forward at a constant speed, and the final stage with droplet front decelerating gradually. Moreover, a parametric study of track divergence angle is further performed, and the influential mechanism of track divergence angle on droplet spreading is further revealed. It is demonstrated that track divergence is responsible for the Laplace pressure gradient and capillary force inside the droplet, which drives the droplet bulge to move forward on the diverging track. With an increase in divergence angle, the capillary force increases linearly, which increases the droplet spreading speed at the initial and intermediate stages, while the peak capillary force comes earlier, and consequently the final decelerating stage comes earlier. On the basis of the parametric study and droplet volume conservation rule, a power law relation between track divergence angle and droplet spreading is proposed, which helps to identify the start of final decelerating stage. Finally, the droplet spreading over various inclined tracks is explored, which can be achieved only when the capillary force at the beginning is larger than the droplet gravity component along the inclined track surface.