[論文レビュー] Long-Time Reynolds Averaging of Reduced Order Models for Fluid Flows: Preliminary Results

Luigi C. Berselli, Traian Iliescu+2 more|arXiv (Cornell University)|Jan 15, 2019

Model Reduction and Neural Networks参考文献 39被引用数 12

ひとこと要約

本稿は、流体流れの低次元モデル（ROM）におけるモード間エネルギー交換の長時間平均について検討し、ストークス固有関数と固有正交分解（POD）モードを比較する。1次元バーガース方程式を用いて、PODモードにおける勾配の非直交性に起因するエネルギー移動の違いを解析的に証明し、十分に長い時間区間において、低インデックスから高インデックスのPODモードへの時間平均エネルギー交換が正であることを数値的に確認した。これは、統計的平衡状態におけるエネルギーの前向き移動を支持する。

ABSTRACT

We perform a theoretical and numerical investigation of the time-average of energy exchange among modes of Reduced Order Models (ROMs) of fluid flows. We are interested in the statistical equilibrium problem, and especially in the possible forward and backward average transfer of energy among ROM basis functions (modes). We consider two types of ROM modes: eigenfunctions of the Stokes operator and Proper Orthogonal Decomposition (POD) modes. We prove analytical results for both types of ROM modes and we highlight the differences between them. We also investigate numerically whether the time-average energy exchange between POD modes is positive. To this end, we utilize the one-dimensional Burgers equation as a simplified mathematical model, which is commonly used in ROM tests. The main conclusion of our numerical study is that, for long enough time intervals, the time-average energy exchange from low index POD modes to high index POD modes is positive, as predicted by our theoretical results.

研究の動機と目的

• 流体流れの低次元モデル（ROM）におけるモード間エネルギー交換の統計的平衡状態を調査すること。
• ROMにおけるストークス固有関数とPODモードの間でのエネルギー移動行動を比較すること。
• 十分に長い時間区間において、PODモード間の時間平均エネルギー交換が正であるかどうかを特定すること。
• POD基底関数の勾配における非直交性がエネルギー移動ダイナミクスに与える影響を評価すること。
• エールゴディシティ仮説に依存せずに、長時間平均化をROMに適用する理論的・数値的基盤を提供すること。

提案手法

• ROMにレインズオーダー平均フレームワークを適用し、時間平均エネルギー交換に関する理論的分析を実施する。
• ストークス固有関数とPODモードの両方について解析的結果を導出し、勾配の直交性の違いに起因する差異を強調する。
• 分片線形有限要素空間離散化と後退オイラー時間積分を用いた1次元バーガース方程式の数値シミュレーション。
• 時間区間と時間刻みを変化させた上で、合成台形則を用いて時間平均を計算する。
• エネルギー交換項 $\mathrm{e}_m(u)$ および $\mathcal{E}_m(u)$ の評価を行い、前向き／後向き移動を評価する。
• 収束性と数値アーチファクトの回避を確認するため、複数の時間区間、時間刻み、部分区間数を用いた体系的検証を実施する。

実験結果

リサーチクエスチョン

• RQ1 十分に長い時間区間において、PODモード間の時間平均エネルギー交換は、低インデックスから高インデックスモードへのネット前向き移動を示すか？
• RQ2 POD基底関数の勾配における直交性の欠如が、ストークス固有関数と比較してエネルギー移動ダイナミクスに与える影響は何か？
• RQ3 1次元バーガース方程式をテストモデルとして用いることで、ROMにおけるエネルギー交換に関する解析的結果を数値的に検証可能か？
• RQ4 数値的時間平均は、異なる時間区間および離散化パラメータにおいて、どの程度収束し、安定性を保っているか？
• RQ5 長時間平均エネルギー交換は、統計的平衡状態における前向きエネルギー移動の理論的予測と整合的か？

主な発見

• 長時間区間（例：$[0,100]$）において、低インデックスから高インデックスのPODモードへの時間平均エネルギー交換が正であることが確認され、理論的予測を裏付けた。
• 時間平均エネルギー交換項 $\mathrm{e}_m(u)$ の大きさは $\mathcal{E}_m(u)$ より顕著に大きく、非線形項によるエネルギー移動が支配的であることを示唆している。
• 短時間区間（$[0,0.1]$）では、一部のモードについて時間平均エネルギー交換が負になることが判明し、一時的で平衡でない挙動を示した。
• 異なる時間刻みと部分区間数においても数値結果が安定しており、時間平均化手順の収束性と信頼性を確認した。
• 理論的分析により、PODモードは勾配の非直交性のため、ストークス固有関数に比べてエネルギー移動の特徴づけがやや不正確であることが示された。
• エールゴディシティを仮定せず、直接的な時間平均化に依拠することで、ROMにおける長時間平均化の基盤を確立した。