Alfvén wave solar model (AWSoM): proton temperature anisotropy and solar wind acceleration
|
|
- Arthur Caldwell
- 6 years ago
- Views:
Transcription
1 doi: /mnras/stv2249 Alfvén wave solar model (AWSoM): proton temrature anisotropy and solar wind acceleration X. Meng, B. van der Holst, G. Tóth and T. I. Gombosi Center for Space Environment Modeling, University of Michigan, Ann Arbor, MI 48109, USA Accepted 2015 September 28. Received 2015 September 25; in original form 2015 April 24 ABSTRACT Temrature anisotropy has been frequently observed in the solar corona and the solar wind, yet poorly represented in computational models of the solar wind. Therefore, we have included proton temrature anisotropy in our Alfvén wave solar model (AWSoM). This model solves the magnetohydrodynamic equations augmented with low-frequency Alfvén wave turbulence. The wave reflection due to Alfvén sed gradient and field-aligned vorticity results in turbulent cascade. At the gyroradius scales, the apportioning of the turbulence dissipation into coronal heating of the protons and electrons is through stochastic heating. This par focuses on the impacts of the proton temrature anisotropy on the solar wind. We apply AWSoM to simulate the steady solar wind from the corona to 1 AU using synoptic magnetograms. The Alfvén wave energy density at the inner boundary is prescribed with a uniform Poynting flux r field strength. We present the proton temrature anisotropy distribution, and investigate the firehose instability in the heliosphere from our simulations. In particular, the comparisons between the simulated and observed solar wind prorties at 1 AU during the ramping-up phase and the maximum of solar cycle 24 imply the importance of addressing the proton temrature anisotropy in solar wind modelling to capture the fast solar wind sed. Key words: methods: numerical solar wind. 1 INTRODUCTION In magnetized plasmas, the plasma temrature along and rndicular to the magnetic field can be different due to insufficient heat exchange between the parallel and rndicular thermal motions. This feature is referred to as temrature anisotropy or thermal pressure anisotropy, and commonly exists in space plasmas. A theoretical foundation is provided by the Chew Goldberger Low approximation (Chew, Goldberger & Low 1956). A large amount of in situ and remote observational data have shown that the plasma temrature is anisotropic in the solar corona (SC; Kohl et al. 1998, 2006; Lietal.1998; Antonucci, Dodero & Giordano 2000; Telloni, Antonucci & Dodero 2007) and the solar wind (Hundhausen 1968; Feldman et al. 1974; Marschetal. 1982; Gary, Skoug & Steinberg 2001; Gary, Goldstein & Neugebauer 2002; Kasr, Lazarus & Gary 2003; Marsch, Ao & Tu 2004; Hellinger et al. 2006). Since electrons transport much more quickly than protons, the proton temrature anisotropy is often much larger than the electron temrature anisotropy in the solar wind. Though being observed frequently, the xingm@umich.edu Present address: Jet Propulsion, 4800 Oak Grove Dr., Mail Stop , Pasadena, CA role of the temrature anisotropy, particularly the proton temrature anisotropy, in the solar wind is not fully understood. Solar wind modelling has become a very active research area (Riley, Linker & Mikić 2001; van der Holst et al. 2005; Feng, Zhou & Wu 2007; Kleimann et al. 2009; Watermann et al. 2009; Gressl et al. 2014). In particular, to address the proton temrature anisotropy in the solar wind, a lot of efforts have been made. The earliest work can be dated back to about 40 years ago. Leer & Axford (1972) develod a one-dimensional (1D) steady-state solar wind model with isotropic electron temrature and anisotropic proton temrature that allows for extended coronal proton heating. The model itself does not solve for the magnetic field, but the solar wind solutions are obtained for the cases of a purely radial and a spiral magnetic field, and the solutions at 1 AU are reasonable. At about the same time, Whang (1972) also develod a 1D steadystate solar wind model including proton thermal anisotropy, which assumes the same isotropic temrature for electrons and protons within a radial distance of 0.4 AU, while different temratures for electrons and protons, as well as proton temrature anisotropy beyond 0.4 AU. The calculated solar wind solution at 1 AU agrees well with observations. Both models are pioneers in simulating proton temrature anisotropy in the solar wind. More work has been done after 1990s. Hu, Esser & Habbal (1997) presented a 1D time-dendent fast solar wind model with isotropic electron temrature and anisotropic proton temrature. C 2015 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society
2 3698 X. Meng et al. The model includes momentum and heat input to the solar wind by Alfvén waves, additional momentum input to the protons and additional heat input to both electrons and protons. The high-sed solar wind solutions they obtained match most of the empirical constraints from observations. This model was further extended to solve the 16-moment bi-maxwellian equations by Li (1999), who has found that the inclusion of proton parallel and rndicular heat flux densities greatly affects proton temrature anisotropy. More recently, Li et al. (2004) reported the first two-dimensional (2D) Alfvén wave turbulence-driven solar wind model with proton temrature anisotropy. They obtained solutions for both fast and slow solar wind, and found that the average proton temrature in the anisotropic case is lower than in the isotropic case. Chandran et al. (2011) develod a 1D solar wind model including proton temrature anisotropy, pitch-angle scattering from mirror and firehose instabilities, and kinetic Alfvén wave turbulence, based on theories of linear wave damping and non-linear stochastic heating. They have found consistency between their model results and a number of measurements. All these studies provide valuable understanding to the modelling of solar wind with anisotropic proton temrature. However, these 1D or 2D models still cannot reveal the complete picture of the spatially and temporally varying solar wind, which could only be provided by time-dendent 3D global models of the corona and heliosphere. For the first time, we present the extension of our 3D global magnetohydrodynamic (MHD) Alfvén wave solar wind model AWSoM to include proton temrature anisotropy. The model assumes isotropic electron temrature; thus, it is a three-temrature (3T) model. Van der Holst et al. (2014) have presented the lower coronal modelling with the AWSoM, and found good agreement between simulated and observed multi-wavelength EUV images. The continuing work, being presented in this par, focuses on the solar wind simulations up to 1 AU. More scifically, we apply the model to simulate the corona and the inner heliosphere (IH) for two Carrington rotations (CRs) during the solar cycle 24. We look into the proton temrature anisotropy and plasma instabilities in the simulations, and validate the simulated solar wind prorties at 1 AU against the observed values made by the STEREO and WIND satellites, as well as from the OMNI data. By comparing the simulation results with the results from the two-temrature (2T) model, i.e. the model with isotropic proton and electron temratures, we find that including proton temrature anisotropy in the model leads to slower and more realistic solar wind seds. The following content of the par is divided into four sections. In Section 2 we describe the equations implemented, with emphasis on the treatment of proton pressure anisotropy limiting through instabilities. Sections 3 and 4 present the solar wind simulation results for CR2107 and CR2123, resctively. Section 5 concludes the par and discusses possible future work. The par also contains apndices about the plasma instabilities we considered in the model and on the electron pressure equation. 2 MODEL DESCRIPTION The governing equations for the AWSoM are ρ + (ρu) = 0, (1) t ρu + [ρuu + p I + p e I + (p p )bb t 1 μ 0 ( BB B2 2 I )] + ( ) w+ + w 2 = ρ GM r 3 r, (2) B t p t + ( u B) = 0, (3) + (p u) + 2p b (b )u = (γ 1) p e p τ ei + (γ 1)Q p + δp δt, (4) ( ) ( ) p p + t γ 1 γ 1 u + p ( u) + (p p )b (b )u = p e p + Q p, (5) τ ei ( ) ( ) + t γ e 1 γ e 1 u + p e ( u) w ± t = q e + p p e τ ei Q rad + Q e, (6) + (uw ± ± V A w ± ) w ± u = R w w + Ɣ ± w ±, (7) where ρ and u are the mass density and velocity, B is the magnetic field, b = B/ B is defined as the unit vector along the magnetic field and μ 0 is the rmeability of vacuum. G is the gravitational constant, M is the solar mass and r is the position vector relative to the centre of the Sun, with r = r. p represents the proton pressure component parallel with the magnetic field, while p is the rndicular component. p = (2p + p )/3 is the average scalar proton pressure. For simplicity, we do not have a subscript that stands for proton. p e represents the electron pressure that is assumed to be isotropic. w ± are the Alfvén wave energy densities, with the subscript + for the Alfvén waves propagating in the direction of B, and the subscript for the Alfvén waves propagating antiparallel to B. V A = B/ μ 0 ρ is the Alfvén velocity. τ ei represents the relaxation time for the heat exchanges between the electrons and the protons due to their collisions. The polytropic index γ is taken to be 5/3. We use a spatially varying electron polytropic index γ e in the electron pressure equation (6): γ e = γf S + γ H (1 f S ), (8) where 1 f S = (9) 1 + (r/r H ) 2 with r H = 5R describes the transition between the regions dominated by collisional and collisionless heat conduction, resctively, similar to the heat flux formulation of Chandran et al. (2011). In the collisionless region, r r H, the polytropic index becomes γ H = γ + α 1 + α, (10) which corresponds to the collisionless heat flux of Hollweg (1978) in the steady-state limit. We set α = 1.05 (Cranmer et al. 2009), resulting in γ H = See Apndix B for details. The collisional heat flux, q e = f S q e,s, is also weighted by the radius-dendent function f S because it is not valid in the
3 collisionless region (Chandran et al. 2011). The Spitzer & Härm (1953) collisional heat conduction is q e,s = κ e T 5/2 e b (b ) T e, (11) where κ e Wm 1 K 7/2 and T e is the electron temrature. The source term Q rad in equation (6) represents the optically thin radiative energy loss, primarily in the lower corona: Q rad = N e N h (T e ), (12) where N e and N h are the electron and hydrogen number densities, resctively, and (T e ) is the radiative cooling curve obtained from the CHIANTI version 7.1 tables (Landi et al. 2013). The reflection and dissipation of the Alfvén waves are addressed on the right-hand side of the wave energy density equation (7). The wave reflection rate R is written as R = min [R imb, max(ɣ ± )] ( ) 1 2 w w + if 4w w + 0 if (1/4)w <w + < 4w, (13) ) if 4w + w ( 2 w+ w 1 R imb = [(V A )logv A ] 2 + (b [ u]) 2, (14) where the subscript imb stands for imbalanced turbulence. The derivation of the reflection rate uses the incompressibility condition for the turbulent fluctuations. The wave dissipation rate can be expressed as Ɣ ± = 2 L w ρ, (15) where L is the transverse correlation length of the Alfvén waves and L B 1/2.WeuseL B = 150 km T 1/2 following Hollweg (1986) and Hollweg & Johnson (1988). The wave dissipation is caused by the non-linear interaction between the oppositely propagating waves, which triggers the transverse energy cascading to small scales around the proton gyroradius and below. The dissipated wave energy is separated to three parts that contribute to proton parallel temrature, proton average temrature and electron temrature, resctively. These are represented by Q p, Q p and Q e in equations (4), (5) and (6) based on kinetic Alfvén wave theory. The expressions for Q p, Q p and Q e can be found in Chandran et al. (2011) and van der Holst et al. (2014). The derivation of the wave energy equation (7), model consistency and overall energy conservation are presented in van der Holst et al. (2014). Here we focus on an important term related to temrature anisotropy. δp /δt in equation (4) is the pressure anisotropy relaxation term that limits the proton pressure anisotropy by the firehose (f), mirror (m) and ion cyclotron (ic) instability constraints in the unstable regions, with an option of globally limiting the anisotropy when necessary: δp δt ( p f p = max ( p m p min τ m τ f, p p ) or τ g, p ic p, p p τ ic τ g ). (16) p f, p m and p ic are the proton parallel pressures for the marginally unstable states, and τ f, τ m and τ ic are the relaxation times based on the instability growth rates (Apndix A). The three instabilities are the major mechanisms that reduce the pressure anisotropy in AWSoM: temrature anisotropy 3699 the corona and heliosphere. The global pressure anisotropy limiting term, when present, is applied everywhere in the computational domain to mimic the effect of some minor mechanisms that also drive the proton pressure towards isotropy. τ g is a constant parameter from input, which is normally set to be much larger than the local values of the instability-based relaxation times. This is to minimize the impact of this term in the regions with instabilities. The model described by equations (1) (7) has been implemented into the BATS-R-US MHD code(powell et al.1999; Meng et al. 2012b), and used to represent both the SC and IH components of the Space Weather Modeling Framework (SWMF; Tóth et al. 2005, 2012). For the model to be used in a rotating frame, we have also implemented additional terms due to centrifugal and Coriolis effects, i.e. we added ρ[ ( r) + 2 u] to the right-hand side of the momentum equation (2), where is the angular velocity of solar rotation. The inner boundary is at the top of the chromosphere, where the ion and electron temratures are assumed to be the same constant and uniform value. The ion temrature is isotropic at the boundary; in other words, the ion parallel temrature is the same as the ion average temrature and the electron temrature. This assumption is approximately valid due to the abundant particle collisions in the relatively dense plasma at the top of the chromosphere. 3 SIMULATION OF CR2107 The coronal simulation results for CR2107 (2011 February 16 through March 16) have been presented and analysed in van der Holst et al. (2014). In this section, we report the modelling results in the heliosphere and at 1 AU for the same CR. The simulation is computed with coupled SC and IH components of the SWMF to propagate the solar wind solution from the corona to 1 AU. Both components use the heliographic rotating frame. For the SC component, the computational domain is a sphere centred at the Sun with a radius of 24 solar radii (R S ). We use a stretched spherical grid with the radial resolution varying from about R S very close to the surface of the Sun to about 1 R S further away from the Sun. The angular cell size is about 1.4 inside 1.7 R S and 2.8 outside, in both meridional and azimuthal directions. The total number of cells is about 3 million. The inner boundary conditions are set identically as described in van der Holst et al. (2014). For the IH component, the computational domain is a square box surrounding the spherical domain of SC. The box extends from 250 to 250 R S along each of the X, Y and Z directions. We use an adaptive Cartesian grid with resolution varying from less than 0.5 R S to about 4 R S. The total number of cells is about 8 million. The global pressure anisotropy relaxation time is τ g = s in both SC and IH components. For comparisons, we also rform an isotropic MHD simulation, which is set up identically as the anisotropic MHD simulation except that the proton temrature is isotropic in the former. We will refer to the isotropic MHD simulation as 2T modelling, while the anisotropic MHD one as 3T modelling in the result analysis. Both simulations are first computed in the local time stepping mode for iterations with the SC component only. At the th step, we refine the grid near the current sheet and a shell with radial distance less than 1.2 R S to the centre of the Sun. We couple the IH component to the SC component once, and then switch off SC, and advance IH from the steady-state SC solution for 4000 iterations in the local time stepping mode. These numbers of iterations are chosen for the solutions to become approximately stationary. The results presented below are taken from the end of the simulations.
4 3700 X. Meng et al. The magnetic field information near the inner boundary of SC is obtained from the Solar Dynamic Observatory (SDO)/Helioseismic and Magnetic Imager, shown as the input magnetogram in Fig. 1. There are some active regions that can be identified by several pairs of large red and blue dots representing the strong magnetic field regions with bipolarity. 3.1 Proton temrature anisotropy and firehose instability First, we examine the spatial distribution of the proton temrature anisotropy in our model. The top panel in Fig. 2 shows the proton pressure anisotropy distribution in the Y = 0 and Z = 0 planes, from both SC and IH solutions. Overall, we find highly rndicular pressure close to the inner boundary, and the pressure anisotropy ratio p /p gradually reduces away from the Sun. This feature is more prominent in the two polar regions than close to the equatorial region as seen in the Y = 0 plane. The polar regions are dominated with on field lines, which can result in adiabatic focusing, i.e. a more rndicular pitch-angle distribution towards the stronger magnetic field region along the magnetic field lines. In the Z = 0 plane, proton pressure anisotropy exhibits corotational structure. The pressure anisotropy distribution can be related to the solar wind sed distribution in the heliosphere beyond several tens of solar radii. The bottom panel in Fig. 2 displays the flow sed from the same 3T simulation. Comparing the anisotropy and the flow sed, we find that the fast solar wind originating from high-latitude regions is characterized with highly parallel proton pressure, while the slow solar wind coming from low-latitude regions is featured with almost isotropic or slightly rndicular proton pressure. In the Z = 0 plane, the highly rndicular pressure region away from the Sun is near the interface between regions with slow and fast solar wind, i.e. the corotating interaction region (CIR). To better understand the proton pressure anisotropy variation in the simulation, we extract the proton pressures and density along the positive Z-axis, which are used to calculate the parallel and rndicular temratures along this direction. The profiles along the Z-axis provide information about the fast solar wind. The temrature and temrature anisotropy profiles are shown in Fig. 3. For the temrature profiles shown in the left plot, outside 10 RS, the parallel temrature decreases slower than the rndicular temrature with increased distance from the Sun along both axes, so that the parallel temrature exceeds the rndicular temrature at around 40 RS. The temrature profiles of the fast solar wind along the Z direction can be compared to the results obtained from the 2D MHD model with proton anisotropy (fig. 2 c in Li et al. 2004) and 1D MHD model with proton anisotropy (fig. 4 in Chandran et al. 2011). Our temrature profiles are very close to theirs in terms of the shas and magnitudes, as well as the measured values between 1 and 3 RS shown in Chandran et al. (2011). Furthermore, Chandran et al. (2011) obtain isotropic proton temrature at a radial distance of about 40 RS, which is the same as our results for CR2107. The right-hand panel of Fig. 3 displays the proton temrature anisotropy ratio T /T along the Z-axis, together with the firehose, mirror and ion cyclotron instability criteria. The firehose instability limits T /T from below, while the other two instability limits the ratio from above. The anisotropy is bounded by the ion cyclotron instability between 2 and 20 RS, and by the firehose instability beyond 200 RS. A global view of the firehose unstable regions in the simulated heliosphere is shown in Fig. 4. The coloured regions are firehose unstable, and the colour scale represents the value of the instability growth rate. The slower the instability grows, the longer τ f is, and the less the anisotropy is limited. The value of τ f also reflects the strength of the instability, or by how much T exceeds the marginal firehose unstable state. Comparing the figure with the flow sed plots in Fig. 2, we find that the firehose instability mostly exists in the simulated interface between the fast and slow solar wind in the Y = 0 plane, and its strength gradually reduces away from the Sun. 3.2 Solar wind prorties at 1 AU As the model validation, we compare the simulated solar wind prorties at 1 AU to satellite observations. Fig. 5 shows the comparison of the simulations with the WIND measurements. We extract the magnetic field strength, flow sed and number density along the WIND trajectory, which is very close to the Earth orbit shown in Fig. 2, from both 3T and 2T simulations. In this way, we can not only evaluate the 3T model, but also compare the 3T and 2T modelling results. Although WIND is at the L1 point instead of the Earth, we neglect this difference since the distance between the L1 point and the Earth is small compared to the heliospheric scale. Also, the grid resolution in our simulation is about 4 RS near 1 AU, which is even larger than the distance between the L1 point and Earth (less than 3 RS ). The comparison indicates that the 3T simulation is similar to the 2T simulation in matching the WIND observations. However, the 3T simulation seems to reproduce the observed fast solar wind sed during March 1 and March 5 better than the 2T simulation. The bottom panel in Fig. 5 compares Figure 1. The input radial magnetic field for CR2107 simulations.
5 AWSoM: temrature anisotropy 3701 Figure 3. Simulated profiles of proton temratures and temrature anisotropy along the positive Z-axis for CR2107. The dotted lines represent the instability thresholds. the simulated and observed proton pressure anisotropy. Similar to the observation, the model shows relatively small deviation from isotropy, which is consistent with the high similarity between 3T and 2T simulation results. A few aks with large p /p in the observation are missed by the model. Here we point out that the sudden jump in the pressure anisotropy during February 18 and 20, together with the simultaneous jumps in the magnetic field strength and number density, is because of an interplanetary coronal mass ejection (ICME), which is not modelled in our simulations presented here. We also compare the 3T and 2T simulations with the STEREO-A and STEREO-B data in Figs 6 and 7, resctively. As shown in the bottom-right panel of Fig. 2, the STEREO-A satellite orbits a little Figure 2. The simulated proton pressure anisotropy ratio p /p in the Y = 0 and Z = 0 planes (top panel) for CR2107. The bottom panel displays the simulated solar wind sed in the two planes, and the plot on the right also shows the trajectories of the Earth, STEREO-A and B satellites projected to the Z = 0 plane. Note that these trajectories are not closed loops in the 3D space. The asterisks mark the locations of the Earth, STEREO-A and B satellites at the beginning of CR2107.
6 3702 X. Meng et al. Figure 4. The simulated firehose unstable region coloured by the growth-rate-based relaxation time τ f in the Y = 0andZ = 0 planes from IH for CR2107. The unit for the axes is solar radii. Figure 5. The simulated solar wind prorties along the WIND orbit and the WIND data during CR2107. The ICME (not modelled) time riod is highlighted in grey. bit closer to the Sun, and the STEREO-B satellite orbits further away from the Sun, compared to the trajectory of the Earth. Similar to the comparison with the WIND data, comparisons with the STEREO-A and STEREO-B measurements show that the 3T simulation captures the fast solar wind sed slightly better than the 2T simulation does. Other than that, the 3T and 2T simulation results are very close. The 1 AU validation for CR2107 implies that the 3T model produces a slower fast solar wind than the 2T model. Since the solar wind accelerates close to its terminal sed within tens of solar
7 AWSoM: temrature anisotropy 3703 Figure 6. The simulated solar wind prorties along the STEREO-A trajectory compared with the actual data from STEREO-A for CR2107. The grey areas mark ICMEs observed by STEREO-A. Figure 7. The simulated solar wind prorties along the STEREO-B trajectory compared with the actual data from STEREO-B for CR2107. The grey areas mark ICMEs observed by STEREO-B.
8 3704 X. Meng et al. Figure 9. The input radial magnetic field for CR2123 simulations. radii, the difference in the fast solar wind sed from the 3T and 2T simulations can be seen clearly in the SC component, shown in the top panel of Fig. 8. The 3T model simulated flow in the high-latitude regions, which forms the fast solar wind, is considerably slower than the 2T model simulated one, escially beyond 10 RS. The corresponding profiles of the solar wind along the positive Z-axis, with extension to the IH component, are shown in the bottom-left plot. Close to 1 AU, the solar wind sed in the 3T simulation is about 200 km s 1 less than that in the 2T simulation. This is most likely due to the p I + (p p )bb term in the momentum equation of the 3T model, while the corresponding term in the 2T model is p I. The separation between the parallel and rndicular pressures in the 3T model allows the two pressure components to evolve differently. The bottom-middle plot of Fig. 8 shows the proton pressure anisotropy in the Y = 0 plane from the 3T simulation. Note that the colour scale is saturated, and the highest value for p /p is more than 10 near the Sun. The plot indicates p > p in the high-latitude regions within the SC domain. Thus, the protons are heated less in the direction parallel to the magnetic field than in the direction rndicular to the magnetic field. This feature Figure 8. The top panel shows the solar wind sed in the Y = 0 plane from the SC component in the 2T and 3T simulations. The bottom-left plot shows the solar wind sed profiles extracted along the Z-axis from the top plots, with extension to the IH domain. The bottom-middle plot displays the proton pressure anisotropy ratio in the Y = 0 plane from SC in the 3T simulation. The colour scale of p /p is saturated. The bottom-right plot compares the scalar pressure and parallel pressure along the Z-axis from the 2T and 3T simulations, resctively.
9 AWSoM: temrature anisotropy 3705 cannot be described by the 2T model, which assumes isotropic proton pressure, i.e. p = p. In the 3T simulation, the ratio p /p decreases with increased radial distance in the polar regions outside a few solar radii. The scial case along the positive Z-axis is shown in the bottom-right plot. Beyond about Z = 2 RS, the parallel pressure gradient in the 3T simulation is smaller than the scalar pressure gradient in the 2T simulation. The smaller pressure gradient in the 3T simulation provides less force to drive the fast solar wind propagating outwards along the on field lines, which leads to a smaller flow sed eventually. 4 S I M U L AT I O N O F C R In this section, we apply the model to simulate CR2123 from 2012 April 28 through May 25, which is during solar maximum and has more active regions with stronger magnetic fields on the solar surface compared to CR2107. The simulation is set up in an identical way as for CR2107, except with a different input magnetogram, which is taken from the National Solar Observatory s Global Oscillation Network Group, shown in Fig. 9. The simulated proton pressure anisotropy and flow sed with the 3T model are shown in Fig. 10. Similar to the results of the CR2107 simulation, highly parallel proton pressure is found in the fast solar wind region, while large p /p occurs close to the Sun Figure 11. The simulated proton pressure anisotropy along the WIND trajectory against the WIND data during CR2123. and the CIR. The slow solar wind is featured with nearly isotropic or slightly rndicular pressure. We extract the proton pressure anisotropy ratio along the Earth trajectory (shown in Fig. 10) from the 3T simulation, and compare it with the WIND satellite measurement. The comparison is shown Figure 10. The simulated proton pressure anisotropy ratio p /p and solar wind sed in the Y = 0 and Z = 0 planes for CR2123. The trajectories of the Earth, STEREO-A and B satellites are projected to the Z = 0 plane in the bottom-right plot. Again, these trajectories are not closed loops in the 3D space. The asterisks mark the locations of the Earth, STEREO-A and B satellites at the beginning of CR2123.
10 3706 X. Meng et al. Figure 12. The simulated solar wind prorties along the Earth orbit and the OMNI data during CR2123. The grey area marks an ICME occurred during this CR. Figure 13. The simulated solar wind prorties along the STEREO-A trajectory compared with the actual data from STEREO-A for CR2123. The grey areas mark ICMEs observed by STEREO-A.
11 AWSoM: temrature anisotropy 3707 Figure 14. The simulated solar wind prorties along the STEREO-B trajectory compared with the actual data from STEREO-B for CR2123. The grey areas mark ICMEs observed by STEREO-B. in Fig. 11, from which we can tell that the simulation roughly captures the observed anisotropy ratio during this CR, though with less variations. Fig. 12 displays the magnetic field magnitude, the solar wind sed, the number density and the proton average temrature along the Earth orbit from both 3T and 2T simulations, as well as the OMNI data, which contains the solar wind and initial mass function (IMF) conditions at 1 AU. Overall, both 3T and 2T simulations capture the general trends in the observed flow sed, density and temrature during most of the time. The poor rformance in reproducing the observed magnetic field is largely due to the inaccuracy of the synoptic magnetograms near the polar regions that are the main source of the magnetic flux in the solar wind. The difference between the 3T and 2T simulations is small. In particular, the fast solar wind sed from the 3T simulation seems slightly lower than the 2T simulation, which fits our discussion at the end of Section 3.2. In the last two data-model comparisons for CR2123 shown in Figs 13 and 14, we compare the solar wind prorties along the STEREO-A and STEREO-B trajectories with the actual data. For both 3T and 2T simulations, the comparisons with STEREO observations are not as good as the comparison with the OMNI data. This is because that the synoptic magnetogram is centred for the Earth, yielding the most realistic simulation results at the Earth. 5 CONCLUSION We have successfully develod a 3D 3T Alfvén wave-driven solar wind model AWSoM. The model assumes anisotropic proton pressure and isotropic electron pressure. The Alfvén wave energy dissipates into three parts: proton parallel, proton rndicular and electron temratures. At the top of the chromosphere where the inner boundary of the model lies, we assume isotropic proton temrature. To limit the proton pressure anisotropy in the model, we have introduced the pressure anisotropy relaxation term in the parallel proton pressure equation, which constrains the anisotropy based on the growth rates of the firehose, mirror and proton cyclotron instabilities, as well as an optional global anisotropy relaxation term. The model is applied to simulate the steady states of two CRs during the solar cycle 24. In the simulations of both CR2107 and CR2123, we find highly rndicular proton pressure close to the Sun, and the anisotropy ratio p /p reduces away from the Sun. There are also corotational structures in the proton pressure anisotropy distributions. In the heliosphere beyond about R S, the fast solar wind is characterized with p /p less than 1, while the slow solar wind is with p /p close to or slightly over 1. Near the CIRs, p /p is typically the highest and larger than 1. We have also obtained reasonable agreement in the comparisons with in situ observational data for the two CRs. At 1 AU, the difference between the 3T and 2T simulations is very small. Interestingly, the 3T simulations predict slightly smaller and better fast solar wind sed than the 2T simulations do, which is primarily due to smaller pressure gradient along the magnetic field to accelerate the fast solar wind in the 3T simulations. Note that we have validated the IMF magnitude instead of its individual components at 1 AU, because the current generation of predictive solar wind models, AWSoM included, are incapable of providing reliable IMF components due to the missing physics in the models and the limited accuracy of magnetograms. In the future when AWSoM is sufficiently advanced, we could validate IMF components, escially B z, which is the most relevant to space weather predictions.
12 3708 X. Meng et al. This is the first time that the proton temrature anisotropy is included in a 3D solar wind model. The model does not offer much improvement over the isotropic pressure 2T AWSoM for the solar wind at 1 AU, unlike the notable improvement over the 2T model for the SC presented in van der Holst et al. (2014). This is excted, given that the temrature anisotropy is more prominent near the Sun than further away. In the heliosphere where the temrature anisotropy is small or even negligible, the 3T model naturally reduces to the 2T model. Nevertheless, our model extends the current capability of 3D solar wind modelling. As next steps, the model can be further improved by incorporating some missing physics. One factor, which also impacts the solar wind acceleration, is that the momentum equation needs to be modified by Alfvénic turbulence with proton temrature anisotropy. More scifically, the Alfvén wave pressure gradient term [the last term on the left of equation (2)] will contain the proton temrature anisotropy after the modification. This will be addressed in the future. ACKNOWLEDGEMENTS This work was supported by the NSF grant AGS We would like to acknowledge high-rformance computing support from Yellowstone (ark:/85065/d7wd3xhc) provided by NCAR s Computational and Information Systems Laboratory, sponsored by the NSF. We also thank J. H. King, N. Papatashvilli, K. Ogilvie at AdnetSystems, NASA GSFC and CDAWeb for providing the OMNI, WIND, STEREO-A and B data. The near- Earth ICMEs were identified based on the WIND ICME list ( and the Richardson/Cane ICMElist. The STEREO-observed ICMEs were obtained from Dr Lan Jian s ICME list. REFERENCES Antonucci E., Dodero M. A., Giordano S., 2000, Sol. Phys., 197, 115 Barnes A., 1966, Phys. Fluids, 9, 8 Chandran B. D. G., Dennis T. J., Quataert E., Bale S. D., 2011, ApJ, 743, 197 Chandrasekhar S., Kaufman A. N., Watson K. M., 1958, Proc. R. Soc. Lond. Ser. A, 245, 435 Chew G. F., Goldberger M. L., Low F. E., 1956, Proc. R. Soc. Lond. Ser. A, 236, 112 Cranmer S. R., Matthaeus W. H., Breech B. A., Kasr J. C., 2009, ApJ, 702, 1604 Feldman W. C., Asbridge J. R., Bame S. J., Montgomery M. D., 1974, Rev. Geophys. Space Phys., 12, 4 Feng X., Zhou Y., Wu S. T., 2007, ApJ, 655, 2 Gary S. P., Montgomery M. D., Feldman W. C., Forslund D. W., 1976, J. Geophys. Res., 81, 1241 Gary S. P., Skoug R. M., Steinberg J. T., 2001, Geophys. Res. Lett., 28, 2759 Gary S. P., Goldstein B. E., Neugebauer M., 2002, J. Geophys. Res., 107, 1169 Gressl C., Veronig A. M., Temmer M., Odstrc Zil D., Linker J. A., MikićZ., Riley P., 2014, Sol. Phys., 289, 5 Hellinger P., Trávníček O., Kasr J. C., Lazarus A. J., 2006, Geophys. Res. Lett., 33, L09101 Hollweg J. V., 1978, Rev. Geophys. Space Phys., 16, 689 Hollweg J. V., 1986, J. Geophys. Res., 91, 4111 Hollweg J. V., Johnson W., 1988, J. Geophys. Res., 93, 9547 Hu Y. Q., Esser R., Habbal S. R., 1997, J. Geophys. Res., 102, Hundhausen A. J., 1968, Space Sci. Rev., 8, 690 Kasr J. C., Lazarus A. J., Gary S. P., Szabo A., 2003, in Velli M., Bruno R., Malara F., eds, AIP Conf. Proc. Vol. 679, Proc. 10th Inter. Solar Wind Conf., Solar Wind Ten. Am. Inst. Phys., New York, p. 538 Kennel C. F., Petschek H. E., 1966, J. Geophys. Res., 71, 1 Kleimann J., Kopp A., Fichtner H., Grauer R., 2009, Ann. Geophys., 27, 989 Kohl J. L. et al., 1998, ApJ, 501, L127 Kohl J. L., Noci G., Cranmer S. R., Raymond J. C., 2006, A&AS, 13, 31 Landi E., Young P. R., Dere K. P., Del Zanna G., Mason H. E., 2013, ApJ, 763, 86 Leer E., Axford W. I., 1972, Sol. Phys., 23, 238 Li X., 1999, J. Geophys. Res., 104, Li X., Habbal S. R., Kohl J. L., Noci G., 1998, ApJ, 501, L133 Li B., Li X., Hu Y.-Q., Habbal S. R., 2004, J. Geophys. Res., 109, A07103 Marsch E., Mühlhäuser K.-H., Schwenn R., Rosenbauer H., Pilipp W., Neubauer F. M., 1982, J. Geophys. Res., 87, 52 Marsch E., Ao X.-Z., Tu C.-Y., 2004, J. Geophys. Res., 109, A04102 Meng X., Tóth G., Liemohn M. W., Gombosi T. I., Runov A., 2012a, J. Geophys. Res., 117, A08216 Meng X., Tóth G., Sokolov I. V., Gombosi T. I., 2012b, J. Comput. Phys., 231, 3610 Powell K., Roe P., Linde T., Gombosi T., De Zeeuw D. L., 1999, J. Geophys. Res., 109, A9 Riley P., Linker J. A., Mikić Z., 2001, J. Geophys. Res., 106, A8 Spitzer L., Härm R., 1953, Phys. Rev., 89, 977 Telloni D., Antonucci E. M., Dodero M. A., 2007, A&A, 476, 1341 Tóth G. et al., 2005, J. Geophys. Res., 110, A12226 Tóth G. et al., 2012, J. Comput. Phys., 231, 870 van der Holst B., Poedts S., Chané E., Jacobs C., Dubey G., Kim D., 2005, Space Sci. Rev., 121, 91 van der Holst B., Meng X., Jin M., Manchester W. B., Tóth G., Gombosi T. I., 2014, ApJ, 782, 81 Watermann J., Vainio R., Lilensten J., Belehaki A., Messerotti M., 2009, Space Sci. Rev., 147, 111 Whang Y. C., 1972, ApJ, 178, 221 APPENDIX A: FIREHOSE, MIRROR AND ION CYCLOTRON INSTABILITIES The firehose, mirror and ion cyclotron instabilities are three tys of instabilities resulted from, and thus constrain, proton pressure anisotropy in magnetized plasmas (Chandrasekhar, Kaufman & Watson 1958; Barnes1966; Kennel & Petschek 1966; Gary et al. 1976). These instabilities, escially the firehose instability, have been frequently observed in the solar wind. Below we briefly introduce each instability and list their criteria, p at marginally unstable states, as well as the relaxation times applied in our model. More detail can be found in Meng et al. (2012a). The firehose instability arises when the parallel pressure is sufficiently large that a small rturbation rndicular to the magnetic field grows. This can be viewed as an analogue of the violent motion of a firehose with water flowing through rapidly. The instability criterion is p > 1 + B2. p μ 0 p (A1) Hence, the parallel pressure in marginally firehose unstable plasmas is given by p f = p + 2B2. (A2) 3μ 0 The maximum growth rate of the firehose instability γ f, m determines τ f in equation (16): τ f = 1 = 2 p (p p /4), γ f,m i p p B 2 /μ 0 where i is proton gyrofrequency. (A3)
13 The mirror and ion cyclotron instabilities occur when the rndicular pressure is large enough. The mirror instability is excited when part of the particles are trapd in small magnetic mirrors formed by the rturbed magnetic field. A typical observational feature is the anti-phased variations in the magnetic field strength and the plasma number density. The criterion is p p > 1 + B2 2μ 0 p. When marginally unstable, ( p m = 1 B2 B 2 ) 2 + 6p + 12B2 p + 9p 3 μ 0 μ 2. 0 μ 0 (A4) (A5) The mirror instability involves kinetic effects, and thus cannot be fully described by anisotropic MHD. The MHD analysis may still be employed to find τ m, because the wavelength corresponding to the maximum growth rate γ m,m is much larger than the ion gyroradius. τ m is approximated as τ m = 1 = 3 5 γ m,m 4 i p 2(p p ) B 2 p /(2μ 0 p ). (A6) The ion cyclotron instability is a kinetic and resonant instability that occurs when the electric field vector of the ion cyclotron wave rotates synchronously with the particle gyration. Its criterion can be written as p p > 1 + C 1 ( B 2 2μ 0 p ) C2, (A7) where C1 andc2are constants that vary with different maximum growth rates. In our model, we use C 1 = 0.3 and C 2 = 0.5. The marginally unstable p then becomes p ic = 1 B B2 + 2p. (A8) 2 μ 0 μ 0 The analytical expression for the maximum growth rate of the ion cyclotron instability can only be obtained through kinetic analysis and is very complicated. We use the following approximate expression for τ ic : τ ic = 1 γ ic,m = 102 i. (A9) APPENDIX B: ELECTRON PRESSURE EQUATION In this apndix, we show how the modified polytropic index defined in equation (8) can be chosen so that it reproduces both the collisional heat conduction (Spitzer & Härm 1953) near the Sun and the collisionless heat flux formulation of Hollweg (1978) at large distances from the Sun. The latter limit is strictly true for steady-state conditions. We start from equations (6) (9): ( ) ( ) + t γ e 1 γ e 1 u + p e ( u) where = q e + p p e τ ei Q rad + Q e (B1) γ e = γf S + γ H (1 f S ) (B2) AWSoM: temrature anisotropy f S = 1 + (r/r H ). (B3) 2 Note that this form of the electron pressure equation describes the time evolution of the electron internal energy p e /(γ e 1) for an arbitrary function of the polytropic index γ e. We choose f S and γ e so that the above electron pressure equation reproduces the effects of the collisional and collisionless heat conductions in the appropriate limits. In the collisional limit r r H,sothatf S = 1, γ e = γ = 5/3 and q e = q e, S is the Spitzer heat flux (equation 11), which results in the usual form of the electron pressure equation with collisional heat conduction. In the collisionless limit r r H,sothatf S = 0, γ = γ H and q e = 0 and we get t ( γ H 1 ) ( ) + γ H 1 u + p e ( u) = p p e Q rad + Q e. (B4) τ ei Taking the steady-state limit and dropping the electron ion heat exchange and radiative cooling terms that are negligible and not essential in the following derivation, the equation reduces to ( ) γ H 1 u + p e ( u) = Q e. (B5) This can be rewritten as ( ) γ 1 u + p e ( u) [ ( 1 = p e u γ H 1 1 )] + Q e. (B6) γ 1 Substituting γ H = (γ + α)/(1 + α) from equation (10) gives ( ) [ ] γ 1 u 1 + p e u = γ 1 αp eu + Q e. (B7) The first term on the right-hand side is the divergence of the collisionless heat flux αp e u of Hollweg (1978). Note that the electron pressure equation with a modified polytropic index is indendent of the frame of reference; therefore, the steady-state solution in the corotating frame is also a valid (timedendent) solution in the inertial frame. As long as the solution is approximately steady state in the corotating frame, the modified polytropic index γ H gives the same result as Hollweg s collisionless heat conduction. This relationship was also pointed out by Hollweg (private communication). This par has been tyset from a TEX/LATEX file prepared by the author.
SOLAR WIND ION AND ELECTRON DISTRIBUTION FUNCTIONS AND THE TRANSITION FROM FLUID TO KINETIC BEHAVIOR
SOLAR WIND ION AND ELECTRON DISTRIBUTION FUNCTIONS AND THE TRANSITION FROM FLUID TO KINETIC BEHAVIOR JUSTIN C. KASPER HARVARD-SMITHSONIAN CENTER FOR ASTROPHYSICS GYPW01, Isaac Newton Institute, July 2010
More informationSolar Wind Turbulent Heating by Interstellar Pickup Protons: 2-Component Model
Solar Wind Turbulent Heating by Interstellar Pickup Protons: 2-Component Model Philip A. Isenberg a, Sean Oughton b, Charles W. Smith a and William H. Matthaeus c a Inst. for Study of Earth, Oceans and
More informationProtons and alpha particles in the expanding solar wind: Hybrid simulations
JOURNAL OF GEOPHYSICAL RESEARCH: SPACE PHYSICS, VOL. 118, 5421 5430, doi:10.1002/jgra.50540, 2013 Protons and alpha particles in the expanding solar wind: Hybrid simulations Petr Hellinger 1,2 and Pavel
More informationHeating of ions by low-frequency Alfven waves
PHYSICS OF PLASMAS 14, 433 7 Heating of ions by low-frequency Alfven waves Quanming Lu School of Earth and Space Sciences, University of Science and Technology of China, Hefei 36, People s Republic of
More informationMagnetic Effects Change Our View of the Heliosheath
Magnetic Effects Change Our View of the Heliosheath M. Opher Λ, P. C. Liewer Λ, M. Velli, T. I. Gombosi ΛΛ, W.Manchester ΛΛ,D. L. DeZeeuw ΛΛ,G.Toth ΛΛ and I. Sokolov ΛΛ Λ Jet Propulsion Laboratory, MS
More informationHybrid simulation of ion cyclotron resonance in the solar wind: Evolution of velocity distribution functions
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110,, doi:10.1029/2005ja011030, 2005 Hybrid simulation of ion cyclotron resonance in the solar wind: Evolution of velocity distribution functions Xing Li Institute
More informationarxiv: v1 [physics.space-ph] 13 Dec 2018
On the generation mechanism of electromagnetic cyclotron waves in the solar wind: statistical results from Wind observations arxiv:1812.05323v1 [physics.space-ph] 13 Dec 2018 G. Q. Zhao 1,2, H. Q. Feng
More informationLecture 5 The Formation and Evolution of CIRS
Lecture 5 The Formation and Evolution of CIRS Fast and Slow Solar Wind Fast solar wind (>600 km/s) is known to come from large coronal holes which have open magnetic field structure. The origin of slow
More informationTHEORETICAL RESEARCH ON SOLAR WIND TURBULENCE
THEORETICAL RESEARCH ON SOLAR WIND TURBULENCE A White Paper Submitted to the NRC Decadal Survey of Solar and Space Physics Benjamin D. G. Chandran, Eliot Quataert, Jean Perez, Aveek Sarkar, Steve Cranmer,
More informationTurbulent Origins of the Sun s Hot Corona and the Solar Wind
Turbulent Origins of the Sun s Hot Corona and the Solar Wind Steven R. Cranmer Harvard-Smithsonian Center for Astrophysics Turbulent Origins of the Sun s Hot Corona and the Solar Wind Outline: 1. Solar
More informationOn MHD Modeling of Coronal Mass Ejections
NUMERICAL MODELING OF SPACE PLASMA FLOWS / ASTRONUM-2007 ASP Conference Series, Vol. 385, c 2008 Nikolai V. Pogorelov, Edouard Audit, and Gary P. Zank, eds. On MHD Modeling of Coronal Mass Ejections Horst
More informationO 5+ at a heliocentric distance of about 2.5 R.
EFFECT OF THE LINE-OF-SIGHT INTEGRATION ON THE PROFILES OF CORONAL LINES N.-E. Raouafi and S. K. Solanki Max-Planck-Institut für Aeronomie, 37191 Katlenburg-Lindau, Germany E-mail: Raouafi@linmpi.mpg.de;
More informationMHD MODELING FOR HMI JON A. LINKER SCIENCE APPLICATIONS INTL. CORP. SAN DIEGO
MHD MODELING FOR HMI ZORAN MIKIĆ JON A. LINKER SCIENCE APPLICATIONS INTL. CORP. SAN DIEGO Presented at the HMI Team Meeting Stanford University, Palo Alto, May 1 2, 23 USEFULNESS OF MHD MODELS A global
More informationarxiv: v2 [astro-ph.sr] 3 Aug 2010
Correlations between the proton temperature anisotropy and transverse high-frequency waves in the solar wind Sofiane Bourouaine 1, Eckart Marsch 1 and Fritz M. Neubauer 2 arxiv:1003.2299v2 [astro-ph.sr]
More informationA NEW MODEL FOR REALISTIC 3-D SIMULATIONS OF SOLAR ENERGETIC PARTICLE EVENTS
A NEW MODEL FOR REALISTIC 3-D SIMULATIONS OF SOLAR ENERGETIC PARTICLE EVENTS Nicolas Wijsen KU Leuven In collaboration with: A. Aran (University of Barcelona) S. Poedts (KU Leuven) J. Pomoell (University
More informationComparison of the effects of wave-particle interactions and the kinetic suprathermal electron population on the acceleration of the solar wind
A&A 395, 1001 1012 (2002) DOI: 10.1051/0004-6361:20021348 c ESO 2002 Astronomy & Astrophysics Comparison of the effects of wave-particle interactions and the kinetic suprathermal electron population on
More informationMagnetic Drivers of CME Defection in the Low Corona
Magnetic Drivers of CME Defection in the Low Corona C. Kay (Boston University) M. Opher (Boston University) R. M. Evans (NASA GSFC/ORAU T. I. Gombosi (University of Michigan) B. van der Holst (University
More informationOpen magnetic structures - Coronal holes and fast solar wind
Open magnetic structures - Coronal holes and fast solar wind The solar corona over the solar cycle Coronal and interplanetary temperatures Coronal holes and fast solar wind Origin of solar wind in magnetic
More informationInverse and normal coronal mass ejections: evolution up to 1 AU. E. Chané, B. Van der Holst, C. Jacobs, S. Poedts, and D.
A&A 447, 727 733 (2006) DOI: 10.1051/0004-6361:20053802 c ESO 2006 Astronomy & Astrophysics Inverse and normal coronal mass ejections: evolution up to 1 AU E. Chané, B. Van der Holst, C. Jacobs, S. Poedts,
More informationLab #2: Activity 5 Exploring the Structure of the Solar Magnetic Field Using the MAS Model
Lab #2: Activity 5 Exploring the Structure of the Solar Magnetic Field Using the MAS Model In this lab activity we will use results from the MAS (Magnetohydrodynamics Around a Sphere) model of the solar
More informationParallel and oblique proton fire hose instabilities in the presence of alpha/proton drift: Hybrid simulations
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 111,, doi:10.1029/2005ja011318, 2006 Parallel and oblique proton fire hose instabilities in the presence of alpha/proton drift: Hybrid simulations Petr Hellinger and
More informationFluid modeling of anisotropic heating and microinstabilities
Fluid modeling of anisotropic heating and microinstabilities in space plasmas Thierry Passot UNS, CNRS, Observatoire de la Côte d Azur, Nice, France Collaborators: D. Laveder, L. Marradi, and P.L. Sulem
More informationHeliophysics Shocks. Merav Opher, George Mason University,
Heliophysics Shocks QuickTime and a TIFF (Uncompressed) decompressor are needed to see this picture. Merav Opher, George Mason University, mopher@gmu.edu Heliophysics Summer School, July 25, 2008 Outline
More informationHeating of Test Particles in Numerical Simulations of MHD Turbulence and the Solar Wind
Heating of Test Particles in Numerical Simulations of MHD Turbulence and the Solar Wind Ian Parrish UC Berkeley Collaborators: Rémi Lehe (ENS), Eliot Quataert (UCB) Einstein Fellows Symposium October 27,
More informationNONLINEAR MHD WAVES THE INTERESTING INFLUENCE OF FIREHOSE AND MIRROR IN ASTROPHYSICAL PLASMAS. Jono Squire (Caltech) UCLA April 2017
NONLINEAR MHD WAVES THE INTERESTING INFLUENCE OF FIREHOSE AND MIRROR IN ASTROPHYSICAL PLASMAS Jono Squire (Caltech) UCLA April 2017 Along with: E. Quataert, A. Schekochihin, M. Kunz, S. Bale, C. Chen,
More informationCoronal Heating versus Solar Wind Acceleration
SOHO 15: Coronal Heating, 6 9 September 2004, University of St. Andrews, Scotland Coronal Heating versus Solar Wind Acceleration Steven R. Cranmer Harvard-Smithsonian Center for Astrophysics, Cambridge,
More informationAlfvén wave heating of heavy ions in the expanding solar wind: Hybrid simulations
JOURNAL OF GEOPHYSICAL RESEARCH, VOL.???, XXXX, DOI:10.1029/, Alfvén wave heating of heavy ions in the expanding solar wind: Hybrid simulations Petr Hellinger, 1 Marco Velli, 2,4 Pavel Trávníček, 1 S.
More informationThe plasma microstate of the solar wind
The plasma microstate of the solar wind Radial gradients of kinetic temperatures Velocity distribution functions Ion composition and suprathermal electrons Coulomb collisions in the solar wind Waves and
More informationThe Physics of Collisionless Accretion Flows. Eliot Quataert (UC Berkeley)
The Physics of Collisionless Accretion Flows Eliot Quataert (UC Berkeley) Accretion Disks: Physical Picture Simple Consequences of Mass, Momentum, & Energy Conservation Matter Inspirals on Approximately
More informationSelf-Consistent Ion Cyclotron Anisotropy-Beta Relation for Solar Wind Protons
1 Self-Consistent Ion Cyclotron Anisotropy-Beta Relation for Solar Wind Protons Philip A. Isenberg 1, Bennett A. Maruca 2, and Justin C. Kasper 3 1 Institute for the Study of Earth, Oceans and Space, University
More informationMagnetic Reconnection in ICME Sheath
WDS'11 Proceedings of Contributed Papers, Part II, 14 18, 2011. ISBN 978-80-7378-185-9 MATFYZPRESS Magnetic Reconnection in ICME Sheath J. Enzl, L. Prech, K. Grygorov, A. Lynnyk Charles University, Faculty
More informationObservations of an interplanetary slow shock associated with magnetic cloud boundary layer
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 33, L15107, doi:10.1029/2006gl026419, 2006 Observations of an interplanetary slow shock associated with magnetic cloud boundary layer P. B.
More informationSTAGNATED OUTFLOW OF O +5 IONS IN THE SOURCE REGION OF THE SLOW SOLAR WIND AT SOLAR MINIMUM
The Astrophysical Journal, 602:415 421, 2004 February 10 # 2004. The American Astronomical Society. All rights reserved. Printed in U.S.A. STAGNATED OUTFLOW OF O +5 IONS IN THE SOURCE REGION OF THE SLOW
More informationAlfvén wave heating of heavy ions in the expanding solar wind: Hybrid simulations
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110, NO. A12, A12109, DOI:129/2005JA011244, 2005 Alfvén wave heating of heavy ions in the expanding solar wind: Hybrid simulations Petr Hellinger, 1 Marco Velli, 2,4
More informationMesoscale Variations in the Heliospheric Magnetic Field and their Consequences in the Outer Heliosphere
Mesoscale Variations in the Heliospheric Magnetic Field and their Consequences in the Outer Heliosphere L. A. Fisk Department of Atmospheric, Oceanic, and Space Sciences, University of Michigan, Ann Arbor,
More informationAlfvénic Turbulence in the Fast Solar Wind: from cradle to grave
Alfvénic Turbulence in the Fast Solar Wind: from cradle to grave, A. A. van Ballegooijen, and the UVCS/SOHO Team Harvard-Smithsonian Center for Astrophysics Alfvénic Turbulence in the Fast Solar Wind:
More informationSolar Wind Turbulence
Solar Wind Turbulence Presentation to the Solar and Heliospheric Survey Panel W H Matthaeus Bartol Research Institute, University of Delaware 2 June 2001 Overview Context and SH Themes Scientific status
More informationThe new Heliospheric Magnetic Field: Observational Implications
The new Heliospheric Magnetic Field: Observational Implications T. H. Zurbuchen, L. A. Fisk, S. Hefti and N. A. Schwa&on Department of Atmospheric, Oceanic and Space Sciences, University of Michigan, Ann
More informationThe MRI in a Collisionless Plasma
The MRI in a Collisionless Plasma Eliot Quataert (UC Berkeley) Collaborators: Prateek Sharma, Greg Hammett, Jim Stone Modes of Accretion thin disk: energy radiated away (relevant to star & planet formation,
More informationKinetic Alfvén waves in space plasmas
Kinetic Alfvén waves in space plasmas Yuriy Voitenko Belgian Institute for Space Aeronomy, Brussels, Belgium Solar-Terrestrial Center of Excellence, Space Pole, Belgium Recent results obtained in collaboration
More informationTurbulent heating and acceleration of He ++ ions by spectra of Alfvén-cyclotron waves in the expanding solar wind: 1.5-D hybrid simulations
JOURNAL OF GEOPHYSICAL RESEARCH: SPACE PHYSICS, VOL. 118, 2842 2853, doi:10.1002/jgra.50363, 2013 Turbulent heating and acceleration of He ++ ions by spectra of Alfvén-cyclotron waves in the expanding
More informationHeating the magnetically open ambient background corona of the Sun by Alfvén waves
A&A 411, L481 L485 (23) DOI: 1.151/4-6361:231587 c ESO 23 Astronomy & Astrophysics Heating the magnetically open ambient background corona of the Sun by Alfvén waves H. Peter 1 and C. Vocks 2 1 Kiepenheuer-Institut
More informationOn the effect of the initial magnetic polarity and of the background wind on the evolution of CME shocks
A&A 432, 331 339 (2005) DOI: 10.1051/0004-6361:20042005 c ESO 2005 Astronomy & Astrophysics On the effect of the initial magnetic polarity and of the background wind on the evolution of CME shocks E. Chané,
More informationMHD Simulation of Solar Chromospheric Evaporation Jets in the Oblique Coronal Magnetic Field
MHD Simulation of Solar Chromospheric Evaporation Jets in the Oblique Coronal Magnetic Field Y. Matsui, T. Yokoyama, H. Hotta and T. Saito Department of Earth and Planetary Science, University of Tokyo,
More informationEmpirical Testing of Solar Coronal and Solar Wind Models
Empirical Testing of Solar Coronal and Solar Wind Models Lauren Woolsey University of Maryland - College Park (2011) Mentor: Dr. Leonard Strachan Introduction What is the Solar Wind? * Outflow of particles
More informationSpace Physics. An Introduction to Plasmas and Particles in the Heliosphere and Magnetospheres. May-Britt Kallenrode. Springer
May-Britt Kallenrode Space Physics An Introduction to Plasmas and Particles in the Heliosphere and Magnetospheres With 170 Figures, 9 Tables, Numerous Exercises and Problems Springer Contents 1. Introduction
More informationMacroscopic plasma description
Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion
More informationThe Magnetic Field at the Inner Boundary of the Heliosphere Around Solar Minimum
The Magnetic Field at the Inner Boundary of the Heliosphere Around Solar Minimum X. P. Zhao and J. T. Hoeksema W. W. Hansen Experimental Physics Laboratory, Stanford University, Stanford, CA 94305-4085
More informationEnergy Analysis During the Collision of Two Successive CMEs
Numerical Modeling of Space Plasma Flows: ASTRONUM-2013 ASP Conference Series, Vol. 488 N.V.Pogorelov, E.Audit,and G.P.Zank,eds. c 2014 Astronomical Society of the Pacific Energy Analysis During the Collision
More informationPickup Proton Instabilities and Scattering in the Distant Solar Wind and the Outer Heliosheath: Hybrid Simulations
Pickup Proton Instabilities and Scattering in the Distant Solar Wind and the Outer Heliosheath: Hybrid Simulations Kaijun Liu 1,2, Eberhard Möbius 2,3, S. P. Gary 2,4, Dan Winske 2 1 Auburn University,
More informationAmbient solar wind s effect on ICME transit times
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L15105, doi:10.1029/2008gl034493, 2008 Ambient solar wind s effect on ICME transit times A. W. Case, 1 H. E. Spence, 1 M. J. Owens, 1
More informationThe solar wind. or rather: The energy budget of the corona-solar wind system. Øystein Lie-Svendsen
Introduction Facts The solar wind or rather: The energy budget of the corona-solar wind system Øystein Lie-Svendsen Norwegian Defence Research Establishment and Institute of theoretical astrophysics, University
More informationKinetic Effects in Coronal Holes & High-Speed Streams: A Roundup of Observational Constraints
Kinetic Effects in Coronal Holes & High-Speed Streams: A Roundup of Observational Constraints Steven R. Cranmer steven.cranmer@colorado.edu University of Colorado Boulder http://lasp.colorado.edu/~cranmer/
More informationPlasma properties at the Voyager 1 crossing of the heliopause
Journal of Physics: Conference Series PAPER Plasma properties at the Voyager 1 crossing of the heliopause Recent citations - Reconnection at the Heliopause: Predictions for Voyager 2 S. A. Fuselier and
More informationHow did the solar wind structure change around the solar maximum? From interplanetary scintillation observation
Annales Geophysicae (2003) 21: 1257 1261 c European Geosciences Union 2003 Annales Geophysicae How did the solar wind structure change around the solar maximum? From interplanetary scintillation observation
More informationMHD simulation of solar wind using solar photospheric magnetic field data
6-16P, LWS workshop 2004 March, Boulder MHD simulation of solar wind using solar photospheric magnetic field data Keiji Hayashi (Stanford University) keiji@quake.stanford.edu Introduction Time-dependent
More informationPlasma and Magnetic Field Observations of Stream Interaction Regions near 1 AU
Plasma and Magnetic Field Observations of Stream Interaction Regions near 1 AU Lan K. Jian 1, C.T. Russell 1, J.G. Luhmann 2, A.B. Gavin 3, D. Odstrcil 4, P.J. MacNeice 5 1 Inst. of Geophysics & Planetary
More informationExtended Coronal Heating and Solar Wind Acceleration over the Solar Cycle
Extended Coronal Heating and Solar Wind Acceleration over the Solar Cycle S. R. Cranmer, J. L. Kohl, M. P. Miralles, & A. A. van Ballegooijen Harvard-Smithsonian Center for Astrophysics Extended Coronal
More informationDynamics of Polar Jets from the Chromosphere to the Corona: Mass, Momentum and Energy Transfer
Dynamics of Polar Jets from the Chromosphere to the Corona: Mass, Momentum and Energy Transfer Judit Szente 1, Gabor Toth 1, Ward Manchester 1, Bart van der Holst 1, Enrico Landi 1, C Richard DeVore 2,
More informationA new mechanism to account for acceleration of the solar wind
A new mechanism to account for acceleration of the solar wind Henry D. May Email: hankmay@earthlink.net Abstract An enormous amount of effort has been expended over the past sixty years in attempts to
More informationWaves & Turbulence in the Solar Wind: Disputed Origins & Predictions for PSP
Waves & Turbulence in the Solar Wind: Disputed Origins & Predictions for PSP Steven R. Cranmer University of Colorado Boulder, LASP A. Schiff, S. Van Kooten, C. Gilbert, L. N. Woolsey, A. A. van Ballegooijen,
More informationA PARALLEL ADAPTIVE 3D MHD SCHEME FOR MODELING CORONAL AND SOLAR WIND PLASMA FLOWS
A PARALLEL ADAPTIVE 3D MHD SCHEME FOR MODELING CORONAL AND SOLAR WIND PLASMA FLOWS C. P. T. GROTH, D. L. DE ZEEUW and T. I. GOMBOSI Space Physics Research Laboratory, Department of Atmospheric, Oceanic
More informationThe Effects of Anisotropic Transport on Dilute Astrophysical Plasmas Eliot Quataert (UC Berkeley)
The Effects of Anisotropic Transport on Dilute Astrophysical Plasmas Eliot Quataert (UC Berkeley) in collaboration with Ian Parrish, Prateek Sharma, Jim Stone, Greg Hammett Hydra A w/ Chandra Galactic
More informationMAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT
MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT ABSTRACT A. G. Tarditi and J. V. Shebalin Advanced Space Propulsion Laboratory NASA Johnson Space Center Houston, TX
More informationSolar Energetic Particles in the Inner Heliosphere
Author: Mariona Adillón Corbera Advisor: Neus Agueda Costafreda Facultat de Física, Universitat de Barcelona, Diagonal 645, 08028 Barcelona, Spain. Abstract: The upcoming missions Solar Orbiter (SolO)
More informationTHE INTERACTION OF TURBULENCE WITH THE HELIOSPHERIC SHOCK
THE INTERACTION OF TURBULENCE WITH THE HELIOSPHERIC SHOCK G.P. Zank, I. Kryukov, N. Pogorelov, S. Borovikov, Dastgeer Shaikh, and X. Ao CSPAR, University of Alabama in Huntsville Heliospheric observations
More informationIntroduction to Plasma Physics
Introduction to Plasma Physics Hartmut Zohm Max-Planck-Institut für Plasmaphysik 85748 Garching DPG Advanced Physics School The Physics of ITER Bad Honnef, 22.09.2014 A simplistic view on a Fusion Power
More informationRadial decay law for large-scale velocity and magnetic field fluctuations in the solar wind
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 111,, doi:10.1029/2005ja011528, 2006 Radial decay law for large-scale velocity and magnetic field fluctuations in the solar wind J. J. Podesta 1 Received 10 November
More informationB.V. Gudiksen. 1. Introduction. Mem. S.A.It. Vol. 75, 282 c SAIt 2007 Memorie della
Mem. S.A.It. Vol. 75, 282 c SAIt 2007 Memorie della À Ø Ò Ø ËÓÐ Ö ÓÖÓÒ B.V. Gudiksen Institute of Theoretical Astrophysics, University of Oslo, Norway e-mail:boris@astro.uio.no Abstract. The heating mechanism
More informationInferring the Structure of the Solar Corona and Inner Heliosphere during the Maunder Minimum using MHD simulations
Inferring the Structure of the Solar Corona and Inner Heliosphere during the Maunder Minimum using MHD simulations Pete Riley, Roberto Lionello, Jon Linker, and Zoran Mikic Predictive Science, Inc. (PSI),
More informationSingle particle motion and trapped particles
Single particle motion and trapped particles Gyromotion of ions and electrons Drifts in electric fields Inhomogeneous magnetic fields Magnetic and general drift motions Trapped magnetospheric particles
More informationJet Stability: A computational survey
Jet Stability Galway 2008-1 Jet Stability: A computational survey Rony Keppens Centre for Plasma-Astrophysics, K.U.Leuven (Belgium) & FOM-Institute for Plasma Physics Rijnhuizen & Astronomical Institute,
More informationcos 6 λ m sin 2 λ m Mirror Point latitude Equatorial Pitch Angle Figure 5.1: Mirror point latitude as function of equatorial pitch angle.
Chapter 5 The Inner Magnetosphere 5.1 Trapped Particles The motion of trapped particles in the inner magnetosphere is a combination of gyro motion, bounce motion, and gradient and curvature drifts. In
More informationBeta-dependent upper bound on ion temperature anisotropy in a laboratory plasma
PHYSICS OF PLASMAS VOLUME 7, NUMBER 3 MARCH 2000 LETTERS The purpose of this Letters section is to provide rapid dissemination of important new results in the fields regularly covered by Physics of Plasmas.
More informationTurbulent dissipation in the solar wind and corona
Turbulent dissipation in the solar wind and corona W. H. Matthaeus, P. Dmitruk, S. Oughton and D. Mullan Bartol Research Institute, University of Delaware, Newark, DE 19716 USA Department of Mathematics,
More informationIon Cyclotron Resonance Frequency Interval Dependence on the O VI Ion Number Density in the North Polar Coronal Hole 1.
J. Astrophys. Astr. (2013) 34, 349 355 c Indian Academy of Sciences Ion Cyclotron Resonance Frequency Interval Dependence on the O VI Ion Number Density in the North Polar Coronal Hole 1.5R 3R Region Özgür
More informationhermally driven winds Viggo H. Hansteen Insititute of Theoretical Astrophysics, University of Oslo
hermally driven winds Viggo H. Hansteen Insititute of Theoretical Astrophysics, University of Oslo Introduction last decade has seen a large change in our understand e solar wind, due both theoretical,
More informationKinetic and Small Scale Solar Wind Physics
Chapter 11 Kinetic and Small Scale Solar Wind Physics Thus far the origin, evolution, and large scale characteristics of the solar wind have been addressed using MHD theory and observations. In this lecture
More informationAnisotropy and Alfvénicity of hourly fluctuations in the fast polar solar wind
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 109,, doi:10.1029/2003ja009947, 2004 Anisotropy and Alfvénicity of hourly fluctuations in the fast polar solar wind M. Neugebauer 1 Lunar and Planetary Laboratory,
More informationMHD turbulence in the solar corona and solar wind
MHD turbulence in the solar corona and solar wind Pablo Dmitruk Departamento de Física, FCEN, Universidad de Buenos Aires Turbulence, magnetic reconnection, particle acceleration Understand the mechanisms
More informationarxiv: v2 [astro-ph.sr] 27 Nov 2016
Data Constrained Coronal Mass Ejections in A Global Magnetohydrodynamics Model M. Jin 1,2, W. B. Manchester 3, B. van der Holst 3, I. Sokolov 3, G. Tóth 3, R. E. Mullinix 4, A. arxiv:1605.05360v2 [astro-ph.sr]
More informationPHYSICS BASIS FOR THE GASDYNAMIC MIRROR (GDM) FUSION ROCKET. Abstract
PHYSICS BASIS FOR THE GASDYNAMIC MIRROR (GDM) FUSION ROCKET T. Kammash Department of Nuclear Engineering and Radiological Sciences University of Michigan Ann Arbor, MI 48109 (313) 764-0205 W. Emrich Jr.
More informationCoronal Jets Simulated with the Alfvén Wave Solar Model
Coronal Jets Simulated with the Alfvén Wave Solar Model J. Szente 1, G. Toth 1, M. Ward 1, B. van der Holst 1, E. Landi 1, C.R. DeVore 2, T.I. Gombosi 1 judithsz@umich.edu Received ; accepted to be submitted
More informationarxiv: v1 [astro-ph.sr] 27 Nov 2016
Chromosphere to 1 AU Simulation of the 2011 March 7th Event: A Comprehensive Study of Coronal Mass Ejection Propagation M. Jin 1,2, W. B. Manchester 3, B. van der Holst 3, I. Sokolov 3, G. Tóth 3, A. Vourlidas
More informationMHD Flow Field and Momentum Transfer Process of Magneto-Plasma Sail
J. Plasma Fusion Res. SERIES, Vol. 8 (2009) MHD Flow Field and Momentum Transfer Process of Magneto-Plasma Sail Hiroyuki NISHIDA, Ikkoh FUNAKI, Yoshifumi INATANI 1) and Kanya KUSANO 2) University of Tokyo,
More informationMagnetic Reconnection in Laboratory, Astrophysical, and Space Plasmas
Magnetic Reconnection in Laboratory, Astrophysical, and Space Plasmas Nick Murphy Harvard-Smithsonian Center for Astrophysics namurphy@cfa.harvard.edu http://www.cfa.harvard.edu/ namurphy/ November 18,
More informationA Comparative Study of Different Approaches and Potential Improvement to Modeling the Solar Wind
A Comparative Study of Different Approaches and Potential Improvement to Modeling the Solar Wind Sun, X. and Hoeksema, J. T. W.W. Hansen Experimental Physics Laboratory (HEPL), Stanford University Abstract:
More informationWeek 8: Stellar winds So far, we have been discussing stars as though they have constant masses throughout their lifetimes. On the other hand, toward
Week 8: Stellar winds So far, we have been discussing stars as though they have constant masses throughout their lifetimes. On the other hand, toward the end of the discussion of what happens for post-main
More informationSolar Flare. A solar flare is a sudden brightening of solar atmosphere (photosphere, chromosphere and corona)
Solar Flares Solar Flare A solar flare is a sudden brightening of solar atmosphere (photosphere, chromosphere and corona) Flares release 1027-1032 ergs energy in tens of minutes. (Note: one H-bomb: 10
More informationSW103: Lecture 2. Magnetohydrodynamics and MHD models
SW103: Lecture 2 Magnetohydrodynamics and MHD models Scale sizes in the Solar Terrestrial System: or why we use MagnetoHydroDynamics Sun-Earth distance = 1 Astronomical Unit (AU) 200 R Sun 20,000 R E 1
More informationThe evolution of solar wind turbulence at kinetic scales
International Association of Geomagnetism and Aeronomy (IAGA) 2 nd Symposium: Solar Wind Space Environment Interaction c 2010 Cairo University Press December 4 th 8 th, 2009, Cairo, Egypt L.Damé & A.Hady
More informationHow is Earth s Radiation Belt Variability Controlled by Solar Wind Changes
How is Earth s Radiation Belt Variability Controlled by Solar Wind Changes Richard M. Thorne Department of Atmospheric and Oceanic Sciences, UCLA Electron (left) and Proton (right) Radiation Belt Models
More informationThe first super geomagnetic storm of solar cycle 24: The St. Patrick day (17 March 2015) event
The first super geomagnetic storm of solar cycle 24: The St. Patrick day (17 March 2015) event Chin Chun Wu 1, Kan Liou 2, Bernard Jackson 3, Hsiu Shan Yu 3, Lynn Hutting 1, R. P. Lepping 4, Simon Plunkett
More informationLandau damping in space plasmas with generalized r, q distribution function
PHYSICS OF PLASMAS 12, 122902 2005 Landau damping in space plasmas with generalized r, q distribution function M. N. S. Qureshi Key Laboratory of Space Weather, CSSAR, Chinese Academy of Sciences, Beijing
More informationObservations and Modeling of Turbulence in the Solar Wind
Observations and Modeling of Turbulence in the Solar Wind Melvyn L. Goldstein NASA Goddard Space Flight Center, USA E-mail: melvyn.l.goldstein@nasa.gov Summary. Alfvénic fluctuations are a ubiquitous component
More informationStudy of Wave-Particle Interaction Using Wind/ACE Data
Study of Wave-Particle Interaction Using Wind/ACE Data Lan Jian (lan.jian@nasa.gov) 1. 2. University of Maryland, College Park NASA Goddard Space Flight Center Collaborators: M. Stevens, S. P. Gary, A.
More informationDifferential velocity between solar wind protons and alpha particles in pressure balance structures
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 109,, doi:10.1029/2003ja010274, 2004 Differential velocity between solar wind protons and alpha particles in pressure balance structures Yohei Yamauchi and Steven
More informationCoronal Modeling and Synchronic Maps*
Coronal Modeling and Synchronic Maps* Jon A. Linker, Roberto Lionello, Zoran Mikic, Pete Riley, and Cooper Downs Predictive Science, Inc. (PSI), San Diego, CA 92121 http://www.predsci.com Carl Henney and
More informationSOLAR WIND MAGNETIC FLUCTUATIONS AND ELECTRON NON-THERMAL TEMPERATURE ANISOTROPY: SURVEY OF WIND-SWE-VEIS OBSERVATIONS
2016. The American Astronomical Society. All rights reserved. doi:10.3847/1538-4357/833/1/49 SOLAR WIND MAGNETIC FLUCTUATIONS AND ELECTRON NON-THERMAL TEMPERATURE ANISOTROPY: SURVEY OF WIND-SWE-VEIS OBSERVATIONS
More informationSolar wind modeling: a computational tool for the classroom. Lauren N. Woolsey. Harvard University. 60 Garden St, M.S. 10, Cambridge, MA 02138
Solar wind modeling: a computational tool for the classroom Lauren N. Woolsey Harvard University 60 Garden St, M.S. 10, Cambridge, MA 02138 lwoolsey@cfa.harvard.edu ABSTRACT: This article presents a Python
More information