[論文レビュー] Iterative solution of piecewise linear systems for the numerical solution of obstacle problems

Luigi Brugnano, Alessandra Sestini|arXiv (Cornell University)|Dec 16, 2009

Matrix Theory and Algorithms参考文献 17被引用数 24

ひとこと要約

本稿では、線形障害問題の数値解法およびその放物型対応から生じる分岐線形系（PLS）を解くための半反復的ニュートン型手法を提案する。この手法は、メッシュの細分化に依存せずに、有限ステップ内で正確な解へのグローバルな単調収束を保証する。有限差分法および有限要素法による離散化を用いた楕円型および放物型障害問題における数値実験により、その有効性が示された。

ABSTRACT

We investigate the use of piecewise linear systems, whose coefficient matrix is a piecewise constant function of the solution itself. Such systems arise, for example, from the numerical solution of linear complementarity problems and in the numerical solution of free-surface problems. In particular, we here study their application to the numerical solution of both the (linear) parabolic obstacle problem and the obstacle problem. We propose a class of effective semi-iterative Newton-type methods to find the exact solution of such piecewise linear systems. We prove that the semiiterative Newton-type methods have a global monotonic convergence property, i.e., the iterates converge monotonically to the exact solution in a finite number of steps. Numerical examples are presented to demonstrate the effectiveness of the proposed methods.

研究の動機と目的

• 障害問題の数値解法に生じる分岐線形系（PLS）を対象とする、効率的な反復型ソルバーの開発を目的とする。
• 反復スキームが正確な解へ有限ステップ内でグローバルに単調収束することを保証することを目的とする。
• 線形楕円型（古典的）障害問題およびその放物型時間発展対応への応用を目的とする。
• 数値実験を通じて、メッシュの細分化に依存しない収束特性を示すこと。

提案手法

• 本手法は、解に依存する係数行列が分岐定数関数である分岐線形系に適用する半反復的ニュートン型アルゴリズムに基づく。
• 各反復ステップで、現在の反復解に基づいて活性集合（解が障害物に接触する領域）を更新し、線形部分問題の系列を生成する。
• 障害問題の補完性形式とM行列の構造を活用することで、反復解が明確に定義されることを保証する。
• 解法プロセスは双対活性集合戦略に相当するが、明確さと効率性を向上させるために、直接的にPLSの観点から定式化されている。
• 放物型の場合には、時間方向に陰的エイラー離散化を用い、各時間ステップでPLSの系列を解く。
• 収束が単調的かつ有限ステップで成立することが証明されており、初期推定値やメッシュ細分化に依存しない。

実験結果

リサーチクエスチョン

• RQ1 障害問題から導かれる分岐線形系に、半反復的ニュートン型手法を効果的に適用可能か？
• RQ2 提案手法が有限ステップ内で正確な解へのグローバルな単調収束を保証するか？
• RQ3 楕円型および放物型障害問題において、収束速度が空間的・時間的メッシュ細分化に依存しないか？
• RQ4 実装上、問題のサイズおよび障害パラメータに応じて反復回数はどのように変化するか？

主な発見

• 提案された半反復的ニュートン型手法は、分岐線形系の正確な解へ有限ステップ内でグローバルに単調収束を達成する。
• 線形楕円型障害問題において、境界条件にかかわらず、全テストメッシュサイズ（N=25から100）で5〜8ステップの収束が達成された。
• τ=10⁴およびν=20の時間ステップを有する放物型障害問題では、各時間ステップにおける反復回数が5〜8に一定であり、空間解像度に依存しなかった。
• 弾塑性ねじり問題の変種では、障害パラメータCに応じて各時間ステップの反復回数が9〜32に変動したが、メッシュサイズに依存せずに安定した。
• 境界条件がディリクレ型とノイマン型の間で、一致集合（解が障害物に接触する領域）はほぼ同一であり、境界付近でのわずかな差異を除いて一致した。
• 収束特性はメッシュ細分化に伴い劣化せず、適応的および高解像度のシミュレーションに適していることが示された。