[论文解读] A Posteriori Error Estimation and Adaptive Algorithm for the Atomistic/Continuum Coupling in 2D

Hao Wang, Mingjie Liao+2 more|arXiv (Cornell University)|Feb 9, 2017

Numerical methods in engineering参考文献 34被引用 4

一句话总结

本文针对二维一致的原子/连续（a/c）耦合方法，提出了一套严格的后验误差估计与自适应算法，采用负Sobolev范数中的残差型误差界。该方法在自由度方面实现了最优收敛率，与先验估计一致，并通过结合残差与能量型误差指标，实现了高效的目标导向性网格加密。

ABSTRACT

Atomistic/continuum coupling methods aim to achieve optimal balance between accuracy and efficiency. Adaptivity is the key for the efficient implementation of such methods. In this paper, we carry out a rigorous a posteriori analysis of the residual, the stability constant, and the error bound, for a consistent atomistic/continuum coupling method in 2D. We design and implement the corresponding adaptive mesh refinement algorithm, and the convergence rate with respect to degrees of freedom is optimal compare with a priori error estimates.

研究动机与目标

• 为基于能量的二维一致原子/连续耦合方法建立数学上严谨的后验误差估计框架。
• 设计并实现一种自适应网格加密算法，通过针对性地优化原子区与连续区的误差源，提升计算效率。
• 建立考虑建模误差以及原子-连续界面附近网格非均匀性的后验误差界。
• 通过数值实验验证自适应算法，并展示其收敛率与先验最优估计相当。

提出的方法

• 在负Sobolev范数中采用残差型方法推导后验误差估计，将一维方法的技术推广至二维。
• 引入一个无散度张量场，以表征一维与二维a/c耦合分析之间的结构差异。
• 基于残差泛函与稳定性常数构造误差估计量，纳入原子到连续过渡区域的建模误差。
• 设计一种自适应算法，根据组合的残差与能量误差指标，动态调整原子区域并加密连续区域网格。
• 采用基于区域化与次样条插值的应力张量公式，实现对一般短程相互作用的一致误差估计。
• 实现两种加密策略：残差驱动与能量驱动，通过可调参数（如τ₃）控制区域扩展与网格加密程度。

实验结果

研究问题

• RQ1 能否在超越一维情形的二维一致原子/连续耦合方法中，建立严格的后验误差估计？
• RQ2 如何将界面附近的建模误差与网格非均匀性纳入后验误差界？
• RQ3 何种自适应加密策略可使a/c耦合问题在自由度方面实现最优收敛率？
• RQ4 残差型与能量型误差指标在指导a/c耦合网格自适应方面有何异同？
• RQ5 所提出的框架能否推广至一般短程相互作用及复杂缺陷（如位错或空位）？

主要发现

• 后验误差估计量在负Sobolev范数中对能量范数误差提供了可靠的上界，具有严格的数学依据。
• 自适应算法在H¹误差与能量差方面均实现了最优收敛率，与理论先验估计一致。
• 数值实验中τ₃ = 0.3与τ₃ = 0.7的对比表明，残差驱动与能量驱动策略的收敛行为相近。
• 该方法能有效处理原子-连续界面附近的建模误差与非均匀网格加密，避免出现虚假力。
• 通过定义良好的延拓算子，成功将方法扩展至空位缺陷，保持能量等价性与梯度界。
• 该框架可推广至其他一致a/c耦合格式（如BQCE、BQCF与BGFC），通过应力张量公式实现。