Molecular Dynamics Simulation of Nanodroplet Evaporation

2000 ◽  
Vol 123 (4) ◽  
pp. 741-748 ◽  
Author(s):  
J. H. Walther ◽  
P. Koumoutsakos

Molecular dynamics simulations are used to study the sub-critical evaporation of a nanometer-size droplet at 300 K and 3 MPa. Classical molecular dynamics techniques are combined with an adaptive tree data structure for the construction of the neighbor lists, allowing efficient simulations using hundreds of thousands of molecules. We present a systematic convergence study of the method demonstrating its convergence for heat conduction problems in submicron scales. These high resolution simulations compute values of the evaporation coefficient that are in excellent agreement with theoretical predictions.

2001 ◽  
Vol 703 ◽  
Author(s):  
Sebastian von Alfthan ◽  
Antti Kuronen ◽  
Kimmo Kaski

ABSTRACTEffect of a crystalline-amorphous interface on heat conduction has been studied using atom-istic simulations of a silicon system. System with amorphous silicon was created using the bond-switching Monte Carlo simulation method and heat conduction near room temperature was studied by molecular dynamics simulations of this system.


Author(s):  
Koji Miyazaki ◽  
Yoshizumi Iida ◽  
Daisuke Nagai ◽  
Hiroshi Tsukamoto

We carried out molecular dynamics simulations (MD) of heat conduction in Si thin film and Si films with a nano-hole to represent the nano-structure, in order to investigate the mechanism of the thermal conductivity reduction of nano-structured materials. The Stillinger-Weber potential is used in this study. Different temperatures are applied at the both sides of boundaries of the calculation domain in the z-direction, and periodic boundary conditions are applied in the x and y directions. The calculated temperature profile of a Si thin film of 10.86nm thickness is compared to that calculated by using the phonon Boltzmann transport equation (BTE). These agreed reasonably well with each other, and the phonon mean free path of Si is estimated to be several tens of nanometers. Molecular dynamics simulation of Si at the uniform temperature of 800K is also carried out. Phonon dispersion curves are calculated by using the time-space 2D Fourier transform. The phonon modes at high frequency are not present in nano-structures of Si. We discuss the mechanism of the reduction of the thermal conductivity of nano-structured material on the atomic scale.


Author(s):  
Koji Miyazaki ◽  
Daisuke Nagai ◽  
Hiroshi Tsukamoto

We carried out molecular dynamics simulations (MD) of heat conduction in Si thin film and Si films with a nano-hole to represent the nano-structure, in order to investigate the mechanism of the thermal conductivity reduction of nano-structured materials. The Stillinger-Weber potential is used in this study. Different temperatures are applied at the both sides of boundaries of the calculation domain in the z-direction, and periodic boundary conditions are applied in the x and y directions. The calculated temperature profile of a Si thin film of 10.86nm thickness is compared to that calculated by using the phonon Boltzmann transport equation (BTE). These agreed reasonably well with each other, and the phonon mean free path of Si is estimated to be several tens of nanometers. Molecular dynamics simulation of Si at the uniform temperature of 800K is also carried out. Phonon dispersion curves are calculated by using the time-space 2D Fourier transform. The phonon modes at high frequency are not present in nano-structures of Si. We discuss the mechanism of the reduction of the thermal conductivity of nano-structured material on the atomic scale.


Author(s):  
Toshihiro Kaneko ◽  
Kenji Yasuoka ◽  
Ayori Mitsutake ◽  
Xiao Cheng Zeng

Multicanonical molecular dynamics simulations are applied, for the first time, to study the liquid-solid and solid-solid transitions in Lennard-Jones (LJ) clusters. The transition temperatures are estimated based on the peak position in the heat capacity versus temperature curve. For LJ31, LJ58 and LJ98, our results on the solid-solid transition temperature are in good agreement with previous ones. For LJ309, the predicted liquid-solid transition temperature is also in agreement with previous result.


RSC Advances ◽  
2020 ◽  
Vol 10 (9) ◽  
pp. 5507-5515
Author(s):  
Liang Song ◽  
Feng-Qi Zhao ◽  
Si-Yu Xu ◽  
Xue-Hai Ju

The bimolecular and fused ring compounds are found in the high-temperature pyrolysis of NONA using ReaxFF molecular dynamics simulations.


2017 ◽  
Vol 890 ◽  
pp. 252-259
Author(s):  
Le Wang ◽  
Guan Cheng Jiang ◽  
Xin Lin ◽  
Xian Min Zhang ◽  
Qi Hui Jiang

Molecular dynamics simulations are used to study the dissociation inhibiting mechanism of lecithin for structure I hydrates. Adsorption characteristics of lecithin and PVP (poly (N-vinylpyrrolidine)) on the hydrate surfaces were performed in the NVT ensemble at temperatures of 277K and the hydrate dissociation process were simulated in the NPT ensemble at same temperature. The results show that hydrate surfaces with lecithin is more stable than the ones with PVP for the lower potential energy. The conformation of lecithin changes constantly after the balanced state is reached while the PVP molecular dose not. Lecithin molecule has interaction with lecithin nearby and hydrocarbon-chains of lecithin molecules will form a network to prevent the diffusion of water and methane molecules, which will narrow the available space for hydrate methane and water movement. Compared with PVP-hydrate simulation, analysis results (snapshots and mass density profile) of the dissociation simulations show that lecithin-hydrate dissociates more slowly.


1994 ◽  
Vol 373 ◽  
Author(s):  
Roger E. Stoller

AbstractA series of high-energy, up to 20 keV, displacement cascades in iron have been investigated for times up to 200 ps at 100 K using the method of molecular dynamics simulation. Thesimulations were carried out using the MOLDY code and a modified version of the many-bodyinteratomic potential developed by Finnis and Sinclair. The paper focuses on those results obtained at the highest energies, 10 and 20 keV. The results indicate that the fraction of the Frenkel pairs surviving in-cascade recombination remains fairly high in iron and that the fraction of the surviving point defects that cluster is lower than in materials such as copper. In particular, vacancy clustering appears to be inhibited in iron. Some of the interstitial clusters were observed to exhibit an unexpectedly complex, three-dimensional morphology. The observations are discussed in terms of their relevance to microstructural evolution and mechanical property changes in irradiated iron-based alloys.


2021 ◽  
Author(s):  
Kazushi Fujimoto ◽  
Tetsuro Nagai ◽  
Tsuyoshi Yamaguchi

<div>The position-dependent diffusion coefficient along with free energy profile are important parameters needed to study mass transport in heterogeneous systems such as biological and polymer membranes, and molecular dynamics (MD) calculation is a popular tool to obtain them. Among many methodologies, the Marrink-Berendsen (MB) method is often employed to calculate the position-dependent diffusion coefficient, in which the autocorrelation function of the force on a fixed molecule is related to the friction on the molecule. However, the diffusion coefficient is shown to be affected by the period of the removal of the center-of-mass velocity, which is necessary when performing MD calculations using the Ewald method for Coulombic interaction. We have clarified theoretically in this study how this operation affects the diffusion coefficient calculated by the MB method, and the theoretical predictions are proven by MD calculations. Therefore, we succeeded in providing guidance on how to select an appropriate the period of the removal of the center-of-mass velocity in estimating the position-dependent diffusion coefficient by the MB method. This guideline is applicable also to the Woolf-Roux method.</div>


2014 ◽  
Vol 10 (1) ◽  
pp. 42-56 ◽  
Author(s):  
Zailani Abdullah ◽  
Tutut Herawan ◽  
A. Noraziah ◽  
Mustafa Mat Deris

Frequent Pattern Tree (FP-Tree) is a compact data structure of representing frequent itemsets. The construction of FP-Tree is very important prior to frequent patterns mining. However, there have been too limited efforts specifically focused on constructing FP-Tree data structure beyond from its original database. In typical FP-Tree construction, besides the prior knowledge on support threshold, it also requires two database scans; first to build and sort the frequent patterns and second to build its prefix paths. Thus, twice database scanning is a key and major limitation in completing the construction of FP-Tree. Therefore, this paper suggests scalable Trie Transformation Technique Algorithm (T3A) to convert our predefined tree data structure, Disorder Support Trie Itemset (DOSTrieIT) into FP-Tree. Experiment results through two UCI benchmark datasets show that the proposed T3A generates FP-Tree up to 3 magnitudes faster than that the benchmarked FP-Growth.


Sign in / Sign up

Export Citation Format

Share Document