[論文レビュー] An Artificial Compressibility Ensemble Timestepping Algorithm for Flow Problems

Joseph A. Fiordilino, Michael McLaughlin|arXiv (Cornell University)|Dec 18, 2017

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

ひとこと要約

本稿は、人工圧縮性とIMEX分割を用いて速度、圧力、温度の解を人工的に分離することで、CFL型の条件のもとで1次収束性と非線形エネルギー安定性を達成する、人工圧縮性アンサンブル（ACE）時間積分法を提案する。この手法により、アンサンブルメンバー間で定数係数行列を再利用することで計算コストと記憶容量を削減し、層流領域におけるアンサンブル流れシミュレーションの効率を著しく向上させる。

ABSTRACT

Ensemble calculations are essential for systems with uncertain data but require substantial increase in computational resources. This increase severely limits ensemble size. To reach beyond current limits, we present a first-order artificial compressibility ensemble algorithm. This algorithm effectively decouples the velocity and pressure solve via artificial compression, thereby reducing computational complexity and execution time. Further reductions in storage and computation time are achieved via a splitting of the convective term. Nonlinear energy stability and first-order convergence of the method are proven under a CFL-type condition involving fluctuations of the velocity. Numerical tests are provided which confirm the theoretical analyses and illustrate the value of ensemble calculations.

研究の動機と目的

• 流体力学における不確実パrameterを伴う大規模アンサンブルシミュレーションの高い計算コストとメモリ要件を軽減すること。
• アンサンブル流れ問題の精度を維持したまま、複雑さと実行時間を低減する時間積分法を開発すること。
• 速度、圧力、温度の解を人工圧縮性によって分離することで、メッシュ解像度を損なわずにより大きなアンサンブルサイズを可能にすること。
• 速度フラクチュエーションに依存するCFL型条件のもとで、非線形エネルギー安定性と1次収束性を証明すること。
• 理論的分析を検証する数値実験を通じて、本手法の有効性を示すこと。

提案手法

• 非線形対流項を分離するためにIMEX（陰・陽）時間離散化を採用し、アンサンブルメンバー間で共通の係数行列を確保する。
• 連続の式に罰則項 $ \epsilon p_t $ を追加することで人工圧縮性を導入し、鞍点問題を勾配-発散安定化付きの対流拡散問題に変換する。
• 速度、圧力、温度の式は完全に分離される：速度と温度は対流拡散系を解き、圧力は速度の発散から代数的に更新される。
• 対流項はアンサンブル平均とフラクチュエーション成分に分解される：$ u \cdot \nabla u = \langle u \rangle \cdot \nabla u + u' \cdot \nabla u $、これによりアンサンブルメンバー間での行列不変性が保証される。
• 速度場のフラクチュエーションに依存するCFL型安定性条件が導出され、主に層流流れに限って適用可能となる。
• 空間離散化が得られた線形系に適用される半陰的1次時間積分フレームワークでアルゴリズムが実装される。

実験結果

リサーチクエスチョン

• RQ1 人工圧縮性をアンサンブル時間積分と組み合わせることで、流れシミュレーションにおける計算複雑性と記憶容量を効果的に低減できるか。
• RQ2 提案されたACEアルゴリズムは、速度フラクチュエーションに依存するCFL条件のもとで1次収束性と非線形エネルギー安定性を達成するか。
• RQ3 速度、圧力、温度の解の分離が、計算コストの低減とより大きなアンサンブルサイズの実現にどの程度寄与するか。
• RQ4 レイノルズ数や初期条件の摂動が異なる場合、本手法の精度と予測可能性はどのように変化するか。
• RQ5 圧力の平均効果的リャプノフ指数が負になるのはなぜか、またその影響が予測可能性にどのように現れるか。

主な発見

• 速度、温度、圧力の1次収束性が数値的に確認され、観測された収束率は理論的予測と一致する。
• 圧力誤差は理論的予測よりも半次の高いレートで収束し、圧力近似の精度が向上していることが示唆される。
• アンサンブル平均は摂動のない解に近く保たれるが、摂動付きの解はレイノルズ数が増加するにつれて著しく分離する。特に速度と圧力で顕著に現れる。
• レイノルズ数が増加するにつれ、分散が増加し、速度と圧力の予測における不確実性が増大するが、温度は比較的信頼性が高い。
• 平均効果的リャプノフ指数は時間とレイノルズ数とともに増加し、速度と温度の予測可能性が低下することを示す。一方、圧力では指数が負になるため、時間経過とともに予測可能性が向上する。
• 予測可能性の $ \delta $-ホライズンはレイノルズ数が増加するにつれて短くなり、乱流領域における予測可能性の低下を確認する。特に $ Ra = 10^4 $ で最小のホライズンが観測された。