[論文レビュー] A second-order time-stepping scheme for simulating ensembles of parameterized flow problems

Max Gunzburger, Nan Jiang+1 more|arXiv (Cornell University)|Jun 9, 2017

Numerical methods for differential equations参考文献 7被引用数 3

ひとこと要約

本稿では、異なるパラメータ、初期条件、境界条件、および外力関数を有する複数の非圧縮性ナビエ＝ストークス流れのシミュレーションを効率的に行うために、2次精度で、アンサンブルに基づく時間積分スキームを提案する。半implicitなアンサンブル平均化定式化を用いることで、各タイムステップで複数の右辺を有する1つの線形系を解くことができ、計算コストを顕著に低減する一方で、2次精度を維持し、安定性および誤差推定を証明する。

ABSTRACT

We consider settings for which one needs to perform multiple flow simulations based on the Navier-Stokes equations, each having different values for the physical parameters and/or different initial condition data, boundary conditions data, and/or forcing functions. For such settings, we propose a second-order time accurate ensemble-based method that to simulate the whole set of solutions, requires, at each time step, the solution of only a single linear system with multiple right-hand-side vectors. Rigorous analyses are given proving the conditional stability and error estimates for the proposed algorithm. Numerical experiments are provided that illustrate the analyses.

研究の動機と目的

• 物理的パラメータ、初期条件、境界条件、および外力関数が異なる複数の流れシミュレーションの高コスト問題に対処すること。
• 多様な流れアンサンブルにわたって効率を維持する、時間的に正確で2次精度の数値法を開発すること。
• 提案されたアンサンブルスキームの安定性を保証し、厳密な誤差推定を導出すること。
• データ同化、低次元モデル化、感度分析などの応用分野における大規模なシミュレーションを可能にすること。

提案手法

• 速度と圧力の更新に、3段階の後退差分公式を用いた2次精度の半implicit時間離散化を採用する。
• アンサンブル平均速度 $ \overline{u}^n $ と共通の粘性係数 $ \overline{\nu} $ を導入することで、すべてのアンサンブルメンバーに同一の係数行列を適用可能にする。
• 非線形対流項をアンサンブル平均項とフラクチュエーション項の組み合わせで分解することで、一貫性と安定性を確保する。
• 離散系は、アンサンブルメンバー数 $ J $ の右辺を有する1つの線形系として解かれ、LU法やブロック反復法による効率的解法が可能になる。
• 各タイムステップで非圧縮性制約 $ \nabla \cdot u_j^{n+1} = 0 $ が保持される。
• エネルギー推定と離散グローヴァルの不等式を用いた安定性および誤差解析により、$ L^2 $ および $ H^1 $ ノルムにおける最適収束率が得られる。

実験結果

リサーチクエスチョン

• RQ1 異なるパラメータと入力を持つ複数のナビエ＝ストークスシミュレーションにわたって、高い効率を維持する2次精度の時間的に正確なアンサンブルスキームを構築可能か？
• RQ2 複数の右辺を有する1つの線形系を用いて、精度や安定性を損なわずに複数の流れアンサンブルをシミュレートする方法は何か？
• RQ3 2次収束を保証するための、アンサンブル平均と個々のメンバー補正の最適なバランスは何か？
• RQ4 粘性係数や外力関数が変化する条件下で、提案スキームの理論的誤差境界と安定性条件は何か？
• RQ5 各アンサンブルメンバーに対して個別にシステムを解く場合と比較して、本手法の計算コストはどのようにスケーリングされるか？

主な発見

• 適切な正則性仮定の下で、速度と圧力の近似において、時間方向に2次精度が達成される。
• アンサンブル粘性係数と解の正則性に依存する時間ステップ制限の下で、条件付き安定性が保証される。
• 誤差推定では最適収束率が得られる：速度に関して $ H^1 $ ノルムで $ \mathcal{O}(h^{k+1}) $、圧力に関して $ L^2 $ ノルムで $ \mathcal{O}(h^{s+1}) $、ここで $ k $ および $ s $ は有限要素多項式の次数に関連する。
• 各タイムステップで $ J $ 個の右辺を有する1つの線形系を解くことで、$ J $ 個の別個のシステムを解く場合と比較して、計算コストが低減される。
• 誤差境界には $ |\nu_j - \overline{\nu}| $ に依存する項が含まれ、粘性係数の変動に敏感であるが、これは $ \overline{\nu}^{-1} $ で有界であり、相対差が小さい場合には安定性が保証される。
• 数値実験により理論的誤差推定が確認され、特に $ J $ が大きい場合に著しい効率向上が実証されている。