Investigating the validity of Schrage relationships for water using molecular dynamics simulations

Effects of Organic Surface Contamination on the Mass Accommodation Coefficient of Water: A Molecular Dynamics Study

The mass accommodation coefficient (MAC), a parameter that quantifies the possibility of a phase change to occur at a liquid-vapor interface, can strongly affect the evaporation and condensation rates at a liquid surface. Due to the various challenges in experimental determination of the MAC, molecular dynamics (MD) simulations have been widely used to study the MAC on liquid surfaces with no impurities or contaminations. However, experimental studies show that airborne hydrocarbons from various sources can adsorb on liquid surfaces and alter the liquid surface properties. In this work, therefore, we study the effects of organic surface contamination, which is immiscible with water, on the MAC of water by equilibrium and nonequilibrium MD simulations. The equilibrium MD simulation results show that the MAC decreases almost linearly with increasing surface coverage of the organic contaminants. With the MAC determined from EMD simulations, the nonequilibrium MD simulation results show that the Schrage equation, which has been proven to be accurate in predicting the evaporation/condensation rates on clean liquid surfaces, is also accurate in predicting the condensation rate at contaminated water surfaces. The key assumption about the molecular velocity distribution in the Schrage analysis is still valid for condensing vapor molecules near contaminated water surfaces. We also find that under nonequilibrium conditions the adsorption of the water vapor molecules on the organic surface results in an adsorption vapor flux near the contaminated water surface. When the water surface is almost fully covered by the model organic contaminants, the adsorption flux dominates over the water condensation flux and leads to a false prediction of the MAC from the Schrage equation.

Nanobubble-Induced Aggregation of Ultrafine Particles: A Molecular Dynamics Study

Nanobubble-induced aggregation (NBIA) of fine and ultrafine particles in liquid is a promising method for enhancing floatation rates in mineral processing, cleaning contaminants from water, and reviving marine ecosystems. Although the current experimental techniques can measure the nanobubble capillary force between two surfaces with controlled approach speed, they are not capable of imaging NBIA dynamics of fine/ultrafine particles by real-time observation with nanoscale spatial resolution. In this work, we use molecular dynamics (MD) simulations to study dynamics of NBIA of Ag particles in a Lennard-Jones fluid system. The molecular-level modeling allows us to study microscopic details of NBIA dynamics that are inaccessible by current experimental means. Using MD simulations, we investigated the effects of NB size, surface wettability, surface roughness, and contact line pinning on NBIA dynamics. Our modeling results show that both concave NB bridges between two hydrophobic surfaces and convex NB bridges between two hydrophilic surfaces can result in an attractive nanobubble capillary force (NBCF) that causes the aggregation of Ag particles in liquids. The equilibrium separation between two fully aggregated particles can be well predicted by the improved capillary force model. We also observe that the change of contact angle after the contact line pinning occurs at the sharp edge of a particle, which slows the aggregation process. Our thermodynamics analysis shows that there is a critical contact angle below which the merged surface NBs will detach from the surface instead of causing aggregation. The prediction of the critical contact angle is corroborated by our MD simulation results.

Functionalized carbon nanocones performance in water harvesting

In this work, we investigate the water capture process for functionalized carbon nanocones (CNCs) through molecular dynamic simulations in the following three scenarios: a single CNC in contact with a reservoir containing liquid water, a single CNC in contact with a water vapor reservoir, and a combination of more than one CNC in contact with vapor. We found that water flows through the nanocones when in contact with the liquid reservoir if the nanocone tip presents hydrophilic functionalization. In contact with steam, we observed the formation of droplets at the base of the nanocone only when hydrophilic functionalization is present. Then, water flows through in a linear manner, a process that is more efficient than that in the liquid reservoir regime. The scalability of the process is tested by analyzing the water flow through more than one nanocone. The results suggest that the distance between the nanocones is a fundamental ingredient for the efficiency of water harvesting.

Atmospheric water harvesting using functionalized carbon nanocones

In this work, we propose a method to harvest liquid water from water vapor using carbon nanocones. The condensation occurs due to the presence of hydrophilic sites at the nanocone entrance. The functionalization, together with the high mobility of water inside nanostructures, leads to a fast water flow through the nanostructure. We show using molecular dynamics simulations that this device is able to collect water if the surface functionalization is properly selected.

Ultrahigh evaporative heat transfer measured locally in submicron water films

Thin film evaporation is a widely-used thermal management solution for micro/nano-devices with high energy densities. Local measurements of the evaporation rate at a liquid-vapor interface, however, are limited. We present a continuous profile of the evaporation heat transfer coefficient ([Formula: see text]) in the submicron thin film region of a water meniscus obtained through local measurements interpreted by a machine learned surrogate of the physical system. Frequency domain thermoreflectance (FDTR), a non-contact laser-based method with micrometer lateral resolution, is used to induce and measure the meniscus evaporation. A neural network is then trained using finite element simulations to extract the [Formula: see text] profile from the FDTR data. For a substrate superheat of 20 K, the maximum [Formula: see text] is [Formula: see text] MW/[Formula: see text]-K at a film thickness of [Formula: see text] nm. This ultrahigh [Formula: see text] value is two orders of magnitude larger than the heat transfer coefficient for single-phase forced convection or evaporation from a bulk liquid. Under the assumption of constant wall temperature, our profiles of [Formula: see text] and meniscus thickness suggest that 62% of the heat transfer comes from the region lying 0.1-1 μm from the meniscus edge, whereas just 29% comes from the next 100 μm.

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.

Probing the temperature profile across a liquid-vapor interface upon phase change

Understanding the temperature profile across a liquid-vapor interface in the presence of phase change is essential for the accurate prediction of evaporation, boiling, and condensation. It has been shown experimentally, from non-equilibrium thermodynamics and using molecular dynamics simulations, the existence of an inverted temperature profile across an evaporating liquid-vapor interface, where the vapor-side interface temperature observes the lowest value and the vapor temperature increases away from the interface, opposite to the direction of heat flow. It is worth noting, however, that an inverted temperature profile is not always the case from other experiments and simulations. In this study, we apply non-equilibrium molecular dynamics simulations to systematically study the temperature profile across a liquid-vapor interface during phase change under various heat fluxes in a two-interface setting consisting of both an evaporating and a condensing interface. The calculated vapor temperature shows different characteristics inside the Knudsen layer and in the bulk vapor. In addition, both the direction and magnitude of the vapor temperature gradient, as well as the temperature jump at the liquid-vapor interface, are functions of the applied heat flux. The interfacial entropy generation rate calculated from the vibrational density of state of the interfacial liquid and vapor molecules shows a positive production during evaporation, and the results qualitatively agree with the predictions from non-equilibrium thermodynamics.