Search NASA⌕ Search

SEARCH · Search NASA

Results for “MHD simulation”

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 217 records · Page 12

The Distant Tail at 200 R(sub E): Comparison Between Geotail Observations and the Results from a Global Magnetohydrodynamic Simulation

This paper reports a comparison between Geotail observations of plasmas and magnetic fields at 200 R(sub E) in the Earth's magnetotail with results from a time-dependent, global magnetohydrodynamic simulation of the interaction of the solar wind with the magnetosphere. The study focuses on observations from July 7, 1993, during which the Geotail spacecraft crossed the distant tail magnetospheric boundary several times while the interplanetary magnetic field (IMF) was predominantly northward and was marked by slow rotations of its clock angle. Simultaneous IMP 8 observations of solar wind ions and the IMF were used as driving input for the MHD simulation, and the resulting time series were compared directly with those from the Geotail spacecraft. The very good agreement found provided the basis for an investigation of the response of the distant tail associated with the clock angle of the IMF. Results from the simulation show that the stresses imposed by the draping of magnetosheath field lines and the asymmetric removal of magnetic flux tailward of the cusps altered considerably the shape of the distant tail as the solar wind discontinuities convected downstream of Earth. As a result, the cross section of the distant tail was considerably flattened along the direction perpendicular to the IMF clock angle, the direction of the neutral sheet following that of the IMF. The simulation also revealed that the combined action of magnetic reconnection and the slow rotation of the IMF clock angle led to a braiding of the distant tail's magnetic field lines along the axis of the tail, with the plane of the braid lying in the direction of the IMF.

Berchem, J.↗

Magnetic Field Orientation Effects on the Standoff Distance of Earth's Bow Shock

Three-dimensional, global MHD simulations of solar wind flow onto a prescribed magnetopause obstacle are used to show that a bow shock's nose location a(sub s), and the relative subsolar magnetosheath thickness Delta(sub ms)/a(sub mp) are strong functions of the IMF cone angle theta (between v(sub sw) and B(sub sw)) and the Alfven Mach number M(sub A). For a given M(sub A) the shock is more distant for higher theta (restricted to the interval 0-90deg by symmetries), while a(sub s)/a(sub mp) and Delta(sub ms/a(sub mp) increase with decreasing M(sub A) for theta greater than or approximately 20deg but decrease with decreasing M(sub A) for theta approximately Odeg. Large differences in Delta(sub ms/a(sub mp) are predicted between theta = Odeg and 90deg at low M(sub A), with smaller differences remaining even at M(sub A) approximately 10. The theta = Odeg results confirm and extend the previous work of Spreiter and Rizzi [1974]. The simulations show that successful models for the subsolar shock location cannot subsume the dependences on M(sub A) and theta into a sole dependence on M(ms). Instead, they confirm a recent prediction [Cairns and Grabbe, 1994] that a(sub s)/a(sub mp) and Delta(sub ms)/a(sub mp) should depend strongly on theta and M(sub A) for M(sub A) less than or approximately 10 (as well as on other MHD variables). Detailed comparisons between theory and data remain to be done. However, preliminary comparisons show good agreement, with distant shock locations found for low M(sub A) and large theta greater than or approximately 45deg and closer locations found for theta less than or approximately 20deg even at M/A approximately 8.

Cairns, Iver H.↗

The Production of Jets From Magnetic Accretion Disks: Simulation of the Blandford-Payne Mechanism

We have performed magnetohydrodynamic (MRD) simulations of the production of jets from magnetized accretion disks with a factor of 5 greater extent in space and time, and with more models, than any study published so far. We find that jets are produced by such disks in a broad range of parameter space, and by at least two different mechanisms. We also are able to follow the propagation of the jet well beyond the accretion disk into the region of hydrodynamic collimation. The code used is our MHD simulation code FLOW (K. Lind, D. Payne, D. Meier, and R. Blandford, 1989), converted to run on Caltech's massively parallel Intel Touchstone Delta supercomputer. Some of these models may be directly applicable to observed radio sources.

Blandford-Payne↗

Observation and simulation of the radial evolution and stream structure of solar wind turbulence

This paper presents a brief overview of the observed evolution in a variety of quantities describing the turbulent evolution of the interplanetary plasma and describes simulation results consistent with many features of the evolution. The turbulence is manifested through a dissipation at small scales in the inner heliosphere with a corresponding evolution in the breakpoint between a relatively flat and a Kolmogoroff spectrum; an evolution from kinetically to (slightly) magnetically dominated energy of the plasma fluctuations; a general decrease in the cross helicity or 'Alfvenicity'; changes in the anisotropy of the fluctuations; and the increasing predominance of quasi-pressure-balanced structures in the compressive component of the fluctuations. MHD simulations with shear layers either side of a central current sheet show that even in the absence of compressibility the lack of a mean field along the direction of the main flow in the current sheet leads to rapid nonlinear evolution and the observed characteristics of 'Elsasser spectra' of the fields in the inner heliosphere. Adding compressibility to the simulations does not greatly change the 'incompressive' quantities but leads in addition to observed correlations between a measure of compression and other quantities.

Roberts, D. A.↗

How a Realistic Magnetosphere Alters the Polarizations of Surface, Fast Magnetosonic, and Alfvén Waves

System-scale magnetohydrodynamic (MHD) waves within Earth's magnetosphere are often understood theoretically using box models. While these have been highly instructive in understanding many fundamental features of the various wave modes present, they neglect the complexities of geospace such as the inhomogeneities and curvilinear geometries present. Here, we show global MHD simulations of resonant waves impulsively excited by a solar wind pressure pulse. Although many aspects of the surface, fast magnetosonic (cavity/waveguide), and Alfvén modes present agree with the box and axially symmetric dipole models, we find some predictions for large-scale waves are significantly altered in a realistic magnetosphere. The radial ordering of fast mode turning points and Alfvén resonant locations may be reversed even with monotonic wave speeds. Additional nodes along field lines that are not present in the displacement/velocity occur in both the perpendicular and compressional components of the magnetic field. Close to the magnetopause, the perpendicular oscillations of the magnetic field have the opposite handedness to the velocity. Finally, widely used detection techniques for standing waves, both across and along the field, can fail to identify their presence. We explain how all these features arise from the MHD equations when accounting for a non-uniform background field and propose modified methods that might be applied to spacecraft observations.

ULF Waves↗

The tilts of corotating interaction regions at mid heliographic latitudes

Ulysses observations of corotating interaction regions (CIRs) at mid heliographic latitudes have shown that the flow downstream of the forward shock (or wave) on the leading edge of a CIR generally turns northward and westward, while the flow downstream of the reverse wave on the trailing edge generally turns southward and eastward. These systematic flow deflections are a natural consequence of large scale pressure gradients associated with the CIRs, and indicate that the forward waves tend to propagate toward and across the equator with increasing heliocentric distance, while the reverse waves tend to propagate toward the pole. Recent determinations of CIR shock normals using the Ulysses magnetic field data appear to confirm these plasma results (Burton, private communication). Numerical simulations indicate that these effects (which imply that CIRs are systematically tilted in the north-south direction at mid latitudes) are a natural consequence of the tilt of the solar magnetic dipole axis relative to the solar rotation axis. The present work utilizes a variety of techniques to analyze the flow deflections observed within CIRs from which we can infer the overall orientations of the CIRs and the speeds and directions of propagation of the waves. Where possible, the observations are quantitatively compared with the results of 3-dimensional MHD simulations.

Riley, Pete↗

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

Shear-induced opening of the coronal magnetic field

This work describes the evolution of a model solar corona in response to motions of the footpoints of its magnetic field. The mathematics involved is semianalytic, with the only numerical solution being that of an ordinary differential equation. This approach, while lacking the flexibility and physical details of full MHD simulations, allows for very rapid computation along with complete and rigorous exploration of the model's implications. We find that the model coronal field bulges upward, at first slowly and then more dramatically, in response to footpoint displacements. The energy in the field rises monotonically from that of the initial potential state, and the field configuration and energy appraoch asymptotically that of a fully open field. Concurrently, electric currents develop and concentrate into a current sheet as the limiting case of the open field is approached. Examination of the equations shows rigorously that in the asymptotic limit of the fully open field, the current layer becomes a true ideal MHD singularity.

Wolfson, Richard↗

Turbulent magnetic reconnection

The results of two-dimensional spectral-method MHD simulations of the effect of turbulence on magnetic-reconnection processes are presented in extensive graphs and tables and characterized in detail. In a periodic sheet-pinch configuration, early-onset, persistent turbulence is found to occur when a low level of broadband fluctuations is introduced into the initial data. Features observed include small-scale unsteady coherent electric-current and vorticity structures, enhanced viscous and resistive dissipation, unsteady and spatially asymmetric fluid flow, large-scale and small-scale magnetic islands, multiple X points, large-amplitude electric-field fluctuations, and spectral characteristics indicating that the turbulence is broadband in both spatial directions and similar to homogeneous MHD turbulence. The importance of reconnection for studies of solar, magnetospheric, and laboratory plasmas in inhomogeneous magnetic fields is indicated.

Matthaeus, W. H.↗

3D Realistic Modeling of the Interaction of Quiet-Sun Magnetic Fields with the Chromosphere

High-resolution observations and 3D simulations suggest that a local dynamo operates near the surface and produces ubiquitous small-scale magnetic elements, thus contributing to the magnetic carpet in the photosphere and to the magnetic structure and dynamics of the solar atmosphere. It appears that the traditional mechanisms of chromospheric energy and mass transport by acoustic waves and shocks are likely to play a secondary role; instead, the primary drivers in the energetics and dynamics of the chromosphere and transition region are small-scale, previously unresolved, quiet-Sun magnetic fields. These fields appear as ubiquitous, rapidly changing (on the scale of a few seconds), tiny magnetic loops and magnetized vortex tubes. Questions then arise about their origin and dynamics in the chromosphere, their links to magnetic fields in the photosphere, and their role in the energy storage and exchange between subsurface layers and the chromosphere. In the talk we will present results of 3D radiative MHD simulations obtained with the StellarBox code and discuss the energetics and dynamical interlinks between the subphotospheric layers and low chromosphere, their effects on the structure of the chromosphere, and signatures of the fine-scale magnetic features in high-resolution spectro-polarimetric observations.

Quie-Sun↗

The Distant Magnetotail Under Long Duration, Very Northward IMF Conditions: October 22-24, 2003

A unique 32 hour interval of very northward Interplanetary Magnetic Field (IMF) on October 22-24, 2003 created a exceptionally thick cold dense magnetotail plasma sheet, a small polar cap and accompanying small tail lobe. These features were detected by the Cluster DMSP and FAST spacecraft and modeled by a global simulation as described in papers by Oieroset et al. (2005) and Li et al. (2005). During the same interval the Wind spacecraft was passing through the center of the magnetotail about 130 Re downstream of Earth. Wind results will be described that reveal a very unusual magnetotail characterized by (1) continual tailward flow of 200-400 km/s with densities in the range 0.2-3/cc, both of whch are clearly less than those expected in the magnetosheath, (2) a mostly northward Bz but with a predominant Bx field component with sign reversals indicating crossings between the two hemispheres of the tail, and (3) velocity waves superposed on the downstream flow with peak-to-peak amplitudes of 100 to 200 km/s, periods of 10 to 20 minutes and clockwise polarization. Low altitude DMSP and Fast measurements reveal an auroral oval with enhanced latitudinal thickness and a small polar cap filled with structured precipitzting electrons and few ions. A new global MHD simulation of the event exhibits a highly elliptical tail of diminished cross-section at 130 Re with major axis aligned with the northward IMF. The tail current sheet also tends to be aligned in a north-south direction with the two tail hemispheres to the east and west with their polarities depending on prior history of the IMF. The simulation appears to be consistent with many, but not all, of the observations. High latitude cusp reconnection and subsequent downtail flow of closed field lines may explain the tail structure, but the waves are more likely due to the Kelvin-Helmholtz instability often thought to occur during northward IMF conditions.

Fairfield, Donald H.↗

Testing Secondary Models for the Origin of Radio Mini-Halos in Galaxy Clusters

We present an MHD simulation of the emergence of a radio minihalo in a galaxy cluster core in a "secondary" model, where the source of the synchrotron-emitting electrons is hadronic interactions between cosmic-ray protons with the thermal intracluster gas, an alternative to the "reacceleration model" where the cosmic ray electrons are reaccelerated by turbulence induced by core sloshing, which we discussed in an earlier work. We follow the evolution of cosmic-ray electron spectra and their radio emission using passive tracer particles, taking into account the time-dependent injection of electrons from hadronic interactions and their energy losses. We find that secondary electrons in a sloshing cluster core can generate diffuse synchrotron emission with luminosity and extent similar to observed radio minihalos. However, we also find important differences with our previous work. We find that the drop in radio emission at cold fronts is less prominent than that in our reacceleration-based simulations, indicating that in this flavor of the secondary model the emission is more spatially extended than in some observed minihalos. We also explore the effect of rapid changes in the magnetic field on the radio spectrum. While the resulting spectra in some regions are steeper than expected from stationary conditions, the change is marginal, with differences in the synchrotron spectral index of delta(spectral index) less than approximately 0.15-0.25, depending on the frequency band. This is a much narrower range than claimed in the best-observed minihalos and produced in the reacceleration model. Our results provide important suggestions to constrain these models with future observations.

clusters: general aEuro" magnetohydrodynamics (MHD↗

Particle acceleration by turbulent magnetohydro-dynamic reconnection

Test particles in a two dimensional, turbulent MHD simulation are found to undergo significant acceleration. The magnetic field configuration is a periodic sheet pinch which undergoes reconnection. The test particles are trapped in the reconnection region for times of order an Alfven transit time in the large electric fields that characterize the turbulent reconnection process at the relatively large magnetic Reynolds number used in the simulation. The maximum speed attained by these particles is consistent with an analytic estimate which depends on the reconnection electric field, the Alfven speed, and the ratio of Larmor period to the Alfven transit time.

Matthaeus, W. H.↗

Anomalous transport by magnetohydrodynamic Kelvin-Helmholtz instabilities in the solar wind-magnetosphere interaction

The high latitude, or downstream flank, and dayside low latitude magnetospheric boundaries are modeled in an MHD simulation of Kelvin-Helmholtz instablities in a compressible plasma for parallel and transverse configurations. Detailed attention is given to the nonlinear consequences of the instabilities for several different Alfven and sound Mach numbers in both configurations. Emphasis is given to the anomalous transport of momentum and energy by the Kelvin-Helmholtz instabilities across the magnetospheric boundary, which is crucial in evaluating the instabilities' contribution to magnetospheric convection. It is concluded that the anomalous tangential stress at the magnetospheric boundary caused by the instability may be of the order of 1 percent of the magnetosheath flow momentum adjacent to the boundary, and gives a contribution to the convection potential drop over the polar cap of 10-30 kV.

Miura, A.↗

Particle-acceleration by turbulent magnetohydrodynamic reconnection

Test particles in a two dimensional, turbulent MHD simulation are found to undergo significant acceleration. The magnetic field configuration is a periodic sheet pinch which undergoes reconnection. The test particles are trapped in the reconnection region for times of order an Alfven transit time in the large electric fields that characterize the turbulent reconnection process at the relatively large magnetic Reynolds number used in the simulation. The maximum speed attained by these particles is consistent with an analytic estimate which depends on the reconnection electric field, the Alfven speed, and the ratio of Larmor period to the Alfven transit time.

Matthaeus, W. H.↗

On the configuration of the polar cusps in earth's magnetosphere

The interaction between the solar wind and the earth's vacuum dipole field leads to the formation of a discontinuity called the magnetopause. In the standard picture, the magnetopause confines the magnetic field in such a manner that the polar cusp field lines originate from high latitudes in the dayside ionosphere and end at the two magnetic neutral points. Wu (1983, 1984) has questioned this standard picture of the polar cusp. MHD simulations indicate the existence of a current sheet above the polar cusp region, called 'the cusp current sheet' by Wu. Wu (1983) concluded that the difference between his cusp configuration and the standard picture is due to the fact that his geometry results from a plasma model, whereas the standard picture is based on a vacuum concept. In the present investigation, Wu's conclusion is questioned, and it is demonstrated that the standard cusp configuration is not restricted to the vacuum magnetosphere.

Voigt, G.-H.↗

The effect of different resistivity models on magnetotail dynamics

Two-dimensional, time-dependent MHD simulations of the interaction between the solar wind and the earth's magnetosphere have been performed to study magnetotail dynamics with varying forms of anomalous resistivity. In general, the resulting models conform to the neutral line model proposed for substorms with near-earth x points and high-speed tailward flows occurring in the magnetotail. However, in the case where Joule heating is included in the physical description of the system, the near-earth x point never moves far down tail, and high-speed tailward flows never stop. Only in the case where there is no Joule heating does the x point move down tail. Simultaneously, the high-speed tailward flows cease. These results indicate that the mechanism of energy dissipation can have an important effect on reconnection processes and the global magnetospheric dynamics.

Lyon, J. G.↗

Dynamic evolution of coronal magnetic fields

The response of coronal magnetic fields to photospheric motion is investigated using a time-dependent, two-dimensional MHD simulation. Starting with an initially uniform field, a circular section of the loop base is slowly rotated to represent the photospheric motion. The field lines at the base move with this flow in a manner consistent with the generated electric fields. The subsequent evolution of the field and flow can be characterized as passing through several distinct configurations. In the earliest phase the kinetic energy is negligible, and the current and field are parallel throughout most of the cylinder. This is followed by a period in which the field rotation increases, the axial field at and near the axis increases, and the acial field decreases in two cylindrical regions away from the axis. When the field in an appreciable portion of the cylinder has undergone one complete rotation, a rapid change in field configuration occurs with a large portion of the field making several rotations at large radii and a corresponding large reduction in the axial field.

Steinolfson, Richard S.↗