Magnetohydrodynamic Modeling of Coronal Evolution and Disruption
In this report we describe progress in our investigation of CMEs using magnetohydrodynamic (MHD) simulations.
Linker, Jon A.↗
SEARCH · Search NASA
Results for “MHD simulation”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Design of a High-Energy, Two-Stage Pulsed Plasma Thruster
Design details of a proposed high-energy (approx. 50 kJ/pulse), two-stage pulsed plasma thruster are presented. The long-term goal of this project is to develop a high-power (approx. 500 kW), high specific impulse (approx. 7500 s), highly efficient (approx. 50%),and mechanically simple thruster for use as primary propulsion in a high-power nuclear electric propulsion system. The proposed thruster (PRC-PPT1) utilizes a valveless, liquid lithium-fed thermal plasma injector (first stage) followed by a high-energy pulsed electromagnetic accelerator (second stage). A numerical circuit model coupled with one-dimensional current sheet dynamics, as well as a numerical MHD simulation, are used to qualitatively predict the thermal plasma injection and current sheet dynamics, as well as to estimate the projected performance of the thruster. A set of further modelling efforts, and the experimental testing of a prototype thruster, is suggested to determine the feasibility of demonstrating a full scale high-power thruster.
Markusic, T. E.↗
Interplanetary Magnetic Field Control of the Entry of Solar Energetic Particles into the Magnetosphere
We have investigated the entry of energetic ions of solar origin into the magnetosphere as a function of the interplanetary magnetic field orientation. We have modeled this entry by following high energy particles (protons and 3 He ions) ranging from 0.1 to 50 MeV in electric and magnetic fields from a global magnetohydrodynamic (MHD) model of the magnetosphere and its interaction with the solar wind. For the most part these particles entered the magnetosphere on or near open field lines except for some above 10 MeV that could enter directly by crossing field lines due to their large gyroradii. The MHD simulation was driven by a series of idealized solar wind and interplanetary magnetic field (IMF) conditions. It was found that the flux of particles in the magnetosphere and transport into the inner magnetosphere varied widely according to the IMF orientation for a constant upstream particle source, with the most efficient entry occurring under southward IMF conditions. The flux inside the magnetosphere could approach that in the solar wind implying that SEPs can contribute significantly to the magnetospheric energetic particle population during typical SEP events depending on the state of the magnetosphere.
Richard, R. L.↗
The Interplanetary Magnetic Field and Magnetospheric Current Systems
We have performed systematic global magnetohydrodynamic (MHD) simulation studies driven by an idealized time series of solar wind parameters to establish basic cause and effect relationships between the solar wind variations and the ionosphere parameters. We studied six cases in which the interplanetary magnetic field (IMF) rotated from southward to northward in one minute. In three cases (cases A, B, and C) we ran five hours of southward IMF with Beta(sub Zeta) = 5 nT, followed by five hours of northward IMF with Beta(sub Zeta) = 5 nT. In the other three cases (cases D, E, and F) the magnetic field magnitude was increased to 10 nT. The solar wind parameters were: For cases A and D a density of 5 cm(exp -3), a thermal pressure of 3.3 nPa, and a solar wind speed 375 km/s, for cases B and E a density of 10 cm(exp -3), a thermal pressure of 9.9 nPa, and a solar wind speed 420 km/s, while for cases C and F a density of 15 cm(exp -3), a thermal pressure of 14.9 nPa, and a solar wind speed of 600 km/s.
El-Alaoui, Mostafa↗
Distinct Magnetospheric Responses to Southward IMF in Two Substorms
Solar wind plasma parameters and the Interplanetary Magnetic Field (IMF) observed by the WIND spacecraft upstream of the bow shock were used as input to magnetohydrodynamic (MHD) simulations of two substorm events. The power deposited into the ionosphere due to electron precipitation was calculated both from VIS observations and from the simulations.
El-Alaoui, Mostafa↗
Interaction of Strong Transient Interplanetary Disturbances with the Dayside Magnetosphere
The objective of the investigation was to gain an understanding of the complex response of the magnetosphere to strong transient interplanetary disturbances. Because the project was only funded for a year, the investigation focused on only one of the three topics proposed in the original three year proposal. We investigated the response of the dayside auroral region to strong transient interplanetary disturbances. The method of the investigation was to use three-dimensional magnetohydrodynamic (MHD) simulations that employ measurements of the solar wind conditions upstream the bow shock to model actual events and then to compare the simulation results with observations from spacecraft located at downstream locations. We modeled an event that occurred on July 14,2000, for which both CLUSTER and IMAGE simultaneous observations were available. The event was marked by high solar wind dynamic pressure and a strong IMF By component. Comparisons showed a very good agreement between intensifications in the auroral emissions measured by IMAGE
Berchem, Jean↗
Kinematic Treatment of CME Evolution in the Solar Wind
We present a kinematic study of the evolution of coronal mass ejections (CMEs) in the solar wind. Specifically, we consider the effects of: (1) spherical expansion; and (2) uniform expansion due to pressure gradients between the Interplanetary CME (ICME) and the ambient solar wind. We compare these results with an MHD model, which allows us to isolate these effects from the combined kinematic and dynamical effects, which are included in MHD models. They also provide compelling evidence that the fundamental cross section of so-called "force-free" flux ropes (or magnetic clouds) is neither circular or elliptical, but rather a convex-outward, "pancake" shape. We apply a force-free fitting to the magnetic vectors from the MHD simulation to assess how the distortion of the flux rope affects the fitting. In spite of these limitations, force-free fittings, which are straightforward to apply, do provide an important description of a number of parameters, including the radial dimension, orientation and chirality of the ICME.
Riley, Pete↗
Temporal-Spatial Structure of Magnetic Merging at the Magnetopause Inferred from 557.7-nm All-Sky Images
We demonstrate that high-resolution 557.7-nm all-sky images are useful tools for investigating the spatial and temporal evolution of merging on the dayside magnetopause. Analysis of ground and satellite measurements leads us to conclude that high-latitude merging events can occur at multiple sites simultaneously and vary asynchronously on time scales of 30 s to 3 min. Variations of 557.7 nm emissions were observed at a 10 s cadence at Ny-Alesund on 19 December 2001, while significant changes in the IMF clock angle were reaching the magnetopause. The optical patterns are consistent with a scenario in which merging occurs around the rim of the high-latitude cusp at positions dictated by the IMF clock angle. Electrons energized at merging sites represent plausible sources for 557.7 nm emissions in the cusp. Polar observations at the magnetopause have directly linked enhanced fluxes of > or = 0.5 keV electrons with merging. Spectra of electrons responsible for some of the emissions, measured during a DMSP F15 overflight, exhibit "inverted-V" features, indicating further acceleration above the ionosphere. SuperDARN spectral width boundaries, characteristic of open-closed field line transitions, are located at the equatorward edge of the 557.7nm emissions. Optical data suggest that with IMF B(sub Y) > 0, the Northern Hemisphere cusp divides into three source regions. When the IMF clock angle was approx. 150 deg structured 557.7-nm emissions came from east of the 13:00 MLT meridian. At larger clock angles the emissions appeared between 12:00 and 13:00 MLT. No significant 557.7-nm emissions were detected in the prenoon MLT sector. MHD simulations corroborate our scenario, showing that with the observed large dipole-tilt and IMF clock angles, merging sites develop near the front and eastern portions of the high-altitude cusp rim in the Northern Hemisphere and near the western part of the cusp rim in the Southern Hemisphere.
Maynard, N. C.↗
Bondi Accretion and the Problem of the Missing Isolated Neutron Stars
A large number of neutron stars (NSs), approximately 10(exp 9), populate the Galaxy, but only a tiny fraction of them is observable during the short radio pulsar lifetime. The majority of these isolated NSs, too cold to be detectable by their own thermal emission, should be visible in X-rays as a result of accretion from the interstellar medium. The ROSAT All-Sky Survey has, however, shown that such accreting isolated NSs are very elusive: only a few tentative candidates have been identified, contrary to theoretical predictions that up to several thousand should be seen. We suggest that the fundamental reason for this discrepancy lies in the use of the standard Bondi formula to estimate the accretion rates. We compute the expected source counts using updated estimates of the pulsar velocity distribution, realistic hydrogen atmosphere spectra, and a modified expression for the Bondi accretion rate, as suggested by recent MHD simulations and supported by direct observations in the case of accretion around supermassive black holes in nearby galaxies and in our own. We find that, whereas the inclusion of atmospheric spectra partly compensates for the reduction in the counts due to the higher mean velocities of the new distribution, the modified Bondi formula dramatically suppresses the source counts. The new predictions are consistent with a null detection at the ROSAT sensitivity.
Perna, Rosalba↗
Dynamical Evolution of the Inner Heliosphere Approaching Solar Activity Maximum: Interpreting Ulysses Observations Using a Global MHD Model
In this study we describe a series of MHD simulations covering the time period from 12 January 1999 to 19 September 2001 (Carrington Rotation 1945 to 1980). This interval coincided with: (1) the Sun s approach toward solar maximum; and (2) Ulysses second descent to the southern polar regions, rapid latitude scan, and arrival into the northern polar regions. We focus on the evolution of several key parameters during this time, including the photospheric magnetic field, the computed coronal hole boundaries, the computed velocity profile near the Sun, and the plasma and magnetic field parameters at the location of Ulysses. The model results provide a global context for interpreting the often complex in situ measurements. We also present a heuristic explanation of stream dynamics to describe the morphology of interaction regions at solar maximum and contrast it with the picture that resulted from Ulysses first orbit, which occurred during more quiescent solar conditions. The simulation results described here are available at: http://sun.saic.com.
Riley, Pete↗
Topological Evolution of a Fast Magnetic Breakout CME in 3-Dimensions
W present the extension of the magnetic breakout model for CME initiation to a fully 3-dimensional, spherical geometry. Given the increased complexity of the dynamic magnetic field interactions in 3-dimensions, we first present a summary of the well known axisymmetric breakout scenario in terms of the topological evolution associated with the various phases of the eruptive process. In this context, we discuss the completely analogous topological evolution during the magnetic breakout CME initiation process in the simplest 3-dimensional multipolar system. We show that an extended bipolar active region embedded in an oppositely directed background dipole field has all the necessary topological features required for magnetic breakout, i.e. a fan separatrix surface between the two distinct flux systems, a pair of spine fieldlines, and a true 3-dimensional coronal null point at their intersection. We then present the results of a numerical MHD simulation of this 3-dimensional system where boundary shearing flows introduce free magnetic energy, eventually leading to a fast magnetic breakout CME. The eruptive flare reconnection facilitates the rapid conversion of this stored free magnetic energy into kinetic energy and the associated acceleration causes the erupting field and plasma structure to reach an asymptotic eruption velocity of greater than or approx. equal to 1100 km/s over an approx.15 minute time period. The simulation results are discussed using the topological insight developed to interpret the various phases of the eruption and the complex, dynamic, and interacting magnetic field structures.
Lynch, B. J.↗
A Model for Solar Polar Jets
We propose a model for the jetting activity that is commonly observed in the Sun's corona, especially in the open-field regions of polar coronal holes. Magnetic reconnection is the process driving the jets and a relevant magnetic configuration is the well-known null point and fan separatrix topology. The primary challenge in explaining the observations is that reconnection must occur in a short-duration energetic burst rather than quasi-continuously as is implied by the observations of long-lived structures in coronal holes, such as polar plumes, for example. The key idea underlying our model for jets is that reconnection is forbidden for an axisymmetric null-point topology. Consequently, by imposing a twisting motion that maintains the axisymmetry, magnetic stress can be built up to large levels until an ideal instability breaks the symmetry and leads to an explosive release of energy via reconnection. Using 3D MHD simulations we demonstrate that this mechanism does produce jets with high speed and mass, driven by nonlinear Alfven waves. We discuss the implications of our results for observations of the solar corona.
Pariat, E.↗
STEREO/SECCHI Stereoscopic Observations Constraining the Initiation of Polar Coronal Jets
We report on the first stereoscopic observations of polar coronal jets made by the EUVI/SECCHI imagers on board the twin STEREO spacecraft. The significantly separated viewpoints (approximately 11 degrees ) allowed us to infer the 3D dynamics and morphology of a well-defined EUV coronal jet for the first time. Triangulations of the jet's location in simultaneous image pairs led to the true 3D position and thereby its kinematics. Initially the jet ascends slowly at approximately equal to 10-20 kilometers per second and then, after an apparent 'jump' takes place, it accelerates impulsively to velocities exceeding 300 kilometers per second with accelerations exceeding the solar gravity. Helical structure is the most important geometrical feature of the jet which shows evidence of untwisting. The jet structure appears strikingly different from each of the two STEREO viewpoints: face-on in the one viewpoint and edge-on in the other. This provides conclusive evidence that the observed helical structure is real and is not resulting from possible projection effects of single viewpoint observations. The clear demonstration of twisted structure in polar jets compares favorably with synthetic images from a recent MHD simulation of jets invoking magnetic untwisting as their driving mechanism. Therefore, the latter can be considered as a viable mechanism for the initiation of polar jets.
Patsourakos, S.↗
The Role of Flux Transfer Events in Solar Wind-Magnetosphere Coupling
Recent global MHD simulations of Flux Transfer Events (FTE's)[Dorelli and Bhattachar ee, JGR, 114, 2009] demonstrate that they are topologically complex flux ropes which extend large distances away from the subsolar magnetopause. Thus, FTE's represent a significant perturbation to the magnetopause magnetic field topology. Specifically, prior to the formation of the first FTE of the simulation, the dayside magnetopause has a relatively simple topology consisting of a single magnetic separator draping northward (even for southward IMF) over the dayside magnetopause (in a manner similar to the "overdraping" which occurs under northward IMF conditions). After the first FTE forms, however, the magnetopause becomes topologically more complex, consisting of multiple separators and braided flux domains. One consequence of this topological complexity is that open flux tubes have direct access to the magnetosphere near the subsolar region (such direct access is not possible during steady separator reconnection). Since the magnetopause topology affects the mapping of the solar wind electric field to the magnetosphere, an interesting question arises: What is the effect of FTE generation on the global magnetospheric convection pattern? In this talk, we address this question by exploring the differences between magnetospheric and ionospheric convection before and after the formation of the first FTE of the simulation.
Dorelli, John C.↗
Understanding Solar Eruptions with SDO/HMI Measuring Photospheric Flows, Testing Models, and Steps Towards Forecasting Solar Eruptions
The Solar Dynamics Observatory (SDO) is carrying the first full-disk imaging vector magnetograph, the Helioseismic and Magnetic Imager (HMI), into an inclined geosynchronous orbit. This magnetograph will provide nearly continuous measurements of photospheric vector magnetic fields at cadences of 90 seconds to 12 minutes with 1" resolution, precise pointing, and unfettered by atmospheric seeing. The enormous data stream of 1.5 Terabytes per day from SAO will provide an unprecedented opportunity to understand the mysteries of solar eruptions. These ground-breaking observations will permit the application of a new technique, the differential affine velocity estimator for vector magnetograms (DAVE4VM), to measure photospheric plasma flows in active regions. These measurements will permit, for the first time, accurate assessments of the coronal free energy available for driving CMEs and flares. The details of photospheric plasma flows, particularly along magnetic neutral-lines, are critical to testing models for initiating coronal mass ejections (CMEs) and flares. Assimilating flows and fields into state-of-the art 3D MHD simulations that model the highly stratified solar atmosphere from the convection zone to the corona represents the next step towards achieving NASA's Living with a Star forecasting goals of predicting "when a solar eruption leading to a CME will occur." Our presentation will describe these major science and predictive advances that will be delivered by SDO/HMI.
Schuck, Peter W.↗
The Evolution of Open Magnetic Flux Driven by Photospheric Dynamics
The coronal magnetic field is of paramount importance in solar and heliospheric physics. Two profoundly different views of the coronal magnetic field have emerged. In quasi-steady models, the predominant source of open magnetic field is in coronal holes. In contrast, in the interchange model, the open magnetic flux is conserved, and the coronal magnetic field can only respond to the photospheric evolution via interchange reconnection. In this view the open magnetic flux diffuses through the closed, streamer belt fields, and substantial open flux is present in the streamer belt during solar minimum. However, Antiochos and co-workers, in the form of a conjecture, argued that truly isolated open flux cannot exist in a configuration with one heliospheric current sheet (HCS) - it will connect via narrow corridors to the polar coronal hole of the same polarity. This contradicts the requirements of the interchange model. We have performed an MHD simulation of the solar corona up to 20R solar to test both the interchange model and the Antiochos conjecture. We use a synoptic map for Carrington Rotation 1913 as the boundary condition for the model, with two small bipoles introduced into the region where a positive polarity extended coronal hole forms. We introduce flows at the photospheric boundary surface to see if open flux associated with the bipoles can be moved into the closed-field region. Interchange reconnection does occur in response to these motions. However, we find that the open magnetic flux cannot be simply injected into closed-field regions - the flux eventually closes down and disconnected flux is created. Flux either opens or closes, as required, to maintain topologically distinct open and closed field regions, with no indiscriminate mixing of the two. The early evolution conforms to the Antiochos conjecture in that a narrow corridor of open flux connects the portion of the coronal hole that is nearly detached by one of the bipoles. In the later evolution, a detached coronal hole forms, in apparent violation of the Antiochos conjecture. Further investigation reveals that this detached coronal hole is actually linked to the extended coronal hole by a separatrix footprint on the photosphere of zero width. Therefore, the essential idea of the conjecture is preserved, if we modify it to state that coronal holes in the same polarity region are always linked, either by finite width corridors or separatrix footprints. The implications of these results for interchange reconnection and the sources of the slow solar wind are briefly discussed.
Linker, Jon A.↗
Reconnection Onset in the Breakout Model for CME Initiation
Fast coronal mass ejections (CMEs) are the most massive explosions in the heliosphere, and the primary drivers of geoeffective space weather. Although it is generally agreed that magnetic reconnection is the key to fast CME initiation, different models incorporate reconnection in different ways. One promising model --- the breakout scenario --- involves reconnection in two distinct yet interconnected locations: breakout reconnection ahead of the CME, and flare reconnection behind it. We will discuss what we have learned about the early evolution of breakout and flare reconnection from recent high-resolution 2.5D adaptively refined MHD simulations of CME initiation, including the evolving properties of the breakout and flare current sheets, the conditions that trigger reconnection onset in each sheet, the ensuing positive feedback between breakout and flare reconnections, and implications for electron acceleration in flares.
Karpen, Judy T.↗
CME Initiation and Reconnection
Coronal mass ejections (CMEs) are the most massive explosions in the heliosphere, and the primary drivers of geoeffective space weather. This talk will be focused on fast CMEs, which travel at Alfvenic speeds as high as 2500 km/s. These ejections are associated with solar flares, prominence eruptions, and energetic particles accelerated near the Sun and in interplanetary space. CMEs require sufficient energy storage, in the form of magnetic stress, and rapid release of this energy. Although it is generally agreed that magnetic reconnection is the key to fast CME initiation, different models incorporate reconnection in different ways. One promising model --- the breakout scenario --- involves reconnection in two distinct yet interconnected locations: breakout reconnection ahead of the CME, and flare reconnect ion behind it. This model has been validated through 2D and 3D MHD simulations and favorable comparison with the observed properties of many fast CMEs. I will discuss what we have learned about the onset and evolution of breakout and flare reconnect ion from recent high-resolution 2D simulations of CME initiation with adaptive mesh refinement and numerical resistivity.
Karpen, Judy T.↗