[论文解读] A Posteriori Error Estimates for Energy-Based Quasicontinuum Approximations of a Periodic Chain

Hao Wang|arXiv (Cornell University)|Dec 22, 2011

Microstructure and mechanical properties参考文献 8被引用 4

一句话总结

本文针对具有最近邻和次近邻相互作用的一维周期性链中的相容能量型非连续元（CEB-QC）方法，提出了后验误差估计。通过将残差分解为能量逼近误差和力逼近误差，并在负Sobolev范数下进行有界性分析，作者推导出变形梯度在$L^2$-范数下的后验误差估计器，并进一步扩展至总能量差，从而实现自适应网格加密，且通过数值验证了其有效性。

ABSTRACT

We present a posteriori error estimates for a recently developed atomistic/continuum coupling method, the Consistent Energy-Based QC Coupling method. The error estimate of the deformation gradient combines a residual estimate and an a posteriori stability analysis. The residual is decomposed into the residual due to the approximation of the stored energy and that due to the approximation of the external force, and are bounded in negative Sobolev norms. In addition, the error estimate of the total energy using the error estimate of the deformation gradient is also presented. Finally, numerical experiments are provided to illustrate our analysis.

研究动机与目标

• 为相容能量型非连续元（CEB-QC）耦合方法建立严格的后验误差分析，以确保一致性和能量守恒。
• 解决基于能量的QC方法中缺乏后验误差估计器的问题，特别是针对连续区域中非均匀有限元网格的问题。
• 通过基于残差和稳定性分析的逐单元可分解误差估计器，实现自适应网格加密。
• 提供一种误差估计框架，结合负Sobolev范数下的残差界与依赖于QC解的后验稳定性常数。
• 通过在二维类似设置下的空位缺陷数值实验，验证理论误差估计器的有效性。

提出的方法

• 该方法采用连续分段线性函数和离散向量表示原子模型，从而实现位移和形变的混合连续-离散表示。
• 将残差分解为两部分：一部分来自存储能量的逼近，另一部分来自外力的逼近，两者均在负Sobolev范数下有界。
• 进行后验稳定性分析，以推导依赖于QC解的稳定性常数，从而确保误差估计的鲁棒性。
• 通过结合残差估计和稳定性常数，构建变形梯度在$L^2$-范数下的误差估计器。
• 基于变形梯度误差估计器，推导出总能量差的误差估计器，从而实现对物理量关键指标的误差控制。
• 开发了基于全局误差估计器中各单元贡献的自适应网格加密算法，并通过数值方法进行了验证。

实验结果

研究问题

• RQ1 如何为具有非均匀有限元网格的一维周期性链中的相容能量型非连续元方法推导后验误差估计？
• RQ2 存储能量逼近误差和外力逼近误差对整体残差的贡献是什么？如何在负Sobolev范数下对它们进行有界性分析？
• RQ3 能否构建一种后验稳定性分析，使得稳定性常数依赖于QC解，从而实现可靠的误差估计？
• RQ4 如何利用残差和稳定性分量，在$L^2$-范数下估计变形梯度的误差？
• RQ5 能否将变形梯度的误差估计器扩展为估计总能量的误差？其在引导自适应加密中的有效性如何？

主要发现

• 存储能量逼近引起的残差和外力逼近引起的残差均在负Sobolev范数下有界，从而实现了严格的基于残差的误差估计。
• 后验稳定性常数被推导为QC解的函数，提供了非恒定的、依赖于解的有界性，这对于精确的误差估计至关重要。
• 通过结合残差估计和稳定性常数，构建了变形梯度在$L^2$-范数下的误差估计器，从而得到全局有效的误差界。
• 总能量差的误差估计器由变形梯度误差估计器推导而来，从而实现了对能量这一关键物理量的误差控制。
• 数值实验表明，误差估计器在引导自适应网格加密方面具有有效性，尤其在捕捉二维类似设置下的局域缺陷（如空位）方面表现突出。
• 所推导的误差估计器具有逐单元可分解性，从而在连续区域中实现了高效且局部化的自适应加密策略。