Verification of Onsager's reciprocal relations for evaporation and condensation using non-equilibrium molecular dynamics

Non-equilibrium molecular dynamic (NEMD) simulations have been used to study heat and mass transfer across a vapor-liquid interface for a one-component system using a Lennard-Jones spline potential. It was confirmed that the relation between the surface tension and the surface temperature in the non-equilibrium system was the same as in equilibrium (local equilibrium). Interfacial transfer coefficients were evaluated for the surface, which expressed the heat and mass fluxes in temperature and chemical potential differences across the interfacial region (film). In this analysis it was assumed that the Onsager reciprocal relations were valid. In this paper we extend the number of simulations such that we can calculate all four interface film transfer coefficients along the whole liquid-vapor coexistence curve. We do this analysis both for the case where we use the measurable heat flux on the vapor side and for the case where we use the measurable heat flux on the liquid side. The most important result we found is that the coupling coefficients within the accuracy of the calculation are equal. This is the first verification of the validity of the Onsager relations for transport through a surface using molecular dynamics. The interfacial film transfer coefficients are found to be a function of the surface temperature alone. New expressions are given for the kinetic theory values of these coefficients which only depend on the surface temperature. The NEMD values were found to be in good agreement with these expressions.     

Dynamic behavior of interfaces: Modeling with nonequilibrium thermodynamics

In multiphase systems the transfer of mass, heat, and momentum, both along and across phase interfaces, has an important impact on the overall dynamics of the system. Familiar examples are the effects of surface diffusion on foam drainage (Marangoni effect), or the effect of surface elasticities on the deformation of vesicles or red blood cells in an arterial flow. In this paper we will review recent work on modeling transfer processes associated with interfaces in the context of nonequilibrium thermodynamics (NET). The focus will be on NET frameworks employing the Gibbs dividing surface model, in which the interface is modeled as a two-dimensional plane. This plane has excess variables associated with it, such as a surface mass density, a surface momentum density, a surface energy density, and a surface entropy density. We will review a number of NET frameworks which can be used to derive balance equations and constitutive models for the time rate of change of these excess variables, as a result of in-plane (tangential) transfer processes, and exchange with the adjoining bulk phases. These balance equations must be solved together with mass, momentum, and energy balances for the bulk phases, and a set of boundary conditions coupling the set of bulk and interface equations. This entire set of equations constitutes a comprehensive continuum model for a multiphase system, and allows us to examine the role of the interfacial dynamics on the overall dynamics of the system. With respect to the constitutive equations we will focus primarily on equations for the surface extra stress tensor.     

Resistances for heat and mass transfer through a liquid-vapor interface in a binary mixture

In this paper we calculate the interfacial resistances to heat and mass transfer through a liquid-vapor interface in a binary mixture. We use two methods, the direct calculation from the actual nonequilibrium solution and integral relations, derived earlier. We verify, that integral relations, being a relatively faster and cheaper method, indeed gives the same results as the direct processing of a nonequilibrium solution. Furthermore we compare the absolute values of the interfacial resistances with the ones obtained from kinetic theory. Matching the diagonal resistances for the binary mixture we find that kinetic theory underestimates the cross coefficients. The heat of transfer is, as a consequence, correspondingly larger.     

Surface conservation laws at microscopically diffuse interfaces

In studies of interfaces with dynamic chemical composition, bulk and interfacial quantities are often coupled via surface conservation laws of excess surface quantities. While this approach is easily justified for microscopically sharp interfaces, its applicability in the context of microscopically diffuse interfaces is less theoretically well-established. Furthermore, surface conservation laws (and interfacial models in general) are often derived phenomenologically rather than systematically. In this article, we first provide a mathematically rigorous justification for surface conservation laws at diffuse interfaces based on an asymptotic analysis of transport processes in the boundary layer and derive general formulae for the surface and normal fluxes that appear in surface conservation laws. Next, we use nonequilibrium thermodynamics to formulate surface conservation laws in terms of chemical potentials and provide a method for systematically deriving the structure of the interfacial layer. Finally, we derive surface conservation laws for a few examples from diffusive and electrochemical transport.     

Evaporation-induced particle microseparations inside droplets floating on a chip

We describe phenomena of colloidal particle transport and separation inside single microdroplets of water floating on the surface of dense fluorinated oil. The experiments were performed on microfluidic chips, where single droplets were manipulated with alternating electric fields applied to arrays of electrodes below the oil. The particles suspended in the droplets were collected in their top region during the evaporation process. Experimental results and numerical simulations show that this microsepration occurs as a result of a series of processes driven by mass and heat transfer. An interfacial tension gradient develops on the surface of the droplet as a result of the nonuniform temperature distribution during the evaporation. This gradient generates an internal convective Marangoni flow. The colloidal particles transported by the flow are collected in the top of the droplets by the hydrodynamic flux, compensating for evaporation through the exposed top surface. The internal flow pattern and temperature distribution within evaporating droplets were simulated using finite element calculations. The results of the simulation were consistent with experiments using tracer particles. Such microseparation processes can be used for on-chip synthesis of advanced particles and innovative microbioassays.     

A Heterogeneous Multiscale Dynamic Model for Simulation of Catalytic Reforming Reactors

The present work aims to establish a generic reforming reaction scheme to evaluate the performance of catalytic reforming systems with the aid of a one‐dimensional heterogeneous dynamic model. The novelty of the numerical model stems from the direct inclusion of interphase (fluid‐to‐particle surface), intraparticle (within particle), and intrareactor heat and mass transport resistances under transient conditions. The developed model accounts for the multicomponent gas mixture physicochemical properties and correlations for calculating mass and heat transfer coefficients. Effective macroscopic properties within the particle are calculated by incorporating diffusivities and conductivities of the porous network characteristics accounting for Knudsen and molecular transport as well as tortuosity and porosity of the overall porous path. The industrial case of a steam‐methane reforming multitubular reactor was studied as the most representative case of the generic reaction scheme, with all mass/energy resistances present under severe pressure and temperature conditions. It was shown that there are notable diffusional limitations within the particle, whereas there are also temperature and partial pressure gradients due to the heat and mass transport resistances in the particle film layer. It is further demonstrated that the proposed model can be utilized as a versatile design tool for catalytic reactor development and optimization.     

The Prediction of Surface Tension and Thermodynamic Analysis of the Surface in Mixtures of Cryogenic Liquids

The surface tensions of 42 binary cryogenic mixtures at low temperature are correlated using the Shereshefsky model and excellent results are obtained. The average percent deviation is about ～ 1.08%. The Gibbs energy change in the surface region is calculated and is used to obtain the excess number of molecular layers in the surface region. Furthermore, the model is used to derive an equation for the standard Gibbs energy of adsorption. The experimental standard Gibbs energy of adsorption is obtained from surface tension data and compared with calculated data. The agreement between experimental and calculated data is found to be very good. The magnitude of the Gibbs energy change in the surface region and the standard Gibbs energy of adsorption are discussed in terms of nature and type of intermolecular interactions in binary mixtures.     

Protein diffusion across the interface in aqueous two-phase systems

We present a detailed study of the diffusive transport of proteins across a fluid phase boundary within aqueous two-phase systems. The aim of the work is to investigate whether local effects at the phase boundary cause a retardation of the diffusive transport between the phases. Possible modifications of interfacial mass transfer could be due to protein adsorption at the phase boundary or local electric fields from electric double layers. Experiments with a microfluidic system have been performed in which protein diffusion (bovine serum albumin and ovalbumin) within a bilaminated configuration of two phases containing polyethylene glycol and dextran is analyzed. A one-dimensional model incorporating phase-specific diffusion constants and the difference in chemical potential between the phases has been formulated. A comparison of experimental and simulation data shows a good overall agreement and suggests that a potential local influence of the phase boundary on protein transport is insignificant for the systems under investigation.     

Local non-equilibrium distribution of charge carriers in a phase-coherent conductor

We use the scattering matrix approach to derive generalized Bardeen-like formulae for the conductances between the contacts of a phase-coherent multiprobe conductor and a tunneling tip which probes its surface. These conductances are proportional to local partial densities of states, called injectivities and emissivities. The current and the current fluctuations measured at the tip are related to and effective local non-equilibrium distribution function. This distribution function contains the quantum-mechanical phase-coherence of the charge carriers in the conductor and is given as products of injectivities and the Fermi distribution functions in the electron reservoirs. The results are illustrated for measurements on ballistic conductors with barriers and for diffusive conductors.     

Indications for a surface temperature excess in heterogeneous catalysis

We explored the surface temperature excess DeltaT = T(r) - T(m) (real reaction temperature T(r), measured catalyst temperature T(m)) on the basis of experimental data, a gradually curved Arrhenius plot for CO oxidation reactions over Pd/gamma-Al(2)O(3) catalysts. Such a plot could be an indication of the surface temperature excess in the 2-dimensional reaction surfaces of catalysts. The positive or negative surface temperature excess could be developed to be a general explanation for a gradually curved Arrhenius plot of a gas-solid catalytic system. This is a new insight into solving the puzzle on such common phenomena in heterogeneous catalysis. By using the reciprocal of the real reaction temperature T(r) in the hypothetical 2-D reaction surface, instead of the experimentally determined catalyst temperature T(m) or the gas temperature T(g), the gradually curved Arrhenius plot becomes linear. We investigated the implications of such a difference among T(r), T(m), and T(g). The surface temperature excess could be the effect of coupling between the fluxes of a chemical reaction and heat transport in the 2-D reaction surface. Its order of magnitude is 10 K for the present model system.(1) The surface temperature excess increases exponentially with the reaction temperature.     

Molecular dynamics simulations of mass transfer resistance in grain boundaries of twinned zeolite membranes

The net mass transfer resistance for gas molecules permeating through zeolite membranes includes contributions from intracrystalline diffusion and contributions from interfacial effects. These interfacial effects can arise either from gas-zeolite interfaces or from interfaces that exist within zeolite crystals due to grain boundaries. We present the first atomically detailed simulations that examine interfacial mass transfer resistance due to internal grain boundaries in zeolites that are relevant for membrane applications. Our calculations examine twinned silicalite crystals in crystallographic configurations that have been identified in previous experiments. We used the dual control volume grand canonical molecular dynamics method to simulate the permeance of CH(4) and CF(4) through thin twinned silicalite crystals. The magnitudes of the grain boundary resistances are quite substantial, at least for the thin crystals that are accessible in our simulations.     

Molecular dynamics of phenol at the liquid-vapor interface of water

Molecular dynamics results are presented for phenol at the water liquid-vapor interface at 300 K. The calculated excess free energy of phenol at the interface is -2.8 +/- 0.4 kcal/mol, in good agreement with the recent experimental results of Eisenthal and co-workers. The most probable orientation of the phenol molecule at the surface is such that the aromatic ring is perpendicular to the interface and the OH group is fully immersed in water. The hydroxyl substituent has a preferred orientation which is similar to the orientation of OH bonds of water at the pure water liquid-vapor interface. The transition between interfacial and bulk-like behavior of phenol is abrupt and occurs when the center of mass of the solute is located about 6 angstroms from the Gibbs surface of water. In this region the para carbon atom of the hydrophobic benzene ring can reach the interface and become partially dehydrated. This result suggests that the width of the interfacial region in which the behavior of a simple amphiphilic solute in water is influenced by the presence of the surface depends primarily on the size of its hydrophobic part. The role of the OH substituent was investigated by comparing phenol at the interface with two model systems: benzene with and without partial charges on carbon and hydrogen atoms. It is shown that in the absence of the hydrophilic substituent the solute is located further away from the liquid phase and is more likely to be oriented parallel to the interface. However, when the center of mass of the solute is moved into the interfacial region where the density of water approaches that of the bulk solvent, all three molecules become oriented perpendicularly to the surface. In this orientation the work of cavity formation needed to accommodate the hydrophobic ring in aqueous solvent is minimized.     

Vapor-liquid interfacial properties of mutually saturated water/1-butanol solutions

Adsorption and ordering at the vapor-liquid interfaces of mutually saturated water/1-butanol solutions at a temperature of 298.15 K were investigated using configurational-bias Monte Carlo simulations in the Gibbs ensemble and compared to the surface properties of neat water and 1-butanol liquids. A dense 1-butanol monolayer is observed at the surface of the water-rich phase, which results in a substantial decrease of its surface tension. In contrast, there is no enrichment of water molecules at the surface of the butanol-rich phase, and its surface tension is not significantly changed. Analysis of the interfacial structures reveals that these systems exhibit orientational ordering and composition heterogeneity. Analysis of the hydrogen-bonding distributions suggests that the formation of the 1-butanol monolayer is driven by an excellent match between water and the primary alcohol; that is, additional hydrogen bonds are formed between the excess free hydrogens of surface water and the excess hydrogen-bond acceptor sites of 1-butanol.