[論文レビュー] On the analysis of block smoothers for saddle point problems

Lorenz John, Ulrich Rüde+2 more|arXiv (Cornell University)|Dec 5, 2016

Advanced Numerical Methods in Computational Mathematics参考文献 21被引用数 3

ひとこと要約

本稿では、鞍点問題に適用された多重グリッド法におけるさまざまなUzawa型ブロック smoother の滑らかさ特性について統一的な分析を提示する。特に、ストークス系を対象として、新しい対称的Uzawa smoother を導入し、不正確およびブロック因子分解型の変種について収束性を厳密に証明している。有限要素離散化における数値的検証を通じて、多重グリッドサイクル全体にわたる頑健な性能を示している。

ABSTRACT

In this article, we discuss several classes of Uzawa smoothers for the application in multigrid methods in the context of saddle point problems. Beside commonly used variants, such as the inexact and block factorization version, we also introduce a new symmetric method, belonging to the class of Uzawa smoothers. For these variants we unify the analysis of the smoothing properties, which is an important part in the multigrid convergence theory. These methods are applied to the Stokes problem for which all smoothers are implemented as pointwise relaxation methods. Several numerical examples illustrate the theoretical results.

研究の動機と目的

• 鞍点問題に対する多重グリッド法における複数のUzawa型 smoother の滑らかさ特性を統一的に分析すること。
• 広く用いられているが理論的分析が少ない不正確Uzawa法の滑らかさ特性を証明すること。
• 不正確Uzawa法の対称化変種としての新規対称的Uzawa smoother を提案し、その解析を行うこと。
• 点緩和スキームを用いてストークス問題に対して理論的結果を数値的に検証すること。
• VサイクルおよびWサイクルの両方の多重グリッド設定における、さまざまなsmoother の計算効率と頑健性を比較すること。

提案手法

• 著者たちは、Uzawa反復に基づくブロック smoother を分析し、多重グリッド緩和における逆演算を近似するためのブロック条件付け行列を用いる。
• 複数のUzawa変種（不正確およびブロック因子分解型を含む）について、多重グリッド収束理論の主要な構成要素である滑らかさ特性を統一的な枠組みで導出し、適用する。
• 反復処理における安定性と対称性を向上させる、新しい対称的Uzawa smoother を提案する。
• すべてのsmoother は、分散メモリシステムとの互換性を保つために、点緩和法（例：ヤコビ法や Gauss-Seidel 法）として実装される。
• 解析は、スレッダー補完子のスペクトル同等性と性質に依存しており、行列ブロックの対称性および半正定値性を仮定する。
• 数値実験は、複数のメッシュレベルにわたる有限要素離散化を用いてストークス問題に対して実施された。

実験結果

リサーチクエスチョン

• RQ1 鞍点問題に対する多重グリッド法におけるさまざまなUzawa型smoother の滑らかさ特性を統一的に分析することは可能か？
• RQ2 新たに提案された対称的Uzawa smoother は、多重グリッド収束のための必要な滑らかさ特性を示すか？
• RQ3 圧力スレッダー補完子の近似法の選択（質量行列対スケーリング済みラプラシアン）が、ソルバー性能にどのように影響するか？
• RQ4 不正確Uzawa smoother はVサイクルおよびWサイクルの両方の多重グリッドスキームにおいて、頑健で効率的か？
• RQ5 理論的滑らかさ特性は、複数のメッシュ細分化レベルにわたって数値的に検証可能か？

主な発見

• 不正確Uzawa smoother が滑らかさ特性を有することを証明し、多重グリッド収束における理論的妥当性を確認した。
• 新規対称的Uzawa smoother は、特にWサイクル多重グリッド応用において、より優れた安定性と頑健性を示した。
• 数値結果から、スレッダー補完子へのスペクトル同等性が欠けるにもかかわらず、質量行列の代わりにスケーリング済み離散圧力ラプラシアンを用いることで、収束特性が著しく向上することが示された。
• Vサイクルでは、不正確Uzawa smoother の収束にはより多くの緩和ステップが必要であり、十分な反復回数が得られない場合には発散することが判明した。これはパrameter 選択への感受性を示している。
• Wサイクルは、大規模問題に対してVサイクルよりも優れた反復回数とソリューション時間を得ており、特に並列環境下では順次処理のボトルネックが生じるが、全体として優れた性能を示した。
• 粗いレベルで緩和ステップ数を増加させる可変Vサイクルは、大規模並列システムにおける頑健性と性能を向上させ、高性能コンputィング応用において実用的利点を示唆した。