Multicomponent Droplet Evaporation on Chemical Micro-Patterned Surfaces

Simulations of contact angle hysteresis effects in droplet sliding on inclined isothermal surfaces and droplet evaporating on heated horizontal surfaces by a lattice Boltzmann method with a variable solid-fluid interaction strength scheme

Evaporation of fluid from a pinned drop placed on solid surface proceeds via constant contact radius (CCR) mode, with a continuous reduction in the contact angle. The reduction of contact angle leads to an imbalance of interfacial tensions at the three-phase contact line. When the unbalanced force is sufficiently strong, the drop slips from the pinned contact line and slides inward. Depinning of the drop alters the mode of evaporation to constant contact angle (CCA) mode till it repins onto the surface. The change in evaporation mode from CCR to CCA is usually achieved by tuning the pinning energy barrier by controlling the surface properties of the substrate. Here, we demonstrate that the evaporation mode can be controlled by solely tailoring the surface tension of the drop, which is achieved in microgel particle-laden sessile drops that show spontaneous adsorption of microgels to the air/water interface, leading to a decrease in the interfacial tension. We show that droplets containing a sufficient number of microgels evaporate predominantly in CCR mode even on a hydrophobic surface, and the contact line remains pinned throughout the evaporation of the drop. Interestingly, the contact line dynamics can be controlled by tuning the softness of the microgels and the particle concentration in the drops.

Numerical simulation of bubble formation on orifice plates with a moving contact line

Asymmetric droplet evaporation on inclined surfaces

This work focuses on the numerical study of electrohydrodynamic multiphase fluid flow problems. Liquid bodies experience deformations caused by strong external electric fields. These deformations can be used in industrial applications to accurately control the motion of single droplets. In the limiting case, strong electric fields can force atomisation of fine microdroplets from much larger liquid bulks. Due to differences in the electrical properties of the different fluid phases, the flow affects in turn the electric field distribution. The analysis of such coupled flows requires numerical modelling of the interdependent hydrodynamic and electrodynamic problems. Since most liquids exhibit some conductivity due to intrinsic ionic species and dissolved impurities, the electrodynamic problem must be modelled by an electroquasistatic model taking into account capacitive, resistive and convective electrical currents. This electroquasistatic problem is coupled to an incompressible fluid flow problem described by the Navier-Stokes equations. Both problems are solved on the same computational grid, using the finite volume method. The fluid-fluid interface is modelled using the volume of fluid method. The resulting diffuse interface captures the motion of the phase boundaries while efficiently handling topology changes. The motion of contact lines is furthermore modelled using a dynamic contact angle model including hysteresis effects. Pinned contact lines and stick-slip contact line motion in transient problems can therefore be represented. The developed solver is readily applicable to a large range of electrohydrodynamic flow problems. This work investigates electrohydrodynamic flows occurring in three technical applications. First, the dynamics of sessile droplets subjected to an AC electric field on the surface of an insulator are considered. The dynamics of the oscillating droplets are compared with experimental data. Partial discharge inception fields are then estimated for similar configurations. Secondly, the detachment dynamics of two liquids in an on-demand droplet generator, where droplet detachment is enforced by electric fields, are considered and compared to experimental data. The different conductivities of the liquids are shown to lead to substantially different detachment dynamics. Relevant parameters in the detachment dynamics are extracted from the simulations. Finally, transient electrosprays in the cone-jet mode are characterised for a number of liquids with different electromechanical properties. The charge-radius correlations of the first ejections are found to obey the scaling laws reported in the literature. Moreover, additional scaling laws are found for the subsequent ejections.

An immersed boundary-lattice Boltzmann model for simulation of deposited particle patterns in an evaporating sessile droplet with dispersed particles

Evaporation behavior of small droplet on wafer with various film was investigated. The droplet evaporation process was categorized into two types, i.e. constant contact radius (CCR) and constant contact angle (CCA) types. In the previous research, the evaporation behavior of relatively large droplets have reported and it was confirmed that the behavior was able to predict by a simple diffusion model. In this study, the same model was applied to the small droplet. As a result, the developed model was agreed with the experimental data in the case of relatively small droplets.

Dynamic contact angles in CFD simulations

For several decades, evaporation phenomena have been intensively investigated for a broad range of applications. However, the dynamics of contact line depinning during droplet evaporation has only been inductively inferred on the basis of experimental data and remains unclear. This study focuses on the dynamics of contact line depinning during droplet evaporation based on thermodynamics. Considering the decrease in the Gibbs free energy of a system with different evaporation modes, a theoretical model was developed to estimate the receding contact angle during contact line depinning as a function of surface conditions. Comparison of experimentally measured and theoretically modeled receding contact angles indicated that the dynamics of contact line depinning during droplet evaporation was caused by the most favorable thermodynamic process encountered during constant contact radius (CCR mode) and constant contact angle (CCA mode) evaporation to rapidly reach an equilibrium state during droplet evaporation.

Exploring the impact of pillar edge effects on water drop evaporation

Liquid bridges with moving contact lines are found in a variety of settings such as capillary feeders and high-speed printing. Although it is often assumed that the length scale for these flows is small enough that inertial effects can be neglected, this is not the case in certain applications. To address this issue, we solve the Navier-Stokes equations with the finite element method for the stretching of a liquid drop between two surfaces for non-zero Reynolds numbers. We consider an axisymmetric liquid bridge between a moving flat plate and either a stationary flat plate or a cavity. The contact lines are allowed to slip, and we evaluate the effect of the Reynolds number and contact angles on the transfer of liquid to the moving plate. In the case of two flat plates, we find that inertia forces the interface to map onto a similarity solution in a manner that shifts the breakup point toward the more wettable surface. Inertia and wettability are thus competing effects, with inertia driving fluid toward the surface with the higher contact angle and wettability driving fluid toward the surface with the lower contact angle. When a cavity is present, contact line pinning on the cavity wall biases breakup toward the cavity as the Reynolds number is increased, leading to improved cavity emptying. As the flat plate is made more wettable, a second pinch-off point can form near that plate, leading to a satellite drop and a reduction in liquid transfer to the plate. Therefore, higher liquid transfer is not always obtained with a more wettable substrate when inertia is present, in contrast to Stokes flow. We also compare our results to those obtained using a model based on the long-wave approximation and find good qualitative agreement, with the long-wave model overpredicting the amount of liquid transferred after breakup.

The movement of the contact line in liquid-liquid-solid systems is a major phenomenon in natural and industrial processes. In particular, n-alkanes are widely occurring in the oil, soil pollution, and chemical industries, yet there is little knowledge on the effects of molecular chain length on the contact line movement. Here, we studied the effects of molecular chain length on the contact line movement in water/n-alkane/solid systems with different surface wettabilities. We used n-heptane (C7), n-decane (C10), and n-hexadecane (C16) as alkanes and α-quartz as the solid surface. We calculated the time-variation contact line moving velocity and also analyzed the jump frequency and the mean distance of the molecular displacement occurring within the contact line zone by molecular-kinetic theory. Molecular dynamics simulation results show that the contact line velocity decreases with increasing the chain length, originally caused by the decreasing the jump frequency and mean distance. These variations with the molecular chain length are related to the more torsions and deformations of the molecules with a longer chain length. In addition, the moving mechanism of the contact line on the same solid surface does not change at different molecular chain lengths, implying that the moving mechanism mainly depends on the three-phase wettability.

When a solid surface makes contact with a fluid-fluid interface, a three-phase contact line is formed. This contact line may move along the surface, driven by both the external fluid flow and the thermodynamically-determined contact angle that the fluid-fluid interface makes with the solid surface. Although the steady motion of contact lines has been studied extensively, the formation and unsteady motion of contact lines has not. The process of contact-line formation is of obvious interest in any physical situation where contact lines exist. In order to study the formation and motion of contact lines we have analysed a particular model problem: the interaction of an initially flat interface between two inviscid fluids with a thin, semi-infinite, flat plate. We extend the results obtained by Billingham and King (J. Fluid Mech 296 1995), where this problem was studied for the case of normal incidence, to incidence at an arbitrary initial angle.

Dynamic wetting characteristics during droplet formation in a microfluidic T-junction

Droplet evaporation is widespread in natural and industrial application, and the rapid and efficient evaporation can significantly improve energy efficiency. However, the fundamental mechanism of contact line dynamics and the microscopic characteristics of evaporating nanodroplets are not well understood. Moreover, how to design a nanostructure surface to enhance nanodroplet evaporation remains unclear. Here, through molecular dynamics simulation, we investigated the evaporation dynamics of nanodroplets on various nanoring surfaces with different geometric parameters and wettability. By measuring the changes of contact radius and contact angle, the results showed that nanodroplets successively exhibit constant contact angle (CCA), constant contact radius (CCR), and mix mode during evaporation, and the evaporation-induced CCA-CCR transition, in essence, is a Cassie-Wenzel wetting transition, whose onset time is remarkably dependent on the surface roughness and wettability. We found that this evaporation-induced wetting transition is postponed on the surface with small nanostructure spacing and weak hydrophilicity, and the evaporation rate of nanodroplets improves accordingly. The dense and hydrophobic nanostructures can not only restrain the Cassie-Wenzel transition, but also enhance the evaporation rate of nanodroplets. Last, through the potential energy field analysis of nanoring substrates, we revealed that the Cassie-Wenzel wetting transition of nanodroplets is a process of molecule migration to low potential energy regions. Our work provides guidance for designing nanostructure surfaces to effectively control the droplet wetting state and enhance its mass transfer performance of phase change.