[论文解读] Convergence of a force-based hybrid method for atomistic and continuum models in three dimension

Jianfeng Lu, Pingbing Ming|arXiv (Cornell University)|Feb 12, 2011

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

一句话总结

该论文首次在三维空间中建立了基于力的原子-连续耦合方法的收敛性结果，证明了当晶格间距趋于零时，混合解相对于原子解具有二次收敛性。该分析通过伪差分算子结合一致性和线性稳定性，确保在一般短程势和简单晶格下实现最优收敛。

ABSTRACT

We study a force-based hybrid method that couples atomistic models with nonlinear Cauchy-Born elasticity models. We show that the proposed scheme converges quadratically to the solution of the atomistic model, as the ratio between lattice parameter and the characteristic length scale of the deformation tends to zero. Convergence is established for general short-ranged atomistic potential and for simple lattices in three dimension. The convergence is based on consistency and stability analysis. General tools are developed in the framework of pseudo-difference operators for stability analysis in arbitrary dimension of the multiscale atomistic and continuum coupling methods.

研究动机与目标

• 建立三维空间中基于力的混合方法，耦合原子模型与Cauchy-Born连续模型的收敛性。
• 解决多维多尺度耦合方法中长期存在的稳定性分析难题，超越一维情形。
• 基于伪差分算子发展通用工具，用于分析任意维度下的稳定性。
• 在对势和晶格结构的最小假设下，证明混合解相对于原子解的二次收敛性。

提出的方法

• 提出一种基于力的混合方法，在原子区域与Cauchy-Born连续区域的界面处耦合原子力与连续力。
• 采用基于力平衡方程的变分形式，而非能量最小化，从而避免一致耦合中的鬼力现象。
• 利用伪差分算子框架分析线性化算子，以处理三维空间中的多尺度结构。
• 通过伪差分算子的正则性估计、线性化格式的一致性以及连续问题的稳定性，综合建立稳定性。
• 采用Strang型收敛性分析，界定混合解与真实原子解之间的误差。
• 应用Brouwer不动点定理，证明解在原子解附近小邻域内的存在性与局部唯一性。

实验结果

研究问题

• RQ1 基于力的原子-连续耦合方法是否能在三维空间中实现二次收敛？
• RQ2 分析高于一维的多尺度耦合方案稳定性，需要哪些通用工具？
• RQ3 如何严格证明基于力的混合方法在一般短程势下的一致性和稳定性？
• RQ4 能否通过精心设计的耦合方式与分析，克服三维系统中的鬼力问题？
• RQ5 当晶格间距趋于零时，混合解相对于原子解的收敛速率是多少？

主要发现

• 当晶格间距 ε → 0 时，基于力的混合方法相对于原子解实现二次收敛，误差被控制在 Cε² 以内。
• 该收敛性结果适用于三维空间中的一般短程两体势和简单布拉维晶格。
• 通过伪差分算子进行稳定性分析，使结果可推广至任意维度。
• 通过证明力平衡方程残差的二次衰减，验证了线性化算子的一致性。
• 由于解点处Hessian矩阵的非退化性，混合方法的解具有局部唯一性。
• 该分析提供了一个通用的多尺度方法稳定性与收敛性框架，可推广至特定混合模型之外。