Evolution of Magnetohydrodynamic Turbulence in the Expanding Solar Wind: Residual Energy and Intermittency
We conduct 3D magnetohydrodynamic simulations of decaying turbulence in the context of the solar wind. To account for the spherical expansion of the solar wind, we implement the expanding box model. The initial turbulence comprises uncorrelated counterpropagating Alfvén waves and exhibits an isotropic power spectrum. Our findings reveal the consistent generation of negative residual energy whenever nonlinear interactions are present, independent of the normalized cross helicity σc and compressibility. The spherical expansion facilitates this process. The resulting residual energy is primarily distributed in the perpendicular direction, with S 2 (b) − S 2 (u) ∝ l ⊥ or equivalently –E r ∝ k$^{–2}_{⟂}$. Here S 2 (b) and S 2 (u) are second-order structure functions of magnetic field and velocity respectively. In most runs, S 2 (b) develops a scaling relation S 2 (b) ∝ $l$$^{1/2}_{⊥}$ (E b ∝ k$^{–3/2}_{⟂}$). In contrast, S 2 (u) is consistently shallower than S 2 (b), which aligns with in situ observations of the solar wind. We observe that the higher-order statistics of the turbulence, which act as a proxy for intermittency, depend on the initial σc and are strongly affected by the expansion effect. Generally, the intermittency is more pronounced when the expansion effect is present. Finally, we find that in our simulations, although the negative residual energy and intermittency grow simultaneously as the turbulence evolves, the causal relation between them seems to be weak, possibly because they are generated on different scales.