• Tribology and Lubricants
• /
• v.32 no.1
• /
• pp.9-17
• /
• 2016

Contact between solid surfaces is one of the most important factors that influence dynamic behavior in micro/nanoscale. Although numerous theories and experimental results on contact behavior have been proposed, a thorough investigation for nanomaterials is still not available owing to technical difficulties. Therefore, molecular dynamics simulation was performed to investigate the contact behavior of nanomaterials, and the application of conventional contact theories to nanoscale was assessed in this work. Particularly, the contact characteristics of cylindrical nanowires were examined via simulation and contact theories. For theoretical analysis, various contact models were utilized and work of adhesion, Hamaker constant and elastic modulus those are required for calculation of the models were obtained from both indentation simulation and tensile simulation. The contact area of the cylindrical nanowire was assessed directly through molecular dynamics simulation and compared with the results obtained from the theories. Determination of the contact area of the nanowires was carried out via simulation by counting each atom, which is within the equilibrium length. The results of the simulation and theoretical calculations were compared, and it was estimated that the discrepancy in the results calculated between the simulation and the theories was less than 10 except in the case of the smallest nanowires. As the result, it was revealed that contact models can be effectively utilized to assess the contact area of nanomaterials.

• Korean Chemical Engineering Research
• /
• v.46 no.6
• /
• pp.1087-1094
• /
• 2008

This article presents a multi-scale simulation approach starting from the molecular level for the adsorption process development, specifically, in n-hexane adsorption on activated carbon. A grand canonical Monte-Carlo(GCMC) method is used for the prediction of adsorption isotherms of n-hexane on activated carbon at the molecular level. Geometric effects and hydrodynamic properties of the adsorption column are examined by means of the two dimensional CFD(computational fluid dynamics) simulation. The adsorption isotherms from the molecular simulation and the axial diffusivity from the CFD simulation are exploited for the process simulation where the elution curve of n-hexane is obtained. For the first moment(mean residence time) of the pulse-response with respect to temperature and flowrate, the process simulation results obtained from this three-steps multiscale simulation approach show a good agreement with experimental data within 20% of maximum difference. The multi-scale simulation approach addressed in this study will be useful to accelerate the adsorption process development, while reducing the number of experiments required.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.8
• /
• pp.559-564
• /
• 2009

The variation in the shape of water droplet moving on atomistically smooth solid surface in the presence of a constant body force is simulated using molecular dynamics simulation. We investigated how the advancing and receding contact angle of the moving water droplet changes on a solid surface having various characteristic energies. From the MD simulation results, we obtained the density profile defined as the number of water molecules at a given position. Then, assuming the water droplet periphery to be a circle, we calculated the contact angles by using a nonlinear fitting of the half-density contour line. The present simulation clearly shows the different profile of the advancing and receding contact angle for these three different interaction potential between the water droplet and the solid surface.

• Korean Chemical Engineering Research
• /
• v.49 no.3
• /
• pp.361-366
• /
• 2011

By using Monte Carlo simulation method we developed a new molecular simulation software which can be used to predict the thermodynamic properties of organic compounds. Starting from molecular structure and intermolecular potential function, rigorous statistical mechanical principles give a probability distribution for the behavior of a system containing many molecules, which enables us to calculate macroscopic thermodynamic properties of the system. The software developed in this work, cheMC, is based on Windows platform providing with easy access. One can efficiently administrate simulations by using an intuitive interface equipped with visualization tool and chart generation. It is expected that molecular simulations supplement the equation of state approach and will play a more important role in the study of thermodynamic properties.

• Applied Chemistry for Engineering
• /
• v.33 no.6
• /
• pp.551-556
• /
• 2022

The desalinated water obtained by the water treatment process based on the membrane is attracting a lot of attention as a promising technology that can solve the global water shortage problem. Reverse osmosis membrane-based desalination, one of the most widely used desalination processes, is a technology that desalinates abundant seawater on Earth, thus having great potential in the desalination industry. To improve the performance of the desalination process, it is necessary to understand the reverse osmosis mechanism of the membrane at the atomic/molecular level. In this review, we introduce molecular simulation, which plays an important role in material research today, and the roles of computational simulation at the atomic/molecular level in the development of reverse osmosis membranes.

• Bulletin of the Korean Chemical Society
• /
• v.34 no.8
• /
• pp.2466-2472
• /
• 2013

Human immunodeficiency virus type-1 (HIV-1) is a causative agent of Acquired immunodeficiency syndrome (AIDS), which has affected a large population of the world. Viral envelope glycoprotein (gp120) is an intrinsic protein for HIV-1 to enter into human host cells. Molecular docking guided molecular dynamics (MD) simulation was performed to explore the interaction mechanism of heterobiaryl derivative with gp120. MD simulation result of inhibitor-gp120 complex demonstrated stability. Our MD simulation results are consistent with most of the previous mutational and modeling studies. Inhibitor has an interaction with the CD4 binding region. Van der Waals interaction between inhibitor and Val255, Thr257, Asn425, Met426 and Trp427 were important. This preliminary MD model could be useful in exploiting heterobiaryl-gp120 interaction in greater detail, and will likely to shed lights for further utilization in the development of more potent inhibitors.

• Journal of Mechanical Science and Technology
• /
• v.16 no.11
• /
• pp.1440-1447
• /
• 2002

Adhesion of a hemispherical tip to the flat surface in nano-structures is simulated using the molecular dynamics technique. The tip and plates are modeled with the Lennard-Jones molecules. The simulation focuses on the deformation of the tip. Detailed descriptions on the evolution of interaction force, the energy dissipation due to adhesion hysteresis, the forma- tion-growth-breakage of adhesive junction as well as the evolution of molecular distribution during the process are presented. The effects of the tip size, the maximum tip approach, the tip temperature, and the affinity between the tip and the mating plate are also discussed.

• Interaction and multiscale mechanics
• /
• v.2 no.4
• /
• pp.353-374
• /
• 2009

A new seamless multiscale simulation was developed for coupling the continuum model with its molecular dynamics. Kriging-based Finite Element Method (K-FEM) is employed to model the continuum base of the entire domain, while the molecular dynamics (MD) is confined in a localized domain of interest. In the coupling zone, where the MD domain overlaps the continuum model, the overall Hamiltonian is postulated by contributions from the continuum and the molecular overlays, based on a quartic spline scaling parameter. The displacement compatibility in this coupling zone is then enforced by the Lagrange multiplier technique. A multiple-time-step velocity Verlet algorithm is adopted for its time integration. The validation of the present method is reported through numerical tests of one dimensional atomic lattice. The results reveal that at the continuum/MD interface, the commonly reported spurious waves in the literature are effectively eliminated in this study. In addition, the smoothness of the transition from MD to the continuum can be significantly improved by either increasing the size of the coupling zone or expanding the nodal domain of influence associated with K-FEM.

• Korean Journal of Materials Research
• /
• v.10 no.8
• /
• pp.581-588
• /
• 2000

Phase transition of the Lennard-Jonesian binary system with a free surface was studied by employing molecular dynamics simulation. The main focus of this study was to investigate the effect of size misfit and solute concentration on phase transition during heating and quenching. For a binary system with a free surface, the melting point and the critical quenching rate decrease as size misfit and solute concentration increase.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2007.10a
• /
• pp.348-351
• /
• 2007

Simulations of the buckling behavior of a single wall carbon nanotube(SWCNT) was carried out using molecular dynamics simulation. Molecular dynamics simulations were done with 1fs of time step. Tersoff's potential function was used as the interatomic potential function since it has been proved to be reliable to describe the C-C bonds in carbon nanotubes. Compressive force was applied by moving the top end of the nanotube at a constant velocity. Buckling behavior under compressive load was observed for (15,15) armchair SWCNTs with 2nm of diameter and 24.9nm of length. Buckling load and critical strain is obtained from the MD simulation. Deformation occurred on the top region of the CNT because of fast downward velocity.