Abstract. A parallel adaptive mesh refinement (AMR) finite-volume scheme for predictingideal MHD flows is used to simulate the initiation, structure, and evolution of a coronal mass ejection (CME) and its interaction with the magnetosphere-ionosphere system. The simulated CME is driven by a local plasma density enhancement on the solar surface with the background initial state of the corona and solar wind represented by a newly devised "steady state" solution. The initial solution has been constructed to provide a reasonable description of the time-averaged solar wind for conditions near solar minimum: (1) the computed magnetic field near the Sun possesses high-latitude polar coronal holes, closed magnetic field flux tubes at low latitudes, and a helmet streamer structure with a neutral line and current sheet; (2) the Archimedean spiral topology of the interplanetary magnetic field is reproduced; (3) the observed two-state nature of the solar wind is also reproduced with the simulation yielding fast and slow solar wind streams at high and low latitudes, respectively; and (4) the predicted solar wind plasma properties at 1 AU are consistent with observations. Starting with the generation of a CME at the Sun, the simulation follows the evolution of the solar wind disturbance as it evolves into a magnetic cloud and travels through interplanetary space and subsequently interacts with the terrestrial magnetosphere-ionosphere system. The density-driven CME exhibits a two-step release process, with the front of the CME rapidly accelerating following the disruption of the near-Sun closed magnetic field line structure and then moving at a nearly constant speed of •560 km/s through interplanetary space. The CME also produces a large magnetic cloud (> 100/•s across) characterized by a magnetic field that smoothly rotates northward and then back again over a period of •2 days at 1 AU. The cloud does not contain a sustained period with a strong southward component of the magnetic field, and, as a consequence, the simulated CME is somewhat ineffective in generating strong geo-magnetic activity at Earth. Nevertheless, the simulation results illustrate the potential, as well as current limitations, of the MHD-based space weather model for enhancing the understanding of coronal physics, solar wind plasma processes, magnetospheric physics, and space weather phenomena. Such models will provide the foundation for future, more comprehensive space weather prediction tools.

Abstract. The solar wind-magnetosphere interaction for northward interplanetary magnetic field (IMF) is studied using a newly developed three-dimensional adaptive mesh refinement (AMR) global MHD simulation model. The simulations show that for northward IMF the magnetosphere is essentially closed. Reconnection between the IMF and magnetospheric field is limited to finite regions near the cusps. When the reconnection process forms newly closed magnetic field lines on the dayside, the solar wind plasma trapped on these reconnected magnetic field lines becomes part of the low-latitude boundary layer (LLBL) plasma and it convects to the nightside along the magnetopause. The last closed magnetic field line marks the topological boundary of the magnetospheric domain. When the last closed magnetic field line disconnects at the cusps and reconnects to the IMF, its plasma content becomes part of the solar wind. Plasma convection in the outer magnetosphere does not directly contribute to the reconnection process. On the dayside the topological boundary between the solar wind and the magnetosphere is located at the inner edge of the magnetopause current layer. At the same time, multiple current layers are observed in the high-altitude cusp region. Our convergence study and diagnostic analysis indicate that the details of the diffusion and the viscous interaction do not play a significant role in controlling the large-scale configuration of the simulated magnetosphere. It is sufficient that these dissipation mechanisms exist in the simulations. In our series of simulations the length of the magnetotail is primarily determined by the balance between the boundary layer driving forces and the drag forces. With a parametric study, we find that the tail length is proportional to the magnetosheath plasma beta near the magnetopause at local noon. A higher solar wind density, weaker IMF, and larger solar wind Mach number results in a longer tail. On the nightside downstream of the last closed magnetic field line the plasma characteristics are similar to that in the magnetotail, posing an observational challenge for identification of the topological status of the corresponding field lines.

No abstract

A 3-D parallel adaptive mesh refinement (AMR) scheme is described for solving the partial-differential equations governing ideal magnetohydrodynamic (MHD) flows. This new algorithm adopts a cell-centered upwind finite-volume discretization procedure and uses limited solution reconstruction, approximate Riemann solvers, and explicit multi-stage time stepping to solve the MHD equations in divergence form, providing a combination of high solution accuracy and computational robustness across a large range in the plasma ( is the ratio of thermal and magnetic pressures). The data structure naturally lends itself to domain decomposition, thereby enabling efficient and scalable implementations on massively parallel supercomputers. Numerical results for MHD simulations of magnetospheric plasma flows are described to demonstrate the validity and capabilities of the approach for space weather applications.