Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magnetohydrodynamical simulations”

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 415 records · Page 23

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

Studies of plasma flow past Jupiter's satellite Io

We have investigated the interaction of Io, Jupiter's innermost Galilean satellite, with the Io plasma torus. The interaction of Io with the plasma surrounding it has been a subject of interest for almost 30 years, dating from the discovery by Bigg (1964) that radio emissions from the Jovian magnetosphere are controlled by Io's position. Since that time, both ground-based and spacecraft observations have shown that Io is a unique satellite that influences the Jovian magnetosphere in important ways. In particular, material from Io is a major source of plasma for the magnetosphere, and the energy that this plasma harnesses from Jupiter's co-rotating magnetic field is an important power source for the magnetosphere. It is apparent that the local interaction of the torus plasma with Io plays a key role in the formation, composition, and energetics of the Io torus; the interaction is also highly nonlinear. We have modeled this interaction using time-dependent three-dimensional magnetohydrodynamic (MHD) simulations. During this past year, we have used NASA support to develop a new MHD code to study the interaction. As part of the Galileo spacecraft's recent successful insertion into orbit around Jupiter, the spacecraft passed within 900 km of Io's surface. Our calculations have focused on using Galileo particles and fields data to examine a question that was not resolved by the Voyager observations: Does Io have an intrinsic magnetic field? In this progress summary, we describe our efforts on this problem to date.

Linker, Jon↗

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

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

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↗

Magnetotail Structure and its Internal Particle Dynamics During Northward IMF

This study uses Global magnetohydrodynamic (MHD) simulations driven by solar wind data along with Geotail observations of the magnetotail to investigate the magnetotail's response to changes in the interplanetary magnetic field (IMF); observed events used in the study occurred on March 29, 1993 and February 9, 1995. For events from February 9, 1995, we also use the time-dependent MHD magnetic and electric fields and the large-scale kinetic (LSK) technique to examine changes in the Geotail ion velocity distributions. Our MHD simulation shows that on March 29, 1993, during a long period of steady northward IMF, the tail was strongly squeezed and twisted around the Sun-Earth axis in response to variations in the IMF B(sub y) component. The mixed (magnetotail and magnetosheath) plasma observed by Geotail results from the spacecraft's close proximity to the magnetopause and its frequent crossings of this boundary. In our second example (February 9, 1995) the IMF was also steady and northward, and in addition had a significant B(sub y) component. Again the magnetotail was twisted, but not as strongly as on March 29, 1993. The Geotail spacecraft, located approximately 30 R(sub E) downtail, observed highly structured ion distribution functions. Using the time-dependent LSK technique, we investigate the ion sources and acceleration mechanisms affecting the Geotail distribution functions during this interval. At 1325 UT most ions are found to enter the magnetosphere on the dusk side earthward of Geotail with a secondary source on the dawn side in the low latitude boundary layer (LLBL). A small percentage come from the ionosphere. By 1347 UT the majority of the ions come from the dawn side LLBL. The distribution functions measured during the later time interval are much warmer, mainly because particles reaching the spacecraft from the dawnside are affected by nonadiabatic scattering and acceleration in the neutral sheet.

Ashour-Abdalia, M.↗

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

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

Magnetotail Structure and its Internal Particle Dynamics During Northward IMF

This study uses Global magnetohydrodynamic (MHD) simulations driven by solar wind data along with Geotail observations of the magnetotail to investigate the magnetotail's response to changes in the interplanetary magnetic field (IMF); observed events used in the study occurred on March 29, 1993 and February 9, 1995. For events from February 9, 1995, we also use the time-dependent MHD magnetic and electric fields and the large-scale kinetic (LSK) technique to examine changes in the Geotail ion velocity distributions. Our MHD simulation shows that on March 29, 1993, during a long period of steady northward IMF, the tail was strongly squeezed and twisted around the Sun-Earth axis in response to variations in the IMF B(sub y) component. The mixed (magnetotail and magnetosheath) plasma observed by Geotail results from the spacecraft's close proximity to the magnetopause and its frequent crossings of this boundary. In our second example (February 9, 1995) the IMF was also steady and northward, and in addition had a significant B(sub y) component. Again the magnetotail was twisted, but not as strongly as on March 29, 1993. The Geotail spacecraft, located approximately 30 R(sub E) downtail, observed highly structured ion distribution functions. Using the time-dependent LSK technique, we investigate the ion sources and acceleration mechanisms affecting the Geotail distribution functions during this interval. At 1325 UT most ions are found to enter the magnetosphere on the dusk side earthward of Geotail with a secondary source on the dawn side in the low latitude boundary layer (LLBL). A small percentage come from the ionosphere. By 1347 UT the majority of the ions come from the dawn side LLBL. The distribution functions measured during the later time interval are much warmer, mainly because particles reaching the spacecraft from the dawn side are affected by nonadiabatic scattering and acceleration in the neutral sheet.

Ashour-Abdalla, M.↗

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

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↗

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

ISTP Solar Maximum Extended Science Program

We have studied the entry of solar energetic particles (SEPS) into the magnetosphere by following particles in the time dependent magnetic and electric field from global magnetohydrodynamic (MHD) simulations of the magnetosphere. The MHD simulations can either be for idealized interplanetary magnetic field (IMF) conditions, or for upstream conditions measured by spacecraft. An important part of the analysis is understanding the response of the magnetosphere to the IMF conditions. In the idealized case, the MHD simulation included a steady interplanetary magnetic field (IMF) B(sub x), velocity and density, while the B(sub y) and B(sub z) components were varied from southward IMF to dawnward and finally to northward IMF. We launched more than ten million protons, as well as about 1 million He-3 ions and a few thousand electrons upstream of the magnetosphere into the solar wind. They were initialized using a kappa distribution, which is a power law distribution with a power law coefficient of 1.5 at high energies. The particles had energies between 0. 1 and 50 MeV. The particles were run in time dependent MHD fields that were advanced in time as the particles moved through the system.

Ashour-Abdalla, Maha↗

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↗

The Formation of Relativistic Jets from Kerr Black Holes

We have performed the first fully three-dimensional general relativistic magnetohydrodynamics (GRMHD) simulation for Schwarzschild and Kerr black holes with a free falling corona and thin accretion disk. The initial simulation results with a Schwarzschild metric show that a jet is created as in the previous axisymmetric simulations with mirror symmetry at the equator. However, the time to form the jet is slightly longer than in the 2-D axisymmetric simulation. We expect that the dynamics of jet formation are modified due to the additional freedom in the azimuth dimension without axisymmetry with respect to the Z axis and reflection symmetry respect to the equatorial plane. The jet which is initially formed due to the twisted magnetic fields and shocks becomes a wind at the later time. The wind flows out with a much wider angle than the initial jet. The twisted magnetic fields at the earlier time were untwisted and less pinched. The accretion disk became thicker than the initial condition. Further simulations with initial perturbations will provide insights for accretion dynamics with instabilities such as magneto-rotational instability (MRI) and accretion-eject instability (AEI). These instabilities may contribute to variabilities observed in microquasars and AGN jets.

Nishikawa, K.-I.↗

Evolutionary Models of Cold, Magnetized, Interstellar Clouds

We modeled the long-term and small-scale evolution of molecular clouds using direct 2D and 3D magnetohydrodynamic (MHD) simulations. This work followed up on previous research by our group under auspices of the ATP in which we studied the energetics of turbulent, magnetized clouds and their internal structure on intermediate scales. Our new work focused on both global and smallscale aspects of the evolution of turbulent, magnetized clouds, and in particular studied the response of turbulent proto-cloud material to passage through the Galactic spiral potential, and the dynamical collapse of turbulent, magnetized (supercritical) clouds into fragments to initiate the formation of a stellar cluster. Technical advances under this program include developing an adaptive-mesh MHD code as a successor to ZEUS (ATHENA) in order to follow cloud fragmentation, developing a shearing-sheet MHD code which includes self-gravity and externally-imposed gravity to follow the evolution of clouds in the Galactic potential, and developing radiative transfer models to evaluate the internal ionization of clumpy clouds exposed to external photoionizing UV and CR radiation. Gammie's work at UIUC focused on the radiative transfer aspects of this program.

Gammie, Charles F.↗