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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Fundamentals of Plasma Sails Propulsion Concept

The Mini-Magnetospheric Plasma Propulsion (M2P2), originally proposed by Winglee et al. [2000], is based on the two-fluid plasma model and requires a 15-km frontal standoff distance (or 20-km cross-sectional diameter) in order for the magnetic bubble to absorb sufficient momentum from the SW to accelerate a spacecraft to the unprecedented speeds of 50-80 km/s after an acceleration period of about three months. Winglee et al. [2000] derived the above size requirement based on an extrapolation of their simulated results in which a system much smaller than a M2P2 was used (p. 21,074 of their study). We submit, however, that a fluid model has no validity for such a small scale size-even in the region near the plasma source! It is assumed in the MHD fluid model, normally applied to the magnetosphere, that the characteristic scale-size is much greater than the Larmor radius and ion skin depth of the SW. In the case of the M2P2, however, the size of the magnetic bubble is actually less than or, at best, comparable to, the scale of these characteristic parameters and, therefore, a kinetic approach, which addresses the smallscale physical mechanisms involved, must be used. A fully three-dimensional version of the hybrid code is used in our M2P2 (Plasma Sails) studies was originally developed by Delamere et al. [1999]. The M2P2 plasma sail is an excellent application for this hybrid code. The primary advantage of this code is the seamless interface between fluid and kinetic descriptions of the ion populations. A kinetic description is not necessary for the dense inner regions of the magnetic bubble and tremendous computational savings can be realized by treating this dense, magnetized ion population with the fluid description. It is essential, however, that the outer bubble regions be treated kinetically as well as the SW protons. Comparison of full size M2P2 simulation based on 3D MHD and kinetic models show that kinetic treatment introduces much more asymmetry to the considering problem and the possibility of kinetic instabilities development.

Khazanov, G. V.↗
NASA NTRS2010

Interchange Reconnection and Coronal Hole Dynamics

We investigate the effect of magnetic reconnection between open and closed fields, often referred to as interchange reconnection, on the dynamics and topology of coronal hole boundaries. The most important and most prevalent three-dimensional topology of the interchange process is that of a small-scale bipolar magnetic field interacting with a large-scale background field. We determine the evolution of such a magnetic topology by numerical solution of the fully three-dimensional MHD equations in spherical coordinates. First, we calculate the evolution of a small-scale bipole that initially is completely inside an open field region and then is driven across a coronal hole boundary by photospheric motions. Next the reverse situation is calculated in which the bipole is initially inside the closed region and driven toward the coronal hole boundary. In both cases, we find that the stress imparted by the photospheric motions results in deformation of the separatrix surface between the closed field of the bipole and the background field, leading to rapid current sheet formation and to efficient reconnection. When the bipole is inside the open field region, the reconnection is of the interchange type in that it exchanges open and closed fields. We examine, in detail, the topology of the field as the bipole moves across the coronal hole boundary and find that the field remains well connected throughout this process. Our results, therefore, provide essential support for the quasi-steady models of the open field, because in these models the open and closed flux are assumed to remain topologically distinct as the photosphere evolves. Our results also support the uniqueness hypothesis for open field regions as postulated by Antiochos et al. On the other hand, the results argue against models in which open flux is assumed to diffusively penetrate deeply inside the closed field region under a helmet streamer. We discuss the implications of this work for coronal observations.

MHD (MAGNETOHYDRODYNAMICS)↗
NASA NTRS1985

Evolution and interaction of large interplanetary streams

A computer simulation for the evolution and interaction of large interplanetary streams based on multi-spacecraft observations and an unsteady, one-dimensional MHD model is presented. Two events, each observed by two or more spacecraft separated by a distance of the order of 10 AU, were studied. The first simulation is based on the plasma and magnetic field observations made by two radially-aligned spacecraft. The second simulation is based on an event observed first by Helios-1 in May 1980 near 0.6 AU and later by Voyager-1 in June 1980 at 8.1 AU. These examples show that the dynamical evolution of large-scale solar wind structures is dominated by the shock process, including the formation, collision, and merging of shocks. The interaction of shocks with stream structures also causes a drastic decrease in the amplitude of the solar wind speed variation with increasing heliocentric distance, and as a result of interactions there is a large variation of shock-strengths and shock-speeds. The simulation results shed light on the interpretation for the interaction and evolution of large interplanetary streams. Observations were made along a few limited trajectories, but simulation results can supplement these by providing the detailed evolution process for large-scale solar wind structures in the vast region not directly observed. The use of a quantitative nonlinear simulation model including shock merging process is crucial in the interpretation of data obtained in the outer heliosphere.

Whang, Y. C.↗
NASA NTRS1987

Macroscale particle simulation of kinetic Alfven waves

Two types of simulations of the kinetic Alfven wave are presented using a macroscale particle simulation code (Tanaka and Sato, 1986) which enables individual particle dynamics to be followed in the MHD scales. In this code, low frequency electromagnetic fields are solved by eliminating high frequency oscillations such as the light modes, and the scalar potential electric field is solved by eliminating Lagrangian oscillations. The dependences of the frequency and the Landau damping on the perpendicular wavenumber were studied, and good agreement was found between simulation and theoretical predictions. Some fundamental nonlinear interactions of the kinetic Alfven wave with the particles (parallel acceleration of the electrons) were also noted.

Tanaka, Motohiko↗
NASA NTRS1987

A magnetohydrodynamic model of the interaction of the solar wind with the Jovian magnetosphere and a magnetohydrodynamic simulation of the interaction of the solar wind with the out flowing plasma from a comet

A three-dimensional code for a rapidly rotating magnetosphere in which the MHD equations and the Maxwell equations were solved by using the two step Lax Endroff scheme, was developed. Preliminary results were presented at the Fall AGU meeting in San Francisco. The basic simulation model to study the solar wind interactions was adapted to other bodies in addition to Jupiter. Because of the recent comet flybys, a comet was chosen as the first model. The aim was to model the formation of the contact surface and the plasma tail. Later, work was begun on a three-dimensional model which would include the effects of mass loading. This model was designed to study the weak cometary bow shocks observed by the probes to comets Halley and Giacobini-Zinner. The model was successful in reproducing the position and shape of the bow shock which was determined by using observations from the Suisei spacecraft.

Walker, R. J.↗
NASA NTRS1987

Effect of a localized minimum in equatorial field strength on resistive tearing instability in the geomagnetotail

A two-dimensional, resistive-MHD computer code is used to investigate the spontaneous reconnection of magnetotaillike configurations. The initial conditions adopted in the simulations are of two types: (1) in which the equatorial normal magnetic field component B(ze) declines monotonically down the tail, and (2) in which B(ze) exhibits a deep minimum in the near-earth plasma sheet. To represent the case where the earthward convection stops before the X line forms, zero-flow boundary conditions are imposed at the edges of the computational box. The initial configurations are in equilibrium and table within ideal MHD. The dynamic evolution of the system starts after the resistivity is turned on. The main results of these simulations basically support the neutral-line model of substorms and confirm Birn's (1980) computer studies. Spontaneous formation of an X-type neutral point and a single O-type plasmoid with strong tailward flow on the tailward side of the X point is found. The time interval from the turning on of the resistivity to the formation of a plasmoid is much shorter in the case where there is an initial deep minimum. A simple analytic calculation is also carried out to demonstrate why the configuration with a deep minimum is more susceptible to the development of the neutral point.

Hau, L.-N.↗
NASA NTRS1987

Numerical simulation of the generation of turbulence from cometary ion pick-up

Observations of magnetic field fluctuations near Comet Halley have revealed a rapid development of a Kolmogoroff-like turbulence spectrum extending from below 0.01 Hz to above 0.1 Hz. Spectra obtained far from the comet have a strong peak in power near the Doppler-shifted ion-cyclotron frequency of singly ionized water. Closer to the comet, the spectrum at higher frequencies is enhanced in power level over the background solar wind spectrum by approximately an order of magnitude. The equations of incompressible MHD are solved using a two-dimensional 256 x 256 mode spectral method code to simulate this spectral evolution as an inertial range turbulent cascade. The initial conditions contained a constant magnetic field and a single coherent wave mode at a low wave number. The solar wind turbulence was modeled by a background noise spectrum having a Kolmogoroff spectral index. The coherent mode decayed into an inertial range spectrum with Kolmogoroff slope within a few eddy-turnover times. Both the time scale and the increase in power level of the turbulence seen in the simulation are in accord with the Giotto observations.

Goldstein, M. L.↗
NASA NTRS2010

Impact of Magnetic Draping, Convection, and Field Line Tying on Magnetopause Reconnection Under Northward IMF

We simulate a northward IMF cusp reconnection event at the magnetopause using the OpenGGCM resistive MHD code. The ACE input data, solar wind parameters, and dipole tilt belong to a 2002 reconnection event observed by IMAGE and Cluster. Based on a fully three-dimensional skeleton separators, nulls, and parallel electric fields, we show magnetic draping, convection, ionospheric field line tying play a role in producing a series of locally reconnecting nulls with flux ropes. The flux ropes in the cusp along the global separator line of symmetry. In 2D projection, the flux ropes the appearance of a tearing mode with a series of 'x's' and 'o's' but bearing a kind of 'guide field' that exists only within the magnetopause. The reconnecting field lines in the string of ropes involve IMF and both open and closed Earth magnetic field lines. The observed magnetic geometry reproduces the findings of a superposed epoch impact parameter study derived from the Cluster magnetometer data for the same event. The observed geometry has repercussions for spacecraft observations of cusp reconnection and for the imposed boundary conditions reconnection simulations.

Wendel, Deirdre E.↗
NASA NTRS1989

Development and simulation study of a new inverse-pinch high Coulomb transfer switch

The inverse-pinch plasma switch was studied using a computer simulation code. The code was based on a 2-D, 2-temperature magnetohydrodynamic (MHD) model. The application of this code was limited to the disk-type inverse-pinch plasma switch. The results of the computer analysis appear to be in agreement with the experimental results when the same parameters are used. An inverse-pinch plasma switch for closing has been designed and tested for high-power switching requirements. An azimuthally uniform initiation of breakdown is a key factor in achieving an inverse-pinch current path in the switch. Thus, various types of triggers, such as trigger pins, wire-brush, ring trigger, and hypocycloidal-pinch (HCP) devices have been tested for uniform breakdown. Recently, triggering was achieved by injection of a plasma-ring (plasma puff) that is produced separately with hypocycloidal-pinch electrodes placed under the cathode of the main gap. The current paths at switch closing, initiated by the injection of a plasma-ring from the HCP trigger are azimuthally uniform, and the local current density is significantly reduced, so that damage to the electrodes and the insulator surfaces is minimized. The test results indicate that electron bombardment on the electrodes and the insulator surfaces is minimized. The test results indicate that electron bombardment on the electrodes is four orders of magnitude less than that of a spark-gap switch for the same switching power. Indeed, a few thousand shots with peak current exceeding a mega-ampere and with hold-off voltage up to 20 kV have been conducted without showing measurable damage to the electrodes and insulators.

Choi, Sang H.↗
NASA NTRS1993

Magnetic effects in Venus/solar wind interaction

The overall objective of this research program is to better understand the interaction of a magnetized solar wind with the Venus atmosphere through the use of numerical solutions of the time-dependent, 2-D and 3-D magnetohydrodynamic (MHD) equations. Due to the more modest CPU requirements for the 2-D simulations, they were used for studies in which useful information not dependent on the third dimension could be obtained. The 2-D simulations served several purposes in addition to providing useful physical insight. They were used to determine the numerical parameters required in the 3-D studies, such as the grid spacing required to resolve particular features, and the damping that must be included to remove high-frequency oscillations. Among the specific studies performed with support from this grant that are discussed in this report are the following: comparison with other available models for purposes of testing the code and obtaining a baseline with which to evaluate the effects of additional physical processes, effects of a finite planet conductivity, bow shock standoff distance, and formation of the magnetic barrier and slippage of the magnetic field around the planet. A brief description of the methodology is presented before discussing the results.

Steinolfson, Richard S.↗
NASA NTRS2019

Initial Results from the GEM Challenge on the Spacecraft Surface Charging Environment

Spacecraft surface charging during geomagnetically disturbed times is one of the most important causes of satellite anomalies. Predicting the surface charging environment is one prevalent task of the geospace environment models. Therefore, the Geospace Environment Modeling (GEM) Focus Group "Inner Magnetosphere Crossenergy/Population Interactions" initiated a communitywide challenge study to assess the capability of several inner magnetosphere ring current models in determining surface charging environment for the Van Allen Probes orbits during the 17 March 2013 storm event. The integrated electron flux between 10 and 50 keV is used as the metrics. Various skill scores are applied to quantitatively measure the modeling performance against observations. Results indicate that no model consistently perform the best in all of the skill scores or for both satellites. We find that from these simulations the ring current model with observational flux boundary condition and Weimer electric potential driver generally reproduces the most realistic flux level around the spacecraft. A simple and weaker VollandStern electric field is not capable of effectively transporting the same plasma at the boundary toward the Earth. On the other hand, if the ring current model solves the electric field selfconsistently and obtains similar strength and pattern in the equatorial plane as the Weimer model, the boundary condition plays another crucial role in determining the electron flux level in the inner region. When the boundary flux spectra based on magnetohydrodynamics (MHD) model/empirical model deviate from the shape or magnitude of the observed distribution function, the simulation produces poor skill scores along Van Allen Probes orbits.

Yu, Yiqun↗
NASA NTRS1992

Numerical simulation of geometric scale effects in cylindrical self-field MPD thrusters

A 2D, two-temperature, single fluid MHD code which incorporates classical plasma transport coefficients and Hall effects has been developed to predict steady-state, self-field MPD thruster performance. The governing equations and numerical methods of solution are outlined and discussed. Experimental comparisons are used to validate model predictions. The model accurately predicts thrust and reproduces trends in the discharge voltage for discharge currents below experimentally measured onset values. However, because the model does not include electrode effects the calculated voltage drops are significantly lower than experimentally measured values. Predictions of thrust and flow efficiency are made for a matrix of fifteen cylindrical thruster geometries assuming a fully ionized argon propellant.

Lapointe, M.↗
NASA NTRS1995

A global 3-D MHD model of the solar wind with Alfven waves

A fully three-dimensional solar wind model that incorporates momentum and heat addition from Alfven waves is developed. The proposed model upgrades the previous one by considering self-consistently the total system consisting of Alfven waves propagating outward from the Sun and the mean polytropic solar wind flow. The simulation region extends from the coronal base (1 R(sub s) out to beyond 1 AU. The fully 3-D MHD equations written in spherical coordinates are solved in the frame of reference corotating with the Sun. At the inner boundary, the photospheric magnetic field observations are taken as boundary condition and wave energy influx is prescribed to be proportional to the magnetic field strength. The results of the model application for several time intervals are presented.

Usmanov, A. V.↗
NASA NTRS1972

Current Sheet Magnetic Model for the Solar Corona

A new magnetic model is developed and compared with previous models and the observed solar corona. An attempt is made to more accurately compute the three-dimensional currents flowing in the solar corona. Physical reasons are given that require most of the large scale currents flowing in the solar corona to lie near thin sheets. The current sheets are not constrained into any particular geometry or symmetry as in the previous models of Altschuler and Newkirk (1969) and Schatten et al. (1969). A comparison with the axisymmetric, isothermal MHD solution of Pneuman and Kopp (1970) suggests that the model is able to simulate to high accuracy an isothermal corona. A comparison of the model with the May 30, 1965, solar eclipse and the November 12, 1966, solar eclipse shows the model is capable of computing many features including the polar plume orientations as well as radial and nonradial streamers in the solar corona.

Schatten, K. H.↗
NASA NTRS1985

Forced MHD turbulence in a uniform external magnetic field

Two-dimensional dissipative MHD turbulence is randomly driven at small spatial scales and is studied by numerical simulation in the presence of a strong uniform external magnetic field. A behavior is observed which is apparently distinct from the inverse cascade which prevails in the absence of an external magnetic field. The magnetic spectrum becomes dominated by the three longest wavelength Alfven waves in the system allowed by the boundary conditions: those which, in a box size of edge 2 pi, have wave numbers (kx' ky) = (1, 1), and (1, -1), where the external magnetic field is in the x direction. At any given instant, one of these three modes dominates the vector potential spectrum, but they do not constitute a resonantly coupled triad. Rather, they are apparently coupled by the smaller-scale turbulence.

Hossain, M.↗
NASA NTRS1985

Forced magnetohydrodynamic turbulence in a uniform external magnetic field

Two-dimensional dissipative MHD turbulence is randomly driven at small spatial scales and is studied by numerical simulation in the presence of a strong uniform external magnetic field. A behavior is observed which is apparently distinct from the inverse cascade which prevails in the absence of an external magnetic field. The magnetic spectrum becomes dominated by the three longest wavelength Alfven waves in the system allowed by the boundary conditions: those which, in a box size of edge 2 pi, have wave numbers (kx, ky) = (1, 1), and (1, -1), where the external magnetic field is in the x direction. At any given instant, one of these three modes dominates the vector potential spectrum, but they do not constitute a resonantly coupled triad. Rather, they are apparently coupled by the smaller-scale turbulence.

Hossain, M.↗
NASA NTRS1988

Numerical simulation of mass injection for the formation of prominence magnetic field configurations. I - Asymmetric injection

A two-dimensional ideal MHD numerical model is used to investigate chromospheric mass injection into an overlying coronal dipole magnetic field. Such injection is shown to produce magnetic field deformations conducive to the formation of active region prominences. The results support a model in which an absorptive strand is formed by chromospheric asymmetric mass injection into the overlying coronal magnetic field. A necessary condition for the accumulation of the strands is that the mass injection forms a Kippenhahn-Schluter-type (1957) field configuration.

An, Chang-Hyuk↗
NASA NTRS1991

Mass-loading and parallel magnetized shocks

Recent observations at comets Giacobini-Zinner and Halley suggest that simple nonreacting gas dynamics or MHD is an inappropriate description for the bow shock. The thickness of the observed (sub)shock implies that mass-loading is an important dynamical process within the shock itself, thereby requiring that the Rankine-Hugoniot conditions possess source terms. This leads to shocks with properties similar to those of combustion shocks. The paper considers parallel magnetized shocks subjected to mass-loading, describes some properties which distinguish them from classical MHD parallel shocks, and establishes the existence of a new kind of MHD compound shock. These results will be of importance both to observations and numerical simulations of the comet-solar wind interaction.

Zank, G. P.↗