[論文レビュー] Variational Multiscale Proper Orthogonal Decomposition: Navier-Stokes Equations

Traian Iliescu, Zhu Wang|arXiv (Cornell University)|Oct 28, 2012

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

ひとこと要約

本稿では、乱流非圧縮性ナビエ＝ストークス方程式のための変分マルチスケール固有直交分解（VMS-POD-ROM）低次元モデルを提案する。この手法は、物理的精度を向上させるために、渦粘性項を小さな解像スケールにのみ統合する。計算効率が高く、Re=1000における3次元流れのシミュレーションで、標準的なPOD-G-ROMおよび混合長ROMと比較して優れた精度を示し、数値的に確認された誤差推定値を有する。

ABSTRACT

We develop a variational multiscale proper orthogonal decomposition reduced-order model for turbulent incompressible Navier-Stokes equations. The error analysis of the full discretization of the model is presented. All error contributions are considered: the spatial discretization error (due to the finite element discretization), the temporal discretization error (due to the backward Euler method), and the proper orthogonal decomposition truncation error. Numerical tests for a three-dimensional turbulent flow past a cylinder at Reynolds number Re=1000 show the improved physical accuracy of the new model over the standard Galerkin and mixing-length proper orthogonal decomposition reduced-order models. The high computational efficiency of the new model is also showcased. Finally, the theoretical error estimates are confirmed by numerical simulations of a two-dimensional Navier-Stokes problem.

研究の動機と目的

• 標準的なPOD-G-ROMが高次モードのダイナミクスを無視するため、乱流流れにおける物理的精度が低いという問題に対処する。
• 捨て去られたPODモードの効果を変分マルチスケールフレームワークを用いて捉える低次元モデルの開発。
• 大きなスケールのダイナミクスを損なわず、小スケールの解像スケールにのみ渦粘性を適用して解を安定化する。
• 完全な有限要素法および後退オイラー離散化における空間的・時間的・POD切断誤差を考慮した包括的な誤差解析を提供する。
• 2次元および3次元ナビエ＝ストークス問題における数値的妥当性の検証を通じて、理論的誤差推定値の確認と計算効率の向上を示す。

提案手法

• 本手法は、フルオーダー解のスナップショットデータセットから主要モードを抽出するために固有直交分解（POD）を用いる。
• 変分マルチスケール（VMS）フレームワークを適用し、射影作用素を介して解を解像スケールと未解像スケールに分解する。
• 小スケールの解像スケール（モード R+1 から r）にのみ渦粘性（EV）項を追加し、大スケールのダイナミクスを汚さずに、小スケール効果をモデル化する。
• 空間的離散化には有限要素法（FE）を、時間的離散化には後退オイラー法を用い、速度と圧力には発散安定FEペア（P^m/P^{m-1}）を用いる。
• VMS-POD-ROMは、低次元POD基底へのガラーキン射影として定式化され、EV項は解像スケール閾値R未満のモードにのみ追加される。
• 誤差解析では、空間的（FE）、時間的（後退オイラー）、POD切断誤差を考慮し、包括的な理論的境界を導出する。

実験結果

リサーチクエスチョン

• RQ1 小スケールの解像スケールにのみ選択的に渦粘性を適用するVMS-POD-ROMは、標準的なPOD-G-ROMと比較して、乱流ナビエ＝ストークスシミュレーションにおける物理的精度を向上させることができるか？
• RQ2 時間刻み幅および解像モード数Rに関するVMS-POD-ROM誤差の収束率は何か？
• RQ3 空間的・時間的・POD切断誤差の各寄与が、完全離散化における総誤差にどのように寄与するか？
• RQ4 数値的実装は、VMS-POD-ROMに対して理論的に導出された誤差推定値を確認できるか？
• RQ5 3次元乱流流れにおいて、VMS-POD-ROMは著しく精度を向上させつつも、高い計算効率を維持できるか？

主な発見

• P_R-VMS-POD-ROMは、L^2ノルムにおいて時間的収束率が\mathcal{O}((\Delta t)^{0.96})を示し、期待される1次収束率に近く、良好な性能を示す。
• 全誤差の二乗 \|\mathbf{u}^M - \mathbf{u}^M_r\|^2 は、\sum_{j=R+1}^d \|\bm{\varphi}_j\|_1^2 \lambda_j の和に関して2次収束を示し、収束率は1.94であり、理論的予測値1.0を上回る。
• Re=1000における円柱の周りの3次元流れにおいて、P_R-VMS-POD-ROMは標準的なPOD-G-ROMおよび混合長ROMと比較して物理的精度が向上している。
• 本モデルは高い計算効率を維持しており、時間刻みを5\times10^{-3}から3.13\times10^{-4}に減少させた際、L^2誤差は2.31\times10^{-2}から1.45\times10^{-3}に減少した。
• 数値結果は理論的誤差推定値を確認しており、主な誤差寄与はEV項によって導入されるPOD切断誤差である。
• 誤差推定値は、異なるメッシュサイズおよび時間刻みに対して安定しており、h^{2m}および\Delta t^2は10^{-8}のオーダーであり、漸近的にPOD寄与が誤差の主因となる。