MOLECULAR DYNAMICS SIMULATIONS OF GRANULAR MATERIALS

1993 ◽  
Vol 04 (02) ◽  
pp. 309-316 ◽  
Author(s):  
HANS J. HERRMANN

When sand or other granular materials are shaken, poured or sheared many intriguing phenomena can be observed. We will model the granular medium by a packing of elastic spheres and simulate it via Molecular Dynamics. Dissipation of energy and shear friction at collisions are included. The onset of fluidization can be determined and is in good agreement with experiments. On a vibrating plate we observe the formation of convection cells due to walls or amplitude modulations. Density and velocity profiles on conveyor belts are measured and the influence of an obstacle discussed. We mention various types of rheology for flow down an inclined chute or through a pipe and outflowing containers.

1993 ◽  
Vol 07 (09n10) ◽  
pp. 1779-1788 ◽  
Author(s):  
JASON A.C. GALLAS ◽  
HANS J. HERRMANN ◽  
STEFAN SOKOLOWSKI

When sand or other granular materials are shaken, poured or sheared many intriguing phenomena can be observed. We will model the granular medium by a packing of elastic spheres and simulate it via Molecular Dynamics. Dissipation of energy and shear friction at collisions are included. The onset of fluidization can be determined and is in good agreement with experiments. On a vibrating plate we observe the formation of convection cells due to walls or amplitude modulations. Density and velocity profiles on conveyor belts are measured and the influence of an obstacle discussed. We mention various types of rheology for flow down an inclined chute or through a pipe and outflowing containers.


2009 ◽  
Vol 60-61 ◽  
pp. 315-319 ◽  
Author(s):  
W.W. Zhang ◽  
Qing An Huang ◽  
H. Yu ◽  
L.B. Lu

Molecular dynamics simulations are carried out to characterize the mechanical properties of [001] and [110] oriented silicon nanowires, with the thickness ranging from 1.05nm to 3.24 nm. The nanowires are taken to have ideal surfaces and (2×1) reconstructed surfaces, respectively. A series of simulations for square cross-section Si nanowires have been performed and Young’s modulus is calculated from energy–strain relationship. The results show that the elasticity of Si nanowires is strongly depended on size and surface reconstruction. Furthermore, the physical origin of above results is analyzed, consistent with the bond loss and saturation concept. The results obtained from the molecular dynamics simulations are in good agreement with the values of first-principles. The molecular dynamics simulations combine the accuracy and efficiency.


2012 ◽  
Vol 501 ◽  
pp. 64-69 ◽  
Author(s):  
Yan He ◽  
Yuan Zheng Tang ◽  
Man Ding ◽  
Lian Xiang Ma

Normal thermal conductivity of amorphous and crystalline SiO2nano-films is calculated by nonequilibrium molecular dynamics (NEMD) simulations in the temperature range from 100 to 700K and thicknesses from 2 to 6nm. The calculated temperature and thickness dependences of thermal conductivity are in good agreement with previous literatures. In the same thickness, higher thermal conductivity is obtained for crystalline SiO2nano-films. And more importantly, for amorphous SiO2nano-films, thickness can be any direction of x, y, z-axis without effect on the normal thermal conductivity, for crystalline SiO2nano-films, the different thickness directions obtain different thermal conductivity results. The different results of amorphous and crystalline SiO2nano-films simply show that film thickness and grain morphology will cause different effects on thermal conductivity.


Author(s):  
Kun Xie ◽  
Chong Qiao ◽  
Hong Shen ◽  
Riyi Yang ◽  
Ming Xu ◽  
...  

Abstract Zr-Rh metallic glass has enabled its many applications in vehicle parts, sports equipment and so on due to its outstanding performance in mechanical property, but the knowledge of the microstructure determining the superb mechanical property remains yet insufficient. Here, we develop a deep neural network potential of Zr-Rh system by using machine learning, which breaks the dilemma between the accuracy and efficiency in molecular dynamics simulations, and greatly improves the simulation scale in both space and time. The results show that the structural features obtained from the neural network method are in good agreement with the cases in ab initio molecular dynamics simulations. Furthermore, we build a large model of 5400 atoms to explore the influences of simulated size and cooling rate on the melt-quenching process of Zr77Rh23. Our study lays a foundation for exploring the complex structures in amorphous Zr77Rh23, which is of great significance for the design and practical application.


2021 ◽  
Author(s):  
JOSH KEMPPAINEN ◽  
IVAN GALLEGOS ◽  
PRATHAMESH DESHPANDE ◽  
JACOB GISSINGER ◽  
GREGORY ODEGARD

Furan resins can be used as precursor resin for Carbon-Carbon Composites but has also been used in adhesives, acid/corrosion resistant materials, and as an alternative fuel precursor [15]. This paper contains the most current understanding of the structure of furan resin and a Molecular Dynamics workflow for computationally simulating its polymerization with the 'fix bond/react' command implemented in LAMMPS. The predicted mechanical properties of the polymerized resin are in good agreement with the literature values.


We report molecular dynamics simulations of n-butane and n-hexane adsorbed in the zeolite silicalite. These systems have been modelled within a rigid framework approximation and a Ryckaert-Bellemans model was adopted to describe the adsorbed molecules. The parametrization due to June, Bell and Theodorou has been used to describe the host-guest molecule interactions. Long simulations (1000 ps) have been performed, modelling these systems at a variety of sorbate loadings and temperatures. These have allowed us to investigate how the presence of the zeolite framework affects the thermodynamic properties, confomational distributions and diffusion related properties of the adsorbed molecules, and their response to changes in the loading and temperature. We have obtained estimates of the diffusion coefficients and activation energies which are in good agreement with experimentally measured data.


2011 ◽  
Vol 133 (7) ◽  
Author(s):  
John C. Duda ◽  
Pamela M. Norris ◽  
Patrick E. Hopkins

We present a new model for predicting thermal boundary conductance in the classical limit. This model takes a different form than those of the traditionally used mismatch theories in the fact that the temperature dependence of thermal boundary conductance is driven by the phononic scattering mechanisms of the materials comprising the interface as opposed to the heat capacities of those materials. The model developed in this work assumes that a phonon on one side of an interface may not scatter at the interface itself but instead scatter with phonons in the adjacent material via the scattering processes intrinsic in the adjacent material. We find that this model is in good agreement with classical molecular dynamics simulations of phonon transport across a Si/Ge interface.


2001 ◽  
Vol 703 ◽  
Author(s):  
Laurent Van Brutzel ◽  
Cindy L. Rountree ◽  
Rajiv K. Kalia ◽  
Aiichiro Nakano ◽  
Priya Vashishta

ABSTRACTParallel molecular dynamics simulations are performed to investigate dynamic fracture in bulk and nanostructured silica glasses at room temperature and 1000 K. In bulk silica the crack front develops multiple branches and nanoscale pores open up ahead of the crack tip. Pores coalesce and then they merge with the advancing crack-front to cause cleavage fracture. The calculated fracture toughness is in good agreement with experiments. In nanostrucutred silica the crack-front meanders along intercluster boundaries, merging with nanoscale pores in these regions to cause intergranular fracture. The failure strain in nanostructured silica is significantly larger than in the bulk systems.


1992 ◽  
Vol 03 (06) ◽  
pp. 1281-1293 ◽  
Author(s):  
GERALD H. RISTOW

In this paper we present an efficient algorithm to perform Molecular Dynamics simulations on a distributed memory parallel computer, the Intel iPSC/860. The proposed model describes the flow properties of granular materials in two dimensions. The specific implementation on a 32 node iPSC/860, especially the message passing and load balancing algorithms, are discussed in detail. Performance data are shown for different computers and varying node numbers of the iPSC/860. As a physical example we calculate some properties of the outflow behavior from a two-dimensional hopper and we discuss possible extensions of our model to three dimensions. Our simulations show that Molecular Dynamics simulations can be implemented quite efficiently on a distributed memory parallel computer if one assures load balancing and optimizes the internode communications.


Sign in / Sign up

Export Citation Format

Share Document