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 469 records · Page 26

Coronal evolution due to shear motion

Numerical solutions of the compressible MHD equations are used here to simulate the evolution of an initially force-free magnetic field in a static corona as a result of slow photospheric motion of the magnetic field footpoints. Simulations have been completed for values of plasma beta from 0.1 to 0.5, maximum shear velocities from 0.5 to 10.3 km/s, and with various amounts of resistive and viscous dissipation. In all cases the evolution proceeds in two qualitatively different stages. In the earlier stage, the field evolves gradually with the field lines, expanding outward at a velocity not unlike the shear velocity. Then, the field begins to expand much more rapidly until it reaches velocities exceeding a characteristic Alfven velocity. Inclusion of the thermodynamics, gravity, and compressibility is shown to have only a quantitative effect on the onset of the eruptive phase, illustrating that the primary interactions are between the dynamics and the magnetic field evolution.

Steinolfson, R. S.↗

Effect of high-latitude ionospheric convection on Sun-aligned polar caps

A coupled magnetospheric-ionospheric (M-I) magnetohydrodynamic (MHD) model has been used to simulate the formation of Sun-aligned polar cap arcs for a variety of interplanetary magnetic field (IMF) dependent polar cap convection fields. The formation process involves launching an Alfven shear wave from the magnetosphere to the ionosphere where the ionospheric conductance can react self-consistently to changes in the upward currents. We assume that the initial Alfven shear wave is the result of solar wind-magnetosphere interactions. The simulations show how the E region density is affected by the changes in the electron precipitation that are associated with the upward currents. These changes in conductance lead to both a modified Alfven wave reflection at the ionosphere and the generation of secondary Alfven waves in the ionosphere. The ensuing bouncing of the Alfven waves between the ionosphere and magnetosphere is followed until an asymptotic solution is obtained. At the magnetosphere the Alfven waves reflect at a fixed boundary. The coupled M-I Sun-aligned polar cap arc model of Zhu et al.(1993a) is used to carry out the simulations. This study focuses on the dependence of the polar cap arc formation on the background (global) convection pattern. Since the polar cap arcs occur for northward and strong B(sub y) IMF conditions, a variety of background convection patterns can exist when the arcs are present. The study shows that polar cap arcs can be formed for all these convection patterns; however, the arc features are dramatically different for the different patterns. For weak sunward convection a relatively confined single pair of current sheets is associated with the imposed Alfven shear wave structure. However, when the electric field exceeds a threshold, the arc structure intensifies, and the conductance increases as does the local Joule heating rate. These increases are faster than a linear dependence on the background electric field strength. Furthermore, above the threshold, the single current sheet pair splits into multiple current sheet pairs. For the fixed initial ionospheric and magnetospheric conditions used in this study, the separation distance between the current pairs was found to be almost independent of the background electric field strength. For either three-cell or distorted two-cell background convection patterns the arc formation favored the positive B(sub y) case in the northern hemisphere.

Sojka, J. J.↗

Resonant absorption of phase-mixed Alfven surface waves in ideal and resistive magnetohydrodynamics - Initial-value problem

An initial-value approach to the study of Alfven surface waves, in which the linearized MHD equations are solved numerically in time and space, is considered. A disturbance of a particular wavelength is excited in a nonuniform plasma and its temporal evolution is simulated with numerical solutions of the ideal and resistive MHD equations in slab geometry. This procedure generates a continuous spectrum and therefore simulates phase-mixing and subsequent resonant absorption. Energy in the phase-mixed surface waves resistively dissipates with the absorption time and width scaling as resistivity to the -1/3 and 1/6 powers, respectively.

Steinolfson, R. S.↗

Steepening of parallel propagating hydromagnetic waves into magnetic pulsations - A simulation study

The steepening mechanism of parallel propagating low-frequency MHD-like waves observed upstream of the earth's quasi-parallel bow shock has been investigated by means of electromagnetic hybrid simulations. It is shown that an ion beam through the resonant electromagnetic ion/ion instability excites large-amplitude waves, which consequently pitch angle scatter, decelerate, and eventually magnetically trap beam ions in regions where the wave amplitudes are largest. As a result, the beam ions become bunched in both space and gyrophase. As these higher-density, nongyrotropic beam segments are formed, the hydromagnetic waves rapidly steepen, resulting in magnetic pulsations, with properties generally in agreement with observations. This steepening process operates on the scale of the linear growth time of the resonant ion/ion instability. Many of the pulsations generated by this mechanism are left-hand polarized in the spacecraft frame.

Akimoto, K.↗

Magnetohydrodynamic Power Generation in the Laboratory Simulated Martian Entry Plasma

This paper addresses the magnetohydrodynamic (MHD) conversion of the energy released during the planetary entry phase of an interplanetary vehicle trajectory. The effect of MHD conversion is multi-fold. It reduces and redirects heat transferred to the vehicle, and regenerates the dissipated energy in reusable and transportable form. A vehicle on an interplanetary mission carries about 10,000 kWh of kinetic energy per ton of its mass. This energy is dissipated into heat during the planetary atmospheric entry phase. For instance, the kinetic energy of Mars Pathfinder was about 4220 kWh. Based on the loss in velocity, Mars Pathfinder lost about 92.5% of that energy during the plasma-sustaining entry phase that is approximately 3900 kWh. An ideal MHD generator, distributed over the probe surface of Mars Pathfinder could convert more than 2000 kWh of this energy loss into electrical energy, which correspond to more than 50% of the kinetic energy loss. That means that the heat transferred to the probe surface can be reduced by at least 50% if the converted energy is adequately stored, or re-radiated, or directly used. Therefore, MHD conversion could act not only as the power generating, but also as the cooling process. In this paper we describe results of preliminary experiments with light and microwave emitters powered by model magnetohydrodynamic generators and discuss method for direct use of converted energy.

Vuskovic, L.↗

The solar-flare induced earth's environment

A composite numerical simulation model developed from a series of MHD models was used to compute the solar-flare-generated disturbances of physical parameters, such as density, temperature, velocity, and magnetic field from the solar surface (i.e., the photospheric level) to the earth's environment. It is shown that the disturbed earth's environment at high latitudes can be approximated by starting with the knowledge of the occurrence and the strength of a solar flare, then simulating the evolutionary consequences of the solar disturbance through interplanetary space up to and through the magnetosphere.

Wu, S. T.↗

Modeling Soft X‐Ray Emissions at the Dayside Magnetopause

In this study, we simulate the Solar Wind Charge Exchange (SWCX) soft X-ray emissions at dayside magnetosheath and cusps by using magnetohydrodynamic (MHD) and LAtmos TEst Particle (LaTeP) models. MHD models are unable to resolve the particle kinetic effects, such as the different behaviors of ions with different q/m, or distinguish the magnetospheric plasma from the solar wind plasma. We investigate these effects with the LaTeP model. As the LaTeP model does not self-compute magnetic and electric field, the magnetic and electric field data obtained from Open Geospace General Circulation Model (OpenGGCM) and Lagrangian version of the piecewise parabolic method (PPMLR) MHD model are used as the input to LaTeP model. The soft X-ray emissivity maps simulated from pure OpenGGCM and PPMLR MHD approaches and from LaTeP-OpenGGCM and LaTeP-PPMLR approaches are presented and compared. The results indicate that the LaTeP model can well resolve the kinetic effects and can be used to investigate the individual spectral characteristics. Therefore, the LaTeP model is a complementary approach for simulating the X-ray emissions near the dayside magnetopause. We also calculate the ratio of integrated OVII/OVIII line intensities, produced by charge exchange of O7+ ions and O8+ ions, respectively. We find a relatively higher ratio at the bow shock compared to the surrounding areas, suggesting that this ratio can be an effective parameter to identify the bow shock location.

Qiuyu Xu↗

Coronal Heating, Weak MHD Turbulence, and Scaling Laws

Long-time high-resolution simulations of the dynamics of a coronal loop in Cartesian geometry are carried out, within the framework of reduced magnetohydrodynamics (RMHD), to understand coronal heating driven by the motion of field lines anchored in the photosphere. We unambiguously identify MHD anisotropic turbulence as the physical mechanism responsible for the transport of energy from the large scales, where energy is injected by photospheric motions, to the small scales, where it is dissipated. As the loop parameters vary, different regimes of turbulence develop: strong turbulence is found for weak axial magnetic fields and long loops, leading to Kolmogorov-like spectra in the perpendicular direction, while weaker and weaker regimes (steeper spectral slopes of total energy) are found for strong axial magnetic fields and short loops. As a consequence we predict that the scaling of the heating rate with axial magnetic field intensity B, which depends on the spectral index of total energy for given loop parameters, must vary from B3/2 for weak fields to B2 for strong fields at a given aspect ratio. The predicted heating rate is within the lower range of observed active region and quiet-Sun coronal energy losses.

corona↗

Evolution and interaction of interplanetary shocks

Whang's unsteady, one-dimensional, one-fluid MHD model is used to carry out a simulation study of the evolution of the solar wind based on two shock events over a distance of the order of 10 AU in the outer heliosphere. The study is based on the observation of two events, each observed by two or more spacecraft. The results show that the shock process, including the formation, collision, and merging of shocks, dominates the dynamical evolution of large-scale solar wind structures. In the outer heliosphere, the large-scale solar wind and magnetic field evolve into a much simpler structure, and MHD shocks are present as a principal component of the solar wind. The simulation results shed new light on the interaction and evolution of large interplanetary streams.

Whang, Y. C.↗

Space Weather Modeling and Prediction for Intermediate Time-scales

Much of solar activity within a sunspot cycle occurs as bursts, or 'seasons' of strong activity over several months, separated by periods of much less activity. The most important space weather effects occur during these bursts. Previous modeling and forecasting efforts have focused on time-scales of hours-to-days and decades-to-centuries. The recent discovery of Rossby waves in the Sun, together with recently developed global models of solar MHD Rossby waves and their interactions with differential rotation and spot-producing magnetic fields, reveal the opportunity to simulate and predict the occurrence, strength and location of enhanced activity bursts a few weeks up to several months in advance. We now have a golden opportunity to fill in this gap in time-scales of forecasting space weather. This requires a) continuous observations of solar Rossby waves by various techniques; b) development of coupled nonlinear MHD models that simulate both global Rossby waves and the much smaller spatial scale emergence of new active regions; c) application of advanced data assimilation techniques to couple surface observations to update the model-system to integrate forward in time for creating forecasts months ahead. Then it will be possible to build operational prediction models to meet the needs of customers and stakeholders, including support of future NASA missions, regarding what kind and level of space weather to expect a few weeks to several months ahead.

SMD↗

The effects on the earth's magnetotail from shocks in the solar wind

A new, minimally diffusive MHD computer code has been used to simulate the time-dependent, 2D interaction of the solar wind with the earth's magnetosphere. The equilibrium with a simulated unmagnetized solar wind shows the basic features of the magnetosphere: bow shock, magnetopause, and long magnetotail. A simulated shock in the solar wind shows the presence of shock focussing in the midplane of the magnetotail. This focussing can be important in the energy budget of the plasma sheet and may have an important role in the initiation of substorms.

Lyon, J.↗

Global simulations of the three-dimensional magnetosphere

Global three-dimensional computer simulations of the magnetosphere using a particle MHD code, reproduce the steady-state Dungey magnetospheric topology in three dimensions. The formation of a compression zone downstream of the tail neutral line that is probably bounded by wake shocks is observed. This compression zone changes its cross-section with distance downstream.

Leboeuf, J. N.↗

Theory and simulation of a high-frequency magnetic drift wave

The equilibrium of a nonneutral plasma in a toroidal vessel with a toroidal magnetic field is analyzed. In the zero inertia limit it is heuristically shown from force balance considerations that there is an electrostatic hoop force and a force due to diamagnetism along the major radius. The problem of equilibrium is formulated in terms of solutions of a 2D partial difference equation. This equation is solved in the large-aspect-ratio limit and a general expression for the shift of the potential axis is obtained which shows that the shift is approximately epsilon and that it depends solely on the internal capacitance of the cloud. The simulation study is based upon the modified MHD equations and the nonlocal nature of the mode is investigated. Applications to sub-Alfvenic plasma expansions, electromagnetic waves in the earth's magnetosphere, and plasma switches are discussed.

Huba, J. D.↗

Anisotropic turbulence in the solar wind

Solar wind turbulence has been viewed traditionally as composed of parallel propagating ('slab' fluctuations) or otherwise as isotropic turbulence. A variety of recent investigations, reviewed here, indicate that the spectrum may contain a significant admixture of two dimensional fluctuations, having variations mainly perpendicular to the local magnetic field. These indications come from simulations, from the theory of nearly incompressible MHD, from cosmic ray transport studies and from transport theory for solar wind turbulence, as well as from interpretations of direct observations. Thus, solar wind turbulence may be more like bundles of spaghetti than like parallel phase fronts.

Matthaeus, W. H.↗

First Lunar Wake Passage of ARTEMIS: Discrimination of Wake Effects and Solar Wind Fluctuations by 3D Hybrid Simulations

The spacecraft P1 of the new ARTEMIS (Acceleration, Reconnection, Turbulence, and Electrodynamics of the Moon's Interaction with the Sun) mission passed the lunar wake for the first time on February 13, 2010. We present magnetic field and plasma data of this event and results of 3D hybrid simulations. As the solar wind magnetic field was highly dynamic during the passage, a simulation with stationary solar wind input cannot distinguish whether distortions were caused by these solar wind variations or by the lunar wake; therefore, a dynamic real-time simulation of the flyby has been performed. The input values of this simulation are taken from NASA OMNI data and adapted to the P1 data, resulting in a good agreement between simulation and measurements. Combined with the stationary simulation showing non-transient lunar wake structures, a separation of solar wind and wake effects is achieved. An anisotropy in the magnitude of the plasma bulk flow velocity caused by a non-vanishing magnetic field component parallel to the solar wind flow and perturbations created by counterstreaming ions in the lunar wake are observed in data and simulations. The simulations help to interpret the data granting us the opportunity to examine the entire lunar plasma environment and, thus, extending the possibilities of measurements alone: A comparison of a simulation cross section to theoretical predictions of MHD wave propagation shows that all three basic MHD modes are present in the lunar wake and that their expansion governs the lunar wake refilling process.

Wiehle, S.↗

Finding the magnetopause location using soft X-ray observations and a statistical inverse method

Variability in the location and shape of the dayside magnetopause is attributed to magnetic reconnection, a fundamental process that enables the transfer of mass, energy, and momentum from the solar wind into the magnetosphere. The spatial and temporal properties of the magnetopause, under varying solar and magnetospheric conditions, remain largely unknown because empirical studies using in-situ observations are challenging to interpret. Global wide field-of-view (FOV) imaging is the only means to simultaneously observe the spatial distribution of the plasma properties over the vast dayside magnetospheric region and, subsequently, quantify the energy transport from the interplanetary medium into the terrestrial magnetosphere. Two upcoming missions, ESA/CAS SMILE and NASA’s LEXI will provide wide-field imagery of the dayside magnetosheath in soft X-rays, an emission generated by charge exchange interactions between high charge-state heavy ions of solar wind origin and exospheric neutral atoms. High-cadence two-dimensional observations of the magnetosheath will allow the estimation of dynamic properties of its inner boundary, the magnetopause, and enable studies of its response to changes in the solar wind dynamic pressure and interplanetary magnetic field orientation. This work introduces a statistically-based estimation approach based on inverse theory to estimate the spatial distribution of magnetosheath soft X-ray emissivities and, with this, identify the location of the magnetopause over the Sun−Earth line. To do so, we simulate the magnetosheath structure using the MHD-based OpenGGCM model and generate synthetic soft X-ray images using LEXI’s orbit and attitude information. Our results show that 3-D estimations using the described statistically-based technique are robust against Poisson-distributed shot noise inherent to soft X-ray images. Also, our proposed methodology shows that the accuracy of both three-dimensional (3-D) estimation and the magnetopause standoff distance calculation highly depends on the observational point.

Gonzalo Cucho-Padin↗

Expected IPS variations due to a disturbance described by a 3-D MHD model

The variations of interplanetary scintillation due to a disturbance described by a three-dimensional, time-dependent, MHD model of the interplanetary medium are calculated. The resulting simulated IPS maps are compared with observations of real disturbances and it is found that there is some qualitative agreement. It is concluded that the MHD model with a more realistic choice of input conditions would probably provide a useful description of many interplanetary disturbances.

Tappin, S. J.↗