[論文レビュー] An efficient algorithm for simulating ensembles of parameterized flow problems

Max Gunzburger, Nan Jiang+1 more|arXiv (Cornell University)|May 25, 2017

Fluid Dynamics and Turbulent Flows参考文献 16被引用数 3

ひとこと要約

本稿では、異なる粘性係数、初期条件、体積力を持つ複数のパラメータ化された非圧縮性ナビエ＝ストークス流れをシミュレートするための新規な半陰的アンサンブルアルゴリズムを提案する。共通の係数行列を用いて複数の右辺をもつ1つの線形系を構築することで、時間積分の1次精度と最適な空間収束性を維持しつつ、計算コストを顕著に削減する。時間刻みとパラメータのずれに関する制約のもとで、厳密な安定性と誤差推定が得られる。

ABSTRACT

Many applications of computational fluid dynamics require multiple simulations of a flow under different input conditions. In this paper, a numerical algorithm is developed to efficiently determine a set of such simulations in which the individually independent members of the set are subject to different viscosity coefficients, initial conditions, and/or body forces. The proposed scheme applied to the flow ensemble leads to need to solve a single linear system with multiple right-hand sides, and thus is computationally more efficient than solving for all the simulations separately. We show that the scheme is nonlinearly and long-term stable under certain conditions on the time-step size and a parameter deviation ratio. Rigorous numerical error estimate shows the scheme is of first-order accuracy in time and optimally accurate in space. Several numerical experiments are presented to illustrate the theoretical results.

研究の動機と目的

• 異なる粘性係数、初期条件、体積力を有するパラメータ化された非圧縮性流れ問題のアンサンブルをシミュレートするための効率的な数値アルゴリズムの開発。
• 各アンサンブルメンバーに対して個別に解くのではなく、複数の右辺をもつ1つの線形系を解くことで、計算コストの低減。
• 時間刻みとパラメータのずれ比に関する制約のもとで、非線形的かつ長期的安定性を確立。
• 包括的な誤差推定を通じて、時間方向の1次精度と空間方向の最適精度を証明。
• ベンチマーク流れ問題に対する数値実験を通じて理論的結果の妥当性を検証。

提案手法

• 非線形な対流項を陽的に、粘性項および圧力項を陰的に扱う、半陰的時間離散化をナビエ＝ストークス方程式に適用。
• 異なる粘性係数、初期条件、および外力関数に対しても、すべてのアンサンブルメンバーで同一の係数行列が使用されることを保証する、新規な空間離散化を導入。
• 各時間ステップで得られる線形系は、1つの係数行列とJ個の右辺（各アンサンブルメンバーごとに1つ）を有し、ブロックソルバーを用いて効率的に解ける。
• 離散エネルギー法を用いた安定性解析を行い、誤差伝搬を制御するための修正された離散グローヴァル不等式を組み込む。
• 複数のノルム（$L^2$、$H^1$、$L^\rho$型ノルムなど）における誤差推定を導出。その境界は、粘性係数のずれ、時間刻み、メッシュサイズに依存する。
• 可変粘性係数を扱うために、基準粘性係数 $\overline{\nu}$ を導入し、安定性および誤差境界における $|\nu_j - \overline{\nu}|$ のずれを分析。

実験結果

リサーチクエスチョン

• RQ1 異なる粘性係数、初期条件、体積力を有するナビエ＝ストークス流れのアンサンブルシミュレーションにおいて、複数の右辺をもつ1つの線形系を構築できるか？
• RQ2 時間刻みと粘性係数のずれ比にどのような条件を課すと、アンサンブルスキームの非線形的かつ長期的安定性が保証されるか？
• RQ3 提案されたアンサンブルアルゴリズムの時間的および空間的収束速度は何か？また、それが最適か？
• RQ4 誤差境界は、粘性係数のずれ、時間刻み、メッシュサイズに対してどのようにスケーリングされるか？
• RQ5 独立に各シミュレーションを解くのと比較して、提案手法が著しく計算コストを削減できるか？

主な発見

• 提案アルゴリズムは、時間方向に1次精度を達成し、速度および圧力近似において最適な $O(h^k)$ 収束性を示す。ここで $k$ は有限要素空間の多項式次数を表す。
• 時間刻みと粘性係数のずれ比が、$\overline{\nu}$ および $|\nu_j - \overline{\nu}|$ を含む特定の境界を満たす限り、非線形的かつ長期的に安定である。
• 誤差境界は $O(\Delta t^2 + h^{2k})$ のオーダーに比例し、定数は $\overline{\nu}^{-1}$、$\|\nabla u_{j,t}\|_{L^2}$、および解の高階微分に依存する。
• 計算コストは、各時間ステップで $J$ 個の別個の線形系を解くのではなく、1つの系に $J$ 個の右辺をもつものに削減され、大幅に効率が向上する。
• 数値実験により理論的誤差推定が確認され、さまざまなパrameter変動において本手法のロバスト性と効率性が示された。
• 粘性係数がアンサンブルメンバー間で大きく異なる場合でも、ずれ比の条件を満たしていれば、安定性と精度を維持する。