Observations of stellar rotation show that low-mass stars lose angular momentum during the main sequence. We simulate the winds of Sun-like stars with a range of rotation rates, covering the fast and slow magneto-rotator regimes, including the transition between the two. We generalize an Alfv\'en-wave driven solar wind model that builds on previous works by including the magneto-centrifugal force explicitly. In this model, the surface-averaged open magnetic flux is assumed to scale as $B_\ast f^{\rm open}_\ast \propto {\rm Ro}^{-1.2}$, where $f^{\rm open}_\ast$ and ${\rm Ro}$ are the surface open-flux filling factor and Rossby number, respectively. We find that, 1. the angular momentum loss rate (torque) of the wind is described as $\tau_w \approx 2.59 \times 10^{30} {\rm \ erg} \ \left( \Omega_\ast / \Omega_\odot \right)^{2.82}$, yielding a spin-down law $\Omega_\ast \propto t^{-0.55}$. 2. the mass-loss rate saturates at $\dot{M}_w \sim 3.4 \times 10^{-14} M_\odot {\rm \ yr^{-1}}$, due to the strong reflection and dissipation of Alfv\'en waves in the chromosphere. This indicates that the chromosphere has a strong impact in connecting the stellar surface and stellar wind. Meanwhile, the wind ram pressure scales as $P_w \propto \Omega_\ast^{0.57}$, which is able to explain the lower-envelope of the observed stellar winds by Wood et al. 3. the location of the Alfv\'en radius is shown to scale in a way that is consistent with 1D analytic theory. Additionally, the precise scaling of the Alfv\'en radius matches previous works which used thermally-driven winds. Our results suggest that the Alfv\'en-wave driven magnetic rotator wind plays a dominant role in the stellar spin-down during the main-sequence.

Using a three-dimensional compressible magnetohydrodynamic (MHD) simulation, we have reproduced the fast solar wind in a direct and self-consistent manner, based on the wave/turbulence driven scenario. As a natural consequence of Alfvenic perturbations at its base, highly compressional and turbulent fluctuations are generated, leading to heating and acceleration of the solar wind. The analysis of power spectra and structure functions reveals that the turbulence is characterized by its imbalanced (in the sense of outward Alfvenic fluctuations) and anisotropic nature. The density fluctuation originates from the parametric decay instability of outwardly propagating Alfven waves and plays a significant role in the Alfven wave reflection that triggers turbulence. Our conclusion is that the fast solar wind is heated and accelerated by compressible MHD turbulence driven by parametric decay instability and resultant Alfven wave reflection.

We conduct a three-dimensional magnetohydrodynamic (MHD) simulation of the parametric decay instability of Alfv\'en waves and resultant compressible MHD turbulence, which is likely to develop in the solar wind acceleration region. Because of the presence of the mean magnetic field, the nonlinear stage is characterized by filament-like structuring and anisotropic cascading. By calculating the timescales of phase mixing and the evolution of Alfv\'en wave turbulence, we have found that the early nonlinear stage is dominated by phase mixing, while the later phase is dominated by imbalanced Alfv\'en wave turbulence. Our results indicate that the regions in the solar atmosphere with large density fluctuation, such as the coronal bottom and wind acceleration region, are heated by phase-mixed Alfv\'en waves, while the other regions are heated by Alfv\'en wave turbulence.

Using numerical simulations, we investigate the onset and suppression of parametric decay instability (PDI) in the solar wind, focusing on the suppression effect by the wind acceleration and expansion. We also discuss the role of PDI in generating the density fluctuation in the solar wind and the cross-helicity evolution up to $1 {\rm \ au}$. Numerical simulations from the coronal base to 1 au are conducted in a self-consistent manner; we take into account the feedback of wave energy and pressure in the background. Monochromatic waves with various injection frequencies $f_0$ are injected to discuss the suppression of PDI, while broadband waves are applied to compare the numerical results with observation. We find, via the simulations with monochromatic waves, that high-frequency ($f_0 \gtrsim10^{-3} {\rm \ Hz}$) Alfv\'en waves are subject to PDI. Meanwhile, the maximum growth rate of the PDI of low-frequency ($f_0 \lesssim10^{-4} {\rm \ Hz}$) Alfv\'en waves becomes negative due to acceleration and expansion effects. Medium-frequency ($f_0 \approx 10^{-3.5} {\rm \ Hz}$) Alfv\'en waves have a positive growth rate but do not show the signature of PDI up to 1 au because the growth rate is too small. The medium-frequency waves experience neither PDI nor reflection so they propagate through the solar wind most efficiently. This is called the frequency-filtering mechanism of the solar wind. The simulations with broadband waves indicate that the observed trend of the density fluctuation is well explained by the evolution of PDI while the observed cross-helicity evolution is in agreement with low-frequency wave propagation.

We propose a novel one-dimensional model that includes both shock and turbulence heating and qualify how these processes contribute to heating the corona and driving the solar wind. Compressible MHD simulations allow us to automatically consider shock formation and dissipation, while turbulent dissipation is modeled via a one-point closure based on Alfv\'en wave turbulence. Numerical simulations were conducted with different photospheric perpendicular correlation lengths $\lambda_0$, which is a critical parameter of Alfv\'en wave turbulence, and different root-mean-square photospheric transverse-wave amplitudes $\delta v_0$. For the various $\lambda_0$, we obtain a low-temperature chromosphere, high-temperature corona, and supersonic solar wind. Our analysis shows that turbulence heating is always dominant when $\lambda_0 \lesssim 1{\rm \ Mm}$. This result does not mean that we can ignore the compressibility because the analysis indicates that the compressible waves and their associated density fluctuations enhance the Alfv\'en wave reflection and therefore the turbulence heating. The density fluctuation and the cross helicity are strongly affected by $\lambda_0$, while the coronal temperature and mass loss rate depend weakly on $\lambda_0$.

Spicule oscillations involve high-frequency components with a typical period approximately corresponding to $40-50$ s. The typical time scale of the photospheric oscillation is a few minutes, and thus, the origin of this high-frequency component is not trivial. In this study, a one-dimensional numerical simulation is performed to demonstrate that the observed spicule oscillations originate from longitudinal-to-transverse mode conversion that occurs around the equipartition layer in the chromosphere. Calculations are conducted in a self-consistent manner with the exception of additional heating to maintain coronal temperature. The analyses indicate the following features: (1) mode conversion efficiently excites high-frequency transverse waves; (2) the typical period of the high-frequency waves corresponds to the sound-crossing time of the mode conversion region; and (3) simulated root-mean-square velocity of the high-frequency component is consistent with the observed value. These results indicate that the observation of spicule oscillation provides direct evidence of mode conversion in the chromosphere.

Using one-dimensional numerical simulations, we study the elementary process of Alfv\'{e}n wave reflection in a uniform medium, including nonlinear effects. In the linear regime, Alfv\'{e}n wave reflection is triggered only by the inhomogeneity of the medium, whereas in the nonlinear regime, it can occur via nonlinear wave-wave interactions. Such nonlinear reflection (backscattering) is typified by decay instability. In most studies of decay instabilities, the initial condition has been a circularly polarized Alfv\'{e}n wave. In this study we consider a linearly polarized Alfv\'en wave, which drives density fluctuations by its magnetic pressure force. For generality, we also assume a broadband wave with a red-noise spectrum. In the data analysis, we decompose the fluctuations into characteristic variables using local eigenvectors, thus revealing the behaviors of the individual modes. Different from circular-polarization case, we find that the wave steepening produces a new energy channel from the parent Alfv\'en wave to the backscattered one. Such nonlinear reflection explains the observed increasing energy ratio of the sunward to the anti-sunward Alfv\'{e}nic fluctuations in the solar wind with distance against the dynamical alignment effect.