Thermal singularity and droplet motion in one-component fluids on solid substrates with thermal gradients

Thermal singularity and contact line motion in pool boiling: Effects of substrate wettability

The dynamic van der Waals theory [Phys. Rev. E 75, 036304 (2007)] is employed to model the growth of a single vapor bubble in a superheated liquid on a flat homogeneous substrate. The bubble spreading dynamics in the pool boiling regime has been numerically investigated for one-component van der Waals fluids close to the critical point, with a focus on the effect of the substrate wettability on bubble growth and contact line motion. The substrate wettability is found to control the apparent contact angle and the rate of bubble growth (the rate of total evaporation), through which the contact line speed is determined. An approximate expression is derived for the contact line speed, showing good agreement with the simulation results. This demonstrates that the contact line speed is primarily governed by (1) the circular shape of interface (for slow bubble growth), (2) the constant apparent contact angle, and (3) the constant bubble growth rate. It follows that the contact line speed has a sensitive dependence on the substrate wettability via the apparent contact angle which also determines the bubble growth rate. Compared to hydrophilic surfaces, hydrophobic surfaces give rise to a thinner shape of bubble and a higher rate of total evaporation, which combine to result in a much faster contact line speed. This can be linked to the earlier formation of a vapor film and hence the onset of boiling crisis.

Single-bubble dynamics in pool boiling of one-component fluids

We numerically investigate the pool boiling of one-component fluids with a focus on the effects of surface wettability on the single-bubble dynamics. We employed the dynamic van der Waals theory [Phys. Rev. E 75, 036304 (2007)], a diffuse-interface model for liquid-vapor flows involving liquid-vapor transition in nonuniform temperature fields. We first perform simulations for bubbles on homogeneous surfaces. We find that an increase in either the contact angle or the surface superheating can enhance the bubble spreading over the heating surface and increase the bubble departure diameter as well and therefore facilitate the transition into film boiling. We then examine the dynamics of bubbles on patterned surfaces, which incorporate the advantages of both hydrophobic and hydrophilic surfaces. The central hydrophobic region increases the thermodynamic probability of bubble nucleation while the surrounding hydrophilic region hinders the continuous bubble spreading by pinning the contact line at the hydrophobic-hydrophilic intersection. This leads to a small bubble departure diameter and therefore prevents the transition from nucleate boiling into film boiling. With the bubble nucleation probability increased and the bubble departure facilitated, the efficiency of heat transfer on such patterned surfaces is highly enhanced, as observed experimentally [Int. J. Heat Mass Transfer 57, 733 (2013)]. In addition, the stick-slip motion of contact line on patterned surfaces is demonstrated in one-component fluids, with the effect weakened by surface superheating.

Molecular dynamics simulations for the motion of evaporative droplets driven by thermal gradients along nanochannels

For a one-component fluid on a solid substrate, a thermal singularity may occur at the contact line where the liquid-vapor interface intersects the solid surface. Physically, the liquid-vapor interface is almost isothermal at the liquid-vapor coexistence temperature in one-component fluids while the solid surface is almost isothermal for solids of high thermal conductivity. Therefore, a temperature discontinuity is formed if the two isothermal interfaces are of different temperatures and intersect at the contact line. This leads to the so-called thermal singularity. The localized hydrodynamics involving evaporation/condensation near the contact line leads to a contact angle depending on the underlying substrate temperature. This dependence has been shown to lead to the motion of liquid droplets on solid substrates with thermal gradients (Xu and Qian 2012 Phys. Rev. E 85 061603). In the present work, we carry out molecular dynamics (MD) simulations as numerical experiments to further confirm the predictions made from our previous continuum hydrodynamic modeling and simulations, which are actually semi-quantitatively accurate down to the small length scales in the problem. Using MD simulations, we investigate the motion of evaporative droplets in one-component Lennard-Jones fluids confined in nanochannels with thermal gradients. The droplet is found to migrate in the direction of decreasing temperature of solid walls, with a migration velocity linearly proportional to the temperature gradient. This agrees with the prediction of our continuum model. We then measure the effect of droplet size on the droplet motion. It is found that the droplet mobility is inversely proportional to a dimensionless coefficient associated with the total rate of dissipation due to droplet movement. Our results show that this coefficient is of order unity and increases with the droplet size for the small droplets (~10 nm) simulated in the present work. These findings are in semi-quantitative agreement with the predictions of our continuum model. Finally, we measure the effect of liquid-vapor coexistence temperature on the droplet motion. Through a theoretical analysis on the size of the thermal singularity, it can be shown that the droplet mobility decreases with decreasing coexistence temperature. This is observed in our MD simulations.

Hydrodynamics of Leidenfrost droplets in one-component fluids

Using the dynamic van der Waals theory [Phys. Rev. E 75, 036304 (2007)], we numerically investigate the hydrodynamics of Leidenfrost droplets under gravity in two dimensions. Some recent theoretical predictions and experimental observations are confirmed in our simulations. A Leidenfrost droplet larger than a critical size is shown to be unstable and break up into smaller droplets due to the Rayleigh-Taylor instability of the bottom surface of the droplet. Our simulations demonstrate that an evaporating Leidenfrost droplet changes continuously from a puddle to a circular droplet, with the droplet shape controlled by its size in comparison with a few characteristic length scales. The geometry of the vapor layer under the droplet is found to mainly depend on the droplet size and is nearly independent of the substrate temperature, as reported in a recent experimental study [Phys. Rev. Lett. 109, 074301 (2012)]. Finally, our simulations demonstrate that a Leidenfrost droplet smaller than a characteristic size takes off from the hot substrate because the levitating force due to evaporation can no longer be balanced by the weight of the droplet, as observed in a recent experimental study [Phys. Rev. Lett. 109, 034501 (2012)].

Enhanced droplet control by transition boiling

A droplet of water on a heated surface can levitate over a film of gas produced by its own evaporation in the Leidenfrost effect. When the surface is prepared with ratchet-like saw-teeth topography, these droplets can self-propel and can even climb uphill. However, the extent to which the droplets can be controlled is limited by the physics of the Leidenfrost effect. Here, we show that transition boiling can be induced even at very high surface temperatures and provide additional control over the droplets. Ratchets with acute protrusions enable droplets to climb steeper inclines while ratchets with sub-structures enable their direction of motion to be controlled by varying the temperature of the surface. The droplets' departure from the Leidenfrost regime is assessed by analysing the sound produced by their boiling. We anticipate these techniques will enable the development of more sophisticated methods for controlling small droplets and heat transfer.