Circulation Statistics in Rayleigh-Bénard Convection

技術要約：レイリー・ベナール対流における循環統計

問題提起
近年の進展により、均質等方乱流（HIT）における速度循環の重要な統計的性質、特に、循環のゆらぎを小規模な初等渦の系として記述する「渦ガスモデル（VGM）」に関する知見が明らかになっている。中心的な問いは、これらの統計的な現象論が、他のクラスの乱流、具体的には浮力駆動型流においても維持されるのか、あるいは修正を受けるのかという点である。本研究では、レイリー数（$Ra$）が$10^9$、プラントル数（$Pr$）が$O(1)$のレイリー・ベナール対流（RBC）における速度循環の統計を調査する。主な目的は、循環のゆらぎを特徴付け、それらの基礎となる初等渦を特定することで、HITで確立されたVGMによる記述がどの程度RBCに拡張可能かを評価することである。特に、エネルギー散逸によって駆動されるHITの渦分布とは対照的に、熱的なプルームがコヒーレントな渦運動を駆動する熱境界層（TBL）に細心の注意を払う。

手法
本研究では、無次元ブシネスク方程式の直接数値シミュレーション（DNS）データを利用する。シミュレーション領域は、横方向のアスペクト比 $L/H = 4$ のボックスであり、9次レジェンドル多項式を用いた $150 \times 150 \times 96$ 要素のスペクトルメッシュを採用し、約 $1.57 \times 10^9$ 個のコロケーションポイントを持つ。流れの解析は、加熱壁から様々な距離にある2D水平断面平面（$\delta_T/4$、$\delta_T/2$、$\delta_T$、$4\delta_T$、および $H/2$）において行う。ここで $\delta_T$ は、温度変動プロファイルから推定された熱境界層の厚さである（$\delta_T/H \approx 7.47 \times 10^{-3}$）。

主要な観測量は以下の通りである：

1. 初等渦： 2D速度勾配テンソルの旋回強度基準（$\lambda_I$）を用いて、2Dスライシング平面上で特定される。構造は $|\lambda_I| \geq \sigma_\lambda/8$（$\sigma_\lambda$ は虚数固有値の標準偏差）を満たす場合に保持される。
2. 循環 ($\Gamma$)： スライシング平面上で定義された、任意の向きを持つ閉曲線（矩形および正方形）上の線積分によって計算される。
3. 統計分析： 渦スポットと温度場の相関、渦の形態（アスペクト比）、および循環の確率密度関数（PDF）を調査する。これらは、VGMからの解析的予測および（ジョンズ・ホプキンス・タービュランス・データベースから得られた）HITのデータと比較される。

主要な貢献および結果

• TBLにおける渦と温度の相関： 熱境界層内において、初等渦の分布は温度場と強く相関している。「ホット渦」（瞬時温度が局所平均を超える領域に位置する渦）の割合は、壁付近で最も高く（$\delta_T/4$ で67.6%）、距離とともに減少する。しかし、ホット渦の相対的な数密度は $\delta_T$ 付近でピークに達しており、これは熱境界層の端に渦が集中していることを示唆している。これは、渦分布がエネルギー散逸場によって決定されるHITとは対照的である。
• 渦の形態： TBL内の渦スポットは、HITよりも著しく伸長している（アスペクト比が高い）。アスペクト比の確率密度関数（PDF）は指数関数的な裾を持ち、平面が壁に近づくにつれて、その裾が緩やかになる（より伸長することを示す）。TBLを超えると、アスペクト比の分布は徐々にHITで観察される統計量へと回復する。
• 初等渦の支配性： 本研究は、RBCにおける循環ゆらぎへの主要な寄与が、特定された初明渦スポットによって運ばれる循環に由来することを裏付けている。個々の渦の寄与から再構成された標準化循環PDFは、渦度フラックスから計算された正確な循環PDFと密接に一致している。
• 循環PDFと面積則：
  • 初等渦： RBCにおける個々の初等渦の循環PDFは、HITとは異なり、より軽い裾を持つ。これは、異なる間欠性メカニズム（熱プルーム vs エネルギー散逸）を反映している。
  • 平面輪郭： 壁に平行な輪郭については、循環PDFはHITに似た特徴を示す。具体的には、「面積則」——循環統計と最小曲面の間の関係——がよく満たされている。同じ面積を囲む輪郭は同じ標準化PDFを共有し、その分散は慣性領域の外側においてのみ幾何学形状に依存する。
  • 自己相似性： TBL内において、循環のゆらぎは近似的な自己相似性を示す。異なる輪郭サイズに対するPDFは、正規化されると崩壊（重なり）し、主に分散の違いによってのみ区別される。
• スケーリング指数： 循環のスケーリング指数 ($\zeta_p$) は、低次においてはモーメントの次数に対して線形依存性を示すが、HITの値からは逸脱する。サンプリング平面が加熱壁に近づくにつれ、循環のヘルダー指数は段階的に減少する。HITの指数からの顕著な逸脱は、バルク（$H/2$）においても観察され、これはチャネル流と同様である。
• 遷移領域： 壁から $4\delta_T$ の距離において、異常な非単調挙動が観察される。この領域では、ホット渦の割合および渦スポットの直径がバルクの値よりも低下する。興味深いことに、ここでのアスペクト比のPDFはHITに酷似しており、面積則も厳密には満たされていない。これは、TBLとバルクの間の遷移的な流れのレジームを示唆している。

意義および主張
本論文は、渦形成を駆動する具体的なメカニズムは異なる（RBCでは熱プルーム、HITではエネルギー散逸）ものの、循環統計の統計的枠組みとしてのVGMは、RBCに対しても概ね適用可能であると主張している。具体的には、循環のゆらぎが初等渦によって支配されていること、および壁に平行な平面輪郭に対して面積則が成立するという発見は、異なる乱流レジーム間における循環統計の強固な普遍性を示唆している。

しかし、著者らはこれらのゆらぎの理論的モデリングに関して慎重な立場をとっている。彼らは、VGMがHITにおける循環の「運動学」を記述することには成功しているものの、RBCのTBLにおける変動する温度場と、その渦分布との相関に関する類似の記述は、「未解決かつ困難な問題」として残っていると述べている。本研究は、RBCにおける渦の間欠的な空間分布が、エネルギー散逸場ではなく熱プルームによって決定されることを強調しており、温度と渦の結合を完全にモデル化するためには、新たな理論的アプローチが必要であることを示している。また、$4\delta_T$ における、異なる統計的特性を持つ遷移領域の特定は、より大きな統計アンサンブルによるさらなる調査を要する知見として提示されている。