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 181 records · Page 10

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

Interaction of Strong Transient Interplanetary Disturbances with the Dayside Magnetosphere

The objective of the investigation was to gain an understanding of the complex response of the magnetosphere to strong transient interplanetary disturbances. Because the project was only funded for a year, the investigation focused on only one of the three topics proposed in the original three year proposal. We investigated the response of the dayside auroral region to strong transient interplanetary disturbances. The method of the investigation was to use three-dimensional magnetohydrodynamic (MHD) simulations that employ measurements of the solar wind conditions upstream the bow shock to model actual events and then to compare the simulation results with observations from spacecraft located at downstream locations. We modeled an event that occurred on July 14,2000, for which both CLUSTER and IMAGE simultaneous observations were available. The event was marked by high solar wind dynamic pressure and a strong IMF By component. Comparisons showed a very good agreement between intensifications in the auroral emissions measured by IMAGE

Berchem, Jean↗

Polar, Cluster and SuperDARN Evidence for High-Latitude Merging during Southward IMF: Temporal/Spatial Evolution

Magnetic merging on the dayside magnetopause often occurs at high latitudes. Polar measured fluxes of accelerated ions and wave Poynting vectors while skimming the subsolar magnetopause. The measurements indicate that their source was located to the north of the spacecraft, well removed from expected component merging sites. This represents the first use of wave Poynting flux as a merging discriminator at the magnetopause. We argue that wave Poynting vectors, like accelerated particle fluxes and the Walen tests, are necessary, but not sufficient, conditions, for identifying merging events. The Polar data are complemented with nearly simultaneous measurements from Cluster in the northern cusp, with correlated observations from the SuperDARN radar, to show that the locations and rates of merging vary. Magnetohydrodynamic (MHD) simulations are used to place the measurements into a global context. The MHD simulations confirm the existence of a high-latitude merging site and suggest that Polar and SuperDARN observed effects are attributable to both exhaust regions of a temporally varying X-line. A survey of 13 merging events places the location at high latitudes whenever the interplanetary magnetic field (IMF) clock angle is less than approximately 150 degrees. While inferred high-latitude merging sites favor the antiparallel merging hypothesis, our data alone cannot exclude the possible existence of a guide field. Merging can even move away from equatorial latitudes when the IMF has a strong southward component. MHD simulations suggest that this happens when the dipole tilt angle increases or when IMF B(sub X) increases the effective dipole tilt.

Maynard, N. C.↗

Magnetohydrodynamic Modeling of the Jovian Magnetosphere

Under this grant we have undertaken a series of magnetohydrodynamic (MHD) simulation and data analysis studies to help better understand the configuration and dynamics of Jupiter's magnetosphere. We approached our studies of Jupiter's magnetosphere in two ways. First we carried out a number of studies using our existing MHD code. We carried out simulation studies of Jupiter s magnetospheric boundaries and their dependence on solar wind parameters, we studied the current systems which give the Jovian magnetosphere its unique configuration and we modeled the dynamics of Jupiter s magnetosphere following a northward turning of the interplanetary magnetic field (IMF). Second we worked to develop a new simulation code for studies of outer planet magnetospheres.

Walker, Raymond↗

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗

Comparison Between Vortices Created and Evolving During Fixed and Dynamic Solar Wind Conditions

We employ Magnetohydrodynamic (MHD) simulations to examine the creation and evolution of plasma vortices within the Earth's magnetosphere for steady solar wind plasma conditions. Very few vortices form during intervals of such solar wind conditions. Those that do remain in fixed positions for long periods (often hours) and exhibit rotation axes that point primarily in the x or y direction, parallel (or antiparallel) to the local magnetospheric magnetic field direction. Occasionally, the orientation of the axes rotates from the x direction to another direction. We compare our results with simulations previously done for unsteady solar wind conditions. By contrast, these vortices that form during intervals of varying solar wind conditions exhibit durations ranging from seconds (in the case of those with axes in the x or y direction) to minutes (in the case of those with axes in the z direction) and convect antisunward. The local-time dependent sense of rotation seen in these previously reported vortices suggests an interpretation in terms of the Kelvin-Helmholtz instability. For steady conditions, the biggest vortices developed on the dayside (about 6R(E) in diameter), had their rotation axes aligned with the y direction and had the longest periods of duration. We attribute these vortices to the flows set up by reconnection on the high latitude magnetopause during intervals of northward Interplanetary Magnetic Field (IMF) orientation. This is the first time that vortices due to high-latitude reconnection have been visualized. The model also successfully predicts the principal characteristics of previously reported plasma vortices within the magnetosphere, namely their dimension, flow velocities, and durations.

solar wind↗

Integrated Real-Time Modeling System for Heliospheric Space Weather Forecasting

The Sun lies at the center of space weather and is the source of its variability. The primary input to coronal and solar wind models is the activity of the magnetic field in the solar photosphere. What is needed is to develop physics-based models for the dynamics of the magnetic field from the deep convection zone of the Sun to the corona with the goal of providing robust near real-time boundary conditions at the base of space weather forecast models. NASA strategic capabilities efforts toward this goal will be summarized during the talk. These have the objective to enable characterizing and predicting the magnetic field structure and flow dynamics of the Sun by assimilating data from helioseismology and magnetic field observations into physics-based realistic magnetohydrodynamics (MHD) simulations. It is argued that an effort in integrating first-principle modeling of solar magnetism and flow dynamics with real-time observational data via advanced data assimilation methods is needed for transformative steps in space weather research and prediction.

Mansour, Nagi N.↗

Modeling a Coronal Mass Ejection from an Extended Filament Channel. I. Eruption and Early Evolution

We present observations and modeling of the magnetic field configuration, morphology, and dynamics of a large-scale, high-latitude filament eruption observed by the Solar Dynamics Observatory. We analyze the 2015 July 9–10 filament eruption and the evolution of the resulting coronal mass ejection (CME) through the solar corona. The slow streamer-blowout CME leaves behind an elongated post-eruption arcade above the extended polarity inversion line that is only poorly visible in extreme ultraviolet (EUV) disk observations and does not resemble a typical bright flare-loop system. Magnetohydrodynamic (MHD) simulation results from our data-inspired modeling of this eruption compare favorably with the EUV and white-light coronagraph observations. We estimate the reconnection flux from the simulation’s flare-arcade growth and examine the magnetic-field orientation and evolution of the erupting prominence, highlighting the transition from an erupting sheared-arcade filament channel into a streamer-blowout flux-rope CME. Our results represent the first numerical modeling of a global-scale filament eruption where multiple ambiguous and complex observational signatures in EUV and white light can be fully understood and explained with the MHD simulation. In this context, our findings also suggest that the so-called “stealth CME” classification, as a driver of unexpected or “problem” geomagnetic storms, belongs more to a continuum of observable/non-observable signatures than to separate or distinct eruption processes.

Benjamin J Lynch↗

Mini-disk accretion onto spinning black hole binaries: quasi-periodicities and outflows

We perform a full 3D general relativistic magnetohydrodynamical (GRMHD) simulation of an equal-mass, spinning, binary black hole approaching merger, surrounded by a circumbinary disk and with mini-disks around each black hole. For this purpose, we evolve the ideal GRMHD equations on top of an approximated spacetime for the binary that is valid in every position of space, including the black hole horizons, during the inspiral regime. We use relaxed initial data for the circumbinary disk from a previous long-term simulation, where the accretion is dominated by an m = 1 overdensity called the lump. We compare our new spinning simulation with a previous non-spinning run, studying how spin influences the mini-disk properties. We analyze the accretion from the inner edge of the lump to the black hole, focusing on the angular momentum budget of the fluid around the mini-disks. We find that mini-disks in the spinning case have more mass over a cycle than the non-spinning case. However, in both cases we find most of the mass received by the black holes is delivered by direct plunging of material from the lump. We also analyze the morphology and variability of the electromagnetic fluxes and we find they share the same periodicities of the accretion rate. In the spinning case, we find that the outflows are stronger than the non-spinning case. Our results will be useful to understand and produce realistic synthetic light curves and spectra, which can be used in future observations.

Astrophysics↗