[论文解读] Homogenization-based Analysis of Quasicontinuum Method for Complex Crystals

Assyr Abdulle, Ping Lin+1 more|arXiv (Cornell University)|Jun 2, 2010

Advanced Mathematical Modeling in Engineering参考文献 31被引用 7

一句话总结

本文提出一种基于均匀化方法的框架，用于分析复杂晶体的非连续元（QC）方法，在一维中建立了收敛性，并证明了离散均匀化与改进的柯西方-玻恩QC公式的等价性。该方法可实现严格的误差估计，并自然地推广至二维情形，数值结果证实了在 $L^2$ 和 $H^1$ 范数下均达到最优收敛率。

ABSTRACT

Among the efficient numerical methods based on atomistic models, the quasicontinuum (QC) method, introduced by Tadmor, Ortiz, and Phillips (1996), has attracted growing interest in recent years. Originally, the QC method was developed for materials with simple crystalline lattice (simple crystals) and later was extended to complex lattice (Tadmor et al, 1999). In the present paper we formulate the QC method for complex lattices in a homogenization framework and perform analysis of such a method in a 1D setting. We also present numerical examples showing that the convergence results are valid in a more general setting.

研究动机与目标

• 开发一种严格的数学框架，用于分析复杂晶态材料中的非连续元（QC）方法，以弥补传统分析的不足。
• 利用离散数值均匀化技术，建立复杂晶格中QC方法的收敛性。
• 通过在均匀化框架内构建QC方法，弥合原子模型与连续力学之间的鸿沟。
• 提供均匀化QC解相对于完整原子解的误差估计。
• 通过稳健的均匀化基础，将QC方法的适用性扩展至有限温度和非晶态材料。

提出的方法

• 在离散均匀化框架内，为复杂晶格中的QC方法建模，将原子系统视为多尺度问题。
• 采用改进的柯西方-玻恩法则，描述复杂晶体中亚晶格之间的相对位移。
• 应用数值均匀化技术，实时计算有效势能，避免预先计算。
• 构建宏观-微观数值格式，利用粗尺度自由度捕捉复杂材料的有效行为。
• 实施后处理重构，从均匀化解中近似得到完整的原子解。
• 通过有限元离散化和自适应采样区域，在一维和二维设置中验证该方法。

实验结果

研究问题

• RQ1 能否使用均匀化框架，对具有多个亚晶格的复杂晶体QC方法进行严格分析？
• RQ2 对于线性、周期性的原子系统，均匀化QC方法在一维中的收敛速率是多少？
• RQ3 离散均匀化方法是否与基于改进柯西方-玻恩法则的标准QC公式等价？
• RQ4 在一维中推导出的误差估计如何推广至非线性和多维问题？
• RQ5 基于均匀化方法的QC方法能否有效应用于具有不同键刚度的二维复杂晶体体系？

主要发现

• 在1D中，均匀化QC方法在 $H^1$-范数下具有第一阶精度，在 $L^2$-范数下具有第二阶精度，与理论预测一致。
• 二维中的数值实验证实，$L^2$ 和 $H^1$ 误差均保持相同的收敛速率，验证了一维分析在更高维中的适用性。
• 在细密网格尺寸下，$L^2$-误差出现平台区，表明具有最优收敛行为。
• 离散均匀化框架与基于改进柯西方-玻恩法则的QC方法等价，确立了理论一致性。
• 该方法可通过后处理实现对完整原子解的精确重构，且误差界由均匀化过程导出。
• 该方法可扩展至有限温度模拟和非晶态材料，为未来多尺度建模提供了可行路径。