Models for the origin of the slow solar wind must account for two seemingly contradictory observations: the slow wind has the composition of the closed-field corona, implying that it originates from ...the continuous opening and closing of flux at the boundary between open and closed field. On the other hand, the slow wind also has large angular width, up to ~60?, suggesting that its source extends far from the open-closed boundary. We propose a model that can explain both observations. The key idea is that the source of the slow wind at the Sun is a network of narrow (possibly singular) open-field corridors that map to a web of separatrices and quasi-separatrix layers in the heliosphere. We compute analytically the topology of an open-field corridor and show that it produces a quasi-separatrix layer in the heliosphere that extends to angles far from the heliospheric current sheet. We then use an MHD code and MDI/SOHO observations of the photospheric magnetic field to calculate numerically, with high spatial resolution, the quasi-steady solar wind, and magnetic field for a time period preceding the 2008 August 1 total solar eclipse. Our numerical results imply that, at least for this time period, a web of separatrices (which we term an S-web) forms with sufficient density and extent in the heliosphere to account for the observed properties of the slow wind. We discuss the implications of our S-web model for the structure and dynamics of the corona and heliosphere and propose further tests of the model.
Sympathetic eruptions on the Sun have been observed for several decades, but the mechanisms by which one eruption can trigger another remain poorly understood. We present a three-dimensional MHD ...simulation that suggests two possible magnetic trigger mechanisms for sympathetic eruptions. We consider a configuration that contains two coronal flux ropes located within a pseudo-streamer and one rope located next to it. A sequence of eruptions is initiated by triggering the eruption of the flux rope next to the streamer. The expansion of the rope leads to two consecutive reconnection events, each of which triggers the eruption of a flux rope by removing a sufficient amount of overlying flux. The simulation qualitatively reproduces important aspects of the global sympathetic event on 2010 August 1 and provides a scenario for the so-called twin filament eruptions. The suggested mechanisms are also applicable for sympathetic eruptions occurring in other magnetic configurations.
In solar coronal loops, thermal non-equilibrium (TNE) is a phenomenon that can occur when the heating is both highly stratified and quasi-constant. Unambiguous observational identification of TNE ...would thus permit us to strongly constrain heating scenarios. While TNE is currently the standard interpretation of coronal rain, the long-term periodic evolution predicted by simulations has never been observed. However, the detection of long-period intensity pulsations (periods of several hours) has been recently reported with the Solar and Heliospheric Observatory/EIT, and this phenomenon appears to be very common in loops. Moreover, the three intensity-pulsation events that we recently studied with the Solar Dynamics Observatory/Atmospheric Imaging Assembly (AIA) show strong evidence for TNE in warm loops. In this paper, a realistic loop geometry from linear force-free field (LFFF) extrapolations is used as input to 1D hydrodynamic simulations. Our simulations show that, for the present loop geometry, the heating has to be asymmetrical to produce TNE. We analyze in detail one particular simulation that reproduces the average thermal behavior of one of the pulsating loop bundle observed with AIA. We compare the properties of this simulation with those deduced from the observations. The magnetic topology of the LFFF extrapolations points to the presence of sites of preferred reconnection at one footpoint, supporting the presence of asymmetric heating. In addition, we can reproduce the temporal large-scale intensity properties of the pulsating loops. This simulation further strengthens the interpretation of the observed pulsations as signatures of TNE. This consequently provides important information on the heating localization and timescale for these loops.
Abstract Observations have shown a clear association of filament/prominence eruptions with the emergence of magnetic flux in or near filament channels. Magnetohydrodynamic (MHD) simulations have been ...employed to systematically study the conditions under which such eruptions occur. These simulations to date have modeled filament channels as 2D flux ropes or 3D uniformly sheared arcades. Here we present MHD simulations of flux emergence into a more realistic configuration consisting of a bipolar active region containing a line-tied 3D flux rope. We use the coronal flux-rope model of Titov et al. as the initial condition and drive our simulations by imposing boundary conditions extracted from a flux emergence simulation by Leake et al. We identify three mechanisms that determine the evolution of the system: (i) reconnection displacing footpoints of field lines overlying the coronal flux rope, (ii) changes of the ambient field due to the intrusion of new flux at the boundary, and (iii) interaction of the (axial) electric currents in the preexisting and newly emerging flux systems. The relative contributions and effects of these mechanisms depend on the properties of the preexisting and emerging flux systems. Here we focus on the location and orientation of the emerging flux relative to the coronal flux rope. Varying these parameters, we investigate under which conditions an eruption of the latter is triggered.
The large-scale, steady-state magnetic field configuration of the solar corona is typically computed using boundary conditions derived from photospheric observations. Two approaches are typically ...used: (1) potential field source surface (PFSS) models, and (2) the magnetohydrodynamic (MHD) models. The former have the advantage that they are simple to develop and implement, require relatively modest computer resources, and can resolve structure on scales beyond those that can be handled by current MHD models. However, they have been criticized because their basic assumptions are seldom met. Moreover, PFSS models cannot directly incorporate time-dependent phenomena, such as magnetic reconnection, and do not include plasma or its effects. In this study, we assess how well PFSS models can reproduce the large-scale magnetic structure of the corona by making detailed comparisons with MHD solutions at different phases in the solar activity cycle. In particular, we (1) compute the shape of the source surface as inferred from the MHD solutions to assess deviations from sphericity, (2) compare the coronal hole boundaries as determined from the two models, and (3) estimate the effects of nonpotentiality. Our results demonstrate that PFSS solutions often closely match MHD results for configurations based on untwisted coronal fields (i.e., when driven by line-of-sight magnetograms). It remains an open question whether MHD solutions will differ more substantially from PFSS solutions when vector magnetograms are used as boundary conditions. This will be addressed in the near future when vector data from SOLIS, the Solar Dynamics Observatory, and Solar-B become incorporated into the MHD models.
In recent work, Antiochos and coworkers argued that the boundary between the open and closed field regions on the Sun can be extremely complex with narrow corridors of open flux connecting seemingly ...disconnected coronal holes from the main polar holes and that these corridors may be the sources of the slow solar wind. We examine, in detail, the topology of such magnetic configurations using an analytical source surface model that allows for analysis of the field with arbitrary resolution. Our analysis reveals three new important results. First, a coronal hole boundary can join stably to the separatrix boundary of a parasitic polarity region. Second, a single parasitic polarity region can produce multiple null points in the corona and, more important, separator lines connecting these points. It is known that such topologies are extremely favorable for magnetic reconnection, because they allow this process to occur over the entire length of the separators rather than being confined to a small region around the nulls. Finally, the coronal holes are not connected by an open-field corridor of finite width, but instead are linked by a singular line that coincides with the separatrix footprint of the parasitic polarity. We investigate how the topological features described above evolve in response to the motion of the parasitic polarity region. The implications of our results for the sources of the slow solar wind and for coronal and heliospheric observations are discussed.
Long-period EUV pulsations, recently discovered to be common in active regions, are understood to be the coronal manifestation of thermal nonequilibrium (TNE). The active regions previously studied ...with EIT/Solar and Heliospheric Observatory and AIA/SDO indicated that long-period intensity pulsations are localized in only one or two loop bundles. The basic idea of this study is to understand why. For this purpose, we tested the response of different loop systems, using different magnetic configurations, to different stratifications and strengths of the heating. We present an extensive parameter-space study using 1D hydrodynamic simulations (1020 in total) and conclude that the occurrence of TNE requires specific combinations of parameters. Our study shows that the TNE cycles are confined to specific ranges in parameter space. This naturally explains why only some loops undergo constant periodic pulsations over several days: since the loop geometry and the heating properties generally vary from one loop to another in an active region, only the ones in which these parameters are compatible exhibit TNE cycles. Furthermore, these parameters (heating and geometry) are likely to vary significantly over the duration of a cycle, which potentially limits the possibilities of periodic behavior. This study also confirms that long-period intensity pulsations and coronal rain are two aspects of the same phenomenon: both phenomena can occur for similar heating conditions and can appear simultaneously in the simulations.
GPU accelerators have had a notable impact on high-performance computing across many disciplines. They provide high performance with low cost/power, and therefore have become a primary compute ...resource on many of the largest supercomputers. Here, we implement multi-GPU acceleration into our Solar MHD code (MAS) using OpenACC in a fully portable, single-source manner. Our preliminary implementation is focused on MAS running in a reduced physics "zero-beta" mode. While valuable on its own, our main goal is to pave the way for a full physics, thermodynamic MHD implementation. We describe the OpenACC implementation methodology and challenges. "Time-to-solution" performance results of a production-level flux rope eruption simulation on multi-CPU and multi-GPU systems are shown. We find that the GPU-accelerated MAS code has the ability to run "zero-beta" simulations on a single multi-GPU server at speeds previously requiring multiple CPU server-nodes of a supercomputer.
The solar magnetic field is the primary agent that drives solar activity and couples the Sun to the hellosphere. Although the details of this coupling depend on the quantitative properties of the ...field, many important aspects of the corona-solar wind connection can be understood by considering only the general topological properties of those regions on the Sun where the field extends from the photosphere out to interplanetary space, the so-called open field regions that are usually observed as coronal holes. From the simple assumptions that underlie the standard quasi-steady corona-wind theoretical models, and that are likely to hold for the Sun as well, we derive two conjectures as to the possible structure and dynamics of coronal holes: (1) coronal holes are unique In that every unipolar region on the photosphere can contain at most one coronal hole, and (2) coronal holes of nested polarity regions must themselves be nested. Magnetic reconnection plays the central role in enforcing these constraints on the field topology. From these conjectures we derive additional properties for the topology of open field regions, and propose several observational predictions for both the slowly varying and transient corona/solar wind.
In this study we describe a three‐dimensional MHD model of the solar corona and heliosphere. We split the modeling region into two distinct parts: the solar corona (1 solar radius, RS, to 30 RS) and ...the inner heliosphere (30 RS to 5 AU). This combined model is driven solely by the observed line‐of‐sight photospheric magnetic field and can thus provide a realistic global picture of the corona and heliosphere for specific time periods of interest. We use the model to illustrate heliospheric structure during three different phases of the solar cycle: (1) Carrington rotation (CR) 1913 (August 22, 1996, to September 18, 1996), which occurred near solar minimum and overlapped the “Whole Sun Month” campaign; (2) CR 1892 (January 27, 1995, to February 23, 1995), which occurred during the declining phase of cycle 22 and coincided with the so‐called “Ulysses rapid latitude scan” and (3) CR 1947 (March 7, 1999, to April 4, 1999), which occurred approximately 15 months before the predicted maximum of solar cycle 23. We compare Ulysses and Wind observations with the simulation for CR 1913 and compare Ulysses observations during its traversal from pole to pole with CR 1892. We find that the simulations reproduce the overall large‐scale features of the observations. We use the near‐solar‐maximum results to speculate on the structure of the high‐latitude solar wind that Ulysses will encounter during its traversal of the southern and northern solar poles in 2000 and 2001, respectively. In particular, the results suggest that because of the presence of equatorial coronal holes the ordered pattern of corotating interaction region tilts and their associated shocks, which was observed during Ulysses' initial southward excursion in 1992, will likely disappear completely as Ulysses moves toward the south pole. We anticipate that Ulysses will encounter fast streams but will not remain within them for more than a fraction of a solar rotaton. Finally, the simulations suggest that crossings of the heliospheric current sheet will persist up to at least ∼70° heliographic latitude.