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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Thermally conductive magnetohydrodynamic flows in helmet-streamer coronal structures

The behavior of thermally conductive plasma flows in helmet-streamer coronal structures is investigated within the framework of the axisymmetric nonrotating one-fluid MHD model. Continuous subsonic-supersonic solutions satisfying observed boundary conditions at the sun as well as the vanishing of the temperature at infinity are obtained and presented. Special attention is paid to the combined effects of conductive flow (and corresponding thermal force) and rapidly diverging magnetic field on the critical points. In this, the heliocentric distance of the neutral point determining the separation between closed and open field lines (cusp) is treated as a free parameter. These thermally conductive solutions are contrasted with those provided by corresponding isothermal models.

Cuperman, S.↗
NASA NTRS1990

Plasma density enhancement at the Comet Halley diamagnetic cavity boundary

A one-dimensional single-fluid MHD model was used by Cravens (1989) to predict the existence of a narrow (50 km) layer of enhanced plasma density at the boundary of the Comet Halley diamagnetic cavity. The existence of such a layer was confirmed by measurements made by the Giotto ion mass spectrometer (Goldstein et al., 1989). Here, the time-dependent coupled continuity equations are solved for several species including H3O(+), H2O(+), OH(+), O(+), NH4(+), NH3(+), CH4(+), and CH3(+). For several ion species, factors affecting the magnitude of this enhancement and its relationship to the thickness of the transition layer are investigated. For example, ion species with short chemical lifetimes are shown to have smaller density enhancements than species with longer lifetimes. The cometocentric distance of the cavity boundary also strongly affects the magnitude of the density enhancement; the enhancement increases with increasing distance.

Keller, C. N.↗
NASA NTRS1991

Numerical simulation of extended corona

A three-dimensional, time-dependent MHD model is presented for the study of coronal dynamics. The model, written in spherical coordinates, extends from the solar surface and was developed with two major issues in mind, namely for interpretation of various steady state and evolutionary dynamical structures in the corona. In order to achieve these objectives, two different numerical techniques are used to seek solutions for these two different, but related, problems; steady state structures and evolutionary structures. These two numerical techniques are: (1) relaxation technique for steady state structures; and (2) FICE (Full-Implicit-Continuous-Eulerian) technique for evolutionary structures. To illustrate this model, numerical results are given for examples of both the steady state and evolutionary structure of the corona which show the additional physical features which cannot be shown by a two-dimensional model.

Wu, S. T.↗
NASA NTRS1991

On dipolarization and its relation to the substorm current wedge

The near-earth properties of magnetospheric substorms are investigated using the Hesse and Birn (1991) self-consistent 3D MHD model of magnetotail activity. The simulation shows that an increase of Bz, a decrease of the cross-tail current, and an increase of region-1 field-aligned currents occurs roughly simultaneously, consistent with the empirical picture described by McPherron et al. (1973). It was also found that the region where the cross-tail current is reduced and diverted is located tailward from the dipolarization region. The results indicate that the apparent tailward propagation of the dipolarization effects does not support an earlier suggestion that the underlying process is initiated further earthward.

Hesse, Michael↗
NASA NTRS1992

A case study of plasma structure in the dusk sector associated with enhanced magnetospheric convection

Consideration is given to a case study based on a combination of ground whistler and satellite measurements of thermal plasma density which provides additional evidence that the abrupt western edge of the bulge region of the magnetosphere, reported earlier from whistlers, is a real phenomenon. The present data and previous MHD modeling work suggest that this distinctive feature develops during periods of steady or declining substorm activity, when dense plasma previously carried sunward under the influence of enhanced convection activity begins to rotate with the earth at angular velocities that decrease with increasing L value and becomes spirallike in form. Whistler data are used to identify a narrow dense plasma feature, separated from the main plasmasphere and extending sunward into the late afternoon sector at L values near the outer observed limits of the main plasmasphere and extending sunward into the edge of the main bulge, found by both whistler stations to be at about 1800 MLT, appeared to be quasi-stationary in sun-earth coordinates during the prevailing conditions of gradually declining geomagnetic agitation.

Carpenter, D. L.↗
NASA NTRS1990

Formation of slow shock pairs associated with coronal mass ejections

The formation of a forward-reverse slow shock pair in the solar corona is presently simulated by an MHD model that uses the Rankine-Hugoniot solution to calculate the flow-property jumps at all shock crossings. The shocks divide the solution-domain into several continuous flow regions whose respective governing characteristics are solved by the method of characteristics. The plasma impact compresses the plasma near the front of the coronal mass ejection (CME); as the CME-associated slow shock pair moves outwards in interplanetary space, it evolves into a pair of fast shocks. All three phenomena are eventually manifested in interplanetary space as a magnetic cloud accompanied by a fast shock pair, with a forward shock preceding the cloud and a reverse shock appearing either within or behind the cloud.

Whang, Y. V.↗
NASA NTRS1992

Validity of the field line resonance expansion

The field line resonance expansion is studied with a magnetohydrodynamic (MHD) model applicable to earth's magnetotail. Using a Frobenius series solution, this expansion is found to be inconsistent. The inconsistency is seen by evaluating the terms in the model equations that the expansion neglects as small. Their effect is to provide a solution more singular than that determined by the expansion procedure. A tentative alternative expansion procedure, retaining these neglected terms approximately, yields a nonsingular solution, suggesting that resonant mode conversion rather than field line resonance is the appropriate phenomenon.

Hansen, P. J.↗
NASA NTRS1993

Termination shock - Solar cycle variations of location and speed

Plasma and magnetic field data from Voyagers 1 and 2 over an 11-year period (1978-1988) together with a one-fluid MHD model are used to study the motion of the termination shock. The observed large fluctuations in solar wind parameters cause variations in the speed, in the jump conditions, and in the location of the termination shock. If the interstellar stagnation pressure is equivalent to a magnetic field of 0.5 nT plus draping enhancements, the termination shock moves between 88 AU and 102 AU. The location of the termination shock is anticorrelated with the sunspot number. The termination shock is very strong with pressure ratios of a few thousands, the fast Mach number is about 12, and the plasma beta ratio jumps from about 0.1 to about 14. During the declining phase of the solar cycle, the shock speed is often greater than or comparable to the spacecraft speed. If the crossing of the termination shock by a spacecraft occurs during the declining phase, the spacecraft will cross the termination shock 3 or more times over a period of about 4 years.

Whang, Y. C.↗
NASA NTRS1994

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

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

Sojka, J. J.↗
NASA NTRS1994

Growing 'Alfvenic' modes in the upstream region of Saturn

Recent studies of low-frequency electromagnetic waves upstream of the Saturn bow shock have shown that these waves, in contrast to those at Earth, are observed not in one, but in at least two, distinct frequency bands. The results of wave mode identification based on the Hall-magnetohydrodynamic (MHD) model of plasma and observed wave polarization suggested that these waves propagate in the high beta intermediate mode. However, the underlying instability was not unambiguously determined. In the present paper we use the full electromagnetic dispersion relation derived from linear Vlasov theory in order to examine which of the plasma modes, with observed properties, are unstable in an isotropic Maxwellian plasma in the presence of backstreaming proton beams consistent with Voyager 2 observations at Saturn. As a result we find that the unstable 'Alfvenic' beam mode, as well as resonant and non-resonant fast magnetosonic modes have properties consistent with the data. Moreover, we find that in contrast to the Earth's upstream waves, at Saturn no 'kinetic' normal mode can account for the observed magnetic polarization.

Orlowski, D. S.↗
NASA NTRS1994

Chaotic dynamics of corotating magnetospheric convection

The corotating plasma convection system of the Jovian magnetosphere is analyzed. The macroscopic (mhd) model introduced by Summers and Mu, (1992) that incorporates the effects of microdiffusion is extended by including previously neglected density effects. We reduce the governing partial differential equations to a third-order ordinary differential system by the Galerkin technique of mode truncation. We carry out such a severe truncation partly in the interests of tractability, and leave open the question of the efficacy of adding additional modes. Exhaustive numerical integrations are carried out to calculate the long-term solutions, and we discover that a rich array of plasma motions is possible, dependent on the value of the height-integrated ionospheric Pederson conductivity Sigma. If Sigma is less than a certain critical value Sigma(sub c), then plasma motion can be expected to be chaotic (or periodic), while if Sigma is greater than Sigma(sub c), then steady state convection is expected. In the former case, whether the plasma motion is chaotic or periodic (and, if periodic, the magnitude of the period) can be very sensitive to the value of Sigma. The value of Sigma(sub c), which is a function of a parameter q that occurs in the assumed form of the stationary radial profile (varies as L(exp -q) of the plasma mass per unit magnetic flux, lies well within the accepted range of values of Sigma for Jupiter, i.e. Sigma greater than or equal to 0.1 mho and less than or equal to 10 mho.

Summers, Danny↗
NASA NTRS1995

Reflection and transmission of GMIR shock at the heliopause and their relation to the 2- and 3-kHz radio emissions

We use Voyager 2 plasma and magnetic field data together with a one-fluid MHD model to study the interactions of the 1991 Global Merged Interaction Region (GMIR) shock with the heliopause. The 1991 GMIR is an extraordinarily large global solar wind structure in radial, longitudinal and latitudinal extents. It has a strong shock at the leading edge. After its penetration through the termination shock, the GMIR shock first propagates through the subsonic solar wind, then interacts with the heliopause. The interaction produces a transmitted shock propagating outward in the interstellar medium, and a reflected shock propagating backward in the subsonic solar wind. We identify the reflected shock and the transmitted shock as the possible source of the radio noise detected at Voyagers. The plasma frequency behind the reflected and the transmitted shock can be, respectively, responsible for the 2- and 3-kHz radio emissions. The two bands of radio noise are emitted from sources on both sides of the heliopause starting at about the same time. If the emission is generated by f(sub p)-radiation then the heliopause is located at R approximately 130 AU. If the emission is generated by 2f(sub p)-radiation the n R approximately 150 AU. Because the relative speed of the interstellar plasma with respect to the sun appears to be sub-Alfvenic, it is very unlikely there is a fast-bow shock of the heliosphere.

Whang, Y. C.↗
NASA NTRS1995

Heating of the solar middle chromosphere by large-scale electric currents

A global resistive, two-dimensional, time-dependent magnetohydrodynamic (MHD) model is used to introduce and support the hypothesis that the quiet solar middle chromosphere is heated by resistive dissipation of large-scale electric currents which fill most of its volume. The scale height and maximum magnitude of the current density are 400 km and 31.3 m/sq m, respectively. The associated magnetic field is almost horizontal, has the same scale height as the current density, and has a maximum magnitude of 153 G. The current is carried by electrons flowing across magnetic field lines at 1 m/s. The resistivity is the electron contribution to the Pedersen resitivity for a weakly ionized, strongly magnetized, hydrogen gas. The model does not include a driving mechanism. Most of the physical quantities in the model decrease exponentially with time on a resistive timescale of 41.3 minutes. However, the initial values and spatial; dependence of these quantities are expected to be essentially the same as they would be if the correct driving mechanism were included in a more general model. The heating rate per unit mass is found to be 4.5 x 10(exp 9) ergs/g/s, independent of height and latitude. The electron density scale height is found to be 800 km. The model predicts that 90% of the thermal energy required to heat the middle chromosphere is deposited in the height range 300-760 km above the temperature minimum. It is shown to be consistent to assume that the radiation rate per unit volume is proportional to the magnetic energy density, and then it follows that the heating rate per unit volume is also proportional to the energy from the photosphere into the overlying chromosphere are briefly discussed as possible driving mechanisms for establishing and maintaining the current system. The case in which part of or all of the current is carried by protons and metal ions, and the contribution of electron-proton scattering to the current are also considered, with the conclusion that these effects do not change the qualitative prediction of the model, but probably change the quantitative predictions slightly, mainly by increasing the maximum magntiude of the current density and magnetic field to at most approximately 100 mA/m and approximately 484 G, respectively. The heating rate per unit mass, current density scale height, magnetic field scale height, temperatures, and pressures are unchanged or are only slightly changed by including these additional effects due to protons and ions.

Goodman, M. L.↗
NASA NTRS1993

Plasmasphere dynamics in the duskside bulge region: A new look at old topic

Data acquired during several multiday periods in 1982 at ground stations Siple, Halley, and Kerguelen and on satellites Dynamics Explorer 1, International Sun Earth Explorer 1, and GEOS 2 have been used to investigate thermal plasma structure and dynamics in the duskside plasmasphere bulge region of the Earth. The distribution of thermal plasma in the dusk bulge sector is difficult to describe realistically, in part because of the time integral manner in which the thermal plasma distribution depends upon on the effects of bulk cross-B flow and interchange plasma flows along B. While relatively simple MHD models can be useful for qualitatively predicting certain effects of enhanced convection on a quiet plasmasphere, such as an initial sunward entrainment of the outer regions, they are of limited value in predicting the duskside thermal plasma structures that are observed. Furthermore, use of such models can be misleading if one fails to realize that they do not address the question of the formation of the steep plasmapause profile or provide for a possible role of instabilities or other irreversible processes in plasmapause formation. Our specific findings, which are based both upon the present case studies and upon earlier work, include the following: (1) during active periods the plasmasphere appears to become divided into two entities, a main plasmasphere and a duskside bulge region. (2) in the aftermath of an increase in convection activity, the main plasmasphere tends (from a statistical point of view) to become roughly circular in equatorial cross section, with only a slight bulge at dusk; (3) the abrupt westward edge of the duskside bulge observed from whistlers represents a state in the evolution of sunward extending streamers; (4) in the aftermath of a weak magnetic storm, 10 to 30% of the plasma 'removed' from the outer plasmasphere appears to remain in the afternoon-dusk sector beyond the main plasmasphere. (5) outlying dense plasma structures may circulate in the outer duskside magetosphere for many days following an increase in convection, unless there is extremely deep quieting; (6) a day-night plasmatrough boundary may be identified in equatorial satellite data; (7) factor-of-2-to-10 density irregularities appear near the plasmatrough from the ionosphere at L = 4.6, predominantly bidirectional field aligned and equatorially trapped light ion pitch angle distributions give away to a predominantly isotropic distribution (as seen by DE 1) when the plasma density reaches a level a factor of about 3 below the satured plasmasphere level; (9) some outlying dense plasma structures are effectively detached from the main plasmasphere, while others appear to be connected to that body.

Carpenter, D. L.↗
NASA NTRS1995

Radial alignment simulation of solar wind streams observed by Pioneers 10 and 11 in 1974

A particularly favorable radial lineup between spacecraft in Earth orbit and Pioneers 10 and 11 (near the ecliptic plane at approximately 5.5 AU and approximately 4.5 AU, respectively) occurred in mid-1974, when the solar corona was in a stable and well-defined warped-dipole configuration. The radial alignment study reported here differs from previous applications of the technique in two respects: (1) It is the first time a two-dimensional (2-D) MHD model has been tested over such a lengthy propagation interval; the 2-D capability is crucial for treatment of the nonradial shearing motions occurring across the stream interface. (2) The three-dimensional (3-D) structure observed in the white light corona at that time is related to the systematic patterns of nonradial flow deflections appearing at the Pioneer corotating interaction region (CIR) fronts. Comparison of predicted and observed flows for pairs of streams in two successive rotations reveals that when the parent coronal hole projects far across the solar equator in a predominantly north-south orientation (i.e., nearest the 2-D idealization), the mapping is accurate down to details of the flow structures. But where the spacecraft tracks along a latitudinal boundary of a hole or the associated stream front is inclined at a shallow angle to the equator, the numerical projections deviate systematically from the observations. Among the sources of error are 3-D dynamical interactions neglected in the model, differential rotation effects, and slow temporal evolution of the coronal structures. A characteristic pattern of north-south and east-west deflections indicative of the 3-D geometry appears across the CIRs, but these patterns evidently reflect primarily the local, as opposed to global, orientation of the 3-D stream fronts. Such patterns appear common in CIRs observed by Pioneer during this period. These findings thus hold special relevance for the analysis of Ulysses observations, since the present coronal configuration is similar to that of 1974.

Pizzo, V. J.↗
NASA NTRS1995

Forward and reverse shocks in the outer heliosphere: Observations from Voyager 2

Observations from Voyager 2 as it moved from 10 to 14 deg S heliographic latitude in the period from 1992 through 1994 were used to gather statistics on the relative number of forward and reverse shocks. These results can be used to compare with observations from the Ulysses spacecraft which moved from 6 deg S to 70 deg S heliographic latitude during that time period. The Ulysses observations are in agreement with a 3-D, MHD model of the evolution of a steady tilted-dipole solar wind flow configuration prevalent in 1993. The model predicts and the Ulysses observations confirm a preponderance of reverse shocks at Ulysses latitudes poleward of streamer-belt latitudes. A preliminary scan of the Voyager data supports the complementary prediction of the model that forward fronts should dominate at large heliocentric distances near the heliographic equatorial plane during the same time period.

Lazarus, A. J.↗
NASA NTRS1995

Transition of slow shocks to fast shocks

As a CME pushes its way through preceding slower solar wind, large disturbances in the interaction region may evolve to form transient MHD shocks. The shocks can be slow shocks in the coronal space but they appear as fast shocks near 1 AU. We use a polytropic MHD model to study the transition of slow shocks to fast shocks in an interaction region sandwiched between a faster solar wind and a slower solar wind near the equatorial plane. The polytropic index controls the radial increase of Beta outside the interaction region. The transition process is attributed chiefly to the increase of Beta and theta to a lesser degree to the radial increase of the shock angle 0. Under the initial condition of small Beta and theta near 0.1 AU, the interaction region evolves to form a pair of slow shocks inside 0.15 AU. As the interaction region convects outward. the increases of B and theta cause a transition of the shock system between 0.15 and 0.3 AU from a slow shock pair to a double shock pair consisting of both slow and fast shocks. As the system moves outward, and Beta continue to increase, the fast shock grows stronger, and the slow shock becomes weaker. Eventually. the slow shocks fade away, and the shock system finally evolves to a pair of fast shocks. Parametric studies of the transition process are carried out for shocks formed over a wide range of disturbances.

Whang, Y. C.↗
NASA NTRS1997

Effects of Neutral Density on Energetic Ions Produced Near High-Current Hollow Cathodes

Energy distributions of ion current from high-current, xenon hollow cathodes, which are essential information to understand erosion phenomena observed in high-power ion thrusters, were obtained using an electrostatic energy analyzer (ESA). The effects of ambient pressure and external flow rate introduced immediately downstream of hollow cathode on ion currents with energies greater than that associated with the cathode-to-anode potential difference were investigated. The results were analyzed to determine the changes in the magnitudes of ion currents to the ESA at various energies. Either increasing the ambient pressure or adding external flow induces an increase in the distribution of ion currents with moderate energies (epsilon less than 25 to 35 eV) and a decrease in the distribution for high energies (epsilon greater than 25 to 35 eV). The magnitude of the current distribution increase in the moderate energy range is greater for a cathode equipped with a toroidal keeper than for one without a keeper, but the distribution in the high energy range does not seem to be affected by a keeper. An MHD model, which has been proposed to describe energetic-ion production mechanism in hollow cathode at high discharge currents, was developed to describe these effects. The results show, however, that this model involves no mechanism by which a significant increase of ion current could occur at any energy. It was found, on the other hand, that the potential-hill model of energetic ion production, which assumes existence of a local maximum of plasma potential, could explain combined increases in the currents of ions with moderate energies and decreases in high energy ions due to increased neutral atom density using a charge-exchange mechanism. The existing, simplified version of the potential-hill model, however, shows poor quantitative agreement with measured ion-current-energy-distribution changes induced by neutral density changes.

Kameyama, Ikuya↗