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.

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

We present a new version of the Alfvén wave solar model, a global model from the upper chromosphere to the corona and the heliosphere. The coronal heating and solar wind acceleration are addressed with low-frequency Alfvén wave turbulence. The injection of Alfvén wave energy at the inner boundary is such that the Poynting flux is proportional to the magnetic field strength. The three-dimensional magnetic field topology is simulated using data from photospheric magnetic field measurements. This model does not impose open-closed magnetic field boundaries; those develop self-consistently. The physics include the following.(1) The model employs three different temperatures, namely the isotropic electron temperature and the parallel and perpendicular ion temperatures. The firehose, mirror, and ion-cyclotron instabilities due to the developing ion temperature anisotropy are accounted for. (2) The Alfvén waves are partially reflected by the Alfvén speed gradient and the vorticity along the field lines. The resulting counter-propagating waves are responsible for the nonlinear turbulent cascade. The balanced turbulence due to uncorrelated waves near the apex of the closed field lines and the resulting elevated temperatures are addressed. (3) To apportion the wave dissipation to the three temperatures, we employ the results of the theories of linear wave damping and nonlinear stochastic heating. (4) We have incorporated the collisional and collisionless electron heat conduction. We compare the simulated multi-wavelength extreme ultraviolet images of CR2107 with the observations from STEREO/EUVI and the Solar Dynamics Observatory/AIA instruments. We demonstrate that the reflection due to strong magnetic fields in the proximity of active regions sufficiently intensifies the dissipation and observable emission.

[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.