Search NASA⌕ Search

SEARCH · Search NASA

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 91 records · Page 5

Geometry of interplanetary magnetic clouds

Two dimensional magnetohydrodynamic simulations are presented of the distortion of a magnetic flux rope that is being accelerated through ambient solar wind plasma. The flux rope magnetic field has an axial component parallel to the solar wind field and an azimuthal component, which lies in the simulation plane. As the flux rope moves through the solar wind plasma, vortices form on its trailing edge and couple strongly to its interior. If the flux rope azimuthal field is weak, it deforms into an elongated banana-like shape a few Alfven transit times. A strong azimuthal field component tends to inhibit this distortion. If the flux rope is taken to model a magnetic cloud, it is suggested that the shape of the cloud at 1 AU is determined by its distortion in the inner solar wind. Distortion timescales beyond 1 AU are estimated as many days. It is estimated that effective drag coefficients somewhat greater than unity are appropriate for modelling flux rope propagation.

Cargill, P. J.↗

Magnetic flux ropes at the high-latitude magnetopause

We examine the consequences of magnetic reconnection at the high-latitude magnetopause using a three-dimensional global magnetohydrodynamic simulation of the solar wind interaction with the Earth's magnetosphere. Magnetic field lines from the simulation reveal the formation of magnetic flux ropes during periods with northward interplanetary magnetic field. These flux ropes result from multiple reconnection processes between the lobes field lines and draped magnetosheath field lines that are convected around the flank of the magnetosphere. The flux ropes identified in the simulation are consistent with features observed in the magnetic field measured by Hawkeye-1 during some high-latitude magnetopause crossings.

Berchem, Jean↗

Forced three-dimensional magnetic reconnection due to linkage of magnetic flux tubes

During periods of southward interplanetary magnetic field (IMF) orientation the magnetic field geometry at the dayside magnetopause is susceptible to magnetic reconnection. It has been suggested that reconnection may occur in a localized manner at several patches on the magnetopause. A major problem with this picture is the interaction of magnetic flux ropes which are generated by different reconnection processes. An individual flux rope is bent elbowlike where it intersects the magnetopause and the magnetic field changes from magnetospheric to interplanetary magnetic field orientation. Multiple patches of reconnection can lead to the formation of interlinked magnetic flux tubes. Although the corresponding flux is connected to the IMF the northward and southward connected branches are hooked into each other and cannot develop independently. We have studied this problem in the framework of three-dimensional magnetohydrodynamic simulations. The results indicate that a singular current sheet forms at the interface of two interlinked flux tubes if no resistivity is present in the simulation. This current sheet is strongly tilted compared to the original current sheet. In the presence of resistivity the interaction of the two flux tubes forces a fast reconnection process which generates helically twisted closed magnetospheric flux. This linkage induced reconnection generates a boundary layer with layers of open and closed magnetospheric flux and may account for the brightening of auroral arcs poleward of the boundary between open and closed magnetic flux.

Otto, A.↗

Substorm associated ion acceleration in the dynamic magnetotail

Ion acceleration and flux increase associated with substorm energetic particle injections are investigated on the basis of geosynchronous observations and test proton orbits in the dynamic fields of a three dimensional magnetohydrodynamic simulation of neutral line formation and dipolarization in the magnetotail. The energetic particle flux changes obtained from the test particle orbits agree well with observations that demonstrate rapid ion flux increases at energies of above 20 keV. The injection region inferred from the test particles has a sharp earthward boundary and a sharp ragged tailward boundary. The earthward portion of the enhanced ion flux can be traced to the enhanced cross-tail electric field associated with the near-earth x-type neutral line. Due to the rapid earthward motion of accelerated ions away from the neutral line, this boundary is displaced earthward to where the energetic ions become more adiabatic in the stronger dipolar field.

Birn, J.↗

Magnetic reconnection during magnetospheric substorms

The near earth reconnection model of substorms represents an attempt to place a broad range of observations into a consistent framework. The roles and requirements of reconnection are discussed. High speed plasma sheet flows, thin current sheet instability, substorm triggering, plasmoids and flux ropes in the distant tail, and magnetohydrodynamic simulations are discussed. Substorms are global, coherent sequences of processes involving solar wind/magnetosphere/ionosphere interaction. Magnetic reconnection is required to explain different dayside and polar cap phenomena, which required nightside reconnection. The modification and expansion of the standard near earth neutral line (NENL) model can integrate breakup arcs, current disruption, current wedge features, and localized plasma flows into the magnetic reconnection framework.

Baker, Daniel N.↗

The Geometric Spreading of Coronal Plumes and Coronal Holes

The geometric spreading in plumes and in the interplume region in coronal holes is calculated, using analytic and numerical theoretical models, between 1.0 and 5.0 solar radius. We apply a two-scale approximation that permits the rapid local spreading at the base of plumes (f(sub t)) to be evaluated separately from the global spreading (f(sub g)) imposed by coronal hole geometry. We show that f(sub t) can be computed from a potential-field model and f(sub g) can be computed from global magnetohydrodynamic simulations of coronal structure. The approximations are valid when the plasma beta is mail with respect to unity and for a plume separation small with respect to a solar radius.

Suess, S. T.↗

The Geometric Spreading of Coronal Plumes and Coronal Holes

The geometric spreading in plumes and in the interplume region in coronal holes is calculated, using analytic and numerical theoretical models, between 1.0 and 5.0 solar radius. We apply two scale approximation that permits the rapid local spreading at the base of plumes to be evaluated separately from the global spreading imposed by coronal hole geometry. We show that fl can be computed from a potential field model and fg can be computed from global magnetohydrodynamic simulations of coronal structure. The approximations are valid when the plasma beta is small with respect to unity and for a plume separation small with respect to a solar radius.

Suess, S. T.↗

Determination of Particle Sources for a Geotail Distribution Function Observed on May 23, 1995

On May 23, 1995, the Comprehensive Plasma Instrumentation (CPI) onboard the Geotail spacecraft observed a complex and structured ion distribution function near the magnetotail midplane at x approximately -10 R(sub E). On the same day, the Wind spacecraft observed a very high density (approximately 40/cubic cm) solar wind and an interplanetary magnetic field (IMF) that was predominantly northward but had several southward turnings. We have inferred the sources of the ions in this distribution function by following approximately 90,000 ion trajectories backward in time using time-dependent electric and magnetic fields obtained from a global MHD (magnetohydrodynamic) simulation. Wind data were used as input for the MHD model. We found that three sources contributed to this distribution: the ionosphere, the plasma mantle which had near-Earth and distant tail components, and the low latitude boundary layer (LLBL). Moreover, distinct structures in the low energy part of the distribution function were found to be associated with individual sources. Structures near 0 deg pitch angle were made up of either ionospheric or plasma mantle ions, while structures near 90 deg pitch angle were dominated by ions from the LLBL source. Particles that underwent nonadiabatic acceleration were numerous in the higher energy part of the ion distribution function, whereas ionospheric and LLBL ions were mostly adiabatic. A large proportion of the near-Earth mantle ions underwent adiabatic acceleration, while most of the distant mantle ions experienced nonadiabatic acceleration.

Ashour-Abdalla, M.↗

Review of Recent Results in Global MHD Modeling: ISTP Project Scientist for Theory and Ground-Based Observations

Global MHD (magnetohydrodynamic) simulations have shown a remarkable ability to describe the global dynamics of geospace. The limitations of the physical approximations underlying MHD would seem to limit the effectiveness of these codes, since kinetic and hybrid effects should manifest themselves by cross-scale coupling from microscales to mesoscales to global scales. However three effects appear to allow the codes to operate much more successfully than one would at first believe. They are:(l) the globally self-consistent nature of the codes with very well defined exterior boundary conditions (the solar wind) which allows the proper intercommunication between magnetospheric regions on MHD scales, (2) the control by global dynamics of the boundary layer locations where micro and meso scale processes operate, and (3) the critical role of numerical diffusion and with a sufficiently high resolution grid, the use of an empirical resistivity term, which if set at a level where the major magnetosphere boundaries properly calibrate against their observed locations, appear to well represent the effects of kinetic and hybrid processes on the global dynamics. The effectiveness of the global MHD codes, which have been developed under the ISTP mission, in describing Wind, Polar and Geotail observations, as well as ground-based observations are described. Particular emphasis is placed upon the Polar imaging data which when combined with ground-based data and global MHD-based synthetic aurora and convection patterns provide a powerful tool in understanding the final link in the solar-terrestrial chain: coupling into the atmosphere and ionosphere.

Curtis, Steven↗

An Investigation of the Large Scale Evolution and Topology of Coronal Mass Ejections in the Solar Wind

This investigation is concerned with the large-scale evolution and topology of coronal mass ejections (CMEs) in the solar wind. During this reporting period we have focused on several aspects of CME properties, their identification and their evolution in the solar wind. The work included both analysis of Ulysses and ACE observations as well as fluid and magnetohydrodynamic simulations. In addition, we analyzed a series of "density holes" observed in the solar wind, that bear many similarities with CMEs. Finally, this work was communicated to the scientific community at three meetings and has led to three scientific papers that are in various stages of review.

Riley, Peter↗

A Global Magnetohydrodynamic Model of Jovian Magnetosphere

The goal of this project was to develop a new global magnetohydrodynamic model of the interaction of the Jovian magnetosphere with the solar wind. Observations from 28 orbits of Jupiter by Galileo along with those from previous spacecraft at Jupiter, Pioneer 10 and 11, Voyager I and 2 and Ulysses, have revealed that the Jovian magnetosphere is a vast, complicated system. The Jovian aurora also has been monitored for several years. Like auroral observations at Earth, these measurements provide us with a global picture of magnetospheric dynamics. Despite this wide range of observations, we have limited quantitative understanding of the Jovian magnetosphere and how it interacts with the solar wind. For the past several years we have been working toward a quantitative understanding of the Jovian magnetosphere and its interaction with the solar wind by employing global magnetohydrodynamic simulations to model the magnetosphere. Our model has been an explicit MHD code (previously used to model the Earth's magnetosphere) to study Jupiter's magnetosphere. We continue to obtain important insights with this code, but it suffers from some severe limitations. In particular with this code we are limited to considering the region outside of 15RJ, with cell sizes of about 1.5R(sub J). The problem arises because of the presence of widely separated time scales throughout the magnetosphere. The numerical stability criterion for explicit MHD codes is the CFL limit and is given by C(sub max)(Delta)t/(Delta)x less than 1 where C(sub max) is the maximum group velocity in a given cell, (Delta)x is the grid spacing and (Delta)t is the time step. If the maximum wave velocity is C(sub w) and the flow speed is C(sub f), C(sub max) = C(sub w) + C(sub f). Near Jupiter the Alfven wave speed becomes very large (it approaches the speed of light at one Jovian radius). Operating with this time step makes the calculation essentially intractable. Therefore under this funding we have been designing a new MHD model that will be able to compute solutions in the wide parameter regime of the Jovian magnetosphere.

Walker, Raymond J.↗

Modeling Magnetospheric Sources

We have used global magnetohydrodynamic, simulations of the interaction between the solar wind and magnetosphere together with single particle trajectory calculations to investigate the sources of plasma entering the magnetosphere. In all of our calculations solar wind plasma primarily enters the magnetosphere when the field line on which it is convecting reconnects. When the interplanetary magnetic field has a northward component the reconnection is in the polar cusp region. In the simulations plasma in the low latitude boundary layer (LLBL) can be on either open or closed field lines. Open field lines occur when the high latitude reconnection occurs in only one cusp. In the MHD calculations the ionosphere does not contribute significantly to the LLBL for northward IMF. The particle trajectory calculations show that ions preferentially enter in the cusp region where they can be accelerated by non-adiabatic motion across the high latitude electric field. For southward IMF in the MHD simulations the plasma in the middle and inner magnetosphere comes from the inner (ionospheric) boundary of the simulation. Solar wind plasma on open field lines is confined to high latitudes and exits the tailward boundary of the simulation without reaching the plasma sheet. The LLBL is populated by both ionospheric and solar wind plasma. When the particle trajectories are included solar wind ions can enter the middle magnetosphere. We have used both the MHD simulations and the particle calculations to estimate source rates for the magnetosphere which are consistent with those inferred from observations.

Walker, Raymond J.↗

Magnetohydrodynamic models of bipolar knotty jet in henize 2-90

A remarkably linear, bipolar, knotty jet was recently discovered in Hen 2-90, an object classified as a young planetary nebula. Using two-dimensional, magnetohydrodynamic simulations, we investigate periodic variations in jet density and velocity as the mechanism for producing the jet and its knotty structures.

bipolar stars winds outflows jet density↗

A Magnetohydrodynamic Boost for Relativistic Jets

We have performed relativistic magnetohydrodynamic simulations of the hydrodynamic boosting mechanism for relativistic jets explored by Aloy & Rezzolla (2006) using the RAISHIN code. Simulation results show that the presence of a magnetic field may change the properties of the shock interface between the tenuous, overpressured jet (V(sub j) (sup z)) flowing tangentially to a dense external medium. Magnetic fields can lead to more efficient acceleration of the jet, in comparison to the pure-hydrodynamic case. A poloidal magnetic field (B(sup z)), tangent to the interface and parallel to the jet flow, produces both a stronger outward moving shock and inward moving rarefaction wave. This leads to a large velocity component normal to the interface in addition to acceleration tangent to the interface, and the jet is thus accelerated to a larger Lorentz factors than those obtained in the pure-hydrodynamic case. In contrast, a strong toroidal magnetic field (B(sup y)), tangent to the interface but perpendicular to the jet flow, also leads to stronger acceleration tangent to the shock interface relative to the pure-hydrodynamic case, but to a lesser extent than found for the poloidal case due to the fact that the velocity component normal to the shock interface is now much smaller. Overall, the acceleration efficiency in the toroidal case is less than that of the poloidal case but both geometries still result in higher Lorentz factors than the pure-hydrodynamic case. Thus, the presence and relative orientation of a magnetic field in relativistic jets can have a significant influence on the hydrodynamic boost mechanism studied by Aloy & Rezzolla (2006).

Mizuno, Yosuke↗

A Magnetohydrodynamic Boost for Relativistic Jets

We performed relativistic magnetohydrodynamic simulations of the hydrodynamic boosting mechanism for relativistic jets explored by Aloy & Rezzolla (2006) using the RAISHIN code. Simulation results show that the presence of a magnetic field changes the properties of the shock interface between the tenuous, overpressured jet (V^z j) flowing tangentially to a dense external medium. We find that magnetic fields can lead to more efficient acceleration of the jet, in comparison to the pure-hydrodynamic case. A "poloidal" magnetic field (B^z), tangent to the interface and parallel to the jet flow, produces both a stronger outward moving shock and a stronger inward moving rarefaction wave. This leads to a large velocity component normal to the interface in addition to acceleration tangent to the interface, and the jet is thus accelerated to larger Lorentz factors than those obtained in the pure-hydrodynamic case. Likewise, a strong "toroidal" magnetic field (B^y), tangent to the interface but perpendicular to the jet flow, also leads to stronger acceleration tangent to the shock interface relative to the pure-hydrodynamic case. Thus. the presence and relative orientation of a magnetic field in relativistic jets can significant modify the hydrodynamic boost mechanism studied by Aloy & Rezzolla (2006).

Mizuno, Yosuke↗

Simulations of Solar Wind Turbulence

Recently we have restructured our approach to simulating magnetohydrodynamic (MHD) turbulence in the solar wind. Previously, we had defined a 'virtual' heliosphere that contained, for example, a tilted rotating current sheet, microstreams, quasi-two-dimensional fluctuations as well as Alfven waves. In this new version of the code, we use the global, time-stationary, WKB Alfven wave-driven solar wind model developed by Usmanov and described in Usmanov and Goldstein [2003] to define the initial state of the system. Consequently, current sheets, and fast and slow streams are computed self-consistently from an inner, photospheric, boundary. To this steady-state configuration, we add fluctuations close to, but above, the surface where the flow become super-Alfvenic. The time-dependent MHD equations are then solved using a semi-discrete third-order Central Weighted Essentially Non-Oscillatory (CWENO) numerical scheme. The computational domain now includes the entire sphere; the geometrical singularity at the poles is removed using the multiple grid approach described in Usmanov [1996]. Wave packets are introduced at the inner boundary such as to satisfy Faraday's Law [Yeh and Dryer, 1985] and their nonlinear evolution are followed in time.

Goldstein, Melvyn L.↗

Searching for Turbulent Spectra in Mutliple Spacecraft Observations

The subject of magnetohydrodynamic simulations in the 1980's, magnetic turbulence has reemerged as a favored explanation for the production of fast magnetic reconnection. Recent three-dimensional simulations reveal the formation of secondary islands that serve to shorten the current sheet and increase the accelerating electric field, while both simulations and observations witness electron holes whose collapse energizes electrons. However, little work has been done to investigate the presence of turbulent magnetic and electric field spectra in association with in situ reconnection observations in space. In this study we address this deficiency by analyzing the wavelength spectra of two reconnection events established by multiple spacecraft observations in previous publications. We apply a k-filtering analysis to determine the wave vectors. To optimize the quality of the k-filtering analysis and minimize aliasing, we select multi-spacecraft observations having optimal spacecraft tetrahedral geometry and dimensions that resolve wavelengths within the ion diffusion region. We have chosen one dayside and one nightside reconnection event each observed by the Cluster spacecraft. Spatial correlation between the reconnecting current sheet and a turbulent magnetic spectrum supports an association between turbulence and reconnection in space. The k-filtering method, a wavelet time-frequency spectrum, and a polarization analysis serve to identify wave modes associated with the turbulence. A calculation of the reconnection electric field based on the normal component of the reconnecting magnetic field provides an estimate of the reconnection rate.

Wendel, Deirdre E.↗

Scientific Visualization to Study Flux Transfer Events at the Community Coordinated Modeling Center

In this paper we present results of modeling of reconnection at the dayside magnetopause with subsequent development of flux transfer event signatures. The tools used include new methods that have been added to the suite of visualization methods that are used at the Community Coordinated Modeling Center (CCMC). Flux transfer events result from localized reconnection that connect magnetosheath magnetic field and plasma with magnetospheric fields and plasma and results in flux rope structures that span the dayside magnetopause. The onset of flux rope formation and the three-dimensional structure of flux ropes are studied as they have been modeled by high-resolution magnetohydrodynamic simulations of the dayside magnetosphere of the Earth. We show that flux transfer events are complex three-dimensional structures that require modern visualization and analysis techniques. Two suites of visualization methods are presented and we demonstrate the usefulness of those methods through the CCMC web site to the general science user.

Rastatter, Lutz↗