[論文レビュー] High-Resolution Numerical Simulations of a Large-Scale Helium Plume Using Adaptive Mesh Refinement

Nicholas T. Wimer, Marc Day+6 more|arXiv (Cornell University)|Jan 11, 2019

Aerodynamics and Acoustics in Jet Flows参考文献 57被引用数 5

ひとこと要約

本研究では、適応メッシュ細分化（AMR）を用いて1 mの直径を有するヘリウム噴流の高解像度数値シミュレーションを実施し、3.91 mmから1.95 mm未塔の解像度が、噴流中心線近くにレイノルズ・ターレイ・不安定性の発生を可能にするため、噴流の動的挙動を顕著に変化させることを明らかにした。これにより、ヘリウムと空気の混合領域が形成され、潜在的コア構造が破壊され、実験的予測を上回る高頻度のパフリングが生じる。本研究は、浮力流の大渦シミュレーションにおけるサブグリッドスケールモデルの検証に向けた高精度なデータセットを提供する。

ABSTRACT

The physical characteristics and evolution of a large-scale helium plume are examined through a series of numerical simulations with increasing physical resolution using adaptive mesh refinement (AMR). The five simulations each model a 1~m diameter circular helium plume exiting into a (4~m)$^3$ domain, and differ solely with respect to the smallest scales resolved using the AMR, spanning resolutions from 15.6~mm down to 0.976~mm. As the physical resolution becomes finer, the helium-air shear layer and subsequent Kelvin-Helmholtz instability are better resolved, leading to a shift in the observed plume structure and dynamics. In particular, a critical resolution is found between 3.91~mm and 1.95~mm, below which the mean statistics and frequency content of the plume are altered by the development of a Rayleigh-Taylor instability near the centerline in close proximity to the base of the plume. This shift corresponds to a plume "puffing" frequency that is slightly higher than would be predicted using empirical relationships developed for buoyant jets. Ultimately, the high-fidelity simulations performed here are intended as a new validation dataset for the development of subgrid-scale models used in large eddy simulations of real-world buoyancy-driven flows.

研究の動機と目的

• 高精度なシミュレーションを用いて、大規模なヘリウム噴流の動的挙動に及ぼす数値解像度の影響を調査すること。
• 浮力噴流におけるケルビン＝ヘルムホルツ不安定性およびレイノルズ・ターレイ不安定性といった主要な不安定性を正確に捉えるために必要な最小解像度を特定すること。
• 将来の大渦シミュレーション（LES）におけるサブグリッドスケール（SGS）モデルの検証に適した、グリッド収束が確認された高精度なデータセットを提供すること。
• 解像度不足が平均統計量および一時的挙動（特にパフリング周波数とコア構造）に与える影響を検討すること。
• 観察された噴流動的挙動の遷移が、高レイノルズ数の浮力ジェットにおけるスケーリング則の変化に対応しているかどうかを評価すること。

提案手法

• 1 mの直径を有するヘリウム噴流を4 m³の領域で、解像度15.6 mmから0.976 mmの範囲で変化させた5つのシミュレーションを、適応メッシュ細分化（AMR）を用いて実施した。
• 2次精度の空間的および時間的離散化を有する有限体積法を採用し、2次統計量までグリッド収束を確保した。
• サブグリッドスケールモデルを一切使用しない明確に定義された初期および境界条件設定により、異なる解像度間での直接比較が可能となった。
• ヘリウム質量分率、速度場、渦度の時間的変化を追跡し、混合、せん断層の発展、不安定性のダイナミクスを分析した。
• シミュレーション全体を通じてパフリング周波数および噴流構造を分析し、流れの挙動に解像度依存の遷移が生じるかを特定した。
• レイノルズ・ターレイ不安定性の発生が噴流動的挙動を根本的に変化させる、臨界解像度範囲（3.91 mmから1.95 mm）を同定した。

実験結果

リサーチクエスチョン

• RQ1 大規模なヘリウム噴流において、ケルビン＝ヘルムホルツ不安定性およびレイノルズ・ターレイ不安定性を正確に解像するための最小空間解像度は何か？
• RQ2 解像度を向上させることで、潜在的コア領域の形成および噴流基部付近のヘリウムと空気の混合領域の発達にどのような影響が生じるか？
• RQ3 実験的スケーリング則と比較して、解像度がパフリング周波数に与える影響はどの程度か？
• RQ4 平均速度や乱流運動エネルギーといった2次統計量が、解像度を向上させることでどのように収束するか？
• RQ5 パフリング周波数の二次的スケーリング則への遷移が、レイノルズ・ターレイ不安定性を解像できる解像度で発生するか？

主な発見

• 3.91 mmから1.95 mmの間の臨界解像度範囲が特定され、この範囲を下回ると、中心線近くにレイノルズ・ターレイ不安定性が発生し、噴流の動的挙動が顕著に変化することが判明した。
• 臨界閾値未満の解像度では、持続的な潜在的コア領域が形成され、Cetegen1996で報告されたRi < 100のための実験的予測と一致するパフリング周波数が観測された。
• より高い解像度では、レイノルズ・ターレイ不安定性によりヘリウムと空気の混合が著しく進行し、混合領域が形成され、主な実験的スケーリング則よりも高いパフリング周波数が生じた。
• 高解像度での観察されたパフリング周波数は、Cetegen1998aの二次的スケーリングよりもわずかに低かったが、異なる力学的状態への遷移を示した。
• 平均速度や乱流運動エネルギーを含む2次統計量は、0.976 mmの最高解像度でグリッド収束が確認され、シミュレーションの信頼性が裏付けられた。
• 初期および境界条件が明確に定義された高精度なデータセットは、実世界の浮力駆動流れにおける大渦シミュレーションのサブグリッドスケールモデルの開発および検証に適していると検証された。