Evolution of MHD turbulence in the expanding solar wind: residual energy and intermittency

We conduct 3D MHD simulations of decaying turbulence in the solar wind context. To account for the spherical expansion of the solar wind, we implement the expanding box model. The initial turbulence comprises uncorrelated counter-propagating 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. The spherical expansion facilitates this process. The resulting residual energy is primarily distributed in the perpendicular direction, with Er∝l⊥ or equivalently Er∝k⊥-2. The magnetic field develops a second-order structure function S2(b)∝l⊥1/2 (Eb∝k⊥-3/2) in most runs. In contrast, the second-order structure function of the velocity S2(v) is consistently shallower than S2(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, has strong dependence on the initial σc and the expansion effect. Generally, the intermittency is more pronounced when the expansion effect is present and when the initial fluctuation is imbalanced (σc ≠0). Finally, we find that in our simulations, although the negative residual energy and intermittency grow simultaneously as the turbulence evolves, there is no obvious causal relation between them as they are generated on different scales.