• Proceedings of the KSME Conference
• /
• 2005.11a
• /
• pp.56.2-56.2
• /
• 2005

• Bulletin of the Korean Chemical Society
• /
• v.25 no.2
• /
• pp.321-324
• /
• 2004

• Bulletin of the Korean Chemical Society
• /
• v.25 no.5
• /
• pp.737-741
• /
• 2004

We presents new results for transport properties of dumbbell fluids by equilibrium molecular dynamics (EMD) simulations using Green-Kubo and Einstein formulas. It is evident that the interaction between dumbbell molecules is less attractive than that between spherical molecules which leads to higher diffusion and to lower friction. The calculated viscosity, however, is almost independent on the molecular elongation within statistical error bar, which is contradicted to the Stokes' law. The calculated thermal conductivity increases and then decreases as molecular elongation increases. These results of viscosity and thermal conductivity for dumbbell molecules by EMD simulations are inconsistent with the earlier results of those by non-equilibrium molecular dynamics (NEMD) simulations. The possible limitation of the Green-Kubo and Einstein formulas with regard to the calculations of viscosity and thermal conductivity for molecular fluids such as the missing rotational degree of freedom is pointed out.

• Bulletin of the Korean Chemical Society
• /
• v.28 no.10
• /
• pp.1697-1704
• /
• 2007

We present results for transport properties of diatomic fluids by isothermal-isobaric (NpT) equilibrium molecular dynamics (EMD) simulations using Green-Kubo and Einstein formulas. As the molecular elongation of diatomic molecules increases from the spherical monatomic molecule, the diffusion coefficient increases, indicating that longish shape molecules diffuse more than spherical molecules, and the rotational diffusion coefficients are almost the same in the statistical error since random rotation decreases. The calculated translational viscosity decreases with the molecular elongation of diatomic molecule within statistical error bar, while the rotational viscosity increases. The total thermal conductivity decreases as the molecular elongation increases. This result of thermal conductivity for diatomic molecules by EMD simulations is again inconsistent with the earlier results of those by non-equilibrium molecular dynamics (NEMD) simulations even though the missing terms related to rotational degree of freedom into the Green-Kubo and Einstein formulas with regard to the calculation of thermal conductivity for molecular fluids are included.

• Bulletin of the Korean Chemical Society
• /
• v.12 no.3
• /
• pp.315-322
• /
• 1991

In a recent $paper^1$ we reported equilibrium (EMD) and non-equilibrium (NEMD) molecular dynamics simulations of liquid argon using the Green-Kubo relations and NEMD algorithms to calculate the thermal transport coefficients-the self-diffusion coefficient, shear viscosity, and thermal conductivity. The overall agreement with experimental data is quite good. In this paper the same technique is applied to calculate the thermal transport coefficients of liquid water at 298.15 K and 1 atm using TIP4P model for the interaction between water molecules. The EMD results show difficulty to apply the Green-Kubo relations since the time-correlation functions of liquid water are oscillating and not decaying rapidly enough except the velocity auto-correlation function. The NEMD results are found to be within approximately ${\pm}$30-40% error bars, which makes it possible to apply the NEMD technique to other molecular liquids.

• Bulletin of the Korean Chemical Society
• /
• v.12 no.3
• /
• pp.309-315
• /
• 1991

The thermal transport coefficients-the self-diffusion coefficient, shear viscosity, and thermal conductivity-of liquid argon at 94.4 K and 1 atm are calculated by non-equilibrium molecular dynamics (NEMD) simulations of a Lennard-Jones potential and compared with those obtained from Green-Kubo relations using equilibrium molecular dynamics (EMD) simulations and with experimental data. The time-correlation functions-the velocity, pressure, and heat flux auto-correlation functions-of liquid argon obtained from the EMD simulations show well-behaved smooth curves which are not oscillating and decaying fast around 1.5 ps. The calculated self-diffusion coefficient from our NEMD simulation is found to be approximately 40% higher than the experimental result. The Lagrange extrapolated shear viscosity is in good agreement with the experimental result and the asymptotic formula of the calculated shear viscosities seems to be an exponential form rather than the square-root form predicted by other NEMD studies of shear viscosity. The agreement for thermal conductivity between the simulation results (NEMD and EMD) and the experimental result is within statistical error. In conclusion, through our NEMD and EMD simulations, the overall agreement is quite good, which means that the Green-Kubo relations and the NEMD algorithms of thermal transport coefficients for simple liquids are valid.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.1 s.232
• /
• pp.159-168
• /
• 2005

In recent studies, it was reported that there existed the temperature discontinuity at a liquid-vapor interface in an equilibrium state. However, from the viewpoint of the classical thermodynamics, it is highly questionable result although considering that the experiments related with a boundary layer is very difficult due to the extremely thin thickness of it. To clarify whether the temperature discontinuity over a liquid-vapor interface really exists, the computer simulations were performed. From the simulation results, it could be concluded that the misconception in a temperature calculation might result in non-uniform temperature distributions over an interface under an equilibrium state.

• Journal of Advanced Marine Engineering and Technology
• /
• v.29 no.1
• /
• pp.60-67
• /
• 2005

Generally. in the molecular dynamics simulations. the Verlet neighboring list algorithm is used for the reduction of a simulation time On the other hand. the application of the Verlet neighboring list forces the time evolution of a simulation system to follow an unrealistic path in a phase space. In equilibrium state, it does not matter with the simulation results because the individual molecule's motion is originally random and any effect due to a small deviation from a real time evolution can be completely ignored. However, if an unsteady state is involved. such a deviation may significantly affect to the results. That is, there is a Possibility that the simulation results Provide ones with any misleading data In this study we evaluated the effect due to the Verlet neighboring list in performing the simulation of a non-equilibrium state and suggested the method to avoid it.

• Composites Research
• /
• v.30 no.6
• /
• pp.371-378
• /
• 2017

In this study, hygroelastic behavior of thermosetting epoxy is predicted by molecular dynamics simulations. Since consistent exposures to humid environments lead to macroscopic degradation of polymer composite, computational simulation study of the hygroscopically aged epoxy cell is essential for long-time durability. Therefore, we modeled amorphous epoxy molecular unit cell structures at a crosslinking ratio of 30, 90% and with the moisture weight fraction of 0, 4 wt% respectively. Diglycidyl ether of bisphenol F (EPON862) and triethylenetetramine (TETA) are chosen as resin and curing agent respectively. Incorporating equilibrium and non-equilibrium ensemble simulation with a classical interatomic potential, various hygroelastic properties including diffusion coefficient of water, coefficient of moisture expansion (CME), stress-strain curve and elastic modulus are predicted. To establish the structural property relationship of pure epoxy, free volume and internal non-bond potential energy of epoxy are examined.

• Proceeding of EDISON Challenge
• /
• 2016.03a
• /
• pp.37-54
• /
• 2016

The methods for estimating the change of free energy in dissipative particle dynamics (DPD) are discussed on the basis of fluctuation theorems. Fluctuation theorems are tactics to evaluate free energy changes from non-equilibrium work distributions and have several forms, as proposed by Jarzynski, Crooks, and Bennett. The validity of these methods however, has been shown merely with the molecular dynamics or Langevin dynamics. In this study, the appropriate forms of fluctuation theorems for dissipative particle dynamics, which has similar structure to that of Langevin dynamics, are suggested using Liouville's theorem, and they are proved equivalent to original fluctuation theorems. Work distribution functions, which are probability distribution functions of works exerted on the system within the systematic change, are the basics of fluctuation theorems and their shapes are turned out to be dependent on the phase space trajectory of the change of the system. The reliability of Jarzynski and Crooks methods is highly dependent on the number of simulations to measure works and the shapes of the work distribution functions. Bennett method, however, can evaluate free energy changes even when Jarzynski and Crooks methods fail to do so.