Coupling of atomistic and continuum simulations using a bridging scale decomposition

We present a new method for coupling molecular dynamics (MD) and continuum mechanics simulations that is based on the projection of the MD solution onto the coarse scale shape functions. This projection, or "bridging scale", represents that part of the solution that is obtainable by both solution methods. By subtracting the bridging scale from the total solution, we arrive at a coarse-fine decomposition that, by a proper choice of projection operator, decouples the kinetic energy of the two simulations. The resulting decomposition can be used in a finite-temperature simulation method in which MD is used only in a localized region, while the continuum simulation covers the entire domain, including the MD region to which it is coupled. One major advantage of this approach is that separate time step sizes can be used in the two simulations, so that the coarse scale time step is not limited to the time scale of the atomic vibrations present in the fine scale. Example problems are demonstrated on a 1D lattice, for which the method is shown to be accurate both for harmonic and anharmonic interatomic potentials.

[1]  Ted Belytschko,et al.  Mixed-time implicit-explicit finite elements for transient analysis , 1982 .

[2]  J. D. Doll,et al.  Generalized Langevin equation approach for atom/solid‐surface scattering: Collinear atom/harmonic chain model , 1974 .

[3]  W. E,et al.  Matching conditions in atomistic-continuum modeling of materials. , 2001, Physical review letters.

[4]  Vivek B. Shenoy,et al.  Finite Temperature Quasicontinuum Methods , 1998 .

[5]  M. Ortiz,et al.  Quasicontinuum analysis of defects in solids , 1996 .

[6]  J. Q. Broughton,et al.  Concurrent coupling of length scales: Methodology and application , 1999 .

[7]  E. B. Tadmor,et al.  Quasicontinuum models of interfacial structure and deformation , 1998 .

[8]  Gregory J. Wagner,et al.  Hierarchical enrichment for bridging scales and mesh-free boundary conditions , 2001 .

[9]  Michael Ortiz,et al.  Mixed Atomistic and Continuum Models of Deformation in Solids , 1996 .

[10]  Michael Ortiz,et al.  Quasicontinuum models of fracture and plasticity , 1998 .

[11]  Steven J. Plimpton,et al.  Computational nanoscale plasticity simulations using embedded atom potentials , 2001 .

[12]  Noam Bernstein,et al.  Spanning the continuum to quantum length scales in a dynamic simulation of brittle fracture , 1998 .

[13]  J. M. Haile,et al.  Molecular dynamics simulation : elementary methods / J.M. Haile , 1992 .

[14]  T. Hughes,et al.  The variational multiscale method—a paradigm for computational mechanics , 1998 .

[15]  Michael Ortiz,et al.  Mixed Atomistic-Continuum Models of Material Behavior: The Art of Transcending Atomistics and Informing Continua , 2001 .

[16]  Thomas J. R. Hughes,et al.  A variational multiscale approach to strain localization – formulation for multidimensional problems , 2000 .

[17]  Wing Kam Liu,et al.  Meshfree and particle methods and their applications , 2002 .

[18]  Ted Belytschko,et al.  Mixed methods for time integration , 1979 .

[19]  F. Reif,et al.  Fundamentals of Statistical and Thermal Physics , 1965 .

[20]  Eduard G. Karpov,et al.  Molecular dynamics boundary conditions for regular crystal lattices , 2004 .

[21]  W. Cai,et al.  Minimizing boundary reflections in coupled-domain simulations. , 2000, Physical review letters.

[22]  C. Farhat,et al.  Two efficient staggered algorithms for the serial and parallel solution of three-dimensional nonlinear transient aeroelastic problems , 2000 .

[23]  Michael Ortiz,et al.  Quasicontinuum simulation of fracture at the atomic scale , 1998 .

[24]  J. D. Doll,et al.  Generalized Langevin equation approach for atom/solid-surface scattering: General formulation for classical scattering off harmonic solids , 1976 .

[25]  Peter Gumbsch,et al.  Atomistic Aspects of Brittle Fracture , 2000 .

[26]  Wing Kam Liu,et al.  Development of mixed time partition procedures for thermal analysis of structures , 1982 .

[27]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[28]  Weimin Han,et al.  Convergence analysis of a hierarchical enrichment of Dirichlet boundary conditions in a mesh‐free method , 2002 .

[29]  Martin H. Müser,et al.  Towards an atomistic understanding of solid friction by computer simulations , 2002 .

[30]  R. A. Uras,et al.  Enrichment of the Finite Element Method With the Reproducing Kernel Particle Method , 1995 .