[論文レビュー] Space-time discretization of the heat equation. A concise Matlab implementation

Roman Andreev|arXiv (Cornell University)|Dec 25, 2012

Advanced Numerical Methods in Computational Mathematics参考文献 12被引用数 6

ひとこと要約

本稿では、熱方程式に対する空間時間ガラーキン有限要素法の、簡潔で安定的かつ並列化可能なMatlab実装を提示する。この手法により、反復的に解く1つの線形系を通じて、空間時間の同時離散化が可能となる。アプローチは、Kronecker積構造と効率的な前処理を利用しており、非一様な時間的メッシュや異なる空間離散化においても最適収束性と頑健な性能を達成する。

ABSTRACT

A concise Matlab implementation of a stable parallelizable space-time Petrov-Galerkin discretization for parabolic evolution equations is given. Emphasis is on reusability of spatial finite element codes.

研究の動機と目的

• 放物型発展方程式、特に熱方程式のためのモジュラーで安定的かつ並列化可能な空間時間離散化手法を開発すること。
• 空間と時間の離散化を統一された変分枠組み内で分離することで、既存の空間有限要素コードの再利用を可能にすること。
• 非一様な時間グリッド、時間依存入力、グローバルな精度制御を備えた反復的解法をサポートする、コンactなMatlab実装を提供すること。
• 構造的および非構造的空間領域における数値実験を通じて、手法の効率性と頑健性を示すこと。

提案手法

• 空間時間領域上で熱方程式の弱形式をテストすることで、空間時間変分形式を導出する。初期条件と境界条件を1つの関数的設定に統合する。
• 試行関数と検証関数の空間時間におけるテンソル積有限要素基底を用いたペトロフ・ガラーキン枠組みを採用し、安定性と準最適性を保証する。
• ノルム誘導作用素を導入して離散解のノルムを定義し、その行列表現をKronecker積構造を用いて計算効率よく構築する。
• 得られた線形系は、一般化されたLSQRアルゴリズムを用いて反復的に解き、ノルム誘導作用素行列の逆行列を前処理として用いる。
• 前処理は、システム行列のKronecker構造を活用し、効率的な行列・ベクトル積を実現するとともに、条件数を低減する。
• 実装はモジュラーである。空間FEMコードは空間離散化を組み込むことで再利用可能であり、時間離散化は時間的有限要素を通じて独立して処理される。

実験結果

リサーチクエスチョン

• RQ1 既存の空間FEMコードとモジュラーに統合可能な、数行のMatlabコードで実装可能な安定的かつ並列化可能な空間時間有限要素法を構築できるか？
• RQ2 ペトロフ・ガラーキン枠組みにおける、時間的メッシュの細分化と検証関数空間の選択が、反復解法の条件数および反復回数に与える影響はいかほどか？
• RQ3 マルチスレッドアーキテクチャ上でのスケーラビリティを考慮した場合、Crank-Nicolsonのような従来の時間ステッピングスキームに比べ、本手法は解法時間とスケーラビリティにおいてどの程度優れているか？
• RQ4 非一様な時間的メッシュと時間依存入力を扱えるが、安定性や収束性を損なわずに行えるか？

主な発見

• タイプ1の時間的部分空間（検証関数空間の細分化なし）では、一般化LSQRの反復回数は時間的メッシュサイズの増加に伴い一時的に増加するが、条件数の低下により最終的に減少する。
• タイプ2の時間的部分空間（検証関数空間の細分化あり）では、メッシュサイズにかかわらず、前処理済みシステムの条件数は約2でほぼ一定となり、反復回数は一貫して低い水準を維持する。
• 本手法は空間時間ノルムにおいて最適収束性を達成し、低次の時間的有限要素を用いても、高次の時間ステッピングスキームと同等の精度を実現する。
• 実行時間の観点では競争力があるが、マルチグリッドやウェーブレットベースの前処理と組み合わせることで、大規模並列アーキテクチャを最大限に活用できるように設計されている。
• 非一様メッシュに対して非局所変換を用いることで、並列環境での通信オーバーヘッドが生じるため、スケーラビリティの観点から、均一またはダイアディックメッシュの使用が好ましいと示唆される。