SEARCH · Search NASA
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.
Spatio-temporal evolution, of large-amplitude Alfven waves, is investigated by means of their characteristic evolution equations as well as by MHD simulations.
We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of collapsars including a rotating black hole. Initially, we assume that the core collapse has failed in this star. A rotating black hole of a few solar masses is inserted by hand into the calculation. The simulation results show the formation of a disklike structure and the generation of a jetlike outflow near the central black hole. The jetlike outflow propagates and accelerated mainly by the magnetic field. The total jet velocity is approximately 0.3c. When the rotation of the black hole is faster, the magnetic field is twisted strongly owing to the frame-dragging effect. The magnetic energy stored by the twisting magnetic field is directly converted to kinetic energy of the jet rather than propagating as an Alfven wave. Thus, as the rotation of the black hole becomes faster, the poloidal velocity of the jet becomes faster.
New tools for data access and visualization promise to make the analysis of space plasma data both more efficient and more powerful, especially for answering questions about the global structure and dynamics of the Sun-Earth system. We will show how new existing tools (particularly the Virtual Space Physics Observatory-VSPO-and the Visual System for Browsing, Analysis and Retrieval of Data-ViSBARD; look for the acronyms in Google) already provide rapid access to such information as spacecraft orbits, browse plots, and detailed data, as well as visualizations that can quickly unite our view of multispacecraft observations. We will show movies illustrating multispacecraft observations of the solar wind and magnetosphere during a magnetic storm, and of simulations of 3 0-spacecraft observations derived from MHD simulations of the magnetosphere sampled along likely trajectories of the spacecraft for the MagCon mission. An important issue remaining to be solved is how best to integrate simulation data and services into the Virtual Observatory environment, and this talk will hopefully stimulate further discussion along these lines.
We examine how the initial state (pre-event corona) affects the numerical MHD simulation for a coronal mass ejection (CME). Earlier simulations based on a pre-event corona with a homogeneous density and temperature distribution at lower boundary (i.e. solar surface) have been used to analyze the role of streamer properties in determining the characteristics of loop-like transients. The present paper extends these studies to show how a broader class of global coronal properties leads not only to different types of CME's, but also modifies the adjacent quiet corona and/or coronal holes. We consider four pre-event coronal cases: (1) Constant boundary conditions and a polytropic gas with gamma = 1.05; (2) Non-constant (latitude dependent) boundary conditions and a polytropic gas with gamma = 1.05; (3) Constant boundary conditions with a volumetric energy source and gamma = 1.67; (4) Non-constant (latitude dependent) boundary conditions with a volumetric energy source and gamma = 1.67. In all models, the pre-event magnetic fields separate the corona into closed field regions (streamers) and open field regions. The CME's initiation is simulated by introducing at the base of the corona, within the streamer region, a standard pressure pulse and velocity change. Boundary values are determined using MHD characteristic theory. The simulations show how different CME's, including loop-like transients, clouds, and bright rays, might occur. There are significant new features in comparison to published results. We conclude that the pre-event corona is a crucial factor in dictating CME's properties.
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. Recent advancements in solar observations and numerical simulations provide a basis for developing 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. The goal is to develop new strategic capabilities that enable characterization and prediction of 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. The integration of first-principle modeling of solar magnetism and flow dynamics with real-time observational data via advanced data assimilation methods is a new, transformative step in space weather research and prediction. This approach will substantially enhance an existing model of magnetic flux distribution and transport developed by the Air Force Research Lab. The development plan is to use the Space Weather Modeling Framework (SWMF) to develop Coupled Models for Emerging flux Simulations (CMES) that couples three existing models: (1) an MHD formulation with the anelastic approximation to simulate the deep convection zone (FSAM code), (2) an MHD formulation with full compressible Navier-Stokes equations and a detailed description of radiative transfer and thermodynamics to simulate near-surface convection and the photosphere (Stagger code), and (3) an MHD formulation with full, compressible Navier-Stokes equations and an approximate description of radiative transfer and heating to simulate the corona (Module in BATS-R-US). CMES will enable simulations of the emergence of magnetic structures from the deep convection zone to the corona. Finally, a plan will be summarized on the development of a Flux Emergence Prediction Tool (FEPT) in which helioseismology-derived data and vector magnetic maps are assimilated into CMES that couples the dynamics of magnetic flux from the deep interior to the corona.
We examine how the initial state (pre-event corona) affects the numerical MHD simulation for a coronal mass ejection (CME). Earlier simulations based on a pre-event corona with a homogeneous density and temperature distribution, at the lower boundary (i.e., solar surface) have been used to analyze the role of streamer properties in determining the characteristics of loop-like transients. The present paper extends these studies to show how a broader class of global coronal properties leads not only to different types of CME's, but also modifies the adjacent quiet corona and/or coronal holes. We consider four pre-event coronal cases: (1) constant boundary conditions and a polytropic gas with gamma = 1.05; (2) non-constant (latitude dependent) boundary conditions and a polytropic gas with gamma = 1.05; (3) constant boundary conditions with a volumetric energy source and gamma = 1.67; (4) non-constant (latitude dependent) boundary conditions with a volumetric energy source and gamma = 1.67. In all models, the pre-event magnetic fields separate the corona into closed field regions (streamers) and open field regions. The CME's initiation is simulated by introducing at the base of the corona, within the streamer region, a standard pressure pulse and velocity change. Boundary values are determined using magnetohydrodynamic (MHD) characteristic theory. The simulations show how different CME's, including loop-like transients, clouds and bright rays, might occur. There are significant new features in comparison to published results. We conclude that the pre-event corona is a crucial factor in dictating CME's properties.
On February 9, 1995, the Comprehensive Plasma Instrumentation (CPI) on the Geotail spacecraft observed a complex, structured ion distribution function near the magnetotail midplane at x approximately -30 R(sub E). On this same day the Wind spacecraft observed a quiet solar wind and an interplanetary magnetic field (IMF) that was northward for more than five hours, and an IMF B(sub y) component with a magnitude comparable to that of the RAF B(sub z) component. In this study, we determined the sources of the ions in this distribution function by following approximately 90,000 ion trajectories backward in time, using the time-dependent electric and magnetic fields obtained from a global MHD simulation. The Wind observations were used as input for the MHD model. The ion distribution function observed by Geotail at 1347 UT was found to consist primarily of particles from the dawn side low latitude boundary layer (LLBL) and from the dusk side LLBL; fewer than 2% of the particles originated in the ionosphere.
We study the ion acceleration in a mesoscale, spontaneously fragmenting flare current sheet (SFCS) characterized by the presence of a plasmoid cascade. The main subject of our investigation is to determine whether and how plasmoid cascades at intermediate scales in a fragmented current sheet of a solar flare can impact the (preferential) acceleration of specific ions. The time evolution of the SFCS is obtained from high-resolution 2.5D MHD simulations. The ion trajectories (in the background fields resulting from the MHD model), energies, and pitch angles are calculated using a relativistic test-particle code based on the half-acceleration–rotation–half-acceleration method. For light ions, the main acceleration effects of electromagnetic fields within the SFCS are analyzed using the guiding center approximation. We identify regions with the most-efficient ion acceleration within the SFCS, the accelerator efficiency, and spectra of the accelerated ions. The influence of the charge-to-mass ratio on ion behavior is also studied and resulting ion abundances are compared with observational data. The main ion acceleration takes place in the regions with a strong polarization term, which is part of the first-order Fermi acceleration. Because the term is mass dependent, heavier ions undergo preferential acceleration. The ion energy spectra, abundance-enhancement factors, and differential fluxes, obtained from the model, exhibit power-law profiles, in agreement with observed solar energetic particle events. Nonetheless, the obtained slopes for the abundance-enhancement factor do not exactly match the observed data. The computed slopes and profiles are not sensitive to changes in the initial plasma temperature.
We describe the results of over 40 MHD simulations, performed on Caltech/JPL parallel supercomputers, of the coronae of magnetized accretion disks around compact objects.
We report the results of MHD simulations of Earth's magnetosphere for idealized steady solar wind plasma and interplanetary magnetic field (IMF) conditions. The simulations feature purely northward and southward magnetic fields and were designed to study turbulence in the magnetotail plasma sheet. We found that the power spectral densities (PSDs) for both northward and southward IMF had the characteristics of turbulent flow. In both cases, the PSDs showed the three scale ranges expected from theory: the energy-containing scale, the inertial range, and the dissipative range. The results were generally consistent with in-situ observations and theoretical predictions. While the two cases studied, northward and southward IMF, had some similar characteristics, there were significant differences as well. For southward IMF, localized reconnection was the main energy source for the turbulence. For northward IMF, remnant reconnection contributed to driving the turbulence. Boundary waves may also have contributed. In both cases, the PSD slopes had spatial distributions in the dissipative range that reflected the pattern of resistive dissipation. For southward IMF there was a trend toward steeper slopes in the dissipative range with distance down the tail. For northward IMF there was a marked dusk-dawn asymmetry with steeper slopes on the dusk side of the tail. The inertial scale PSDs had a dusk-dawn symmetry during the northward IMF interval with steeper slopes on the dawn side. This asymmetry was not found in the distribution of inertial range slopes for southward IMF. The inertial range PSD slopes were clustered around values close to the theoretical expectation for both northward and southward IMF. In the dissipative range, however, the slopes were broadly distributed and the median values were significantly different, consistent with a different distribution of resistivity.
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.
In situ observations indicate that the dynamical processes in the geoplasma environment generally entail localized intermittent processes and anomalous global transports. It was suggested by T. Chang that instead of considering the turbulence as a mixture of interacting waves, such type of patchy intermittency could be more easily understood in terms of the development, interaction, merging, preferential acceleration and evolution of coherent magnetic structures. In this three-year project, we have used direct numerical MHD simulations to study some aspects of the MHD dynamics in Chang's model. Our large-scale numerical calculations and simulations have been supplemented by and coordinated with theoretical studies conducted by Chang and his colleagues.
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.
Context. The generation of the slow solar wind remains an open problem in heliophysics. One of the current theories among those aimed at explaining the injection of coronal plasma in the interplanetary medium is based on interchange reconnection. It assumes that the exchange of magnetic connectivity between closed and open fields allows the injection of coronal plasma in the interplanetary medium to travel along the newly reconnected open field. However, the exact mechanism underlying this effect is still poorly understood. Aims. Our objective is to study this scenario in a particular magnetic structure of the solar corona: a pseudo-streamer. This topological structure lies at the interface between open and closed magnetic field and is thought to be involved in the generation of the slow solar wind. Methods. We performed innovative 3D magnetohydrodynamic (MHD) simulations of the solar corona with a pseudo-streamer, using the Adaptively Refined MHD Solver (ARMS). By perturbing the quasi-steady ambient state with a simple photospheric, large-scale velocity flow, we were able to generate a complex dynamics of the open-and-closed boundary of the pseudo-streamer. We studied the evolution of the connectivity of numerous field lines to understand its precise dynamics. Results. We witnessed different scenarios of opening of the magnetic field initially closed under the pseudo-streamer: one-step interchange reconnection dynamics, along with more complex scenarios, including a coupling between pseudo-streamer and helmet streamer, as well as back-and-forth reconnections between open and closed connectivity domains. Finally, our analysis revealed large-scale motions of a newly opened magnetic field high in the corona that may be explained by slipping reconnection. Conclusions. By introducing a new analysis method for the magnetic connectivity evolution based on distinct closed-field domains, this study provides an understanding of the precise dynamics underway during the opening of a closed field, which enables the injection of closed-field, coronal plasma in the interplanetary medium. Further studies shall provide synthetic observations for these diverse outgoing flows, which could be measured by Parker Solar Probe and Solar Orbiter.
Realistic 3D radiative MHD simulations open new perspectives for understanding the turbulent dynamics of the solar surface, its coupling to the atmosphere, and the physical mechanisms of generation and transport of non-thermal energy. Traditionally, plasma eruptions and wave phenomena in the solar atmosphere are modeled by prescribing artificial driving mechanisms using magnetic or gas pressure forces that might arise from magnetic field emergence or reconnection instabilities. In contrast, our 'ab initio' simulations provide a realistic description of solar dynamics naturally driven by solar energy flow. By simulating the upper convection zone and the solar atmosphere, we can investigate in detail the physical processes of turbulent magnetoconvection, generation and amplification of magnetic fields, excitation of MHD waves, and plasma eruptions. We present recent simulation results of the multi-scale dynamics of quiet-Sun regions, and energetic effects in the atmosphere and compare with observations. For the comparisons we calculate synthetic spectro-polarimetric data to model observational data of SDO, Hinode, and New Solar Telescope.
The mini-magnetospheric plasma propulsion (M2P2) device, originally proposed by Winglee et al., predicts that a 15-km standoff distance (or 20-km cross-sectional dimension) of the magnetic bubble will provide for sufficient momentum transfer from the solar wind to accelerate a spacecraft to unprecedented speeds of 50 C80 km/s after an acceleration period of 3 mo. Such velocities will enable travel out of the solar system in period of 7 yr almost an order of magnitude improvement over present chemical-based propulsion systems. However, for the parameters of the simulation of Winglee et al., a fluid model for the interaction of M2P2 with the solar wind is not valid. It is assumed in the magnetohydrodynamic (MHD) fluid model, normally applied to planetary magnetospheres, that the characteristic scale size is much greater than the Larmor radius and ion skin depth of the solar wind. In the case of M2P2, the size of the magnetic bubble is actually less than or comparable to the scale of these characteristic parameters. Therefore, a kinetic approach, which addresses the small-scale physical mechanisms, must be used. A two-component approach to determining a preliminary estimate of the momentum transfer to the plasma sail has been adopted. The first component is a self-consistent MHD simulation of the small-scale expansion phase of the magnetic bubble. The fluid treatment is valid to roughly 5 km from the source and the steady-state MHD solution at the 5 km boundary was then used as initial conditions for the hybrid simulation. The hybrid simulations showed that the forces delivered to the innermost regions of the plasma sail are considerably ( 10 times) smaller than the MHD counterpart, are dominated by the magnetic field pressure gradient, and are directed primarily in the transverse direction.
A forced magnetic reconnection process with a temporal evolution of resistivity is studied for a plasma sheet with a nonuniform resistivity profile based on the nonlocal mode structure of the lower hybrid drift type instability. The growth rate of the mode found is almost independent of the resistivity at the neutral sheet, but depends on the resistivity of the region of maximum density gradient away from the neutral sheet. This is studied by using both a nonlinear numerical MHD simulation and a linear theory. The mode may be relevant to the prevalent theoretical concept of MHD reconnection and the localized anomalous resistivity profile based on the lower hybrid drift instability.
Based on satellite observations of a recurrent solar wind structure conducted in 1974, an MHD simulation model, and input functions generated from plasma and magnetic field data, the continuing evolution of the solar wind structure outside 5 AU is studied. The model uses the Rankine-Hugoniot relations to describe the jumps in flow properties across the shocks, and it treats shocks as surfaces of discontinuity with zero thickness. Two interaction processes (the collision and the merging of shocks) play important roles in restructuring the solar wind in the outer heliosphere. The simulation result shows that period doubling occurs between 5 and 10 AU. The recurrent solar wind appears to be a persistent new structure between 10 and 20 AU, and it consists of one merged interaction region per solar rotation.