A posteriori analysis and adaptive algorithms for blended type atomistic-to-continuum coupling with higher-order finite elements

• Published in: Computer physics communications Vol. 310; p. 109533
• Main Author: Wang, Yangshuai
• Format: Journal Article
• Language: English
• Published: Elsevier B.V 01.05.2025
• Subjects: A posteriori error estimate; Adaptivity; Atomistic-to-continuum coupling; Crystalline defects; Higher-order finite element methods; Adaptivity; Atomistic-to-continuum coupling; Higher-order finite element methods; A posteriori error estimate; Crystalline defects
• ISSN: 0010-4655
• DOI: 10.1016/j.cpc.2025.109533

The accurate and efficient simulation of material systems with defects using atomistic-to-continuum (a/c) coupling methods is a significant focus in computational materials science. Achieving a balance between accuracy and computational cost requires the application of a posteriori error analysis and adaptive algorithms. In this paper, we provide a rigorous a posteriori error analysis for three common blended a/c methods: the blended energy-based quasi-continuum (BQCE) method, the blended force-based quasi-continuum (BQCF) method, and the atomistic/continuum blending with ghost force correction (BGFC) method. We discretize the Cauchy-Born model in the continuum region using first- and second-order finite element methods, with the potential for extending to higher-order schemes. The resulting error estimator provides both an upper bound on the true error and a reliable lower bound, subject to a controllable truncation term. Furthermore, we offer an a posteriori analysis of the energy error. We develop and implement an adaptive mesh refinement algorithm applied to two typical defect scenarios: a micro-crack and a Frenkel defect. In both cases, our numerical experiments demonstrate optimal convergence rates with respect to degrees of freedom, in agreement with a priori error estimates. •An innovative a posteriori error estimator for blended type a/c coupling is developed.•The estimator is rigorously shown to provide both upper and lower bounds for the actual approximation error.•The estimator is easily generalizable to higher-order finite element methods and energy error.•Adaptive algorithms guiding mesh refinement and region allocation are developed.