• Journal of Biomedical Engineering Research
• /
• v.24 no.5
• /
• pp.447-451
• /
• 2003

A new method for analyzing the bone-remodeling process using quasi-molecular dynamics was proposed in this study. The effect of pressure due to bone marrow, which could not be considered in previous methods, was also considered in this method. Bone-remodeling response of the 2D vertebral body of lumbar spine to a uniaxial compressive displacement of 1.8564mm. corresponding to approximately 2kN of compressive load, was studied. Converged shape change of the cortical shell and rearrangement of cancellous bone structure matched well with a normal shape of the vertebral body. The calculated responses in the spinal elements also shows closed results compared with experimental results.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2008.11a
• /
• pp.668-669
• /
• 2008

In this paper, we introduce various coarse-grained elastic network modeling (ENM) techniques as a novel computational method for simulating atomic scale dynamics in macromolecules including DNA, RNA, protein, and polymer. In ENM, a system is modeled as a spring network among representative atoms in which each linear elastic spring is well designed to replace both bonded and nonbonded interactions among atoms in the sense of quantum mechanics. Based on this simplified system, a harmonic Hookean potential is defined and used for not only calculating intrinsic vibration modes of a given system, but also predicting its anharmonic conformational change, both of which are strongly related with its functional features. Various nano and bio applications of ENM such as fracture mechanics of nanocomposite and protein dynamics show that ENM is one of promising tools for simulating atomic scale dynamics in a more effective and efficient way comparing to the traditional molecular dynamics simulation.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.433-438
• /
• 2007

Mechanical behavior of copper Nanowire is investigated, An FCC Nanowire model composed of 1,408 atoms is used for NID simulation, Simulations are performed within NVT ensemble setting without periodic boundary conditions, Nose-Poincare MD algorithm is employed to guarantee preservation of Hamiltonian. Numerical tensile tests are carried out with constant strain rate, Stress-strain curve is constructed from the calculated Cauchy stresses and specified strain values, Non-linear behavior appears around $\varepsilon$=0.064, At this instance, starting of structural reorientations are observed.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2011.04a
• /
• pp.253-256
• /
• 2011

본 논문에서는 표면효과와 비선형 탄성효과를 고려한 FCC 나노박막의 순차적 멀티스케일 해석 모델을 제시한다. 표면에서의 구성방정식은 표면응력과 표면탄성계수를 이용하여 선형으로 표시되며, 표면효과를 나타내기 위한 표면물성들은 EAM 포텐셜을 이용한 원자적 계산 방법으로 계산된다. 두께가 얇은 나노박막은 표면응력으로 인하여 면내 방향으로 수축 또는 인장의 변형이 발생하게 된다. 나노박막의 평형상태에서의 변형율은 두께가 얇은 박막의 경우 재료가 선형 탄성 영역을 벗어나는 값을 가지는 경우가 많으므로 나노박막의 해석시 벌크 영역의 비선형 탄성 효과를 고려해야 한다. 이러한 비선형 탄성 효과를 고려하기 위해 본 연구에서는 FCC 구조를 가지는 금속의 비선형 탄성 모델을 제시하고, EAM 포텐셜로 계산된 응력과 탄성 계수를 이용하여 매칭 기법을 통하여 비선형 탄성 모델의 계수들을 결정한다. 또한 Cauchy-Born Rule 모델과 분자동역학 전산모사를 통하여 본 연구에서 제안된 비선형 탄성 모델에 대한 검증을 수행한다.

• Proceedings of the Korean Society Of Semiconductor Equipment Technology
• /
• 2007.06a
• /
• pp.203-206
• /
• 2007

연속 전자 모델과 결합된 종래의 분자 동역학 방법은 원자 사이의 힘과 원자의 전기용량에 의해 야기되는 탄소 나노튜브의 구부러지는 성질의 특성을 해석하였다. 탄소 원자의 전기 용량은 탄소 원자의 길이에 따라 변하였다. 본 연구는 11.567nm($L_{CNT}$)의 길이와 $0.9{\sim}1.5nm(H)$의 안쪽 깊이를 가진 (5,5) 탄소 나노튜브 브리지로 MD 시뮬레이션을 수행하였다. 탄소 나노튜브는 금 표면에 부딪힌 후 탄소 나노튜브 브리지는 약 ${\sim}1{\AA}$의 크기로 금 표면에서 진동하며, 크기는 차츰 감소하였다. $H{\leq}1.3nm$일 때, 탄소 나노튜브 브리지는 첫 번째 충돌 후에 금 표면과 계속 접촉해 있었고, $H{\leq}1.4nm$일 때, 탄소 나노튜브 브리지는 몇 번의 충돌 후에 금 표면과 안정한 접촉상태가 되었다. $H/L_{CNT}$가 0.13보다 작을 때, 탄소 나노튜브 초소형 전자기기 메모리는 반영구적인 비활성의 메모리 장치가 되는 반면에 $H/L_{CNT}$가 0.14보다 클 때 탄소 나노튜브 초소형 전자기기 메모리는 휘발성이거나 스위치 장치로 동작할 수 있다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.2
• /
• pp.103-108
• /
• 2019

Non-equilibrium molecular dynamics simulation on the thermal boundary resistance(TBR) of an aluminum(Al)/silicon(Si) interface was performed in the present study. The constant heat flux across the Si/Al interface was simulated by adding the kinetic energy in hot Si region and removing the same amount of the energy from the cold Al region. The TBR estimated from the sharp temperature drop at the interface was independent of heat flux and equal to $5.13{\pm}0.17K{\cdot}m^2/GW$ at 300K. The simulation result was experimentally confirmed by the time-domain thermoreflectance technique. A 90nm thick Al film was deposited on a Si(100) wafer using an e-beam evaporator and the TBR on the film/substrate interface was measured using the time-domain thermoreflectance technique based on a femtosecond laser system. A numerical solution of the transient heat conduction equation was obtained using the finite difference method to estimate the TBR value. Experimental results were compared to the prediction and discussions on the nanoscale thermal transport phenomena were made.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.4
• /
• pp.315-321
• /
• 2012

In this study, a sequential multiscale homogenization method to characterize the effective thermal conductivity of nano particulate polymer nanocomposites is proposed through a molecular dynamics(MD) simulations and a finite element-based homogenization method. The thermal conductivity of the nanocomposites embedding different-sized nanoparticles at a fixed volume fraction of 5.8% are obtained from MD simulations. Due to the Kapitza thermal resistance, the thermal conductivity of the nanocomposites decreases as the size of the embedded nanoparticle decreases. In order to describe the nanoparticle size effect using the homogenization method with accuracy, the Kapitza interface in which the temperature discontinuity condition appears and the effective interphase zone formed by highly densified matrix polymer are modeled as independent phases that constitutes the nanocomposites microstructure, thus, the overall nanocomposites domain is modeled as a four-phase structure consists of the nanoparticle, Kapitza interface, effective interphase, and polymer matrix. The thermal conductivity of the effective interphase is inversely predicted from the thermal conductivity of the nanocomposites through the multiscale homogenization method, then, exponentially fitted to a function of the particle radius. Using the multiscale homogenization method, the thermal conductivities of the nanocomposites at various particle radii and volume fractions are obtained, and parametric studies are conducted to examine the effect of the effective interphase on the overall thermal conductivity of the nanocomposites.

• Composites Research
• /
• v.37 no.2
• /
• pp.132-138
• /
• 2024

In this study, a molecular dynamics simulation study was performed to investigate the size-dependent electroelastic properties of single-walled boron nitride nanotubes(BNNT). To describe the elasticity and polarization of BNNT under mechanical loading, the Tersoff potential model and rigid ion approximation were adopted. For the prediction of piezoelectric constants and Young's modulus of BNNTs, piezoelectric constitutive equations based on the Maxwell's equation were used to calculate the strain-electric displacement and strain-stress relationships. It was found that the piezoelectric constants of BNNTs gradually decreases as the radius of the tubes increases showing a nonnegligible size effect. On the other hand, the elastic constants of the BNNTs showed opposites trends according to the equivalent geometrical assumption of the tubular structures. To establish the structure-property relationships, localized configurational change of the primarily bonded B-N bonded topology was investigated in detail to elucidate the BNNT curvature dependent elasticity.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.14 no.3
• /
• pp.115-122
• /
• 2004

When porous materials are dried, the particles flocculate into fish-net structure in gel phase. In order to exactly analyze the stress distribution of porous materials during drying process, the elastic tensor of microscopic gel structures has to be predicted considering pore shapes as well as porosities of porous materials. The elastic characteristics of porous materials associated with porosities were predicted analyzing microscopic gel structures with circular and cross pores via homogenization method and the drying processes of the electric porous ceramic insulator were simulated using finite element method (FEM). Comparing analysis results between consideration and negligence of pores, the deformed shape and distributions of temperature and moisture were similar but the residual stress was significantly different.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.7 no.6
• /
• pp.93-99
• /
• 2003

In this study, new shock absorbing system is proposed by using nano-technology based on the theoretical analysis. The new shock absorbing system is complementary to the hydraulic damper, having a cylinder-piston-orifice construction. Particularly for new shock absorbing system, the hydraulic oil is replaced by a colloidal suspension, which is composed of a porous matrix and a lyophobic fluid. The matrix of the suspension is consisted of porous micro-grains with a special architecture: they present nano-pores serially connected to micro-cavities. Until now, only experimentally qualitative studies of new shock absorbing system have been performed, but the mechanism of energy dissipation has not been clarified. This paper presents a modeling and theoretical analysis of the new shock absorbing system thermodynamics, nono-flows and energy dissipation. Compared with hydraulic system, the new shock absorbing system behaves more efficiently, which absorb a large amount of mechanical energy, without heating. The theoretical computations agree reasonably well with the experimental results. As a result. the proposed new shock absorbing system was proved to be an effective one, which can replace with the conventional one.