[論文レビュー] A Posteriori Error Estimation and Adaptive Algorithm for the Atomistic/Continuum Coupling in 2D

Hao Wang, Mingjie Liao+2 more|arXiv (Cornell University)|Feb 9, 2017

Numerical methods in engineering参考文献 34被引用数 4

ひとこと要約

本稿では、2次元における一貫性のある原子論的/連続体（a/c）結合法について、負のソボレフノルムにおける残差に基づく誤差境界を用いて、きめ細やかな事後誤差推定と適応型アルゴリズムを提示する。この手法は自由度に関して最適収束率を達成し、事前推定と一致する。また、残差とエネルギーに基づく誤差指標を組み合わせることで、効率的かつ目的指向のメッシュ細分化を可能にする。

ABSTRACT

Atomistic/continuum coupling methods aim to achieve optimal balance between accuracy and efficiency. Adaptivity is the key for the efficient implementation of such methods. In this paper, we carry out a rigorous a posteriori analysis of the residual, the stability constant, and the error bound, for a consistent atomistic/continuum coupling method in 2D. We design and implement the corresponding adaptive mesh refinement algorithm, and the convergence rate with respect to degrees of freedom is optimal compare with a priori error estimates.

研究の動機と目的

• 2次元におけるエネルギーに基づく一貫性のある原子論的/連続体結合法の数学的に厳密な事後誤差推定フレームワークの構築。
• 原子論的領域と連続体領域の両方の誤差源に焦点を当て、計算効率を最適化するための適応的メッシュ細分化アルゴリズムの設計および実装。
• 原子論的-連続体界面近傍における非一様メッシュ細分化とモデリング誤差を考慮した事後誤差境界の確立。
• 数値的検証を通じて、事前推定と同等の収束率を示す適応型アルゴリズムの有効性を検証。

提案手法

• 1次元から2次元への拡張により、負のソボレフノルムにおける残差に基づく事後誤差推定を導出。
• 1次元と2次元のa/c結合解析における構造的差を特徴付けるため、発散なしのテンソル場を導入。
• 残差関数と安定性定数に基づく誤差推定器を構築し、原子論的から連続体への遷移に起因するモデリング誤差を組み込む。
• 残差とエネルギーに基づく誤差指標を統合した適応型アルゴリズムを設計し、動的に原子論的領域を調整し、連続体メッシュを細分化。
• 局所化と準補間に基づく応力テンソル定式化を採用し、一般の短距離相互作用に対して一貫性のある誤差推定を可能にする。
• 2種類の細分化戦略（残差駆動型とエネルギー駆動型）を実装し、領域拡張とメッシュ細分化を制御する調整可能なパラメータ（例：τ₃）を導入。

実験結果

リサーチクエスチョン

• RQ1 1次元の範囲を越えて、一貫性のある2次元原子論的/連続体結合法に対して数学的に厳密な事後誤差推定を確立できるか？
• RQ2 界面近傍におけるモデリング誤差と非一様メッシュ細分化を、事後誤差境界に組み込む方法は何か？
• RQ3 a/c結合問題において、自由度に関して最適収束率を達成するための適応的細分化戦略は何か？
• RQ4 残差に基づく誤差指標とエネルギーに基づく誤差指標は、a/c結合におけるメッシュ適応の誘導に比べてどのように異なるか？
• RQ5 提案されたフレームワークは、一般の短距離相互作用およびディスlocationやバケーションのような複雑な欠陥に対しても拡張可能か？

主な発見

• 事後誤差推定器はエネルギーノルム誤差の信頼性の高い境界を提供し、負のソボレフノルムにおける数学的根拠を有する。
• 適応型アルゴリズムは、H¹誤差とエネルギー差の両方において、理論的事前推定と一致する最適収束率を達成する。
• τ₃ = 0.3 および τ₃ = 0.7 の数値実験において、残差駆動型とエネルギー駆動型の両戦略が同等の収束行動を示す。
• 本手法は、原子論的-連続体界面近傍におけるモデリング誤差と非一様メッシュ細分化を効果的に処理し、誤った力（ゴースト力）の発生を回避する。
• バケーション欠陥への拡張は、エネルギー同等性と勾配バインディングを保持する明確な拡張作用素を介して達成される。
• 応力テンソル定式化を介して、BQCE、BQCF、BGFCなどの他の一貫性のあるa/c結合スキームに対しても、本フレームワークは一般化可能である。