SEARCH · Search NASA

Results for “ideal MHD”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Measurement of the rugged invariants of magnetohydrodynamic turbulence in the solar wind

Measurements of the total energy, cross helicity, and magnetic helicity of the solar wind at 1, 2.8, and 5 AU are presented. These quantities are the three rugged invariants of three-dimensional ideal incompressible MHD turbulence theory. The theoretical technique for measuring the magnetic helicity from the matrix of two-point correlations is shown. The length scales characterizing the magnetic helicity are found to be equal to or greater than those which characterize the magnetic energy. The magnetic helicity typically lies at scales larger than the magnetic correlation length, consistent with the expectations of the inverse cascade and selective decay hypotheses of three-dimensional MHD turbulence. At smaller scales, the magnetic helicity oscillates in sign. Measurements of the cross helicity are not fully consistent with the usual interpretation in terms of outward propagating Alfvenic functions. Especially during the interval at 5 AU the cross helicity is found to oscillate in sign indicating fluctuations propagating both outward and inward.

Matthaeus, W. H.↗

NASA NTRS1986

Filament formation due to photospheric shear

Recently, Wu et al. (1983, 1986) developed a three-dimensional, time-dependent axisymmetric ideal magnetohydrodynamic (MHD) model which has been used to study lower solar atmospheric response due to photospheric motion. They have shown that various observed solar features can be formed from these numerical simulation experiments. Some work currently in progress which has resulted from a specific numerical simulation experiment is reported. These results may elucidate the possible physical mechanism of formation of filament.

Wu, S. T.↗

NASA NTRS1986

Formation of active region and quiescent prominence magnetic field configurations

To investigate the formation of prominences, researchers studied chromospheric mass injection into an overlying coronal dipole magnetic field using a 2-D ideal magnetohydrodynamic (MHD) numerical model. Researchers propose that active region prominences are formed by chromospheric plasmas injected directly into the overlying coronal magnetic field and that quiescent prominences are formed by plasmas evaporated at the interface between spicules and corona. Hence, for the simulation of an active region prominence magnetic field we inject the mass from one side, but use a symmetric mass injection to form a quiescent prominence field configuration. Researchers try to find optimum conditions for the formation of Kippenhahn-Schuluter(K-S)type field configuration for stable support of the injection plasmas. They find that the formation of K-S type field configuration by mass injection requires a delicate balance between injection velocity, density, and overlying magnetic fields. These results may explain why a prominence does not form on every neutral line.

An, C.-H.↗

NASA NTRS1990

Self-similar evolution of the nonlinear magnetic buoyancy instability

A new type of self-similar solution of ideal magnetohydrodynamics (MHD) in the nonlinear stage of the undular model (k parallel to B) of the magnetic buoyancy instability (the ballooning instability in fusion plasma physics or the Parker instability in astrophysics) is found through MHD simulation and theory. The linear theory developed agrees well with the simulation in the early (linear) stage. The nonlinear stages of the instability in the simulation show the self-similar evolution. One of the solutions obtained from the nonlinear analysis has the characteristics of nonlinear instability in Lagrangian coordinates; the fluid velocity and the Alfven speed on each magnetic loop increase exponentially with time, because the loop is evacuated by the field-aligned motion of matter resulting from gravitational acceleration. In the later stage of the nonlinear evolution, the solution property changes from exponential to power-law time dependence. The latter corresponds to a force-free expansion solution. The later saturation of the velocity increment is also discussed.

Shibata, K.↗

NASA NTRS1992

The effects of radiation drag on radial, relativistic hydromagnetic winds

The effects of drag on an idealized relativistic MHD wind of radial geometry are studied. The astrophysical motivation is to understand the effects of radiation drag on the dynamics of a jet or wind passing through the intense radiation field of an accreting compact object. From a critical point analysis, it is found that a slow magnetosonic point can appear in a dragged flow even in the absence of gravitational force, as a result of a balance between the drag force and the combination of thermal pressure and centrifugal forces. As in the undragged case, the Alfven point does not impose any constraints on the flow. Although it is formally possible for a dragged flow to possess more than one fast magnetosonic point, it is shown that this is unlikely in practice. In the limit of a 'cold', centrifugally driven flow, it is shown that the fast magnetosonic point moves to infinite radius, just as in the drag-free case. For a given mass flux, the total energy output carried to infinity, and the final partition between the kinetic energy and the Poynting flux, are the same for the dragged and the drag-free flows. The main effects of radiation drag are to increase the amount of energy and angular momentum extracted from the source and to redistribute the regions where acceleration occurs in the flow. This is accomplished through the storage and release of magnetic energy, as a result of additional winding and compression of the field caused by the action of the drag. For a relativistic wind, the dissipated energy can exceed the final kinetic energy of the flow and may be comparable to the total flow energy (which is dominated by Poynting flux). The energy lost to radiation drag will appear as a Doppler-boosted beam of scattered radiation, which could dominate the background radiation if the flow is well-collimated.

Li, Zhi-Yun↗

NASA NTRS1994

An approximate Riemann solver for magnetohydrodynamics (that works in more than one dimension)

An approximate Riemann solver is developed for the governing equations of ideal magnetohydrodynamics (MHD). The Riemann solver has an eight-wave structure, where seven of the waves are those used in previous work on upwind schemes for MHD, and the eighth wave is related to the divergence of the magnetic field. The structure of the eighth wave is not immediately obvious from the governing equations as they are usually written, but arises from a modification of the equations that is presented in this paper. The addition of the eighth wave allows multidimensional MHD problems to be solved without the use of staggered grids or a projection scheme, one or the other of which was necessary in previous work on upwind schemes for MHD. A test problem made up of a shock tube with rotated initial conditions is solved to show that the two-dimensional code yields answers consistent with the one-dimensional methods developed previously.

Powell, Kenneth G.↗

NASA NTRS1994

On the maximum energy release in flux-rope models of eruptive flares

We determine the photospheric boundary conditions which maximize the magnetic energy released by a loss of ideal-magnetohydrodynamic (MHD) equilibrium in two-dimensional flux-rope models. In these models a loss of equilibrium causes a transition of the flux rope to a lower magnetic energy state at a higher altitude. During the transition a vertical current sheet forms below the flux rope, and reconnection in this current sheet releases additional energy. Here we compute how much energy is released by the loss of equilibrium relative to the total energy release. When the flux-rope radius is small compared to its height, it is possible to obtain general solutions of the Grad-Shafranov equation for a wide range of boundary conditions. Variational principles can then be used to find the particular boundary condition which maximizes the magnetic energy released for a given class of conditions. We apply this procedure to a class of models known as cusp-type catastrophes, and we find that the maximum energy released by the loss of equilibrium is 20.8% of the total energy release for any model in this class. If the additional restriction is imposed that the photospheric magnetic field forms a simple arcade in the absence of coronal currents, then the maximum energy release reduces to 8.6%

Forbes, T. G.↗

NASA NTRS1994

Local shear instabilities in weakly ionized, weakly magnetized disks

We extend the analysis of axisymmetric magnetic shear instabilities from ideal magnetohydrodynamic (MHD) flows to weakly ionized plasmas with coupling between ions and neutrals caused by collisions, ionization, and recombination. As part of the analysis, we derive the single-fluid MHD dispersion relation without invoking the Boussinesq approximation. This work expands the range of applications of these instabilities from fully ionized accretion disks to molecular disks in galaxies and, with somewhat more uncertainty, to protostellar disks. Instability generally requires the angular velocity to decrease outward, the magnetic field strengths to be subthermal, and the ions and neutrals to be sufficiently well coupled. If ionization and recombination processes can be neglected on an orbital timescale, adequate coupling is achieved when the collision frequency of a given neutral with the ions exceeds the local epicyclic freqency. When ionization equilibrium is maintained on an orbital timescale, a new feature is present in the disk dynamics: in contrast to a single-fluid system, subthermal azimuthal fields can affect the axisymmetric stability of weakly ionized two-fluid systems. We discuss the underlying causes for this behavior. Azimuthal fields tend to be stabilizing under these circumstances, and good coupling between the neutrals and ions requires the collision frequency to exceed the epicyclic frequency by a potentially large secant factor related to the magnetic field geometry. When the instability is present, subthermal azimuthal fields may also reduce the growth rate unless the collision frequency is high, but this is important only if the field strengths are very subthermal and/or the azimuthal field is the dominant field component. We briefly discuss our results in the context of the Galactic center circumnuclear disk, and suggest that the shear instability might be present there, and be responsible for the observed turbulent motions.

Blaes, Omer M.↗

NASA NTRS1998

Gas-Kinetic Theory Based Flux Splitting Method for Ideal Magnetohydrodynamics

A gas-kinetic solver is developed for the ideal magnetohydrodynamics (MHD) equations. The new scheme is based on the direct splitting of the flux function of the MHD equations with the inclusion of "particle" collisions in the transport process. Consequently, the artificial dissipation in the new scheme is much reduced in comparison with the MHD Flux Vector Splitting Scheme. At the same time, the new scheme is compared with the well-developed Roe-type MHD solver. It is concluded that the kinetic MHD scheme is more robust and efficient than the Roe- type method, and the accuracy is competitive. In this paper the general principle of splitting the macroscopic flux function based on the gas-kinetic theory is presented. The flux construction strategy may shed some light on the possible modification of AUSM- and CUSP-type schemes for the compressible Euler equations, as well as to the development of new schemes for a non-strictly hyperbolic system.

Xu, Kun↗

NASA NTRS1998

Streamer Evaporation

Evaporation is the consequence of heating near the top of streamers in ideal Magnetohydrodynamics (MHD) models, where the plasma is weakly contained by the magnetic field. Heating causes slow opening of field lines and release of new solar wind. It was discovered in simulations and, due to the absence of loss mechanisms, the ultimate end point is the complete evaporation of the streamer. Of course streamers do not behave in this way because there are losses by thermal conduction and radiation. Physically, heating is also expected to depend on ambient conditions. We use our global MHD model with thermal conduction to examine the effect of changing the heating scale height. We also apply and extend an analytic model of streamers developed by Pneuman (1968) to show that steady streamers are unable to contain plasma for temperatures near the cusp greater than approximately 2 x 10(exp 6) K.

Suess, S. T.↗

NASA NTRS2002

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↗

NASA NTRS2007

Broken Symmetries and Magnetic Dynamos

Phase space symmetries inherent in the statistical theory of ideal magnetohydrodynamic (MHD) turbulence are known to be broken dynamically to produce large-scale coherent magnetic structure. Here, results of a numerical study of decaying MHD turbulence are presented that show large-scale coherent structure also arises and persists in the presence of dissipation. Dynamically broken symmetries in MHD turbulence may thus play a fundamental role in the dynamo process.

Shebalin, John V.↗

NASA NTRS2013

A Simple GPU-Accelerated Two-Dimensional MUSCL-Hancock Solver for Ideal Magnetohydrodynamics

We describe our experience using NVIDIA's CUDA (Compute Unified Device Architecture) C programming environment to implement a two-dimensional second-order MUSCL-Hancock ideal magnetohydrodynamics (MHD) solver on a GTX 480 Graphics Processing Unit (GPU). Taking a simple approach in which the MHD variables are stored exclusively in the global memory of the GTX 480 and accessed in a cache-friendly manner (without further optimizing memory access by, for example, staging data in the GPU's faster shared memory), we achieved a maximum speed-up of approx. = 126 for a sq 1024 grid relative to the sequential C code running on a single Intel Nehalem (2.8 GHz) core. This speedup is consistent with simple estimates based on the known floating point performance, memory throughput and parallel processing capacity of the GTX 480.

graphics processing units↗

NASA NTRS2007

Multi-Scale Modeling of Magnetospheric Reconnection

One of the major challenges in modeling the magnetospheric magnetic reconnection is to quantify the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions near reconnection sites. There is still considerable debate as to what degree microphysical processes on kinetic scales affect the global evolution and how important it is to substitute numerical dissipation and/or ad hoc anomalous resistivity by a physically motivated model of dissipation. Comparative studies of magnetic reconnection in small scale geometries demonstrated that MHD simulations that included non-ideal processes in terms of a resistive term $\eta J$ did not produce the fast reconnection rates observed in kinetic simulations. For a broad range of physical parameters in collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. We utilize the global MHD code BATSRUS and incorporate nongyrotropic effects in diffusion regions in terms of corrections to the induction equation. We developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated selfconsistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites is updated during the simulations. To clarify the role of nongyrotropic effects in diffusion region on the global magnetospheric dynamic we perform simulations with steady southward IMF driving of the magnetosphere. Ideal MHD simulations with magnetic reconnection supported by numerical resistivity produce steady configuration with almost stationary near-earth neutral line (NENL). Simulations with non-gyrotropic corrections demonstrate dynamic quasi-periodic response to the steady driving condition. The loading/unloading cycle in non-gyrotropic MHD results has a non-stationary reconnection site in the magnetotail, with the retreating during the stretching phase and then a new NENL forming in the resulting thin plasma sheet. We expect that this model will lead to improved representations of space weather event in the magnetosphere.

Kuznetsova, M. M.↗

NASA NTRS1997

Design Study: Rocket Based MHD Generator

This report addresses the technical feasibility and design of a rocket based MHD generator using a sub-scale LOx/RP rocket motor. The design study was constrained by assuming the generator must function within the performance and structural limits of an existing magnet and by assuming realistic limits on (1) the axial electric field, (2) the Hall parameter, (3) current density, and (4) heat flux (given the criteria of heat sink operation). The major results of the work are summarized as follows: (1) A Faraday type of generator with rectangular cross section is designed to operate with a combustor pressure of 300 psi. Based on a magnetic field strength of 1.5 Tesla, the electrical power output from this generator is estimated to be 54.2 KW with potassium seed (weight fraction 3.74%) and 92 KW with cesium seed (weight fraction 9.66%). The former corresponds to a enthalpy extraction ratio of 2.36% while that for the latter is 4.16%; (2) A conceptual design of the Faraday MHD channel is proposed, based on a maximum operating time of 10 to 15 seconds. This concept utilizes a phenolic back wall for inserting the electrodes and inter-electrode insulators. Copper electrode and aluminum oxide insulator are suggested for this channel; and (3) A testing configuration for the sub-scale rocket based MHD system is proposed. An estimate of performance of an ideal rocket based MHD accelerator is performed. With a current density constraint of 5 Amps/cm(exp 2) and a conductivity of 30 Siemens/m, the push power density can be 250, 431, and 750 MW/m(sup 3) when the induced voltage uB have values of 5, 10, and 15 KV/m, respectively.

Source record↗

NASA NTRS2004

Theory and Simulation of Real and Ideal Magnetohydrodynamic Turbulence

Incompressible, homogeneous magnetohydrodynamic (MHD) turbulence consists of fluctuating vorticity and magnetic fields, which are represented in terms of their Fourier coefficients. Here, a set of five Fourier spectral transform method numerical simulations of two-dimensional (2-D) MHD turbulence on a 512(sup 2) grid is described. Each simulation is a numerically realized dynamical system consisting of Fourier modes associated with wave vectors k, with integer components, such that k = |k| less than or equal to k(sub max). The simulation set consists of one ideal (non-dissipative) case and four real (dissipative) cases. All five runs had equivalent initial conditions. The dimensions of the dynamical systems associated with these cases are the numbers of independent real and imaginary parts of the Fourier modes. The ideal simulation has a dimension of 366104, while each real simulation has a dimension of 411712. The real runs vary in magnetic Prandtl number P(sub M), with P(sub M) is a member of {0.1, 0.25, 1, 4}. In the results presented here, all runs have been taken to a simulation time of t = 25. Although ideal and real Fourier spectra are quite different at high k, they are similar at low values of k. Their low k behavior indicates the existence of broken symmetry and coherent structure in real MHD turbulence, similar to what exists in ideal MHD turbulence. The value of PM strongly affects the ratio of kinetic to magnetic energy and energy dissipation (which is mostly ohmic). The relevance of these results to 3-D Navier-Stokes and MHD turbulence is discussed.

Shebalin, John V.↗

NASA NTRS2002

Analysis of the Magneto-Hydrodynamic (MHD) Energy Bypass Engine for High-Speed Air-Breathing Propulsion

The performance of the MHD energy bypass air-breathing engine for high-speed propulsion is analyzed in this investigation. This engine is a specific type of the general class of inverse cycle engines. In this paper, the general relationship between engine performance (specific impulse and specific thrust) and the overall total pressure ratio through an engine (from inlet plane to exit plane) is first developed and illustrated. Engines with large total pressure decreases, regardless of cause or source, are seen to have exponentially decreasing performance. The ideal inverse cycle engine (of which the MHD engine is a sub-set) is then demonstrated to have a significant total pressure decrease across the engine; this total pressure decrease is cycle-driven, degrades rapidly with energy bypass ratio, and is independent of any irreversibility. The ideal MHD engine (inverse cycle engine with no irreversibility other than that inherent in the MHD work interaction processes) is next examined and is seen to have an additional large total pressure decrease due to MHD-generated irreversibility in the decelerator and the accelerator. This irreversibility mainly occurs in the deceleration process. Both inherent total pressure losses (inverse cycle and MHD irreversibility) result in a significant narrowing of the performance capability of the MHD bypass engine. The fundamental characteristics of MHD flow acceleration and flow deceleration from the standpoint of irreversibility and second-law constraints are next examined in order to clarify issues regarding flow losses and parameter selection in the MM modules. Severe constraints are seen to exist in the decelerator in terms of allowable deceleration Mach numbers and volumetric (length) required for meaningful energy bypass (work interaction). Considerable difficulties are also encountered and discussed due to thermal/work choking phenomena associated with the deceleration process. Lastly, full engine simulations utilizing inlet shock systems, finite-rate chemistry, wall cooling with thermally balanced engine (fuel heat sink), fuel injection and mixing, friction, etc. are shown and discussed for both the MHD engine and the conventional scramjet. The MHD bypass engine has significantly lower performance in all categories across the Mach number range (8 to 12.2). The lower performance is attributed to the combined effects of 1) additional irreversibility and cooling requirements associated with the MHD components and 2) the total pressure decrease associated with the inverse cycle itself.

Riggins, David W.↗

NASA NTRS2010

Coherent Eigenmodes in Homogeneous MHD Turbulence

The statistical mechanics of Fourier models of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence is discussed, along with their relevance for dissipative magnetofluids. Although statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation, i.e., we have coherent structure. We use eigenanalysis of the modal covariance matrices in the probability density function to 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 provide examples from 2-D and 3-D magnetohydrodynamic simulations of homogeneous turbulence, and show new results from long-time simulations of MHD turbulence with and without a mean magnetic field

Shebalin, John V.↗