Massively parallel finite difference elasticity using block-structured adaptive mesh refinement with a geometric multigrid solver
arXiv:2001.04789 · doi:10.1016/j.jcp.2020.110065
Abstract
Computationally solving the equations of elasticity is a key component in many materials science and mechanics simulations. Phenomena such as deformation-induced microstructure evolution, microfracture, and microvoid nucleation are examples of applications for which accurate stress and strain fields are required. A characteristic feature of these simulations is that the problem domain is simple (typically a rectilinear representative volume element (RVE)), but the evolution of internal topological features is extremely complex. Traditionally, the finite element method (FEM) is used for elasticity calculations; FEM is nearly ubiquituous due to (1) its ability to handle meshes of complex geometry using isoparametric elements, and (2) the weak formulation which eschews the need for computation of second derivatives. However, variable topology problems (e.g. microstructure evolution) require either remeshing, or adaptive mesh refinement (AMR) - both of which can cause extensive overhead and limited scaling. Block-structured AMR (BSAMR) is a method for adaptive mesh refinement that exhibits good scaling and is well-suited for many problems in materials science. Here, it is shown that the equations of elasticity can be efficiently solved using BSAMR using the finite difference method. The boundary operator method is used to treat different types of boundary conditions, and the "reflux-free" method is introduced to efficiently and easily treat the coarse-fine boundaries that arise in BSAMR. Examples are presented that demonstrate the use of this method in a variety of cases relevant to materials science: Eshelby inclusions, fracture, and microstructure evolution. Reasonable scaling is demonstrated up to 4000 processors with tens of millions of grid points, and good AMR efficiency is observed.
References in corpus (1)
Cited by in corpus (10)
- Block structured adaptive mesh refinement and strong form elasticity approach to phase field fracture with applications to delamination, crack branching and crack deflection
- Multiphase field modeling of grain boundary migration mediated by emergent disconnections
- TAC: Optimizing Error-Bounded Lossy Compression for Three-Dimensional Adaptive Mesh Refinement Simulations
- TAC+: Optimizing Error-Bounded Lossy Compression for 3D AMR Simulations
- Atomistic-informed phase field modeling of magnesium twin growth by disconnections
- Robust, strong form mechanics on an adaptive structured grid: efficiently solving variable-geometry near-singular problems with diffuse interfaces
- The Alamo multiphysics solver for phase field simulations with strong-form mechanics and block structured adaptive mesh refinement
- Comparison of evolving interfaces, triple points, and quadruple points for discrete and diffuse interface methods
- A ghost-point smoothing strategy for geometric multigrid on curved boundaries
- Mechanics of incompatible asymmetric grain boundary migration