[論文レビュー] Lyapunov stable Galerkin models of post-transient incompressible flows

Maciej Balajewicz|arXiv (Cornell University)|Dec 1, 2013

Model Reduction and Neural Networks参考文献 9被引用数 7

ひとこと要約

本稿では、一般化されたガラーキン近似を用いて反復的に展開モードを最適化することで、非圧縮性流れの定常状態以降の、証明可能にリャプノフ安定なガラーキンモデルを構築する新規手法を提案する。この手法により、線形化されたシステム行列のエルミート部が消去され、固定点の近傍で有界な軌道が保証され、混合層およびリードド・キャビティの両方で検証され、人工的散逸を用いずにナビエ＝ストークス吸引子のダイナミクスを高精度で再現する。

ABSTRACT

A method for deriving provably stable low-dimensional Galerkin models of post-transient incompressible flows is introduced. The proposed approach involves an iterative procedure for expansion modes that satisfy Lyapunov stability in the neighborhood of a fixed point. The approach is demonstrated using two prototypical flow configurations: a two-dimensional mixing layer, and two-dimensional flow inside a lid-driven cavity. This new methodology can be a building block in an effort to develop more accurate and more robust low-dimensional models of incompressible flows.

研究の動機と目的

• 従来の低次元POD-ガラーキンモデルが長期間にわたり有界性を保つために経験的修正を要するという不安定性を解消すること。
• 線形化されたシステム行列のエルミート部が消去されることで、ガラーキンモデルの固定点のリャプノフ安定性を保証する体系的な手法を開発すること。
• このような安定モデルが、人工的散逸を用いずにナビエ＝ストークス吸引子の長期的ダイナミクスを正確に捉えられることを示すこと。
• 軌道が吸引子から逸脱した場合にのみ作動する経験的散逸項を将来的に統合可能にし、内在的な安定性を維持すること。
• モード導出アルゴリズムに依存しないフレームワークを提供することで、POD以外の手法（DMDなど）への拡張を可能にすること。

提案手法

• 展開モードをより多くのPODモード（N > n）の線形結合として表現し、直交変換行列 X ∈ ℝ^{N×n} を用いる一般化ガラーキン近似を採用する。
• ガラーキンテンソル（C, L, Q）をXの関数として再表現することで、システムの線形化ダイナミクスの最適化を可能にする。
• リャプノフ候補関数 V(z) = ½∑zᵢ² を用い、線形化行列Aのエルミート部 A_H = (A + Aᵀ)/2 が恒等的にゼロである場合に安定性が保証される。
• Xの直交性を満たす制約のもとで、A_Hのフロベニウスノルムのゼロからの逸脱を最小化する最適化問題を定式化し、MATLABの非線形最小二乗ソルバを用いる。
• 有限差分によるヤコビアン近似を用いて、fsolveを介したXの制約付き最適化により、A_H = 0 を強制する。
• 得られたモデルは初期条件のエネルギーを引き継ぎ、吸引子上では非作用となる経験的散逸項を追加可能であり、安定性を保持する。

実験結果

リサーチクエスチョン

• RQ1 定常状態以降の非圧縮性流れのガラーキンモデルを、固定点がリャプノフ安定であり、すべての時間で有界な軌道を保証できるように構築できるか？
• RQ2 一般化ガラーキン近似を用いて、線形化システム行列のエルミート部が消えるように流体モードを最適化する方法は何か？
• RQ3 リャプノフ安定なガラーキンモデルは、人工的散逸を用いずにナビエ＝ストークス吸引子の長期的ダイナミクスをどの程度正確に再現できるか？
• RQ4 提案手法は、軌道が吸引子から逸脱した場合にのみ作動する経験的散逸項を統合できるか？
• RQ5 本手法は、混合層やリードド・キャビティなどの異なる流れ状態においても、頑健かつ一般化可能か？

主な発見

• 提案手法により、エルミート部がゼロ（A_H = 0）となるガラーキンモデルが成功裏に生成され、固定点のリャプノフ安定性と吸引子近傍での有界な軌道が保証された。
• システムの内在的エネルギー保存性のおかげで、人工的散逸を用いなくても、ナビエ＝ストークス吸引子のダイナミクスが高精度で再現された。
• 2次元混合層およびリードド・キャビティという2つの代表的流れに対して、本手法の異なる流れ領域への適用可能性が示された。
• モデルが初期条件のエネルギーを引き継ぎ、A_H = 0 であることがエネルギー保存を直接的に示している。
• 吸引子上では非作用となる経験的散逸項を追加可能であり、コアモデルの安定性を損なわず長期シミュレーションを正確に実行可能である。
• 最適化フレームワークはモード導出アルゴリズムに依存しないため、POD、DMD、あるいはその他の低次元モデル化手法へも適用可能である。