Search NASA⌕ Search

SEARCH · Search NASA

Results for “MHD modeling”

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 19 records
All sourcesNASA NTRSDOE OSTI
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 NTRS2023

Is the Global MHD Modeling of the Magnetosphere Adequate for GIC Prediction: the May 27–28, 2017 Storm

Practical steps taken by the international community to reduce the damage to technological systems from space weather include the development of numerical models capable of real-time predictions of electromagnetic disturbances at the Earth’s surface. Here we examine the feasibility of a version of the Space Weather Modeling Framework (SWMF) global MHD simulation code similar to that used by the NOAA Space Weather Prediction Center to predict the level of geomagnetic field variability, and consequently geomagnetically induced currents (GICs). We consider the contribution of geomagnetic disturbances to the bursts of GIC in the electric power line of the Kola Peninsula during the May 27–28, 2017 storm and compare the observations with results of the global MHD model. During the maximal disturbance magnetic field variations at East Scandinavian stations become more irregular, as intense Pi3 pulsations are superposed on the magnetic bay. These pulsations are not quasi-sinusoidal waves like typical Pc5 pulsations, but they are rather a quasi-periodic sequence of magnetic impulses with time scales ~5–15 min. During this period with elevated Pi3 activity very high values of GIC were recorded (variations >100 A) in the electric power transmission line. The SWMF modeling reasonably well reproduces the global magnetospheric parameters, such as SYM-H index or cross-polar potential. However, the magnetic field variability dB/dt in the East Scandinavia predicted by the modeling has turned out to be more than order of magnitude less than that observed. Thus, the version of SWMF with the grid used by NOAA SWPC still cannot adequately predict for the May 27–28 event the fine structure of the storm/substorm—Pi3 geomagnetic disturbances, and consequently the magnitude of the GIC that they drive.

geomagnetically induced currents↗
NASA NTRS1986

Magnetohydrodynamics (MHD) modelling of flare energy buildup, the energy release phase, and its propagation into heliospheric space

Solar flare energy buildup at the photospheric level and energy release and transport into heliospheric space are examined using a composite MHD model. A four phase composite MHD model is described. An example demonstrating the applicability of the model is presented; the model was applied to the active region AR 2372. The limitations of this composite MHD model approach to analyzing solar flare energy buildup are discussed.

Wu, S. T.↗
NASA NTRS2004

3D MHD Models of Active Region Loops

Present imaging and spectroscopic observations of active region loops allow to determine many physical parameters of the coronal loops, such as the density, temperature, velocity of flows in loops, and the magnetic field. However, due to projection effects many of these parameters remain ambiguous. Three dimensional imaging in EUV by the STEREO spacecraft will help to resolve the projection ambiguities, and the observations could be used to setup 3D MHD models of active region loops to study the dynamics and stability of active regions. Here the results of 3D MHD models of active region loops are presented, and the progress towards more realistic 3D MHD models of active regions. In particular the effects of impulsive events on the excitation of active region loop oscillations, and the generation, propagations and reflection of EIT waves are shown. It is shown how 3D MHD models together with 3D EUV observations can be used as a diagnostic tool for active region loop physical parameters, and to advance the science of the sources of solar coronal activity.

Ofman, Leon↗
NASA NTRSDate not supplied

Observations and 3D MHD Modeling of Small-scale Solar Prominence Oscillations

Small scaleoscillations were recently observed in prominences with high resolution Hinode/SOT and IRIS SJI in unprecedented detailand identified as fast magnetosonic waves. The high-resolution observations provide evidence of flows of cool material, fluid instabilities such as Kelvin-Helmholtz and Rayleigh-Taylor, and the spectroscopic information provides evidence of Doppler shifts, heating, and cooling of the prominence material. Recently nonlinear fast magnetosonic waves were observed withperiods in the range of 5–11 minutes and wavelengths in the plane of the sky (POS) of ~2000kmand flows, and the flows in narrow threads withPOS speed ∼16–46 km/s. Thenonlinear fast magnetosonic waves and flows were modeled using 2.5D MHD (Ofman & Kucera 2020). Here, we extend the studyof the nonlinear fast magnetosonic waves using more realistic 3D MHD model. The model includes the 3D structure of the density and the magnetic fieldas well as mode coupling anddemonstrate the propagating density compressions associated with the nonlinear fast magnetosonicwaves, as well as the guided modesin the pillar. The 3DMHD modelingimprove understanding ofthe POSeffects in limb prominence oscillations observations.The results of the modeling are useful for coronal seismology of the prominence structure.

Leon Ofman↗
NASA NTRS2017

The Substorm Cycle as Reproduced by Global MHD Models

Recently, Gordeev et al. (2015) suggested a method to test global MHD models against statistical empirical data. They showed that four community-available global MHD models supported by the Community Coordinated Modeling Center (CCMC) produce a reasonable agreement with reality for those key parameters (the magnetospheric size, magnetic field, and pressure) that are directly related to the large-scale equilibria in the outer magnetosphere. Based on the same set of simulation runs, here we investigate how the models reproduce the global loading-unloading cycle. We found that in terms of global magnetic flux transport, three examined CCMC models display systematically different response to idealized2 h north then 2 h south interplanetary magnetic field (IMF) Bz variation. The LFM model shows a depressed return convection and high loading rate during the growth phase as well as enhanced return convection and high unloading rate during the expansion phase, with the amount of loaded unloaded magnetotail flux and the growth phase duration being the closest to their observed empirical values during isolated substorms. Two other models exhibit drastically different behavior. In the BATS-R-US model the plasma sheet convection shows a smooth transition to the steady convection regime after the IMF southward turning. In the Open GGCM a weak plasma sheet convection has comparable intensities during both the growth phase and the following slow unloading phase. We also demonstrate potential technical problem in the publicly available simulations which is related to post processing interpolation and could affect the accuracy of magnetic field tracing and of other related procedures.

Gordeev, E.↗
NASA NTRS2005

A Numerical Study of Resistivity and Hall Effects for a Compressible MHD Model

The effect of resistive, Hall, and viscous terms on the flow structure compared with compressible ideal MHD is studied numerically for a one-fluid non-ideal MHD model. The goal of the present study is to shed some light on the emerging area of non-ideal MHD modeling and simulation. Numerical experiments are performed on a hypersonic blunt body flow with future application to plasma aerodynamics flow control in reentry vehicles. Numerical experiments are also performed on a magnetized time-developing mixing layer with possible application to magnetic/turbulence mixing.

Yee, H. C.↗
NASA NTRS2001

Global MHD Models of the Solar Corona

Global magnetohydrodynamic (MHD) models of the solar corona are computationally intensive, numerically complex simulations that have produced important new results over the past few years. After a brief overview of how these models usually work, I will address three topics: (1) How these models are now routinely used to predict the morphology of the corona and analyze Earth and space-based remote observations of the Sun; (2) The direct application of these models to the analysis of physical processes in the corona and chromosphere and to the interpretation of in situ solar wind observations; and (3) The use of results from global models to validate the approximations used to make detailed studies of physical processes in the corona that are not otherwise possible using the global models themselves.

Suess, S. T.↗
NASA NTRS1983

Global MHD model of the earth's magnetosphere

A global MHD model of the earth's magnetosphere is defined. An introduction to numerical methods for solving the MHD equations is given with emphasis on the shock-capturing technique. Finally, results concerning the shape of the magnetosphere and the plasma flows inside the magnetosphere are presented.

Wu, C. C.↗
NASA NTRS1988

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

A New Axisymmetric MHD Model of the Interaction of the Solar Wind with Venus

A new two-dimensional axisymmetric MHD model is used to study the interaction of the solar wind with Venus under conditions where the interplanetary field is approximately aligned with the solar wind velocity. This numerical model solves the MHD transport equations for density, velocity, pressure, and magnetic field on an adaptively refined, unstructured grid system. This use of an adaptive grid allows high spatial resolution in regions of large density/velocity gradients and yet can be run on a workstation. The actual grid sizes vary from about 0.06 R(sub v) near the bowshock to 2 R(sub v) in the unperturbed solar wind. The results of the calculations are compared with observed magnetic field values obtained from the magnetometer on the Pioneer Venus Orbiter, at a time when the angle between the solar wind velocity vector and the interplanetary magnetic field (IMF) was only 7.6 deg. Good qualitative agreement between the observed and calculated field behavior is found. The overall results suggest that the induced magnetotail disappears when the IMF is radial for an extended time period and implies that it weakens when the field rotated through a near-radial orientation.

DeZeeuw, Darren L.↗
NASA NTRS1988

Magnetohydrodynamic (MHD) modelling of solar active phenomena via numerical methods

Numerical ideal MHD models for the study of solar active phenomena are summarized. Particular attention is given to the following physical phenomena: (1) local heating of a coronal loop in an isothermal and stratified atmosphere, and (2) the coronal dynamic responses due to magnetic field movement. The results suggest that local heating of a magnetic loop will lead to the enhancement of the density of the neighboring loops through MHD wave compression. It is noted that field lines can be pinched off and may form a self-contained magnetized plasma blob that may move outward into interplanetary space.

Wu, S. T.↗
DOE OSTI2022

A resistive MHD model and simulation on plasma flow evolution in the presence of resonant magnetic perturbation in a tokamak

Nonaxisymmetric magnetic fields such as the intrinsic error field and the externally applied resonant magnetic perturbation (RMP) in a tokamak are known to influence the plasma momentum transport and flow evolution through plasma response, which itself strongly depends on the plasma flow as well. The nonlinear interaction between plasma response and flow has been previously modeled in the conventional error field theory with the “no-slip” condition, which has been recently extended to allow the “free-slip” condition. In this work, we further target this specific process and numerically simulate the nonlinear plasma response and flow evolution in the presence of a single-helicity RMP in a circular-shaped model tokamak configuration, based on the full resistive MHD model in the initial-value code NIMROD. Time evolution of the parallel (to k) flow or “slip frequency” profile and its asymptotic steady state obtained from the NIMROD simulations are compared with both conventional and extended nonlinear response theories. Here, k is the wave vector of the propagating island. Good agreement with the extended theory with free-slip condition has been achieved for the parallel flow profile evolution in response to RMP in all resistive regimes, whereas the difference from the conventional theory with the no-slip condition tends to diminish as the plasma resistivity approaches zero.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

Three-Dimensional Multiscale MHD Model of Cometary Plasma Environments

First results of a three-dimensional multiscale MHD model of the interaction of an expanding cometary atmosphere with the magnetized solar wind are presented. The model starts with a supersonic and super-Alfvenic solar wind far upstream of the comet (25 Gm upstream of the nucleus) with arbitrary interplanetary magnetic field orientation. The solar wind is continuously mass loaded with cometary ions originating from a 10-km size nucleus. The effects of photoionization, electron impact ionization, recombination, and ion-neutral frictional drag are taken into account in the model. The governing equations are solved on an adaptively refined unstructured Cartesian grid using our new multiscale upwind scalar conservation laws-type numerical technique (MUSCL). We have named this the multiscale adaptive upwind scheme for MHD (MAUS-MHD). The combination of the adaptive refinement with the MUSCL-scheme allows the entire cometary atmosphere to be modeled, while still resolving both the shock and the diamagnetic cavity of the comet. The main findings are the following: (1) Mass loading decelerates the solar wind flow upstream of the weak cometary shock wave (M approximately equals 2, M(sub A) approximately equals 2), which forms at a subsolar standoff distance of about 0.35 Gm. (2) A cometary plasma cavity is formed at around 3 x 10(exp 3) km from the nucleus. Inside this cavity the plasma expands outward due to the frictional interaction between ions and neutrals. On the nightside this plasma cavity considerably narrows and a relatively fast and dense cometary plasma beam is ejected into the tail. (3) Inside the plasma cavity a teardrop-shaped inner shock is formed, which is terminated by a Mach disk on the nightside. Only the region inside the inner shock is the 'true' diamagnetic cavity. (4) The model predicts four distinct current systems in the inner coma: the density peak current, the cavity boundary current, the inner shock current, and finally the cross-tail current. (5) The calculated plasma parameters (magnetic field, plasma density, speed, and temperature) are in very good agreement with published Giotto observations.

Gombosi, Tamas I.↗
NASA NTRS2003

Three-Dimensional MHD Modeling of The Solar Corona and Solar Wind: Comparison with The Wang-Sheeley Model

We present simulation results from a tilted-dipole steady-state MHD model of the solar corona and solar wind and compare the output from our model with the Wang-Sheeley model which relates the divergence rate of magnetic flux tubes near the Sun (inferred from solar magnetograms) to the solar wind speed observed near Earth and at Ulysses. The boundary conditions in our model specified at the coronal base and our simulation region extends out to 10 AU. We assumed that a flux of Alfven waves with amplitude of 35 km per second emanates from the Sun and provides additional heating and acceleration for the coronal outflow in the open field regions. The waves are treated in the WKB approximation. The incorporation of wave acceleration allows us to reproduce the fast wind measurements obtained by Ulysses, while preserving reasonable agreement with plasma densities typically found at the coronal base. We find that our simulation results agree well with Wang and Sheeley's empirical model.

Usmanov, A. V.↗
NASA NTRS2023

Comparing Magnetopause Predictions From Two MHD Models During A Geomagnetic Storm and A Quiet Period

Magnetopause location is an important prediction of numerical simulations of the magnetosphere, yet the models can err, either under-predicting or over-predicting the motion of the boundary. This study compares results from two of the most widely used magnetohydrodynamic (MHD) models, the Lyon–Fedder–Mobarry (LFM) model and the Space Weather Modeling Framework (SWMF), to data from the GOES 13 and 15 satellites during the geomagnetic storm on 22 June 2015, and to THEMIS A, D, and E during a quiet period on 31 January 2013. The models not only reproduce the magnetopause crossings of the spacecraft during the storm, but they also predict spurious magnetopause motion after the crossings seen in the GOES data. We investigate the possible causes of the over-predictions during the storm and find the following. First, using different ionospheric conductance models does not significantly alter predictions of the magnetopause location. Second, coupling the Rice Convection Model (RCM) to the MHD codes improves the SWMF magnetopause predictions more than it does for the LFM predictions. Third, the SWMF produces a stronger ring current than LFM, both with and without the RCM and regardless of the LFM spatial resolution. During the non-storm event, LFM predicts the THEMIS magnetopause crossings due to the southward interplanetary magnetic field better than the SWMF. Additionally, increasing the LFM spatial grid resolution improves the THEMIS predictions, while increasing the SWMF grid resolutions does not.

magnetohydrodynamics↗
NASA NTRS1998

Io's Plasma Environment During the Galileo Flyby: Global Three-Dimensional MHD Modeling with Adaptive Mesh Refinement

The first results for applying a three-dimensional multimedia ideal MHD model for the mass-loaded flow of Jupiter's corotating magnetospheric plasma past Io are presented. The model is able to consider simultaneously physically realistic conditions for ion mass loading, ion-neutral drag, and intrinsic magnetic field in a full global calculation without imposing artificial dissipation. Io is modeled with an extended neutral atmosphere which loads the corotating plasma torus flow with mass, momentum, and energy. The governing equations are solved using adaptive mesh refinement on an unstructured Cartesian grid using an upwind scheme for AHMED. For the work described in this paper we explored a range of models without an intrinsic magnetic field for Io. We compare our results with particle and field measurements made during the December 7, 1995, flyby of to, as published by the Galileo Orbiter experiment teams. For two extreme cases of lower boundary conditions at Io, our model can quantitatively explain the variation of density along the spacecraft trajectory and can reproduce the general appearance of the variations of magnetic field and ion pressure and temperature. The net fresh ion mass-loading rates are in the range of approximately 300-650 kg/s, and equivalent charge exchange mass-loading rates are in the range approximately 540-1150 kg/s in the vicinity of Io.

Combi, M. R.↗
NASA NTRS1986

Three-dimensional, time-dependent, MHD model of a solar flare-generated interplanetary shock wave

A three-dimensional time-dependent MHD model of the propagation of an interplanetary shock wave into an ambient three-dimensional heliospheric solar wind is initialized with a peak velocity of 1000 km/s at the center of a right circular cone of 18 deg included angle at 18 solar radii. Differences from a previous 2-1/2 simulation (Wu et al., 1983; Gislason et al., 1984; Dryer et al., 1984) include diminuation of the solar peak velocity and concentration of the peak density at each radius. The IMF magnitude starts with high-latitude peaks, and helical-like IMF rotation is noted due to a large-amplitude nonlinear Alfven wave in the shocked plasma.

Dryer, M.↗