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.

Study Of Mechanical Properties Of Amorphous Copper With Molecular Dynamics Simulation

The formation and mechanical properties of amorphous copper are studied using molecular dynamics simulation. The simulations of tension and shearing show that more pronounced plasticity is found under shearing, compared to tension. Apparent strain hardening and strain rate effect are observed. Interestingly, the variations of number density of atoms during deformation indicate free volume creation, especially under higher strain rate. In particular, it is found that shear induced dilatation does appear in the amorphous metal.

Molecular dynamics simulations of DFZ

Dislocation emission from the crack tip in copper under mode II loading is simulated with molecular dynamics method. After 26 partial dislocations are emitted and then relaxed to reach the equilibrium under the constant displacement, the double pile-ups (including an inverse pile-up and a pile-up) are formed. i.e., the first dislocation is piled up before the obstruction, and the last dislocation is piled up ahead of the crack tip. These results conform to the TEM observations.

Molecular dynamics studies on dislocations in crystallites of nanocrystalline alpha-iron

Molecular dynamics simulations have been carried our to study the atomic structure of the crystalline component of nanocrystalline alpha-iron. A two-dimensional computational block is used to simulate the consolidation process. It is found that dislocations are generated in the crystallites during consolidation when the grain size is large enough. The critical value of the grain size for dislocation generation appears to be about 9 nm. This result agrees with experiment qualitatively. AN dislocations that are preset in the original grains glide out during consolidation. It shows that dislocations in the crystallites we generated in consolidation process, but not in the original grains. Higher consolidation pressure results in more dislocations. Furthermore, new interfaces are found within crystallites. These interfaces might result from the special environment of nanomaterial. (C) 1998 Acta Metallurgica Inc.

Control dynamics of severe acute respiratory syndrome transmission

Severe acute respiratory syndrome (SARS) is a serious disease with many puzzling features. We present a simple, dynamic model to assess the epidemic potential of SARS and the effectiveness of control measures. With this model, we analysed the SARS epidemic data in Beijing. The data fitting gives the basic case reproduction number of 2.16 leading to the outbreak, and the variation of the effective reproduction number reflecting the control effect. Noticeably, our study shows that the response time and the strength of control measures have significant effects on the scale of the outbreak and the lasting time of the epidemic.

Streamline topology and dilute particle dynamics in a Karman vortex street flow

Three types of streamline topology in a Karman vortex street flow are shown under the variation of spatial parameters. For the motion of dilute particles in the Karman vortex street flow, there exist a route of bifurcation to a chaotic orbit and more attractors in a bifurcation diagram for the proportion of particle density to fluid density. Along with the increase of spatial parameters in the flow field, the bifurcation process is suspended, as well as more and more attractors emerge. In the motion of dilute particles, a drag term and gravity term dominate and result in the bifurcation phenomenon.

Size effects on the melting of nickel nanowires: a molecular dynamics study

The melting process of nickel nanowires are simulated by using molecular dynamics with the quantum Sutten-Chen many-body force field. The wires studied were approximately cylindrical in cross-section and periodic boundary conditions were applied along their length; the atoms were arranged initially in a face-centred cubic structure with the [0 0 1] direction parallel to the long axis of the wire. The size effects of the nanowires on the melting temperatures are investigated. We find that for the nanoscale regime, the melting temperatures of Ni nanowires are much lower than that of the bulk and are linear with the reciprocal of the diameter of the nanowire. When a nanowire is heated up above the melting temperature, the neck of the nanowire begins to arise and the diameter of neck decreases rapidly with the equilibrated running time. Finally, the breaking of nanowire arises, which leads to the formation of the spherical clusters. (C) 2004 Elsevier B.V. All rights reserved.

Molecular dynamics simulation of plastic deformation of nanotwinned copper

The plastic deformation of polycrystalline Cu with ultrathin lamella twins has been studied using molecular dynamics simulations. The results of uniaxial tensile deformation simulation show that the abundance of twin boundaries provides obstacles to dislocation motion, which in consequence leads to a high strain hardening rate in the nanotwinned Cu. We also show that the twin lamellar spacing plays a vital role in controlling the strengthening effects, i.e., the thinner the thickness of the twin lamella, the harder the material. Additionally, twin boundaries can act as dislocation nucleation sites as they gradually lose coherency at large strain. These results indicate that controlled introduction of nanosized twins into metals can be an effective way of improving strength without suppression tensile ductility. (C) 2007 American Institute of Physics.

Nonlinear dynamics of atomic force microscopy with intermittent contact

When the atomic force microscopy (AFM) in tapping mode is in intermittent contact with a soft substrate, the contact time can be a significant portion of a cycle, resulting in invalidity of the impact oscillator model, where the contact time is assumed to be infinitely small. Furthermore, we demonstrate that the AFM intermittent contact with soft substrate can induce the motion of higher modes in the AFM dynamic response. Traditional ways of modeling AFM (one degree of freedom (DOF) system or single mode analysis) are shown to have serious mistakes when applied to this kind of problem. A more reasonable displacement criterion on contact is proposed, where the contact time is a function of the mechanical properties of AFM and substrate, driving frequencies/amplitude, initial conditions, etc. Multi-modal analysis is presented and mode coupling is also shown. (c) 2006 Published by Elsevier Ltd.

Adsorption of His-tagged peptide to Ni, Cu and Au (100) surfaces: Molecular dynamics simulation

Molecular dynamics (MD) simulations are performed to study the interaction of His-tagged peptide with three different metal surfaces in explicit water. The equilibrium properties are analyzed by using pair correlation functions (PCF) to give an insight into the behavior of the peptide adsorption to metal surfaces in water solvent. The intermolecular interactions between peptide residues and the metal surfaces are evaluated. By pulling the peptide away from the peptide in the presence of solvent water, peeling forces are obtained and reveal the binding strength of peptide adsorption on nickel, copper and gold. From the analysis of the dynamics properties of the peptide interaction with the metal surfaces, it is shown that the affinity of peptide to Ni surface is the strongest, while on Cu and An the affinity is a little weaker. In MD simulations including metals, the His-tagged region interacts with the substrate to an extent greater than the other regions. The work presented here reveals various interactions between His-tagged peptide and Ni/Cu/Au surfaces. The interesting affinities and dynamical properties of the peptide are also derived. The results give predictions for the structure of His-tagged peptide adsorbing on three different metal surfaces and show the different affinities between them, which assist the understanding of how peptides behave on metal surfaces and of how designers select amino sequences in molecule devices design. (c) 2007 Elsevier Ltd. All rights reserved.

Computer simulation on the collision-sticking dynamics of two colloidal particles in an optical trap

Collisions of a particle pair induced by optical tweezers have been employed to study colloidal stability. In order to deepen insights regarding the collision-sticking dynamics of a particle pair in the optical trap that were observed in experimental approaches at the particle level, the authors carry out a Brownian dynamics simulation. In the simulation, various contributing factors, including the Derjaguin-Landau-Verwey-Overbeek interaction of particles, hydrodynamic interactions, optical trapping forces on the two particles, and the Brownian motion, were all taken into account. The simulation reproduces the tendencies of the accumulated sticking probability during the trapping duration for the trapped particle pair described in our previous study and provides an explanation for why the two entangled particles in the trap experience two different statuses. (c) 2007 American Institute of Physics.

Vortex dynamics in the studies of looping in tropical cyclone tracks

This paper proposes criteria for predicting the tendency of looping in tropical cyclone tracks using the approach of vortex dynamics. We model the asymmetric structure of a cyclone by a system of vortex patches. The evolution of such system of vortices is simulated by the method of contour dynamics. A new set of exact analytic formulas for contour dynamics calculations is derived, which is shown to be more computationally effective. Based on point-vortex models, we derive analytic formulas for the criteria of looping in a cyclone track. From numerical experiments, the simulated trajectories obtained from the point-vortex system and vortex patch system agree quite well. Hence, the looping criteria obtained from the point-vortex system can be applied by forecasters to stay alert for tendency of looping in a cyclone track. To demonstrate the applicability of the proposed criteria, the trajectory of Typhoon Yancy (9012), whose field data are available from ''TCM-90'', is simulated. The case study shows that the asymmetric structure similar to the pattern of a beta gyre is responsible for its recurvature when Yancy landed Fujian Province, China on 20 August 1990.

Dislocations in the crystallites of nanocrystalline alpha-Fe - A molecular dynamics study

Molecular dynamics simulations are carried out in order to study the atomic structure of crystalline component, of nanocrystalline alpha-Fe when it is consolidated from small grains. A two-dimensional computational block is used to simulate the consolidation process. All the preset dislocations in the original grains glide out of them in the consolidation process, but new dislocations can generate when the grain size is large enough. It shows that dislocations exist in the consolidated material rather than in the original grains. Whether dislocations exist in the crystalline component of the resultant model nana-material depends upon grain size. The critical value of grain size for dislocation generation appears to be about 9 nm. This result agrees with experiments qualitatively.