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

• Journal of the Korean Vacuum Society
• /
• v.19 no.2
• /
• pp.121-127
• /
• 2010

Neutral gas in a Hall-effect thruster in a small satellite is simulated using a molecular dynamics code. Investigated are neutral density, pressure, axial average velocity, and temperature for the variation of diffusive reflection ratio, initial gas temperature, and channel length. Expected through this research are improving of discharge simulation through the neutral simulation and understanding of real system.

• Nuclear Engineering and Technology
• /
• v.50 no.6
• /
• pp.907-914
• /
• 2018

In this article, we conducted molecular dynamics simulations to investigate the effect of applied strain and temperature on irradiation-induced damage in alpha-zirconium. Cascade simulations were performed with primary knock-on atom energies ranging between 1 and 20 KeV, hydrostatic and uniaxial strain values ranging from -2% (compression) to 2% (tensile), and temperatures ranging from 100 to 1000 K. Results demonstrated that the number of defects increased when the displacement cascade proceeded under tensile uniaxial hydrostatic strain. In contrast, compressive strain states tended to decrease the defect production rate as compared with the reference no-strain condition. The proportions of vacancy and interstitial clustering increased by approximately 45% and 55% and 25% and 32% for 2% hydrostatic and uniaxial strain systems, respectively, as compared with the unstrained system, whereas both strain fields resulted in a 15-30% decrease in vacancy and interstitial clustering under compressive conditions. Tensile strains, specifically hydrostatic strain, tended to produce larger sized vacancy and interstitial clusters, whereas compressive strain systems did not significantly affect the size of defect clusters as compared with the reference no-strain condition. The influence of the strain system on radiation damage became more significant at lower temperatures because of less annealing than in higher temperature systems.

• Coupled systems mechanics
• /
• v.6 no.1
• /
• pp.41-54
• /
• 2017

Shape-memory alloys (SMA) have interesting behaviors and important mechanical properties due to the solid-solid phase transformation. These phenomena are dominated by the evolution of microstructures. In recent years, the microstructures in SMAs have been studied extensively and modeled using molecular dynamics (MD) simulations. However, it remains difficult to identify the crystal variants in the simulation results, which consist of large numbers of atoms. In the present work, a method is developed to identify the austenite phase and the monoclinic martensite crystal variants in MD results. The transformation matrix of each lattice is calculated to determine the corresponding crystal variant. Evolution of the volume fraction of the crystal variants and the microstructure in Ni-Ti SMAs under thermal and mechanical boundary conditions are examined. The method is validated by comparing MD-simulated interface normals with theoretical solutions. In addition, the results show that, in certain cases, the interatomic potential used in the current study leads to inconsistent monoclinic lattices compared with crystallographic theory. Thus, a specific modification is applied and the applicability of the potential is discussed.

• Coupled systems mechanics
• /
• v.3 no.4
• /
• pp.329-344
• /
• 2014

In this work, quantum molecular dynamics simulations (QMD) are preformed to study the hydrogen molecules in three types of carbon nanostructures, $C_{60}$ fullerene, (5,5) and (9,0) carbon nanotubes and graphene layers. Interactions between hydrogen and the nanostructures is of importance to understand hydrogen storage for the development of hydrogen economy. The QMD method overcomes the difficulties with empirical interatomic potentials to model the interaction among hydrogen and carbon atoms in the confined geometry. In QMD, the interatomic forces are calculated by solving the Schrodinger's equation with the density functional theory (DFT) formulation, and the positions of the atomic nucleus are calculated with the Newton's second law in accordance with the Born-Oppenheimer approximation. It is found that the number of hydrogen atoms that is less than 58 can be stored in the $C_{60}$ fullerene. With larger carbon fullerenes, more hydrogen may be stored. For hydrogen molecules passing though the fullerene, a particular orientation is required to obtain least energy barrier. For carbon nanotubes and graphene, adsorption may adhere hydrogen atoms to carbon atoms. In addition, hydrogen molecules can also be stored inside the nanotubes or between the adjacent layers in graphite, multi-layer graphene.

• The Journal of the Korea institute of electronic communication sciences
• /
• v.14 no.1
• /
• pp.147-154
• /
• 2019

In this research, a molecular dynamics system was designed and developed to calculate trajectories of molecules in nozzles and turbin blades. The Lennard-Jones potential model was used to approximate the interaction between a pair of molecules and the Verlet integration is used as a numerical method to integrate Newton's equations of motion. To compute Lennard-Jones potential functions, for the number of molecules N, the computation complexity $O(N^2)$ for interactions of all pairs of molecules is reduced to O(N) by using cutoff radius $r_c$. This was implemented to save CPU times.

• The Journal of the Korea institute of electronic communication sciences
• /
• v.14 no.2
• /
• pp.425-432
• /
• 2019

Molecular dynamics simulations were carried out to analyze the impact on turbine rotational power based on nozzle shapes. y varying two parameters, turbine rotational velocity w and initial velocities ($V_x,V_y$) for generating molecules, turbine rotation velocities of molecules in turbine entrance arc region and turbine exit arc region were calculated. Average momentums were calculated over two regions, respectively and the differences (AMD) were compared and analyzed. The optimal range of w values to enhance AMD was investigated over 4 different nozzle shapes used, and an excellent nozzle shape to achieve more turbine rotational power was drawn.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.18 no.5
• /
• pp.60-65
• /
• 2019

Since nanofluids (NFs), which are a mixture of a small amount of nanoparticles and a bulk liquid solvent, were first proposed by Stephen Choi at the Argonne National Lab in 1995, they have been considered for use in many technical studies of power cooling systems and their practical application due to their high thermal conductivity and heat transfer coefficients compared to conventional coolants. Although nanofluids are a well-known form of engineering fluid that show great promise for use in future cooling systems, their underlying physics as demonstrated in experiments remain unclear. One proven method of determining the heat transfer performance of nanofluids is measuring the concentration of nanoparticles in a mixture. However, it is experimentally inefficient to build testbeds to systematically observe particle distributions on a nanoscale. In this paper, we demonstrate the distribution of nanoparticles under a temperature gradient in a solution using molecular dynamics simulations. First, temperature profiles based on substrate temperature are introduced. Following this, the radial pair distribution functions of pairs of nanoparticles, solvents, and substrates are calculated. Finally, the distribution of nanoparticles in different heating regions is determined.

• Nuclear Engineering and Technology
• /
• v.53 no.5
• /
• pp.1587-1592
• /
• 2021

The effects of the symmetrical tilt Σ5 (310) <001> grain boundary (GB) on the evolution of radiation-induced point defects in pure δ-plutonium (Pu) were studied by Molecular dynamics (MD) simulations. The evolution of radiation-induced point defects was obtained when primary knock-on atom (PKA) was respectively set as -15 Å and 15 Å far from the GB and the number of residual defects was obtained as the distance from PKA to GB was changed. According to the results, compared with vacancies, interstitial atoms were more easily absorbed by GB. In addition, the formation energy of point defects was also calculated. The results showed that there was almost no difference for the formation energy of vacancies in the all matrix. However, the formation energy of interstitial atoms close to the GB was lower than that in the other bulk regions.

• Bulletin of the Korean Chemical Society
• /
• v.12 no.3
• /
• pp.309-315
• /
• 1991

The thermal transport coefficients-the self-diffusion coefficient, shear viscosity, and thermal conductivity-of liquid argon at 94.4 K and 1 atm are calculated by non-equilibrium molecular dynamics (NEMD) simulations of a Lennard-Jones potential and compared with those obtained from Green-Kubo relations using equilibrium molecular dynamics (EMD) simulations and with experimental data. The time-correlation functions-the velocity, pressure, and heat flux auto-correlation functions-of liquid argon obtained from the EMD simulations show well-behaved smooth curves which are not oscillating and decaying fast around 1.5 ps. The calculated self-diffusion coefficient from our NEMD simulation is found to be approximately 40% higher than the experimental result. The Lagrange extrapolated shear viscosity is in good agreement with the experimental result and the asymptotic formula of the calculated shear viscosities seems to be an exponential form rather than the square-root form predicted by other NEMD studies of shear viscosity. The agreement for thermal conductivity between the simulation results (NEMD and EMD) and the experimental result is within statistical error. In conclusion, through our NEMD and EMD simulations, the overall agreement is quite good, which means that the Green-Kubo relations and the NEMD algorithms of thermal transport coefficients for simple liquids are valid.