Large-scale simulations of microstructural evolution in elastically stressed solids

Norio Akaiwa, K. Thornton*, P. W. Voorhees

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

33 Scopus citations

Abstract

We present a set of numerical methods for simulations of microstructural evolution in elastically stressed solids. We combine three powerful tools to achieve computational efficiency: the boundary integral method to provide excellent resolution, the fast multipole method to reduce computational cost, and the small-scale decomposition technique with a two-level time-stepping scheme to remove the stiffness from the time advance. Although we apply these methods to study the topic of our interest, the details of how they are implemented can be useful in many other situations. We extend the fast multipole method to calculate the anisotropic stress field in periodic two-dimensional domains and in doing so address issues associated with the conditional convergence of the summations. In addition, we introduce a new formula for the potential in periodically arranged two-dimensional cells in the absence of an applied field through a summation in physical space without using the Ewald sum. Furthermore, we implement a time-stepping scheme that enables us to speed up the calculation by an additional factor of 100 over a straightforward implementation of the small-scale decomposition technique. The computational complexity scales as the number of mesh points N, and thus we are able to employ N ∼ 500, 000 in a typical calculation.

Original languageEnglish (US)
Pages (from-to)61-86
Number of pages26
JournalJournal of Computational Physics
Volume173
Issue number1
DOIs
StatePublished - Oct 10 2001

Funding

KT is grateful for the hospitality of, and stimulating interaction with, W. C. Carter and his group at the Massachusetts Institute of Technology. We also thank National Science Foundation for its support through Grant DMR-9707073.

Keywords

  • Anisotropic elasticity
  • Boundary integral method
  • Fast multiple method
  • Ostwald ripening
  • Periodic potential

ASJC Scopus subject areas

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • General Physics and Astronomy
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'Large-scale simulations of microstructural evolution in elastically stressed solids'. Together they form a unique fingerprint.

Cite this