THREE-FLUID, THREE-DIMENSIONAL MAGNETOHYDRODYNAMIC SOLAR WIND MODEL WITH EDDY VISCOSITY AND TURBULENT RESISTIVITY

Abstract

We have developed a three-fluid, three-dimensional magnetohydrodynamic solar wind model that incorporates turbulence transport, eddy viscosity, turbulent resistivity, and turbulent heating. The solar wind plasma is described as a system of co-moving solar wind protons, electrons, and interstellar pickup protons, with separate energy equations for each species. Numerical steady-state solutions of Reynolds-averaged solar wind equations coupled with turbulence transport equations for turbulence energy, cross helicity, and correlation length are obtained by the time relaxation method in the corotating with the Sun frame of reference in the region from 0.3 to 100 AU (but still inside the termination shock). The model equations include the effects of electron heat conduction, Coulomb collisions, photoionization of interstellar hydrogen atoms and their charge exchange with the solar wind protons, turbulence energy generation by pickup protons, and turbulent heating of solar wind protons and electrons. The turbulence transport model is based on the Reynolds decomposition and turbulence phenomenologies that describe the conversion of fluctuation energy into heat due to a turbulent cascade. In addition to using separate energy equations for the solar wind protons and electrons, a significant improvement over our previous work is that the turbulence model now uses an eddy viscosity approximation for the Reynolds stress tensor and the mean turbulent electric field. The approximation allows the turbulence model to account for driving of turbulence by large-scale velocity gradients. Using either a dipole approximation for the solar magnetic field or synoptic solar magnetograms from the Wilcox Solar Observatory for assigning boundary conditions at the coronal base, we apply the model to study the global structure of the solar wind and its three-dimensional properties, including embedded turbulence, heating, and acceleration throughout the heliosphere. The model results are compared with plasma and magnetic field observations on WIND, Ulysses, and Voyager 2 spacecraft.