No abstract

Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by di fferent physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and e fficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solarwind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a blockadaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional blockadaptive grid. We describe the algorithms of BATL and demonstrate its e fficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that di fferent time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics based space weather modeling and even forecasting.https://ntrs.nasa.gov/search.jsp?R=20110005631 2018-05-12T08:50:45+00:00Z

[1] The Space Weather Modeling Framework (SWMF) provides a high-performance flexible framework for physics-based space weather simulations, as well as for various space physics applications. The SWMF integrates numerical models of the Solar Corona, Eruptive Event Generator, Inner Heliosphere, Solar Energetic Particles, Global Magnetosphere, Inner Magnetosphere, Radiation Belt, Ionosphere Electrodynamics, and Upper Atmosphere into a high-performance coupled model. The components can be represented with alternative physics models, and any physically meaningful subset of the components can be used. The components are coupled to the control module via standardized interfaces, and an efficient parallel coupling toolkit is used for the pairwise coupling of the components. The execution and parallel layout of the components is controlled by the SWMF. Both sequential and concurrent execution models are supported. The SWMF enables simulations that were not possible with the individual physics models. Using reasonably high spatial and temporal resolutions in all of the coupled components, the SWMF runs significantly faster than real time on massively parallel supercomputers. This paper presents the design and implementation of the SWMF and some demonstrative tests. Future papers will describe validation (comparison of model results with measurements) and applications to challenging space weather events. The SWMF is publicly available to the scientific community for doing geophysical research. We also intend to expand the SWMF in collaboration with other model developers.

Abstract.It is well known that the ionosphere plays a role in determining the global state of the magnetosphere. The ionosphere allows magnetospheric currents to close, thereby allowing magnetospheric convection to occur. The amount of current which can be carried through the ionosphere is mainly determined by the ionospheric conductivity. This paper starts to quantify the nonlinear relationship between the ionospheric conductivity and the global state of the magnetosphere. It is found that the steady-state magnetosphere acts neither as a current nor as a voltage generator; a uniform Hall conductance can influence the potential pattern at low latitudes, but not at high latitude; the EUV generated conductance forces the currents to close in the sunlight, while the potential is large on the nightside; the solar generated Hall conductances cause a large asymmetry between the dawn and dusk potential, which effects the pressure distribution in the magnetosphere; a uniform polar cap potential removes some of this asymmetry; the potential difference between solar minimum and maximum is ∼11%; and the auroral precipitation can be related to the local field-aligned current through an exponential function.

We present a three-dimensional numerical magnetohydrodynamic simulation designed to model the emergence of a magnetic flux rope passing from below the photosphere into the corona. For the initial state, we prescribe a plane-parallel atmosphere that comprises a polytropic convection zone, photosphere, transition region, and corona. Embedded in this system is an isolated horizontal magnetic flux rope located 10 photospheric pressure scale heights below the photosphere. The flux rope is uniformly twisted, with the plasma temperature inside the rope reduced to compensate for the magnetic pressure. Density is reduced in the middle of the rope, so that this section buoyantly rises. The early evolution proceeds with the middle of the rope rising to the photosphere and expanding into the corona. Just as it seems the system might approach equilibrium, the upper part of the flux rope begins to separate from the lower, mass-laden part. The separation occurs through stretching of the field, which forms a current sheet, where reconnection severs the field lines to form a new system of closed flux. This flux then erupts into the corona. Essential to the eruption process are shearing motions driven by the Lorentz force, which naturally occur as the rope expands in the pressure-stratified atmosphere. The shearing motions transport axial flux and energy to the expanding portion of the magnetic field, driving the eruption. Subject headingg s: MHD -Sun: magnetic fields -Sun: photosphere

[1] We present a three-dimensional (3-D) numerical ideal magnetohydrodynamics (MHD) model describing the time-dependent expulsion of a coronal mass ejection (CME) from the solar corona propagating to 1 astronomical unit (AU). The simulations are performed using the Block Adaptive Tree Solar-Wind Roe Upwind Scheme (BATS-R-US) code. We begin by developing a global steady-state model of the corona that possesses high-latitude coronal holes and a helmet streamer structure with a current sheet at the equator. The Archimedean spiral topology of the interplanetary magnetic field is reproduced along with fast and slow speed solar wind. Within this model system, we drive a CME to erupt by the introduction of a Gibson-Low magnetic flux rope that is anchored at both ends in the photosphere and embedded in the helmet streamer in an initial state of force imbalance. The flux rope rapidly expands and is ejected from the corona with maximum speeds in excess of 1000 km/s. Physics-based adaptive mesh refinement (AMR) allows us to capture the structure of the CME focused on a particular Sun-Earth line with high spatial resolution given to the bow shock ahead of the flux rope as well as to the current sheet behind. The CME produces a large magnetic cloud at 1 AU (>100 R ) in which Bz undergoes a full rotation from north to south with an amplitude of 20 nT. In a companion paper, we find that the CME is very effective in generating strong geomagnetic activity at the Earth in two ways. First, through the strong sustained southward Bz (lasting more than 10 hours) and, second, by a pressure increase associated with the CMEdriven shock that compresses the magnetosphere.

[1] This paper describes the coupling of BATS-R-US (Block Adaptive Tree Solar-wind Roe-type Upwind Scheme), a magnetohydrodynamics (MHD) code representing the Earth's global magnetosphere and its coupling to the ionosphere and solar wind, and the Rice Convection Model (RCM), which represents the inner magnetosphere and its coupling to the ionosphere. The MHD code provides a time-evolving magnetic field model for the RCM as well as continuously updated boundary conditions for the electric potential and plasma. The RCM computes the distribution functions of inner magnetospheric particles, including transport of inner plasmasheet and ring current particles by gradient/curvature drifts; it thus calculates more accurate inner magnetospheric pressures, which are frequently passed to the MHD model and used to nudge the MHD values. Results are presented for an initial run with the coupled code for the case of uniform ionospheric conductance with steady solar wind and southward interplanetary magnetic field (IMF). The results are compared with those for a run of the MHD code alone. The coupled-code run shows significantly higher inner magnetospheric particle pressures. It also exhibits several well-established characteristics of inner magnetospheric electrodynamics, including strong region-2 Birkeland currents and partial shielding of the inner magnetosphere from the main force of the convection electric field. A sudden northward turning of the IMF causes the ring current to become more nearly symmetric. The inner magnetosphere exhibits an overshielding (dusk-to-dawn) electric field that begins about 10 min after the northward turning reaches the magnetopause and lasts just over an hour.

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.