Search

Steady state pressure driven flow of liquid argon through a finite length cylindrical nanopore was investigated numerically by classical Navier-Stokes (NS) hydrodynamic models and nonequilibrium molecular dynamics (MD) simulations. In both approaches, the nanopore was nominally 2.2 nm in diameter and 6 nm long. For the MD simulations, the intermolecular properties of the walls were specified independently from the liquid. Comparisons between the approaches were made in terms of the gross feature of total flow rate through the nanopore, as well as the more refined considerations of the spatial distributions of pressure, density, and velocity. The results showed that for the NS equations to predict the same trends in total flow rate with increasing pressure difference as the MD simulation, submodels for variations in density and viscosity with pressure are needed to be included. The classical NS boundary conditions quantitatively agreed with the flow rate predictions from MD simulations only under the condition of having a neutral-like solid-liquid interaction. Under these conditions, the NS and MD models also agreed well in streamwise distributions of pressure, density, and velocity, but not in the radial direction. [ABSTRACT FROM AUTHOR]

Fluid transport through a nanopore in a membrane was investigated by using a novel molecular dynamics approach proposed in this study. The advantages of this method, relative to dual-control-volume grand-canonical molecular dynamics method, are that it eliminates disruptions to the system dynamics that are normally created by inserting or deleting particles from control volumes, and that it functions well for dense systems due to the number of particles being fixed in the system. Using the proposed method, we examined liquid argon transport through a nanopore by performing nonequilibrium molecular dynamics (NEMD) simulations under different back pressures. Validation of the code was performed by comparing simulation results to published experimental data obtained under equilibrium conditions. NEMD results show that constant pressure difference across the membrane was readily achieved. [ABSTRACT FROM AUTHOR]

For a model system of polyethylene of chain lengths 40 and 100 carbon atoms, we calculated the pressure at different densities and compared them with the experimental values. The simulation was conducted on the second nearest neighbor diamond lattice, and the pressure was calculated using the virtual-volume-variation method after the system was reverse mapped to its fully atomistic form in continuous space and energy minimized. In addition, the pressure was also calculated from the virial route by conducting a short molecular dynamics simulation starting from the energy minimized structure. We show that the pressure obtained from our simulations is quite reasonable in the length of simulation time (in Monte Carlo steps) normally employed in our group. These results provide additional evidence for the equilibration of our model systems, and methodology to calculate the pressure in our lattice models. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

An amorphous 50/50 atactic polypropylene (aPP)/isotactic polypropylene (iPP) mixture at 125 °C was simulated using a second nearest neighbor diamond lattice and a three states rotational isometric state model. The result suggests that at the liquidlike density that corresponds to the atmospheric pressure, aPP prefers to interact with other aPP chains rather than with iPP chains. The result is consistent with the inference of Keith and Padden [J. Appl. Phys. 35, 1286 (1964)] that aPP and iPP will tend to separate from one another in their melt at 125 °C, before the onset of crystallization of iPP. The tendency for immiscibility of the amorphous aPP/iPP blend is likely attributed to the presence of short syndiotactic sequences in the aPP chains adopting all-trans conformations. The attractive intermolecular interaction of pairs of such subchains at 125 °C promotes the separation of aPP from iPP. This interaction is weakened at higher temperature, where aPP and iPP become miscible. The result also shows that miscibility of the blend increases with increasing pressure. However, the origin of the pressure effect is not clear. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations using a modified Drieding 2.21 force field were carried out to study the coalescence behavior of nanometer-sized water droplets in vacuum and in n-heptane. The coalescence mechanisms of the water droplets in the above-noted environments are fairly similar in a sense that the water droplets form a bridge linking the droplets before they merge. However, in the latter situation, due to the presence of n-heptane molecules in between the water droplets, the coalescence was observed to be slowed down considerably, especially in the first 10 ps of the process. However, once the bridge is formed, the water droplets, in both situations, spend about the same amount of time to form a single droplet. The maximum distance between the droplets above which coalescence does not occur was found to be 10 Å. In terms of the dynamics, the diffusion coefficient of n-heptane in the emulsion system was very close to its value in the pure liquid form. This may be because n-heptane is the continuous phase. Nonetheless, the dynamic behavior of water in n-heptane is different from that of pure water during and after the coalescence. In particular, the self-diffusion coefficient of water molecules in n-heptane is about 20% higher than the experimental value of pure water. Due to the lack of strong attraction forces between water and n-heptane molecules, the n-heptane molecules were observed to orient themselves perpendicularly to the water/n-heptane interfaces so that the contacting area is minimized. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations were carried out on single chain models of linear low-density polyethylene in vacuum to study the effects of branch length, branch content, and branch distribution on the polymer’s crystalline structure at 300 K. The trans/gauche (t/g) ratios of the backbones of the modeled molecules were calculated and utilized to characterize their degree of crystallinity. The results show that the t/g ratio decreases with increasing branch content regardless of branch length and branch distribution, indicating that branch content is the key molecular parameter that controls the degree of crystallinity. Although t/g ratios of the models with the same branch content vary, they are of secondary importance. However, our data suggests that branch distribution (regular or random) has a significant effect on the degree of crystallinity for models containing 10 hexyl branches/1,000 backbone carbons. The fractions of branches that resided in the equilibrium crystalline structures of the models were also calculated. On average, 9.8% and 2.5% of the branches were found in the crystallites of the molecules with ethyl and hexyl branches while C[sup 13] NMR experiments showed that the respective probabilities of branch inclusion for ethyl and hexyl branches are 10% and 6% [Hosoda et al., Polymer 1990, 31, 1999–2005]. However, the degree of branch inclusion seems to be insensitive to the branch content and branch distribution. © 2001 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Steady state pressure driven flow of liquid argon through a finite length cylindrical nanopore was investigated numerically by classical Navier-Stokes (NS) hydrodynamic models and nonequilibrium molecular dynamics (MD) simulations. In both approaches, the nanopore was nominally 2.2 nm in diameter and 6 nm long. For the MD simulations, the intermolecular properties of the walls were specified independently from the liquid. Comparisons between the approaches were made in terms of the gross feature of total flow rate through the nanopore, as well as the more refined considerations of the spatial distributions of pressure, density, and velocity. The results showed that for the NS equations to predict the same trends in total flow rate with increasing pressure difference as the MD simulation, submodels for variations in density and viscosity with pressure are needed to be included. The classical NS boundary conditions quantitatively agreed with the flow rate predictions from MD simulations only under the condition of having a neutral-like solid-liquid interaction. Under these conditions, the NS and MD models also agreed well in streamwise distributions of pressure, density, and velocity, but not in the radial direction.

Fluid transport through a nanopore in a membrane was investigated by using a novel molecular dynamics approach proposed in this study. The advantages of this method, relative to dual-control-volume grand-canonical molecular dynamics method, are that it eliminates disruptions to the system dynamics that are normally created by inserting or deleting particles from control volumes, and that it functions well for dense systems due to the number of particles being fixed in the system. Using the proposed method, we examined liquid argon transport through a nanopore by performing nonequilibrium molecular dynamics (NEMD) simulations under different back pressures. Validation of the code was performed by comparing simulation results to published experimental data obtained under equilibrium conditions. NEMD results show that constant pressure difference across the membrane was readily achieved.