[論文レビュー] Finite Element Technique for Solving the Stream Function Form of a Linearized Navier-Stokes Equations Using Argyris Element

Faisal Fairag, Noah Al-Mulla|ArXiv.org|Jun 3, 2004

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

ひとこと要約

本稿では、線形化されたナビエ＝ストークス方程式のストリーム関数形式を解くために、Argyris要素を用いた有限要素法を提示する。行列の帯域幅を最小化するノード再順序付けを通じて計算効率を向上させることに注力している。この手法は最適な収束性と安定性を達成しており、Re=1の単位正方形領域におけるマニュフェイクド正確解を用いた数値結果は、適応的台形則を用いて高精度な解を得ていることを示している。

ABSTRACT

The numerical implementation of finite element discretization method for the stream function formulation of a linearized Navier-Stokes equations is considered. Algorithm 1 is applied using Argyris element. Three global orderings of nodes are selected and registered in order to conclude the best banded structure of matrix and a fluid flow calculation is considered to test a problem which has a known solution. Visualization of global node orderings, matrix sparsity patterns and stream function contours are displayed showing the main features of the flow.

研究の動機と目的

• 線形化されたナビエ＝ストークス方程式のストリーム関数形式に対する安定的かつ効率的な有限要素離散化を開発すること。
• 得られる剛性行列の帯域幅とスパarsity構造に及ぼす異なるグローバルノード順序付けの影響を評価すること。
• 有限要素離散化から生じる非対称かつスパースな線形方程式系を解くために、共役グレディエント安定化法（BiCGSTAB）を実装・テストすること。
• 単位正方形領域におけるマニュフェイクド正確解を用いて、数値スキームの妥当性を検証すること。
• メッシュの細分化や台形則の変更に伴う計算コスト、収束速度、反復ソルバーの性能を分析すること。

提案手法

• ストリーム関数形式により圧力の項が消去され、非圧縮性が自動的に満たされるため、問題は4階の双調和型偏微分方程式に簡略化される。
• 弱形式は、Sobolev空間 $ H^2_0(\Omega) $ を用いて導出され、双線形形式 $ a(\psi,\varphi) $ と $ b(\psi^*;\psi,\varphi) $ がそれぞれラプラシアン項と対流項を表す。
• Argyris要素が $ H^2_0(\Omega) $ 内の整合的有限要素として用いられ、4階問題に必要な $ C^1 $ 連続性を保証する。
• 行列の帯域幅を最小化し、剛性行列の条件数を改善するために、3つのグローバルノード再順序付け戦略をテストした。
• 得られたスパースかつ非対称な線形方程式系に対して、BiCGSTAB反復ソルバーを適用し、残差および反復解のノルムの許容誤差を用いて収束をモニタリングした。
• 数値積分には4点および6点のガウス台形則を用い、要素行列の評価における精度と計算コストを評価した。

実験結果

リサーチクエスチョン

• RQ1 ノード再順序付けは、$ H^2 $-整合離散化における有限要素行列の帯域幅とスパarsityパターンにどのように影響を与えるか？
• RQ2 4点と6点の台形則選択が、Argyris有限要素解の精度と計算コストに及ぼす影響は何か？
• RQ3 線形化問題において、メッシュ細分化（$ h = 1/3, 1/5, 1/9 $）に伴い、BiCGSTAB反復回数とCPU時間はどのように変化するか？
• RQ4 マニュフェイクド正確解を用いたストリーム関数形式において、Argyris有限要素法が最適な収束率を達成できるか？
• RQ5 バイハーモニック方程式を解いて得られる初期推定値 $ \psi^* $ の選択が、反復ソルバーの収束に与える影響は何か？

主な発見

• ノード再順序付け1が最も狭い行列帯域幅と最小のCPU時間（28.45秒）を達成し、計算効率において他の手法を上回った。
• $ h = 1/9 $ の場合、$ L^2 $-ノルムで相対誤差 $ 1.423 \times 10^{-4} $、$ H^1 $-ノルムで $ 6.986 \times 10^{-4} $ を達成し、BiCGSTAB反復回数は567.5回であった。
• 算術演算回数はメッシュ細分化に伴い著しく増加し、$ h=1/3 $ の約1750万回から $ h=1/9 $ の約35億回にまで増加した。これは、$ H^2 $-整合要素の $ O(h^{-4}) $ スケーリングを反映している。
• 6点台形則の使用は4点と比較してわずかに誤差と演算回数を増加させたが、安定性と収束行動の改善に寄与した。
• 流れのストリームライン可視化により、すべてのメッシュサイズで物理的に整合した流れのパターンが確認され、この手法が渦構造を正しく捉えられることを裏付けた。
• 全テストされたメッシュサイズおよび台形則において、収束が安定しており、残差および反復解のノルムが許容誤差 $ 10^{-5} $ 未満に低下した。