Search NASA⌕ Search

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.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

MHD Simulations of Earth's Bow Shock at low Mach Numbers: Standoff Distances

Global, three-dimensional, ideal MHD simulations of Earth's bow shock are reported for low Alfven Mach numbers M(sub A) and quasi-perpendicular magnetic field orientations. The simulations use a hard, infinitely conducting magnetopauause obstacle, with axisymmetric three-dimensional location given by scaled standard model, to directly address previous gasdynamic (GD) and field-aligned MHD (FA-MHD) work. Tests of the simulated shocks' density jumps X for 1.4 approx. less than MA approx. less than 10 and the high M(sub A) shock location, and reproduction of the GD relation between magnetosheath thickness and X for quasi-gasdynamic MHD runs with M(sub A) much greater than M(sub s), confirm that the MHD code is working correctly. The MHD simulations show the standoff distance a(sub s), increasing monotonically with decreasing M(sub A). Significantly larger a(sub s), are found at low M(sub A) than predicted by GD and phenomenological MHD models and FA-MHD simulations, as required qualitatively by observations. The GD and FA-MHD predictions err qualitatively, predicting either constant or decreasing a(sub s), with decreasing M(sub A). This qualitative difference between quasi- perpendicular MHD and FA-MHD simulations is direct evidence for a(sub s), depending on the magnetic field orientation Theta. The enhancement factor over the phenomenological MHD predictions at MA approx. 2.4 agrees quantitatively with one observatiorial estimate. A linear relationship is found between the magnetosheath thickness and X, modified both quantitatively and intrinsically by MHD effects from the GD result. The MHD and GD results agree in the high M(sub A) limit. An MHD theory is developed for a(sub s), restricted to sufficiently perpendicular Theta and high sonic Mach numbers M(sub s). It explains the simulation results with excellent accuracy. Observational and further simulation testing of this MHD theory, and of its predicted M(sub A), Theta, and M(sub s) effects, is desirable.

Cairns, Iver H.↗
NASA NTRS1995

MHD simulation of a propagation of loop-like and bubble-like magnetic clouds

Propagation and evolution of magnetic clouds in the ambient solar wind flow is studied self-consistently using ideal MHD equations in three dimensions. Magnetic clouds as ideal force-free objects (cylinders or spheres) are ejected near the Sun and followed beyond the Earth's orbit. We investigate the influence of various initial parameters like the injection velocity, magnetic field strength, magnetic helicity, orientation of the clouds' axis, etc., on their propagation and evolution. We demonstrate that the injection velocity and magnetic field strength have a major influence on propagation. Simulation results are compared with analytical solutions of magnetic cloud evolution.

Vandas, M.↗
NASA NTRS2014

Forced Reconnection in the Near Magnetotail: Onset and Energy Conversion in PIC and MHD Simulations

Using two-dimensional particle-in-cell (PIC) together with magnetohydrodynamic (MHD) Q1 simulations of magnetotail dynamics, we investigate the evolution toward onset of reconnection and the subsequent energy transfer and conversion. In either case, reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. The boundary deformation leads to the formation of thin embedded current sheets, which are bifurcated in the near tail, converging to a single sheet farther out in the MHD simulations. The thin current sheets in the PIC simulation are carried by electrons and are associated with a strong perpendicular electrostatic field, which may provide a connection to parallel potentials and auroral arcs and an ionospheric signal even prior to the onset of reconnection. The PIC simulation very well satisfies integral entropy conservation (intrinsic to ideal MHD) during this phase, supporting ideal ballooning stability. Eventually, the current intensification leads to the onset of reconnection, the formation and ejection of a plasmoid, and a collapse of the inner tail. The earthward flow shows the characteristics of a dipolarization front: enhancement of Bz, associated with a thin vertical electron current sheet in the PIC simulation. Both MHD and PIC simulations show a dominance of energy conversion from incoming Poynting flux to outgoing enthalpy flux, resulting in heating of the inner tail. Localized Joule dissipation plays only a minor role.

two-dimensional↗
NASA NTRS1991

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗
NASA NTRS1994

Oblique magnetohydrodynamic cosmic-ray-modified shocks: Two-fluid numerical simulations

We present the first results of time-dependent two-fluid cosmic-ray (CR) modified magnetohydrodynamic (MHD) shock simulations. The calculations were carried out with a new numerical code for one-dimensional ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks, including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they apppear to have converged to dynamical steady states. We also present results from an initial time-dependent simulation which extends the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescales than parallel shocks.

Frank, Adam↗
NASA NTRS1984

Comments on the MHD stability of coronal plasmas with line-tying

Reasonable boundary conditions and test functions for photospheric line-tying are discussed, and their effects on the ideal MHD stability of coronal loops are considered. It is concluded that the plasma at the footpoints of a loop may safely be assumed to be stationary. When a simple perturbed test function is used for the energy principle, the function should be helical in form. Constraints on the test function should be discarded because they overestimate stability.

An, C.-H.↗
NASA NTRS1994

MHD simulations of the magnetosheath

Ideal magnetohydrodynamic (MHD) simulations of the dayside magnetosheath are reported. The calculations show a depletion layer near the magnetopause and a new stagnation flow structure. The stagnation flow appears to be a generalization of the stagnation line for ideal MHD. For Mach numbers typical of the solar wind at the Earth, there is no evidence of the slow mode waves recently observed, but for low Mach numbers prominent slow mode waves are seen in the magnetosheath.

Lyon, J. G.↗
NASA NTRS2023

Modeling the Formation and Eruption of Coronal Structures by Linking Data-Driven Magnetofrictional and MHD Simulations for AR 12673

Context .The data-driven and time-dependent modeling of coronal magnetic fields is crucial for understanding solar eruptions. These efforts are complicated by the challenges of finding a balance between physical realism and computing efficiency. One possible technique is to couple two modeling approaches. Aims . Our aim here is to showcase our progress in using time-dependent magnetofrictional model (TMFM) results as input to dynamical magnetohydrodynamic (MHD) simulations. However, due to the different evolution processes in these two models, using TMFM snapshots in an MHD simulation is nontrivial. We address these issues, both physically and numerically, discuss the incompatibility of the TMFM output to serve as the initial condition in MHD simulations, and show our methods of mitigating this. The evolution of the flux systems and the cause of the eruption are investigated. Methods . TMFM is a prevalent approach that has proven to be a very useful tool in the study of the formation of unstable structures in the solar corona. In particular, it is capable of incorporating observational data as initial and boundary conditions and requires shorter computational time compared to MHD simulations. To leverage the efficiency of data-driven TMFM and also to simulate eruptive events in the MHD framework, one can apply TMFM up to a certain time before the expected eruption(s) and then proceed with the simulation in the full or ideal MHD regime in order to more accurately capture the eruption process. Results . We show the results of a benchmark test case with a linked TMFM and MHD simulation to study the evolution of NOAA active region 12673. A rise of a twisted flux bundle through the MHD simulation domain is observed, but we find that the rate of the rise and the altitude reached depends on the time of the TMFM snapshot that was used to initialize the MHD simulation and the helicity injected into the system. The analysis suggested that torus instability and slip-running reconnection could play an important role in the eruption. Conclusions . The results show that the linkage of TMFM and zero-β MHD models can be successfully used to model the eruptive coronal magnetic fields.

F. Daei↗
NASA NTRS1985

An MHD model of the earth's magnetosphere

It is pointed out that the earth's magnetosphere arises from the interaction of the solar wind with the earth's geomagnetic field. A global magnetohydrodynamics (MHD) model of the earth's magnetosphere has drawn much attention in recent years. In this model, MHD equations are used to describe the solar wind interaction with the magnetosphere. In the present paper, some numerical aspects of the model are considered. Attention is given to the ideal MHD equations, an equation of state for the plasma, the model as an initial- and boundary-value problem, the shock capturing technique, computational requirements and techniques for global MHD modeling, a three-dimensional mesh system employed in the global MHD model, and some computational results.

Wu, C. C.↗
NASA NTRS1984

Numerical simulation of reconnection in an emerging magnetic flux region

The resistive MHD equations are numerically solved in two dimensions for an initial-boundary-value problem which simulates reconnection between an emerging magnetic flux region and an overlying coronal magnetic field. The solution involves both ideal-MHD and resistive-MHD processes, and the solution shows an evolution which is remarkably suggestive of the preflare, impulsive, and main phases of the flare-cycle.

Forbes, T. C.↗
NASA NTRS1983

MHD stability of incompressible coronal loops with radiative energy loss

Previous studies of the magnetohydrodynamic (MHD) stability of solar coronal loops have not taken into account the effects of radiative energy loss in the energy equation. However, since coronal loops continuously lose energy by radiation and heat conduction, it is important to understand how these energy loss mechanisms affect MHD stability. We investigate the problem assuming that a magnetic loop has cylindrical geometry. As a first step, stability is studied for a localized mode, and the result is applied to a specific equilibrium. We find that the radiative energy loss effect not only changes the growth rate of ideally unstable modes, but also alters the stability boundary predicted by ideal MHD theory.

An, C.-H.↗
NASA NTRS2016

Magnetohydrodynamic Turbulence and the Geodynamo

Recent research results concerning forced, dissipative, rotating magnetohydrodynamic (MHD) turbulence will be discussed. In particular, we present new results from long-time Fourier method (periodic box) simulations in which forcing contains varying amounts of magnetic and kinetic helicity. Numerical results indicate that if MHD turbulence is forced so as to produce a state of relatively constant energy, then the largest-scale components are dominant and quasistationary, and in fact, have an effective dipole moment vector that aligns closely with the rotation axis. The relationship of this work to established results in ideal MHD turbulence, as well as to models of MHD turbulence in a spherical shell will also be presented. These results appear to be very pertinent to understanding the Geodynamo and the origin of its dominant dipole component. Our conclusion is that MHD turbulence, per se, may well contain the origin of the Earth's dipole magnetic field.

Shebalin, John V.↗
NASA NTRS1995

Heating of coronal holes by the resonant absorption and dissipation of Alfven waves and its relation to solar wind acceleration

Coronal hole regions are well known sources of high-speed solar wind, however to account for the observed properties of the solar wind a source of momentum and heat must be included. Alfven waves were suggested as the possible source of heating that accelerates the solar wind. We investigate the propagation of the Alfven waves in coronal holes via numerical solution of the linearized 2-D resistive MHD equations in slab geometry. The Alfven waves are driven at the lower boundary of the coronal hole and propagate into the corona. The waves are reflected at the coronal hole boundary and part of the wave energy leaks out of the coronal hole. We compare the calculated wavelengths and the attenuation rate of the fast mode Alfven waves in the leaky waveguide formed by the coronal hole with the analytical ideal MHD solutions. The formation of resonance heating layers is found to occur when shear Alfven waves propagate in an inhomogeneous coronal hole. The heating is enhanced when fast mode waves couple to the shear Alfven waves. The narrow heating layers are formed near the location of the ideal resonance, which might occur near the coronal hole boundary for a nearly constant density coronal hole, surrounded by a higher density plasma. We investigate the dependence of the heating on the driver frequency, the Lundquist number, and on the heliocentric distance. and find that the low frequency Alfven waves can be an efficient source of heating at large distances from the Sun. We discuss the relation of our results to the observed properties of high-speed solar wind and coronal holes.

Ofman, L.↗
NASA NTRS1989

Nonlinear Parker instability of isolated magnetic flux in a plasma

The nonlinear evolution of the Parker instability in an isolated horizontal magnetic-flux sheet embedded in a two-temperature layer atmosphere is studied by using a two-dimensional MHD code. In the solar case, this two-layer model is regarded as a simplified abstraction of the sun's photosphere/chromosphere and its overlying much hotter (coronal) envelope. The horizontal flux sheet is initially located in the lower temperature atmosphere so as to satisfy magnetostatic equilibrium under a constant gravitational acceleration. Ideal MHD is assumed, and only perturbations with k parallel to the magnetic-field lines are investigated. As the instability develops, the gas slides down the expanding loop, and the evacuated loop rises as a result of enhanced magnetic buoyancy. In the nonlinear regime of the instability, both the rise velocity of a magnetic loop and the local Alfven velocity at the top of the loop increase linearly with height and show self-similar behavior with height as long as the wavelength of the initial perturbation is much smaller than the horizontal size of the computing domain.

Shibata, K.↗
NASA NTRS2012

Particle Acceleration in Dissipative Pulsar Magnetospheres

Pulsar magnetospheres represent unipolar inductor-type electrical circuits at which an EM potential across the polar cap (due to the rotation of their magnetic field) drives currents that run in and out of the polar cap and close at infinity. An estimate ofthe magnitude of this current can be obtained by dividing the potential induced across the polar cap V approx = B(sub O) R(sub O)(Omega R(sub O)/c)(exp 2) by the impedance of free space Z approx eq 4 pi/c; the resulting polar cap current density is close to $n {GJ} c$ where $n_{GJ}$ is the Goldreich-Julian (GJ) charge density. This argument suggests that even at current densities close to the GJ one, pulsar magnetospheres have a significant component of electric field $E_{parallel}$, parallel to the magnetic field, a condition necessary for particle acceleration and the production of radiation. We present the magnetic and electric field structures as well as the currents, charge densities, spin down rates and potential drops along the magnetic field lines of pulsar magnetospheres which do not obey the ideal MHD condition $E cdot B = 0$. By relating the current density along the poloidal field lines to the parallel electric field via a kind of Ohm's law $J = sigma E_{parallel}$ we study the structure of these magnetospheres as a function of the conductivity $sigma$. We find that for $sigma gg OmegaS the solution tends to the (ideal) Force-Free one and to the Vacuum one for $sigma 11 OmegaS. Finally, we present dissipative magnetospheric solutions with spatially variable $sigma$ that supports various microphysical properties and are compatible with the observations.

Kazanas, Z.↗
NASA NTRS1985

Resistive decay of Alfven waves in a non-uniform plasma

The effect of resistive dissipation on the propagation of an MHD disturbance in a nonuniform plasma is examined. The present analysis, based on a boundary-layer technique, shows the existence of resistive normal modes with complex eigenfrequencies. The real part of the eigenfrequency is associated with an oscillatory behavior and defines the location in space of the layer where resistivity is important. The dissipation mechanism is responsible for the damping of the wave, in contrast with previous works in which the ideal MHD theory was used.

Mok, Y.↗
NASA NTRS1990

Global simulation of the solar wind-magnetosphere interaction - The importance of its numerical validity

This paper examines the global interaction of the solar wind with the geomagnetic field, using a newly developed three-dimensional high-precision MHD simulation code which employs the fourth-order Runge-Kutta-Gill scheme in time and the direct finite-difference method in space, and which allows much less numerical error than the conventional code. The simulation reveals in detail the formation processes of the magnetosphere and the bow shock, as well as the plasma sheet formation process. It is shown that, contrary to the conventional understanding, no Kelvin-Helmholtz instability occurs along the magnetopause in an ideal MHD.

Watanabe, K.↗
NASA NTRS1990

Dynamical evolution of twisted magnetic flux tubes. I - Equilibrium and linear stability

The three-dimensional dynamical evolution of twisted magnetic flux tubes is studied using a time-dependent magnetohydrodynamic (MHD) model. The flux tubes are intended to model solar coronal loops, and include the stabilizing effect of photospheric line tying. The model permits the complete evolution of flux tubes to be followed self-consistently, including the formation, equilibrium, linear instability, and nonlinear behavior. Starting from an initial uniform background magnetic field, a twisted flux tube is created by the application of slow, localized photospheric vortex flows. The flux tube evolves quasi-statically through sequences of equilibria with increasing twist, until it becomes linearly unstable to an ideal MHD kink mode. In this paper, the equilibrium properties and the linear stability behavior are discussed. The application of the method to the uniform-twist, Gold-Hoyle field confirms the previous stability threshold for kink instability and provides estimates of the resulting growth rate.

Mikic, Zoran↗