Residual energy in weakly compressible turbulence with a mean guide field

技術要約：平均ガイド場を持つ弱圧縮性乱流における残留エネルギー

問題の定義
磁気流体力学（MHD）乱流は、空間スケールにわたるエネルギー分布によって特徴付けられ、その性質は音速マッハ数（$M_S$）およびアルヴェン・マッハ数（$M_A$）、ならびに駆動メカニズムに強く影響される。この分布における重要な指標の一つが、運動エネルギーと磁気エネルギーの不均衡を定量化する残留エネルギー $E_r = E_{kin} - E_{mag}$ である。先行研究では、異なる領域における明確な挙動が確立されている。すなわち、非圧縮で強磁場な乱流は、通常 $E_r < 0$（磁気過剰）となり、これは太陽風の観測結果と一致する。一方で、高度に圧縮された超音速乱流は、しばしば $E_r > 0$（運動エネルギー過剰）を示す。しかし、中間的な領域である弱圧縮性領域（$M_S \sim 0.1$）かつ強い平均ガイド場が存在する場合の残留エネルギーの性質については、特に、異なる駆動メカニズム（運動学的駆動か磁気的駆動か）がカスケードや残留エネルギーのスペクトル・スケーリングにどのように影響するかという点において、依然として理解が進んでいない。

手法
著者らは、等温・高粘性・高抵抗MHD方程式を解くために、PENCILコードを用いた一連の直接数値シミュレーション（DNS）を実施した。シミュレーションは、周期境界条件を持つ立方体ドメイン内で、$256^3$ の解像度（収束テストには $512^3$ を使用）で行われた。

• 領域: 弱圧縮性の影響を孤立させるため、音速マッハ数 $M_S \approx 0.1$ の準静的な領域を維持した。
• パラメータ: 3つの初期磁場強度（$B_0 = 0.5, 1.0, 2.0$）をテストした。これらは、プラズマベータ値 $\beta = 4.0, 1.0, 0.3$ に対応している。アルヴェン・マッハ数（$M_A$）もこれに従って変化させた。
• 駆動: 2つの異なる駆動メカニズムを採用した。
  1. 運動学的駆動 ($K$): 大スケール（$k_f \sim 1.5$）において速度変動を注入する。
  2. 磁気的駆動 ($M$): 大スケールにおいて磁場変動を注入する。
• 解析: 体積平均されたエネルギー論、運動・磁気・密度・残留エネルギーのパワースペクトル、およびソレノイダルモードと圧縮モードの相関長を解析した。また、正規化されたクロス・ヘリシティ（$\sigma_c$）および残留エネルギー（$\sigma_r$）の分布についても検討した。

主な結果
シミュレーションの結果、グローバルなマッハ数が同程度であっても、エネルギー分配とスペクトル・スケーリングを決定するのは駆動メカニズムが支配的な要因であることが明らかになった。

1. 磁気駆動乱流:

   • エネルギーバランス: これらのシミュレーションは、慣性領域全体を通じて $E_r \approx 0$ となる均衡したカスケードを示す。アルヴェン比（$R_A = E_{kin}/E_{mag}$）は約 $0.7-0.8$（$E_r < 0$）である。
   • スケーリング: 運動および磁気エネルギーの両方のパワースペクトルは、イロシュニコフ・クライチャニ（IK）スケーリング $E(k) \propto k^{-3/2}$ に従う。
   • ダイナミクス: 乱流はアルヴェン的な揺らぎによって支配されている。グローバルなカスケードは均衡しているが、局所的なパッチでは強いクロス・ヘリシティの不均衡が見られ、これはダイナミック・アライメントの理論と一致する。圧縮モードはアルヴェン・カスケードに受動的に従い、スローモードの圧力バランス関係を満たしている。

2. 運動学的駆動乱流:

   • エネルギーバランス: これらのシミュレーションは、慣性領域のすべてのスケールにおいて、磁気エネルギーに対する運動エネルギーの系統的な過剰（$E_r > 0$）を示す。$R_A$ は $2.5 - 3.4$ である。
   • スケーリング: 運動および磁気エネルギーの両方のパワースペクトルは著しく浅くなり、$E(k) \propto k^{-1}$ とスケーリングする。
   • 残留エネルギーのスケーリング: 残留エネルギーのカスケード（$E_r \propto k^\alpha$）は $\beta$ に依存する：
     • $\beta = 4.0$ のとき：$-2 \lesssim \alpha \lesssim -5/3$
     • $\beta = 1.0$ のとき：$-5/3 \lesssim \alpha \lesssim -3/2$
     • $\beta = 0.3$ のとき：$\alpha \approx -1$
   • ダイナミクス: $k^{-1}$ のスケーリングは、高いコヒーレンスを保持する弱く相互作用するモードを示唆している。正規化されたクロス・ヘリシティは弱く断片化されており、強いダイナミック・アライメントの欠如を示している。密度揺らぎは直接注入され、移流されることで、エルサツァー・フィールド間の結合を維持する緩やかに変化する不均質な背景を作り出している。

3. 揺らぎのスケーリング:
   運動学的駆動における $\delta u \sim \sqrt{B_0}\delta B$ という非圧縮性の研究による示唆に反し、本研究では、弱圧縮性領域において両方の駆動タイプで線形なアルヴェン的スケーリング（$\delta u \sim \delta B$）が見出された。

意義および結論
本論文は、弱圧縮性乱流におけるエネルギー分配は、単なるグローバルなマッハ数ではなく、駆動メカニズムによって強く規定されると結論付けている。

• メカニズムへの依存性: 磁気駆動は、ダイナミック・アライメント理論と一致する、残留エネルギーが無視できるほど小さい均衡したアルヴェン的カスケード（$k^{-3/2}$）をもたらす。一方、運動学的駆動は、密度不均質性の持続によってエルサツァー・フィールドが結合されることで、正の残留エネルギーを伴う浅いカスケード（$k^{-1}$）を生成する。
• 太陽風への文脈: これらの結果は、特定の太陽風領域（例：低速太陽風ストリームやシース領域）で見られる正の残留エネルギーと $k^{-1}$ スケーリングが、運動学的駆動に近い駆動、あるいは密度不均質性に起因している可能性を示唆している。一方で、標準的な $k^{-3/2}$ スケーリングと磁気過剰（$E_r < 0$）は、磁気駆動のアルヴェン的乱流の特徴である。
• 理論的含意: 本研究は、非圧縮性と高度に圧縮された領域の間の溝を埋めるものである。たとえ弱圧縮性（$M_S \sim 0.1$）であっても、運動学的に駆動される場合には密度不均質性が乱流カスケードを根本的に変え、純粋な磁気駆動シナリオではデコヒーレンス（脱コヒーレンス）が起こるはずの相関を維持させることを実証している。