[論文レビュー] Controlling the level of the ideal invariant fluxes for MHD turbulence using TURBO spectral solver

Bogdan Teaca, Cristian C. Lalescu+2 more|arXiv (Cornell University)|Aug 12, 2011

Solar and Space Plasma Dynamics参考文献 15被引用数 3

ひとこと要約

本論文は、TURBOスペクトルソルバーにおいて、磁気流体乱流における3つの理想的保存量である全エネルギー、クロスヘリシティ、磁気ヘリシティの注入率を独立に制御できる新しいフォースティング機構を提示する。フーリエ空間におけるヘリカル分解を用いることで、統計的定常状態におけるフラックスレベルを正確に制御でき、$5/3$スケーリングが検証された、非対称かつヘリカルなMHD状態の系統的分析が可能になる。

ABSTRACT

The ideal invariants present in the formalism of magnetohydrodynamics (MHD), i.e. global quantities that are conserved in the absence of sources and dissipative effects, play an important role in various theoretical and numerical studies of MHD turbulence. The fluxes of these ideal invariants represent separate channels that transfer the information across different scales in a turbulent system. Once a statistically stationary state of turbulence is reached, the amount of any ideal invariant quantity introduced in the system by a forcing mechanism equals the amount of the same quantity removed by the dissipative effects from the system. For highly developed turbulence, these two mechanisms act predominantly at different scales that are largely separated. Since the ideal invariant quantities cascade between scales, a constant flux is generated with great implication on the state of the system. Numerically, controlling the ideal invariant fluxes levels for a turbulent MHD system is important for the analysis of fundamental MHD turbulence properties. We propose a forcing mechanism that controls the three ideal invariants of MHD turbulence: the total energy, the cross-helicity and the magnetic helicity. This forcing is implemented in the freely available TURBO solver, that is also briefly presented.

研究の動機と目的

• 磁気流体乱流における3つの理想的保存量である全エネルギー、クロスヘリシティ、磁気ヘリシティの注入率を独立に制御するフォースティング機構の開発。
• 特に非対称およびヘリカルな領域において、これらの保存量のフラックスレベルを制御した数値的解析を可能にする。
• 非圧縮性MHD乱流の高精度直接数値シミュレーション（DNS）を実行可能な、自由に利用可能なTURBOスペクトルソルバーへの実装。
• 制御された注入条件下で疑似エネルギースペクトルに一貫した$5/3$スケーリングが得られることを示すことによる手法の妥当性の検証。
• クロスヘリシティおよび磁気ヘリシティの制御注入を可能にすることで、太陽風乱流およびダイナモ物理学分野の研究を支援する。

提案手法

• フォースティングはエールサッター変数を用いてフーリエ空間で定式化され、速度場と磁場を$\mathbf{z}^{\pm} = \mathbf{u} \pm \mathbf{b}$と分解することで、エネルギーとヘリシティフラックスの別々の制御が可能になる。
• 所望の注入率$\varepsilon^{u}_{\pm}$, $\varepsilon^{b}_{\pm}$およびクロスヘリシティ制御パラメータ$\sigma$から、フォースティング振幅$\alpha^{u}_{\pm}(\mathbf{k})$, $\beta^{u}_{\pm}(\mathbf{k})$, $\alpha^{b}_{\pm}(\mathbf{k})$, $\beta^{b}_{\pm}(\mathbf{k})$が導出される。
• 全エネルギー、クロスヘリシティ、磁気ヘリシティの注入率が$\varepsilon^{u}_{\pm}$, $\varepsilon^{b}_{\pm}$, $\sigma$を介して独立に設定され、コーシー＝シュワルツ不等式による制約により特異性を回避する。
• 不均一性を最小限に抑えて統計的定常性を確保するため、多数の$N_f$個のフーリエモードにフォースティングを適用し、力は$\hat{\mathbf{f}}^{\pm}(\mathbf{k}) = \frac{\varepsilon_{\pm}}{N_f} \frac{\hat{\mathbf{u}}_{\pm}(\mathbf{k})}{E^{u}_{\pm}(\mathbf{k})}$と表現される。
• この手法はTURBOスペクトルコードに実装されており、周期的境界条件とスペクトル空間離散化を用いて非圧縮性MHD方程式を解く。
• 統計的定常状態におけるスペクトルとフラックスの計算により、手法が検証され、非対称条件下でも疑似エネルギースペクトルに$5/3$スケーリングが確認され、フラックスの分離が示された。

実験結果

リサーチクエスチョン

• RQ1 スペクトルシミュレーションにおけるMHD乱流において、全エネルギー、クロスヘリシティ、磁気ヘリシティの3つの理想的保存量の注入率をどのように独立に制御できるか。
• RQ2 非ゼロのクロスヘリシティ注入が、MHD乱流における疑似エネルギーフラックスおよびスペクトルの分離に与える影響は何か。
• RQ3 提案されたフォースティング機構は、非対称乱流条件下でも疑似エネルギースペクトルに$5/3$スケーリングを保持するか。
• RQ4 スペクトルソルバーにこの手法を実装しても、不適切な不均一性や数値的不安定性が生じないか。
• RQ5 磁気ヘリシティの制御注入によって、MHD乱流におけるダイナモ過程の研究がどの程度可能になるか。

主な発見

• 提案されたフォースティング機構は、MHD乱流のスペクトルシミュレーションにおいて、全エネルギー、クロスヘリシティ、磁気ヘリシティの注入率を独立に制御できることを示した。
• 非対称乱流（$\sigma = 0.4$）において、2つの疑似エネルギーフラックス（$E^+$ と $E^-$）が分離し、正および負のヘリシティモードにおけるエネルギー伝達ダイナミクスの相違が示された。
• $E^+$ と $E^-$ のスペクトルはインertial領域で$5/3$スケーリングを示し、コルモゴロフに類似したスケーリングと一致した。
• コーシー＝シュワルツ不等式の違反を回避するための制約$\sigma \varepsilon_{\pm} \leq \sqrt{4\varepsilon^{u}_{\pm}\varepsilon^{b}_{\pm}}$を適用することで、数値的安定性が維持され、特異性が回避された。
• クロスヘリシティの注入により、2つの疑似エネルギーの散逸率が変化し、統計的定常状態でそれらのフラックスレベルに顕著な差が生じた。
• TURBOソルバーへの実装により、保存量フラックスを制御した高精度DNSが可能となり、太陽風およびダイナモ現象の研究を支援する。