Search for: All records

Abstract 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σ_cand compressibility. The spherical expansion facilitates this process. The resulting residual energy is primarily distributed in the perpendicular direction, withS₂(b) − S₂(u) ∝ l_⊥or equivalently $-E_r\propto k_{\perp }^{-2}$. HereS₂(b) andS₂(u) are second-order structure functions of magnetic field and velocity respectively. In most runs,S₂(b) develops a scaling relation $S_2\left(\mathit{b}\right)\propto l_{\perp }^{1/2}$( $E_b\propto k_{\perp }^{-3/2}$). In contrast,S₂(u) is consistently shallower thanS₂(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σ_cand 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. 

Abstract The heliosphere is permeated with highly structured solar wind originating from the Sun. One of the primary science objectives of Parker Solar Probe (PSP) is to determine the structures and dynamics of the plasma and magnetic fields at the sources of the solar wind. However, establishing the connection between in situ measurements and structures and dynamics in the solar atmosphere is challenging: most of the magnetic footpoint mapping techniques have significant uncertainties in the source localization of a plasma parcel observed in situ, and the PSP plasma measurements suffer from a limited field of view. Therefore, it lacks a universal tool to self-contextualize the in situ measurements. Here we develop a novel time series visualization method named Gaussianity Scalogram. Utilizing this method, by analyzing the magnetic magnitude data from both PSP and Ulysses, we successfully identify in situ structures that are possible remnants of solar atmospheric and magnetic structures spanning more than 7 orders of magnitude, from years to seconds, including polar and midlatitude coronal holes, as well as structures compatible with supergranulation, “jetlets” and “picoflares.” Furthermore, computer simulations of Alfvénic turbulence successfully reproduce the Gaussianization of magnetic magnitude, supporting the observed distribution. Building upon these discoveries, the Gaussianity Scalogram can help future studies to reveal the fractal-like fine structures in the solar wind time series from both PSP and a decades-old data archive. 

Abstract Parker Solar Probe observations reveal that the near-Sun space is almost filled with magnetic switchbacks (“switchbacks” hereinafter), which may be a major contributor to the heating and acceleration of solar wind. Here, for the first time, we develop an analytic model of an axisymmetric switchback with uniform magnetic field strength. In this model, three parameters control the geometry of the switchback: height (length along the background magnetic field), width (thickness along radial direction perpendicular to the background field), and the radial distance from the center of switchback to the central axis, which is a proxy of the size of the switchback along the third dimension. We carry out 3D magnetohydrodynamic simulations to investigate the dynamic evolution of the switchback. Comparing simulations conducted with compressible and incompressible codes, we verify that compressibility, i.e., parametric decay instability, is necessary for destabilizing the switchback. Our simulations also reveal that the geometry of the switchback significantly affects how fast the switchback destabilizes. The most stable switchbacks are 2D-like (planar) structures with large aspect ratios (length to width), consistent with the observations. We show that when plasma beta (β) is smaller than one, the switchback is more stable asβincreases. However, whenβis greater than 1, the switchback becomes very unstable as the pattern of the growing compressive fluctuations changes. Our results may explain some of the observational features of switchbacks, including the large aspect ratios and nearly constant occurrence rates in the inner heliosphere. 

Numerical simulations have been an increasingly important tool in space physics. Here, we introduce an open-source three-dimensional compressible Hall-Magnetohydrodynamic (MHD) simulation codeLAPS(UCLA-Pseudo-Spectral,https://github.com/chenshihelio/LAPS). The code adopts a pseudo-spectral method based on Fourier Transform to evaluate spatial derivatives, and third-order explicit Runge-Kutta method for time advancement. It is parallelized using Message-Passing-Interface (MPI) with a “pencil” parallelization strategy and has very high scalability. The Expanding-Box-Model is implemented to incorporate spherical expansion effects of the solar wind. We carry out test simulations based on four classic (Hall)-MHD processes, namely, 1) incompressible Hall-MHD waves, 2) incompressible tearing mode instability, 3) Orszag-Tang vortex, and 4) parametric decay instability. The test results agree perfectly with theory predictions and results of previous studies. Given all its features,LAPSis a powerful tool for large-scale simulations of solar wind turbulence as well as other MHD and Hall-MHD processes happening in space.