Expressions for the Evaporation and Condensation Coefficients in the Hertz-Knudsen Relation

Accelerated Zymonic Acid Formation from Pyruvic Acid at the Interface of Aqueous Nanodroplets

To explore the role of the liquid interface in mediating reactivity in small compartments, the formation kinetics of zymonic acid (ZA) is measured in submicron aerosols (average radius = 240 nm) using mass spectrometry. The formation of ZA, from a condensation reaction of two pyruvic acid (PA) molecules, proceeds over days in bulk solutions, while in submicron aerosols, it occurs in minutes. The experimental results are replicated in a kinetic model using an apparent interfacial reaction rate coefficient of krxn = (0.9 ± 0.2) × 10-3 M-1 s-1. The simulation reveals that surface activity of PA coupled with an enhanced interfacial reaction rate drives accelerated ZA formation in aerosols. Experimental and simulated results provide compelling evidence that the condensation reaction of PA occurs exclusively at the aerosol interface with a reaction rate coefficient that is enhanced by 4 orders of magnitude (∼104) relative to what is estimated for macroscale solutions.

The Problem With National Institute of Standards and Technology Thermodynamics Tables in Continuum Mechanics

Thermodynamics is a fundamental topic of continuum mechanics and biomechanics, with a wide range of applications to physiological and biological processes. This study addresses two fundamental limitations of current thermodynamic treatments. First, thermodynamics tables distributed online by the U.S. National Institute of Standards and Technology (NIST) report properties of fluids as a function of absolute temperature T and absolute pressure P. These properties include mass density ρ, specific internal energy u, enthalpy h=u+P/ρ, and entropy s. However, formulations of jump conditions across phase boundaries derived from Newton's second law of motion and the first law of thermodynamics employ the gauge pressure p=P-Pr, where Pr is an arbitrarily selected referential absolute pressure. Interchanging p with P is not innocuous as it alters tabulated NIST values for u while keeping h and s unchanged. Using p for functions of state and governing equations solves the problem with using NIST entries for the specific internal energy u in standard thermodynamics tables and analyses of phase transformation in continuum mechanics. Second, constitutive models for the free energy of fluids, such as water and air, are not typically provided in standard thermodynamics treatments. This study proposes a set of constitutive models and validates them against suitably modified NIST data.

Chemical composition from photos: Dried solution drops reveal a morphogenetic tree

Under nonequilibrium conditions, inorganic systems can produce a wealth of life-like shapes and patterns which, compared to well-formed crystalline materials, remain widely unexplored. A seemingly simple example is the formation of salt deposits during the evaporation of sessile droplets. These evaporites show great variations in their specific patterns including single rings, creep, small crystals, fractals, and featureless disks. We have explored the patterns of 42 different salts at otherwise constant conditions. Based on 7,500 images, we show that distinct pattern families can be identified and that some salts (e.g., Na2SO4 and NH4NO3) are bifurcated creating two distinct motifs. Family affiliations cannot be predicted a priori from composition alone but rather emerge from the complex interplay of evaporation, crystallization, thermodynamics, capillarity, and fluid flow. Nonetheless, chemical composition can be predicted from the deposit pattern with surprisingly high accuracy even if the set of reference images is small. These findings suggest possible applications including smartphone-based analyses and lightweight tools for space missions.

Unified theoretical framework for temperature regulation via phase transition

Phase transition usually consumes or releases energy to produce cooling or heating within different materials, providing a generalized framework for temperature regulation in practical applications. Because of the strong coupling between the enthalpy change in thermodynamics and heat-mass transfer kinetics, unveiling the mechanism of temperature regulation via the phase transition remains a great challenge. Here, we develop a new theoretical method by establishing a connection of enthalpy change from thermodynamics to phase transition dynamics to study evaporation-induced cooling as an example. Our new approach can spontaneously generate evaporative cooling at interfaces, and the predicted results are consistent with recent experiments. The evaporation-induced steady vapor is dictated by an anomalous cold-to-hot mass transfer through temperature-dependent chemical potentials, which enables temperature regulation inside liquids via a thermodynamic-kinetic interplay. Moreover, we show that a simple prohibition of heat exchange between liquids and reservoir can greatly enhance the cooling magnitude by a factor of 2∼4, which is highly dependent on the thermodynamics and kinetic coefficients of liquids. Our new method paves the way for exploration of cooling or heating induced by different phase transitions, such as evaporation, sublimation, or condensation, in a unified framework, which can significantly promote the development of temperature regulation by phase transitions.

Unexpected concentration dependence of the mass accommodation coefficient of water on aqueous triethylene glycol droplets

The mass accommodation coefficient αM of water on aqueous triethylene glycol droplets was determined for water mole fractions in the range xmol = 0.1-0.93 and temperatures between 21 and 26 °C from modulated Mie scattering measurement on single optically-trapped droplets in combination with a kinetic multilayer model. αM reaches minimum values around 0.005 at a critical water concentration of xmol = 0.38, and increases with decreasing water content to a value of ≈0.1 for almost pure triethylene glycol droplets, essentially independent of the temperature. Above xmol = 0.38, αM first increases with increasing water content and then stabilises at a value of ≈0.1 at the lowest temperatures, while at the highest temperature its value remains around 0.005. We analysed the unexpected concentration and temperature dependence with a previously proposed two-step model for mass accommodation which provides concentration and temperature-dependent activation enthalpies and entropies. We suggest that the unexpected minimum in αM at intermediate water concentrations might arise from a more or less saturated hydrogen-bond network that forms at the droplet surface.

A new setup for measurements of absolute saturation vapor pressures using a dynamical method: Experimental concept and validation

We describe a new experimental system for direct measurements of the absolute saturation vapor pressures of liquid or solid samples. The setup allows the isolation of the sample under steady conditions in an ultra-high vacuum chamber, where the measurement of the sample's vapor pressure as a function of its temperature can be performed in a range around room temperature and in a pressure range defined only by the applied absolute pressure sensor. We characterize the setup and illustrate its capability to measure saturation vapor pressures as well as enthalpies of evaporation around room temperature with explicit measurements on four liquid compounds (diethyl phthalate, 1-decanol, 1-heptanol, and 1-hexanol) for which accurate vapor pressures have previously been reported.

Solar-driven abnormal evaporation of nanoconfined water

Intrinsic water evaporation demands a high energy input, which limits the efficacy of conventional interfacial solar evaporators. Here, we propose a nanoconfinement strategy altering inherent properties of water for solar-driven water evaporation using a highly uniform composite of vertically aligned Janus carbon nanotubes (CNTs). The water evaporation from the CNT shows the unexpected diameter-dependent evaporation rate, increasing abnormally with decreasing nanochannel diameter. The evaporation rate of CNT10@AAO evaporator thermodynamically exceeds the theoretical limit (1.47 kg m-2 hour-1 under one sun). A hybrid experimental, theoretical, and molecular simulation approach provided fundamental evidence of different nanoconfined water properties. The decreased number of H-bonds and lower interaction energy barrier of water molecules within CNT and formed water clusters may be one of the reasons for the less evaporative energy activating rapid nanoconfined water vaporization.

Evaluation of Evaporation Fluxes for Pesticides and Low Volatile Hazardous Materials Based on Evaporation Kinetics of Net Liquids

Evaporation is the phase transition process that plays a significant role in many spheres of life and science. Volatilization of hazardous materials, pesticides, petroleum spills, etc., impacts the environment and biosphere. Predicting evaporation fluxes under specific environmental conditions is challenging from theoretical and empirical points of view. A new practical method for estimating fluxes is proposed based on our experimental results and previously published data. It is demonstrated that some parameters in theoretical equations for near-equilibrium evaporation can be estimated from experiments, and these formulas can be exploited to predict steady-state evaporation fluxes in the air in a range of 8 orders of magnitude based on a single experiment carried out for nontoxic volatile compounds.

Highly efficient molecular film for inhibiting volatilization of hazardous nitric acid

Nitric acid, an important basic chemical raw material, plays an important role in promoting the development of national economy. However, such liquid hazardous chemicals are easy to cause accidental leakage during production, transportation, storage and use. The high concentration and corrosive toxic gas generated from decomposition shows tremendous harm to the surrounding environment and human life safety. Therefore, how to inhibit the volatilization of nitric acid and effectively control and block the generation of the toxic gas in the first time are the key to deal with the nitric acid leakage accident. Herein, a new method of molecular film obstruction is proposed to inhibit the nitric acid volatilization. The molecular film inhibitor spontaneously spread and form an insoluble molecular film on the gas-liquid interface, changing the state of nitric acid liquid surface and inhibiting the volatilization on the molecular scale. The inhibition rate up to 96% can be achieved below 45 °C within 400 min. Cluster structure simulation and energy barrier calculation is performed to elucidate the inhibition mechanism. Theoretical analysis of energy barrier shows that the specific resistance of the inhibitor significantly increased to 460 s·cm-1 at 45 °C, and the generated energy barrier is about 17,000 kJ·mol-1, which is much higher than the maximum energy required for nitric acid volatilization of 107.97 kJ·mol-1. The molecular film obstruction strategy can effectively inhibit the volatilization of nitric acid. This strategy paves the way for preventing the volatilization of liquid hazardous chemicals in accidental leakage treatment.

Handling Nanoparticle Content in Nanocomposite Thin Films Deposited by Misty Plasma Processes through Controlled Flash Boiling Atomization

Misty plasma processes based on colloidal solutions sprayed into low-pressure plasmas have recently shown great potential for multifunctional thin film deposition. In such processes, nanoparticle accumulation in ring-shaped structures remains the main obstacle to the synthesis of high-quality coatings containing abundant, small-scale, and evenly dispersed nanoparticles. These local buildups appear after a colloidal droplet evaporates from a substrate. Accordingly, controlling the droplets' size in the spray is of key importance to ensure a uniform nanoparticle content in the plasma-deposited nanocomposite film. In this work, it is shown that the use of more volatile solvents produces finer droplets on the substrate, thereby improving nanoparticle dispersion in the matrix. A one-dimensional evaporation model is further developed and used to show that, contrary to what one might expect, this result cannot be attributed to faster evaporation during droplet transport in the low-pressure plasma. Instead, a so-called "flash" boiling atomization mechanism is discussed to support the experimental findings.

Structural segregation in GSR from mercuric primers. A preliminary study

Since the second half of the XX century, primer mixtures based on mercury fulminate have become a rare occurrence on small ammunition markets in Western Europe and North America. As a consequence, Hg-containing gunshot residue (GSR) particles have not been as deeply investigated as residues from lead-based primer mixtures. As a matter of fact, no mention of GSR particles from mercuric primers is made by the current ASTM standard procedure for gunshot residue analysis. However, those laboratories dealing with ammunition and firearms produced in Eastern Europe or Asia still have a forensic interest in Hg-containing GSR. In this paper, a brief description of chemical composition and inner morphology of GSR particles from three different mercuric primers is reported. Regarding composition, arguments are given to promote SbSnHg residues to Characteristic of GSR particles when mercuric primers are discharged. From a morphological point of view, presence of inner nodules and other inhomogeneities were shown in GSR particles milled in a FIB/SEM. Moreover, mercury vaporization under the electron beam was observed for a particle reduced to a lamella. Mercury evanescence in GSR was interpreted in terms of mercury segregation during particle formation and higher mobility of Hg atoms in presence of defects (vacancies) in a strained lattice.

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.

Nonequilibrium liquid-vapor interfaces: Linear and nonlinear descriptions

While it is often assumed that liquid-vapor interfaces in nonequilibrium processes are in states of local thermodynamic equilibrium, this might not be the case for strong deviations from equilibrium. Clausius-Clapeyron equations for bulk properties yield a consistently defined temperature of the interface that is close to the liquid bulk temperature. The alternative interface temperature defined through the surface tension will be different for stronger nonequilibrium processes. Structural variables are introduced to extend the thermodynamic description of interfaces to a wider range of processes. Interfacial resistivities will depend on interface temperature as well as mass and heat flux through the interface.

On paradoxical phenomena during evaporation and condensation between two parallel plates

Kinetic theory has long predicted that temperature inversion may happen in the vapor-phase for evaporation and condensation between two parallel plates, i.e., the vapor temperature at the condensation interface is higher than that at the evaporation interface. However, past studies have neglected transport in the liquid phases, which usually determine the evaporation and condensation rates. This disconnect has limited the acceptance of the kinetic theory in practical heat transfer models. In this paper, we combine interfacial conditions for mass and heat fluxes with continuum descriptions in the bulk regions of the vapor and the liquid phases to obtain a complete picture for the classical problem of evaporation and condensation between two parallel plates. The criterion for temperature inversion is rederived analytically. We also prove that the temperature jump at each interface is in the same direction as externally applied temperature difference, i.e., liquid surface is at a higher temperature than its adjacent vapor on the evaporating interface and at a lower temperature than its adjacent vapor on the condensing interface. We explain the interfacial temperature jump and temperature inversion using the interfacial cooling and heating processes, and we predict that this process can lead to a vapor phase temperature much lower than the lowest wall temperatures and much higher than the highest wall temperature imposed. When the latent heat of evaporation is small, we found that evaporation can happen at the low temperature side while condensation occurs at the high temperature side, opposing the temperature gradient.

Mesoscopic Model for Disjoining Pressure Effects in Nanoscale Thin Liquid Films and Evaporating Extended Meniscuses

Disjoining pressure effect is the key to describe contact line dynamics, micro/nanoscale liquid-vapor phase change heat transfer, and liquid transport in nanopores. In this paper, by combining a mesoscopic approach for nanoscale liquid-vapor interfacial transport and a mean-field approximation of the long-range solid-fluid molecular interaction, a mesoscopic model for the disjoining pressure effect in nanoscale thin liquid films is proposed. The capability of this model to delineate the disjoining pressure effect is validated. We demonstrate that the Hamaker constant determined from our model agrees very well with molecular dynamics (MD) simulation and that the transient evaporation/condensation mass flux predicted by this mesoscopic model is also consistent with the kinetic theory. Using this model, we investigate the characteristics of the evaporating extended meniscus in a nanochannel. The nonevaporating film region, the evaporating thin-film region, and the intrinsic meniscus region are successfully captured by our model. Our results suggest that the apparent contact angle and thickness of the nonevaporating liquid film are self-tuned according to the evaporation rate, and a higher evaporation rate results a in larger apparent contact angle and thinner nonevaporating liquid film. We also show that disjoining pressure plays a dominant role in the nonevaporating film region and suppresses the evaporation in this region, while capillary pressure dominates the intrinsic meniscus region. Strong evaporation takes place in the thin-film region, and both the disjoining pressure and capillary pressure contribute to the total pressure difference that delivers the liquid from the intrinsic meniscus region to the evaporating thin-film region, compensating for the liquid mass loss due to strong evaporation. Our work provides a new avenue for investigating thin liquid film spreading, liquid transport in nanopores, and microscopic liquid-vapor phase change heat/mass transfer mechanisms near the three-phase contact line region.

Beginner's Guide to Micro- and Nanoscale Electrochemical Additive Manufacturing

Electrochemical additive manufacturing is an advanced microfabrication technology capable of producing features of almost unlimited geometrical complexity. A unique combination of the capacity to process conductive materials, design freedom, and micro- to nanoscale resolution offered by these electrochemical techniques promises tremendous opportunities for a multitude of future applications spanning microelectronics, sensing, robotics, and energy storage. This review aims to equip readers with the basic principles of electrochemical 3D printing at the small length scale. By describing the basic principles of electrochemical additive manufacturing technology and using the recent advances in the field, this beginner's guide illustrates how controlling the fundamental phenomena that underpin the print process can be used to vary dimensions, morphology, and microstructure of printed structures.

Mechanism, Kinetics, and Potential of Mean Force of Evaporation of Water from Aqueous Sodium Chloride Solutions of Varying Concentrations

The mechanism, kinetics, and potential of mean force of evaporation of water from aqueous NaCl solutions are investigated through both unbiased molecular dynamics simulations and also biased simulations using the umbrella sampling method. The results are obtained for aqueous solutions of three different NaCl concentrations ranging from 0.6 to 6.0 m and also for pure water. The rate of evaporation is found to decrease in the presence of ions. It is found that the process of evaporation of a surface water molecule from ionic solutions can be triggered through its collision with another water or chloride ion. Such collisions provide the additional kinetic energy that is required for evaporation. However, when the collision takes place with a Cl^- ion, the evaporation of the escaping water also involves a collision with water in the vicinity of the ion at the same time along with the ion-water collision. These two collisions together provide the required kinetic energy for escape of the evaporating water molecule. Thus, the mechanism of evaporation process of ionic solutions can be more complex than that of pure water. The potential of mean force (PMF) of evaporation is found to be positive and it increases with increasing ion concentration. Also, no barrier in the PMF is found to be present for the condensation of water from vapor phase to the surfaces of the solutions. A detailed analysis of the unsuccessful evaporation attempts by surface water molecules is also made in the current study.

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.

Pore-Level Multiphase Simulations of Realistic Distillation Membranes for Water Desalination

Membrane distillation (MD) is a thermally driven separation process that is operated below boiling point. Since the performance of MD modules is still comparatively low, current research aims to improve the understanding of the membrane structure and its underlying mechanisms at the pore level. Based on existing realistic 3D membrane geometries (up to 0.5 billion voxels with 39nm resolution) obtained from ptychographic X-ray computed tomography, the D3Q27 lattice Boltzmann (LB) method was used to investigate the interaction of the liquid and gaseous phase with the porous membrane material. In particular, the Shan and Chen multi-phase model was used to simulate multi-phase flow at the pore level. We investigated the liquid entry pressure of different membrane samples and analysed the influence of different micropillar structures on the Wenzel and Cassie-Baxter state of water droplets on rough hydrophobic surfaces. Moreover, we calculated the liquid entry pressure required for entering the membrane pores and extracted realistic water contact surfaces for different membrane samples. The influence of the micropillars and flow on the water-membrane contact surface was investigated. Finally, we determined the air-water interface within a partially saturated membrane, finding that the droplet size and distribution correlated with the porosity of the membrane.

Formation of Crystalline Bulk Heterojunctions in Organic Solar Cells: Insights from Phase-Field Simulations

The performance of organic solar cells strongly depends on the bulk-heterojunction (BHJ) morphology of the photoactive layer. This BHJ forms during the drying of the wet-deposited solution, because of physical processes such as crystallization and/or liquid-liquid phase separation (LLPS). However, the process-structure relationship remains insufficiently understood. In this work, a recently developed, coupled phase-field-fluid mechanics framework is used to simulate the BHJ formation upon drying. For the first time, this allows to investigate the interplay between all the relevant physical processes (evaporation, crystal nucleation and growth, liquid demixing, composition-dependent kinetic properties), within a single coherent theoretical framework. Simulations for the model system P3HT-PCBM are presented. The comparison with previously reported in situ characterization of the drying structure is very convincing: The morphology formation pathways, crystallization kinetics, and final morphology are in line with experimental results. The final BHJ morphology is a subtle mixture of pure crystalline donor and acceptor phases, pure and mixed amorphous domains, which depends on the process parameters and material properties. The expected benefit of such an approach is to identify physical design rules for ink formulation and processing conditions to optimize the cell's performance. It could be applied to recent organic material systems in the future.

Enhanced Water Evaporation from Å-Scale Graphene Nanopores

Enhancing the kinetics of liquid-vapor transition from nanoscale confinements is an attractive strategy for developing evaporation and separation applications. The ultimate limit of confinement for evaporation is an atom thick interface hosting angstrom-scale nanopores. Herein, using a combined experimental/computational approach, we report highly enhanced water evaporation rates when angstrom sized oxygen-functionalized graphene nanopores are placed at the liquid-vapor interface. The evaporation flux increases for the smaller nanopores with an enhancement up to 35-fold with respect to the bare liquid-vapor interface. Molecular dynamics simulations reveal that oxygen-functionalized nanopores render rapid rotational and translational dynamics to the water molecules due to a reduced and short-lived water-water hydrogen bonding. The potential of mean force (PMF) reveals that the free energy barrier for water evaporation decreases in the presence of nanopores at the atomically thin interface, which further explains the enhancement in evaporation flux. These findings can enable the development of energy-efficient technologies relying on water evaporation.

Evaporation of Water Nanodroplets on Heated Surfaces: Does Nano Matter?

While experiments and continuum models have provided a relatively good understanding of the evaporation of macroscopic water droplets, elucidating how sessile nanodroplets evaporate is an open question critical for advancing nanotechnological applications where nanodroplets can play an essential role. Here, using molecular dynamics simulations, we find that evaporating nanodroplets, in contrast to their macroscopic counterparts, are not always in thermal equilibrium with the substrate and that the vapor concentration on the nanodroplet surface does not reach a steady state. As a result, the evaporative behavior of nanodroplets is significantly different. Regardless of hydrophobicity, nanodroplets do not follow conventional evaporation modes but instead exhibit dynamic wetting behavior characterized by huge, non-equilibrium, isovolumetric fluctuations in the contact angle and contact radius. For hydrophilic nanodroplets, the evaporation rate, controlled by the vapor concentration, decays exponentially over time. Hydrophobic nanodroplets follow stretched exponential kinetics arising from the slower thermalization with the substrate. The evaporative half-lifetime of the nanodroplets is directly related to the thermalization time scale and therefore increases monotonically with the hydrophobicity of the substrate. Finally, the evaporative flux profile along the nanodroplet surface is highly nonuniform but does not diverge at the contact line as the macroscopic continuum models predict.

Temperature measurements of liquid flat jets in vacuum

Sub-μm thin samples are essential for spectroscopic purposes. The development of flat micro-jets enabled novel spectroscopic and scattering methods for investigating molecular systems in the liquid phase. However, the temperature of these ultra-thin liquid sheets in vacuum has not been systematically investigated. Here, we present a comprehensive temperature characterization using optical Raman spectroscopy of sub-micron flatjets produced by two different methods: colliding of two cylindrical jets and a cylindrical jet compressed by a high pressure gas. Our results reveal the dependence of the cooling rate on the material properties and the source characteristics, i.e., nozzle-orifice size, flow rate, and pressure. We show that materials with higher vapor pressures exhibit faster cooling rates, which is illustrated by comparing the temperature profiles of water and ethanol flatjets. In a sub-μm liquid sheet, the temperature of the water sample reaches around 268 K and the ethanol around 253 K close to the flatjet's terminus.

Phase-Field Simulation of Liquid-Vapor Equilibrium and Evaporation of Fluid Mixtures

In solution processing of thin films, the material layer is deposited from a solution composed of several solutes and solvents. The final morphology and hence the properties of the film often depend on the time needed for the evaporation of the solvents. This is typically the case for organic photoactive or electronic layers. Therefore, it is important to be able to predict the evaporation kinetics of such mixtures. We propose here a new phase-field model for the simulation of evaporating fluid mixtures and simulate their evaporation kinetics. Similar to the Hertz-Knudsen theory, the local liquid-vapor (LV) equilibrium is assumed to be reached at the film surface and evaporation is driven by diffusion away from this gas layer. In the situation where the evaporation is purely driven by the LV equilibrium, the simulations match the behavior expected theoretically from the free energy: for evaporation of pure solvents, the evaporation rate is constant and proportional to the vapor pressure. For mixtures, the evaporation rate is in general strongly time-dependent because of the changing composition of the film. Nevertheless, for highly nonideal mixtures, such as poorly compatible fluids or polymer solutions, the evaporation rate becomes almost constant in the limit of low Biot numbers. The results of the simulation have been successfully compared to experiments on a polystyrene-toluene mixture. The model allows to take into account deformations of the liquid-vapor interface and, therefore, to simulate film roughness or dewetting.

Accidental Combustion Phenomena at Cryogenic Conditions

The presented state of the art can be intended as an overview of the current understandings and the remaining challenges on the phenomenological aspects involving systems operating at ultra-low temperature, which typically characterize the cryogenic fuels, i.e., liquefied natural gas and liquefied hydrogen. To this aim, thermodynamic, kinetic, and technological aspects were included and integrated. Either experimental or numerical techniques currently available for the evaluation of safety parameters and the overall reactivity of systems at cryogenic temperatures were discussed. The main advantages and disadvantages of different alternatives were compared. Theoretical background and suitable models were reported given possible implementation to the analyzed conditions. Attention was paid to models describing peculiar phenomena mainly relevant at cryogenic temperatures (e.g., para-to-ortho transformation and thermal stratification in case of accidental release) as well as critical aspects involving standard phenomena (e.g., ultra-low temperature combustion and evaporation rate).

Remarkable Rate of Water Evaporation through Naked Veins of Natural Tree Leaves

In the form of leaves, nature designs the finest photothermal evaporators, and the tremendous evaporation efficiency of leaves is supported by a precisely designed network of veins. Here, we have demonstrated that the vein network of a natural leaf can be extracted through a simple water-assisted digestion process and exploited for low-energy steam generation. The naked leaf veins exhibit a remarkable flux (evaporation rate, 1.5 kg·m-2·h-1) of capillary evaporation under ambient conditions (25 °C and 30% RH), close to the photothermal material-based evaporators reported in the recent literature. Even inside a dark box, naked veins exhibit an evaporation rate up to 4.5 kg·m-2·h-1 (at 30% relative humidity (RH) and a wind speed of 22 km·h-1). The mechanistic studies performed with variable atmospheric conditions (temperature, humidity, and wind speed) suggest the evaporation process through the naked veins to be a kinetic-limited process. Naked veins with remarkable evaporation efficiency are found to be suitable for applications like water desalination and streaming potential harvesting. Experiments with the naked veins also unveiled that the biofluidic channels in leaves not only exhibit the characteristics of surface charge-governed ionic transport but also support an exceptional water transport velocity of 1444 μm·s-1.

Evaporation and propagation of liquid drop streams at vacuum pressures: Experiments and modeling

Evaporation of streams of liquid droplets in environments at vacuum pressures below the vapor pressure has not been widely studied. Here, experiments and simulations are reported that quantify the change in droplet diameter when a steady stream of ≈100 μm diameter drops is injected into a chamber initially evacuated to <10^{-8}bar. In experiments, droplets fall through the center of a 0.8 m long liquid nitrogen cooled shroud, simulating infinity radiation and vapor mass flux boundary conditions. Experimentally measured changes in drop diameters vary from ≈0 to 6 μm when the initial vapor pressure is increased from 10^{-6} to 10^{-3} bar by heating the liquid. Measured diameter changes are predicted by a model based on the Hertz-Knudsen equation. One uncertainty in the calculation is the "sticking coefficient" β. Assuming a constant β for all conditions studied here, predicted diameter changes best match measurements with β≈0.3. This value falls within the range of β reported in the literature for organic liquids. Finally, at the higher vapor pressure conditions considered here, the drop stream disperses transverse to the main flow direction. This spread is attributed to forces imparted by an absolute pressure gradient produced by the evaporating stream.

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.

Shining New Light on the Kinetics of Water Uptake by Organic Aerosol Particles

The uptake of water vapor by various organic aerosols is important in a number of applications ranging from medical delivery of pharmaceutical aerosols to cloud formation in the atmosphere. The coefficient that describes the probability that the impinging gas-phase molecule sticks to the surface of interest is called the mass accommodation coefficient, α_M. Despite the importance of this coefficient for the description of water uptake kinetics, accurate values are still lacking for many systems. In this Feature Article, we present various experimental techniques that have been evoked in the literature to study the interfacial transport of water and discuss the corresponding strengths and limitations. This includes our recently developed technique called photothermal single-particle spectroscopy (PSPS). The PSPS technique allows for a retrieval of α_M values from three independent, yet simultaneous measurements operating close to equilibrium, providing a robust assessment of interfacial mass transport. We review the currently available data for α_M for water on various organics and discuss the few studies that address the temperature and relative humidity dependence of α_M for water on organics. The knowledge of the latter, for example, is crucial to assess the water uptake kinetics of organic aerosols in the Earth's atmosphere. Finally, we argue that PSPS might also be a viable method to better restrict the α_M value for water on liquid water.

Numerical Investigation of Spray Collapse in GDI with OpenFOAM

During certain operating conditions in spark-ignited direct injection engines (GDI), the injected fuel will be superheated and begin to rapidly vaporize. Fast vaporization can be beneficial for fuel–oxidizer mixing and subsequent combustion, but it poses the risk of spray collapse. In this work, spray collapse is numerically investigated for a single hole and the spray G eight-hole injector of an engine combustion network (ECN). Results from a new OpenFOAM solver are first compared against results of the commercial CONVERGE software for single-hole injectors and validated. The results corroborate the perception that the superheat ratio Rp, which is typically used for the classification of flashing regimes, cannot describe spray collapse behavior. Three cases using the eight-hole spray G injector geometry are compared with experimental data. The first case is the standard G2 test case, with iso-octane as an injected fluid, which is only slightly superheated, whereas the two other cases use propane and show spray collapse behavior in the experiment. The numerical results support the assumption that the interaction of shocks due to the underexpanded vapor jet causes spray collapse. Further, the spray structures match well with experimental data, and shock interactions that provide an explanation for the observed phenomenon are discussed.

A purely kinetic description of the evaporation of water droplets

The process of water evaporation, although deeply studied, does not enjoy a kinetic description that captures known physics and can be integrated with other detailed processes such as drying of catalytic membranes embedded in vapor-fed devices and chemical reactions in aerosol whose volumes are changing dynamically. In this work, we present a simple, three-step kinetic model for water evaporation that is based on theory and validated by using well-established thermodynamic models of droplet size as a function of time, temperature, and relative humidity as well as data from time-resolved measurements of evaporating droplet size. The kinetic mechanism for evaporation is a combination of two limiting processes occurring in the highly dynamic liquid-vapor interfacial region: direct first order desorption of a single water molecule and desorption resulting from a local fluctuation, described using third order kinetics. The model reproduces data over a range of relative humidities and temperatures only if the interface that separates bulk water from gas phase water has a finite width, consistent with previous experimental and theoretical studies. The influence of droplet cooling during rapid evaporation on the kinetics is discussed; discrepancies between the various models point to the need for additional experimental data to identify their origin.

Dynamical Density Functional Theory for the Drying and Stratification of Binary Colloidal Dispersions

We develop a dynamical density functional theory based model for the drying of colloidal films on planar surfaces. We consider mixtures of two different sizes of hard-sphere colloids. Depending on the solvent evaporation rate and the initial concentrations of the two species, we observe varying degrees of stratification in the final dried films. Our model predicts the various structures described in the literature previously from experiments and computer simulations, in particular the small-on-top stratified films. Our model also includes the influence of adsorption of particles to the interfaces.

Distributions of Two Atoms Collisions over the Surface of the Condensed Phase

The processes of heat and mass transfer are closely related to the evaporation of a substance from the surface of the condensed phase. The interaction of outgoing molecules from the surface of the condensed phase with condensed phase molecules plays a fundamental role. A simpler case of evaporation is the departure of atoms from the surface of the condensed phase, i.e. the atoms overcome the potential barrier on the surface of the condensed phase. Depending on the evaporation rate, a Knudsen layer appears above the surface of the condensed phase. In this paper, based on the model of rigid spheres, the density distributions of the collision distances and the average values of the collision distances of two atoms emitted simultaneously from the surface of the condensed phase above the surface are analyzed. Distributions of the collision distance depending on the surface temperature, the size of the potential barrier, and the size of the evaporation area are obtained. Computer experiments were performed using the Monte Carlo method. To obtain the results of numerical simulation, a parallel algorithm adapted to calculations on graphics processors with CUDA technology was developed.

Maximum evaporating flux of molecular fluids from a planar liquid surface

In this work, we use the kinetic theory of gases (KTG) to develop a theoretical model to understand the role of internal motions of molecules on the maximum evaporation flux from a planar liquid surface. The kinetic theory is applied to study the evaporation of molecular fluids into a vacuum and predict the dimensionless maximum evaporation flux (J_{R,max}, i.e., the ratio of the maximum evaporation flux to the molar flux emitted from a liquid surface). The key assumptions regarding the velocity distribution function (VDF) of polyatomic molecules in the highly nonequilibrium vapor near the evaporating surface are validated by the VDF obtained directly from molecular dynamics (MD) simulations. Our KTG-based analysis shows that J_{R,max} is affected by the specific heat (c_{V,int}) associated with internal degrees of freedom of fluid molecules. When the maximum evaporation flux is reached, the isotropic evaporating vapor far from the liquid surface moves at its speed of sound regardless of whether it is a monatomic vapor or polyatomic vapor. To fundamentally understand the evaporation of a molecular fluid into a vacuum, we solve the Boltzmann transport equation (BTE) to obtain the temperature, density, and flow speed distributions in the highly nonequilibrium evaporating vapor flow. Our BTE solutions indicate that there are several universal features of the evaporating vapor when the maximum evaporation flux occurs. In particular, we find that the evaporating vapor flow speed reaches the maximum value of sqrt[1.5] times the most probable thermal speed in the vapor flow direction at the vacuum boundary, and this maximum value is independent of fluid properties. All theoretical predictions in this work are verified by the MD simulation results of the evaporation of the model liquid Ar and the model liquid n-dodecane into a vacuum, and existing experimental data.

How Adsorption of Pheromones on Aerosols Controls Their Transport

We propose a general transport theory for pheromone molecules in an atmosphere containing aerosols. Many pheromones are hydrophobic molecules containing polar groups. They are low volatile and have some properties similar to those of hydrotropes. They therefore form a nonsoluble film at the water-air interface of aerosols. The fate of a small pheromone puff in air is computed through reaction-diffusion equations. Partitioning of pheromones between the gas and the aerosol surface over time is studied for various climate conditions (available aerosol surface) and adsorption affinities (energy of adsorption). We show that, for adsorption energy above 30 k B T per molecule, transport of pheromones on aerosols dominates over molecular transport typically 10 s after pheromone emission, even when few adsorbing aerosols are present. This new communication path for airborne chemicals leads to distinctive features including enhanced signal sensibility and increased persistence of pheromone concentration in the air due to slow diffusion of aerosols. Each aerosol droplet has the ability to adsorb thousands of pheromones to the surface, keeping a "history" of the atmospheric content between emission and reception. This new mechanism of pheromone transport leads to dramatic consequences on insect sensing revisiting the way we figure the capture of chemical signals.

Thermal transport across the interface between liquid n-dodecane and its own vapor: A molecular dynamics study

There are two possible thermal transport mechanisms at liquid-gas interfaces, namely, evaporation/condensation (i.e., heat transfer by liquid-vapor phase change at liquid surfaces) and heat conduction (i.e., heat exchange by collisions between gas molecules and liquid surfaces). Using molecular dynamics (MD) simulations, we study thermal transport across the liquid-vapor interface of a model n-dodecane (C₁₂H₂₆) under various driving force conditions. In each MD simulation, we restrict the thermal energy to be transferred across the liquid-vapor interface by only one mechanism. In spite of the complex intramolecular interactions in n-dodecane molecules, our modeling results indicate that the Schrage relationships, which were shown to give accurate predictions of evaporation and condensation rates of monatomic fluids, are also valid in the prediction of evaporation and condensation rates of n-dodecane. In the case of heat conduction at the liquid-vapor interface of n-dodecane, the interfacial thermal conductance obtained from MD simulations is consistent with the prediction from the kinetic theory of gases. The fundamental understanding of thermal transport mechanisms at liquid-gas interfaces will allow us to formulate appropriate boundary conditions for continuum modeling of heating and evaporation of small fuel droplets.