Molecular-dynamics simulation of evaporation processes of fluid bridges confined in slitlike pores

Mass Transfer through Vapor-Liquid Interfaces Studied by Non-Stationary Molecular Dynamics Simulations

Molecular dynamics (MD) simulations are highly attractive for studying the influence of interfacial effects, such as the enrichment of components, on the mass transfer through the interface. In a recent work, we have presented a steady-state MD simulation method for investigating this phenomenon and tested it using model mixtures with and without interfacial enrichment. The present study extends this work by introducing a non-stationary MD simulation method. A rectangular simulation box that contains a mixture of two components 1 + 2 with a vapor phase in the middle and two liquid phases on both sides is used. Starting from a vapor-liquid equilibrium state, a non-stationary molar flux of component 2 is induced by inserting particles of component 2 into the center of the vapor phase in a pulse-like manner. During the isothermal relaxation process, particles of component 2 pass through the vapor phase, cross the vapor-liquid interface, and enter the liquid phase. The system thereby relaxes into a new vapor-liquid equilibrium state. During the relaxation process, spatially resolved responses for the component densities, fluxes, and pressure are sampled. To reduce the noise and provide measures for the uncertainty of the observables, a set of replicas of simulations is carried out. The new simulation method was applied to study mass transfer in two binary Lennard-Jones mixtures: one that exhibits a strong enrichment of the low-boiling component 2 at the vapor-liquid interface and one that shows no enrichment. Even though both mixtures have similar transport coefficients in the bulk phases, the results for mass transfer differ significantly, indicating that the interfacial enrichment influences the mass transfer.

Prediction of the structure and mechanical properties of polycaprolactone-silica nanocomposites and the interphase region by molecular dynamics simulations: the effect of PEGylation

Polymer/silica (PS) nanocomposites are, among numerous combinations of inorganic/organic nanocomposites, one of the most important materials reported in the literature and have been employed in a wide variety of applications. Due to this great interest in the scientific and industry community, knowledge about their physiochemistry allows for a better understanding of their development and improvement. One area of interest found in biopolymers is silica, where silica nanoparticles can be used to increase their mechanical properties and give them higher opportunities to replace synthetic plastics. With this aim in mind, molecular dynamics (MD) simulations were used to predict the structure and mechanical properties of the interphase region and nanocomposite systems of polycaprolactone (PCL), a common poly(hydroxy acid) type biopolymer, reinforced with silica nanoparticles. Two types of nanoparticles were studied to assess the effect of PEGylation: hydroxyl (ungrafted) and polyethylene glycol (PEG) (grafted or PEGylated) functionalized silica. The interaction energy between the nanoparticle and the polymeric matrix was determined, showing an increase of the affinity between each component due to the PEGylation of the nanoparticle. Through the analysis of polymer density profiles, the structure and thickness of the interphase region were determined, and it was observed that PEGylation increased the interphase thickness from 10.80 Å to 13.04 Å while it decreased the peak and average polymer density of the interphase region. Using compressed and expanded molecular models of the neat PCL polymer, the mechanical properties of the interphase region were related to its density through an interpolation model, and mechanical property profiles were obtained, from which the average values of the Young's modulus, Poisson's ratio and shear modulus of the interphase region were calculated. Finally, the mechanical properties of the nanocomposites were determined by molecular mechanics simulations, showing that the silica nanoparticles increased the stiffness of the composite system to about 7-8% with respect to that of the neat polymer, having a 2.09% weight of bare silica or 2.82% weight of PEGylated silica. PEGylation did not show an additional effect on the overall mechanical properties. A mean field micromechanics model (Mori-Tanaka) corroborated the properties calculated for the interphase region using MD simulations. It was concluded that the PEGylation of the nanoparticle improved the affinity, and thus the dispersion, of the silica nanoparticles towards the PCL matrix, but with no further increase in the mechanical properties of the composite.

Interfacial phenomena at a surface of individual and complex fumed nanooxides

Investigations of interfacial and temperature behaviors of nonpolar and polar adsorbates interacting with individual and complex fumed metal or metalloid oxides (FMO), initial and subjected to various treatments or chemical functionalization and compared to such porous adsorbents as silica gels, precipitated silica, mesoporous ordered silicas, filled polymeric composites, were analyzed. Complex nanooxides include core-shell nanoparticles, CSNP (50-200nm in size) with titania or alumina cores and silica or alumina shells in contrast to simple and smaller nanoparticles of individual FMO. CSNP could be destroyed under high-pressure cryogelation (HPCG) or mechanochemical activation (MCA). These treatments affect the structure of aggregates of nanoparticles and agglomerates of aggregates, resulting in their becoming more compacted. The analysis shows that complex FMO could be more sensitive to external actions than simple nanooxides such as fumed silica. Any treatment of 'soft' FMO affects the interfacial and temperature behaviors of polar and nonpolar adsorbates. Rearrangement of secondary particles and surface functionalization affects the freezing-melting point depression of adsorbates. For some adsorbates, open hysteresis loops became readily apparent in adsorption-desorption isotherms. Clustering of adsorbates bound in textural pores in aggregates of nanoparticles (i.e., voids between nanoparticles in secondary structures) causes reduced changes in enthalpy during phase transitions (freezing, fusion, evaporation). Freezing point depression and melting point elevation cause significant hysteresis freezing-melting effects for adsorbates bound to FMO in the textural pores. Relaxation phenomena for both low- and high-molecular weight adsorbates or filled polymeric composites are affected by the morphology of primary particles, structural organization of secondary particles of differently treated or functionalized FMO, content of adsorbates, co-adsorption order, and temperature.

Contact angle of sessile drops in Lennard-Jones systems

Molecular dynamics simulations are used for studying the contact angle of nanoscale sessile drops on a planar solid wall in a system interacting via the truncated and shifted Lennard-Jones potential. The entire range between total wetting and dewetting is investigated by varying the solid-fluid dispersive interaction energy. The temperature is varied between the triple point and the critical temperature. A correlation is obtained for the contact angle in dependence of the temperature and the dispersive interaction energy. Size effects are studied by varying the number of fluid particles at otherwise constant conditions, using up to 150,000 particles. For particle numbers below 10,000, a decrease of the contact angle is found. This is attributed to a dependence of the solid-liquid surface tension on the droplet size. A convergence to a constant contact angle is observed for larger system sizes. The influence of the wall model is studied by varying the density of the wall. The effective solid-fluid dispersive interaction energy at a contact angle of θ = 90° is found to be independent of temperature and to decrease linearly with the solid density. A correlation is developed that describes the contact angle as a function of the dispersive interaction, the temperature, and the solid density. The density profile of the sessile drop and the surrounding vapor phase is described by a correlation combining a sigmoidal function and an oscillation term.

Effects of viscosity, surface tension, and evaporation rate of solvent on dry colloidal structures: a lattice Boltzmann study

Understanding the mechanisms of how colloidal solution properties and drying processes result in dry colloidal structures is essential for industrial applications such as paint, ceramics, and electrodes. In this study, we develop a computational method to simulate the drying process of colloidal suspensions containing solid particles and polymers. The method consists of a solvent evaporation model, a fluid particle dynamics method, and a two-phase lattice Boltzmann method. We determine that a high-viscosity solvent, small surface tension, and a high evaporation rate of the solvent lead to a structure with dispersed particles and interconnected pores. When these conditions are not present, the particles agglomerate and the pores are disconnected.

Molecular dynamics simulations for the motion of evaporative droplets driven by thermal gradients along nanochannels

For a one-component fluid on a solid substrate, a thermal singularity may occur at the contact line where the liquid-vapor interface intersects the solid surface. Physically, the liquid-vapor interface is almost isothermal at the liquid-vapor coexistence temperature in one-component fluids while the solid surface is almost isothermal for solids of high thermal conductivity. Therefore, a temperature discontinuity is formed if the two isothermal interfaces are of different temperatures and intersect at the contact line. This leads to the so-called thermal singularity. The localized hydrodynamics involving evaporation/condensation near the contact line leads to a contact angle depending on the underlying substrate temperature. This dependence has been shown to lead to the motion of liquid droplets on solid substrates with thermal gradients (Xu and Qian 2012 Phys. Rev. E 85 061603). In the present work, we carry out molecular dynamics (MD) simulations as numerical experiments to further confirm the predictions made from our previous continuum hydrodynamic modeling and simulations, which are actually semi-quantitatively accurate down to the small length scales in the problem. Using MD simulations, we investigate the motion of evaporative droplets in one-component Lennard-Jones fluids confined in nanochannels with thermal gradients. The droplet is found to migrate in the direction of decreasing temperature of solid walls, with a migration velocity linearly proportional to the temperature gradient. This agrees with the prediction of our continuum model. We then measure the effect of droplet size on the droplet motion. It is found that the droplet mobility is inversely proportional to a dimensionless coefficient associated with the total rate of dissipation due to droplet movement. Our results show that this coefficient is of order unity and increases with the droplet size for the small droplets (~10 nm) simulated in the present work. These findings are in semi-quantitative agreement with the predictions of our continuum model. Finally, we measure the effect of liquid-vapor coexistence temperature on the droplet motion. Through a theoretical analysis on the size of the thermal singularity, it can be shown that the droplet mobility decreases with decreasing coexistence temperature. This is observed in our MD simulations.

From concave to convex: capillary bridges in slit pore geometry

We investigate the morphological evolution of nonaxisymmetric capillary bridges in slit-pore geometry as the height of the pore and aspect ratio of the bridge are varied. The liquid bridges are formed between two hydrophobic surfaces patterned with hydrophilic strips. The aspect ratio of the capillary bridges (length/width) is varied from 2.5 to 120 by changing the separation between the surfaces, the width of the strips, or the fluid volume. As the bridge height is increased, the aspect ratio decreases and we observe a large increase in the mean curvature of the bridge. More specifically, the following counterintuitive result is observed: the mean curvature of the bridges changes sign and goes from negative (concave bridge) to positive (convex bridge) when the height is increased at constant volume. These experimental observations are in quantitative agreement with Surface Evolver simulations. Scaling shows a collapse of the data indicating that this transition in the sign of the Laplace pressure is universal for capillary bridges with high aspect ratios. Finally, we show that the morphology diagrams obtained from our 3D analysis are considerably different from those expected from a 2D analysis.

Diffusive domain coarsening: early time dynamics and finite-size effects

We study the diffusive dynamics of phase separation in a symmetric binary (A + B) mixture with a 50:50 composition of A and B particles, following a quench below the demixing critical temperature, both in spatial dimensions d=2 and d=3. The particular focus of this work is to obtain information about the effects of system size and correction to the growth law via the appropriate application of the finite-size scaling method to the results obtained from the Kawasaki exchange Monte Carlo simulation of the Ising model. Observations of only weak size effects and a very small correction to scaling in the growth law are significant. The methods used in this work and information thus gathered will be useful in the study of the kinetics of phase separation in fluids and other problems of growing length scale. We also provide a detailed discussion of the standard methods of understanding simulation results which may lead to inappropriate conclusions.

Contact angle dependence on the fluid-wall dispersive energy

Menisci of the truncated and shifted Lennard-Jones fluid between parallel planar walls are investigated by molecular dynamics simulation. Thereby, the characteristic energy of the unlike dispersive interaction between fluid molecules and wall atoms is systematically varied to determine its influence on the contact angle. The temperature is varied as well, covering most of the range between the triple-point temperature and the critical temperature of the bulk fluid. The transition between obtuse and acute angles is found to occur at a temperature-independent magnitude of the fluid-wall dispersive interaction energy. On the basis of the present simulation results, fluid-wall interaction potentials can be adjusted to contact angle measurements.