Search NASA⌕ Search

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.

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

MHD modeling of magnetotail instability for localized resistivity

We present results of a three-dimensional magnetohydrodynamic (MHD) simulation of magnetotail evolution initiated by a sudden occurrence or increase of spatially localized resistivity as the major expected concequence of some localized microinstability. Because of the absence of a quantitative model, possible variations of resistivity levels with current density, or the reduction thereof, are not incorporated in the present investigation. The emphasis of the study is on an investigation of the changes to the overall evolution brought about by this localization, in particular, on the disruption and diversion of the cross-tail current and the nonlinear evolution of the magnetotail instability. The immediate consequences of the occurrence of the localized resistance and the resulting electric field are a reduction and diversion of the electric current around the region of high resistivity, associated with an increase of B(sub z) ('dipolarization') at the earthward edge and a decrease of B(sub z) at the tailward edge of this region. These effects, however, are localized and do not involve a reduction of the total cross-tail current and hence do not lead to the global development of a 'substorm current wedge,' which includes not only the reduction of the cross-tail current but also the buildup of a global field-aligned current system of 'regional 1' type (toward the Earth on the dawnside and away on the duskside of the tail). Such signatures develop at a later time, as consequences of a three-dimensional tearing instability, which is triggered by the occurrence of the resistivity. These features are found in combination with plasmoid formation and ejection, quite similar to results of earlier simulations with uniform resistivity. Differences are found in the timescale of the evolution, which tends to be shorter for localized resistivity, and in the propagation of the dipolarization effects in the equatorial plane. Whereas for uniform resistivity the temporal increase in B(sub z) tends to propagate tailward, apparently due to a pileup effect in the near-Earth region, an earthward propagation is found for the localized resistivity. This propagation results from an earthward convection of the gradient in B(sub z), which is set up by the reconnection process further tailward.

Hesse, Michael↗
NASA NTRS1994

Tearing mode instability in a multiple current sheet system

The tearing mode and magnetic reconnection are studied for multiple current sheet systems by two-dimensional magnetohydrodynamic (MHD) simulations. Both the linear and nonlinear evolution of this process are anaylsed for laminar perturbations. The results illustrate the existence of a linear regime with a symmetric and antisymmetric mode and agree with previous analytic results (Otto and Birk, 1992). The nonlinear evolution shows a number of interesting new features and may explain some properties in corresponding studies of turbulent reconnection. For wavelengths larger than twice the current sheet separation the evolution of antisymmetric modes leads to an entire reconfiguration of the magnetic field and converts a major portion of the magnetic energy into kinetic energy. Antisymmetric modes with smaller wavelengths and symmetric modes are found to saturate. The influence of the value of the resistivity on the reconnection rate decreases in the nonlinear evolution, and the ratio of current sheet separation to wavelength seems to be of major importance. A comparion of the dynamics of periodic current sheets with the evolution of only two current sheets indicates that some of the results for the periodic system also apply to the evolution of only two interacting current sheets. The results are discussed with respect to observations of large-scale plasma and magnetic field reconfigurations in the magnetosheath and near the Earth's bow shock.

Yan, M.↗
NASA NTRS1994

The population of the magnetosphere by solar winds ions when the interplanetary magnetic field is northward

We have examined some possible entry mechanisms of solar wind ions into the magnetosphere by calculating the trajectories of thousands of non-interacting ions in the magnetic and electric fields from a three dimensional global magnetohydrodynamic (MHD) simulation of the magnetosphere and the magnetosheath, under northward interplanetary magnetic field (IMF) conditions. Particles, launched in the solar wind, entered the magnetosphere and formed the low latitude boundary layer (LLBL), plasma sheet and a region of trapped particles near the Earth. The densities and temperatures we obtained in these regions were realistic, with the exception of trapped particle densities. The dominant entry mechanism was convection into the magnetosphere on reconnecting field lines.

Richard, Robert L.↗
NASA NTRS1995

Magnetic reconnection in the presence of sheared flow and density asymmetry: Applications to the Earth's magnetopause

Classical models of magnetic reconnection consist of a small diffusion region bounded by two symmetric slow shocks, across which the plasma is accelerated. Asymmetries often present in space plasmas are sheared plasma flow and dissimilar plasma densities on the two sides of current sheets. In this paper, we investigate magnetic reconnection in the presence of a shear flow and an asymmetric density across the current sheet using two-dimensional magnetohydrodynamic (MHD) simulations. The results demonstrate that magnetic reconnection can occur only for a plasma flow velocity (in the frame of the X line) which is below the Alfven speed in each inflow region. This limits the velocity of the X line to a certain range for a given flow shear and provides an upper limit to the total velocity shear at which reconnection ceases to operate. Depending on the direction of the flow in the adjacent inflow region, the effects from the sheared flow and from the density asymmetry will compete with or enhance each other in respect to the magnitude and location of the currents which bound the outflow regions. The results are applied to the dayside and flank regions of the magnetosphere. For the dayside region where the magnetosheath flow is slow, the magnetic field transition region is thin and the accelerated flow is earthward of the sharp current layer (magnetopause). At the flanks tailward of the X line, shear flow and density asymmetry effects compete making the magnetic field transition layer broad with the high-speed flow contained within the transition region which explains corresponding observations. At the flanks sunward of the X line, shear flow and density asymmetry effects enhance each other and lead to a strong current sheet on the magnetosheath side of the accelerated flow. The total volume affected by magnetic reconnection is much larger than the steady state region. A large bulge region precedes the steady state region. Qualitatively, the bulge and the steady state region have similar signatures and both can explain observations. We provide criteria in order to distinguish between the bulge and the steady state region in observations.

La Belle-Hamer, A. L.↗
NASA NTRS1995

Flow properties along field lines in a 3-D tilted-dipole geometry

A 3-D MHD simulation of a global, tilted-dipole solar wind flow pattern is analyzed to determine flow properties along individual magnetic field lines. In the model, flow conditions near the Sun are chosen to provide a reasonable match to the interplanetary configuration prevailing during the recent south polar passage by Ulysses, i.e., a streamer belt inclined approximately 30 deg to the solar equator and speeds ranging from 325-800 km/s. Field lines all across the stream pattern are traced from 1 to 10 AU by following the motion of marker particles embedded in the flow. It is found that those field lines threading the core of the interaction region are subject to significant latitudinal and relative longitudinal displacement over this range of heliocentric distance. Thus, observations taken at a fixed latitude in the inner solar system sample, over the course of a solar rotation, field lines which connect to a range of latitudes in the outer heliosphere. Maps of the field line displacements are presented to help visualize these connections. In addition, it is found that depending upon the location relative to the CIR structure, the radial evolution of fluid elements frozen to different field lines can deviate considerably from that of the canonical solar wind. That is, for selected subsets of field lines, large speed changes (not just at shocks) can be experienced; the density variation can be far from 1/r(exp 2), and the magnetic field intensity need not decay monotonically with distance.

Pizzo, V. J.↗
NASA NTRS1995

CME stimulated by eruptive prominence

A model of CME arising due to drift motion in the corona in the presence of an eruptive prominence is presented. Magnetic field configuration is in accordance with a model of inverse polarity. In a region where a magnetic pressure of the filament magnetic field is greater than a gas pressure, plasma motion can be assumed as a drift motion. Its characteristic is such that the further one gets from the filament, the higher is drift velocity of a plasma. That sort of motion leads to a rarefaction of plasma and formation of a cavity around the filament. But a current strength in eruptive prominences estimated from observations is such that the region b is less than 1 has a limited size. Near the boundary b = 1 plasma deceleration is occurred and as a result of it coronal density is increasing. Plasma condensation near the surface b = 1 leads to formation of a dense envelope which can be collated with an outer loop of CME. Two dimensional numerical MHD simulation displays a process of cavity and loop formation. If a current is large enough, two compact regions of compressed dense matter arise at both sides of the rising filament and two narrow jets are developed. This scenario, perhaps, corresponds to CMEs in which a top of the loop is faint or it is absent at all.

Filippov, B. P.↗
NASA NTRS1995

Coronal mass ejections and the evolution of the large-scale corona

Coronal mass ejections (CMEs) are dynamic, large-scale events in the solar corona that expel plasma and magnetic fields into the solar wind. The structure of the large-scale corona influences CME occurrence; for example, CMEs are most frequently observed to be disruptions of coronal or helmet streamer configurations on the Sun. We have investigated the evolution of the large-scale corona in the presence of differential rotation, using time-dependent magnetohydrodynamic (MHD) simulations in two and three dimensions. An equilibrium coronal configuration is first computed by superimposing a Parker solar wind flow on an initial potential field, the equations are integrated in time until a steady-state is reached. When differential rotation is introduced, we find that the configuration does not reach a steady state; the closed field regions (helmet streamers) disrupt recurrently, ejecting plasmoids into the solar wind. Our results suggest that differential rotation may be one of the mechanisms by which mass ejections are initiated. We will describe the evolution of the polarization brightness (pB) for the resulting configurations, as well as the propagation of the ejected plasmoids in the solar wind, and we will compare the effects of differential rotation with other possible initiation mechanisms, such as the emergence of new magnetic flux.

Linker, J. A.↗
NASA NTRS1995

Scientific visualization tools for the ISTP project: Mission planning, data analysis and model interpretation

Visualization tools are being developed to meet the challenges of mission planning and data analysis presented by the International Solar-Terrestrial Physics (ISTP) program. ISTP encompasses a large number of spacecraft, multiple ground-based observatories, and several theoretical investigations, with the goal of understanding the global behavior of the solar wind/magnetosphere/ionosphere system. The tools include three-dimensional displays of key boundaries in geospace along with spacecraft trajectories, which can be animated and synchronized to universal time. Magnetic field models and MHD simulation results can be invoked to reveal the magnetic topology or to identify magnetic conjunctions between spacecraft and/or ground-based facilities. Simultaneous displays of satellite trajectories, spacecraft-borne observations, and model predictions are available to facilitate data processing and interpretation efforts. The current status of these tools is described, and their implementation at the ISTP Science Planning and Operations Facility and distribution to the entire ISTP community are discussed.

Peredo, M.↗
NASA NTRS1997

Studies of Plasma Flow Past Jupiters Satellite Io

We have investigated the interaction of Io, Jupiter's innermost Galilean satellite, with the Io plasma torus, and the interaction of Ganymede with the corotating Jovian plasma. With the successful insertion of the Galileo spacecraft into orbit around Jupiter, many new observations have been made of the Jovian magnetosphere. Some of the most exciting results thus far have been in regards to Jupiter's satellites, Io and Ganymede. In both cases the large perturbations to the background (Jovian) magnetic field have been consistent with the satellites' possession of an intrinsic magnetic field. The gravity measurements implying a differentiated core at both Io and Ganymede makes internal generation of a magnetic field by dynamo action in these satellites plausible, and, in the case of Ganymede, the identification of an intrinsic field is apparently unambiguous. For Io the situation is less clear, and further analysis is necessary to answer this important question. During the past year, we have used time-dependent three-dimensional magnetohydrodynamic (MHD) simulations to study these plasma-moon interactions. The results from these simulations have been used directly in the analysis of the Galileo magnetometer data. Our primary emphasis has been on the Io interaction, but we recently presented results on the Ganymede interaction as well. In this progress summary we describe our efforts on these problems to date.

Linker, Jon A.↗
NASA NTRS1996

Latitudinal Dependence of the Radial IMF Component: Coronal Imprint

Measurements by Ulysses have confirmed that there is no significant gradient with respect to heliomagnetic latitude in the radial component, B(sub r,) of the interplanetary magnetic field. In the corona, the plasma, beta is much less than 1, except directly above streamers, so longitudinal and latitudinal gradients in field strength will relax due to the transverse magnetic pressure gradient force as the solar wind carries magnetic flux away from the Sun. This happens quickly enough so that the field is essentially uniform by 5 - 10 solar radius, apparently remaining so as it is carried to beyond 1 AU. Here, we illustrate the coronal relaxation with a qualitative physical argument and by reference to a detailed Magneto HydroDynamics (MHD) simulation.

Suess, S. T.↗
NASA NTRS1997

Cross-field Current Instability for Substorm Expansions

The funding provided by the above-referenced NASA grant has enabled us: (1) to investigate the quasi-linear evolution of the IWI [Lui et al., 1993] and that of the generalized MTSI/IWI [Yoon and Lui, 1993], (2) to carry out the linear analysis of the LHDI to elucidate the difference between it and the MTSI/PM instability [Yoon et al., 1994], (3) to conduct some preliminary nonlocal analyses of the MTSI [Lui et al., 1995] and the IWI [Yoon and Lui, 1996] modes, (4) to study low-frequency shear-driven instability and its nonlinear evolution, which might compete with the CCI [Yoon et al., 1996], and (5) to study the evolution of current sheet during late substorm growth phase by means of 2-D Hall-MHD simulation in order to obtain a better understanding of the current sheet equilibrium crucial for CCI theory [Yoon and Lui, 1997].

Lui, Anthony↗
NASA NTRS1998

Theoretical Investigation of the High-Altitude Cusp Region using Observations from Interball and ISTP Spacecraft

A fundamental goal of magnetospheric physics is to understand the transport of plasma through the solar wind-magnetosphere-ionosphere system. To attain such an understanding, we must determine the sources of the plasma, the trajectories of the particles through the magnetospheric electric and magnetic fields to the point of observation, and the acceleration processes they undergo enroute. This study employed plasma distributions observed in the near-Earth plasma sheet by Interball and Geotail spacecraft together with theoretical techniques to investigate the ion sources and the transport of plasma. We used ion trajectory calculations in magnetic and electric fields from a global Magnetohydrodynamics (MHD) simulation to investigate the transport and to identify common ion sources for ions observed in the near-Earth magnetotail by the Interball and Geotail spacecraft. Our first step was to examine a number of distribution functions and identify distinct boundaries in both configuration and phase space that are indicative of different plasma sources and transport mechanisms. We examined events from October 26, 1995, November 29-30, 1996, and December 22, 1996. During the first event Interball and Geotail were separated by approximately 10 R(sub E) in z, and during the second event the spacecraft were separated by approximately 4(sub RE). Both of these events had a strong IMF By component pointing toward the dawnside. On October 26, 1995, the IMF B(sub Z) component was northward, and on November 1-9-30, 1996, the IMF B sub Z) component was near 0. During the first event, Geotail was located near the equator on the dawn flank, while Interball was for the most part in the lobe region. The distribution function from the Coral instrument on Interball showed less structure and resembled a drifting Maxwellian. The observed distribution on Geotail, on the other hand, included a great number of structures at both low and high energies. During the third event (December 22, 1996) both spacecraft were in the plasma sheet and were separated bY approximately 20 R(sub E) in the y direction. During this event the IMF was southward.

Ashour-Abdalla, Maha↗
NASA NTRS1997

Ion Sources and Acceleration Mechanisms Inferred from Local Distribution Functions

This study investigates the sources of the ions up the complex and nonisotropic H(+) velocity distribution functions observed by the Geotail spacecraft on May 23, 1995, in the near-Earth magnetotail region and recently reported by Frank et al. [1996]. A distribution function observed by Geotail at -10 R(sub E) downtail is used as input for the large scale kinetic (LSK) technique to follow the trajectories of approximately 90,000 H(+) ions backward in time. Time-dependent magnetic and electric fields are taken from a global magnetohydrodynamic (MHD) simulation of the magnetosphere and its interactions with appropriate solar wind and IMF conditions. The ion population described by the Geotail distribution function was found to consist of a mixture of particles originating from three distinct sources: the ionosphere, the low latitude boundary layer (LLBL), and the high latitude plasma mantle. Ionospheric particles had direct access along field lines to Geotail, and LLBL ions convected adiabatically to the Geotail location. Plasma mantle ions, on the other hand, exhibited two distinct types of behavior. Most near-Earth mantle ions reached Geotail on adiabatic orbits, while distant mantle ions interacted with the current sheet tailward of Geotail and had mostly nonadiabatic orbits. Ions from the ionosphere, the LLBL, and the near-Earth mantle were directly responsible for the well-separated, low energy structures easily discernible in the observed and modeled distribution functions. Distant mantle ions formed the higher energy portion of the Geotail distribution. Thus, we have been successful in extracting useful information about particle sources, their relative contribution to the measured distribution and the acceleration processes that affected particle transport during this time.

Ashour-Abdalla, M.↗
NASA NTRS1998

A Study of Transport in the Near-Earth Plasma Sheet During A Substorm Using Time-Dependent Large Scale Kinetics

In this study we investigate the transport of H+ ions that made up the complex ion distribution function observed by the Geotail spacecraft at 0740 UT on November 24, 1996. This ion distribution function, observed by Geotail at approximately 20 R(sub E) downtail, was used to initialize a time-dependent large-scale kinetic (LSK) calculation of the trajectories of 75,000 ions forward in time. Time-dependent magnetic and electric fields were obtained from a global magnetohydrodynamic (MHD) simulation of the magnetosphere and its interaction with the solar wind and the interplanetary magnetic field (IMF) as observed during the interval of the observation of the distribution function. Our calculations indicate that the particles observed by Geotail were scattered across the equatorial plane by the multiple interactions with the current sheet and then convected sunward. They were energized by the dawn-dusk electric field during their transport from Geotail location and ultimately were lost at the ionospheric boundary or into the magnetopause.

El-Alaoui, M.↗
NASA NTRS1999

The Interaction of Io and Europa With Jupiter's Magnetosphere

The Galileo orbiter has provided radio occultation measurements of the electron density profiles of the plasma surrounding Io and Europa. There have been six occultations of Io, providing twelve electron density profiles at various locations relative to the ram direction of the impinging particles of the Jupiter magnetosphere on Io, and eight profiles on Europa. The two satellites were found to have very different plasma environments, with Io having a proper ionosphere produced on top of an endogenous SO2 atmosphere by magnetospheric particle precipitation and solar EUV, while Europa has a tenuous plasma environment produced by the same mechanisms from an oxygen atmosphere itself also produced by sputtering of water ice from its surface by impinging magnetospheric particles. In both cases the observed electron density profiles are highly asymmetrical, with a compressed profile on the ram side, and an extended one on the wake side. The presence of several measurements for each satellite at different ram-to-wake directions provided data for estimating an approximate distribution of ionization from the ram direction to the wake direction, which were compared with the results of MHD simulations.

Kliore, Arvydas J.↗
NASA NTRS2000

Modeling Ultraviolet Emissions Near Io

In this report, we describe work awarded to Science Applications International Corporation, for the period 6/l/99 to 5/31/00. During this time period, we have investigated the interaction of Io, Jupiter's innermost Galilean satellite, with the Io plasma torus, and the role this interaction plays in producing ultraviolet (UV) emissions from neutral oxygen and sulfur. Io, the innermost of Jupiter's Galilean satellites, plays a unique role in the jovian magnetosphere. Neutral material that escapes from Io is ionized to form the lo torus, a dense, heavy-ion plasma that corotates with Jupiter and interacts with Io. Io supplies not only the torus, but is a major source of plasma for the entire magnetosphere. Ionization and charge-exchange of neutrals near lo strongly influences the plasma interaction, and Io's neutral atmosphere plays an important role in the generation of currents that couple Io to Jupiter. There have been no in situ measurements of the neutral density near Io, but remote observations of neutrals near lo have been performed for many years. Recent observations from the Hubble Space Telescope (HST) have shown detailed structure in UV emissions from neutral species near Io. Electron-impact of the neutrals by the Io torus plasma is the primary mechanism responsible for exciting these emissions. Previously, we have modeled the Io plasma environment using three-dimensional magnetohydrodynamic (MHD) simulations, and we have shown that the interaction between Io and the plasma torus plays an important role in producing the morphology of the observed emissions. In the past year, we have extended these studies to use both UV observations and Galileo particle and field measurements to investigate the Io interaction.

Linker, Jon A.↗
NASA NTRS2001

Properties and Radial Trends of Coronal Mass Ejecta and Their Associated Shocks Observed by Ulysses in the Ecliptic Plane

In this paper, magnetic and plasma measurements are used to analyze 17 interplanetary coronal mass ejections (CMEs) identified by Ulysses during its in-ecliptic passage to Jupiter. We focus on the expansion characteristics of these CMEs (as inferred from the time rate of change of the velocity profiles through the CMEs) and the properties of 14 forward shocks unambiguously associated with these CMEs. We highlight radial trends from 1 to 5.4 AU. Our results indicate that the CMEs are generally expanding at all heliocentric distances. With regard to the shocks preceding these ejecta, we note the following: (1) There is a clear tendency for the shock speed (in the upstream frame of reference) to decrease with increasing heliocentric distance as the CMEs transfer momentum to the ambient solar wind and slow down; (2) 86% of the shock fronts are oriented in the ecliptic plane such that their normals point westward (i.e., in the direction of planetary motion about the Sun), (3) 86% of the shocks are propagating toward the heliographic equator; and (4) no clear trend was found in the strength of the shocks versus heliocentric distance. These results are interpreted using simple dynamical arguments and are supported by fluid and magnetohydrodynamic (MHD) simulations.

Riley, Pete↗
NASA NTRS2001

Global Magnetohydrodynamic Modeling of the Solar Corona

This report describes the progress made in the investigation of the solar corona using magnetohydrodynamic (MHD) simulations. Coronal mass ejections (CME) are believed to be the primary cause of nonrecurrent geomagnetic storms and these have been investigated through the use of three-dimensional computer simulation.

Linker, Jon A.↗