[论文解读] 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谱求解器中控制MHD湍流中三个理想不变量——总能量、交叉螺旋度和磁螺旋度——注入速率的新型强迫机制。通过在傅里叶空间中对强迫进行螺旋分解，该方法能够在统计稳态湍流中精确控制通量水平，从而实现对非对称及具有螺旋度的MHD状态的系统研究，并在伪能谱中验证了$5/3$标度。

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.

研究动机与目标

• 开发一种可独立控制MHD湍流中三个理想不变量——总能量、交叉螺旋度和磁螺旋度——注入速率的强迫机制。
• 实现对这些不变量通量水平受控的MHD湍流数值研究，尤其针对非对称和螺旋度区域。
• 在免费提供的TURBO谱求解器中实现该强迫方法，以支持不可压缩MHD湍流的高保真直接数值模拟（DNS）。
• 通过在受控注入条件下展示伪能谱中一致的$5/3$标度，验证该方法的有效性。
• 通过实现交叉螺旋度和磁螺旋度的受控注入，支持太阳风湍流与发电机物理的研究。

提出的方法

• 强迫在傅里叶空间中通过Elsässer变量表述，将速度场和磁场分解为$\mathbf{z}^{\pm} = \mathbf{u} \pm \mathbf{b}$，从而实现对能量与螺旋度通量的独立控制。
• 强迫振幅$\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$推导得出。
• 该方法通过参数$\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湍流的谱模拟中，如何独立控制总能量、交叉螺旋度和磁螺旋度的注入速率？
• RQ2 非零交叉螺旋度注入对MHD湍流中伪能通量与谱分离的影响是什么？
• RQ3 在非对称湍流条件下，所提出的强迫机制是否能保持伪能谱中的$5/3$标度？
• RQ4 该方法是否可在谱求解器中实现而不引入虚假各向异性和数值不稳定性？
• RQ5 受控的磁螺旋度注入在多大程度上能够支持对MHD湍流中发电机过程的研究？

主要发现

• 所提出的强迫机制在MHD湍流的谱模拟中成功实现了对总能量、交叉螺旋度和磁螺旋度注入速率的独立控制。
• 在$\sigma = 0.4$的非对称湍流中，两个伪能通量（$E^+$和$E^-$）发生分离，表明正负螺旋度模态具有不同的能量传递动力学。
• 在惯性区，$E^+$和$E^-$的谱呈现$5/3$标度，与Kolmogorov类标度一致，即使在非对称强迫条件下亦成立。
• 通过施加约束$\sigma \varepsilon_{\pm} \leq \sqrt{4\varepsilon^{u}_{\pm}\varepsilon^{b}_{\pm}}$，该方法保持了数值稳定性并避免了奇点，确保不违反柯西-施瓦茨不等式。
• 交叉螺旋度的注入改变了两种伪能的耗散率，导致在统计稳态下其通量水平出现可测量的差异。
• 在TURBO求解器中的实现支持了具有受控不变量通量的MHD湍流高保真DNS，为太阳风与发电机现象的研究提供了支持。