[論文レビュー] Sharp Stability Estimates for the Force-based Quasicontinuum Method

Matthew Dobson, Mitchell Luskin+1 more|arXiv (Cornell University)|Jul 22, 2009

Microstructure and mechanical properties参考文献 24被引用数 12

ひとこと要約

本稿は、次近接相互作用を伴う1次元原子連続モデルに対する力に基づく準連続（QCF）法の鋭い安定性解析を提供する。QCF作用素は一般に非正規かつ不定的であるが、著者らはその固有値スペクトルがエネルギーに基づく準非局所（QNL）作用素と一致することを示し、臨界破壊荷重まで動的安定性が保証されることを示し、特定の関数空間の選択のもとで逆作用素の一様有界性を確立した。これは欠陥予測における主要な安定性懸念を解消するものである。

ABSTRACT

A sharp stability analysis of atomistic-to-continuum coupling methods is essential for evaluating their capabilities for predicting the formation and motion of lattice defects. We formulate a simple one-dimensional model problem and give a detailed analysis of the stability of the force-based quasicontinuum (QCF) method. The focus of the analysis is the question whether the QCF method is able to predict a critical load at which fracture occurs. Numerical experiments show that the spectrum of a linearized QCF operator is identical to the spectrum of a linearized energy-based quasi-nonlocal quasicontinuum operator (QNL), which we know from our previous analyses to be positive below the critical load. However, the QCF operator is non-normal and it turns out that it is not generally positive definite, even when all of its eigenvalues are positive. Using a combination of rigorous analysis and numerical experiments, we investigate in detail for which choices of "function spaces" the QCF operator is stable, uniformly in the size of the atomistic system. Force-based multi-physics coupling methods are popular techniques to circumvent the difficulties faced in formulating consistent energy-based coupling pproaches. Even though the QCF method is possibly the simplest coupling method of this kind, we anticipate that many of our observations apply more generally.

研究の動機と目的

• 結晶性材料における破壊の臨界荷重を力に基づく準連続（QCF）法が信頼性を持って予測できるかどうかという根本的問題を解決すること。
• 保存則に従わない非保存的・非正規的性質を示すQCF法の安定性を、固有値の正定性を越えて分析すること。
• 系のサイズNが増加する際、線形化QCF作用素の逆作用素が一様に有界となるような離散関数空間の選択を特定すること。
• QCF法のスペクトル的性質と物理的動的安定性との間の関係を、欠陥形成および運動の文脈で確立すること。
• QCF法の格子欠陥動態の予測能力を明確にする、厳密で鋭い安定性推定を提供すること。

提案手法

• 引張ひずみ下での破壊を研究するために、次近接相互作用を含む1次元周期的鎖モデルを定式化する。
• 線形化QCF作用素を導出し、そのスペクトル的性質を解析し、エネルギーに基づく準非局所（QNL）法と、厳密な関数解析的手法を用いて比較する。
• 重み付きℓp空間における Hölder の不等式およびノルム推定を用いて、線形化QCF作用素の作用素ノルムを系のサイズNに依存しない形で有界化する。
• 作用素の挙動を調査するため、微分係数が区分的定数関数であるように巧みに構築されたテスト関数vを導入し、v′のℓpノルムに対する下界を導出する。
• 数値実験を用いて線形化QCF作用素とQNL作用素のスペクトルを比較し、テストされた範囲では両者のスペクトルが完全に一致することを明らかにする。
• 逆作用素の一様有界性が、関数空間の位相の選択に強く依存することを示し、特定のノルムでは不安定性を示す反例を提示する。

実験結果

リサーチクエスチョン

• RQ1 線形化作用素が非正規的かつ不定的であるにもかかわらず、QCF法は臨界破壊荷重に達するまで動的安定性を示すのか？
• RQ2 系のサイズNが増加する際、線形化QCF作用素の逆作用素が一様に有界となるような離散関数空間（位相）の選択は何か？
• RQ3 QCF作用素のスペクトル的性質は、臨界荷重未満で安定であることが知られているエネルギーに基づくQNL作用素とどのように比較できるか？
• RQ4 線形化作用素が強収縮性（正定性）を欠いているにもかかわらず、QCF法は破壊の発生を信頼性高く予測できるのか？
• RQ5 QCF作用素の非正規性が、その安定性および物理的予測可能性に果たす役割は何か？

主な発見

• 線形化QCF作用素のスペクトルは、線形化QNL作用素と同一であり、QCF法が臨界破壊荷重までQNL法の動的安定性を引き継ぐことを示唆する。
• 非正規構造のため、固有値がすべて正であっても、線形化QCF作用素は一般に不定的であり、正定性が安定性の十分条件でないことを示す。
• 適切に関数空間が選ばれていれば、線形化QCF作用素の逆作用素はℓp(ε)ノルム（p > 1）において一様に有界であり、大規模系の極限でも安定性が保証される。
• 反例により、ℓ1(ε)ノルムでは一様有界性が成立しないことが示され、安定性が関数空間の選択に強く依存することを強調する。
• 線形化QCF作用素の作用素ノルムはC N^{-1/p}で有界であることが示され、条件数が系のサイズNに従ってゆっくりと増加することを支持し、数値的安定性を裏付ける。
• 非退化なℓpノルムを有するテスト関数の構築と、作用素への作用の制御により、ノルム推定の鋭さおよび関数空間選択の必要性が確認される。