[論文レビュー] Lagrangian basis method for dimensionality reduction of convection dominated nonlinear flows

Rambod Mojgani, Maciej Balajewicz|arXiv (Cornell University)|Jan 16, 2017

Model Reduction and Neural Networks参考文献 26被引用数 37

ひとこと要約

本稿では、移動する衝撃波や急峻な勾配を示す波動的解を、ラグランジュ座標系で計算領域を変化させることで、わずかな数のグローバル基底関数のみで効率的に圧縮できるラグランジュ基底法を提案する。この手法は、このような解に内在する低ランク構造を捉えることができ、スカラー移流拡散方程式および衝撃波を伴うエーラー方程式の数値実験により、従来のオイラー座標系手法を凌駕することを示している。

ABSTRACT

Foundations of a new projection-based model reduction approach for convection dominated nonlinear fluid flows are summarized. In this method the evolution of the flow is approximated in the Lagrangian frame of reference. Global basis functions are used to approximate both the state and the position of the Lagrangian computational domain. It is demonstrated that in this framework, certain wave-like solutions exhibit low-rank structure and thus, can be efficiently compressed using relatively few global basis. The proposed approach is successfully demonstrated for the reduction of several simple but representative problems.

研究の動機と目的

• 移動する衝撃波、不連続性、急峻な勾配を示す解の圧縮に失敗する従来の射影ベースのモデル次元削減（MOR）の問題に対処すること。
• 波動的ダイナミクスを扱いながらもグローバル基底構造を保持する一般的かつ単純な次元削減フレームワークを構築すること。
• ラグランジュ座標系において波動的解が低ランク構造を示すことを実証し、最小限の基底関数で効率的な圧縮が可能であることを示すこと。
• オイラー座標系手法がスペクトル的非効率性のため、複雑な移流優勢な流れに応用できない状況において、MORを拡張する基盤を提供すること。

提案手法

• 本手法は、状態および計算グリッドの位置の両方をグローバル基底関数で近似することで、ラグランジュ座標系における流れの時間発展を定式化する。
• ラグランジュ座標系における解スナップショットの低ランク近似を用い、式 (3.5) で示される修正されたランク制約付き最適化問題を解くことで基底関数を構築する。
• 状態およびグリッド位置は、それぞれグローバル基底行列 Uw および Ux が張る低次元部分空間に射影され、これにより低次元モデル化が可能になる。
• 本手法はラグランジュ座標系におけるガラーキン射影を用い、移動グリッドに従って進化する低次元系を導出する。
• 本手法は使用するハイパーリダクション技術に依存しないため、さらなる計算速度向上のための既存のハイパーリダクション戦略と互換性を持つ。
• 適切な直交分解（POD）を基底生成に用いるが、DMD やクープマンモードなど他の手法にも適用可能である。

実験結果

リサーチクエスチョン

• RQ1 移動する衝撃波および急峻な勾配を示す波動的解は、ラグランジュフレームワークにおいてグローバル基底関数を用いて効率的に圧縮可能か？
• RQ2 なぜ従来のオイラー座標系射影ベースMOR手法はこのような解に対して失敗するのか？また、座標系の参照フレームを変えることでこの失敗を克服できるか？
• RQ3 ラグランジュ定式化により、オイラー座標系では圧縮不能であった解に低ランク構造が顕在化するか？
• RQ4 少数のグローバル基底関数で、移流優勢な流れにおける複雑な衝撃ダイナミクスを正確に再現できるか？
• RQ5 本手法は、高次元またはより複雑な流体力学的問題へ一般化可能か？

主な発見

• Re = 1000 のスカラーバーガーズ方程式において、k = 2 のラグランジュROMは高精度モデル（HFM）と区別できない解を得たが、オイラー座標系ROMは衝撃近傍で大きな振動を示した。
• 断面積が変化するノズル流れ問題において、µ10 の条件下で k = 2 のラグランジュROMは衝撃位置と密度プロファイルを正確に捉えたが、オイラー座標系ROMは顕著な誤差と振動を示した。
• 移流拡散問題において、ラグランジュROMは k = 5 時に相対誤差が 10−3 未満に抑えられたのに対し、オイラー座標系ROMは同程度の精度を得るためには k = 15 必要だった。
• 本手法は、移動する不連続性を伴う解がラグランジュ座標系において低ランク構造を示すことを実証し、わずかな数のグローバル基底関数で効率的な圧縮が可能であることを示した。
• 誤差減衰曲線（図4）により、衝撃を伴う問題において、従来のオイラー座標系MORに比べ、収束性と安定性の面で著しく優れた性能を示した。
• ラグランジュフレームワークにより、k = 2 の基底関数で µ1 から µ10 までの複数のパrameterインスタンスにおいて、全パrametric変動を正確に再現できた。