Search NASA⌕ Search

SEARCH · Search NASA

Results for “ideal MHD”

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.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Initiation of Coronal Mass Ejections

This paper is a synopsis of the initiation of the strong-field magnetic explosions that produce large, fast coronal mass ejections. Cartoons based on observations are used to describe the inferred basic physical processes and sequences that trigger and drive the explosion. The magnetic field that explodes is a sheared-core bipole that may or may not be embedded in surrounding strong magnetic field, and may or may not contain a flux rope before it starts to explode. We describe three different mechanisms that singly or in combination trigger the explosion: (1) runaway internal tether-cutting reconnection, (2) runaway external tether-cutting reconnection, and (3) ideal MHD instability or loss or equilibrium. For most eruptions, high-resolution, high-cadence magnetograms and chromospheric and coronal movies (such as from TRACE and/or Solar-B) of the pre-eruption region and of the onset of the eruption and flare are needed to tell which one or which combination of these mechanisms is the trigger. Whatever the trigger, it leads to the production of an erupting flux rope. Using a simple model flux rope, we demonstrate that the explosion can be driven by the magnetic pressure of the expanding flux rope, provided the shape of the expansion is "fat" enough.

Moore, Ronald L.↗
NASA NTRS2005

Shape and Reconnection of the Exploding Magnetic Field in the Onset of CMEs

From chromospheric and coronal images and line-of-sight and vector magnetograms of magnetic regions that produce CMEs, and from chromospheric and coronal movies of the onsets of CME eruptions, it appears that the magnetic field that explodes to drive the CME is initially the strongly sheared core of a magnetic arcade encasing a polarity dividing line in the magnetic flux. Before or during the onset of the explosion, the sheared core field becomes a flux rope, often carrying chromospheric material within it. For the erupting flux rope to drive the explosion, that is, for its magnetic energy content to decrease in the explosion, the flux rope's cross-sectional area must increase faster than its length. For instance, for isotropic expansion, the area increases as the square of the length, and the magnetic energy content of the flux rope decreases as the inverse of the length. The instability that initiates the eruption of the flux rope might be an ideal MHD kink instability, or might involve runaway tether-cutting reconnection. The reconnection begins below the flux rope (internal to the arcade) when the overall field configuration of the region is effectively that of a single bipole. When the flux rope resides in a multi- bipolar configuration having a magnetic null above the flux rope, the runaway tether- cutting reconnection might begin either below the flux rope or at the null above (external to) the arcade. We present examples of observed CME onsets that illustrate the above alternatives. In each example, reconnection below the flux rope begins early in the eruption. This indicates that internal tether cutting reconnection (classic tether-cutting reconnection) is important in unleashing the CME explosion in all cases, including those in which the explosion may be triggered by MHD kinking or by external reconnection (classic breakout reconnection).

Moore, R. L.↗
NASA NTRS2006

Initiation of Solar Eruptions: Recent Observations and Implications for Theories

Solar eruptions involve the violent disruption of a system of magnetic field. Just how the field is destabilized and explodes to produce flares and coronal mass ejections (CMEs) is still being debated in the solar community. Here I discuss recent observational work into these questions by ourselves (me and my colleagues) and others. Our work has concentrated mainly on eruptions that include filaments. We use the filament motion early in the event as a tracer of the motion of the general erupting coronal field in and around the filament, since that field itself is hard to distinguish otherwise. Our main data sources are EUV images from SOHO/EIT and TRACE, soft Xray images from Yohkoh, and magnetograms from SOHO/MDI, supplemented with coronagraph images from SOHO/LASCO, hard X-ray data, and ground-based observations. We consider the observational findings in terms of three proposed eruption-initiation mechanisms: (i) runaway internal tether-cutting reconnection, (ii) slow external tether-cutting reconnection ("breakout"), and (iii) ideal MHD instability.

Sterling, A. C.↗
NASA NTRS2008

Simulating Sources of Superstorm Plasmas

We evaluated the contributions to magnetospheric pressure (ring current) of the solar wind, polar wind, auroral wind, and plasmaspheric wind, with the surprising result that the main phase pressure is dominated by plasmaspheric protons. We used global simulation fields from the LFM single fluid ideal MHD model. We embedded the Comprehensive Ring Current Model within it, driven by the LFM transpolar potential, and supplied with plasmas at its boundary including solar wind protons, polar wind protons, auroral wind O+, and plasmaspheric protons. We included auroral outflows and acceleration driven by the LFM ionospheric boundary condition, including parallel ion acceleration driven by upward currents. Our plasmasphere model runs within the CRCM and is driven by it. Ionospheric sources were treated using our Global Ion Kinetics code based on full equations of motion. This treatment neglects inertial loading and pressure exerted by the ionospheric plasmas, and will be superceded by multifluid simulations that include those effects. However, these simulations provide new insights into the respective role of ionospheric sources in storm-time magnetospheric dynamics.

Fok, Mei-Ching↗
NASA NTRS2010

The Writhe of Helical Structures in the Solar Corona

Context. Helicity is a fundamental property of magnetic fields, conserved in ideal MHD. In flux rope topology, it consists of twist and writhe helicity. Despite the common occurrence of helical structures in the solar atmosphere, little is known about how their shape relates to the writhe, which fraction of helicity is contained in writhe, and how much helicity is exchanged between twist and writhe when they erupt. Aims. Here we perform a quantitative investigation of these questions relevant for coronal flux ropes. Methods. The decomposition of the writhe of a curve into local and nonlocal components greatly facilitates its computation. We use it to study the relation between writhe and projected S shape of helical curves and to measure writhe and twist in numerical simulations of flux rope instabilities. The results are discussed with regard to filament eruptions and coronal mass ejections (CMEs). Results. (1) We demonstrate that the relation between writhe and projected S shape is not unique in principle, but that the ambiguity does not affect low-lying structures, thus supporting the established empirical rule which associates stable forward (reverse) S shaped structures low in the corona with positive (negative) helicity. (2) Kink-unstable erupting flux ropes are found to transform a far smaller fraction of their twist helicity into writhe helicity than often assumed. (3) Confined flux rope eruptions tend to show stronger writhe at low heights than ejective eruptions (CMEs). This argues against suggestions that the writhing facilitates the rise of the rope through the overlying field. (4) Erupting filaments which are S shaped already before the eruption and keep the sign of their axis writhe (which is expected if field of one chirality dominates the source volume of the eruption), must reverse their S shape in the course of the rise. Implications for the occurrence of the helical kink instability in such events are discussed.

WRITHE↗
NASA NTRS2012

The Extended Pulsar Magnetosphere

We present the structure of the 3D ideal MHD pulsar magnetosphere to a radius ten times that of the light cylinder, a distance about an order of magnitude larger than any previous such numerical treatment. Its overall structure exhibits a stable, smooth, well-defined undulating current sheet which approaches the kinematic split monopole solution of Bogovalov 1999 only after a careful introduction of diffusivity even in the highest resolution simulations. It also exhibits an intriguing spiral region at the crossing of two zero charge surfaces on the current sheet, which shows a destabilizing behavior more prominent in higher resolution simulations. We discuss the possibility that this region is physically (and not numerically) unstable. Finally, we present the spiral pulsar antenna radiation pattern.

Constantinos, Kalapotharakos↗
NASA NTRS2014

New Aspects of a Lid-Removal Mechanism in the Onset of a SEP-Producing Eruption Sequence

We examine a sequence of two ejective eruptions from a single active region on 2012 January 23, using magnetograms and EUV images from SDO/HMI and SDO/AIA, and EUV images from STEREO. Cheng et al. (2013) showed that the first eruption's ("Eruption 1'') flux rope was apparent only in "hotter'' AIA channels, and that it removed overlying field that allowed the second eruption (``Eruption 2'') to begin via ideal MHD instability; here we say Eruption 2 began via a ``lid removal'' mechanism. We show that during Eruption-1's onset, its flux rope underwent ``tether weakening'' (TW) reconnection with the field of an adjacent active region. Standard flare loops from Eruption 1 developed over Eruption-2's flux rope and enclosed filament, but these overarching new loops were unable to confine that flux rope/filament. Eruption-1's flare loops, from both TW reconnection and standard-flare-model internal reconnection, were much cooler than Eruption-2's flare loops (GOES thermal temperatures of approx. 9 MK compared to approx. 14 MK). This eruption sequence produced a strong solar energetic particle (SEP) event (10 MeV protons, >10(exp 3) pfu for 43 hrs), apparently starting when Eruption-2's CME blasted through Eruption-1's CME at 5-10 R_s. This occurred because the two CMEs originated in close proximity and in close time sequence: Eruption-1's fast rise started soon after the TW reconnection; the lid removal by Eruption-1's ejection triggered the slow onset of Eruption 2; and Eruption-2's CME, which started approx. 1 hr later, was three times faster than Eruption-1's CME.

Sterling, Alphonse C.↗
NASA NTRS1980

Comments on the dissipation of hydromagnetic surface waves

A recent paper by Wentzel, which claims to calculate a plasma heating rate due to dissipation of surface waves in an ideal magnetohydrodynamic (MHD) fluid, is found to be in error in interpretation. A well-established general theorem pertaining to the conservative ideal MHD fluid requires that the normal mode calculated by Wentzel be oscillatory in time. Within ideal MHD, dissipation and plasma heating are therefore impossible.

Lee, M. A.↗
NASA NTRS1995

Numerical magetohydrodynamics in astronphysics: Algorithm and tests for one-dimensional flow`

We describe a numerical code to solve the equations for ideal magnetohydrodynamics (MHD). It is based on an explicit finite difference scheme on an Eulerian grid, called the total variation diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. We also describe a non-linear Riemann solver for ideal MHD, which includes rarefractions as well as shocks. The numerical code and the Riemann solver have been used to test each other. Extensive tests encompassing all the possible ideal MHD structures with planar symmetries (i.e., one-dimensional flows) are presented. These include those for which the field structure is two dimensional (i.e., those flows often called '1 + 1/2 dimensional') as well as those for which the magnetic field plane rotates (i.e., those flows often called '1 + 1/2 + 1/2 dimensional'). Results indicate that the code can resolve strong fast, slow, and magnetosonic shocks within two to four cells, but more cells are required if shocks become weak. With proper steepening, we could resolve rotational discontinuities within three to five cells. However, without successful implementation of steepening, contact discontinuities are resolved with approximately 10 cells and tangential discountinuities are resolved with approximately 15 cells. Out tests confirm that slow compound structures with tow-dimensional magnetic fields are composed of intermediate shocks (so-called 2-4 intermediate shocks) followed by slow rarefaction waves. Finally, tests demostrate that in two-dimensional magnetohydrodynamics, fast compound structures, which are composed of intermediate shocks (so-called 1-3 intermediate shocks) preceeded by fast rarefaction waves, are also possible.

Ryu, Dongsu↗
NASA NTRS2007

Broken Symmetry and Coherent Structure in MHD Turbulence

Absolute equilibrium ensemble theory for ideal homogeneous magnetohydrodynamic (MHD) turbulence is fairly well developed. Theory and Simulation indicate that ideal MHD turbulence non-ergodic and contains coherent structure. The question of applicability real (i.e., dissipative) MHD turbulence is examined. Results from several very long time numerical simulations on a 64(exp 3) grid are presented. It is seen that coherent structure begins to form before decay dominates over nonlinearity. The connection with inverse spectral cascades and selective decay will also be discussed.

Shebalin, John V.↗
NASA NTRS1994

Numerical simulation of torus-driven plasma transport in the Jovian magnetosphere

The Rice convection model has been modified for application to the transport of Io-generated plasma through the Jovian magnetosphere. The new code, called the RCM-J, has been used for several ideal-magnetohydrodynamic (MHD) numerical simulations to study how interchange instability causes an initially assumed torus configuration to break up. In simulations that start from a realistic torus configuration but include no energetic particles, the torus disintegrates too quickly (approximately 50 hours). By adding an impounding distribution of energetic particles to suppress the interchange instability, resonable lifetimes were obtained. For cases in which impoundment is insufficient to produce ideal-MHD stability, the torus breaks up predominantly into long fingers, unless the initial condition strongly favors some other geometrical form. If the initial torus has more mass on one side of the planet than the other, fingers form predominatly on the heavy side (which we associate with the active sector). Coriolis force bends the fingers to lag corotation. The simulation results are consistent with the idea that the fingers are formed with a longitudinal thickness that is roughly equal to the latitudinal distance over which the invariant density declines at the outer edges of the initial torus. Our calculations give an average longitudinal distance between plasma fingers of about 15 deg which corresponds to 20 to 30 minutes of rotation of the torus. We point to some Voyager and Ulysses data that are consistent with this scale of torus longitudinal irregularity.

Yang, Y. S.↗
NASA NTRS1975

Acceleration of charged particles in oblique MHD shocks

The motion of high energy charged particles in ideal oblique MHD shocks has been studied extensively. The orbits of charged particles can be solved exactly from Lorentz force equation. It is found that a charged particle may cross the shock front many times before it leaves the shock and it is energized during crossings. The most energetic particles are those with initial pitch and phase near the boundary between crossing and reflected particles. The reflected particles cross the shock about twice as many times as the passing particles do. The energy gain increases for larger values of shock strength and speed. If reflection occurs, the energy gain is larger for smaller angles between the field and shock surface. The introduction of scattering results in a few more higher energy particles.

Chen, G.↗
NASA NTRS1999

On Europa's Magnetospheric Interaction: A MHD Simulation of the E4 Flyby

The global three-dimensional interaction of Europa with the Jovian magnetosphere is modeled by using a complete set of ideal magnetohydrodynamic (MHD) equations. The model accounts for exospheric mass loading, ion-neutral charge exchange, recombination, and a possible intrinsic dipole magnetic field of Europa. The single-fluid MHD equations are solved by using a modem, finite volume, higher-order, Godunov-type method on an adoptively refined unstructured grid, which allows detailed modeling of the region near Europa while still resolving both the upstream region and the satellite's wake. The magnetic field and plasma density measured during Galileo's EGA flyby of December 19, 1996, are reproduced reasonably well in the simulation. We find the agreement between the data and our model particularly convincing if we assume that the plasma velocity during the EGA flyby deviated from the nominal corotation direction by approximately 20 deg. Evidence from the Galileo energetic particle detector also supports this assumption. In this case, we can fit the data using a dipole with orientation close to that of an induced dipole arising from the interaction of a hypothetical conducting subsurface layer on Europa with the periodically changing magnetic field of Jupiter. However, the magnitude of the dipole in our model is somewhat smaller (70%) than that suggested by Khurana et al. The total mass loading and ion-neutral charge exchange rates are consistent with the estimates of Europa's atmosphere and ionosphere.

Kabin, K.↗
NASA NTRS2000

Two-Species, 3D, MHD Simulation of Europa's Interaction with Jupiter's Magnetosphere

The interaction of Europa with the Jovian a magnetosphere has been studied by using a two species in ideal magnetohydrodynamic (MHD) numerical model. This model considers the upstream plasma in the Jovian magnetosphere and the molecular oxygen ions in the ionosphere of Europa, separately. We present results a from simulation studies, which take into account impact ionization, recombination, and the effect of a possible induced dipole magnetic field of Europa. The total mass loading of the magnetospheric flow and the ionization frequency used in the model are consistent with the estimates of Europa's ionosphere and atmosphere. The multi-species MHD equations are solved by using a finite volume, high-order, Godunov-type method on an adoptively refined unstructured grid, which allows a detailed modeling of the region near Europa's surface, while still resolving both the upstream region and the satellite's wake. We have paid special attention to the wake of Europa, in order to be able to make comparisons with the Galileo's E4 flyby observations, as well as other model calculations. The calculated escape flux of a O2+ down the tail was found to be about 5.6 x 10(exp 25) s(sup -1).

Liu, Yifan↗
NASA NTRS1990

Magnetohydrodynamic modeling of the solar corona

The ideal and resistive magnetohydrodynamic (MHD) model is used to examine the dynamics and structure of the solar corona. When the coronal magnetic field is deformed by photospheric flow it can evolve to states that become unstable to ideal MHD modes. The nonlinear evolution of these instabilities can lead to the generation of current sheets, field line reconnection, and energy release. The disruption of an arcade field and the kinking of coronal loops is described. The braiding of the large-scale coronal field by convective photospheric motions develops fine-scale structure in the magnetic field and leads to the development of intense current filaments. The resistive dissipation of these currents can provide an efficient coronal heating mechanism.

Mikic, Zoran↗
NASA NTRS1997

A Self-Consistent Numerical Magnetohydrodynamic (MHD) Model of Helmet Streamer and Flux-Rope Interactions: Initiation and Propagation of Coronal Mass Ejections (CMEs)

We present results for an investigation of the interaction of a helmet streamer arcade and a helical flux-rope emerging from the sub-photosphere. These results are obtained by using a three-dimensional axisymmetric, time-dependent ideal magnetohydrodynamic (MHD) model. Because of the physical nature of the flux-rope, we investigate two types of flux-ropes; (1) high density flux-rope (i.e. flux-rope without cavity), and (2) low density flux rope (i.e. flux-rope with cavity). When the streamer is disrupted by the flux-rope, it will evolve into a configuration resembling the typical observed loop-like Coronal Mass Ejection (CMES) for both cases. The streamer-flux rope system with cavity is easier to be disrupted and the propagation speed of the CME is faster than the streamer-flux rope system without cavity. Our results demonstrate that magnetic buoyancy force plays an important role in disrupting the streamer.

Wu, S. T.↗
NASA NTRS1997

Self-Consistent and Time-Dependent Solar Wind Models

We describe the first results from a self-consistent study of Alfven waves for the time-dependent, single-fluid magnetohydrodynamic (MHD) solar wind equations, using a modified version of the ZEUS MHD code. The wind models we examine are radially symmetrical and magnetized; the initial outflow is described by the standard Parker wind solution. Our study focuses on the effects of Alfven waves on the outflow and is based on solving the full set of the ideal nonlinear MHD equations. In contrast to previous studies, no assumptions regarding wave linearity, wave damping, and wave-flow interaction are made; thus, the models naturally account for the back-reaction of the wind on the waves, as well as for the nonlinear interaction between different types of MHD waves. Our results clearly demonstrate when momentum deposition by Alfven waves in the solar wind can be sufficient to explain the origin of fast streams in solar coronal holes; we discuss the range of wave amplitudes required to obtained such fast stream solutions.

Ong, K. K.↗
NASA NTRS1982

A mechanism for a class of solar coronal disturbances

A new ideal magnetohydrodynamic (MHD) model is used to examine the dynamical response of the upper solar atmosphere to the injection of cold mass from the photosphere, a surge perturbation. Theoretical results show that mass ejections from the photosphere will form loop structures with internally-moving plasma blobs, which is due to the mass injections from both legs in active regions leading to a reflection of waves at the top of the loop. The mass-loaded loop will remain as long as ejection continues, and will excite MHD waves propagating outward to the upper atmosphere, which may trigger a class of coronal disturbances.

Wu, S. T.↗