• Bulletin of the Korean Chemical Society
• /
• v.4 no.2
• /
• pp.87-91
• /
• 1983

Molecular dynamic method has been applied to a single polymer chain immersed in a solvent. The interactions for the pairs, of two solvent molecules (SS), of a chain element and a solvent molecules (CS), and of two non-neighbor chain elements (CC) are given by the Lennard-Jones potential, and the interaction between two bonded chain elements is given by a harmonic potential. We changed the CS interaction parameter ${\varepsilon}_{CS}$ to 0.5, 1.0 and 2.0 times of the SS interaction ${\varepsilon}_{SS}$. We calculated the pair correlation functions for the SS, CS, and CC pairs, end-to-end distance and radius of gyration with the varying ${\varepsilon}_{CS}$ parameters. The results showed that a phase separation occurs between the polymer and solvent in the 0.5 system where ${\varepsilon}_{CS}$ = 0.5 {\varepsilon}_{SS}$. The autocorrelation functions for end-to-end distance and radius of gyration were also calculated.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1382-1387
• /
• 2004

From the viewpoint of a macro state, there is no thermal boundary resistance (TBR) at an interface if both surfaces at an interface are perfectly contacted. However, recent molecular dynamics (MD) studies reported that there still exists the TDR at the interface in an ideal epitaxial superlttice. Our previous studies suggested the model to predict the TBR not only quantitatively also qualitatively in superlattices. The suggested model was based on the classical theory of a wave reflection, and provided highly satisfactory results for an engineering purpose. However, it was not the complete model because our previous model was derived by considering only the effects from a mass ratio and a potential ratio of two species. The interaction of two species presented by the Lennard-Jones (L-J) potential is governed by the mutual ratio of the masses, the potential well depths, and the diameters. In this study, we performed the preliminary simulations to investigate the effect resulting from the diameter ratio of two species for the completion of our model and confirmed that it was also a ruling factor to the TBR at an interface in superlattices.

• Bulletin of the Korean Chemical Society
• /
• v.33 no.9
• /
• pp.3039-3042
• /
• 2012

Pressure tensors at the planar surface of liquid-vapor argon are evaluated from the virial theorem, Irving-Kirkwood, and Harasima versions using a test-area molecular dynamics simulation method through a Lennard-Jones intermolecular potential at two temperatures. We found that the normal and transverse components of the pressure tensor, $p_N(z)$ and $p_T(z)$, obtained from the virial theorem and Harasima version are essentially the same. The normal component of the pressure tensor from Irving-Kirkwood version, $p_N^{IK}(z)$, is shown to be a nearly constant at the lower temperature, independent of z, as agreed in a previous study, but not for $p_N^H$(z), while the transverse components, $p_T^{IK}(z)$ and $p_T^H(z)$, are almost the same. The values of surface tension for both versions computed from $p_N(z)-p_T(z)$ are also the same and are fully consistent with the experimental data.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2003.10a
• /
• pp.347-354
• /
• 2003

The dynamics of fluid flow through nanomachines is completely different from that of continuum. In this study, molecular dynamics simulations were performed for the flow of helium, neon, argon inside carbon(graphite) nanotubes of several sizes. The fluid was introduced into the nanotube at a given initial velocity according to given temperature. Diffusion coefficients were evaluated by Green-Kubo equation derived from Einstein relationship. The behaviour of the fluid was strongly dependent on the density of fluid and tube diameter, not on the tube length. It was found that the diffusion Coefficients increased With decreasing the density of molecules and increasing the diameter and temperature.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.11
• /
• pp.2436-2441
• /
• 2002

Finite element technique considering adhesive forces is proposed and applied to analyze the behavior of elastic hemispherical asperity adhesively contacting the plane surface of semi -infinite rigid body. It is demonstrated that the finite element model simulates interfacial phenomena such as jump -to-contact and adhesion hysteresis that cannot be simulated with the currently available adhesive contact continuum models. This simulation aiso provides valuable information on contact pressure, contact region and stress distributions. This technique is anticipated to be utilized in designing a low-adhesion surface profile for MEMS/NEMS applications since various contact geometries can be analyzed with this technique.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.337-342
• /
• 2007

In this study, molecular dynamics simulation of nano imprint lithography in which patterned stamp is pressed onto amorphous polyethylene(PE) surface are performed to study the behaviour of polymer. Force fields including bond, angle, torsion, and Lennard Jones potential are used to describe the inter-molecular and intra-molecular force of PE molecules and stamp, substrate. Periodic boundary condition is used in horizontal direction and canonical NVT ensemble is used to control the system temperature. As the simulation results, the behaviour of polymer is investigated during the imprinting process. The mechanism of polymer deformation is studied by means of inspecting the surface shape, volume, density, atom distribution. Deformation of the polymer resist was found for various of the stamp geometry and the alignment state of the polymer molecules.

• Journal of the Korean Society for Precision Engineering
• /
• v.14 no.3
• /
• pp.122-129
• /
• 1997

Machining technology has emerged to the point of performing atomic-scale fabrication. In tail paper atomic-scale machining characteristics are investigated by using Molecular Statics simulation method. The cutting model used in this work simulates machining with tools such as an AFM. It is shown that built-up edge formation and cutting forces depend on tool tip geometry. Also, the material flow during cutting is shown for various cutting conditions such as depth of cut, rake angle, and edge radius of tool.

• Journal of Advanced Marine Engineering and Technology
• /
• v.28 no.3
• /
• pp.441-450
• /
• 2004

The velocity distribution of gas molecules from the experimental results was confirmed as the same with the Maxwell-Boltzmann's theoretical results within the experimental error. This study is on the realization of the Maxwell-Boltzmann's velocity distribution of gas molecules by the molecular dynamics(MD) method. The Maxwell-Boltzmann's velocity distribution of gas molecules is extremely important to confirm the equilibrium state because the properties of a thermodynamic system shall be obtained from the system's equilibrium configuration in the MD method. This study is the first trial in the successive researches to calculate the properties of a thermodynamic system by the computer simulations. We confirmed that the maxwell-boltzmann's velocity distribution is developed in some transient time after starting a simulation and dependent on the size of a system. Also it is found that the velocity distribution has no relation with an initial configuration of gas molecules.

• Bulletin of the Korean Chemical Society
• /
• v.35 no.12
• /
• pp.3527-3531
• /
• 2014

In this paper, we report thermodynamic and transport properties (diffusion coefficient, viscosity, and thermal conductivity) of diatomic gases ($H_2$, $N_2$, $O_2$, and $Cl_2$) at 273.15 K and 1.00 atm by performing molecular dynamics simulations using Lennard-Jones intermolecular potential and modified Green-Kubo formulas. The results of self-diffusion coefficients of diatomic gases obtained from velocity auto-correlation functions by Green-Kubo relation are in good agreement with those obtained from mean square displacements by Einstein relation. While the results for viscosities of diatomic gases obtained from stress auto-correlation functions underestimate the experimental results, those for thermal conductivities obtained from heat flux auto-correlation functions overestimate the experimental data except $H_2$.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.660-661
• /
• 2008

Nanotubes fabricated with diverse materials show different hydrophobic properties. The hydrophobic property is one of key properties for possible applications to ion channels due to their affinity. This study focuses on the structures of water molecules inside nanotubes with varying hydrophobicity using molecular dynamics simulation. Hydrophobicity here is determined by varying the attraction term in Lennard-Jones potential. The number of water molecules inside hydrophilic nanotubes increase, as expected, and their mobilities also increase. This trend is rather discrete with increasing number of water molecules and this discreteness is attributed to hydrogen bond. We plan to perform energy analysis to understand these structural results.