A Molecular Dynamics Simulation: the Effect of Finite Size on the Thermal Conductivity in a Single Crystal Silicon

Non-equilibrium molecular dynamics (NEMD) simulations are performed to calculate thermal conductivity. The environment-dependent interatomic potential (EDIP) potential on crystal silicon is adopted as a model system. The issues are related to nonlinear response, local thermal equilibrium and statistical averaging. The simulation results by non-equilibrium molecular dynamics show that the calculated thermal conductivity decreases almost linearly as the film thickness reduced at the nanometre scale. The effect of size on the thermal conductivity is also obtained by a theoretic analysis of the kinetic theory and formulas of the heat capacity. The analysis reveals that the contributions of phonon mean free path (MFP) and phonon number in a finite cell to thermal conductivity are very important.

A Constrained Particle Dynamics for Continuum-Particle Hybrid Method in Micro- and Nano-Fluidics

A hybrid method of continuum and particle dynamics is developed for micro- and nano-fluidics, where fluids are described by a molecular dynamics (MD) in one domain and by the Navier-Stokes (NS) equations in another domain. In order to ensure the continuity of momentum flux, the continuum and molecular dynamics in the overlap domain are coupled through a constrained particle dynamics. The constrained particle dynamics is constructed with a virtual damping force and a virtual added mass force. The sudden-start Couette flows with either non-slip or slip boundary condition are used to test the hybrid method. It is shown that the results obtained are quantitatively in agreement with the analytical solutions under the non-slip boundary conditions and the full MD simulations under the slip boundary conditions.

Dynamics Of Fine Particles In Liquid-Solid Fluidized Beds

On the basis of the Local Equilibrium Model (LEM), fine particles with large Richardson-Zaki exponent n show, under certain conditions during bed expansion and collapse, different dynamic behavior from particles with small n. For an expansion process there may be a concentration discontinuity propagating upward from the distributor, and, on the contrary, for a collapse process there may be a progressively broadening and upward-propagating continuous transition zone instead of discontinuity. The predictions of the bed height variation and the discontinuity trace have been validated experimentally. (c) 2007 Chinese Society of Particuology and Institute of Process Engineering, Chinese Academy of Sciences. Published by Elsevier B.V. All rights reserved.

Tolman Effect On Fluid Dynamics In Carbon Nanotubes

The numerical solutions of or(R) given by two different methods (Samsonov et al., 2003; and Lu et al., 2005) are compared with the result that they are coincident closely (the difference is within 4%). We conclude that it is necessary to consider the Tolman correction in the calculation of fluid dynamics in carbon nanotubes. Although our conclusion is the same as that of Prylutskyy et al. (2005), the sign of our Tolman correction is opposite to theirs, and the difference can be attributed to the errors appeared in the paper of Prylutskyy et al.

Molecular Dynamics And Density Functional Studies Of A Body-Centered-Tetragonal Polymorph Of Zno

We report a previously unknown body-centered-tetragonal structure for ZnO. This structure results from a phase transformation from wurtzite in [0001]-oriented nanorods during uniaxial tensile loading and is the most stable phase for ZnO when stress is above 7 GPa. The stress-induced phase transformation has important implications for the electronic, piezoelectric, mechanical, and thermal responses of ZnO. The discovery of this polymorph brings about a more complete understanding of the extent and nature of polymorphism in ZnO. A crystalline structure-load triaxiality map is developed to summarize the relationship between structure and loading.

Size-Dependent Elastic Properties Of Ni Nanofilms By Molecular Dynamics Simulation

Size-dependent elastic properties of Ni nanofilms are investigated by molecular dynamics ( MD) simulations with embedded atom method (EAM). The surface effects are considered by calculating the surface relaxation, surface energy, and surface stress. The Young's modulus and yield stress are obtained as functions of thickness and crystallographic orientation. It is shown that the surface relaxation has important effects on the the elastic properties at nanoscale. When the surface relaxation is outward, the Young's modulus decreases with the film thickness decreasing, and vice versa. The results also show that the yield stresses of the films increase with the films becoming thinner. With the thickness of the nanofilms decreasing, the surface effects on the elastic properties become dominant.

Molecular Dynamics Simulation Of Adhesion Forces In A Dipalmitoylphosphatidylcholine Membrane

Adhesion forces of Dipalmitoylphosphatidylcholine ( DPPC) membrane in the gel phase are investigated by molecular dynamics ( MD) simulation. In the simulations, individual DPPC molecules are pulled out of DPPC membranes with different rates and we get the maximum adhesion forces of DPPC membrane. We find that the maximum adhesion forces increase with pull rate, from about 400 to 700 pN when pull rates are from 0.001 to 0.03 nm/ps. We analyze the relationship between pull rate and adhesion forces of different origins using Brownian dynamics and notice that viscosity of solvent plays an important role in adhesion forces. Then we simulate the motion of a single DPPC molecule in solvent and it elucidates that the maximum drag force is almost linear with respect to the pull rate. We use Stokes' relation to describe the motion of a single DPPC molecule and deduce the effective length of a DPPC molecule. Conformational analyses indicate that the free energy variation of a DPPC molecule inside and outside of the DPPC membrane is an essential part of adhesion energy.

Bubble Dynamics In Nucleate Pool Boiling On Thin Wires In Microgravity

A temperature-controlled pool boiling (TCPB) device has been developed to study the bubble behavior and heat transfer in pool boiling phenomenon both in normal gravity and in microgravity. A thin platinum wire of 60 mu m in diameter and 30 mm in length is simultaneously used as heater and thermometer. The fluid is R113 at 0.1 MPa and subcooled by 26 degrees C nominally for all cases. Three modes of heat transfer, namely single-phase natural convection, nucleate boiling, and two-mode transition boiling, are observed in the experiment both in microgravity aboard the 22nd Chinese recoverable satellite and in normal gravity on the ground before and after the space flight. Dynamic behaviors of vapor bubbles observed in these experiments are reported and analyzed in the present paper. In the regime of fully developed nucleate boiling, the interface oscillation due to coalescence of adjacent tiny bubbles is the primary reason of the departure of bubbles in microgravity. On the contrary, in the discrete bubble regime, it's observed that there exist three critical bubble diameters in microgravity, dividing the whole range of the observed bubbles into four regimes. Firstly, tiny bubbles are continually forming and growing on the heating surface before departing slowly from the wire when their sizes exceed some value of the order of 10(-1) mm. The bigger bubbles with about several millimeters in diameter stay on the wire, oscillate along the wire, and coalesce with adjacent bubbles. The biggest bubble with diameter of the order of 10 mm, which was formed immediately after the onset of boiling, stays continuously

Strengthening Mechanisms And Dislocation Dynamics In Twinned Metal Nanowires

Microtwins are frequently observed in face-centered-cubic (fcc) metal nanowires with low stacking fault energy. The authors have previously reported that the tensile Yield strength of copper nanowires can be increased by, the presence of twin boundaries. lit this work, simulations are carried out under both uniaxial tension and compression loading, to demonstrate that the strengthening effects are inherent to these nanowires, independent of the loading condition (tensile/compressive). It appears that the strengthening mechanism of the twinned nanowires can be attributed to stress redistribution due to the change of crystallographic orientations across twin boundaries, which requires larger external stress to make them Yield as compared to the twin-free wire.

Molecular Dynamics Study Of Mechanical Behaviour Of Screw Dislocation During Cutting With Diamond Tip On Silicon

By means of Tersoff and Morse potentials, a three-dimensional molecular dynamics simulation is performed to study atomic force microscopy cutting on silicon monocrystal surface. The interatomic forces between the workpiece and the pin tool and the atoms of workpiece themselves are calculated. A screw dislocation is introduced into workpiece Si. It is found that motion of dislocations does not occur during the atomic force microscopy cutting processing. Simulation results show that the shear stress acting on dislocation is far below the yield strength of Si.

Large Scale Numerical Simulations For Multi-Phase Fluid Dynamics With Moving Interfaces

This short communication presents our recent studies to implement numerical simulations for multi-phase flows on top-ranked supercomputer systems with distributed memory architecture. The numerical model is designed so as to make full use of the capacity of the hardware. Satisfactory scalability in terms of both the parallel speed-up rate and the size of the problem has been obtained on two high rank systems with massively parallel processors, the Earth Simulator (Earth simulator research center, Yokohama Kanagawa, Japan) and the TSUBAME (Tokyo Institute of Technology, Tokyo, Japan) supercomputers.

Grain Boundary Effects On Plastic Deformation And Fracture Mechanisms In Cu Nanowires: Molecular Dynamics Simulations

Metallic nanowires have many attractive properties such as ultra-high yield strength and large tensile elongation. However, recent experiments show that metallic nanowires often contain grain boundaries, which are expected to significantly affect mechanical properties. By using molecular dynamics simulations, here, we demonstrate that polycrystalline Cu nanowires exhibit tensile deformation behavior distinctly different from their single-crystal counterparts. A significantly lowered yield strength was observed as a result of dislocation emission from grain boundaries rather than from free surfaces, despite of the very high surface to volume ratio. Necking starts from the grain boundary followed by fracture, resulting in reduced tensile ductility. The high stresses found in the grain boundary region clearly play a dominant role in controlling both inelastic deformation and fracture processes in nanoscale objects. These findings have implications for designing stronger and more ductile structures and devices on nanoscale.

Molecular Dynamics Study Of The Mechanical Behavior Of Nickel Nanowire: Strain Rate Effects

We present the analysis of uniaxial deformation of nickel nanowires using molecular dynamics simulations, and address the strain rate effects on mechanical responses and deformation behavior. The applied strain rate is ranging from 1 x 10(8) s(-1) to 1.4 x 10(11) s(-1). The results show that two critical strain rates, i.e., 5 x 10(9) s(-1) and 8 x 10(10) s(-1), are observed to play a pivotal role in switching between plastic deformation modes. At strain rate below 5 x 10(9) s(-1), Ni nanowire maintains its crystalline structure with neck occurring at the end of loading, and the plastic deformation is characterized by {111} slippages associated with Shockley partial dislocations and rearrangements of atoms close to necking region. At strain rate above 8x10(10) s(-1), Ni nanowire transforms from a fcc crystal into a completely amorphous state once beyond the yield point, and hereafter it deforms uniformly without obvious necking until the end of simulation. For strain rate between 5 x 10(9) s(-1) and 8 x 10(10) s(-1), only part of the nanowire exhibits amorphous state after yielding while the other part remains crystalline state. Both the {111} slippages in ordered region and homogenous deformation in amorphous region contribute to the plastic deformation. (C) 2007 Published by Elsevier B.V.