• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2005년도 창립60주년 기념 추계 학술대회
• /
• pp.56.2-56.2
• /
• 2005

• Bulletin of the Korean Chemical Society
• /
• 제25권2호
• /
• pp.321-324
• /
• 2004

• Bulletin of the Korean Chemical Society
• /
• 제25권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
• /
• 제28권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
• /
• 제12권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
• /
• 제12권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.

• 대한기계학회논문집B
• /
• 제29권1호
• /
• 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
• /
• 제29권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
• /
• 제30권6호
• /
• pp.371-378
• /
• 2017

본 연구에서는 대표적인 열경화성 재료인 에폭시 기지의 흡습탄성 거동을 예측하기 위해 분자동역학 전산모사를수행하였다. 고분자 복합재가 오랜 시간 동안 흡습환경에 지속적으로 노출될 경우, 거시적 물성의 저하가 발생하기 때문에 복합재의 내구설계 측면에 있어 흡습노화 현상에 대해 분자스케일적으로 접근하는 방법은 매우 중요하다. 본 연구에서는 $EPON862^{(R)}$ 수지와 아민계 Triethylenetetramine (TETA) 경화제로 비정질 에폭시 분자모델을 구성하였으며, 각각 30과 90%의 가교 상태에서 수분 흡수 유무에 따른 물성변화를 관찰하였다. 건조상태의 에폭시와 수분이 4 wt% 포함된 에폭시 단위셀에 대한 평형 및 비평형 앙상블 전산모사 과정을 통해, 에폭시의 수분팽창계수, 응력-변형률 선도 및 탄성계수 그리고 침투된 수분의 수지 내 확산계수를 예측하였다. 또한 흡습된 구조와 그에 따른 물성변화의 상관관계를 규명하기 위해, 자유체적 변화 및 흡습에 따른 에폭시 수지의 비결합 포텐셜 에너지 변화를 관찰하였다.

• EDISON SW 활용 경진대회 논문집
• /
• 제5회(2016년)
• /
• 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.