PNNL

Abstract

The structure of thin water films during the rupture process was investigated by a new approach, which combines molecular dynamics simulation (MDS) with image processing analysis. The analysis procedure was developed to convert MDS trajectories to readable 3D images. The water films were studied at different thicknesses by MDS to determine the critical thickness at which the film ruptures. The potential energy of each specific film thickness during the simulation time was analyzed, and the results showed that the potential energy of stable films remained unchanged while the potential energy kept decreasing for films which ruptured during the simulation time. By applying the new procedure, the molecular porosity, which is defined as the void fraction between the volume of molecular pores in the water film and the total volume of the water film, was calculated. The results of molecular porosity for different film thicknesses during the simulation time suggested a critical molecular porosity as 49%. In other words, stable films have a molecular porosity of less than 49%. If a water film has a molecular porosity greater than 49%, rupture occurs during the simulation. This research was supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences.