[論文レビュー] A Parallel Integrated Computational-Statistical Platform for Turbulent Transport Phenomena

Ali Akhavan-Safaei, Mohsen Zayernouri|arXiv (Cornell University)|Dec 9, 2020

Fluid Dynamics and Turbulent Flows参考文献 32被引用数 6

ひとこと要約

本論文は、分散メモリスーパーコンピュータ上で、三重周期的立方領域における非圧縮性ナビエ-ストークス方程式および拡散輸送方程式のための擬似スペクトル法を用いた、均一かつ等方的乱流およびスカラー輸送の直接数値シミュレーション（DNS）を実行する、オープンソースで並列処理が可能な高精度な計算統計プラットフォームを提示する。本プラットフォームは、統計的量の自動追跡を可能にし、エネルギー、渦度量およびスカラー分散の生成・散逸率の平衡状態により、完全に発達した乱流状態が確認され、25回分の大型渦の回転時間にわたって結果が検証されている。

ABSTRACT

In this paper, we present an open-source, automated, and multi-faceted computational-statistical platform to obtain synthetic homogeneous isotropic turbulent flow and passive scalar transport. A parallel implementation of the well-known pseudo-spectral method in addition to the comprehensive record of the statistical and small-scale quantities of the turbulent transport are offered for executing on distributed memory CPU-based supercomputers. The user-friendly workflow and easy-to-run design of the developed package is disclosed through an extensive and step-by-step example. The resulting low- and high-order statistical records vividly verify well-established and fully-developed turbulent state as well as the seamless statistical balance of conservation laws. Post-processing tools provided in this platform would let the user to readily construct multiple important transport quantities from the primitive turbulent fields.

研究の動機と目的

• 均一かつ等方的乱流の高精度なDNSを実行するための持続可能でオープンソースかつ使いやすい計算統計プラットフォームの開発。
• 分散メモリアーキテクチャ上で、自動的かつ並列的かつ高次の精度で非圧縮性ナビエ-ストークス方程式および移流拡散方程式をシミュレートすることの実現。
• 乱流場の統計的モニタリングを包括的に行い、低次および高次モーメントを含め、完全に発達した統計的に定常な乱流状態の確認を可能にする。
• データ駆動型乱流モデルおよびサブグリッドスケール閉じ込めの訓練および検証に適した高精度なデータ生成を支援すること。
• 流体力学、CFD、乱流輸送分野における教育的ギャップを埋めるために、研究者および学生が統合的かつ容易に利用可能なツールを提供すること。

提案手法

• プラットフォームは、三重周期的立方領域上での非圧縮性ナビエ-ストークス方程式および移流拡散方程式の空間離散化に、フーリエコロケーションを用いた擬似スペクトル法を実装している。
• 時間積分には、乱流の統計的定常性を維持するために動的アーティファクト強制を用いた4次ルンゲ・クッタ（RK4）スキームが採用されている。
• ソフトウェアはMPIおよびPythonを用い、NumPy、SciPy、MPI4Pyとの依存関係を活用して、CPUベースのスーパーコンピュータ上でスケーラブルな並列実行を実現している。
• 初期条件は、所定のエネルギースペクトルのスペクトル的分解により、発散なしで均一かつ等方的な速度場を生成する。
• 各時間ステップで、速度およびスカラー勾配モーメント、分散生成および散逸率といった統計的量が計算され、時系列形式で記録されている。
• 後処理ツールにより、派生的な輸送量の計算が可能であり、統計的平衡状態への収束を統計的均衡を用いて検証できる。

実験結果

リサーチクエスチョン

• RQ1 統計的モニタリングを統合した、スケーラブルでオープンソースのプラットフォームをどのように設計すれば、均一かつ等方的乱流の高精度なDNSを自動化できるか？
• RQ2 DNSにおいて完全に発達した統計的に定常な乱流状態の到来を信頼性高く確認する統計的指標は何か？
• RQ3 完全に発達した乱流速度場において、スカラー勾配統計を適切に解像する条件下で、スカラー輸送をどのように高精度に解像できるか？
• RQ4 プラットフォームの統計的モニタリングは、保存則のバランスおよび確立された乱流スケーリング則の妥当性をどの程度確認できるか？
• RQ5 このプラットフォームは、データ駆動型乱流モデルおよびサブグリッドスケール閉じ込めの訓練および検証に適した高精度なデータを生成できるか？

主な発見

• エネルギーおよび渦度量スペクトルの平衡状態により、プラットフォームは約25回分の大型渦の回転時間経過後に統計的に定常かつ完全に発達した乱流状態に到達したことが確認された。
• スカラー分散の生成率と散逸率の比は、約2回分の大型渦の回転時間後に約1.0に安定化し、スカラー輸送における平衡状態の到来が示された。
• スカラー勾配の高次統計的モーメント、例えばフラットネス係数（Kφ,2）および歪度（Sφ,2）は、t/Te ≥ 5 の時点でそれぞれ約20.8および約1.4の定常値に達した。
• スケール数Sc ≥ 1の範囲で空間解像度が十分であり、バッチォラー波長ηBが速度場の解像度により正しく解像されており、スカラーの動的挙動が良好に解像されていることが保証された。
• プラットフォームの統計的記録は、保存則の滑らかなバランスを確認し、複数の独立した指標を用いて乱流状態の妥当性を検証した。
• オープンソースでモジュール構成された設計により、再現可能でスケーラブルかつ拡張可能なシミュレーションが可能となり、高度な乱流研究およびモデル開発に適している。