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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

A new computational model for the prediction of mass loading phenomena for solar wind interactions with cometary and planetary ionospheres

The modified gasdynamic convected-magnetic-field MHD model developed by Spreiter and Stahara (1980) to simulate the supersonic flow of the solar wind past planetary magnetoionospheres is extended to account for cases (such as Venus and comets) in which significant numbers of neutral atmospheric atoms become ionized in the surrounding flow and add to its momentum, energy, and mass. The mathematical model and the solution procedures for the nose and tail regions are explained; typical computational grids are shown; and numerical results for a comet and for the Venus ionosheath and bow shock are presented graphically. It is found that the bow shock weakens and moves further upstream of the obstacle as mass loading is increased and the flow upstream of the bow shock becomes more compressed.

Stahara, Stephen S.↗
NASA NTRS1995

Large-scale structure of the solar corona and inner heliosphere

The large-scale structure of the solar corona influences solar activity particularly coronal mass ejections (CMEs). The helmet streamers that are observed to dominate the structure of the inner corona are formed by the interaction of the solar wind with coronal magnetic fields. We have simulated this interaction in three dimensions using the magnetohydrodynamic (MHD) equations. In order to create a realistic model, we use the magnetic field that is observed at the Sun's surface (deduced from daily Wilcox Solar Observatory magnetograms) as input, in combination with specified density and temperature profiles at the surface. A self-consistent 3D solar-wind solution is developed by integrating the MHD equations in time to steady state. Such solutions can reproduce the observed structures that are seen in coronagraph images and eclipse photographs of the corona. This model allows us to accurately determine the position of the heliospheric current sheet. We will compare the results obtained from our model with Ulysses observations during the period May-June 1993, and with an eclipse photograph of the corona on November 3, 1994.

Mikic, Z.↗
NASA NTRS2002

A Theoretical Study of the Global Geospace System Simulations

This report describes work undertaken by the principal investigator (PI) and her colleagues in three areas of research: Relativistic Electron Dynamics, Solar Energetic Particle Study, and ULF Wave Statistical Study. The researchers conducted three dimensional modeling of relativistic electron dynamics, simulated Solar Energetic Particles (SEP) trapping using a magnetohydrodynamic (MHD) code, and a study of Pc 5 magnetic pulsations.

Hudson, Mary K.↗
NASA NTRS2010

Dynamics of Ring Current and Electric Fields in the Inner Magnetosphere During Disturbed Periods: CRCM-BATS-R-US Coupled Model

We present simulation results from a one-way coupled global MHD model (Block-Adaptive-Tree Solar-Wind Roe-Type Upwind Scheme, BATS-R-US) and kinetic ring current models (Comprehensive Ring Current Model, CRCM, and Fok Ring Current, FokRC). The BATS-R-US provides the CRCM/FokRC with magnetic field information and plasma density/temperature at the polar CRCM/FokRC boundary. The CRCM uses an electric potential from the BATS-R-US ionospheric solver at the polar CRCM boundary in order to calculate the electric field pattern consistent with the CRCM pressure distribution. The FokRC electric field potential is taken from BATS-R-US ionospheric solver everywhere in the modeled region, and the effect of Region II currents is neglected. We show that for an idealized case with southward-northward-southward Bz IMF turning, CRCM-BATS-R-US reproduces well known features of inner magnetosphere electrodynamics: strong/weak convection under the southward/northward Bz; electric field shielding/overshielding/penetration effects; an injection during the substorm development; Subauroral Ion Drift or Polarization Jet (SAID/PJ) signature in the dusk sector. Furthermore, we find for the idealized case that SAID/PJ forms during the substorm growth phase, and that substorm injection has its own structure of field-aligned currents which resembles a substorm current wedge. For an actual event (12 August 2000 storm), we calculate ENA emissions and compare with Imager for Magnetopause-to-Aurora Global Exploration/High Energy Neutral Atom data. The CRCM-BATS-R-US reproduces both the global morphology of ring current and the fine structure of ring current injection. The FokRC-BATS-R-US shows the effect of a realistic description of Region II currents in ring current-MHD coupled models.

Buzulukova, N.↗
NASA NTRS2018

Simulated Encounters of the Parker Solar Probe with a Coronal-Hole Jet

Solar coronal jets are small, transient, collimated ejections most easily observed in coronal holes (CHs). The upcoming Parker Solar Probe (PSP) mission provides the first opportunity to encounter CH jets in situ near the Sun and examine their internal structure and dynamics. Using projected mission orbital parameters, we have simulated PSP encounters with a fully three-dimensional magnetohydrodynamic (MHD) model of a CH jet. We find that three internal jet regions, featuring different wave modes and levels of compressibility, have distinct identifying signatures detectable by PSP. The leading Alfvén wave front and its immediate wake are characterized by transAlfvénic plasma flows with mild density enhancements. This front exhibits characteristics of a fast switch-on MHD shock, whose arrival is signaled by the sudden onset of large-amplitude transverse velocity and magnetic-field oscillations highly correlated in space and time. The trailing portion is characterized by supersonic but subAlfvénicout flows of dense plasma with uncorrelated velocity and magnetic-field oscillations. This compressible region contains most of the jet's mass. The volume between the immediate wake and dense jet, the remote wake,mixes and transitions the characteristics of the two other regions. In addition to probing each region separately, we also simulate a corotational PSP-jet encounter. In this scenario, the simulated spacecraft hovers over the jet producing CH, as may occur during the mission's corotational phases, sampling each jet region in turn. We estimate that PSP will encounter numerous CH jets over the lifetime of the mission.

Roberts, Merrill A.↗
NASA NTRS2001

MHD Instability and Turbulence in the Tachocline

In this quarter we have begun simulations on the Cray T3E at PSC and we are debugging our code on the TSC. The PSC simulations are examining stratified shear turbulence with a flow-aligned magnetic field and passive tracer particles. We have conducted analysis of neutral simulations to establish a firm basis of comparison. Second-order structure functions have been computed, fit, and compared to theoretical expressions relating the dissipation fields and the structure-function-fit parameters. Agreement with high-Reynolds number observations is excellent, giving us confidence that the lower-Re simulations are relevant to higher-Re flows. We have also evaluated the neutral layer anisotropy.

Werne, Joseph↗
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 NTRS1992

Preliminary study of the CRRES magnetospheric barium releases

Preliminary theoretical and computational analyses of the Combined Release and Radiation Effects Satellite (CRRES) magnetospheric barium releases are presented. The focus of the studies is on the evolution of the diamagnetic cavity which is formed by the barium ions as they expand outward, and on the structuring of the density and magnetic field during the expansion phase of the releases. Two sets of simulation studies are discussed. The first set is based upon a 2D ideal MHD code and provides estimates of the time and length scales associated with the formation and collapse of the diamagnetic cavity. The second set uses a nonideal MHD code; specifically, the Hall term is included. This additional term is critical to the dynamics of sub-Alfvenic plasma expansions, such as the CRRES barium releases, because it leads to instability of the expanding plasma. Detailed simulations of the G4 and G10 releases were performed. In both cases the expanding plasma rapidly structured: the G4 release structured at time t less than about 3 s and developed scale sizes of about 1-2 km, while the G10 release structured at time t less than about 22 s and developed scale sizes of about 10-15 km. It is also found that the diamagnetic cavity size is reduced from those obtained from the ideal MHD results because of the structure. On the other hand, the structuring allows the formation of plasma blobs which appear to free stream across the magnetic field; thus, the barium plasma can propagate to larger distances traverse to the magnetic field than the case where no structuring occurs. Finally, a new normal mode of the system was discovered which may be excited at the leading edge of the expanding barium plasma.

Huba, J. D.↗
NASA NTRS1978

Magnetohydrodynamic simulation of coronal mass ejections into the solar wind

A state of the art fluid continuum technique to describe the MHD transient respose of the corona below 10 solar radii during two well observed events (flares on 21 August 1973 and 5 September 1973) is reviewed. It is concluded that the computer simulation should be subjected to in situ verification of as many of its initial assumptions as possible. Its ability to provide a rational basis for physical understanding of mass ejections suggests its use as one of the tools used in the planning and analysis of such encounter missions.

Dryer, M.↗
NASA NTRS1986

A numerical model of magnetosphere-ionosphere coupling Preliminary results

A three-dimensional simulation model was developed to study magnetosphere-ionosphere coupling in the auroral region. One-fluid MHD equations are adopted to model the magnetosphere, and current density continuity equations were solved consistently to model the ionosphere. In the preliminary simulation runs described here, the electrodynamics of region 1 field-aligned currents were modeled. Initially, the electric field is taken to be in the magnetospheric equatorial plane. Alfven waves then propagate down to the ionosphere, accompanying a field-aligned current and exciting the electrostatic potential in the ionosphere by electrodynamic coupling. The results of the preliminary runs directly correspond with the fundamental characteristics of global magnetosphere-ionosphere coupling. These characteristics include an ionospheric electrostatic potential which varies in its development in time, depending upon the ratio of ionospheric resistance to magnetospheric impedance, the flowing of field-aligned currents into the ionosphere on the dawnside and out on the duskside, and the distribution along the geomagnetic field of the amplitude of the field-aligned current density in proportion to the field intensity.

Watanabe, K.↗
NASA NTRS1983

Reconnection rates, small scale structures and simulations

The study of reconnection in the context of one fluid, two dimensional magnetohydrodynamics (MHD), with spatially uniform constant density, viscosity and resistivity is though to retain most of the physics important in reconnection. Much of the existing reconnection literature makes use of this approach. This discussion focuses on attempts to determine the properties of reconnection solutions to MHD as precisely as possible without regard to the intrinsic limitations of the model.

Matthaeus, W. H.↗
NASA NTRS1991

A numerical simulation of magnetic reconnection and radiative cooling in line-tied current sheets

Radiative MHD equations are used for an optically thin plasma to carry out a numerical experiment related to the formation of 'postflare' loops. The numerical experiment starts with a current sheet that is in mechanical and thermal equilibrium but is unstable to both tearing-mode and thermal-condensation instabilities. The current sheet is line-tied at one end to a photospheric-like boundary and evolves asymmetrically. The effects of thermal conduction, resistivity variation, and gravity are ignored. In general, reconnection in the nonlinear stage of the tearing-mode instability can strongly affect the onset of condensations unless the radiative-cooling time scale is much smaller than the tearing-mode time scale. When the ambient plasma is less than 0.2, the reconnection enters a regime where the outflow from the reconnection region is supermagnetosonic with respect to the fast-mode wave speed. In the supermagnetosonic regime the most rapidly condensing regions occur downstream of a fast-mode shock that forms where the outflow impinges on closed loops attached to the photospheric-like boundary. A similar shock-induced condensation might occur during the formation of 'postflare' loops.

Forbes, T. G.↗
NASA NTRS1997

Gamma-Ray Bursts and Relativistic MHD around Coalescing Neutron Stars

In this proposal we performed the first relativistic hydrodynamic simulations of the binary neutron stars near coalescence. Calculations were performed in three spatial dimensions and were supplemented with magnetohydrodynamic simulations of the magnetic reconnection as the neutron star fluid responds to the relativistic forces of the binary. The ultimate goal of the proposed work was to develop a theoretical model for the temporal and spectral evolution of cosmological gamma-ray bursts produced during the magnetic reconnection. This is the first time such hydrodynamic calculations were performed and some surprising results were obtained. The neutrons stars were observed to compress and heat in the simulations. They could even collapse to black holes many orbits before coalescence. The question then arose as to whether this strange phenomenon could power a gamma-ray burst. The bulk of computing effort in this proposal was therefore devoted to an analysis of the implications of this possible new gamma-ray burst paradigm. The initial focus of the work performed under this proposal was therefore to utilize (3+1) dimensional and spherical numerical general relativistic hydrodynamics to study the origin, evolution, and parametric sensitivity of the collapse instability. We derived modified conditions of hydrostatic equilibrium for stars in the curved space of quasi-static orbits. We followed the magnetic field evolution for approximately one orbit (10 msec). The magnetic field energy was observed to exponentiate with an e-folding time of about 1 msec. Thus, the field should build up very quickly to a magnitude such that reconnection and back reaction of the fluid inhibits further growth.

Mathews, Grant J.↗
NASA NTRS1979

Numerical study of an explosion in a non-homogeneous medium with and without magnetic fields

A version of the two step Lax-Wendroff difference method with second order accuracy is used to seek solutions of the unsteady magnetohydrodynamic (MHD) equations for the study of an explosion in a non-homogeneous medium with and without magnetic fields. The explosion is initiated by introducing a finite amount of energy in a small volume of gas which leads to an instantaneous increase of gas temperature. Two types of magnetic field configurations (i.e. open and closed) are considered to illustrate the dependence and differences of magnetic and a non-magnetic flow motion. Numerical results show the strong dependence of induced MHD flow field on the magnetic field configuration and strength. The development and decay of the fast and slow MHD shock waves, as well as the ordinary gasdynamic shock waves, are also represented well by the present computer simulations. In conclusion, we have demonstrated that the numerical scheme we have used is a reliable one for studying 2-dimensional MHD problems in non-homogeneous medium.

Han, S. M.↗
NASA NTRS1992

Hybrid simulations of intermediate shocks - Coplanar and noncoplanar solutions

A hybrid code is used to investigate the kinetic structure and stability of subfast intermediate shocks (IS) formed dynamically by the interaction between a flowing plasma and a stationary piston. Results of the kinetic simulation of noncoplanar ISs are compared with predictions of the MHD theory. The relevance of the results of the study to observations of the magnetopause is discussed.

Karimabadi, H.↗
NASA NTRS1992

A powerful local shear instability in weakly magnetized disks. III - Long-term evolution in a shearing sheet. IV - Nonaxisymmetric perturbations

The nonlinear evolution of the recently identified accretion disk magnetic shear instability is investigated through a series of numerical simulations. Finite-difference computations of the equations of compressible MHD are carried out on an axisymmetric shearing sheet system with periodic boundary conditions designed to approximate a local region within an accretion disk. Initial field configurations that include some net vertical component evolve into a nonlinear, exponentially growing solution with large poloidal velocities and magnetic fields with energies comparable to the thermal energy density. The stability of a purely azimuthal field configuration is examined, and it is found that nonaxisymmetric instability is present, but with a growth time measured in tens of orbital periods. In general, the most rapid growth occurs for very small radial and azimuthal wavenumbers, leading to coherent magnetic field structure in planes parallel to the disk. It is suggested that this instability is a key ingredient for the generation of magnetic fields in disks.

Hawley, John F.↗
NASA NTRS2003

A Plasmoid Thruster for Space Propulsion

There are a number of possible advantages to using accelerated plasmoids for in-space propulsion. A plasmoid is a compact plasma structure with an integral magnetic field. They have been studied extensively in controlled fusion research and are classified according to the relative strength of the poloidal and toroidal magnetic field (BP and Bt, respectively). An Object with B P t >> 1 is classified as a Field Reverse Configuration (FRC); if B, = Bt, it is called a Spheromak. The plasmoid thruster operates by producing FRC-like plasmoids, and subsequently ejecting them from the device at high velocity. The plasmoid is formed inside of a single turn conical theta-pinch coil. As this process is inductive, there are no electrodes. Similar experiments have yielded plasmoid velocities of at least 50 km/s (l), and calculations indicate that velocities in excess of 100 km/s should be possible. This concept should be capable of producing Isp s in the range of 5,000 - 10,000 s with thrust densities of order 10(exp 5) N/sq m. The current experiment is designed to produce jet powers in the range of 5-10 kW, although the concept should be scalable to several MW s. The plasmoids mass and velocity will be measured with a variety of diagnostics, including internal and external B-dot probes, flux loops, Langmuir probes, high-speed cameras, and a laser interferometer. Also of key importance will be measurements of the efficiency and mass utilization. Simulations of the plasmoid thruster using MOQUI, a time dependent MHD code, will be carried out concurrently with experimental testing.

Koelfgen, Syri J.↗
NASA NTRS2003

The Plasmoid Thruster Experiment (PTX)

A plasmoid is a compact plasma structure with an integral magnetic field. They have been studied extensively in controlled fusion research and are categorized according to the relative strength of the poloidal and toroidal magnetic field (B(phi), and B(tau), respectively). An object with B(phi)/B(tau) >> 1 is classified as a Field Reverse Configuration (FRC); if B(phi) = B(tau), it is called a Spheromak. There are a number of possible advantages to using accelerated plasmoids for in-space propulsion. A thruster based on this concept would operate by repetitively producing plasmoids and ejecting them from the device at high velocity. The plasmoid is formed inside of a single turn conical theta-pinch coil; as this process is inductive, there are no life-limiting electrodes. Similar experiments have yielded plasmoid velocities of at least 50 km/s (l), and calculations indicate that velocities in excess of 100 km/s are possible. A thruster based on this concept would be capable of producing an I(sp) in the range of 5,000 - 10,OOO s, with thrust densities of order 10(exp 5) N/m(exp 2). The current experiment is designed to produce jet powers in the range of 5-10 kW, although the concept should be scalable to higher power. The purpose of this experiment is to determine the feasibility of this plasma propulsion concept. To accomplish this, it will be necessary to determine: a.) specific impulse and thrust, b.) efficiency and mass utilization, c.) which type of plasmoid (FRC-like or Spheromak-like) gives the best performance, and d.) the characteristics required of actual thruster components (i.e., switch and capacitor technology). The plasmoid mass and velocity will be measured with a variety of diagnostics, including internal and external B-dot probes, flux loops, Langmuir probes, high-speed cameras, and an interferometer. Simulations of the plasmoid thruster using MOQUI, a time dependent MHD code, will be carried out concurrently with experimental testing. The PTX device is currently undergoing initial testing and preliminary experimental results are presented.

Eskridge, Richard↗