[論文レビュー] Delineation of the flow and mixing induced by Rayleigh-Taylor instability through tracers

Ge Zhang, Aiguo Xu+4 more|arXiv (Cornell University)|Dec 31, 2019

Fluid Dynamics and Turbulent Flows参考文献 87被引用数 22

ひとこと要約

本研究では、可圧縮レイノルズ＝ターナー不安定性（RTI）における流れ構造および混合ダイナミクスを可視化・定量するために、離散ボルツマンモデル（DBM）にトレーサーを導入する。色分けされたトレーサーを追跡することで、微細構造の流れのパターンと熱力学的非平衡（TNE）行動が明らかになり、トレーサーに基づく混合度指標により、大規模な初期構造と小規模な後期乱流に起因する二段階混合進化が解明された。主な貢献は、粘性係数の飽和効果と、混合に及ぼす圧縮性・粘性の二段階的影響の同定であり、RTI駆動混合メカニズムに関する新たな知見を提供する。

ABSTRACT

Rayleigh-Taylor-instability(RTI) induced flow and mixing are of great importance in both nature and engineering scenarios. To capture the underpinning physics, tracers are introduced to make a supplement to discrete Boltzmann simulation of RTI in compressible flows. Via marking two types of tracers with different colors, the tracer distribution provides a clear boundary of two fluids during the RTI evolution. Fine structures of the flow and thermodynamic nonequilibrium behavior around the interface in a miscible two-fluid system are delineated. Distribution of tracers in its velocity phase space makes a charming pattern showing quite dense information on the flow behavior, which opens a new perspective for analyzing and accessing significantly deep insights into the flow system. RTI mixing is further investigated via tracer defined local mixedness. The appearance of Kelvin-Helmholtz instability is quantitatively captured by mixedness averaged align the direction of the pressure gradient. The role of compressibility and viscosity on mixing are investigated separately, both of which show two-stage effect. The underlying mechanism of the two-stage effect is interpreted as the development of large structures at the initial stage and the generation of small structures at the late stage. At the late stage, for a fixed time, a saturation phenomenon of viscosity is found that further increase of viscosity cannot see an evident decline in mixedness. The mixing statues of heavy and light fluids are not synchronous and the mixing of a RTI system is heterogenous. The results are helpful for understanding the mechanism of flow and mixing induced by RTI.

研究の動機と目的

• 可圧縮で可混和性を示す流体系におけるレイノルズ＝ターナー不安定性（RTI）が誘発する詳細な流れおよび混合ダイナミクスを調査すること。
• 流体界面における熱力学的非平衡（TNE）効果を捉えることが難しい従来のナビエ＝ストークスモデルの限界を克服すること。
• 高時間・空間分解能を実現するトレーサーに基づく手法を開発し、RTIにおける微細構造の進化および界面挙動を可視化すること。
• トレーサー分布から導出した新しい混合度指標を用いて、圧縮性および粘性がRTI混合に与える影響を定量すること。
• RTI系で観察された二段階混合挙動の背後にある物理的メカニズムを解明すること。

提案手法

• 流体界面の進化を追跡し、重い流体と軽い流体の領域を区別するために、多重緩和時間離散ボルツマンモデル（DBM）にトレーサーを導入する。
• 異なる色で着色された2種類のトレーサーを用いることで、RTI発展過程における流体間の明確な視覚的境界を提供する。
• トレーサーの空間的分布に基づいて、領域内各点における混合度の局所的指標χpを定義し、混合度の定量的分析を可能にする。
• 混合度を縦方向に平均化し、横方向に分析することで、ケルビン＝ヘルムホルツ不安定性（KHI）の発生時期と強度を検出する。
• 混合度プロファイルのロバストネスを評価し、過剰な平滑化や過剰なフラクチュエーションを回避するために、異なるサイズの空間セル（n = 1.0および8.0）を用いた統計解析を実施する。
• 最適な指数フィッティング（χp = χ0 + C0·exp(t*/t0)）を用いたスケール解析により、粘性および圧縮性の関数として特徴時間尺度t0および特徴混合度C0を抽出する。

実験結果

リサーチクエスチョン

• RQ1 トレーサーは、RTIにおける流れ構造および界面ダイナミクスの可視化と定量的評価をどのように向上させるか？
• RQ2 圧縮性は、RTI混合の二段階的進化において果たす役割は何か？ また、特徴時間尺度t0および混合度C0にどのように影響を与えるか？
• RQ3 粘性はRTI混合にどのように影響を与えるか？ また、観察された二段階的挙動および粘性係数の飽和効果の原因は何か？
• RQ4 トレーサーに基づく混合度指標は、RTI系におけるケルビン＝ヘルムホルツ不安定性（KHI）の発生時期と強度を信頼性高く検出できるか？
• RQ5 RTI混合における大規模構造から小規模構造への優位性の転移を引き起こす物理的メカニズムは何か？

主な発見

• トレーサー手法により、特に速度位相空間における微細構造の流れのパターンと熱力学的非平衡（TNE）行動が明確に可視化され、RTIダイナミクスに関する豊富な知見が得られた。
• 混合度χpは時間とともにほぼ指数関数的に増加し、圧縮性が増加するにつれて特徴時間尺度t0が単調に減少することが示された。
• 粘性は混合度に二段階的効果を示す：初期段階では大規模構造の浸透を促進して混合度を増加させるが、後期段階では小規模構造の抑制により混合度を低下させ、最終的に粘性係数の飽和効果が現れ、さらなる粘性係数の増加が混合度に影響しなくなる。
• 特徴時間尺度t0は初期段階で粘性の増加に伴い増加し、その後わずかに減少する一方、特徴混合度C0も同様の二段階的傾向を示し、初期段階と後期段階で異なる物理的メカニズムが働くことが示された。
• 横方向に平均化した混合度プロファイルから、ケルビン＝ヘルムホルツ不安定性（KHI）の発生が検出され、圧力勾配方向に整列した平均混合度によりその強度が定量的に評価された。
• 重い流体と軽い流体の混合は非同期的であり、混合プロセスは不均一である。システムは流体の圧縮性および粘性に応じて顕著な異なる挙動を示した。