[論文レビュー] A Posteriori Error Estimates for Energy-Based Quasicontinuum Approximations of a Periodic Chain

Hao Wang|arXiv (Cornell University)|Dec 22, 2011

Microstructure and mechanical properties参考文献 8被引用数 4

ひとこと要約

本稿では、最近接および第二近隣相互作用を有する一次元周期的鎖に対する一貫エネルギー基準準連続体（CEB-QC）法の事後誤差推定を提示する。残差をエネルギー近似誤差と力の近似誤差に分解し、それらを負のソボレフノルムで評価することで、変形勾配の$L^2$ノルムにおける事後誤差推定器を導出し、全エネルギー差へと拡張した。これにより、数値的検証を伴う適応的メッシュ細分化が可能になった。

ABSTRACT

We present a posteriori error estimates for a recently developed atomistic/continuum coupling method, the Consistent Energy-Based QC Coupling method. The error estimate of the deformation gradient combines a residual estimate and an a posteriori stability analysis. The residual is decomposed into the residual due to the approximation of the stored energy and that due to the approximation of the external force, and are bounded in negative Sobolev norms. In addition, the error estimate of the total energy using the error estimate of the deformation gradient is also presented. Finally, numerical experiments are provided to illustrate our analysis.

研究の動機と目的

• 一貫エネルギー基準準連続体（CEB-QC）結合法の厳密な事後誤差解析を構築し、一貫性およびエネルギー保存を保証すること。
• 特に連続領域における非一様有限要素メッシュを有する問題に対して、エネルギー基準QC法における事後誤差推定器の欠如を是正すること。
• 残差と安定性解析に基づき、要素ごとに分解可能な誤差推定器を導出し、適応的メッシュ細分化を可能とすること。
• 残差の負のソボレフノルムにおける境界と、QC解に依存する事後安定性定数を組み合わせた誤差推定のフレームワークを提供すること。
• 2次元的設定におけるバニティ型欠陥を対象とした数値実験を通じて、理論的誤差推定器の妥当性を検証すること。

提案手法

• 連続的区分線形関数と離散的ベクトルを用いて原子的モデルを定式化し、変位および変形のハイブリッド連続的・離散的表現を可能とした。
• 残差を2つの部分に分解した：保存エネルギーの近似に起因する部分と、外部力の近似に起因する部分。両者とも負のソボレフノルムで評価した。
• 事後安定性解析を実施し、QC解に依存する安定性定数を導出し、誤差推定の堅牢性を確保した。
• 変形勾配の$L^2$ノルム誤差推定器は、残差推定と安定性定数を組み合わせることで構築した。
• 全エネルギー差の誤差推定器は、変形勾配誤差推定器を用いて導出し、物理的関心の高い量（エネルギー）の誤差制御を可能とした。
• 要素ごとのグローバル誤差推定器への寄与に基づき、適応的メッシュ細分化アルゴリズムを構築し、数値的に検証した。

実験結果

リサーチクエスチョン

• RQ1 非一様有限要素メッシュを有する1次元周期的鎖における一貫エネルギー基準準連続体法の事後誤差推定をどのように導出できるか？
• RQ2 保存エネルギー近似誤差と外部力近似誤差が全体の残差に与える寄与は何か？ また、それらはどのように負のソボレフノルムで評価できるか？
• RQ3 安定性定数がQC解に依存する形で構築可能な事後安定性解析は可能か？ これにより信頼性の高い誤差推定が可能になるか？
• RQ4 残差と安定性成分を用いて、変形勾配の$L^2$ノルムにおける誤差をどのように推定できるか？
• RQ5 変形勾配の誤差推定器を全エネルギーの誤差推定に拡張可能か？ また、その有効性は適応的細分化の誘導にどの程度寄与するか？

主な発見

• 保存エネルギー近似に起因する残差と、外部力近似に起因する残差の両方が、負のソボレフノルムで評価可能であり、厳密な残差に基づく誤差推定が可能となった。
• 事後安定性定数はQC解に依存する関数として導出され、誤差推定の鋭さを確保するための非定数かつ解に依存する境界を提供した。
• 変形勾配の$L^2$ノルム誤差推定器は、残差推定と安定性定数を組み合わせることで構築され、グローバルに有効な誤差境界を提供した。
• 全エネルギー差の誤差推定器は、変形勾配誤差推定器から導出され、エネルギー（応用分野における主要な物理的量）の誤差制御を可能とした。
• 数値実験により、誤差推定器が適応的メッシュ細分化を効果的に誘導することが示された。特に2次元的設定における局在的欠陥（バニティ）の捉え込みに優れた性能を示した。
• 導出された誤差推定器は要素ごとに分解可能であり、連続領域における効率的かつ局所的な適応的細分化戦略を可能とした。