Thermocapillary migration of a droplet on a substrate with wettability difference: A comparison between slip and precursor film models in three dimensions

Abstract

Moving contact line dynamics calculations include two models: precursor film models and slip models. The lubrication approximation method is used to establish a three-dimensional mathematical model to analyze the droplet thermocapillary migration behavior on a non-uniformly heated solid substrate with a wettability track. The contact line dynamics in the slip model and the disjoining pressure effect in the precursor model are proposed to regulate the substrate wettability. Both models are numerically implemented to investigate droplet spreading for three cases: free spreading on an isothermal substrate, thermocapillary migration on a uniform wettability substrate, and thermocapillary migration on a wettability-confined track. For the case of free spreading on an isothermal substrate, the three-dimensional results of the slip and precursor contact line models are essentially consistent with two-dimensional slip model results. For the case of thermocapillary migration on a uniform wettability substrate, the results of the two models essentially agree with the experimental results. Decreasing the thermal gradient reduces the discrepancies between the two models that result from the coordinate transformation method used in the slip model, which reduces the contact angles measured in the y-direction and enlarges the advancing contact angle in the migration direction. For the case of thermocapillary migration on a wettability-confined track, the slip model gradually shows a “dynamic-pinning” behavior with increasing equilibrium contact angle in the hydrophobic region. By contrast, the precursor film model maintains a stationary pinning behavior but separates a residual liquid outside the track. The precursor film model is preferred over the slip model in lubrication approximations for three-dimensional fluids when calculating complex moving contact dynamics caused by wettability differences. However, the precursor film model must be further optimized to prevent numerical instability.