Abstract:This paper introduces the MeshAC package, which generates three-dimensional adaptive meshes tailored for the efficient and robust implementation of multiscale coupling methods. While Delaunay triangulation is commonly used for mesh generation across the entire computational domain, generating meshes for multiscale coupling methods is more challenging due to intrinsic discrete structures such as defects, and the need to match these structures to the continuum domain at the interface. The MeshAC package tackles these challenges by generating meshes that align with fine-level discrete structures. It also incorporates localized modification and reconstruction operations specifically designed for interfaces. These enhancements improve both the implementation efficiency and the quality of the coupled mesh. Furthermore, MeshAC introduces a novel adaptive feature that utilizes gradient-based a posteriori error estimation, which automatically adjusts the atomistic region and continuum mesh, ensuring an optimal balance between accuracy and efficiency. This package can be directly applied to the geometry optimization problems of a/c coupling in static mechanics, with potential extensions to many other scenarios. Its capabilities are demonstrated for complex material defects, including straight edge dislocation in BCC W and double voids in FCC Cu. These results suggest that MeshAC can be a valuable tool for researchers and practitioners in computational mechanics.

What problem does this paper attempt to address?

The problem this paper attempts to address is: The challenges of generating and adapting 3D meshes in atomistic-to-continuum (a/c coupling) simulations. Specifically, the paper aims to develop an open-source software package named **MeshAC**, which is dedicated to addressing mesh generation and adaptation issues in a/c coupling methods for physically relevant 3D scenarios. These challenges include: 1. **Handling lattice structures and complex defects**: The generated mesh needs to account for crystal structures (such as body-centered cubic BCC and face-centered cubic FCC) as well as complex crystal defects (such as dislocations, voids, grain boundaries, and cracks), leading to anisotropic atomic regions and computational domains. 2. **Interface matching**: To ensure effective a/c coupling, the mesh must align with the underlying atomic grid, especially at the a/c interface. As the interface evolves, complex mesh operations are required at these interfaces to balance accuracy and computational cost. 3. **Adaptive optimization**: To achieve the optimal balance between cost and accuracy, robust a posteriori error estimators need to be utilized to drive the redistribution of the a/c interface, computational domain boundaries, and the mesh. By addressing these issues, MeshAC can improve the efficiency and quality of implementing multiscale coupling methods, particularly suitable for geometric optimization problems in static mechanics, and can be extended to other scenarios such as discrete-continuum coupling, random alloys, truss structures, polymer networks, and magnetization dynamics. The paper demonstrates the capabilities of MeshAC by applying the BGFC method (Blended Ghost Force Correction atomistic/continuum hybrid method) to simulate real material defects (such as straight dislocations and double voids), showing that MeshAC is a valuable tool for studying complex material defect behaviors in computational mechanics.

MeshAC: A 3D Mesh Generation and Adaptation Package for Multiscale Coupling Methods

Geometry-based adaptive mesh generation for continuous and discrete parametric surfaces

Geometry-Based Adaptive Parametric Surface Mesh Generation

Unstructured tetrahedral mesh generation for complex configurations

An Extended Advancing Front Technique For Closed Surfaces Mesh Generation

Advancing Front Mesh Generation on Dirty Composite Surfaces

Geometry-based Adaptive Surface Mesh Generation for STL Models

Mesh Adaptation for Curing the Pathological Behaviors of an Upwind Scheme

An Optimization Model and Algorithm for the Smoothing of Sizing Functions of Three-Dimensional Finite Element Meshes

Adaptive Multiscale Coupling Methods of Molecular Mechanics based on a Unified Framework of a Posteriori Error Estimates

Automatic Mesh Generation Techniques And Adaptive Softwares

Adaptive Multigrid Strategy for Geometry Optimization of Large-Scale Three Dimensional Molecular Mechanics

HybridOctree_Hex: Hybrid Octree-Based Adaptive All-Hexahedral Mesh Generation with Jacobian Control

Hi-performance Decoupled Parallel Mesh Generation Algorithm for Unstructured Tetrahedral Mesh

A Posteriori Error Estimate and Adaptive Mesh Refinement Algorithm for Atomistic/Continuum Coupling with Finite Range Interactions in Two Dimensions

Generation of adaptive refinement tetrahedral meshes over domains of multi-chamber partition

On Multi-Mesh H -Adaptive Methods

Adaptive and Parallel Multiscale Framework for Modeling Cohesive Failure in Engineering Scale Systems

An Adaptive and Highly Accurate FDTD Mesh Generation Technique for Objects with Complex Edge Structures

A novel unresolved/semi‐resolved CFD‐DEM coupling method with dynamic unstructured mesh

An automatic mesh generator for coupled 1D–2D hydrodynamic models