• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1558-1563
• /
• 2004

This work applies the nonequilibrium molecular dynamics simulation method to study a Lennard-Jones liquid thin film suspended in the vapor and calculates the thermal conductivity by linear response function. As a preliminary test, the thermal conductivity of pure argon fluid are calculated by nonequilibrium molecular dynamics simulation. It is found that the thermal conductivity decrease with decreasing the density. When both argon liquid and vapor phase are present, the effects of the system temperature on the thermal conductivity are investigated. It can be seen that the thermal conductivity of liquid-vapor interface is constant with increasing the temperature

• Journal of the Korea Society for Simulation
• /
• v.6 no.1
• /
• pp.53-60
• /
• 1997

Results from a nonequilibrium molecular dynamics (NEMD) simulation are presented for an argon liquid subject to a shear flow. The segmented molecular dynamics method and the subtraction technique used in NEMD program to reduce the thermal fluctuation noise in data are studied with different shear rates. The standard deviation in the shear stress reduced from 0.030 to 0.004 by the segmented molecular dynamics method for 50 repeated segments. On the other hand, the standard deviation of the data remained the same when the subtraction technique was applied, where as the results of shear stress by constant value in a random way.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.13 no.1
• /
• pp.156-163
• /
• 2010

A nonequilibrium molecular dynamics(NEMD) simulation method is developed and applied to study shock waves propagating through argon gas. In this simulation method, shock waves are generated by pushing a piston at a constant speed from one side of a simulation box filled with argon molecules. A linear relationship between piston speeds and shock speeds is observed. Thermodynamic properties including density, temperature, and pressure before and after the shock front are obtained from the simulations and compared with the well-known Rankine-Hugoniot equations based on ideal gases. The comparison shows an excellent agreement, indicating that this NEMD simulation method can be employed to investigate various physical properties of shock waves further.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.14 no.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.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1880-1884
• /
• 2004

To analyze the fluid flow and thermal characteristics in a nanoscale system, the planar Poiseuille flow of a Lennar-Jones liquid through parallel plates formed by fixed atoms is studied using nonequilibrium molecular dynamics simulations. The role of important simulation parameters such as the channel width, the magnitude of external field, the temperatures of the top and bottom plates, and the interaction potential parameter between fluid and wall atoms, which affect flow patterns and heat transfer rate inside the channel, are investigated. Under the various simulation conditions, interesting phenomena deviated from the continuum predictions have found.

• Bulletin of the Korean Chemical Society
• /
• v.21 no.4
• /
• pp.434-440
• /
• 2000

We present the results of computer simulation for the steady shear flows of rodlike molecules using nonequi-librium molecular dynamics simulation (NEMD) method. The model particle is a rigid rod composed of lin-early connected 6-sites and the Lennard-Jones 12-6 potential governs interactions between sites in different molecules. The system of rodlike molecules exhibits the change of orientational structure, that is, isotropic-nematic transition at high shear rates. We elucidate the nature of the ordered system developed from an isotro-pic phase by steady shear through an analysis of various quantities: orientational order parameters, orientational pair correlation functions, orientational distribution function, and snapshots of configurations. The effects of temperature and density on the shear rate dependence of orientational structure are described.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.332-336
• /
• 2007

Temperature is an essential process variable in nanoimprint lithography(NIL) where the temperature varies between room temperature and above the glass transition temperature. To simulate NIL process, we employ both the Nose-Poincare method for temperature controlled molecular dynamics(MD) and force field for polymer material i.e. polymethyl methacrylate(PMMA), which is most widely selected as NIL resist. Nose-Poincare method, which convinces the conservation of Hamiltonian structure and time-reversal symmetry, overcomes the drawbacks inherent in the conventional methods such as Nose thermostat and Nose-Hoover thermostat. Thus, this method exhibits enhanced numerical stability even when the temperature fluctuation is large. To describe PMMA, we adopt the force field which account for bond stretch, bending, torsion, inversion, partial charge, and van der Waals energy.

• Nuclear Engineering and Technology
• /
• v.50 no.5
• /
• pp.731-737
• /
• 2018

In this research paper, the thermal transport in thorium dioxide is investigated by using nonequilibrium molecular dynamics. The thermal conductivity of bulk thorium dioxide was measured to be 20.8 W/m-K, confirming reported values, and the phonon mean free path was estimated to be between 7 and 8.5 nm at 300 K. It was observed that the thermal conductivity of thorium dioxide shows a strong dependency on temperature; the highest thermal conductivity was estimated to be 77.3 W/m-K at 100 K, and the lowest thermal conductivity was estimated to be 4.3 W/m-K at 1200 K. In addition, by simulating thorium dioxide structures with different lengths at different temperatures, it was identified that short wavelength phonons dominate thermal transport in thorium dioxide at high temperatures, resulting in decreased intrinsic phonon mean free paths and minimal effect of boundary scattering while long wavelength phonons dominate the thermal transport in thorium dioxide at low temperatures.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 2002.10b
• /
• pp.73-75
• /
• 2002

Recently, remarkable agreement has been reported between nonequilibrium molecular dynamics simulation and high-pressure Couette rheometry on squalane. We utilized the parameters obtained from this unique collaboration along with high-pressure viscometer measurements to calculate the elastohydrodynamic traction curve. A comparison with measured traction at 1.29 GPa shows excellent agreement, confirming the validity of the measurements and simulations. It should no longer be necessary to invoke a different rheological response to explain film thickness and traction.