CaltechTHESIS
  A Caltech Library Service

A Fully-nonlocal Energy-based Formulation and High-performance Realization of the Quasicontinuum Method

Citation

Amelang, Jeffrey Scott (2016) A Fully-nonlocal Energy-based Formulation and High-performance Realization of the Quasicontinuum Method. Dissertation (Ph.D.), California Institute of Technology. doi:10.7907/Z9SB43PH. http://resolver.caltech.edu/CaltechTHESIS:09152015-212147583

Abstract

The quasicontinuum (QC) method was introduced to coarse-grain crystalline atomic ensembles in order to bridge the scales from individual atoms to the micro- and mesoscales. Though many QC formulations have been proposed with varying characteristics and capabilities, a crucial cornerstone of all QC techniques is the concept of summation rules, which attempt to efficiently approximate the total Hamiltonian of a crystalline atomic ensemble by a weighted sum over a small subset of atoms. In this work we propose a novel, fully-nonlocal, energy-based formulation of the QC method with support for legacy and new summation rules through a general energy-sampling scheme. Our formulation does not conceptually differentiate between atomistic and coarse-grained regions and thus allows for seamless bridging without domain-coupling interfaces. Within this structure, we introduce a new class of summation rules which leverage the affine kinematics of this QC formulation to most accurately integrate thermodynamic quantities of interest. By comparing this new class of summation rules to commonly-employed rules through analysis of energy and spurious force errors, we find that the new rules produce no residual or spurious force artifacts in the large-element limit under arbitrary affine deformation, while allowing us to seamlessly bridge to full atomistics. We verify that the new summation rules exhibit significantly smaller force artifacts and energy approximation errors than all comparable previous summation rules through a comprehensive suite of examples with spatially non-uniform QC discretizations in two and three dimensions. Due to the unique structure of these summation rules, we also use the new formulation to study scenarios with large regions of free surface, a class of problems previously out of reach of the QC method. Lastly, we present the key components of a high-performance, distributed-memory realization of the new method, including a novel algorithm for supporting unparalleled levels of deformation. Overall, this new formulation and implementation allows us to efficiently perform simulations containing an unprecedented number of degrees of freedom with low approximation error.

Item Type:Thesis (Dissertation (Ph.D.))
Subject Keywords:atomistic to continuum transition ; computational science ; multiscale modeling ; nanoscale mechanics
Degree Grantor:California Institute of Technology
Division:Engineering and Applied Science
Major Option:Mechanical Engineering
Awards:Charles D. Babcock Award, 2012
Thesis Availability:Public (worldwide access)
Research Advisor(s):
  • Kochmann, Dennis M.
Thesis Committee:
  • Ortiz, Michael (chair)
  • Ravichandran, Guruswami
  • Desbrun, Mathieu
  • Kochmann, Dennis M.
Defense Date:24 June 2015
Non-Caltech Author Email:jeff.amelang (AT) gmail.com
Funders:
Funding AgencyGrant Number
Department of Energy National Nuclear Security AdministrationDE-FC52-08NA28613
National Science FoundationCMMI-1234364
Record Number:CaltechTHESIS:09152015-212147583
Persistent URL:http://resolver.caltech.edu/CaltechTHESIS:09152015-212147583
DOI:10.7907/Z9SB43PH
Related URLs:
URLURL TypeDescription
http://dx.doi.org/10.1016/j.jmps.2015.03.007DOIArticle adapted for ch. 2
Default Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:9155
Collection:CaltechTHESIS
Deposited By: Jeffrey Amelang
Deposited On:25 Sep 2015 16:11
Last Modified:18 May 2017 17:21

Thesis Files

[img]
Preview
PDF - Final Version
See Usage Policy.

287Mb

Repository Staff Only: item control page