[論文レビュー] A Multilevel Correction Scheme for Nonsymmetric Eigenvalue Problems by Finite Element Methods

Hehu Xie, Zhimin Zhang|arXiv (Cornell University)|May 23, 2015

Advanced Numerical Methods in Computational Mathematics参考文献 22被引用数 8

ひとこと要約

本稿では、有限要素法を用いて非対称固有値問題に対する多層補正スキームを提案し、固有値の上昇が1であることを仮定せず、効率的かつ高精度な近似を可能にする。粗いメッシュ上でのみ小さな固有値問題を解き、細かい空間上での源問題を解くことで、最適収束率を達成し、特に多重グリッドまたは適応的グリッド・リファインメントと組み合わせることで計算効率を著しく向上させる。

ABSTRACT

A multilevel correction scheme is proposed to solve defective and nodefective of nonsymmetric partial differential operators by the finite element method. The method includes multi correction steps in a sequence of finite element spaces. In each correction step, we only need to solve two source problems on a finer finite element space and two eigenvalue problems on the coarsest finite element space. The accuracy of the eigenpair approximation is improved after each correction step. This correction scheme improves overall efficiency of the finite element method in solving nonsymmetric eigenvalue problems.

研究の動機と目的

• 流体動力学および安定性解析に生じる非対称固有値問題を解くための効率的な有限要素法の開発。
• 従来、対称問題にのみ適用されてきた多層補正技術を、固有値の上昇が1であることを仮定しない非対称作用素へと拡張すること。
• 細かい有限要素空間における大規模固有値問題の解法コストを低減することで、計算効率の向上。
• 多重グリッド法および適応的有限要素法とシームレスに統合し、最適収束性と並列処理を実現すること。

提案手法

• 解像度を段階的に高めた入れ子の有限要素空間の系列を用いる。
• 各補正ステップでは、前回の固有関数を荷重ベクトルとして用い、拡張された有限要素空間上で源問題を解く。
• 次に、補正済み固有関数を含む空間上で新しい固有値問題を解き、固有値計算には唯一粗い空間を使用する。
• 繰り返し補正ステップを通じて、固有値および固有関数の両方の成分を段階的に改善する。
• 細かいグリッド上で完全な固有値問題を解く代わりに、細かい空間上での安価な源問題と粗い空間上での小さな固有値問題に置き換えることで、全問題を回避する。
• 境界値問題の多重グリッドソルバーおよびZZ回復に基づく後騒誤差推定子を用いた適応的リファインメントと両立可能である。

実験結果

リサーチクエスチョン

• RQ1 固有値の上昇が1であることを仮定しない非対称固有値問題に適した多層補正スキームを設計できるか？
• RQ2 有限要素法を用いて大規模非対称固有値問題の計算コストをどのように低減できるか？
• RQ3 特異固有関数（例：L字型領域上）を有する問題に適用した場合、多層補正スキームが最適収束率を維持できるか？
• RQ4 微分作用素に複素ベクトル場が含まれる場合、この手法はどのように動作するか？
• RQ5 多重グリッドおよび適応的有限要素技術と効果的に組み合わせられ、計算効率が向上するか？

主な発見

• 特異固有関数を有する非凸領域（例：L字型領域）上の問題に対しても、固有値および固有関数の近似において最適収束率を達成する。
• 数値結果により、準均一メッシュ上での最初の固有値に対して、収束率が理論的期待通りの $h^{4/3}$ のオーダーを維持することが確認された。
• 適応的メッシュ上でも、誤差が期待される割合で減少し、初期メッシュの選択が最終的な精度にほとんど影響しないことが示された。
• 微分作用素に複素ベクトル場が含まれる場合でも、最初の6つの固有値に対して急速な収束が観察された。
• 複数の空間を用いるアプローチは、最初の6つの固有値の和を計算する際、多重グリッド法を上回る高い精度を達成しており、高精度な設定において優位性を示している。
• 本スキームは非常に効率的である：粗い空間上での小さな固有値問題と細かい空間上での源問題を解くだけで、細かい空間上での完全な固有値問題を解いたのと同等の精度が得られる。