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

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1584-1589
• /
• 2003

This work applies the equilibrium molecular dynamics simulation method to study a Lennard-Jones liquid thin film suspended in the vapor and calculates diffusion coefficients by Green-Kubo equation derived from Einstein relationship. As a preliminary test, the diffusion coefficients of the pure argon fluid are calculated by equilibrium molecular dynamics simulation. It is found that the diffusion coefficients increase with decreasing the density and increasing the temperature. When both argon liquid and vapor phases are present, the effects of the system temperature on the diffusion coefficient are investigated. It can be seen that the diffusion coefficient significantly increases with the temperature of the system.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.19 no.6
• /
• pp.54-63
• /
• 2015

We present a microscopic understanding of the chemical erosion due to combustion product on the nozzle throat using molecular dynamics simulations. The present erosion process consists of molecule-addition step and equilibrium step. First, either $CO_2$ or $H_2O$ are introduced into the system with high velocity to provoke the collision with graphite surface. Then, the equilibrium simulation is followed. The collision-included dissociation and its influence on the erosion is emphasized and the present molecular observations are compared with the macroscopic chemical reaction model.

• Journal of the Korean Institute of Telematics and Electronics D
• /
• v.35D no.12
• /
• pp.30-40
• /
• 1998

We have calculated ultra-low energy silicon-self ion implantations and silicon damages through classical molecular dynamics simulation using empirical potentials. We tested whether the recently developed Environment-Dependent Interatomic Potential(EDIP) was suitable for ultra low energy ion implantation simulation, and found that point defects formation energies were in good agreement with other theoretical calculations, but the calculated vacancy migration energy was overestimated. Most of the damages that are produced by collision cascades are concentrated into amorphous-like pockets. Also, We upgraded MDRANGE code for silicon ion implantation process simulation. We simulated ultra-low energy boron ion implantation, 200eV, 500eV, and 1000eV respectively, and calculated boron profiles with silicon substrate temperature and tilt angle. We investigated that below 1000eV, channeling effect must be considered.

• 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.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2009.04a
• /
• pp.34-37
• /
• 2009

본 논문에서는 원자적 계산(atomistic calculation)을 이용한 나노박막의 평형상태(self-equilibrium state)에 대한 해석기법을 제시한다. 두께가 얇은 나노박막은 표면 응력(surface stress)에 의한 영향으로 원자간 거리가 벌크상태의 거리보다 작아진다. 두께가 얇은 나노박막에서의 원자 사이의 거리는 표면 응력과 탄성계수들의 표현식으로 계산이 가능하며, 본 논문에서는 {100}, {111}, {110} 표면을 가지는 나노박막의 평형상태의 해석을 위한 해석적 방법을 제시한다. 원자 사이의 거리를 계산하기 위해서는 보다 정확한 표면 응력의 계산방법이 필요하다. 본 연구에서는 나노박막의 평형상태에 대한 해석을 위해 surface relaxation model을 제시하고, 이 모델을 이용하여 표면응력(surface stress)과 표면강성계수(surface stiffness tensor)와 같은 surface parameter의 계산을 수행한다. 본 논문에서 제시된 surface relaxation model을 검증하기 위하여 분자동역학 전산모사(molecular dynamics simulation)의 수치 결과를 제시하고, 본 연구에서 계산한 equilibrium strain과 비교 검증한다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.4
• /
• pp.301-305
• /
• 2012

In this work the mechanical behaviors of GaN nanowires are analyzed during tension, compression, and unloading deformations. The thermal conductivity of the nanowires at each deformed state is evaluated using an equilibrium Green-Kubo approach. Under tensile loading, the [0001]-oriented nanowires with hexagonal cross-sections undergo a phase transformation from wurtzite to a tetragonal structure. The phase transformation is not observed under compressive loading. The thermal conductivity decreases on going from compressive strains to tensile strains. The strain dependence of the thermal conductivity results from the relaxation time of phonon. A reverse transformation from the tetragonal structure to the wurtzite structure is observed during unloading. The thermal conductivities in the intermediate states are lower than the conductivity in the wurtzite structure at same strain. Such differences in the thermal conductivity between different atomic structures are mainly due to changes in the group velocity of phonon.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.4
• /
• pp.331-334
• /
• 2009

In this paper, we present a surface relaxation model in atomistic calculations for thin nanofilms. This surface relaxation model is very simple model which have only two degrees of freedoms to determine the atomic positions of nanofilms. Whereas in conventional molecular statics simulations, the same number of degrees of freedoms at all atom positions are used as unknown variables. In order to prove the reliability of the presented model, we present the results of self-equilibrium strain calculations with the surface parameters obtained from this model.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.7
• /
• pp.60-68
• /
• 2004

High velocity impacts are accompanied with large deformations, which generate a large amount of heat due to plastic works, resulting in a significant temperature rise of the material. Because the elevated temperature affects the dynamic properties of materials, it is important to predict the temperature rise during high-stram-rate deformations. Both existing vacancies and excess vacancies are credited to the stored energy, yet it is difficult to distinguish one from another in contribution to the stored energy using macroscopic level materials models. In this study, an atomistic material model for fee materials such as copper is set up to calculate the stored energy using molecular dynamics (MD) simulations. It is concluded that excess vacancies play an important role for the stored energy during a high-strain-rate deformation.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.1086-1091
• /
• 2004

The boundary layer is a very important characteristic of a liquid-vapor interface since it governs the heat and mass transfer phenomena across an interface. However, the thickness of a boundary layer is generally micro- or nano-sized, which requires highly accurate measurement devices and, consequently, costs the related experiments very high and time-consuming. Due to these size dependent limitations, the experiments related with a nano-scaled size have suffered from the errors and the reliability of the obtained data. This study is performed to grasp the characteristics of a liquid-vapor interface, by using a molecular dynamics method. The simulation results were compared with other studies if possible. Although other studies reported that there existed a temperature discontinuity over an interface when the system was reduced to micro- or nano-sized, we confirmed that there was no such a temperature discontinuity.