[論文レビュー] Proper orthogonal decomposition closure models for fluid flows: Burgers equation

Omer San, Traian Iliescu|arXiv (Cornell University)|Aug 14, 2013

Model Reduction and Neural Networks参考文献 36被引用数 67

ひとこと要約

本稿では、移動する衝撃波を有する1次元のバーガース方程式をテストベッドとして用い、流体の流れのための固有直交分解低次元モデル（POD-ROM）の10種類の閉じ込めモデルを提案および評価する。新規の渦粘性およびエネルギー保存型モデルを含むこれらのモデルは、標準的なガレルキンPOD-ROM-Gよりも顕著に精度を向上させつつ、計算コストを低く維持しており、特にパrameterチューニング下でPOD-ROM-RおよびPOD-ROM-RQが最も優れた全体的な性能を示している。

ABSTRACT

This paper puts forth several closure models for the proper orthogonal decomposition (POD) reduced order modeling of fluid flows. These new closure models, together with other standard closure models, are investigated in the numerical simulation of the Burgers equation. This simplified setting represents just the first step in the investigation of the new closure models. It allows a thorough assessment of the performance of the new models, including a parameter sensitivity study. Two challenging test problems displaying moving shock waves are chosen in the numerical investigation. The closure models and a standard Galerkin POD reduced order model are benchmarked against the fine resolution numerical simulation. Both numerical accuracy and computational efficiency are used to assess the performance of the models.

研究の動機と目的

• 流体の流れにおけるPOD-ROMのための新しい閉じ込めモデルの開発および評価、特に標準的なPOD-ROM-Gが失敗する乱流または衝撃を伴う系において。
• 複雑な3次元流れへの応用の前に、挑戦的ではあるが簡略化された設定—移動する衝撃波を有する1次元バーガース方程式—におけるこれらのモデルの性能を評価すること。
• 精度、計算効率、パrameter感受性に基づいて、閉じ込めモデルを包括的に比較すること。
• 流体動力学のPOD-ROMにおいて、精度と計算コストの間の最良のトレードオフを提供する閉じ込めモデルを特定すること。
• 将来的な研究の基盤を築くこと—特に、スナップショットトレーニング期間を超えた動的パrameter選択および予測能力について。

提案手法

• 3つのクラスの閉じ込めモデルを提案：(1) PODモードに依存する係数を有する定数渦粘性モデル；(2) スマゴリンスキーLESにインspiredされた動的渦粘性モデル；(3) 自由パrameterを回避するエネルギー保存型モデル。
• 移動する衝撃波を模擬できる1次元バーガース方程式をテスト問題として使用し、閉じ込めモデルの性能を制御された形で評価可能にする。
• 高精度なDNSスナップショットからPODモードを生成し、ガレルキン射影を用いて支配方程式を低次元POD部分空間上に射影する。
• 一貫した比較が可能な統一フレームワーク内に、10の閉じ込めモデル（うち3つが新規）を導入する。
• 自由パrameter（例：渦粘性係数）のパrameter感受性を調査し、ロバストネスおよび最適チューニングを評価する。
• 誤差ノルムおよび計算コストの観点から、すべてのモデルを細分解のDNSおよび標準的なPOD-ROM-Gと比較する。

実験結果

リサーチクエスチョン

• RQ1 移動する衝撃波を有するバーガース方程式のダイナミクスを近似する際、どの閉じ込めモデルが最も高い精度を示すか？
• RQ2 各閉じ込めモデルの計算コストは、標準的なPOD-ROM-GおよびDNSと比較してどの程度か？
• RQ3 自由パrameterの選択にどの程度感受性を示し、最適なパラメータ範囲は何か？
• RQ4 より複雑な閉じ込めモデル（例：動的またはエネルギー保存型）は、単純な定数係数モデルよりも顕著に優れた結果をもたらすか？
• RQ5 1次元バーガース方程式での閉じ込めモデルの性能は、より複雑な3次元乱流流れへ一般化可能か？

主な発見

• 10のすべての閉じ込めモデルが、標準的なPOD-ROM-Gよりも顕著に精度を向上させている。POD-ROM-Gは衝撃ダイナミクスを適切に捉えられていない。
• POD-ROM-RおよびPOD-ROM-RQモデルは、2つのテスト問題における誤差ノルムの最小化において、他のモデルを一貫して上回っている。
• POD係数の時間発展を追跡した場合、POD-ROM-RおよびPOD-ROM-RQが再び優れた性能を示したが、他のいくつかのモデルも競合可能であった。
• すべての閉じ込めモデルが、POD-ROM-Gと同等の計算コスト（1実験あたり約4秒）を達成しており、これはDNS（95〜130秒）と比べて顕著に低い。
• エネルギー保存型モデル（POD-ROM-C）および動的スマゴリンスキー型モデル（POD-ROM-S）は、より高い計算オーバーヘッドにもかかわらず、精度が向上しなかった。これは、このテストケースでは複雑さに利益がないことを示唆している。
• 本研究は、パrameterを丁寧にチューニングした場合、PODモードに依存する定数渦粘性モデルがバーガース方程式に対して最も効果的であると結論づけている。