Wave Modeling of the Solar Wind

The acceleration and heating of the solar wind have been studied for decades using satellite observations and models. However, the exact mechanism that leads to solar wind heating and acceleration is poorly understood. In order to improve the understanding of the physical mechanisms that are involved in these processes a combination of modeling and observational analysis is required. Recent models constrained by satellite observations show that wave heating in the low-frequency (MHD), and high-frequency (ion-cyclotron) range may provide the necessary momentum and heat input to coronal plasma and produce the solar wind. This review is focused on the results of several recent solar modeling studies that include waves explicitly in the MHD and the kinetic regime. The current status of the understanding of the solar wind acceleration and heating by waves is reviewed.


Introduction
The solar wind plasma plays a major role in the physical connection between the Sun and the Earth, and affects the plasma conditions in the heliosphere. The solar wind is an important component of space weather, and forms the background upon which solar disturbances, such as CMEs and energetic particles propagate towards the Earth. However, the exact physics of formation of the solar wind is poorly understood. Remote sensing and in-situ observations indicate that there are two types of solar wind: slow and fast. The fast solar wind, reaching the speed of ∼ 800 km s -1 at 1 AU, is steady and of low density (few particles per cm 3 at 1 AU); the slow solar wind reaches half the above speed with an order of magnitude higher density (see Figure 1). At solar minimum the solar wind speed is a clear function of latitude as measured by Ulysses' SWOOPS instrument with the fast wind emerging from polar coronal holes, and the slow wind confined to the equatorial regions. Near the solar maximum the corona is dominated by streamers and by slow solar wind. The heavy ion composition of the solar wind is correlated with the fast and slow wind, indicating their coronal origin. Figure 1: The solar wind speed as a function of latitude (in km s -1 ) measured by Ulysses' SWOOPS instrument near solar minimum (left panel) and near solar maximum (right panel). The direction of the magnetic field is marked (red-outward; blue-inward). The typical composite solar images near minimum (8/17/96) and maximum (12/07/00) are shown using data from SOHO/LASCO, EIT, and Mauna Loa K-coronameter images (McComas et al., 2003).
The twin Helios spacecraft investigated the heliosphere in the region 0.3 to 1 AU and found evidence for magnetic fluctuations spectrum, and investigated the radial dependence of plasma parameters and the velocity distributions in the solar wind. These measurements provide the basis of many theoretical studies of solar wind acceleration and heating (see the review by Marsch, 2006). Recent in-situ observations by ACE (Stone et al., 1998), Wind (e.g., Lepping et al., 1995), Ulysses, (e.g., Balogh et al., 1992), and other spacecrafts confirm that the solar wind contains magnetic fluctuations that have power law dependence with frequency (see Figure 2 adapted from Smith et al., 2006). The magnetic fluctuations exhibit high correlation with velocity fluctuations indicating their Alfvénic nature. Remote sensing observations show that the solar wind is heated and accelerated close to the Sun within 10 ⊙ (SOHO/UVCS; Kohl et al., 1997). The observed kinetic and compositional properties provide clues on coronal origin (i.e., coronal holes, active regions) (e.g., Ko et al., 2008), and on the acceleration and heating mechanism. However, the interpretation of observations requires theoretical and computational modeling of the global and kinetic properties of the solar wind to understand the physics and the dynamics of the multi-ion solar wind plasma acceleration and heating, and improve the accuracy of space weather forecasting. Year/day of sample (decimal day of year in 1998) and Power Spectral Density (PSD) are shown. The top curve represents the summed power in the two components perpendicular to the mean magnetic field, and the lower curve represents the power in the parallel component. The fit functions with power exponent -1.56 ± 0.3 (upper) and -1.8 ± 0.03 (lower) are shown (see Smith et al., 2006, for further details).
In the region between 1 ⊙ and 1 AU the corona and the heliosphere span many orders of magnitudes in relevant temporal-spatial scales, density, magnetic field strength, and in other physical parameters. As a result, the physical processes in the plasma and the modeling approaches are dramatically different in the collisional lower corona with densities on the order of 10 8-9 cm -3 , magnetic field strength of 1 to 100 G, and the heliospheric plasma near 1 AU with densities of few particles per cm -3 and magnetic field strength of few nano-tesla (few tens of micro-Gauss). Correspondingly, the plasma frequency, collision frequency, and the proton gyrofrequency vary over many orders of magnitude between 1 ⊙ and 1 AU. This disparity of scales and physical regimes leads to difficulty in the physical modeling of the heliospheric plasma with a single 'all-inclusive' model. For example, solving the Boltzmann equations in three spatial dimensions and six degrees of freedom is not practical with current or foreseeable future computational resources. A practical modeling approach that is applicable on large scales (usually, MHD) can not resolve the physics of the small scale (kinetic) phenomena. The global models provide little or no information on the kinetics of the heating and acceleration processes. The solution of this difficulty is a multi-level modeling approach, where global models are used for large scale structures, with small scales pa-rameterized as external (to MHD) diffusion coefficients and heating/loss functions. Kinetic models are used to study the physics of the local heating and dissipation phenomena on small spatialtemporal scales. The output of the kinetic models can serve as a guide of the diffusion processes that need to be included in global models.
As a result, there are two main types of solar wind models: (1) observationally driven global MHD models that provide the overall global structure of the solar wind as shaped by the interaction with the solar magnetic field, rotation, gravity, and heating; (2) local models that deal with the physical processes of heating and acceleration of the solar wind magnetized plasma with the initial state and boundary conditions not necessarily tied to a particular observation. The output of these models provides the solar wind speed as a function of position, density, temperature, and other parameters. Usually, the acceleration of the solar wind is modeled by introducing simplifying assumption into the energy equation, such as taking the plasma to be isothermal at coronal temperature, or assuming polytropic index below the adiabatic = 5/3, that can also vary with heliocentric distance (Cohen et al., 2007). The first isothermal and polytropic solar wind expansion models were developed by Parker (1958Parker ( , 1963 in one spatial dimension. Present global models provide 3D structure of the solar wind and try to match empirically the observations at 1 AU by adjusting the model parameters (e.g., Mikić et al., 1999;Linker et al., 1999;Usmanov et al., 2000;Roussev et al., 2003;Usmanov and Goldstein, 2003;Cohen et al., 2007). Recently, 2D MHD models that use observations to constrain the solar wind heating function and momentum input were developed (Sittler Jr and Guhathakurta, 1999;Vásquez et al., 2003;Guhathakurta et al., 2006;Sittler Jr and Ofman, 2006). The empirical heating function was recently included in 2.5 solar wind MHD models Airapetian et al., 2011).
On the other end of the scale are the kinetic models which provide the description of the small scale interaction in the solar wind plasma between the waves, the ions, and the background magnetic field. These models usually do not provide the information on the global structure of the solar wind in the heliosphere. However, the kinetic models are well suited for the investigation of the kinetic processes and instabilities involved in heating of solar wind-like plasma. Due to the complexity of the kinetic models, and the necessity to resolve the fine scale on the proton, or even electron Larmor radius scale, it is computationally difficult to include global scale structures.
The range in between the above two modeling approaches is occupied by MHD models that include explicit heating and acceleration of the solar wind by MHD waves. Following Osterbrock (1961) study suggesting MHD waves for the heating of the solar chromosphere and corona, the acceleration of the solar wind by Alfvén waves was studied in the past (Barnes, 1969;Alazraki and Couturier, 1971;Belcher, 1971;Belcher and Davis Jr, 1971;Heinemann and Olbert, 1980). Recently, one-dimensional (Cranmer and van Ballegooijen, 2005;Inutsuka, 2005, 2006;Cranmer et al., 2007) and three-dimensional MHD models (e.g., Usmanov et al., 2000;Evans et al., 2009) were developed (see the review by Ofman, 2005). However, due to the requirements of time step and resolution, the waves are not included explicitly in global models. They are modeled by an additional wave-pressure term, and wave-energy equations. Only few models consider the resolved waves explicitly in 2.5D MHD models Davila, 1997, 1998). The next level of plasma approximation of fully resolved wave driven wind is via multi-fluid models, that describe each particle species as separate fluid (Ofman and Davila, 2001;Ofman, 2004a). The fluids interact through momentum and energy exchanges, and through electromagnetic interactions resulting from quasi-neutrality condition, and through their contribution to the net current. These models can be tested directly by comparing to observations that contain heavy ion emission (e.g., Ofman, 2004b;Abbo et al., 2010).
Observations with SOHO Ultraviolet Coronagraph Spectrometer (UVCS) show that heavy ions such as O 5+ , and Mg 9+ undergo preferential perpendicular heating, causing large temperature anisotropy ( ⊥ / ‖ > 10), are hotter and flow faster in coronal holes than protons (Kohl et al., 1997Li et al., 1998;Cranmer et al., 1999). Enhanced perpendicular heating of ions compared Living Reviews in Solar Physics http://www.livingreviews.org/lrsp-2010-4 to protons has also been observed in streamers (Strachan et al., 2002;Uzzo et al., 2007). The magnitude of this effect is significant, but smaller in streamers than in coronal holes. Ulysses and Helios in-situ measurements in the heliosphere have shown that minor ions flow faster than protons by the local Alfvén speed, and are preferentially heated as well, and proton distributions often appear double-peaked with an average relative drift parallel to the background magnetic field (e.g., Marsch et al., 1982b;Feldman et al., 1996;Neugebauer et al., 1996).
The temperature anisotropy of protons deduced from remote sensing and in-situ observations of fast solar wind streams provides indirect evidence for the presence of the ion-cyclotron waves in coronal plasma, since the anisotropy can be produced by the resonant absorption of the ioncyclotron waves. Purely adiabatic expansion of the solar plasma is expected to result in an opposite effect: ⊥ < ‖ due to the conservations of magnetic moment of the expanding ions in the decreasing radial magnetic field (e.g., Marsch, 2006). However, ⊥ / ‖ > 1 is observed in the heliosphere (e.g., Marsch et al., 1982b;Gazis and Lazarus, 1982). In the past, several theories of ion-cyclotron resonance have been developed and applied to the heating of the solar corona and the solar wind (e.g., Axford and McKenzie, 1992;Tu and Marsch, 1997;Li et al., 1999;Hollweg, 2000;Hu et al., 2000;Cranmer, 2000;Hollweg and Isenberg, 2002). However, there are theoretical difficulties with the application of the ion-cyclotron mechanism for coronal heating, and its role is not yet fully understood (e.g., Cranmer, 2000;Isenberg, 2004). Most such theories may be classified as either fluid-like or quasi-linear kinetic models. The limitation of the fluid or quasi-linear kinetic models is the assumption of fixed-shape ion velocity distribution and quasi-linear limits (i.e., small magnetic fluctuation amplitude allowing simplified description of wave-particle interactions). In the hybrid models the electrons are treated as a fluid, and the ions are treated fully kinetically as particles. Hybrid simulations (see Section 2.3 below) allow relaxing many approximations used in the fluid, multi-fluid, and in linear or quasi-linear kinetic theory. The model is nonlinear, and can describe both the brief initial linear evolution of the plasma, as well as the nonlinear saturated state. Recently, new kinetic models of heating and acceleration of solar coronal plasma in inhomogeneous magnetic field by Alfvén waves were developed (Galinsky and Shevchenko, 2013a,b). The generation of the solar wind by parallel (Isenberg, 2012;Mecheri, 2013) and oblique Isenberg and Vasquez, 2011) ion-cyclotron waves were also studied. The possible role of kinetic Alfvén waves (KAW) (e.g., Voitenko and Goossens, 2006;Dwivedi et al., 2012) and Alfvén waves turbulence (Chandran, 2010;Chandran et al., 2011;Li et al., 2011;Cranmer and van Ballegooijen, 2012) on the acceleration and heating of the solar wind was recently considered.
Recently, Cranmer and van Ballegooijen (2005) solved the linearized wave equation for Alfvén waves in the heliosphere. The dependence on heliocentric distance of the frequency-integrated Alfvénic velocity amplitude obtained from the solution of the linearized Alfvén wave equation driven by a spectrum of transverse photospheric fluctuations with an amplitude of 3 km s -1 at the photosphere is compared to Alfvén wave amplitude inferred from spectroscopic observations (SOHO/SUMER and UVCS), IPS measurements, and in-situ data in Figure 3. The study shows that there is generally good qualitative agreement between observations and theoretical prediction of the Alfvén wave evolution in the heliosphere, even for a linearized model.
The turbulence in the solar wind magnetized plasma has been studied for decades in the past (see the review by Velli, 2003), and recently (Verdini et al., 2009;Chandran and Hollweg, 2009;Verdini et al., 2010;Markovskii et al., 2010) as the possible state that leads to cascade of energy from the observed large scale fluctuations and waves to small scale structures, down to dissipation scales that can heat the solar wind plasma. Observations of Alfvénic fluctuations in the solar wind by Helios and Ulysses spacecraft show that the turbulent energy carried by these fluctuations is distributed in frequency according to a power law, at high frequencies going as −5/3 , a Kolmogorov spectrum, while at lower frequencies the spectrum flattens to −1 (where is the fluctuation frequency) (Goldstein et al., 1995). Alfvénic turbulence is predominant in fast wind streams while in slow solar wind the turbulence is of more complex nature with low Alfvénicity (see the reviews by Tu and Marsch, 1995;Bruno and Carbone, 2013). Recent observations by ACE spacecraft of the solar wind protons at 1 AU indicate that the turbulent cascade rate agrees better with Kraichnan ( −3/2 ), rather than with Kolmogorov ( −5/3 ) rate (Vasquez et al., 2007). Similar results were seen by the Wind spacecraft (Podesta et al., 2006). Part of the fluctuating power at low frequencies can be attributed to propagating structures in the solar wind. However, there is strong evidence that the fluctuations are Alfvénic at frequencies of milli-Hertz and higher.
Recent observations by Hinode satellite show that Alfvénic fluctuations are the likely energy source that drives the solar wind (e.g., De Pontieu et al., 2007;Ofman and Wang, 2008;Hahn et al., 2012;Hahn and Savin, 2013). A review of observational evidence for propagating MHD waves in coronal holes that may accelerate the solar wind is found in Banerjee et al. (2011). Recently launched NASA's Solar Dynamics Observatory (SDO) provides unparalleled opportunity to look for the solar coronal wave spectrum over the entire disk of the Sun at high temporal and spatial resolution. The analysis of the data from the Atmospheric Imaging Assembly (AIA) onboard SDO will likely provide the constraint on the input wave spectrum that drives the solar wind. Although the observed spectrum is limited to the MHD frequency range, since the temporal resolution does not allow resolving high frequencies down to the gyroresonant scale (∼ kHz), the form of the spectrum is likely to provide clues on the relevant turbulent cascade processes.
The following level of solar wind plasma approximation is via hybrid models (see Section 2.3), that describe protons and other ions kinetically as particles, and electrons as neutralizing background fluid (Winske and Omidi, 1993). Hybrid simulations can represent more completely (than fluid model) and self-consistently the wave-particle interactions in the multi-ion solar wind magnetized plasma. The models can be used to describe the kinetic processes involved in heating by a spectrum of waves, and can be used to study the nonlinear and resonant interactions of the turbulent spectrum with the ions. Such numerical methods have the potential to model the wave-Living Reviews in Solar Physics http://www.livingreviews.org/lrsp-2010-4 particle interactions, and the corresponding velocity distribution and magnetic fluctuations in the nonlinear saturated state that can be compared to in-situ observations. Recently, one dimensional hybrid simulations of multi-ion solar wind plasma were used to study the heating by wave spectrum, beams, and the stability of solar wind multi-ion plasma (Liewer et al., 2001;Ofman et al., , 2002Xie et al., 2004;Lu and Wang, 2005;Li and Habbal, 2005;Hellinger et al., 2005;Ofman et al., 2005). Due to the local nature of the hybrid models, they require special treatment taking into account the global properties of the solar wind, such as expansion of the solar wind into the heliosphere (Liewer et al., 2001;Hellinger et al., 2005;Ofman et al., 2011). We will review some of the recent results of hybrid simulation models of solar wind plasma heating.
Two-dimensional hybrid models of homogeneous multi-ion plasma heating were also studied recently (Gary et al., , 2003Kaghashvili et al., 2003;Hellinger and Trávníček, 2006;Gary et al., 2006). However, only few studies considered the effect of time dependent wave fluctuations on solar wind plasma heating with 2D hybrid models (Ofman and Viñas, 2007;Ofman, 2010;Markovskii et al., 2010). It was found that the input wave spectrum can heat the ions by resonant interaction, as well as through non-resonant parametric decay instability of Alfvén waves (e.g., Araneda et al., 2007). It was also found that the presence of small scale inhomogeneity in the background plasma can enhance the heating by the high frequency Alfvén wave spectrum (Ofman, 2010).
This review is focused on selected wave-acceleration models (both, MHD and kinetic) of the solar wind. Other reviews of solar wind models were recently published (Hansteen and Velli, 2012;Cranmer, 2012

Model Equations
Below we provide the typical basic equations used in the three classes of models reviewed here: MHD, multi-fluid, and hybrid.

Single fluid MHD
The single-fluid, normalized visco-resistive MHD equations with gravity are (see, e.g., Priest, 1982;Ofman, 2005) where is the fluid density, v is the fluid velocity, B is the magnetic field, = ∇ × B is the current density, and is the plasma pressure. The normalization of the magnetic field is by the typical magnetic field 0 at the base of the corona, distances ( ) are normalization by the solar radius ⊙ , the density normalization is by the typical density 0 at the coronal base, the velocity normalization is by Alfvén speed = 0 /(4 0 ) 1/2 , and the time normalization is by the Alfvén time defined as = / , where is the typical lengthscale of the problem (for convenience, = ⊙ is used). The pressure is normalized by 0 = 0 0 , where is the Boltzmann constant, 0 = / ( ) is the typical number density, where is the proton mass, and is the average mass number density of the coronal plasma, and 0 is the typical temperature at the base of the corona. The Froude number is = 2 ⊙ /( ⊙ ), is the gravitational constant, ⊙ is the solar radius, ⊙ is the solar mass, = 2 2 /( 2 ) is the ratio of thermal to magnetic pressure, where = ( / ) 1/2 is the sound speed. The Lundquist number = / is the ratio of the typical Alfvén time to the resistive diffusion time ( = 4 2 / 2 , where is the resistivity, and is the speed of light), and F is the viscous force term (see Braginskii, 1965).
The heating and loss terms are ℎ , and , respectively, where ℎ is the coronal heating function (assumed, or obtained empirically), and represents the losses due to thermal conduction and radiation (for example, Landi and Landini, 1999;Colgan et al., 2008, for optically thin plasma). The polytropic index = 1 in isothermal plasma, = 1.05 in a commonly used polytropic model of the solar wind, and = 5/3 for solar wind models with explicit heating terms. Variable was used by Cohen et al. (2007) to match solar wind properties at 1 AU. The above set of equations is supplemented by the equation of state = , and by the solenoidality condition ∇ · B = 0. The visco-resistive single fluid MHD model was used recently to reproduce the global emission properties of the solar corona (Lionello et al., 2009;Downs et al., 2010).

Multi-fluid models
Here, we describe the multi-fluid equations and model utilized by Ofman (2004a) to model the fast solar wind in coronal holes accelerated by nonlinear MHD waves. Neglecting electron inertia ( ≪ ), relativistic effects ( ≪ ), assuming quasi-neutrality ( = + ), where is the charge number, the normalized three-fluid MHD equations can be written as Living Reviews in Solar Physics http://www.livingreviews.org/lrsp-2010-4 where the index = , (in Equation (9) = , , ), is the heat conduction term, is the heating term of each fluid, is the energy coupling term between the various fluids (Li et al., 1997;Ofman, 2004a), F = ∇ · Π is the viscous force term due to ions, where Π the viscous stress tensor and the Coulomb friction terms F ,coul given by Braginskii (1965), = 5/3 is the polytropic index of each species, is the mass number, and Ω = 0 is the normalized gyrofrequency.
The three-fluid equations are normalized by → / ⊙ , where ⊙ is the solar radius; → / ; The following parameters enter in the above equations: the Lundquist number; , = ( , ,0 / )/ 2 the electron and proton Euler number; = ( ,0 / )/ 2 the ion Euler number; the Froude number; is the atomic mass of species ; = 0 /(4 0 ⊙ ), ,0 is the normalization temperature, is the mass of the particles, 0 is the normalization magnetic field. The heat conduction term in Equation (9) is normalized by → ( ⊙ / 2.5 0 ). In Ofman (2004a) model the fast solar wind is produced by a broad band spectrum of waves. The linearly polarized Alfvén waves are driven at the base of the corona as follows: where = /2 , = −1 for the −1 spectrum, the discrete frequencies are given by = 1 + ( − 1)Δ , and the range is defined by Δ = ( − 1 )/( −1), where is the number of modes, Γ ( ) is the random phase that depends on the solar latitude . Typically, the frequencies are in the mHz range, the driving amplitude is few percent of , with an order of 100 modes used to model the desired spectrum (see Figure 4). As described by Ofman (2004a) dissipation of the waves occurs through viscous and resistive terms in the momentum and inductance equations, respectively. The dissipation coefficients used in that model are hyper-viscosity and hyper-resistivity, i.e., their values are much larger than the classical resistivity and viscosity, accounting empirically for kinetic and turbulent effects.
In addition to the waves, an empirical heating function is introduced in this model to heat the ions: where 0, ≡ 0, is the amplitude of the heat input in normalized units, and is the length scale of the heating in ⊙ . This is necessary, since the Alfvén wave spectrum constrained by available observations can not account for the observed (e.g., Kohl et al., 1997)  and heating of heavy ions. The asymptotic solar wind parameters (speed of various ion species, mass flux, temperature, etc.), can be matched by fitting the parameters of the heating function, combined with the parameters of the driving Alfvén wave spectrum.

Hybrid models
In the hybrid model the ions are represented as particles, neglecting collisions, while the electrons are described as a finite temperature massless fluid to maintain quasineutrality of the plasma. This method allows one to resolve the ion dynamics and to integrate the equations over many ion-cyclotron periods, while neglecting the small temporal-and spatial-scales of the electron kinetic motions. In the hybrid model, the three components of order million particle velocities are used to calculate the currents, and the fields in the 1D, 2D, or, in some cases, 3D grid. Note that each numerical particle represents large number of real particles, determined by the density normalization. The required number of particles per cell is determined by the required limitation on the overall statistical noise, and could be increased by an order of magnitude as needed. The following equations of motion are solved for each particle of species (k): where is the particle mass, is the charge number, is the electron charge, and is the speed of light. The electron momentum equation is solved by neglecting the electron inertia where = is used for closure, and quasi-neutrality implies = + , where is the number density of electron, protons, and ions, respectively. The above equations are supplemented with Maxwell's equations where the displacement current is neglected in non-relativistic plasma, and The field solutions are obtained on the 1D, 2D, or 3D grid, and the proton and ion equations of motions are solved as the particle motions respond to the fields at each time step. The method has been tested and used successfully in many studies. In Ofman and Viñas (2007) 2D study 128 × 128 grid with 100 particles/cell/species were used. The particle and field equations were integrated in time using a rational Runge-Kutta (RRK) method (Wambecq, 1978) whereas the spatial derivatives were calculated by pseudospectral FFT method. When non-periodic boundary conditions are applied finite difference method is used for the field solver. The hybrid model allows computing the self-consistent evolution of the velocity distribution of the ions that includes the nonlinear effects of wave-particle interactions without additional assumptions. Moreover, the hybrid model is well suited to describe the nonlinear saturated state of the plasma.
Since the hybrid models usually describe a region of several hundred ion inertial length ( = / ) across, the method is limited to model local small scale structures in the corona. For a typical solar plasma density of 10 4 cm -3 at 10 ⊙ and a simulation box with a side of 440 we get about 1000 km for the extent of the simulated region in each dimension. At 1 AU the plasma density is much lower and the modeled region covers about 45 000 km in each dimension. A way to overcome the computational limitation to small scales is to use an 'expanding box' model (e.g., Grappin and Velli, 1996;Liewer et al., 2001;Hellinger et al., 2005). This approach employs transformation of variables to the moving solar wind frame that expands together with the size of the parcel of plasma as it propagates outward from the Sun. In particular, it is assumed that a small packet of plasma of length ≪ 0 and width ( ), where ( )/ 0 ≪ 1 expands in the lateral direction only as it moves away from the Sun at constant speed 0 . The initial distance 0 is ( ⊙ ). Thus, the modeled regions position is ( ) = ⊙ + 0 , and the width ( ) = ( )/ 0 . Using these transformations the coordinates are transformed as ′ = − ( ), ′ = / ( ), and ′ = / ( ), and the equations of motions together with the field equations are transformed to the moving and expanding frame. Although, the method requires several severe simplifying assumptions (i.e., lateral expansion only, constant solar wind speed) and approximations (the original spherical coordinates and the mean magnetic field transformed to new coordinates using second order expansion (see Liewer et al., 2001) to remain tractable, it provides qualitatively good description of the solar wind expansions, thus connecting the disparate scales of the plasma in the various parts of the heliosphere.

Selected Model Results
In this section a brief overview of solar wind model results is given. Here, we concentrate on the results of wave driven 2.5D MHD, 2.5D multi-fluid models, as well as 1D and 2D hybrid models. In the reviewed models the waves are included explicitly, fully resolved, and their damping or resonant absorption was calculated explicitly. This allows more accurate description of the physics and interaction between the waves and the solar wind plasma than in WKB models, or models that parameterize the propagation and dissipation of the waves. The results of a global 3D MHD solar wind model computed with SWMF that incorporates the effects of Alfvén wave heating and acceleration in the WKB approximation is shown in Figure 5. The figure shows a cut of the 3D MHD model results in the meridional plane for an idealized tilted-dipole magnetic field configuration. The formation of the bi-model solar wind with slow wind in the streamer belt and fast at higher latitudes is evident in the radial outflow velocity values. Wave models in the single fluid MHD are limited to frequencies much smaller than the proton gyrofrequency and correspondingly to wavelengths that are much larger than the proton gyroradius. The acceleration of the solar wind plasma by the waves due to momentum transfer (wave reflection) and gradient of the wave pressure is modeled. The dissipation of the waves by Ohmic and viscous dissipation terms is included. The multi-fluid models allow including waves with frequencies in the MHD and in the range of proton and ion gyrosresonant frequency. The multi-ion cyclotron resonant dispersion of waves is reproduced by this model even in the linear regime. It has been shown that multi-fluid dispersion relation is equivalent to Vlasov's dispersion for cold plasma (for example, see Ofman et al., 2005). In addition, the multi-fluid models can include separate heating and dissipation processes for electrons, protons, and each ion species with different dissipation coefficients. The fluids are coupled through collisional energy exchange terms and Coulomb friction, Living Reviews in Solar Physics http://www.livingreviews.org/lrsp-2010-4 and through electromagnetic interactions. The multi-fluid models provide the next level of plasma approximation between the MHD and the kinetic descriptions.
The hybrid simulations extend the modeled physics of the solar wind plasma to even smaller scale in time and space to the kinetic regime, and the wave frequencies in the ion-and protongyroresonant scale are resolved. In addition, ion velocity space instabilities and ion kinetic processes are modeled fully. The hybrid models are limited to waves with frequencies below the electron gyroresonant frequency, since the electrons are treated as fluid in these models. Other limitations of these models are outlined in Winske and Omidi (1993). The 1D hybrid models are limited to parallel propagating waves in one spatial direction, or oblique waves with fixed angle of propagation. The more general 2D hybrid models include the description of waves with arbitrary propagation direction and can be used to model inhomogeneous plasma in two spatial dimensions.

Fast solar wind in coronal holes
It is well known that thermally driven Parker's solar wind model with typical coronal temperature of 1 -2 MK can produce the slow solar wind asymptotic speed of about 400 km s -1 , but can not explain the fast solar wind that is observed to reach 800 km s -1 within 10 ⊙ and is associated with coronal holes with typical temperatures < 1 MK (Aschwanden, 2004). The common approach is to include an additional source of momentum in the MHD equations, such as Alfvén waves as an empirical WKB momentum addition term (e.g., Usmanov et al., 2000). This approach was recently extended to include the effects of turbulence dissipation in a global wave driven solar wind model, and implemented in the Space Weather Modeling Framework (SWMF) (Tóth et al., 2005) coronal 3D MHD code (Evans et al., 2012;Sokolov et al., 2013;Oran et al., 2013). Two-temperature Alfvén wave driven fast solar wind models were also developed in SWMF (van der Holst et al., 2010). Lau and Siregar (1996) studied the acceleration of the solar wind by resolved nonlinear Alfvén waves in 1.5D MHD model. Ofman and Davila (1997) were the first to use the single fluid 2.5D MHD model to study resolved Alfvén wave driven fast solar wind in a coronal hole. In their model the Alfvén waves were launched at the solar boundary of a coronal hole, and were resolved throughout the coronal hole to 40 ⊙ . The acceleration of the solar wind occurs through momentum transfer from the waves to the solar wind plasma. The heating of the solar wind plasma was not included explicitly in this model and an isothermal approximation was used ( = 1). However, wave dissipation does occur through resistive dissipation with finite value of . In Figure 6, a snapshot of the spatial dependence of the solutions in terms of / 1/2 , , , and is shown at = 255 = 32.5 h. The velocities and / 1/2 are in units of = 1527 km s -1 , and the density is in units of 10 8 cm -3 . The monochromatic Alfvén waves launched in this model are evident in and in / 1/2 . The nonlinear longitudinal waves produced by the gradient of the compressions associated with the Alfvén wave, 2 , are evident in and . The large amplitude, long wavelength compressional velocity and density fluctuation propagate in-phase (Ofman and Davila, 1998). Ofman and Davila (1998) found that low-frequency (0.35 mHz) Alfvén waves with amplitude of 46 km s -1 can produce the fast solar wind in coronal holes. Grappin et al. (2002) were the fist to study resolved Alfvén waves driven wind that include both, closed and open field regions using 2.5D MHD model. They found that onset of Alfvén wave flux in one hemisphere generates a stable global circulation pattern in the closed loops region that can lead to global north-south asymmetry of the solar corona.
In Figure 7 a cut through the center of the coronal hole is shown. The and solar wind velocity components are shown for two Alfvén wave driving frequencies, and the green curve shows Parker's isothermal solar wind solution. It is evident that the low frequency waves (f = 0.35 mHz) lead to significant acceleration of the fast solar wind above Parker's solution, and produce the fast solar wind far from the Sun. The higher frequency waves provide acceleration close to the Sun below 10 ⊙ .
Living Reviews in Solar Physics http://www.livingreviews.org/lrsp-2010-4 The 2.5D model discussed above includes only the coronal part of the solar wind, with the driving Alfvén waves applied at the lower coronal boundary. Recently, Suzuki and Inutsuka (2005) modeled the acceleration of the fast solar wind by Alfvén waves from the photosphere to 0.3 AU using a 1.5D model (Figure 8). This approach allowed connecting directly photospheric motions of observationally constrained magnitude to solar wind speed at 0.3 AU. Although the model does not include the effects of cross-field gradients, the model demonstrates that sufficient Alfvén wave energy flux reaches the corona to accelerate the solar wind. The 1.5D model results compare favorably to IPS observation (Grall et al., 1996;Canals et al., 2002) and to SOHO observations (see, Suzuki and Inutsuka, 2005, for the details).
Living Reviews in Solar Physics http://www.livingreviews.org/lrsp-2010-4 Living Reviews in Solar Physics http://www.livingreviews.org/lrsp-2010-4 3.2 Fast solar wind: 2.5D multi-fluid models In Figures 9 and 10 we show the results of the 3-fluid model of the Alfvén wave driven fast solar wind in a coronal hole obtained by Ofman (2004a). In this model a broad band spectrum of Alfvénic fluctuations was applied at the lower coronal hole boundary. The fast solar wind was produced by acceleration and heating with the spectrum of Alfvén waves, that were fully resolved in the model. In Figure 9 the Alfvén waves are evident in and the accelerating solar wind in velocity components for He ++ ions (left panels) and protons (right panels) at = 114 are shown. Note that compressive fluctuations are also seen in due to the local variation of wave pressure gradient. The velocity is in units of 1527 km s -1 . The distance is in units of ⊙ , and the latitude is in radians. The typical form of the magnetic fluctuations obtained with the 3-fluid model at = 18 ⊙ is shown in Figure 10. It is interesting to note that the −1 spectrum launched at the base of the coronal hole results in −2 spectrum at larger distances. The steepening of the magnetic fluctuations spectrum is expected due to turbulence and dissipation that affects shorter wavelengths Living Reviews in Solar Physics http://www.livingreviews.org/lrsp-2010-4 Wave Modeling of the Solar Wind 21 and correspondingly higher frequencies more than the long wavelength (low frequencies) fluctuations. The power law dependence is close to Kolmogorov's turbulent power spectrum of −5/3 . At frequencies higher than 200 −1 the spectrum steepens due to increased dissipation. In Figure 11 the -averaged outflow speeds of protons, He ++ , and O 5+ ion fluid are shown for four sets of model parameters. The parameters are (a) 0 = 0.5, 0He++ = 12, = 0.034, (b) 0 = 0.5, 0O5+ = 10, = 0.034, and (c) 0 = 0.0, 0He++ = 12, = 0.05, where 0 is the heating rate per particle for an empirical heating term used in Ofman (2004a), and is the amplitude of the Alfvén wave spectrum. The corresponding temperatures and densities are shown in Figure 12. In Figure 11a the solutions of the 3-fluid model with empirical heating term [Equation (12)] in addition to Alfvén wave spectrum are shown. The heating term parameters were chosen to match the observed fast solar wind speed, and the faster outflow of He ++ ions compared to protons, observed at 0.3 AU and beyond with Helios and Ulysses spacecraft (Marsch et al., 1982a,b;Feldman et al., 1996;Neugebauer et al., 2001). In Figure 11b the O 5+ ions were included as the third fluid, and the heating function parameters for O 5+ ions were adjusted to get faster than proton outflow. In Figure 11c the same heating per particle was deposited in protons and He ++ ions. Evidently, in this case the He ++ ions outflow speed is slower than the proton outflow speed, contrary to observations. In Figure 11d the solar wind protons are accelerated and heated solely by the Alfvén wave spectrum (i.e., 0 = 0). This was achieved by increasing the input wave amplitude, compared to the values used in Figures 11a-c. Note that the temperature structure of protons and O 5+ ions as seen in Figure 12b is in qualitative agreement with SOHO/UVCS observations (Kohl et al., 1997;Cranmer et al., 1999;Antonucci et al., 2000) (at present there are no observations of He ++ temperature in this region). The model shows that in all cases electron heating can be achieved by thermal coupling between electrons and protons alone (through the thermal coupling term in Equation (9)).  Living Reviews in Solar Physics http://www.livingreviews.org/lrsp-2010-4

Leon Ofman
In spectroscopic observations of emission lines the observed finite line width is the result of broadening by Doppler shift due to the motion of the emitting ions in the line of sight. The motions are usually attributed to two components: (1) thermal or kinetic motions due to the finite width of the ion velocity distribution; (2) non-thermal motions, arising from any unresolved macroscopic motions of the plasma in the line of sight. The effect of observationally unresolved Alfvén waves on the apparent emission line widths was modeled with the 3-fluid model by Ofman and Davila (2001) and Ofman (2004a). These models allow separating the contribution of waves to the observed line profiles. In Figure 13 the Doppler-broadened emission line, resulting from combined thermal and non-thermal motions calculated with the 3-fluid model is shown. The solid line shows the simulated emission line profile of protons at 4 ⊙ at a temperature of 3.5 MK, and the dashes show the emission line broadened by unresolved Alfvénic fluctuation. In Figure 14 the effective temperature calculated from the kinetic temperature and the Alfvénic wave motion contribution is shown for the wave-driven fast solar wind as a function of the heliocentric distance. The effective proton temperature is affected significantly by the non-thermal component, while the relative effect on He ++ is smaller close to the Sun. By comparing the results of the 3-fluid model to observations, it is possible to evaluate better the thermal and non-thermal motions for protons and heavier ions in the observational data.  protons was demonstrated in this model. In Figure 15 the evolution of the temperature anisotropy for protons and O 5+ ions is shown. It is evident that after ∼ 600 Ω −1 the perpendicular heating of the ions saturates at anisotropy level of ∼ 7, and the proton are not heated significantly. The level of saturated anisotropy is determined by the temperature dependent nonlinear balance between ion-cyclotron unstable ion velocity distribution that releases electromagnetic ion-cyclotron waves and the resonant absorption of magnetic fluctuations together with parallel heating of the ions. After inspecting the perpendicular and parallel temperatures of O 5+ at the end of the run it is evident that the heating was predominantly in the perpendicular direction .

1D hybrid models
In Figure 16 the velocity distribution of the protons and O 5+ ions is shown at the end of the run. It is evident that the proton velocity distribution is isotropic and the O 5+ ions are hotter in the perpendicular direction than in the parallel direction. The O 5+ velocity distribution is close to bi-Maxwellian with small non-Maxwellian features in the parallel velocity distribution, likely produced by the small parallel heating due to nonlinear compressive modes driven by the Alfvénic fluctuations spectrum.
The relaxation of O 5+ ion temperature anisotropy due to ion-cyclotron instability for the parameter range relevant to fast solar wind in coronal holes was recently studied using 1D hybrid model ) (see Figure 23). The study was motivated by SOHO/UVCS observations indicating large temperature anisotropy of O 5+ ions (Kohl et al., 1997;Cranmer et al., 1999). It was found that the scaling of the relaxed ⊥ / ‖ − 1 with the final ‖ (full circles) and the scaling of the relaxation time, rel , with the initial ‖ (circles) agree well with the theoretical scaling law −0.41 ‖ (Gary, 1993). The "x"'s mark the values ⊥ / ‖ − 1 at = 0. The enhanced O 5+ abundance relative to protons of 6 × 10 -4 in this model was implemented in oder to shorten the computation times. Similar result was found in 2D hybrid model by Gary et al. (2003)   waves in a multi-ion coronal plasma. The dispersion relation of parallel propagating Alfvén cyclotron waves in the multi-ions coronal plasma was determined using 1D hybrid model (see Figure 17) and compared with multi-fluid and Vlasov dispersion relation. It was found that the three methods are in good qualitative agreement in the weakly damped regime ( /Ω < 1). The ratio of the ion to proton fluid velocities perpendicular to the direction of the magnetic field was calculated for each wave modes for typical coronal parameters (see Figure 18). It was found that the O 6+ perpendicular fluid velocity exhibits strong (factor of 20 -100) enhancement and He ++ perpendicular velocity is enhanced by a factor of 3.5 -5 compared with proton perpendicular fluid velocity, in qualitative agreement with SOHO/UVCS observations of large perpendicular velocity of heavy ions in coronal holes (e.g., Kohl et al., 1997;Cranmer et al., 1999). The study demonstrated how the results of hybrid models can be used to better understand the observations of coronal ion emission. Recently, Araneda et al. (2007Araneda et al. ( , 2008 used Vlasov theory and one-dimensional hybrid simulations to study the effects of compressible fluctuations driven by parametric instabilities of Alfvéncyclotron waves. They found that field-aligned proton beams are generated during the saturation phase of the wave-particle interaction, with a drift speed somewhat above the Alfvén speed. This finding agrees with typically observed velocity distributions of protons in the solar wind that contain a thermal anisotropic core and a beam component (see the review by Marsch, 2006). The expanding box model (Grappin and Velli, 1996;Liewer et al., 2001)  hybrid models of H + -He ++ solar wind plasma heated by a spectrum of turbulent Alfvénic fluctuations and in solar wind plasma with super-Alfvénic ion relative drift Maneva et al., 2013). In particular, Maneva et al. (2013) studied the turbulent heating and acceleration of He ++ ions by an initial self-consistent spectra of Alfvén-cyclotron waves in the expanding solar wind plasma using 1.5D hybrid simulations. They found that the He ++ ions are preferentially heated by the broad-band initial spectrum, resulting in much more than mass-proportional temperature increase (see Figure 19). Maneva et al. (2013) also found that the differential acceleration of protons and He ++ ions depend on the amplitude and spectral index of the magnetic fluctuation, while solar wind expansion suppresses the differential streaming. They also find that the expansion leads in general to perpendicular cooling for protons and alphas. However, the cooling effect of the expansion is small and the waves provide sufficient heating, maintaining significant temperature anisotropy, in agreement with observations. Inspection of the proton and alpha velocity distribution in the ‖ -⊥ plane shows the formation of non-Maxwellian features due to the effects of the broad band spectrum, such as perpendicular broadening (i.e., temperature anisotropy), as well as formation of a population of accelerated particles by the waves (see

2D hybrid models
The 2D hybrid codes solve similar set of equations as in the 1D hybrid codes but in two spatial dimensions. This allows an additional degree of freedom for particle motions, and the wave propagation is not limited to parallel propagating waves, allowing oblique propagation. In addition, the parallel magnetic field component does not have to be constant in order to conserve ∇ · B = 0. As a result, a broader range of possible wave modes, wave-particle, and wave-wave interactions are included in the 2D model compared to 1D model. Obviously, the 2D models are computationally intensive, and may require parallel processing for similar resolution in 2D and similar numbers of particles per cell as in the 1D models that can be run on a desktop workstation.
The 2D hybrid models have been used extensively in the past to model successfully the electromagnetic interactions in magnetized plasmas (McKean et al., 1994;Gary et al., 1997;Daughton et al., 1999;Gary et al., 2000Gary et al., , 2001Gary et al., , 2003Ofman et al., , 2002Xie et al., 2004;Ofman and Viñas, 2007). Comparisons between one-and two-dimensional hybrid simulations often show qualitative agreement in the ion response (Winske and Quest, 1986;Ofman and Viñas, 2007). In addition to allowing oblique waves, the 2D code allows including spatial inhomogeneity of plasma density perpendicular to the magnetic field, as well as divergent magnetic field geometry. These features are needed to describe solar wind acceleration and heating more consistently with corona conditions.
Recently, Ofman and Viñas (2007) studied the heating and the acceleration of protons, and heavy ions by a spectrum of waves in the solar wind, as well as the nonlinear influence of heavy ions on the wave structure using the 2D hybrid model. They considered for the first time the heating and acceleration of protons and heavy ions by a driven input spectrum of Alfvén/cyclotron waves, and by heavy ion beam in the multi-species coronal plasma in two spatial dimensions. They found that in the homogeneous plasma the ion beams heat the ions faster than the driven wave spectrum constrained by solar wind parameters, and produce temperature anisotropy with ⊥ > ‖ in qualitative agreement with observation. The beam-heating model requires that the beam speed is larger than the local Alfvén speed. Since any reconnection process produces Alfvénic beams as an exhaust (e.g., Priest, 1982;Aschwanden, 2004), the beams could readily become super-Alfvénic as the plasma moves to regions of lower local Alfvén speed. Since the threshold of beam stability is the Alfvén speed, it is possible that remnants of this process that takes place close to the Sun in the acceleration region of the solar wind are seen in proton data beyond 0.3 AU (Marsch, 2006). Below, the results obtained recently by Ofman and Viñas (2007) are reviewed. Figure 21: Comparison of 1D and 2D model results. The evolution of O 5+ temperature anisotropy calculated with the 1D hybrid (dashes) and 2D hybrid (solid) models show good agreement (Ofman and Viñas, 2007).
Living Reviews in Solar Physics http://www.livingreviews.org/lrsp-2010-4 Leon Ofman Ofman and Viñas (2007) compared the evolution of O 5+ ion anisotropy relaxation by ion cyclotron instability, by modeling the coronal plasma with both 1D and 2D hybrid codes. In Figure 21 the results of a 1D and 2D hybrid models runs are shown. The initial temperature anisotropy was set to 50 in both cases with parallel temperature of 1.4 × 10 6 K. It is evident that the temperature anisotropy has relaxed to similar values in 1D and 2D runs, close to the marginally stable value of ∼ 10 obtained for the parameters used with the Vlasov stability analysis (Gary, 1993). The agreement that was found between 1D hybrid and 2D hybrid evolution is consistent with the Vlasov dispersion relation that shows maximal growth of the ion cyclotron instability for parallel propagating modes. Figure 22: Results of the parametric study of He ++ anisotropy relaxation obtained with 2D hybrid code by Gary et al. (2003). The parameters were / = 0.05 with initial / ‖ = 1.0, ‖ / ‖ = 4.0, and isotropic protons. The crosses correspond to = 0, the squares indicate plasma parameters at saturation of the fluctuating magnetic fields, and the dots represent later times. The dashed line indicates the best fit of the anisotropies at Ω = 400 (adapted from Gary et al., 2003).
In Figure 22 the results of the parametric study of He ++ anisotropy relaxation obtained with 2D hybrid code by Gary et al. (2003) is shown. In that study the parameters were / = 0.05 with initial = ‖ , ‖ / ‖ = 4.0, and isotropic protons. The initial anisotropy of He ++ was chosen to maximize the linear growth rate of ion-cyclotron instability. The crosses show the anisotropy at = 0, while the squares show the anisotropy at magnetic energy saturation, and the dots represent later times. The dashed line shows the best fit of the anisotropies at Ω = 400, that produces the scaling law ⊥ / ‖ − 1 = 0.71/ 0.45 . Note the qualitative agreement between the 2D hybrid study for He ++ anisotropy relaxation, the 1D hybrid study for the O 5+ anisotropy relaxation shown in Figure 23, and the scaling law obtained analytically (Gary, 1993). Gary et al. (2003) have shown that the model results are consistent with Ulysses in-situ observations of solar wind protons and He ++ ions. Ofman and Viñas (2007) found that the perpendicular heating occurs for the beam-driven instability, which quickly saturates nonlinearly, and due to the driven spectrum of waves. In the driven wave spectrum case the amplitude of the magnetic field fluctuations was / 0 = 0.06, and the frequency range of the driver was below the proton gyroresonance. It was found (see Figure 24) that the O 5+ anisotropy grows quickly (within 400 Ω −1 ) to ⊥ / ‖ ≈ 4, and than saturates nonlinearly remaining in the range 4 -5 throughout rest of the evolution. The "x"'s mark the values ⊥ / ‖ − 1 at = 0. The enhanced O 5+ abundance of 6 × 10 -4 in this parametric study leads to shorter computation times . range of the wave spectrum included the O 5+ ion resonant frequency at rest. The anisotropy of the protons remains close to unity throughout the run. No significant net drift was found between the protons and O 5+ ions in the wave driven case. Similar results were obtained for He ++ ions.
The non-Maxwellian features of the ion velocity distribution are evident in the perpendicular to the magnetic field phase space plane. When the initial distribution is drifting Maxwellian with the drift velocity = 1.5 , the perpendicular velocity distribution of the O 5+ is shell-like, with decreased phase-space density in the central part of the distribution, compared to the perimeter. When the initial drift velocity was increased to 2 , the shell like structure of the phase-space density of O 5+ ions becomes even more apparent (Figure 25). It is interesting to note, that the He ++ perpendicular velocity distribution for the drifting case is nearly bi-Maxwellian, and does not exhibit the shell structure.
Recently, Ofman (2010) expanded this study and considered in a parametric study the effect of inhomogeneous background density on the heating by high frequency circularly polarized Alfvén waves with, and without drift between the protons and heavier ions. Ofman (2010) found that the inhomogeneity, and the drift lead to increased heating of the solar wind ions, compared to the homogeneous case, and the spectrum of magnetic fluctuations steepens beyond Kolmogorov's slope of -5/3. In Figure 26 the magnetic energy fluctuations spectrum obtained in the 2D hybrid simulation with inhomogeneous background density is shown. The dashed-dotted line shows the best fit power law to the spectrum in the regions where the slope did not change considerably.
Living Reviews in Solar Physics http://www.livingreviews.org/lrsp-2010-4 Ofman (2010) found that in the low density region the slope was m = -1.66 for wave driven case, and m = -1.81 for the beam driven case. However, in the high density region the slopes were m = -2.53 for the wave driven case, and m = -2.80 for the beam driven case, indicating enhanced dissipation due to the refraction of Alfvén waves, and the generation of small scale magnetosonic fluctuations that dissipate more effectively than Alfvénic fluctuations. Ofman et al. (2011) explored additional forms of background inhomogeneity on the magnetosonic drift instability and solar wind plasma heating by a spectrum of Alfvénic fluctuations. The expansion of the distant solar wind plasma at 0.3 AU and beyond and the generation of the associated kinetic instabilities and waves was considered in 2D hybrid models by Trávníček (2011, 2013).

Open Questions and Challenges
Although significant progress was made in observing and modeling of the solar wind over the past decades, there are still several important questions that are unanswered. This situation stems from the lack of unambiguous observations, which point to a specific physical mechanism for coronal heating, and solar wind accelerations, as well as due to the limitations of present models and theories. In particular, the following questions remain open: 1. What is the exact physical mechanism that produces the fast and slow solar wind? This question relates directly to the question of coronal heating mechanism.
2. What is the role of waves (in a broad frequency range from kinetic to MHD) in the acceleration and heating of the solar wind?
3. What is the role of density inhomogeneity and small scale turbulence (cascade) in the heating and the acceleration of the fast and slow solar wind?
4. How are the non-Maxwellian velocity distributions of protons and ions in the solar wind formed?
5. What determines the heavy ion composition (i.e., elemental abundance and charge states, see Zurbuchen, 2007) of the fast and slow solar wind?
6. What is the role of electrons in solar wind acceleration and heating? 7. How does Earth's global space environment respond to solar wind variations?
The works reviewed here bear on the first five questions, showing that MHD waves with a given spectrum provide plausible acceleration mechanism of the fast solar wind in coronal holes, and heating of coronal plasma may occur through resonant and non-resonant dissipation of the waves energy. However, the fluid models do not provide the kinetic details of the dissipation processes, and the hybrid models show only limited aspects of the resonant dissipation processes. The formation and the evolution of the ion velocity distribution of the solar wind plasma is not modeled in detail from the Sun to 1 AU. Multi-fluid models address limited aspects (i.e., within the ion-fluids approximation) of the compositional variation of the solar wind in open and closed structures (Ofman, 2000(Ofman, , 2004a. In the reviewed models the electrons are only studied as a fluid and their role in solar wind heating and acceleration only includes basic aspects (i.e., heat conduction, collisions coupling to ions). The last question can be addressed with global models that include the Earth's magnetosphere and ionosphere. However, the study of the kinetic processes that are at the roots of the solar-wind-magnetosphere-ionosphere interactions is far from complete.
The above questions are on the forefront of current research, and the answer can be obtained by combination of improved observations and modeling. A possible way to answer these questions is by obtaining in-situ measurements of the solar wind plasma in the region close to the Sun, where the acceleration and heating processes are still significant (McComas et al., 2007). This is the goal of the future European Solar Orbiter, and NASA's Solar Probe Plus missions.

Summary and Discussion
Satellite observations provide ample evidence for the presence of low-frequency (MHD) waves in the solar corona and the solar wind. The presence of ion-cyclotron waves and beams is evident as well in in-situ observations at 0.3 AU and beyond, and deduced from remote-sensing spectroscopic observations. Motivated by these observations wave-driven solar wind acceleration and heating models were developed with various degrees of approximation. The formation and the effect of beams on solar wind plasma heating were studied with hybrid models. In the present paper we have reviewed several such models of solar wind acceleration and plasma heating. The emphasis in this review is on wave driven models with fully resolved wave spectrum in 2.5D MHD, 2.5D multi-fluid models, and in hybrid kinetic 1D and 2D models.
Thermal conduction alone can not explain the acceleration of the solar wind to fast wind speed for plasma temperatures of 1 -2 MK commonly deduced from observations in open magnetic structures. The 2.5D MHD and multi-fluid models show that Alfvén wave spectrum in the MHD frequency range (millihertz) accelerates the fast solar wind to the observed speed of ∼ 800 km s -1 and provides the necessary energy to heat the solar wind. The advantage of the WKB approximation is that it allows incorporating the effects of Alfvén wave heating and acceleration in global 3D MHD models. However, the models that include fully resolved waves provide more accurate and realistic account of the interaction between the waves and the solar wind plasma than the WKB approximation and the MHD models that use ad-hoc heating function, momentum input, or variation of the polytropic index with distance from the Sun. The main limitations of the wave-driven solar wind MHD models are that the heating is described by Ohmic and viscous dissipation with empirical dissipation coefficients, and the exact kinetic process that underlay the fluid description can only be modeled in detail by kinetic approach.
Multi-fluid models extend beyond MHD by providing insights on the compositional variation of the solar wind plasma, on separate heating processes for electrons, protons, and heavy ions, and on the interactions between the various plasma constituents. The results of multi-fluid models are compared directly with observations of the coronal emission, consisting of ion emission lines, and white light polarized brightness that comes from electron Thompson scattering. These comparisons provide more stringent observational constrains on solar wind models than can be achieved with single fluid MHD, since all modeled particle species must conform to the observed properties (e.g., electron temperature, proton temperature, relative abundance of heavy ions in various magnetic structures, wave signatures in separate fluids, etc.), and the various fluids are coupled through Coulomb and electromagnetic interactions.
The 1D and 2D hybrid models provide the next level of physical modeling, and are reliable tools that have been tested and used for decades to study ion kinetic processes in space plasmas. The reviewed studies concentrate on the resonant dissipation of wave spectrum in the multi-ion solar wind plasma and include the effects of beams. The models show that the high frequency waves in the proton and ion gyrosresonant frequency range can heat the solar wind heavy ions preferentially and anisotropically and produce the anisotropic ion velocity distributions deduced from observations. High-amplitude waves can lead to beam formation, while solar wind expansion can lead to perpendicular cooling of the ions. The hybrid models show that heating can be enhanced further by the instability of super-Alfvénic beams of heavy ions. The reviewed studies show that protons are not heated significantly by these waves due to resonant absorption by heavier ions. Thus, the spectrum of waves that heats and accelerates the solar wind must contain both, lowfrequency (non-resonant) and high-frequency Alfvén waves. The hybrid models do not include the kinetics of electrons, and their possible role in solar wind energy balance and the dissipation of low-frequency waves is not modeled beyond the fluid description.
The planned NASA Solar Probe Plus mission and the European Solar Orbiter missions will provide new measurements in the unexplored region of the inner heliosphere. In particular, in-situ Living Reviews in Solar Physics http://www.livingreviews.org/lrsp-2010-4 Wave Modeling of the Solar Wind 39 measurement of non-Maxwellian features in proton, ion, and electron velocity distributions such as anisotropy and beams, and measurement of magnetic fluctuations spectrum in the acceleration region of the solar wind close to the Sun will provide the necessary information that will improve our understanding of solar wind acceleration and heating. These measurements will provide improved constraints for future theoretical studies and numerical models of solar wind plasma heating and acceleration for all levels of plasma approximations.