Igor V. Sokolov

Space Weather Modeling Framework: A new tool for the space science community

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

Adaptive numerical algorithms in space weather modeling

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 different relevant physics in different domains. A multi-physics system can be modeled by a software framework comprising 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 efficient 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-type Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamic (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 block-adaptive 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 block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency 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 different 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 short-term forecasting.

Three-dimensional, multispecies, high spatial resolution MHD studies of the solar wind interaction with Mars

[1] We present the results of model calculations, using our new, four-species, spherical MHD model. Our results are compared with the relevant and limited available data. The resulting comparisons help us to increase our understanding of the interaction processes between the solar wind and the Martian atmosphere/ionosphere. This new model with a spherical grid structure allowed us to use small (� 10 km) radial grid spacing in the ionospheric region. We found that the calculated bow shock positions agree reasonably well with the observed values. The calculated results vary with interplanetary magnetic field orientation, solar cycle conditions, and subsolar location. We found that our calculated ion densities, with parameters corresponding to solar cycle minimum conditions, reproduced the Viking 1 observed ion densities well. The calculated solar cycle maximum densities, above � 140 km, are also consistent with the appropriate Mars Global Surveyor radio occultation electron densities. Both the calculated solar cycle maximum and solar cycle minimum total transterminator and escape fluxes are significantly smaller than our previously published values. This decrease is due to the improved temperature values used for the recombination rates in this new model, which in turn results in lower ion densities and lower fluxes. INDEX TERMS: 2780 Magnetospheric Physics: Solar wind interactions with unmagnetized bodies; 6026 Planetology: Comets and Small Bodies: Ionospheres— composition and chemistry; 6028 Planetology: Comets and Small Bodies: Ionospheres—structure and dynamics; 2728 Magnetospheric Physics: Magnetosheath; KEYWORDS: Mars, MHD, bow shock, ionosphere, solar wind interaction

Hole Boring in a DT Pellet and Fast-Ion Ignition with Ultraintense Laser Pulses

Recently achieved high intensities of short laser pulses open new prospects in their application to hole boring in inhomogeneous overdense plasmas and for ignition in precompressed DT fusion targets. A simple analytical model and numerical simulations demonstrate that pulses with intensities exceeding 10;{22} W/cm;{2} may penetrate deeply into the plasma as a result of efficient ponderomotive acceleration of ions in the forward direction. The penetration depth as big as hundreds of microns depends on the laser fluence, which has to exceed a few tens of GJ/cm;{2}. The fast ions, accelerated at the bottom of the channel with an efficiency of more than 20%, show a high directionality and may heat the precompressed target core to fusion conditions.

A THREE-DIMENSIONAL FLUX ROPE MODEL FOR CORONAL MASS EJECTIONS BASED ON A LOSS OF EQUILIBRIUM

A series of simulation runs are carried out to investigate the loss of equilibrium of the three-dimensional flux rope configuration of Titov & Demoulin as a suitable mechanism for the initiation of coronal mass ejections. By means of these simulations, we are able to determine the conditions for which stable equilibria no longer exist. Our results imply that it is possible to achieve a loss of equilibrium even though the ends of the flux rope are anchored to the solar surface. However, in order to have the flux rope escape, it is necessary to modify the configuration by eliminating the arcade field.

A Semiempirical Magnetohydrodynamical Model of the Solar Wind

We present a new MHD model for simulating the large-scale structure of the solar corona and solar wind under “steady state” conditions stemming from the Wang-Sheeley-Arge empirical model. The processes of turbulent heating in the solar wind are parameterized using a phenomenological, thermodynamical model with a varied polytropic index. We employ the Bernoulli integral to bridge the asymptotic solar wind speed with the assumed distribution of the polytropic index on the solar surface. We successfully reproduce the mass flux from Sun to Earth, the temperature structure, and the large-scale structure of the magnetic field. We reproduce the solar wind speed bimodal structure in the inner heliosphere. However, the solar wind speed is in a quantitative agreement with observations at 1 AU for solar maximum conditions only. The magnetic field comparison demonstrates that the input magnetogram needs to be multiplied by a scaling factor in order to obtain the correct magnitude at 1 AU.

Alfvén wave solar model (AWSoM): Coronal heating

We present a new version of the Alfven 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 Alfven wave turbulence. The injection of Alfven 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 Alfven waves are partially reflected by the Alfven 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.

Three-dimensional MHD Simulation of the 2003 October 28 Coronal Mass Ejection: Comparison with LASCO Coronagraph Observations

We numerically model the coronal mass ejection (CME) event of 2003 October 28 that erupted from AR 10486 and propagated to Earth in less than 20 hr, causing severe geomagnetic storms. The magnetohydrodynamic (MHD) model is formulated by first arriving at a steady state corona and solar wind employing synoptic magnetograms. We initiate two CMEs from the same active region, one approximately a day earlier that preconditions the solar wind for the much faster CME on the 28th. This second CME travels through the corona at a rate of over 2500 km s−1, driving a strong forward shock. We clearly identify this shock in an image produced by the Large Angle Spectrometric Coronagraph (LASCO) C3 and reproduce the shock and its appearance in synthetic white-light images from the simulation. We find excellent agreement with both the general morphology and the quantitative brightness of the model CME with LASCO observations. These results demonstrate that the CME shape is largely determined by its interaction with the ambient solar wind and may not be sensitive to the initiation process. We then show how the CME would appear as observed by wide-angle coronagraphs on board the Solar Terrestrial Relations Observatory (STEREO) spacecraft. We find complex time evolution of the white-light images as a result of the way in which the density structures pass through the Thomson sphere. The simulation is performed with the Space Weather Modeling Framework (SWMF).

CORONAL MASS EJECTION SHOCK AND SHEATH STRUCTURES RELEVANT TO PARTICLE ACCELERATION

Most high-energy solar energetic particles are believed to be accelerated at shock waves driven by coronal mass ejections (CMEs). The acceleration process strongly depends on the shock geometry and the structure of the sheath that forms behind the shock. In an effort to understand the structure and time evolution of such CME-driven shocks andtheirrelevancetoparticleacceleration,weinvestigatetheinteractionofafastCMEwiththeambientsolarwind by means of a three-dimensional numerical ideal MHD model. Our global steady state coronal model possesses high-latitudecoronalholesandahelmetstreamerstructurewithacurrentsheetneartheequator,reminiscentofnear solar minimum conditions. Fast and slow solar winds flow at high and low latitude, respectively, and the Archimedean spiral geometry of the interplanetary magnetic field is reproduced by solar rotation. Within this model system, we drive a CME to erupt by introducing a Gibson-Low magnetic flux rope that is 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 � 1 , driving a fast-mode shock from the inner corona to a distance of 1 AU. We find that the ambient solar wind structure strongly affects the evolution of the CME-driven shocks, causing deviations of the fast-mode shocks from their expected global configuration. These deflections lead to substantial compressions of the plasma and magnetic field in their associated sheath region. The sudden postshock increase in magneticfieldstrengthonlow-latitudefieldlinesisfoundtobeeffectiveforacceleratingparticlestotheGeVrange. Subject heading gs: acceleration of particles — MHD — shock waves — Sun: coronal mass ejections (CMEs)

Relativistic laser piston model: Ponderomotive ion acceleration in dense plasmas using ultraintense laser pulses

Laser ponderomotive force at superhigh intensities provides an efficient ion acceleration in bulk dense targets and evacuates a channel enabling further laser beam propagation. The developed quasistationary model of a laser piston—a double layer structure supported by the radiation pressure—predicts the general parameters of the acceleration process in homogeneous and inhomogeneous overdense plasmas. Particle-in-cell simulations confirm the estimated characteristics in a wide range of laser intensities and ion densities and show advantages of circularly polarized laser pulses. Two nonstationary effects are identified in the simulations. First, oscillations of the piston velocity and of the thickness of the ion charge separation layer broaden the energy spectrum of accelerated ions. Second, the electrons accelerated toward the incoming laser wave emit strong high-frequency radiation, enabling a cooling effect, which helps to sustain high charge neutrality in the piston and to maintain an efficient ion acce...