Open Access

ARTICLE

Multiscale Simulation Using Generalized Interpolation Material Point (GIMP) Method and Molecular Dynamics (MD)¹

by J. Ma², H. Lu², B. Wang², R. Hornung³, A. Wissink³, R. Kom,uri²

1 This work was performed under the auspices of the U.S. Department of Energy by University of California, Lawrence Livermore National Laboratory under Contract W-7405-Eng-48. UCRL-JRNL-216232.
2 School of Mechanical and Aerospace Engineering, Oklahoma State University, Stillwater, OK 74078. *Correspondence author, Tel: 405-744-5900; Fax: 405-744-7873; e-mail: ranga@ceat.okstate.edu
3 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore, CA 94551.

Computer Modeling in Engineering & Sciences 2006, 14(2), 101-118. https://doi.org/10.3970/cmes.2006.014.101

Abstract

A new method for multiscale simulation bridging two scales, namely, the continuum scale using the generalized interpolation material point (GIMP) method and the atomistic scale using the molecular dynamics (MD), is presented and verified in 2D. The atomistic strain from the molecular dynamics simulation is determined through interpolation of the displacement field into an Eulerian background grid using the same generalized interpolation functions as that in the GIMP method. The atomistic strain is consistent with that determined from the virial theorem for interior points but provides more accurate values at the boundary of the MD region and in the transition region between MD-GIMP. A material point in the continuum is split into smaller material points using multi-level refinement until it has nearly reached the atom size to couple with atoms in the MD region. Consequently, coupling between GIMP and MD is achieved by using compatible deformation, force, and energy fields in the transition region between GIMP and MD. The coupling algorithm is implemented in the Structural Adaptive Mesh Refinement Application Infrastructure (SAMRAI) for parallel processing. Both mode I and mode II crack problems are simulated using the coupling algorithm. The stress field near the crack tip was verified by comparing the results from the coupled simulations with purely GIMP simulations of the same model. Coupled simulation results were also compared with pure MD simulation results. In both cases, a very good agreement was obtained.

---

Keywords

---