[論文レビュー] Truncated Nonsmooth Newton Multigrid for phase-field brittle-fracture problems.

Carsten Gräser, Daniel Kienle+1 more|arXiv (Cornell University)|Jul 23, 2020

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

ひとこと要約

本稿では、相場分野の破壊問題における空間的離散化を効率的に解くために、切り捨てられた滑らかでないニュートン多重グリッド法（TNNMG）を導入する。変分構造を活用し、罰則法や履歴場を用いずに破壊の不可逆制約を直接扱うことで、TNNMGは1回の線形解法程度の時間で収束を達成し、従来のオペレータ分割法よりも10倍以上高速でありながらも、安定性を維持する。

ABSTRACT

We propose the Truncated Nonsmooth Newton Multigrid Method (TNNMG) as a solver for the spatial problems of the small-strain brittle-fracture phase-field equations. TNNMG is a nonsmooth multigrid method that can solve biconvex, block-separably nonsmooth minimization problems in roughly the time of solving one linear system of equations. It exploits the variational structure inherent in the problem, and handles the pointwise irreversibility constraint on the damage variable directly, without penalization or the introduction of a local history field. Memory consumption is significantly lower compared to approaches based on direct solvers. In the paper we introduce the method and show how it can be applied to several established models of phase-field brittle fracture. We then prove convergence of the solver to a solution of the nonsmooth Euler-Lagrange equations of the spatial problem for any load and initial iterate. Numerical comparisons to an operator-splitting algorithm show a speed increase of more than one order of magnitude, without loss of robustness.

研究の動機と目的

• 微小変形の相場分野破壊モデルに生じる空間的問題に対して、高速で、安定的かつメモリ効率の良い解法を開発すること。
• 罰則法や補助履歴場を用いずに、損傷変数の局所的不可逆制約を直接強制すること。
• 特にオペレータ分割に基づく反復解法と比較して、計算コストを低減すること。
• さまざまな荷重条件や初期推定値に対しても安定性を維持すること。
• 任意の荷重および初期推定値に対して、滑らかでないオイラー＝ラグランジュ方程式の解に解法が収束することを証明すること。

提案手法

• TNNMG法は、双凸でブロック別非滑らか最小化問題を対象とした非滑らか多重グリッド解法である。
• 相場分野破壊エネルギー関数の変分構造を活用して、安定的かつ効率的な解法プロセスを構築する。
• 最適化フレームワーク内で損傷変数の局所的不可逆制約を直接処理し、正則化や補助場を避ける。
• 粗いグリッド上で切り捨てられたニュートンステップを用い、多重グリッドサイクルを用いてスケールを跨いで補正を効率的に伝搬する。
• 非滑らか性の影響を受けても収束を保証するため、非単調ラインサーチを採用する。
• TNNMG法は、アムブロシオ＝トルトレッリおよびミヘーらのモデルを含む、既に確立された相場分野破壊モデルに適用される。

実験結果

リサーチクエスチョン

• RQ1 物理的破壊の不可逆性を保ちつつ、空間的相場分野破壊問題を効率的に解く非滑らか多重グリッド法を設計できるか？
• RQ2 TNNMG法は、相場分野破壊問題における標準的なオペレータ分割アルゴリズムと比較して、速度と安定性で優れているか？
• RQ3 さまざまな荷重や初期推定値の下でのTNNMG法の収束挙動はいかなるものか？
• RQ4 メモリ消費量を増加させず、安定性を損なわず、顕著な高速化を達成できるか？
• RQ5 罰則法とは異なり、不可逆性を直接強制することで、数値的安定性や精度が向上するか？

主な発見

• TNNMG法は、任意の荷重および初期推定値に対して、滑らかでないオイラー＝ラグランジュ方程式の解に収束することが保証され、グローバル収束を示す。
• 解法時間は1回の線形システムを解くコストにほぼ等しく、従来の反復解法よりも顕著に高速である。
• 数値的比較により、オペレータ分割法と比較して10倍以上の高速化が確認された。
• 多様な荷重条件や初期推定値に対しても、安定性が維持されている。
• 反復的で多重グリッドに基づく手法であるため、直接解法と比較して、メモリ消費量が著しく低い。
• 不可逆制約の直接処理により、罰則項や履歴場の必要がなくなり、定式化が簡素化され、物理的整合性が向上した。