[論文レビュー] Iterative Methods for the Force-based Quasicontinuum Approximation

Matthew Dobson, Mitchell Luskin+1 more|arXiv (Cornell University)|Oct 12, 2009

Nonlocal and gradient elasticity in micro/nano structures参考文献 28被引用数 22

ひとこと要約

本稿では、力に基づく準連続体（QCF）方程式の線形化系を解くためのプリコンディショニング付きGMRES法を提案する。QCF近似の非保存的性質により、この方程式系は非対称かつ不定型である。本手法はQCLプリコンディショナーと特化した内積を用い、臨界ひずみに達するまで安定かつ信頼性のある収束を達成し、残差が有効な誤差予測子として機能する。

ABSTRACT

Force-based atomistic-continuum hybrid methods are the only known pointwise consistent methods for coupling a general atomistic model to a finite element continuum model. For this reason, and due to their algorithmic simplicity, force-based coupling methods have become a popular class of atomistic-continuum hybrid models as well as other types of multiphysics models. However, the recently discovered unusual stability properties of the linearized force-based quasicontinuum (QCF) approximation, especially its indefiniteness, present a challenge to the development of efficient and reliable iterative methods. We present analytic and computational results for the generalized minimal residual (GMRES) solution of the linearized QCF equilibrium equations. We show that the GMRES method accurately reproduces the stability of the force-based approximation and conclude that an appropriately preconditioned GMRES method results in a reliable and efficient solution method.

研究の動機と目的

• 線形化力に基づく準連続体（QCF）方程式を解くために、従来の反復法（例えば、修正非線形共役勾配法やゴースト力補正（GFC））の数値的不安定性を解消すること。
• 標準的な反復法の収束性と安定性を損なう、線形化QCF作用素の不定性および非対称性という課題を克服すること。
• 臨界ひずみに達するまで、元の原子論的モデルの安定性特性を正確に捉えることができる信頼性があり効率的な反復解法を開発すること。
• 残差が実用的シミュレーションにおける収束モニタリングに有効である、信頼性の高い誤差指標として機能することを保証すること。

提案手法

• 非対称かつ不定型の線形系に強い一般化最小残差（GMRES）法を、線形化QCF系に適用し、その堅牢性を活用する。
• GMRESの反復収束性と安定性を向上させるために、QCL（局所補正付き準連続体）法をプリコンディショナーとして用いる。
• GMRESにおける基本的な内積として、QCL作用素に基づく特化した$υ^{1,2}$-内積を採用し、条件数の改善と収束特性の向上を図る。
• 高重複固有値を扱うために、ブロック構造とデフラーション技術を用いて一般化固有値問題$L^{-1}L^{\text{qcf}}_F$の固有基底を構築し、スペクトル解析における数値的安定性を確保する。
• 残差と誤差ノルムの両方を用いて、異なる系サイズ$N$、細分化レベル$K$、力パラメータ$f$を想定した1次元モデル問題に対して、本手法を実装・検証する。
• 収束性と信頼性を定量的に評価するために、残差の$\mathcal{U}^{-1,2}$-ノルムと誤差の$\mathcal{U}^{1,2}$-ノルムを用いる。

実験結果

リサーチクエスチョン

• RQ1 修正非線形共役勾配法やゴースト力補正（GFC）といった標準的な反復法は、臨界ひずみに達する前でさえ、線形化QCF系を数値的に安定して解くことができるか？
• RQ2 適切にプリコンディショニングされ、適切な内積を備えたGMRES法は、QCF作用素の不定性および非対称性を克服し、安定な収束を保証できるか？
• RQ3 GMRES反復における残差は、QCF系の真の誤差を信頼性高く予測できるか？
• RQ4 プリコンディショニング付きGMRES法の収束速度は、システムの固有値分布に基づく理論的予測と一致するか？
• RQ5 臨界ひずみに達するまで、QCF系が不安定化するが、本手法は原子論的モデルの正しい安定性特性を保持しているか？

主な発見

• 修正非線形共役勾配法は、線形化QCF作用素の不定性に起因する不安定性により、臨界ひずみに達する前でも数値的に失敗する。
• QCEエネルギーをプリコンディショナーとして用いるゴースト力補正（GFC）法は、臨界ひずみに達する前でも不安定化し、欠陥形成の臨界ひずみが不正確に低く予測される。
• QCLプリコンディショナーと$\mathcal{U}^{1,2}$-内積を備えた本手法のプリコンディショニング付きGMRES法は、臨界ひずみに達するまで安定な収束を達成し、システムの正しい安定性挙動を保持する。
• 残差ノルム$\|r^{(m)}\|_{\mathcal{U}^{-1,2}}$は、理論的予測（命題5.3）と一致する線形収束レート$q = \frac{1 - \sqrt{A_F / \phi_F''}}{1 + \sqrt{A_F / \phi_F''}}$で減少する。
• 誤差ノルム$\|e^{(m)}\|_{\mathcal{U}^{1,2}}$は、残差ノルム$\|L^{-1/2}r^{(m)}\|_{\ell^2_\varepsilon}$と密接に一致しており、残差が信頼性のある収束指標であることが確認される。
• 本手法は、異なる系サイズ$N$、細分化レベル$K$、右辺ベクトル$f$に対して、1次元モデル問題においてもスケーラビリティと信頼性を維持し、堅牢である。