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

• Interaction and multiscale mechanics
• /
• v.1 no.4
• /
• pp.397-409
• /
• 2008

The interfacial thermal resistance of pristine and defective carbon nanotubes (CNTs) embedded in low-density polyethylene matrix is studied in this paper. Interface thermal resistance in nanosystems is one of the most important factors that lead to the large variation in thermal conductivities in literature and the novelty of this paper lies in the estimation of the interfacial thermal resistance for defective nanotubes-systems. Thermal properties of CNT nanostructures are estimated using molecular dynamics (MD) simulations and the simulations were carried out for various temperatures by rescaling the velocities of carbon atoms in the nanotube. This paper also deals with the mesoscale thermal conductivities of composite systems, using effective medium theories by considering the size effect in the form of interfacial thermal resistance and also using the conventional micromechanical methods like Hashin-Shtrikman bounds and Wakashima-Tsukamoto estimates.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.7
• /
• pp.60-68
• /
• 2004

High velocity impacts are accompanied with large deformations, which generate a large amount of heat due to plastic works, resulting in a significant temperature rise of the material. Because the elevated temperature affects the dynamic properties of materials, it is important to predict the temperature rise during high-stram-rate deformations. Both existing vacancies and excess vacancies are credited to the stored energy, yet it is difficult to distinguish one from another in contribution to the stored energy using macroscopic level materials models. In this study, an atomistic material model for fee materials such as copper is set up to calculate the stored energy using molecular dynamics (MD) simulations. It is concluded that excess vacancies play an important role for the stored energy during a high-strain-rate deformation.

• Membrane Journal
• /
• v.30 no.4
• /
• pp.242-251
• /
• 2020

Computer simulation tools mainly used for polymer materials and polymeric membranes are divided into various fields depending on the size of the object to be simulated and the time to be simulated. The computer simulations introduced in this review are classified into three categories: Quantum mechanics (QM), molecular dynamics (MD), and mesoscale modeling, which are mainly used in computational material chemistry. The computer simulation used in polymer research has different research target for each kind of computational simulation. Quantum mechanics deals with microscopic phenomena such as molecules, atoms, and electrons to study small-sized phenomena, molecular dynamics calculates the movement of atoms and molecules calculated by Newton's equation of motion when a potential or force of is given, and mesoscale simulation is a study to determine macroscopically by reducing the computation time with large molecules by forming beads by grouping atoms together. In this review, various computer simulation programs mainly used for polymers and polymeric membranes divided into the three types classified above will be introduced according to each feature and field of use.

• Composites Research
• /
• v.31 no.6
• /
• pp.371-378
• /
• 2018

In this study, the thermoelastic properties of poly lactic acid (PLA) based nanocomposites are predicted by molecular dynamics (MD) simulation and a micromechanics model. The stereocomplex mixed with L-lactic acid (PLLA) and D-lactic acid (PDLA) is modeled as matrix phase and a single walled carbon nanotube is embedded as reinforcement. The glass transition temperature, elastic moduli and thermal expansion coefficients of pure matrix and nanocomposites unit cells are predicted though ensemble simulations according to the hydrolysis. In micromechanics model, the double inclusion (D-I) model with a perfect interface condition is adopted to predict the properties of nanocomposites at the same composition. It is found that the stereocomplex nanocomposites show prominent improvement in thermal stability and interfacial adsorption regardless of the hydrolysis. Moreover, it is confirmed from the comparison of MD simulation results with those from the D-I model that the interface between CNT and the stereocomplex matrix is slightly weak in nature.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.4A
• /
• pp.529-536
• /
• 2008

The failure behavior of structures is changed under different loading rates, which might arise from the rate dependency of materials. This phenomenon has been focused in the engineering fields. However, the failure mechanism is not fully understood yet, so that it is hard to be implemented in numerical simulations. In this study, the numerical experiments to a brittle material are simulated by the Molecular Dynamics (MD) for understanding the rate dependent failure behavior. The material specimen with a notch is modeled for the compact tension test simulation. Lennard-Jones potential is used to describe the properties of a brittle material. Several dynamic failure features under 6 different loading rates are achieved from the numerical experiments, where remarkable characteristics such as crack roughness, crack recession/arrest, and crack branching are observed during the crack propagation. These observations are interpreted by the energy inflow-consumption rates. This study will provides insight about the dynamic failure mechanism under different loading rates. In addition, the applicability of the MD to the macroscopic mechanics is estimated by simulating the previous experimental research.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.97-102
• /
• 2009

Translocation of biopolymers such as DNA and RNA through a nano-pore is an important process in biotechnology applications. The translocation process of a biopolymer through an artificial nano-pore in the presence of a fluid solvent is simulated. The polymer motion is simulated by Langevin molecular dynamics (MD) techniques while the solvent dynamics are taken into account by lattice-Boltzmann method (LBM). The hydrodynamic interactions are considered explicitly by coupling the polymer and solvent through the frictional and the random forces. From simulation results we found that the hydrodynamic interactions between polymer and solvent speed-up the translocation process. The translocation time ${\tao}_T$ scales with the chain length N as ${{\tau}_T}^{\propto}N^{\alpha}$. The value of scaling exponents($\alpha$) obtained from our simulations are $1.29{\pm}0.03$ and $1.41{\pm}0.03$, with and without hydrodynamic interactions, respectively. Our simulation results are in good agreement with the experimentally observed value of $\alpha$, which is equal to $1.27{\pm}0.03$, particularly when hydrodynamic interaction effects are taken into account.

• Current Optics and Photonics
• /
• v.1 no.4
• /
• pp.278-283
• /
• 2017

We theoretically study the dynamics of water tetramer in intense femtosecond laser pulses with different frequencies. The simulations are carried out by incorporating the molecular dynamics method non-adiabatically into the time-dependent local-density approximation (TDLDA-MD). Three typical scenarios of water tetramer including the normal vibration with enlarged OH bonds, free OH bonds breaking and the pure Coulomb explosion are presented by investigating the electronic and ionic dynamics. The result indicates that the ionization is enhanced and the corresponding fragmentation effect as well as the damping of the dipole moment are found more notably when increasing the laser frequency especially when the frequency falls in the resonant region of the absorption spectra. The study of the level depletion reveals that the ratio of the emission amount from different levels can be controlled by changing the laser frequency referring to the Keldysh parameter.

• Bulletin of the Korean Chemical Society
• /
• v.29 no.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.

• Proceedings of the Korean Institute of Electrical and Electronic Material Engineers Conference
• /
• 2008.06a
• /
• pp.379-380
• /
• 2008

Molecular Dynamics (MD) simulations have been used to model the dynamics of the charge-compensating sodium ions in the non-stoichiometric hollandite Nax$(Ti_{8-x}Cr_x)O_{16}$. These interstitial ions reside in 'tunnels' in the crystal structure and move under the forces of both the ions making up the cage structure and the many body interactions of the other sodium ions in the tunnel. The Velocity Autocorrelation Function (VAF) of the sodium ions is calculated for a range of temperature from 250K to 1000K and converted into the linear ac-conductivity and ac-susceptibility response via Fourier transformation. A peak is found in the conductivity around $6\times10^{12}$ Hz that has some of the character of a Poley absorption. Here it is shown to be due to an harmonically coupled site vibrations of the sodium atoms, which extend only over a limited range. At frequencies below the peak the conductivity tends towards a constant i.e. dc value corresponding to a constant flow of ions through the simulation cell. At high temperatures the conductivity due to this ion transport process behaves like a metal with an insulator to metal transition occurring around a specific temperature.