• 설비공학논문집
• /
• 제22권2호
• /
• pp.64-69
• /
• 2010

We investigated the variation in contact angle of a nano-sized water droplet on a nano stripe-patterned surface using molecular dynamics simulation. By changing the height and width of the stripe pillar, and the gap width of the stripes, we observed the contact angle of water droplet in equilibrium. When the surface energies were 0.1 and 0.3 kcal/mol, the calculated contact angles were in good agreement with the Cassie and Baxter equation. However, when the surface energy is 0.5 kcal/mol, the contact angles are observed to be perturbed along the Cassie and Baxter equation.

• Bulletin of the Korean Chemical Society
• /
• 제28권10호
• /
• pp.1689-1696
• /
• 2007

Equilibrium molecular dynamics (EMD) simulations are used to evaluate the transport coefficients of argonkrypton mixtures at two liquid states (state A: 94.4 K and 1 atm; state B: 135 K and 39.5 atm) via modified Green-Kubo formulas. The composition dependency of the volume at state A obeys close to the linear model for ideal liquid mixture, while that at state B differs from the linear model probably due to the high pressure. The radial distribution functions for the Ar-Kr mixture (x = 2/3) show a mixing effect: the first peak of g11 is higher than that of g(r) for pure Ar and the first peak of g22 is lower than that of g(r) for pure Kr. An exponential model of engineering correlation for diffusion coefficient (D) and shear viscosity (η) is superior to the simple linear model for ideal liquid mixtures. All three components of thermal conductivity (λpm, λtm, and λti) at state A and hence the total thermal conductivity decrease with the increase of x. At state B, the change in λtm is dominant over those in λpm and λti, and hence the total thermal conductivity decrease with the increase of x.

• 한국전산구조공학회논문집
• /
• 제22권4호
• /
• pp.331-334
• /
• 2009

본 논문에서는 원자적 계산(atomistic calculation)을 위한 해석적 모델로 surface relaxation model을 제시한다. 기존의 분자정역학(molecular statics)이 모든 원자의 위치를 자유도로 선정하여 사용하는데 반하여, 이 모델은 면내방향에 해당하는 두 개의 자유도로 나노박막의 원자 위치를 기술하는 매우 간단한 방법이다. 본 연구에서는 surface relaxation model을 이용하여 표면응력(surface stress)과 표면강성계수(surface stiffness tensor)와 같은 표면인자(surface parameter)의 계산을 수행하고, surface stress model을 이용하여 평형상태에서의 원자의 위치정보를 계산한다. 그리고 surface relaxation model을 검증하기 위하여 분자동역학 전산모사(molecular dynamics simulation)의 수치 결과가 제시되며, 본 연구에서 계산한 equilibrium strain과 비교 검증한다.

• Macromolecular Research
• /
• 제10권6호
• /
• pp.297-303
• /
• 2002

We have investigated the microstructural properties of a weakly charged polyelectrolyte modeled with both Hookean spring and Debye-Huckel potential, by employing a novel hybrid scheme of molecular dynamics (MD) and Monte Carlo (MC) simulations. Although the off-lattice pivot step facilitates the earlier computations stage, it gives rise to oscillations and hinders the stable equilibrium state. In order to overcome this problem, we adopt the MC off-lattice pivot step in early stage only, and then switch the computation to a pure MD step. The result shows that the computational speed-up compared to the previous method is entirely above 10 to 50, without loss of the accuracy. We examined the conformations of polyelectrolyte in solvents in terms of the end-to-end distance, radius of gyration, and structure factor with variations of the screening effects of solvent and the monomer charges. The emphasis can favorably be given on the elongation behavior of a polyelectrolyte chain, with observing the simultaneous snapshots.

• Bulletin of the Korean Chemical Society
• /
• 제29권8호
• /
• pp.1554-1560
• /
• 2008

The probe diffusion and friction constants of methyl yellow (MY) in liquid n-alkanes of increasing chain length were calculated by equilibrium molecular dynamics (MD) simulations at temperatures of 318, 418, 518 and 618 K. Lennard-Jones particles with masses of 225 and 114 g/mol are modeled for MY. We observed that the diffusion constant of the probe molecule follows a power law dependence on the molecular weight of nalkanes, DMY${\sim}M^{-\gamma}$ well. As the molecular weight of n-alkanes increases, the exponent $\gamma$ shows sharp transitions near n-dotriacontane ($C_{32}$) for the large probe molecule (MY2) at low temperatures of 318 and 418 K. For the small probe molecule (MY1) $D_{MY1}$ in $C_{12}$ to C80 at all the temperatures are always larger than Dself of n-alkanes and longer chain n-alkanes offer a reduced friction relative to the shorter chain n-alkanes, but this reduction in the microscopic friction for MY1 is not large enough to cause a transition in the power law exponent in the log-log plot of DMY1 vs M of n-alkane. For the large probe molecule (MY2) at high temperatures, the situation is very similar to that for MY1. At low temperatures and at low molecular weights of n-alkanes, $D_{MY2}$ are smaller than $D_{self}$ of n-alkanes due to the relatively large molecular size of MY2, and MY2 experiences the full shear viscosity of the medium. As the molecular weight of n-alkane increases, $D_{self}$ of n-alkanes decreases much faster than $D_{MY2}$ and at the higher molecular weights of n-alkane, MY2 diffuses faster than the solvent fluctuations. Therefore there is a large reduction of friction in longer chains compared to the shorter chains, which enhances the diffusion of MY2. The calculated friction constants of MY1 and MY2 in liquid n-alkanes supported these observations. We deem that this is the origin of the so-called“solventoligomer”transition.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2009년도 정기 학술대회
• /
• 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과 비교 검증한다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 추계학술대회
• /
• 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.

• 한국재료학회지
• /
• 제16권2호
• /
• pp.92-98
• /
• 2006

Molecular dynamics simulations were performed to study interface structures between an $Al_2O_3$ crystalline phase and a interface phase of $CaAl_2Si_2O_8$. We calculated atomic structures and excess interface energies in systems with different thicknesses of the interface film. It was found that excess interface energies at first readily decreased with increasing film thickness, but increased for larger thicknesses of more than 2 nm. The excess energies of $Al_2O_3/CaAl_2Si_2O_8$ interfaces exhibit a minimum at a thickness around 1 nm. In this range of film thicknesses, the atoms in the interface film show a short-range ordered structure and slow diffusion rather than the random structure and rapid diffusion expected to an observation of an equilibrium thickness for interface films in ceramics.

• 한국군사과학기술학회지
• /
• 제14권1호
• /
• pp.147-153
• /
• 2011

A series of nonequilibrium molecular dynamics(NEMD) simulations are performed to investigate the kinetic energy and velocity distributions of molecules in shock waves. In the simulations, argon molecules are used as a medium gas through which shock waves are propagating. The kinetic energy distribution profiles reveals that as a strong shock wave whose Mach number is 27.1 is applied, 39.6% of argon molecules inside the shock wave have larger kinetic energy than molecular ionization energy. This indicates that an application of a strong shock wave to argon gas can give rise to an intense light. The velocity distribution profiles in z direction along which shock waves propagate clearly represent two Maxwell-Boltzmann distributions of molecular velocities in two equilibrium regions and one bimodal velocity distribution profile that is attributed to a nonequilibrium region. The peak appearing in the directional temperature in z direction is discussed on a basis of the bimodal velocity distribution in the nonequilibrium region.

• 설비공학논문집
• /
• 제19권12호
• /
• pp.883-888
• /
• 2007

As a starting point of investigating what molecular dynamic simulations can reveal about the nature of atomic level of heating and cooling process, argon described by the LJ potential is considered. Stepwise heating and cooling of constant rates are simulated in the NPT (constant number, pressure and temperature) ensemble. Hysteresis is found due to the superheating and supercooling. Drastic change of volume and energy is involved with phase change, but the melting point can not be obtained by simply observing the changes of these quantities. Since liquid and solid phases can co-exist at the same temperature, Gibbs free energy should be calculated to find the temperature where the Gibbs free energy of liquid is equal to that of the solid since the equilibrium state is the state of minimum Gibbs free energy. The obtained melting temperature, $T^*=0.685$, is close to that of the experiment with only 2% error.