We study the effect of "ghost forces" for a quasicontinuum method in three dimension with a planar interface. "Ghost forces" are the inconsistency of the quasicontinuum method across the interface between the atom...We study the effect of "ghost forces" for a quasicontinuum method in three dimension with a planar interface. "Ghost forces" are the inconsistency of the quasicontinuum method across the interface between the atomistic region and the continuum region. Numerical results suggest that "ghost forces" may lead to a negilible error on the solution, while lead to a finite size error on the gradient of the solution. The error has a layer-like profile, and the interfacial layer width is of O(ε). The error in certain component of the displacement gradient decays algebraically from O(1) to O(ε) away from the interface. A surrogate model is proposed and analyzed, which suggests the same scenario for the effect of "ghost forces". Our analysis is based on the explicit solution of the surrogate model.展开更多
基金supported by National Natural Science Foundation of China(Grant Nos.1093201191230203 and 11021101)
文摘We study the effect of "ghost forces" for a quasicontinuum method in three dimension with a planar interface. "Ghost forces" are the inconsistency of the quasicontinuum method across the interface between the atomistic region and the continuum region. Numerical results suggest that "ghost forces" may lead to a negilible error on the solution, while lead to a finite size error on the gradient of the solution. The error has a layer-like profile, and the interfacial layer width is of O(ε). The error in certain component of the displacement gradient decays algebraically from O(1) to O(ε) away from the interface. A surrogate model is proposed and analyzed, which suggests the same scenario for the effect of "ghost forces". Our analysis is based on the explicit solution of the surrogate model.