Development and application of a multiphase Lagrangian structure function model in anisotropic turbulence

技術要約：異方性乱流における多相ラグランジュ構造関数モデルの開発と応用

問題提起
慣性粒子による乱流へのエネルギー応答は、多相ダイナミクスの理解および粒子の分散に関する低次モデルの開発において極めて重要である。これまでに、均質等方性乱流（HIT）における粒子のエネルギー特性に関する多くの研究が行われてきたが、実在する流れは、特に固体境界付近において、空間的な不均一性と異方性を呈することが多い。乱流境界層（TBL）のような壁境界乱流では、平均せん断と境界条件によって大規模な異方性が生成される。主要な未解決の課題は、この大規模な異方性がどのように粒子相へと転送されるのか、そしてそれが形式的に等方的な小スケール領域まで持続するのかどうかである。さらに、これらの環境における粒子の慣性、重力沈降、および流れの異方性の間の複雑な結合関係を正確にモデル化するためには、さらなる調査が必要である。

手法
著者らは、純沈降条件下における乱流境界層内でのスケールごとの粒子エネルギー特性を調査するために、理論と数値計算を組み合わせたアプローチを採用している。

1. 理論的枠組み:

   • 本研究では、キネティック理論を用い、任意のスケールにおける粒子相の二次の構造関数（SoSF）に関する、厳密ではあるが閉じられていない保存則を導出している。
   • 粒子ダイナミクスは、粒子と流体の密度比が非常に大きい極限における、小さく球状で非相互作用な粒子のマックスレイ・ライリー方程式を用いてモデル化されている。
   • 著者らは、縦方向と垂直方向の二次構造関数の差に基づき、粒子相の異方性の指標 $\eta^{(p)}$ を定義している。
   • 相対速度モーメントと動径分布関数に対して冪乗則の挙動を仮定した上で、導出された連続体方程式に対して漸近解析を行い、エネルギーフラックスと異方性を支配する項のバランスを解明している。

2. 数値シミュレーション:

   • 一方向結合の粒子含有壁境界乱流の直接数値シミュレーション（DNS）を、中程度の摩擦レイノルズ数（$Re_\tau = 315$）で行っている。
   • シミュレーションは、TBLの対数領域およびバッファ層を対象としている。
   • 粒子の慣性を評価するためにストークス数（$St^+$）を、重力ドリフトの影響を評価するために沈降パラメータ（$Sv^+$）を変化させている。
   • オイラー的な視点（流体）とラグランジュ的な視点（粒子）の両方が用いられ、特に局所的な垂直ダイナミクスを分離するために、垂直方向の分離距離がゼロである「スラブ方向」の統計に焦点を当てている。

主な貢献

• 保存則の導出: 任意のスケールにおける粒子相のエネルギー特性に関する厳密な保存方程式を導出し、輸送、一次モーメントへの結合、およびエネルギーフラックスのための項を明示的に特定した。これにより、異方性多相流におけるスケールごとのエネルギー転送を分析するための連続体フレームワークが提供される。
• 異方性転送の定量化: 大規模な流れの異方性がどのように粒子相に「寄与」するか、また粒子の慣性と沈降がこの転送をどのように修正するかを系統的に定量化した。
• 漸近解析: 著者らは、DNSデータを解釈するために漸近的な議論を提供し、観測された異方性指標の非単調な挙動を、大規模な流れの異方性と小スケールの間スケールエネルギー転送との間の競争に関連付けて説明している。

結果

• 慣性と異方性: 小さなストークス数では、粒子相は流体と同様に小スケールでほぼ等方的になる。しかし、慣性が増加すると（$St^+ = 10$）、粒子相は小スケールにおいても顕著な異方性を保持しており、これは粒子慣性が大規模な流れの異方性の記憶を保持していることを示している。
• 沈降の影響: 沈降パラメータ（$Sv^+$）の影響は、流体領域に強く依存する。もともと異方性を持つバッファ層（$z^+ = 30$）では、$Sv^+$ の増加に伴い、垂直成分が増加し縦方向成分が減少するという複雑な構造関数の変化が生じる。小スケールが等方的な対数層（$z^+ = 90$）では、高い沈降率（$Sv^+$）は粒子相を小スケールにおいて等方的な状態へと向かわせる傾向がある。
• スケールの非局所性: 粒子の構造関数は、期待される流体の $r^2$ スケーリングと比較して、散逸領域においてより緩やかなスケーリング（$r^{2/3}$ に接近）を示す。これは、主に慣性と沈降によって変調される、慣性粒子のスケール非局所性を反映している。
• 異方性指標の挙動: 粒子の異方性指標 $\eta^{(p)}$ は、分離スケール $r$ が減少するにつれて非単調な挙動を示す。それは極小値に達した後、$r/z \to 0$ に向かって再び増加する。これは、大規模な流れによって課される異方性と、粒子相内の間スケールエネルギー転送によって生成される異方性との間の競争を示唆している。

意義と主張
本論文は、異方性乱流における慣性粒子のスケールごとのエネルギーフラックスに関する重要な知見を提供すると主張している。粒子の慣性が、形式的に等方的な小スケールにおいても大規模な異方性を保持し得ることを示すことで、著者らは粒子相の局所等方性を仮定するモデルの限界を指摘している。導出された保存則と特定された異方性転送のメカニズムは、慣性沈降粒子に関する連続体モデルの開発、特にこれらの非局所的な異方性効果が二方向結合ダイナミクスにどのように影響を与えるかという点において、直接的な示唆を与えるものである。本研究は、正確なモデリングのためには、粒子の慣性と重力ドリフトを同時に考慮することが必要であることを強調している。なぜなら、それらの効果は単に加算的なものではなく、粒子ダイナミクスの中で暗黙的に結合しているからである。