Multiscale modeling of the dynamics of solids at finite temperature

We develop a general multiscale method for coupling atomistic and continuum simulations using the framework of the heterogeneous multiscale method (HMM). Both the atomistic and the continuum models are formulated in the form of conservation laws of mass, momentum and energy. A macroscale solver, here the finite volume scheme, is used everywhere on a macrogrid; whenever necessary the macroscale fluxes are computed using the microscale model, which is in turn constrained by the local macrostate of the system, e.g. the deformation gradient tensor, the mean velocity and the local temperature. We discuss how these constraints can be imposed in the form of boundary conditions. When isolated defects are present, we develop an additional strategy for defect tracking. This method naturally decouples the atomistic time scales from the continuum time scale. Applications to shock propagation, thermal expansion, phase boundary and twin boundary dynamics are presented.     

Effective particle size from molecular dynamics simulations in fluids

We report molecular dynamics simulations designed to investigate the effective size of colloidal particles suspended in a fluid in the vicinity of a rigid wall where all interactions are defined by smooth atomic potential functions. These simulations are used to assess how the behavior of this system at the atomistic length scale compares to continuum mechanics models. In order to determine the effective size of the particles, we calculate the solvent forces on spherical particles of different radii as a function of different positions near and overlapping with the atomistically defined wall and compare them to continuum models. This procedure also then determines the effective position of the wall. Our analysis is based solely on forces that the particles sense, ensuring self-consistency of the method. The simulations were carried out using both Weeks–Chandler–Andersen and modified Lennard-Jones (LJ) potentials to identify the different contributions of simple repulsion and van der Waals attractive forces. Upon correction for behavior arising the discreteness of the atomic system, the underlying continuum physics analysis appeared to be correct down to much less than the particle radius. For both particle types, the effective radius was found to be \(\sim 0.75\sigma \), where \(\sigma \) defines the length scale of the force interaction (the LJ diameter). The effective “hydrodynamic” radii determined by this means are distinct from commonly assumed values of \(0.5\sigma \) and \(1.0\sigma \), but agree with a value developed from the atomistic analysis of the viscosity of such systems.     

An atomistic–continuum hybrid scheme for numerical simulation of droplet spreading on a solid surface

We present an atomistic–continuum hybrid method to investigate spreading dynamics of drops on solid surfaces. The Navier–Stokes equations are solved by the finite-volume method in a continuum domain comprised of the main body of the drop, and atomistic molecular dynamics simulations are used in a particle domain in the vicinity of the contact line. The spatial coupling between the continuum and particle domains is achieved through constrained dynamics of flux continuities in an overlap domain.     

Multiscale plasticity modeling: coupled atomistics and discrete dislocation mechanics

A computational method (CADD) is presented whereby a continuum region containing dislocation defects is coupled to a fully atomistic region. The model is related to previous hybrid models in which continuum finite elements are coupled to a fully atomistic region, with two key advantages: the ability to accomodate discrete dislocations in the continuum region and an algorithm for automatically detecting dislocations as they move from the atomistic region to the continuum region and then correctly “converting” the atomistic dislocations into discrete dislocations, or vice-versa. The resulting CADD model allows for the study of 2d problems involving large numbers of defects where the system size is too big for fully atomistic simulation, and improves upon existing discrete dislocation techniques by preserving accurate atomistic details of dislocation nucleation and other atomic scale phenomena. Applications to nanoindentation, atomic scale void growth under tensile stress, and fracture are used to validate and demonstrate the capabilities of the model.     

An atomistic-based continuum theory for carbon nanotubes: analysis of fracture nucleation

Carbon nanotubes (CNTs) display unique properties and have many potential applications. Prior theoretical studies on CNTs are based on atomistic models such as empirical potential molecular dynamics (MD), tight-binding methods, or first-principles calculations. Here we develop an atomistic-based continuum theory for CNTs. The interatomic potential is directly incorporated into the continuum analysis through constitutive models. Such an approach involves no additional parameter fitting beyond those introduced in the interatomic potential. The atomistic-based continuum theory is then applied to study fracture nucleation in CNTs by modelling it as a bifurcation problem. The results agree well with the MD simulations.     

Nanoscale finite element models for vibrations of single-walled carbon nanotubes: atomistic versus continuum

By the atomistic and continuum finite element models, the free vibration behavior of single-walled carbon nanotubes （SWCNTs） is studied. In the atomistic finite element model, the bonds and atoms are modeled by the beam and point mass elements, respectively. The molecular mechanics is linked to structural mechanics to determine the elastic properties of the mentioned beam elements. In the continuum finite element approach, by neglecting the discrete nature of the atomic structure of the nanotubes, they are modeled with shell elements. By both models, the natural frequencies of SWCNTs are computed, and the effects of the geometrical parameters, the atomic structure, and the boundary conditions are investigated. The accuracy of the utilized methods is verified in comparison with molecular dynamic simulations. The molecular structural model leads to more reliable results, especially for lower aspect ratios. The present analysis provides valuable information about application of continuum models in the investigation of the mechanical behaviors of nanotubes.     

An atomistically validated continuum model for strain relaxation and misfit dislocation formation

In this paper, molecular dynamics (MD) calculations have been used to examine the physics behind continuum models of misfit dislocation formation and to assess the limitations and consequences of approximations made within these models. Without compromising the physics of misfit dislocations below a surface, our MD calculations consider arrays of dislocation dipoles constituting a mirror imaged “surface”. This allows use of periodic boundary conditions to create a direct correspondence between atomistic and continuum representations of dislocations, which would be difficult to achieve with free surfaces. Additionally, by using long-time averages of system properties, we have essentially reduced the errors of atomistic simulations of large systems to “zero”. This enables us to deterministically compare atomistic and continuum calculations. Our work results in a robust approach that uses atomistic simulation to accurately calculate dislocation core radius and energy without the continuum boundary conditions typically assumed in the past, and the novel insight that continuum misfit dislocation models can be inaccurate when incorrect definitions of dislocation spacing and Burgers vector in lattice-mismatched systems are used. We show that when these insights are properly incorporated into the continuum model, the resulting energy density expression of the lattice-mismatched systems is essentially indistinguishable from the MD results.     

Evaluation of coarse-grained CFD-DEM models with the validation of PEPT measurements

Computational Fluid Dynamics coupled with Discrete Element Method (CFD-DEM) is a commonly used numerical method to model gas-solid flow in fluidised beds and other multiphase systems. A significant limitation of CFD-DEM is the feasibility of the realistic simulation of large numbers of particles. Coarse-graining (CG) approaches, through which groups of multiple individual particles are represented by single, larger particles, can substantially reduce the total number of particles while maintaining similar system dynamics. As these three CG models have not previously been compared, there remains some debate, however, about the best practice in the application of CG in CFD-DEM simulations. In this paper, we evaluate the performance of three typical CG methods based on simulations of a bubbling fluidised bed. This is achieved through the use of a numerical validation framework, which makes full use of the high-resolution 3D positron emission particle tracking (PEPT) measurements to rigorously validate the outputs of CFD-DEM simulations conducted using various different coarse-graining models, and various different degrees of coarse-graining. The particle flow behaviours in terms of the particle occupancy field, velocity field, circulation time, and bubble size and velocity, are comprehensively analysed. It is shown that the CG simulation starts to fail when the size ratio between the bed chamber and the particles decreases to approximately 20. It is also observed, somewhat surprisingly, that the specific CG approach applied to interparticle contact parameters does not have a substantial effect on the simulation results for the bubbling bed simulations across a wide range of CG factors.     

A multiscale method for dislocation nucleation and seamlessly passing scale boundaries

In the concurrent multiscale analysis, it is difficult to have truly seamless transition between the atomistic and continuum scale. This situation is even worse when defects pass through the boundary between different scales. For example, there is a lack of effective methods to handle the dislocation passing through scale boundaries which is important to investigate plasticity at the nanoscale. In this paper, the generalized particle (GP) method proposed by the first author is further developed so that a seamless transition and dislocation passing between different scales can be realized. Specifically, the linkage between different scales is through material neighbor-link cells (NLC) with scale duality. This indicates that material elements can be high-scale particles through a lumping process and can also be atoms via decomposition depending on the needs of the simulation. At the interface, the information transfer from bottom scale-up or from top scale-down is through the particles or atoms in the NLC. They are with the same material structure, all possess nonlocal constitutive behavior; thus, the smooth transition at the interface between different scales can be attained and validated to avoid non-physical responses. To save degrees of freedom, atoms are lumped together into a generalized particle in the domain in which the deformation gradient is near homogeneous. On the other hand, when defects such as dislocations in the atomistic domain are near the particle domain, the particles along dislocation propagation path and its surrounding region will be decomposed into atoms so dislocations can freely pass through the scale boundary and propagate inside the model just as it propagates in the deformed atomistic crystal structure. The method is verified first for seamless transition of variables at the scale boundary by a one-dimensional model and then verified for dislocation nucleation and propagation passing through scale boundaries in two cases, one is near the free surface and the other is inside of the copper nanowire. All the validations are through comparisons with fully atomistic analyses under same conditions. The comparison is satisfactory.     

Continuum Mechanics Modeling and Simulation of Carbon Nanotubes

The understanding of the mechanics of atomistic systems greatly benefits from continuum mechanics. One appealing approach aims at deductively constructing continuum theories starting from models of the interatomic interactions. This viewpoint has become extremely popular with the quasicontinuum method. The application of these ideas to carbon nanotubes presents a peculiarity with respect to usual crystalline materials: their structure relies on a two-dimensional curved lattice. This renders the cornerstone of crystal elasticity, the Cauchy–Born rule, insufficient to describe the effect of curvature. We discuss the application of a theory which corrects this deficiency to the mechanics of carbon nanotubes (CNTs). We review recent developments of this theory, which include the study of the convergence characteristics of the proposed continuum models to the parent atomistic models, as well as large scale simulations based on this theory. The latter have unveiled the complex nonlinear elastic response of thick multiwalled carbon nanotubes (MWCNTs), with an anomalous elastic regime following an almost absent harmonic range.     

Atomistic and continuum studies of stress and strain fields near a rapidly propagating crack in a harmonic lattice

Large-scale atomistic simulations of a mode I crack propagating in a harmonic lattice are presented. The objective of this work is to study the stress and strain fields near a rapidly propagating mode I crack. The asymptotic continuum mechanics solutions of the elastic fields are compared quantitatively with molecular-dynamics simulation results for different crack velocities. It is observed that both atomistic stress and atomistic strain can be successfully related to the corresponding continuum quantities. The study reveals that the atomistic simulation results agree well with the continuum theory predictions, which suggests that the continuum theory can be applied for nano-scale dynamic problems.     

Estimation of anisotropic elastic properties of nanocomposites using atomistic-continuum interphase model

We have revised classical micromechanics by accounting for the effect of interface to predict the effective anisotropic elastic properties of heterogeneous materials containing nano-inhomogeneities. In contrast to sharp interface between the matrix and inhomogeneity, we introduce the concept of interphase to account for the interfacial-stress effect at the nano-scale. The interphase’s constitutive properties are derived from atomistic simulations and then incorporated in a micromechanics-based interphase model to compute effective properties of nanocomposites. This scale transition approach bridges the gap between discrete atomic level interactions and continuum mechanics. An advantage of this approach is that it combines atomistic with continuum models that consider inhomogeneity and interphase morphology. It thereby enables us to account simultaneously for both the shape and the anisotropy of a nano-inhomogeneity and interphase at the continuum level when we compute material’s overall properties. In so doing, it frees us from making any assumptions about the interface characteristics between matrix and the nano-inhomogeneity.     

Assessing continuum postulates in simulations of granular flow

Continuum mechanics relies on the fundamental notion of a mesoscopic volume “element” in which properties averaged over discrete particles obey deterministic relationships. Recent work on granular materials suggests that a continuum law may be inapplicable, revealing inhomogeneities at the particle level, such as force chains and slow cage breaking. Here, we analyze large-scale three-dimensional discrete-element method (DEM) simulations of different granular flows and show that an approximate “granular element” defined at the scale of observed dynamical correlations (roughly three to five particle diameters) has a reasonable continuum interpretation. By viewing all the simulations as an ensemble of granular elements which deform and move with the flow, we can track material evolution at a local level. Our results confirm some of the hypotheses of classical plasticity theory while contradicting others and suggest a subtle physical picture of granular failure, combining liquid-like dependence on deformation rate and solid-like dependence on strain. Our computational methods and results can be used to guide the development of more realistic continuum models, based on observed local relationships between average variables.     

Multiscale modeling of crack initiation and propagation at the nanoscale

Fracture occurs on multiple interacting length scales; atoms separate on the atomic scale while plasticity develops on the microscale. A dynamic multiscale approach (CADD: coupled atomistics and discrete dislocations) is employed to investigate an edge-cracked specimen of single-crystal nickel, Ni, (brittle failure) and aluminum, Al, (ductile failure) subjected to mode-I loading. The dynamic model couples continuum finite elements to a fully atomistic region, with key advantages such as the ability to accommodate discrete dislocations in the continuum region and an algorithm for automatically detecting dislocations as they move from the atomistic region to the continuum region and then correctly “converting” the atomistic dislocations into discrete dislocations, or vice-versa. An ad hoc computational technique is also applied to dissipate localized waves formed during crack advance in the atomistic zone, whereby an embedded damping zone at the atomistic/continuum interface effectively eliminates the spurious reflection of high-frequency phonons, while allowing low-frequency phonons to pass into the continuum region.The simulations accurately capture the essential physics of the crack propagation in a Ni specimen at different temperatures, including the formation of nano-voids and the sudden acceleration of the crack tip to a velocity close to the material Rayleigh wave speed. The nanoscale brittle fracture happens through the crack growth in the form of nano-void nucleation, growth and coalescence ahead of the crack tip, and as such resembles fracture at the microscale. When the crack tip behaves in a ductile manner, the crack does not advance rapidly after the pre-opening process but is blunted by dislocation generation from its tip. The effect of temperature on crack speed is found to be perceptible in both ductile and brittle specimens.     

Multiscale boundary conditions in crystalline solids: Theory and application to nanoindentation

This paper presents a systematic approach to treating the interfaces between the localized (fine grain) and peripheral (coarse grain) domains in atomic scale simulations of crystalline solids. Based on Fourier analysis of regular lattices structures, this approach allows elimination of unnecessary atomic degrees of freedom over the coarse grain, without involving an explicit continuum model for the latter. The mathematical formulation involves compact convolution operators that relate displacements of the interface atoms and the adjacent atoms on the coarse grain. These operators are defined by geometry of the lattice structure, and interatomic potentials. Application and performance are illustrated on quasistatic nanoindentation simulations with a crystalline gold substrate. Complete atomistic resolution on the coarse grain is alternatively employed to give the benchmark solutions. The results are found to match well for the multiscale and the full atomistic simulations.     

Computational nanoscale plasticity simulations using embedded atom potentials

In determining structure–property relations for plasticity at different size scales, it is desired to bridge concepts from the continuum to the atom. This raises many questions related to volume averaging, appropriate length scales of focus for an analysis, and postulates in continuum mechanics. In a preliminary effort to evaluate bridging size scales and continuum concepts with descritized phenomena, simple shear molecular dynamics simulations using the Embedded Atom Method (EAM) potentials were performed on single crystals. In order to help evaluate the continuum quantities related to the kinematic and thermodynamic force variables, finite element simulations (with different material models) and macroscale experiments were also performed. In this scoping study, various parametric effects on the stress state and kinematics have been quantified. The parameters included the following: crystal orientation (single slip, double slip, quadruple slip, octal slip), temperature (300 and 500 K), applied strain rate (10⁶–10¹² s⁻¹), specimen size (10 atoms to 2 μm), specimen aspect ratio size (1:8–8:1), deformation path (compression, tension, simple shear, and torsion), and material (nickel, aluminum, and copper). Although many conclusions can be drawn from this work, which has provided fodder for more studies, several major conclusions can be drawn.

• The yield stress is a function of a size scale parameter (volume-per-surface area) that was determined from atomistic simulations coupled with experiments. As the size decreases, the yield stress increases.

• Although the thermodynamic force (stress) varies at different size scales, the kinematics of deformation appears to be very similar based on atomistic simulations, finite element simulations, and physical experiments.

Atomistic simulations, that inherently include extreme strain rates and size scales, give results that agree with the phenomenological attributes of plasticity observed in macroscale experiments. These include strain rate dependence of the flow stress into a rate independent regime; approximate Schmid type behavior; size scale dependence on the flow stress, and kinematic behavior of large deformation plasticity.     

Coarse-grained atomistic simulation of dislocations

This paper presents a new methodology for coarse-grained atomistic simulation of dislocation dynamics. The methodology combines an atomistic formulation of balance equations and a modified finite element method employing rhombohedral-shaped 3D solid elements suitable for fcc crystals. With significantly less degrees of freedom than that of a fully atomistic model and without additional constitutive rules to govern dislocation activities, this new coarse-graining (CG) method is shown to be able to reproduce key phenomena of dislocation dynamics for fcc crystals, including dislocation nucleation and migration, formation of stacking faults and Lomer-Cottrell locks, and splitting of stacking faults, all comparable with fully resolved molecular dynamics simulations. Using a uniform coarse mesh, the CG method is then applied to simulate an initially dislocation-free submicron-sized thin Cu sheet. The results show that the CG simulation has captured the nucleation and migration of large number of dislocations, formation of multiple stacking fault ribbons, and the occurrence of complex dislocation phenomena such as dislocation annihilation, cutting, and passing through the stacking faults. The distinctions of this method from existing coarse-graining or multiscale methods and its potential applications and limitations are also discussed.     

COARSE-GRAINED ATOMISTIC MODELING AND SIMULATION OF INELASTIC MATERIAL BEHAVIOR

This paper presents a new methodology for coarse-grained atomistic simulation of inelastic material behavior including phase transformations in ceramics and dislocation mediated plasticity in metals. The methodology combines an atomistic formulation of balance equations and a modified finite element method. With significantly fewer degrees of freedom than those of a fully atomistic model and without additional constitutive rules but the interatomic force field, the new coarse-grained (CG) method is shown to be feasible in predicting the nonlinear constitutive responses of materials and also reproducing atomic-scale phenomena such as phase transformations (diamond →β-Sn) in silicon and dislocation nucleation and migration, formation of dislocation loops and stacking faults ribbons in single crystal nickel. Direct comparisons between CG and the corresponding full molecular dynamics (MD) simulations show that the present methodology is efficient and promising in modeling and simulation of inelastic material behavior without losing the essential atomistic features. The potential applications and the limitations of the CG method are also discussed.     

Surface effects in non-uniform nanobeams: Continuum vs. atomistic modeling

Nanobeams are expected to be one of the key structural elements in nanotechnology. Contrary to macroscopic structures, surface effects can strongly influence the stress and deformation properties of nano-devices. In addition, at such small scales, material non-uniformity becomes significant and must be considered.In this work, a continuum model for nanobeams, including both surface effects and material heterogeneity is developed. The model treats the surfaces as separate material layers with finite thickness. The continuum solution is compared with atomistic simulations, from which the effective bulk and surface properties are computed independently. A special case of self-deflection due to surface non-uniformities, which is important for design of nanosensors, is studied. A comparison between continuum and atomistic solutions reveals differences, which originate from local transition effects in the neighborhood of strong non-uniformities. This discrepancy is addressed by defining an effective length, found by correlating the beam deflections from both methods.