[論文レビュー] Adaptive methods for PDE's: wavelets or mesh refinement?

Albert Cohen|ArXiv.org|Dec 1, 2002

Advanced Numerical Methods in Computational Mathematics参考文献 28被引用数 15

ひとこと要約

この論文は、偏微分方程式（PDE）を解くための適応的ウェーブレット法と適応的メッシュ・リファインメントを比較し、非線形近似理論に基づく理論的最適収束率を達成するが、幾何的複雑性、データ構造のオーバーヘッド、および非等方的解像度の低さのため、実用的な制限に直面することを示している。強固な理論的基盤を持つにもかかわらず、多くの問題において有限要素法に比して計算効率が劣る。

ABSTRACT

Adaptive mesh refinement techniques are nowadays an established and powerful tool for the numerical discretization of PDE's. In recent years, wavelet bases have been proposed as an alternative to these techniques. The main motivation for the use of such bases in this context is their good performances in data compression and the approximation theoretic foundations which allow to analyze and optimize these performances. We shall discuss these theoretical foundations, as well as one of the approaches which has been followed in developing efficient adaptive wavelet solvers. We shall also discuss the similarities and differences between wavelet methods and adaptive mesh refinement.

研究の動機と目的

• 特異性を有するPDEを解くための適応的ウェーブレット法と適応的メッシュリファインメントの理論的・実用的性能を比較すること。
• 特に関数空間における非線形近似理論を用いて、両手法の収束率を分析すること。
• 実世界の応用においてウェーブレットに基づく適応的ソルバーに内在する根本的制限（「呪い」）を特定すること。
• ウェーブレット法が実際の計算で最適収束率を達成できるかどうか、特に有限要素法と比較して評価すること。
• 幾何的複雑性、データ構造管理、および非等方性がウェーブレットに基づく適応的ソルバーの効率に与える影響を評価すること。

提案手法

• 非線形近似理論を用いて、ウェーブレットと有限要素法の最良近似誤差の減少率を比較し、誤差はソボレフ空間または$L^2$ノルムで測定する。
• 関数に特異性（不連続性や境界層など）が存在する場合のウェーブレット係数の減少を定量的に評価するため、$N$-項近似の概念を適用する。
• 動的基底関数選択（適応的集合$\Lambda_n$）とスパース表現を用いた行列・ベクトル演算に基づく適応的ウェーブレットアルゴリズムを構築する。
• ウェーブレット・ガレルキン行列のスパarsityを活用して、線形方程式系を効率的に解くための行列・ベクトルアルゴリズム（式24）を採用する。
• 理論的収束境界を維持するための粗化ステップをアルゴリズムに組み込むが、実際には収束率に影響を及げないため、しばしば省略される。
• 後験的誤差推定と階層的メッシュリファインメントを用いて、ウェーブレット法と適応的有限要素法を比較する。

実験結果

リサーチクエスチョン

• RQ1 特異性を有するPDEに対して、ウェーブレットに基づく適応的手法は最適収束率を達成できるか。また、適応的有限要素法と比較してどうなるか。
• RQ2 ウェーブレット法が最適収束を達成できる理論的基盤は何か。非線形近似理論とどのように関連しているか。
• RQ3 最適収束率を達成するにもかかわらず、なぜウェーブレットベースのソルバーは実用的に適応的有限要素法に劣るのか。
• RQ4 幾何的複雑性、データ構造管理、非等方性は、ウェーブレットに基づく適応的ソルバーの性能にどのように影響を与えるか。
• RQ5 「呪い」（幾何、データ構造、非等方性）が、PDEソルバーにおけるウェーブレット法の実用的適用性をどの程度制限するか。

主な発見

• ウェーブレットに基づく適応的手法は、特異性を有する特定のクラスのPDEに対して、非線形近似理論に基づく強固な理論的基盤により最適収束率を達成する。
• 滑らかな領域の特性関数に対する$L^2$空間における最良$N$-項近似レートは、$r = 1/(2d-2)$に制限され、次元が増加するにつれて劣化する。これは、幾何的不連続性を扱う際の根本的制限を示している。
• 理論的最適性にもかかわらず、係数およびインデックスの動的データ構造管理に起因する高い計算オーバーヘッドにより、実際には有限要素法に比して遅くなる。
• 幾何の呪いは、非長方形領域におけるウェーブレット基底が、ドメイン分割などの複雑な構成を必要とし、テンソル積系の単純さを失うことに起因する。
• 非等方性の呪いは、ウェーブレット法が等方的リファインメントしか行えないため、ショック波や境界層のような方向的特徴を効率的に解像できないことに起因する。
• 収束の理論的証明に不可欠な粗化ステップは、実際には収束率に影響を与えず省略されがちであり、理論と実装の間にはギャップが存在する。