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 685 records · Page 38

A two-dimensional MHD global coronal model - Steady-state streamers

We describe a two-dimensional time-dependent, numerical, magnetohydrodynamic model for the determination of the physical properties of coronal streamers from the top of the transition zone (solar radius = 1) to 15 solar radii. Four examples are given: for dipole, quadrupole, and hexapole initial field topologies. The computed parameters are density, temperature, velocity, and magnetic field. In addition to the properties of the solutions, their accuracy is discussed. We use the model as the basis for a general discussion of the way boundary conditions are specified in this and similar simulations.

Wang, A.-H.↗

Multiscale Hybrid Modeling of the Impact Response of the Earth's Magnetotail to Ionospheric O+ Outflow

Ionospheric outflow plays an important role in coupling the ionosphere with the solar wind-magnetosphere system. Previous multi-fluid MHD studies explored the global influence of oxygen ions of ionospheric origin (O+) on magnetospheric dynamics. A detailed exploration of the interaction of ionospheric ions with the magnetotail requires kinetic treatment for ions. We perform a self-consistent investigation of these processes with a three-dimensional space-time adaptive hybrid code, HYPERS, powered by an intelligent Event-driven Multi-Agent Planning System (EMAPS). By comparing simulations with and without outflow we conclude that oxygen ions, flowing from the ionosphere through the lobes into the tail or directly entering the inner magnetosphere, are able to significantly modify the magnetotail configuration and induce X-points and current sheet structures at locations where magnetic reconnection does not occur in a simulation without outflow, potentially very close to the Earth. This finding may have implications for interpreting substorms and magnetotail reconnection events observed for southward magnetic field simultaneously with significant contents of oxygen ions of ionospheric origin.

multiscale↗

The Influence of Pickup Protons, from Interstellar Neutral Hydrogen, on the Propagation of Interplanetary Shocks from the Halloween 2003 Solar Events to ACE and Ulysses: A 3-D MHD Modeling Study

We describe our 3-D, time ]dependent, MHD solar wind model that we recently modified to include the physics of pickup protons from interstellar neutral hydrogen. The model has a time-dependent lower boundary condition, at 0.1 AU, that is driven by source surface map files through an empirical interface module. We describe the empirical interface and its parameter tuning to maximize model agreement with background (quiet) solar wind observations at ACE. We then give results of a simulation study of the famous Halloween 2003 series of solar events. We began with shock inputs from the Fearless Forecast real ]time shock arrival prediction study, and then we iteratively adjusted input shock speeds to obtain agreement between observed and simulated shock arrival times at ACE. We then extended the model grid to 5.5 AU and compared those simulation results with Ulysses observations at 5.2 AU. Next we undertook the more difficult tuning of shock speeds and locations to get matching shock arrival times at both ACE and Ulysses. Then we ran this last case again with neutral hydrogen density set to zero, to identify the effect of pickup ions. We show that the speed of interplanetary shocks propagating from the Sun to Ulysses is reduced by the effects of pickup protons. We plan to make further improvements to the model as we continue our benchmarking process to 10 AU, comparing our results with Cassini observations, and eventually on to 100 AU, comparing our results with Voyager 1 and 2 observations.

Detman, T. R.↗

A predicted energetic ion signature of a neutral line in the geomagnetic tail

Under certain assumptions one can predict the formation of a magnetic neutral line on theoretical grounds, generally as a result of an unstable current sheet. Neutral points also seen to form spontaneously in a wide variety of plasma simulatioons, from MHD to fully kinetic. Hence, such magnetic null regions are central to some models of magnetotail energization and magnetic substorms. Yet the observational evidence for such regions remains controversial. This paper presents a new signature of a neutral line: the field-aligned ion beam characteristic of current sheet acceleration becomes dispersed into a 'ridge' in the velocity space distribution function as a neutral line is approached. The results were obtained using single-particle simulations in a fixed magnetic field model. The origin of the ridge is shown to be due to initially earthward going ions which pass through the neutral line region. These ions come from a smaller part of the initially antiearthward flowing distribution, generating the depletion of the distribution function below the ridge.

Martin, R. F., Jr.↗

Hydrodynamic Stability and Magnetic Reconnection in Disks and Stars

The purpose of this grant is to study parametric instability. The simplest example of parametric instability is a harmonic oscillator with a periodic modulation of the spring constant. If the modulation frequency is close to twice the natural frequency of the oscillator, the amplitude of oscillation tends to grow exponentially. The growth rate is proportional to the strength of the modulation, but it also depends upon the closeness to resonance of the two frequencies, and upon natural damping rate or "Q" of the oscillator. Parametric instabilities are very common in physics. A familiar example is a jogger's ponytail--normally a very strongly damped pendulum, it can be destabilized by the variation in effective gravity during the jogger's stride. Observation confirms that the period of the pendulum is half that of the jogger's vertical motion. In astrophysics, parametric instability may occur by external tidal forcing, or by interaction among eigenmodes. In the latter case, an energetic eigenmode may destabilize modes of half its frequency, provided some weak nonlinearity exists to couple them. Under a previous Astrophysical Theory grant (NAGW-2419), the PI discovered a parametric instability of tidally forced disks such as the accretion disks in cataclysmic variables and X ray binaries [2]. The destabilized modes are tightly-wound, incompressible, three-dimensional waves analogous to g-modes and r-modes in stars. Later work has confirmed our analysis [4]. It was hoped that these modes might provide a source of turbulence and angular momentum transport in accretion disks. However, a follow-up investigation of this instability by local numerical simulations, although confirming the analytically estimated growth rates, found negligible angular momentum flux [3]. Other work, partly supported by the ATP, now strongly indicates that the transport mechanism in such disks is magnetohydrodynamic turbulence [6]. Nevertheless, the parametric mechanism may truncate the outer edges of disks in close binaries [2], and it may be important in disks of very low ionization such as protostellar disks, or even cataclysmic-variable disks in quiescence where the MHD mechanism may be ineffective [5]. All analyses up to 1996 were done in a local approximation where the orbital frequency, shear rate, and tidal field were treated as constants. The locally computed growth rate turns out to depend strongly on radius, and it was unclear how to average these local rates to obtain the correct global rate. This is a critical issue for accretion disks in close binaries, because the local growth rate is comparable to the orbital frequency towards the outer edge of the disk but decreases rapidly inwards. Paper #1 examined this issue in a simplified global model where the destabilizing terms vary with position. We found that the global growth rate is essentially equal to the maximum local rate, provided that the latter is smoothed over a radial range equal to the distance that the destabilized wave propagates at its group speed in one growth time. Thus, in an accretion disk, waves would grow rapidly in the outer parts but would propagate both inwards and outwards at a maximum group speed of order the disk thickness divided by the orbital period.

Goodman, Jeremy↗

A magnetospheric magnetic field model with flexible current systems driven by independent physical parameters

A tilt-dependent magnetic field model of the Earth's magnetosphere with variable magnetopause standoff distance is presented. Flexible analytic representations for the ring and cross-tail currents, each composed of the elements derived from the Tsyganenko and Usmanov (1982) model, are combined with the fully shielded vacuum dipole configurations of Voigt (1981). Although the current sheet does not warp in the y-z plane, changes in the shape and position of the neutral sheet with dipole tilt are consistent with both MHD equilibrium theory and observations. In addition, there is good agreement with observed Delta B profiles and the average equatorial contours of magnetic field magnitude. While the dipole field is rigorously shielded within the defined magnetopause, the ring and cross-tails currents are not similarly confined, consequently, the model's region of validity is limited to the inner magnetosphere. The model depends on four independent external parameters. We present a simple but limited method of simulating several substorm related magnetic field changes associated with the disrupion of the near-Earth cross-tail current sheet and collapse of the midnight magnetotail field region. This feature further facilitates the generation of magnetic field configuration time sequences useful in plasma convection simulations of real magnetospheric events.

Hilmer, Robert V.↗

Experimental Plans for Subsystems of a Shock Wave Driven Gas Core Reactor

This Contractor Report proposes a number of plans for experiments on subsystems of a shock wave driven pulsed magnetic induction gas core reactor (PMI-GCR, or PMD-GCR pulsed magnet driven gas core reactor). Computer models of shock generation and collision in a large-scale PMI-GCR shock tube have been performed. Based upon the simulation results a number of issues arose that can only be addressed adequately by capturing experimental data on high pressure (approx.1 atmosphere or greater) partial plasma shock wave effects in large bore shock tubes ( 10 cm radius). There are three main subsystems that are of immediate interest (for appraisal of the concept viability). These are (1) the shock generation in a high pressure gas using either a plasma thruster or pulsed high magnetic field, (2) collision of MHD or gas dynamic shocks, their interaction time, and collision pile-up region thickness, and (3) magnetic flux compression power generation (not included here).

Kazeminezhad, F.↗

Cometary environments; Proceedings of Symposium 5, Workshop IV, and Topical Meeting of the 27th COSPAR Plenary Meeting, Espoo, Finland, July 18-29, 1988

Papers on the environment of comets are presented, covering topics such as constraints on the interstellar dust model of comet dust, dust particles and comet nuclei, models of cometary nuclei, subliming gas in the near-nuclear layer of the comet coma, the nucleus and rotation of Comet Halley, the surface albedo of comet nuclei, observational studies on Comet Halley, comet simulations, comet ion composition, and chemical abundances in comets. Additional topics include the O 1D and H2O production rate from comets, collisional coma models, the gas coma of Comet Giacobini-Zinner, IR properties of rough comet grains, studies of Comet Halley by Giotto, the impact of large dust particles on the Vega spacecraft, and carbonaceous materials as components of comet dust. Also, consideration is given to comet plasma boundaries, the comet ionopause, the solar wind-comet interaction, MHD turbulence and particle acceleration in a mass-loaded solar wind, combined first and second order Fermi acceleration at comets, discrete wave packets upstream from the earth and comets, and the visual appearance of comets under varying solar wind conditions.

Gombosi, T. I.↗

Magnetohydrodynamic models of coronal transients in the meridional plane. I - The effect of the magnetic field

The propagation of coronal transients through the lower corona in the meridional plane is studied on the basis of numerical solutions of the time-dependent MHD equations of motion. The importance of the topology of the initial coronal magnetic field is demonstrated by considering two configurations: one which is essentially radial (open), and another which is essentially parallel to the solar surface (closed). The effect of coronal field magnitude is investigated by using values for the beta parameter of 1 and 0.1 at the coronal base. The consequences of including radiative losses and of varying the polytropic index are also considered. The solar event is simulated by a step-function increase in the pressure at the base of an initially hydrostatic atmosphere by a factor of 5, lasting for a period of 5 min.

Steinolfson, R. S.↗

Aerospace Applications of Non-Equilibrium Plasma

Nonequilibrium plasma/non-thermal plasma/cold plasmas are being used in a wide range of new applications in aeronautics, active flow control, heat transfer reduction, plasma-assisted ignition and combustion, noise suppression, and power generation. Industrial applications may be found in pollution control, materials surface treatment, and water purification. In order for these plasma processes to become practical, efficient means of ionization are necessary. A primary challenge for these applications is to create a desired non-equilibrium plasma in air by preventing the discharge from transitioning into an arc. Of particular interest is the impact on simulations and experimental data with and without detailed consideration of non-equilibrium effects, and the consequences of neglecting non-equilibrium. This presentation will provide an assessment of the presence and influence of non-equilibrium phenomena for various aerospace needs and applications. Specific examples to be considered will include the forward energy deposition of laser-induced non-equilibrium plasmoids for sonic boom mitigation, weakly ionized flows obtained from pulsed nanosecond discharges for an annular Hall type MHD generator duct for turbojet energy bypass, and fundamental mechanisms affecting the design and operation of novel plasma-assisted reactive systems in dielectric liquids (water purification, in-pipe modification of fuels, etc.).

Nonequilibrium Plasma↗

Influence of the Convection Electric Field Models on Predicted Plasmapause Positions During Magnetic Storms

In the present work, we determine how three well documented models of the magnetospheric electric field, and two different mechanisms proposed for the formation of the plasmapause influence the radial distance, the shape and the evolution of the plasmapause during the geomagnetic storms of 28 October 2001 and of 17 April 2002. The convection electric field models considered are: Mcllwain's E51) electric field model, Volland-Stern's model and Weimer's statistical model compiled from low-Earth orbit satellite data. The mechanisms for the formation of the plasmapause to be tested are: (i) the MHD theory where the plasmapause should correspond to the last-closed- equipotential (LCE) or last-closed-streamline (LCS), if the E-field distribution is stationary or time-dependent respectively; (ii) the interchange mechanism where the plasmapause corresponds to streamlines tangent to a Zero-Parallel-Force surface where the field-aligned plasma distribution becomes convectively unstable during enhancements of the E-field intensity in the nightside local time sector. The results of the different time dependent simulations are compared with concomitant EUV observations when available. The plasmatails or plumes observed after both selected geomagnetic storms are predicted in all simulations and for all E-field models. However, their shapes are quite different depending on the E-field models and the mechanisms that are used. Despite the partial success of the simulations to reproduce plumes during magnetic storms and substorms, there remains a long way to go before the detailed structures observed in the EUV observations during periods of geomagnetic activity can be accounted for very precisely by the existing E-field models. Furthermore, it cannot be excluded that the mechanisms currently identified to explain the formation of "Carpenter's knee" during substorm events, will', have to be revised or complemented in the cases of geomagnetic storms.

Pierrard, V.↗

Resistive wave dissipation on magnetic inhomogeneities Normal modes and phase mixing

Numerical solutions of the linearized, resistive MHD equations indicate the presence of two characteristic forms of wave decay on magnetic inhomogeneities. Solutions for relatively long-wavelength disturbances and large values of the resistivity have the behavior of decaying normal modes. In the opposite limits, phase mixing becomes important, with the eventual build-up of large spatial gradients in which the resistive terms dominate. This composite behavior is shown to be conceptually consistent with analytic results, which predict that the complete solution consists of the sum of collective and noncollective contributions. The numerical simulations go beyond analytic theory by defining where each contribution prevails and computing wave decay when phase mixing is effective. The dispersion relation in the normal-mode regime is also determined analytically, with some approximations, and is in good agreement with the numerical solution. The decay time of the normal modes varies with Lundquist number at S1/6, while the phase-mixed decay time scales as S1/3.

Steinolfson, R. S.↗

Hybrid simulations of whistler waves generation and current closure by a pulsed tether in the ionosphere

The dynamic response of a magnetized collisionless plasma to an externally driven, finite size, sudden switch-on current source across the magnetic field has been studied using a two dimensional hybrid code. It was found that the predominant plasma response was the excitation of whistler waves and the formation of current closure by induced currents in the plasma. The results show that the current closure path consists of: (a) two antiparallel field-aligned current channels at the end of the imposed current sheet; and (b) a cross-field current region connecting these channels. The formation of the current closure path occured in the whistler timescale much shorter than that of MHD and the closure region expanded continuously in time. The current closure process was accompanied by significant energy loss due to whistler radiation.

Chang, C. L.↗

Comparison of Computed and Measured Performance of a Pulsed Inductive Thruster Operating on Argon Propellant

Pulsed inductive plasma accelerators are electrodeless space propulsion devices where a capacitor is charged to an initial voltage and then discharged through a coil as a high-current pulse that inductively couples energy into the propellant. The field produced by this pulse ionizes the propellant, producing a plasma near the face of the coil. Once a plasma is formed if can be accelerated and expelled at a high exhaust velocity by the Lorentz force arising from the interaction of an induced plasma current and the magnetic field. A recent review of the developmental history of planar-geometry pulsed inductive thrusters, where the coil take the shape of a flat spiral, can be found in Ref. [1]. Two concepts that have employed this geometry are the Pulsed Inductive Thruster (PIT)[2, 3] and the Faraday Accelerator with Radio-frequency Assisted Discharge (FARAD)[4]. There exists a 1-D pulsed inductive acceleration model that employs a set of circuit equations coupled to a one-dimensional momentum equation. The model was originally developed and used by Lovberg and Dailey[2, 3] and has since been nondimensionalized and used by Polzin et al.[5, 6] to define a set of scaling parameters and gain general insight into their effect on thruster performance. The circuit presented in Fig. 1 provides a description of the electrical coupling between the current flowing in the thruster I1 and the plasma current I2. Recently, the model was upgraded to include an equation governing the deposition of energy into various modes present in a pulsed inductive thruster system (acceleration, magnetic flux generation, resistive heating, etc.)[7]. An MHD description of the plasma energy density evolution was tailored to the thruster geometry by assuming only one-dimensional motion and averaging the plasma properties over the spatial dimensions of the current sheet to obtain an equation for the time-evolution of the total energy. The equation set governing the dynamics of the coupled electrodynamic-current sheet system is composed of first-order, coupled ordinary differential equations that can be easily solved numerically without having to resort to much more complex 2-D finite element plasma simulations.

Polzin, Kurt A.↗

Magnetohydrodynamic simulation of the coronal transient associated with the solar limb flare of 1980, June 29, 18:21 UT

The spatial and temporal behavior of excess and depleted density regions of the coronal transient accompanying the west limb solar flare of 18:21 UT on June 29, 1980, is modeled mathematically on the basis of SMM-X-ray-polychrometer data, and the results are compared to observations made with the radio spectrograph at Harvard Radio Astronomy Station at Fort Davis and with the NCAR/HAO Mark III K-coronameter at Mauna Loa. Input data for the model include T(max) = about 20 x 10 to the 6th K, n(max) = about 4 x 10 to the 11th/cu cm, and an assumed ejection velocity of 200 km/sec. Computations using an improved 2D nonplane MHD model are carried out for locally open and closed magnetic topologies. The spatially wide, large-amplitude, temporarily steepened MHD wave predicted by the model for both magnetic topologies is shown to agree well with the observations, except that the predicted density enhancement exceeded observed values by at least 50 percent for the closed field and by a factor of 3 for the open field. This discrepancy is seen as an indication that the mass emitting soft X-rays was confined in closed-field regions near the sun during the obervation period.

Wu, S. T.↗

State of the AGN: Progress Toward Understanding Black Hole Accretion Processes

Accretion of plasma onto black holes power some of the most powerful systems in the cosmos. Supermassive black holes at the centers of galaxies represent the high-mass limit of these objects and so account for the most luminous accretors. As a result, their influence spans vast spatial and temporal scales of cosmic phenomena: intracluster heating, intergalactic media, galactic feedback and star formation, kiloparsec-scale jets/outflows, variability over time scales of minutes to centuries, and luminous multi-wavelength electromagnetic emission extending all the way down to its event horizon. Their intrigue is heightened by the fact that they lie at the intersection of various physical laws---e.g., general relativistic gravity, magnetohydrodynamics, radiation, high-energy particle physics, thermodynamics, and photo-ionization---which all must be reconciled to arrive at a fundamental understanding and probe for new physics, like tests of general relativity. These physics ingredients must be incorporated into simulations performed somehow on dynamical scales ranging from that of the event horizon to parsec-scales. Fortunately, new computational and theoretical techniques---such as GPU computing, radiation transport, and and novel gridding techniques---are enabling progress to larger scales, more degrees of freedom, and even to binary systems. Some of the topics we will survey include recent progress on simulating the relationship between the disk-jet interaction, how tilted black holes behave, radiation-dominated flow, and how binary AGN affect the standard picture of black hole accretion. Along the way, we will highlight how new technologies have enabled these scientific rewards. Future directions and open questions will be provided to inspire discussion and interaction during the session.

LISA↗

Dynamic MHD modeling of the solar wind disturbances during the August 1972 events

A time-dependent one-dimensional MHD theoretical model is tested by using plasma and magnetic field observations of Pioneer 9 and Pioneer 10 during the August 1972 events on the sun and in the interplanetary medium. These spacecraft were nearly aligned along a common heliocentric radius during these events, considered now to be the most spectacular and best-documented events during solar cycle 20. The observations of Pioneer 9 at 0.78 AU were used as input for the theoretical model. The plasma and magnetic field forcing functions were superimposed upon a preexisting ambient solar wind at this inner boundary, and the response was simulated as far as 8 AU. The simulated output at 2.2 AU is compared directly with the Pioneer 10 observations at 2.2 AU. Qualitative comparison is good, although several limitations of the one-dimensional theory are noted.

Dryer, M.↗

The Generation of Smooth High Speed Solar Wind From Plume-Interplume Mixing

Ulysses has shown that fast solar wind is extremely smooth, with a variance of less than 5%, in contrast to slow wind with a variance of approximately 30%. Now UVCS has produced the surprising result that the flow speed within coronal holes, the source of fast wind, is not at all smooth. Specifically, Giordano et al. (ApJ, v531, L79-L82, 2000) report that at 1.7 R(sub SUN) the interplume flow speed is typically more than twice the plume flow speed. Other less direct evidence supports this same result, with speeds from less than 300 to over 1000 km/s reported at approximately 5 R(sub SUN). This presents the paradox of how strongly differing plume and interplume flow speeds can exist near the Sun and be absent far from the Sun. The only answer is that plume and interplume material or momentum must be strongly mixed and that the mixing must occur mainly inside 0.3 AU to be consistent with Helios observations of smooth fast wind. Pressure balance structures (PBSs) and He abundance anomalies (Reisenfeld, et al., GRL, v26(13), 1805-1808, 1999) have been identified as interplanetary remnants of plumes, implying momentum mixing is the dominant coronal process. One possible source for plume/interplume momentum mixing is MHD Kelvin-Helmholtz (KH) shear instabilities occurring on the velocity shear interfaces. The velocity shear is a source of free energy and KH fluctuations could, through nonlinear cascade, provide the forcing required for the plasma oscillations (Cranmer, ApJ, v532, 1197-1208, 2000) reported to exist in coronal holes. The physical properties in coronal holes are now sufficiently well known that we can show plume/interplume shear interfaces become unstable to the KH instability at 5 - 10 R(sub SUN). The KH dispersion relation can be used to analyze marginal stability, the most unstable wavelengths, and linear growth rates. Numerical simulations can be used to verify results from the linear analysis and study the nonlinear development of KH modes. Here we will describe the marginal stability criterion for the KH instability, how and where this condition is met in the corona, and the possible character of the resultant fluctuations.

Suess, Steve↗