[论文解读] An Optimization-Based Atomistic-to-Continuum Coupling Method

Derek Olson, Pavel Bochev+2 more|arXiv (Cornell University)|Apr 17, 2013

Block Copolymer Self-Assembly参考文献 15被引用 4

一句话总结

该论文提出了一种基于优化的原子-连续体（AtC）耦合方法，将耦合问题表述为一个约束优化问题，以最小化重叠区域中原子与连续体位移之间的不匹配。通过在界面上对原子和连续体模型分别施加虚拟狄利克雷控制，该方法避免了混合处理，并在一维线性化最近邻模型中实现了最优收敛速率。

ABSTRACT

We present a new optimization-based method for atomistic-to-continuum (AtC) coupling. The main idea is to cast the coupling of the atomistic and continuum models as a constrained optimization problem with virtual Dirichlet controls on the interfaces between the atomistic and continuum subdomains. The optimization objective is to minimize the error between the atomistic and continuum solutions on the overlap between the two subdomains, while the atomistic and continuum force balance equations provide the constraints. Splitting of the atomistic and continuum problems instead of blending them and their subsequent use as constraints in the optimization problem distinguishes our approach from the existing AtC formulations. We present and analyze the method in the context of a one-dimensional chain of atoms modeled using a linearized two-body next-nearest neighbor interactions.

研究动机与目标

• 开发一种新型AtC耦合方法，通过独立处理原子和连续体模型，避免能量或力的混合。
• 通过约束优化最小化重叠区域中原子与连续体解之间的差异。
• 利用人工界面上的虚拟狄利克雷控制，建立理论收敛性与稳定性。
• 提供一个框架，使得方法的关键特性可直接从原子和连续体组分的独立性质中推导得出。
• 在一维线性化两体相互作用模型中展示最优收敛速率。

提出的方法

• 该方法将AtC耦合表述为一个带有虚拟狄利克雷控制的约束优化问题，控制对象为原子与连续体子域之间的界面。
• 目标函数最小化重叠区域Ωₒ内原子与连续体位移差异的ℓ²-范数。
• 约束条件分别在原子子域Ωₐ和连续体子域Ωᶜ中强制实现力平衡。
• 该方法将原子与连续体问题分离处理，而非进行混合，从而与现有基于能量或力的混合方法区分开来。
• 该方法基于位移场的变分形式，并通过虚拟控制施加边界条件，以确保问题的适定性。
• 理论分析依赖于求解位移模态的线性系统，并利用矩阵逆范数估计与几何级数进行范数估计。

实验结果

研究问题

• RQ1 能否通过优化方法构建一种非混合的AtC方法，以最小化重叠区域内的位移不匹配？
• RQ2 在人工界面上使用虚拟狄利克雷控制如何影响耦合系统的稳定性与收敛性？
• RQ3 与现有基于混合的公式相比，该基于优化的AtC方法的收敛速率如何？
• RQ4 耦合方法的关键特性能否直接从原子和连续体模型的独立性质中推导得出？
• RQ5 在渐近缩放下，位移解中的系数如何变化，这对误差界意味着什么？

主要发现

• 该方法在位移误差的ℓ²-范数下实现了最优收敛速率，误差被一个常数乘以网格尺寸所界定。
• 原子区域中位移解的范数被矩阵T_L的逆所界定，当L→∞时，T_L收敛于一个固定矩阵T。
• 重叠区域的误差估计依赖于二次型的特征值，最小特征值下界为(1/24)γ²（当γ较小时）。
• 当适当缩放时，收敛速率与重叠尺寸无关，表明对域分解具有鲁棒性。
• 理论分析证实，该方法的稳定性和收敛性可仅从原子和连续体模型的性质中推导得出。
• 该方法通过最小化不匹配而非混合能量或力，避免了传统混合AtC方法中常见的相容性与稳定性问题。