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 163 records · Page 9

A catastrophe mechanism for coronal mass ejections

The ideal-MHD equations are used to show that a coronal current filament can suddenly lose equilibrium if its magnetic energy exceeds a critical value. The loss of equilibrium in the configuration results from an imbalance between magnetic tension and compression, and this imbalance ejects the filament upwards. Near the critical value, the equilibrium configuration develops a vertical current sheet attached to the photosphere at the point directly below the filament. When equilibrium is lost, field lines anchored to the photosphere are stretched upwards, and the current sheet rapidly grows longer. Without reconnection in the current sheet, the filament travels only a short distance before reaching a new equilibrium, and the net magnetic energy released is less than 1 percent of the stored magnetic energy. However, with reconnection, the filament travels upwards indefinitely, and all of the stored energy is released.

Forbes, T. G.↗

Magnetically driven jets and winds

Four equations for the origin and propagation of nonrelativistic jets and winds are derived from the basic conservation laws of ideal MHD. The axial current density is negative in the vicinity of the axis and positive at larger radii; there is no net current because this is energetically favored. The magnetic field is essential for the jet solutions in that the zz-component of the magnetic stress acts, in opposition to gravity, to drive matter through the slow magnetosonic critical point. For a representative self-consistent disk/jet solution relevant to a protostellar system, the reaction of the accreted mass expelled in the jets is 0.1, the ratio of the power carried by the jets to the disk luminosity is 0.66, and the ratio of the boundary layer to disk luminosities is less than about 0.13. The star's rotation rate decreases with time even for rotation rates much less than the breakup rate.

Lovelace, R. V. E.↗

On steady field-aligned double-adiabatic flow

The paper summarizes the four general double-adiabatic invariants governing ideal MHD flow in a narrow tube and renders them in a form suitable for the calculation of flow in tubes of nonconstant cross-sectional area. Results for converging-diverging flux tubes are developed and compared to single-adiabatic results which are identical to those for ordinary gasdynamic flow. It is found that in certain parameter regimes the governing equation for the vector potential is elliptic and can be reduced to Laplace's equation by simple stretching of one coordinate. In other regimes the equation is hyperbolic, indicating the presence of wave patterns. The relevance of these results to flow over magnetopause undulations, including flux transfer events, is discussed.

Sonnerup, B. U. O.↗

Preliminary investigation of power flow and electrode phenomena in a multi-megawatt coaxial plasma thruster

The present report on preliminary results of theoretical and experimental investigations of power flow in a large, unoptimized, multimegawatt coaxial thruster evaluates the significance of these data for the development of efficient, megawatt-class magnetoplasmadynamic (MPD) thrusters. The good agreement obtained between thruster operational performance and model predictions suggests that ideal MHD processes, including those of a magnetic nozzle, play an important role in coaxial plasma thruster dynamics at power levels relevant to advanced space propulsion. An optimized magnetic nozzle design would aid the development of efficient, multimegawatt MPD thrusters.

Schoenberg, Kurt F.↗

Is the earth's magnetotail balloon unstable?

In the past, the onset of magnetospheric substorms has been attributed to the plasma tearing mode instabilities. This paper investigates the ideal MHD ballooning instability of the near-and middle-tail magnetosphere region, as a first step toward determining whether it could trigger the tearing mode, by using the energy principle to investigate whether standard 2D tail models with the 'hard' ionospheric boundary condition are unstable to ballooning instability. Numerical results are presented for compressible ballooning modes that are symmetric about the center of the current sheet. It is shown that, for such a hard boundary condition, no reasonable magnetotail configuration exists that would be unstable to compressible symmetric ballooning but stable against interchange.

Lee, D.-Y.↗

Dynamics of aging magnetic clouds

The dynamics of radially expanding magnetic clouds is rigorously analyzed within the framework of ideal MHD. The cloud is modelled as a cylindrically symmetric magnetic flux rope. In the force balance we include the gas pressure gradient and the Lorentz force. Interaction with the ambient solar wind due to expansion of the magnetic cloud is represented by a drag force proportional to the bulk velocity. We consider the self-similar expansion of a polytrope, and reduce the problem to an ordinary nonlinear differential equation for the evolution function. Analyzing the asymptotic behavior of the evolution function, we formulate theoretical expectations for the long-term behavior of cloud parameters. We focus on the temporal evolution of (1) the magnetic field strength; (2) the twist of the field lines; (3) the asymmetry of the total field profile; and (4) the bulk flow speed. We present data from two magnetic clouds observed at 1 AU and 2 AU, respectively, and find good agreement with theoretical expectations. For a peak magnetic field strength at 1 AU of 25 nT and a polytropic index of 0.5, we find that a magnetic cloud can be distinguished from the background interplanetary field up to a distance of about 5 AU. Taking larger magnetic fields and bigger polytropic indices this distance can double.

Osherovich, V. A.↗

Magnetically driven relativistic jets and winds: Exact solutions

We present self-consistent solutions of the full set of ideal MHD equations which describe steady-state relativistic cold outflows from thin accretion disks. The magnetic field forms a spiral which is anchored in the disk, rotates with it, and accelerates the flow out of the disk plane. The collimation at large distances depends on the total amount of electric current that flows along the jet. We considered various distributions of electric current and derived the result that in straight jets which extend to infinite distances, a strong electric current flows along their axis of symmetry. The asymptotic flow velocities are of the order of the initial rotational velocity at the base of the flow (a few tenths of the speed of light). The solutions are applied to both galactic (small-scale) and extragalactic (large-scale) jets.

Contopoulos, J.↗

On current sheet approximations in models of eruptive flares

We consider an approximation sometimes used for current sheets in flux-rope models of eruptive flares. This approximation is based on a linear expansion of the background field in the vicinity of the current sheet, and it is valid when the length of the current sheet is small compared to the scale length of the coronal magnetic field. However, we find that flux-rope models which use this approximation predict the occurrence of an eruption due to a loss of ideal-MHD equilibrium even when the corresponding exact solution shows that no such eruption occurs. Determination of whether a loss of equilibrium exists can only be obtained by including higher order terms in the expansion of the field or by using the exact solution.

Bungey, T. N.↗

Physics and evolution of constant opening angle jets using a quasi-one-dimensional magnetohydrodynamic model

We discuss the significance of the assumptions of infinite conductivity and time independence in the context of an ideal MHD model for constant opening angle jets. The model is developed by projecting the MHD equations onto the jet axis. We find that for initially sub-Alfvenic flows (i.e., flows emanating from active galactic nuclei and neutron stars) wind-type solutions exist only when the field lines at the origin are wound up in a direction opposite to the direction of rotation of the compact source. We discuss the possibility that the time evolution of these outflows may be a cycle between breeze- and wind-type solutions as a result of continuous changes in the boundary conditions at the origin due to accretion. We propose that such cycles may explain the apparent one-sideness of some jets, especially the ones for which we cannot use arguments of relativistic beaming. We examine the dependence of the wind-type solutions on the following parameters describing the outflow at the origin: the degree of winding of the field lines, the value of the gas pressure, the polytropic index, the strength of the magnetic field, the value of the rotational velocity, the gravitational potential of the compact object, and the injection velocity. We compare results with results obtained previously, and discuss briefly the qualitative features and physical interpretation of the solutions for outflows emanating from neutron stars and protostars.

Koupelis, Theodoros↗

Nonlinear evolution of magnetic flux ropes. 2: Finite beta plasma

In this second paper on the evolution of magnetic flux ropes we study the effects of gas pressure. We assume that the energy transport is described by a polytropic relationship and reduce the set of ideal MHD equations to a single, second-order, nonlinear, ordinary differential equation for the evolution function. For this conservative system we obtain a first integral of motion. To analyze the possible motions, we use a mechanical analogue -- a one-dimensional, nonlinear oscillator. We find that the effective potential for such an oscillator depends on two parameters: the polytropic index gamma and a dimensionless quantity kappa the latter being a function of the plasma beta, the strength of the azimuthal magnetic field relative to the axial field of the flux rope, and gamma. Through a study of this effective potential we classify all possible modes of evolution of the system. In the main body of the paper, we focus on magnetic flux ropes whose field and gas pressure increase steadily towards the symmetry axis. In this case, for gamma greater than 1 and all values of kappa, only oscillations are possible. For gamma less than 1, however, both oscillations and expansion are allowed. For gamma less than 1 and kappa below a critical value, the energy of the nonlinear oscillator determines whether the flux rope will oscillate or expand to infinity. For gamma less than 1 and kappa above critical, however, only expansion occurs. Thus by increasing kappa while keeping gamma fixed (less than 1), a phase transition occurs at kappa = kappa(sub critical) and the oscillatory mode disappears. We illustrate the above theoretical considerations by the example of a flux rope of constant field line twist evolving self-similarly. For this example, we present the full numerical MHD solution. In an appendix to the paper we catalogue all possible evolutions when (1) either the magnetic field or (2) the gas pressure decreases monotonically toward the axis. We find that in these cases critical conditions can occur for gamma greater than 1. While in most cases the flux rope collapses, there are notable exceptions when, for certain ranges of kappa and gamma, collapse may be averted.

Osherovich, V. A.↗

Alfven wave resonances and flow induced by non-linear Alfven waves in a stratified atmosphere

A nonlinear, time-dependent, ideal MHD code has been developed and used to compute the flow induced by nonlinear Alfven waves propagating in an isothermal, stratified, plane-parallel atmosphere. The code is based on characteristic equations solved in a Lagrangian frame and is highly accurate. Results show that resonance behavior of Alfven waves exists in the presence of a continuous density gradient and that the waves with periods corresponding to resonant peaks exert considerably more force on the medium than off-resonance periods; this leads to enhanced flow. If only off-peak periods are considered, the relationship between the wave period and induced longitudinal velocity shows that short period WKB waves push more on the background medium than longer period, non-WKB, waves. The results also show the development of the longitudinal waves produced by the finite amplitude of the Alfven waves. The longitudinal wave becomes strong as the Alfven wave relative amplitude grows above 10 percent and will lead to strong damping of the Alfven waves.

Stark, B. A.↗

Physics of Rotation Powered Pulsars and Their Nebulae

This letter is my progress report for the Astrophysics Theory grant. The first goal of the research supported by this grant is to produce a time dependent theory of the unsteady relativistic collisionless shock wave terminating the relativistic wind from a pulsar, and compare the predicted surface brightness fluctuations to Hubble Space Telescope observations of the wisps in the Crab Nebula. The second goal is to model the production of electron-positron pairs over the polar caps of rotation powered pulsars, and use the results to predict the heating of the surface due to particle trapping and bombardment of the atmosphere at the base of the polar field lines. We have succeeded in creating a one dimensional hybrid code, in which the electron-positron pairs incident on the shock structure are modeled as a relativistic, adiabatic, ideal MHD fluid, while the heavy ions are treated as particles using a particle-in-cell algorithm. The electromagnetic fields are evaluated from the currents and charge densities in the pairs and the ions, while the particles and the fluid accelerate in response to the computed self-consistent electromagnetic fields. The results are promising, in that the underlying ion cyclotron instability generates finite amplitude, propagating magnetosonic waves in the pairs, whose wavelengths and frequencies, when translated into physical units, are comparable to the observed running waves of brightness observed by HST near the Crab pulsar. The code is undergoing a number of tests, to assure us that this preliminary correspondence is not an artifact. In the coming year, the observational appearance of the models will be computed and compared to the HST observations of the Crab now in hand, and used to predict the HST results which will be obtained the year after next. WE also developed a one dimensional cascade theory for pair creation over pulsars' polar caps. A linear integral equation describing the synchrotron cascade has been derived and solved by iterative techniques, in the case when a high energy electron moving parallel to a star centered dipole magnetic field initiates the cascade through curvature gamma ray emission.

Arons, Jonathan↗

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

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

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

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↗

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↗

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↗