Spreading Time of Impacting Nanodroplets

What Controls the Hole Formation of Nanodroplets: Hydrodynamic or Thermodynamic Instability?

Using molecular dynamics simulations, we investigate the air hole formation of water nanodroplets impacting hydrophilic to hydrophobic surfaces in the range of static contact angles from 30° to 140° with different initial surface temperatures ranging from 300 to 1000 K. We show that the hole dynamics of nanodroplets are different from those observed in millimeter-sized droplets. The hole formation can be observed on smooth surfaces for nanodroplets; however, it only occurs on nonsmooth surfaces for millimeter-sized droplets. We clarify that the hole formation of nanodroplets is triggered by a nucleated vapor bubble due to thermodynamic instability, whereas it is initiated by air bubble entrapment during impact due to hydrodynamic instability for millimeter-sized droplets. The hole formation of nanodroplets relies heavily on the surface temperature and surface wettability, because the nucleated vapor bubble more easily occurs and grows on the surface with high initial temperatures and hydrophobic surfaces. Based on the thermal stability analysis, a criterion is developed to predict the hole formation of nanodroplets, which verifies the dependence of hole formation on the surface temperature and wettability. Furthermore, we show that the ring-bouncing of nanodroplets is triggered by the nucleated vapor bubble. We clarify the reasons for the reduced contact time of nanodroplets caused by the ring-bouncing.

Impact of nanodroplets on cone-textured surfaces

Molecular dynamics simulations have been performed to study the dynamics of nanodroplets impacting on a flat superhydrophobic surface and surfaces covered with nanocone structures. We present a panorama of nanodroplet behaviors for a wide range of impact velocities and different cone geometrics, and develop a model to predict whether a nanodroplet impacting onto cone-textured surfaces will touch the underlying substrate during impact. The advantages and disadvantages of applying nanocone structures to the solid surface are revealed by the investigations into restitution coefficient and contact time. The effects of nanocone structures on droplet bouncing dynamics are probed using momentum analysis rather than conventional energy analysis. We further demonstrate that a single Weber number is inadequate for unifying the dynamics of macroscale and nanoscale droplets on cone-textured surfaces, and propose a combined dimensionless number to address it. The extensive findings of this study carry noteworthy implications for engineering applications, such as nanoprinting and nanomedicine on functional patterned surfaces, providing fundamental support for these technologies.

Oblique Impacts of Nanodroplets upon Surfaces

In this work, oblique impacts of nanodroplets impacting surfaces in a wide range of impact angles (α) are investigated in detail via molecular dynamics simulations. Five outcomes are observed, including deposition, prompt splashing, break-up, separation, and shattering. With increasing impact angle, the outcomes of prompt splashing, break-up, separation, and shattering are enlarged but the one of deposition is compressed. By drawing a We_n ∼ α phase diagram, the outcome regimes and corresponding boundaries of them can be successfully identified, and the boundary between the deposition and other outcome regimes is theoretically modeled and shows good agreement with the phase diagram, where We_n is the normal impact Weber number. For further understanding of the oblique impacts, the maximum spreading factor, as the feature parameter of spreading, is investigated. Asymmetry spreading behaviors are observed, noting that β_max,∥ is always larger than β_max,⊥. β_max,⊥ is tested that it only depends on We_n with wide impact angles and could be predicted by the scaling law of β_max,⊥ = 0.7We_n^1/4. However, β_max,∥ depends on not only We_n but also impact angles. A modified model is proposed for predicting β_max,∥ as 0.7We_n^1/4 + 0.001(We_n tan² α)^3/2, which shows good agreement with data on surfaces with θ from 73 to 105° in wide We_n and α ranges.