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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

General Relativistic Radiative Transfer and General Relativistic MHD Simulations of Accretion and Outflows of Black Holes

We have calculated the emission from relativistic flows in black hole systems using a fully general relativistic radiative transfer, with flow structures obtained by general relativistic magnetohydrodynamic simulations. We consider thermal free-free emission and thermal synchrotron emission. Bright filament-like features are found protruding (visually) from the accretion disk surface, which are enhancements of synchrotron emission when the magnetic field is roughly aligned with the line-of-sight in the co-moving frame. The features move back and forth as the accretion flow evolves, but their visibility and morphology are robust. We propose that variations and location drifts of the features are responsible for certain X-ray quasi-periodic oscillations (QPOs) observed in black-hole X-ray binaries.

Fuerst, Steven V.↗
NASA NTRS1995

Physics of magnetospheric boundary layers

This final report was concerned with the ideas that: (1) magnetospheric boundary layers link disparate regions of the magnetosphere-solar wind system together; and (2) global behavior of the magnetosphere can be understood only by understanding its internal linking mechanisms and those with the solar wind. The research project involved simultaneous research on the global-, meso-, and micro-scale physics of the magnetosphere and its boundary layers, which included the bow shock, the magnetosheath, the plasma sheet boundary layer, and the ionosphere. Analytic, numerical, and simulation projects were performed on these subjects, as well as comparisons of theoretical results with observational data. Other related activity included in the research included: (1) prediction of geomagnetic activity; (2) global MHD (magnetohydrodynamic) simulations; (3) Alfven resonance heating; and (4) Critical Ionization Velocity (CIV) effect. In the appendixes are list of personnel involved, list of papers published; and reprints or photocopies of papers produced for this report.

Cairns, Iver H.↗
NASA NTRS2021

Spatially Separated Electron and Proton Beams in a Simulated Solar Coronal Jet

Magnetic reconnection is widely accepted to be a major contributor to non-thermal particle acceleration in the solar atmosphere. In this paper we investigate particle acceleration during the impulsive phase of a coronal jet, which involves bursty reconnection at a magnetic null point. A test-particle approach is employed, using electromagnetic fields from a magnetohydrodynamic simulation of such a jet. Protons and electrons are found to be accelerated non-thermally both downwards toward the domain’s lower boundary and the solar photosphere, and outwards along the axis of the coronal jet and into the heliosphere. A key finding is that a circular ribbon of particle deposition on the photosphere is predicted, with the protons and electrons concentrated in different parts of the ribbon. Furthermore, the outgoing protons and electrons form two spatially separated beams parallel to the axis of the jet, signatures that may be observable inin-situobservations of the heliosphere.

Ross Pallister↗
NASA NTRS2011

CUDA Simulation of Homogeneous, Incompressible Turbulence

We discuss very fast Compute Unified Device Architecture (CUDA) simulations of ideal homogeneous incompressible turbulence based on Fourier models. These models have associated statistical theories that predict that Fourier coefficients of fluid velocity and magnetic fields (if present) are zero-mean random variables. Prior numerical simulations have shown that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation. We review the theoretical basis of this "broken ergodicity" as applied to 2-D and 3-D fluid and magnetohydrodynamic simulations of homogeneous turbulence. Our new simulations examine the phenomenon of broken ergodicity through very long time and large grid size runs performed on a state-of-the-art CUDA platform. Results comparing various CUDA hardware configurations and grid sizes are discussed. NS and MHD results are compared.

Morin, Lee↗
NASA NTRS1999

Kelvin Helmholtz Instability at the Equatorial Magnetotail Boundary: MHD Simulation and Comparison with Geotail Observations

On March 24, 1995 the Geotail spacecraft observed large fluctuations of the magnetic field and plasma properties in the Low Latitude Boundary Layer (LLBL) about 15 R(sub E) tailward of the dusk meridian. Although the magnetospheric and the magnetosheath field were strongly northward, the B(sub z) component showed strong short duration fluctuations in which B(sub z) could even reach negative values. We have used two-dimensional magnetohydrodynamic simulations with magnetospheric and magnetosheath input parameters specifically chosen for this. Geotail event to identify the processes which cause the observed boundary properties. It is shown that these fluctuations can be explained by the Kelvin-Helmholtz instability if the k vector of the instability has a component along the magnetic field direction. The simulation results show many of the characteristic properties of the Geotail observations. In particular, the quasi-periodic strong fluctuations are well explained by satellite crossings through the Kelvin-Helmholtz vortices. It is illustrated how the interior structure of the Kelvin-Helmholtz vortices leads to the rapid fluctuations in the Geotail observations. Our results suggest an average Kelvin-Helmholtz wavelength of about 5 R(sub E) with a vortex size of close to 2 R(sub E) for an average repetition time of 2.5 minutes. The growth time for these waves implies a source region of about 10 to 16 R(sub E) upstream from the location of the Geotail spacecraft (i.e., near the dusk meridian). The results also indicate a considerable mass transport of magnetosheath material into the magnetosphere by magnetic reconnection in the Kelvin-Helmholtz vortices.

Fairfield, Donald H.↗
NASA NTRS1993

Current loop coalescence studied by 3-D electromagnetic particle code

Solar flare plasma data from the Yohkoh satellite is analyzed. The interactions of current loops were observed in the active regions on the Sun. This observation pointed out the importance of the idea that the solar flare is generated by the coalescence of current loops. The three dimensional electromagnetic particle simulations are to help in understanding the global interaction between two current loops including the evolution of the twist of loops due to instabilities. Associated rapid dynamics of current loop coalescence such as reconnection, shock waves and associated kinetic processes such as energy transfer, acceleration of particles, and electromagnetic emissions are to be studied by the code to complement analytical theories and magnetohydrodynamic simulations of the current loop coalescence. The simulation results show the strong interactions between two current loops, beam and whistler instabilities, and associated parallel and perpendicular particle heating.

Nishikawa, Ken-Ichi↗
NASA NTRS2010

Broken Ergodicity in Ideal, Homogeneous, Incompressible Turbulence

We discuss the statistical mechanics of numerical models of ideal homogeneous, incompressible turbulence and their relevance for dissipative fluids and magnetofluids. These numerical models are based on Fourier series and the relevant statistical theory predicts that Fourier coefficients of fluid velocity and magnetic fields (if present) are zero-mean random variables. However, numerical simulations clearly show that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation. We explain this phenomena in terms of broken ergodicity', which is defined to occur when dynamical behavior does not match ensemble predictions on very long time-scales. We review the theoretical basis of broken ergodicity, apply it to 2-D and 3-D fluid and magnetohydrodynamic simulations of homogeneous turbulence, and show new results from simulations using GPU (graphical processing unit) computers.

Morin, Lee↗
NASA NTRS2010

Observations and Simulations of Electron Dynamics Near an Active Neutral Line

Recent observations in the Earth's magnetotail have shown rapid increases in the fluxes of energetic electrons with energies up to 100's of keV associated with dipolarization fronts that propagate into the inner magnetosphere. On August 15, 2001 the four Cluster spacecraft located slightly dawnward of midnight (yGSM approx. -5.4RE) at xGSM approx. -18RE observed a series of earthward propagating dipolarization fronts [Hwang et al., 2010]. At least 6 dipolarization fronts were observed in a 20m interval. Unlike previously reported cases the fluxes of electrons up to 95keV decreased during the passage of the first three fronts over the spacecraft. The energetic electron fluxes increased during the passage of the last three fronts. We have performed a global magnetohydrodynamic simulation of this event using solar wind observations from the ACE satellite to drive the simulation. In the simulation a very complex reconnection system in the near-Earth tail at XGSM approx. -20RE launched a series of earthward propagating dipolarization fronts that are similar to those observed on Cluster. The simulation results indicate that the Cluster spacecraft were just earthward of the reconnection site. In this paper we will present a study of the dynamics of electrons associated with these events by using the large-scale kinetic simulation approach in which we launch a large number of electrons into the electric and magnetic fields from this simulation.

Goldstein, M. L.↗
NASA NTRS2015

Deep Chandra Observation and Numerical Studies of the Nearest Cluster Cold Front in the Sky

We present the results of a very deep (500 ks) Chandra observation, along with tailored numerical simulations, of the nearest, best resolved cluster cold front in the sky, which lies 90 kpc (19 arcmin) to the north-west of M87. The northern part of the front appears the sharpest, with a width smaller than 2.5 kpc (1.5 Coulomb mean free paths; at 99 per cent confidence). Everywhere along the front, the temperature discontinuity is narrower than 4-8 kpc and the metallicity gradient is narrower than 6 kpc, indicating that diffusion, conduction and mixing are suppressed across the interface. Such transport processes can be naturally suppressed by magnetic fields aligned with the cold front. Interestingly, comparison to magnetohydrodynamic simulations indicates that in order to maintain the observed sharp density and temperature discontinuities, conduction must also be suppressed along the magnetic field lines. However, the northwestern part of the cold front is observed to have a non-zero width. While other explanations are possible, the broadening is consistent with the presence of Kelvin-Helmholtz instabilities (KHI) on length-scales of a few kpc. Based on comparison with simulations, the presence of KHI would imply that the effective viscosity of the intracluster medium is suppressed by more than an order of magnitude with respect to the isotropic Spitzer-like temperature dependent viscosity. Underneath the cold front, we observe quasi-linear features that are approximately 10 per cent brighter than the surrounding gas and are separated by approximately 15 kpc from each other in projection. Comparison to tailored numerical simulations suggests that the observed phenomena may be due to the amplification of magnetic fields by gas sloshing in wide layers below the cold front, where the magnetic pressure reaches approximately 5-10 per cent of the thermal pressure, reducing the gas density between the bright features.

hydrodynamics↗
NASA NTRS1987

Evolution of the solar wind structure in the outer heliosphere

Shocks and interaction regions play very important roles in the evolution of large-scale solar wind structure in the outer heliosphere. This study is based on (1) plasma and magnetic field data observed from Voyager and Pioneer spacecraft, and (2) a quantitative magnetohydrodynamic simulation model. Interaction regions bounded by a forward and a reverse shock begin to form near 1 AU at the leading edges of a large-scale stream. The total pressure in the region is greater than the ambient pressure by a factor of ten or more. Large jumps in pressure remain as a prominant feature of the interplanetary structure even as the jumps in flow speed become less visible in the outer heliosphere. The propagation of the forward and reverse shocks widens the dimension of an interaction region. As a result, two interaction regions belonging to neighboring streams coalesce to form a merged interaction region (MIR). Collision and merging of shocks take place during the coalescence process. Two MIRs can themselves merge again at greater heliocentric distances. Simulation results agree well with spacecraft observations, and they explain major restructuring of the solar wind in the outer heliosphere.

Whang, Y. C.↗
NASA NTRS1990

The nonlinear evolution of magnetized solar filaments

Thermal instability driven by optically thin radiation is believed to initiate the formation of plasma filaments in the solar corona. The fact that filaments are observed generally to separate regions of opposite, line-of-sight, magnetic polarity in the underlying photosphere suggests that filament formation requires the presence of a highly sheared, local magnetic field. Two-dimensional, nonlinear, magnetohydrodynamic simulations of the local genesis and growth of solar filaments in a force-free, sheared, magnetic field were performed, and the evolution of generic perturbations possessing broad spatial profiles was traced. It was found that simulations of the evolution of initial random-noise perturbations produce filamentary plasma structures that exhibit densities and temperatures characteristic of observed solar filaments. Furthermore, in each of these simulations, the filament axis lies at a finite angle with respect to the local magnetic field, consistent with solar observations.

Sparks, L.↗
NASA NTRS1991

Magnetic field draping at the low-latitude magnetopause

Magnetohydrodynamic simulations are used to investigate the structure of the low-latitude magnetopause for interplanetary magnetic field conditions with a dominant southward component. The structure is self-consistently calculated as an initial-value problem in which the system is allowed to evolve into a quasi-steady state. All components of the 3D velocity and magnetic field as well as compressibility, resistivity, and viscosity are included in the 2D calculation. The simulation model shows that magnetic field draping can occur at the magnetopause boundary when magnetic merging takes place in the presence of a tangential shear flow. For 'normal' (positive Bx) draping, the higher-latitude portion of the field lines are curved toward the sun on the magnetospheric side of the magnetopause and away from the sun on the magnetosheath side. The thickness of the normal draping structure scales with the viscosity. The field-aligned current system that accompanies normal magnetic draping is consistent with the sense of the region 1 currents that flow into the dayside ionosphere.

Richard, R. L.↗
NASA NTRS1991

A parametric study of the solar wind interaction with comets

The Naval Research Laboratory's magnetohydrodynamic simulation code is used to simulate the solar wind interction with comet Halley for two different outgassing rates and several different solar wind states. The magnetic field is more strongly draped for fast solar wind conditions than slow. For higher mass loading rates, the tail becomes wider and contains more magnetic flux. The visual appearance of the comet differs for the case in which the interplanetary magnetic field lies in the plane of the sky from the case when it lies along the line of sight. The ion tail appears shorter in the latter case. Thus variation in the IMF direction can cause significant changes in the appearance of comets. The comet also creates a large momentum flux deficit in the solar wind with a narrow enhanced region within it corresponding to the ion tail.

Russell, C. T.↗
NASA NTRS1992

Anomalous aspects of magnetosheath flow and of the shape and oscillations of the magnetopause during an interval of strongly northward interplanetary magnetic field

On 15 Feb. 1978, the orientation of the interplanetary magnetic field (IMF) remained steadily northward for more than 12 hours. The ISEE 1 and 2 spacecraft were located near apogee on the dawn side flank of the magnetotail. IMP 8 was almost symmetrically located in the magnetosheath on the dusk flank and IMP 7 was upstream in the solar wind. Using plasma and magnetic field data, we show the following: (1) the magnetosheath flow speed on the flanks of the magnetotail steadily exceeded the solar wind speed by 20 percent; (2) surface waves with approximately a 5-min period and very non-sinusoidal waveform were persistently present on the dawn magnetopause and waves of similar period were present in the dusk magnetosheath; and (3) the magnetotail ceased to flare at an antisunward distance of 15 R(sub E). We propose that the acceleration of the magnetosheath flow is achieved by magnetic tension in the draped field configuration for northward IMF and that the reduction of tail flaring is consistent with a decreased amount of open magnetic flux and a larger standoff distance of the subsolar magnetopause. Results of a three-dimensional magnetohydrodynamic simulation support this phenomenological model.

Chen, Sheng-Hsien↗
NASA NTRS1993

Numerical modeling of the energy storage and release in solar flares

This paper reports on investigation of the photospheric magnetic field-line footpoint motion (usually referred to as shear motion) and magnetic flux emerging from below the surface in relation to energy storage in a solar flare. These causality relationships are demonstrated by using numerical magnetohydrodynamic simulations. From these results, one may conclude that the energy stored in solar flares is in the form of currents. The dynamic process through which these currents reach a critical value is discussed as well as how these currents lead to energy release, such as the explosive events of solar flares.

Wu, S. T.↗
NASA NTRS1994

Overview of magnetospheric research

During the funding period for NASA Grant NAG5-1480 which prior to December, 1, 1990 was known as NASA Grant NAGW-78, the group has made substantial progress on the various topics originally proposed. The research performed has resulted in two Ph.D. Theses, more than 50 refereed papers in various journals and conference proceedings, and 31 invited and 104 contributed talks at conferences and symposia throughout the world. The main results from this work are summarized in each of the sections outlined in the original proposal, followed by a complete list of the group publications associated with this grant, a list of all invited talks given during the last three years, and finally a listing of the contributed talks. Research topics include global magnetohydrodynamic simulations, structure of the dayside magnetopause, merging at the dayside magnetopause, polar wind, plasma waves in the distant magnetotail, slow shocks in the distant magnetotail, collisionless tearing instabilities in the magnetotail, and other problems.

Ashour-Abdalla, Maha↗
NASA NTRS1992

Magnetotail instability

The stability of the geomagnetic tail is investigated on the basis of three dimensional resistive magnetohydrodynamic simulations, using different dynamic constraints and different initial equilibria. Different forms of the energy equation for isotropic pressure are found to have no significant effect on the dynamic growth of a resistive tearing instability, which is responsible for near Earth reconnection, plasmoid formation and ejection, and the generation of fast plasma flows. The constraints of a modified double adiabatic approach, however, can quench the tearing instability through the development of large, mirror type, anisotropies in the boundary regions of the plasma sheet, unless isotropization occurs on fast, nearly Alfvenic, time scales. The presence of a net cross tail magnetic field component B(sub yN) can reduce the growth of the instability without complete stabilization. An increase of B(sub z) from midnight toward the tail flanks, however, by more than a factor of about 3, apparently completely stabilizes the tearing mode. Stabilization and destabilization thus may depend on properties and constraints (and their release) in regions other than the neutral sheet where reconnection is initiated.

Birn, Joachim↗
NASA NTRS1994

The deformation of flux tubes in the solar wind with applications to the structure of magnetic clouds and CMEs

Two dimensional magnetohydrodynamic simulations of the distortion of a magnetic flux tube, accelerated through ambient solar wind plasma, are presented. Vortices form on the trailing edge of the flux tube, and couple strongly to its interior. If the flux tube azimuthal field is weak, it deforms into an elongated banana-like shape after a few Alfven transit times. A significant azimuthal field component inhibits this distortion. In the case of magnetic clouds in the solar wind, it is suggested that the shape observed at 1 AU was determined by distortion of the cloud in the inner heliosphere. Distortion of the cloud beyond 1 AU takes many days. It is estimated that effective drag coefficients slightly greater than unity are appropriate for modeling flux tube propagation. Synthetic magnetic field profiles as would be seen by a spacecraft traversing the cloud are presented.

Cargill, Peter J.↗