Atomic Scale Interfacial Transport at an Extended Evaporating Meniscus

Nanoscale Meniscus Dynamics in Evaporating Thin Films: Insights from Molecular Dynamics Simulations

Evaporation studies are focused on unraveling heat transfer and flow dynamics near the solid-liquid-vapor contact line, particularly focusing on the meniscus, which encompasses the nonevaporating adsorbed layer, thin-film, and bulk meniscus regions. Continuum models assume that there are no evaporating adsorbed layers due to the strong intermolecular forces. However, recent molecular dynamics (MD) simulations have unveiled the significant role of adsorbed layers in thin-film evaporation. Leveraging a recently published energy-based interface detection method, the current study presents nonequilibrium MD simulation results for thin-film evaporation from a phase-change-driven nanopump using liquid argon confined between parallel platinum plates. Notably, unlike the transient simulations often encountered in the literature, the simulation system achieves a statistically steady transport. In this context, we showcase the shapes of the evaporating menisci for two distinct channel heights, 8 and 16 nm, and elucidate the underlying flow physics through velocity vectors and temperature contours. This comprehensive investigation advances our understanding of thin-film evaporation and its mechanisms, offering insights that span from nanoscale phenomena to broader thermal management applications.

The Critical Diameter for Continuous Evaporation Is between 3 and 4 nm for Hydrophilic Nanopores

Evaporation studies of water using classical molecular dynamics simulations are largely limited due to their high computational expense. This study addresses that issue by developing coarse-grained molecular dynamics models based on Morse potential. Models are optimized based on multi-temperature and at room temperature using machine learning techniques like Genetic Algorithm, Nelder-Mead algorithm, and Strength Pareto Evolutionary Algorithm. The multi-temperature-based model named as Morse-D is found to be more accurate than the single temperature model in representing the water properties at higher temperatures. Using this Morse-D water model, evaporation from hydrophilic nanopores with pore diameter varying from 2 to 5 nm is studied. Our results show that the critical diameter to initiate continuous evaporation at nanopores lies between 3 and 4 nm. A maximum heat flux of 21.3 kW/cm² is observed for a pore diameter of 4.5 nm and a maximum mass flow rate of 16.2 ng/s for a pore diameter of 5 nm. The observed heat flux is an order of magnitude times larger than the currently reported values from experiments in the literature for water, which indicates that we need to focus on nanoscale evaporation to enhance the critical heat flux.

Energy-Based Interface Detection for Phase Change Processes of Monatomic Fluids in Nanoconfinements

An energy-based liquid-vapor interface detection method is presented using molecular dynamics simulations of liquid menisci confined between two parallel plates under equilibrium and evaporation/condensation conditions. This method defines the liquid-vapor interface at the location where the average kinetic energy of atoms first exceeds the average potential energy imposed by all neighboring molecules. This definition naturally adapts to the location of the menisci relative to the walls and can properly model the behavior of the liquid adsorbed layers. Unlike the density cutoff methods frequently used in the literature that suffer from density layering effects, this new method gives smooth and continuous liquid-vapor interfaces in nanoconfinements. Surface tension values calculated from the equilibrium MD simulations match the Young-Laplace equation better when using the radius of curvatures calculated from this method. Overall, this energy-based liquid-vapor interface detection method can be used in studies of nanoscale phase change processes and other relevant applications.

Enhanced Evaporation of Ultrathin Water Films on Silicon-Terminated Si3N4 Nanopore Membranes

Water evaporation confined in nanoscale is a ubiquitous phenomenon in nature and has crucial importance in a broad range of technical applications. With the nonequilibrium molecular dynamics simulation, we elucidate nanothin water film evaporation characteristics on a silicon nitride nanopore membrane considering the effects of pore size and pore chemistry. Pore chemistry plays the main role in regulating the evaporation flux. The terminated Si atoms on the pore surface lead to a higher evaporation intensity than the N ones. We attribute this enhancement to the transition of the structural properties of fluid, where liquid molecules are packed loosely and disorderedly under the inducement of terminated silicon atoms. The findings in the present work can contribute to the fundamental understanding of the nanopore-enhanced evaporation process and provide new guidance to the design of advanced nanopore membrane materials.

Evaporation in nano/molecular materials

Evaporation is a physical phenomenon with fundamental significance to both nature and technology ranging from plant transpiration to DNA engineering. Various analytical and empirical relationships have been proposed to characterize evaporation kinetics at macroscopic scales. However, theoretical models to describe the kinetics of evaporation from nano and sub-nanometer (molecular) confinements are absent. On the other hand, the fast advancements in technology concentrated on development of nano/molecular-scale devices demand appropriate models that can accurately predict physics of phase-change in these systems. A thorough understanding of the physics of evaporation in nano/molecular materials is, thus, of critical importance to develop the required models. This understanding is also crucial to explain the intriguing evaporation-related phenomena that only take place when the characteristic length of the system drops to several nanometers. Here, we comprehensively review the underlying physics of evaporation phenomenon and discuss the effects of nano/molecular confinement on evaporation. The role of liquid-wall interface-related phenomena including the effects of disjoining pressure and flow slippage on evaporation from nano/molecular confinements are discussed. Different driving forces that can induce evaporation in small confinements, such as heat transfer, pressure drop, cavitation and density fluctuations are elaborated. Hydrophobic confinement induced evaporation and its potential application for synthetic ion channels are discussed in detail. Evaporation of water as molecular clusters rather than isolated molecules is discussed. Despite the lack of experimental investigations on evaporation at nanoscale, there exist an extensive body of literature that have applied different simulation techniques to predict the phase change behavior of liquids in nanoconfinements. We infer that exploring the effect of electrostatic interactions and flow slippage to enhance evaporation from nanoconduits is an interesting topic for future endeavors. Further future studies could be devoted to developing nano/molecular channels with evaporation-based gating mechanism and utilization of 2D materials to tune energy barrier for evaporation leading to enhanced evaporation.

Stable and Efficient Nanofilm Pure Evaporation on Nanopillar Surfaces

Molecular dynamics simulations were conducted to systematically investigate how to maintain and enhance nanofilm pure evaporation on nanopillar surfaces. First, the dynamics of the evaporation meniscus and the onset and evolution of nanobubbles on nanopillar surfaces were characterized. The meniscus can be pinned at the top surface of the nanopillars during evaporation for perfectly wetting fluid. The curvature of the meniscus close to nanopillars varies dramatically. Nanobubbles do not originate from the solid surface, where there is an ultrathin nonevaporation film due to strong solid-fluid interaction, but originate and evolve from the corner of nanopillars, where there is a quick increase in potential energy of the fluid. Second, according to a parametric study, the smaller pitch between nanopillars (P) and larger diameter of nanopillars (D) are found to enhance evaporation but also raise the possibility of boiling, whereas the smaller height of nanopillars (H) is found to enhance evaporation and suppress boiling. Finally, it is revealed that the nanofilm thickness should be maintained beyond a threshold, which is 20 Å in this work, to avoid the suppression effect of disjoining pressure on evaporation. Moreover, it is revealed that whether the evaporative heat transfer is enhanced on the nanopillar surface compared with the smooth surface is also affected by the nanofilm thickness. The value of nanofilm thickness should be determined by the competition between the suppression effect on evaporation due to the decrease in the volume of supplied fluid and the existence of capillary pressure and the enhancement effect on evaporation due to the increase in the heating area. Our work serves as the guidelines to achieve stable and efficient nanofilm pure evaporative heat transfer on nanopillar surfaces.

Temperature Effect on Capillary Flow Dynamics in 1D Array of Open Nanotextured Microchannels Produced by Femtosecond Laser on Silicon

Capillary flow of water in an array of open nanotextured microgrooves fabricated by femtosecond laser processing of silicon is studied as a function of temperature using high-speed video recording. In a temperature range of 23-80 °C, the produced wicking material provides extremely fast liquid flow with a maximum velocity of 37 cm/s in the initial spreading stage prior to visco-inertial regime. The capillary performance of the material enhances with increasing temperature in the inertial, visco-inertial, and partially in Washburn flow regimes. The classic universal Washburn's regime is observed at all studied temperatures, giving the evidence of its universality at high temperatures as well. The obtained results are of great significance for creating capillary materials for applications in cooling of electronics, energy harvesting, enhancing the critical heat flux of industrial boilers, and Maisotsenko cycle technologies.

Joint effect of advection, diffusion, and capillary attraction on the spatial structure of particle depositions from evaporating droplets

A simplified model is developed, which allows us to perform computer simulations of the particles transport in an evaporating droplet with a contact line pinned to a hydrophilic substrate. The model accounts for advection in the droplet, diffusion, and particle attraction by capillary forces. On the basis of the simulations, we analyze the physical mechanisms of forming of individual chains of particles inside the annular sediment. The parameters chosen correspond to the experiments of Park and Moon [Langmuir 22, 3506 (2006)LANGD50743-746310.1021/la053450j], where an annular deposition and snakelike chains of colloid particles have been identified. The annular sediment is formed by advection and diffusion transport. We find that the close packing of the particles in the sediment is possible if the evaporation time exceeds the characteristic time of diffusion-based ordering. We show that the chains are formed by the end of the evaporation process due to capillary attraction of particles in the region bounded by a fixing radius, where the local droplet height is comparable to the particle size. At the beginning of the evaporation, the annular deposition is shown to expand faster than the fixing radius moves. However, by the end of the process, the fixing radius rapidly outreaches the expanding inner front of the ring. The snakelike chains are formed at this final stage when the fixing radius moves toward the symmetry axis.