[論文レビュー] Ensemble Timestepping Algorithms for Natural Convection

Joseph A. Fiordilino, Khankan, Sarah|arXiv (Cornell University)|Aug 1, 2017

Meteorological Phenomena and Simulations参考文献 10被引用数 6

ひとこと要約

本稿では、乱流でない自然対流のための2つのアンサンブル時刻ステッピング法を提示する。これらの手法は、共有係数行列を用いて複数の右辺を同時に解くことで、記憶容量と計算コストを削減する。手法は、摂動と関連する時間刻み条件の下で安定性と収束性を保証し、アンサンブル平均による予測信頼性の向上を実証する。数値的検証により、予測可能性限界と個々のアンサンブル成員に比べ優れたエネルギー推定が確認された。

ABSTRACT

This paper presents two algorithms for calculating an ensemble of solutions to laminar natural convection problems. The ensemble average is the most likely temperature distribution and its variance gives an estimate of prediction reliability. Solutions are calculated by solving two coupled linear systems, each involving a shared coefficient matrix, for multiple right-hand sides at each timestep. Storage requirements and computational costs to solve the system are thereby reduced. Stability and convergence of the method are proven under a timestep condition involving fluctuations. A series of numerical tests, including predictability horizons, are provided which confirm the theoretical analyses and illustrate uses of ensemble simulations.

研究の動機と目的

• 不確実な初期データを伴う自然対流のシミュレーションにおいて、アンサンブルサイズとメッシュ解像度の間の計算コストのトレードオフを解消すること。
• 複数のアンサンブル解を同時に解くための効率的な時刻ステッピング法を開発し、1つの共有係数行列を用いること。
• アンサンブル分散による予測信頼性の推定と、アンサンブル平均によるエネルギー推定の改善。
• 解の摂動を含む時間刻み条件の下で、アルゴリズムの安定性と収束性を証明すること。
• 予測可能性限界の分析と収束性の検証を通じて、数値的妥当性を検証すること。

提案手法

• 平均速度と摂動項を遅延させることで、流体問題と熱的問題を分離する、陰・陽解法（IMEX）時間離散化を採用する。
• 各時間刻みで、2つの線形系を解く：1つは速度／圧力用、もう1つは温度用。両者とも、すべてのアンサンブルメンバーに共通の係数行列を有する。
• 流体運動量方程式では、対流項にアンサンブル平均速度を用いる。一方、温度方程式には遅延された速度摂動項を含む。
• 薄膜モデルでは、右辺に追加の$u_1$項が含まれ、安定性に影響を与え、解析には離散的グロワール不等式が必要となる。
• アルゴリズムは弱形式で定式化され、空間的離散化に有限要素法が用いられ、適切なノルムにおける収束率が計算される。
• BV（分岐変動）アルゴリズムにより、予測可能性限界の探索とリャプノフ指数の計算のための摂動初期条件が生成される。

実験結果

リサーチクエスチョン

• RQ1 共有係数行列を用いたアンサンブル時刻ステッピングは、計算コストと記憶容量を削減しながらも、自然対流シミュレーションの精度を維持できるか？
• RQ2 アンサンブル平均は、個々のアンサンブルメンバーに比べて、予測信頼性をどのように向上させるか？
• RQ3 薄膜モデルにおける$u_1$項が解の安定性と予測可能性に与える影響は何か？
• RQ4 提案されたアンサンブルアルゴリズムの収束率と安定性条件は何か？
• RQ5 さまざまなレイノルズ数における自然対流の流れの予測可能性限界は何か？

主な発見

• 数値的実験により、アンサンブル平均による解のエネルギー推定が、個々のアンサンブルメンバーに比べて真の系エネルギーをよりよく再現することが示された。
• アルゴリズムは最適な収束率を達成している：速度、温度、圧力の各変数について、時間方向に2次、空間方向に1次（それぞれのノルムにおいて）。
• Ra = 10^4の場合、平均有効リャプノフ指数がt \approx 0.11で正の値を示し、この時点以降に予測不可能性が生じることを示している。
• すべてのテストしたレイノルズ数において、$\delta$-予測可能性限界は負であり、シミュレーション区間$t^* = 0.5$の間、流れは予測可能であることを示している。
• 数値実験により、アンサンブルメンバーの分散が低く保たれていることが確認され、平均解が真の解を強固に予測していることが示された。
• 1つの行列系に複数の右辺を解くことで、$J$個の別個のシステムを解くのではなく、計算コストと記憶容量が削減された。