Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Magnetic reconnection with large separatrix angles

The magnetic reconnection process is studied here using incompressible MHD simulations with different inflow boundary conditions and different magnetic Reynolds numbers R(m). The angle between the magnetic separatrices is in steady state reconnection depends mainly on the normal magnetic field on the inflow boundary. In steady state nonuniform reconnection with large separatrix angles, field-aligned plasma jets appear slightly downstream of the magnetic separatrices. The field-aligned plasma jet are stronger when R(m) is larger. Each field-aligned plasma jet consists of two parts: a slow shock and a fast-mode compressional wave. The slow shock converts the magnetic energy into plasma kinetic energy by acceleration and heating. The fast-mode compressional wave decelerates the plasma to a smaller outflow speed and heats it further. Nearly all the magnetic energy flowing into the diffusion region is converted into other forms. The length and width of the diffusion region depend on the values of the reconnection rate, R(m), and the normal magnetic field on the inflow boundary.

Yan, M.↗
NASA NTRS1992

MHD turbulence, reconnection, and test-particle acceleration

We examine homogeneous MHD turbulence and turbulent magnetic reconnection as possible mechanisms for accelerating cosmic ray particles. Test particle calculations are performed using fields from MHD simulations, and initially Maxwellian particle distributions are shown to evolve into power-law distributions. Simple estimates for both the maximum energy attainable and the mean energies of the accelerated particles are fairly successful and are consistent with timescales for flares and cosmic rays.

Gray, Perry C.↗
NASA NTRS1993

The thermal structure of the magnetized solar transition region

The detailed thermal structure of the magnetized solar transition region, as measured by its differential emission measure DEM(T), is unknown. Proposals have been made that envision a significant lower-temperature contribution to the energy balance from cross-field (ion) heat flux. In this paper, we describe a self-consistent 2D MHD simulation (including the full effects of anisotropic thermal conduction) of a conceptual model due to Athay (1990). We display the detailed irregular thermal and magnetic structure of the transition region and demonstrate that the predicted DEM agrees with observations, particularly in the T less than 10 exp 5 K regime where previous theories had difficulty.

Mok, Y.↗
NASA NTRS1993

The differential emission measure of nested hot and cool magnetic loops

The detailed thermal structure of the magnetized solar transition region, as measured by its differential emission measure (DEM(T)), is poorly known. Building on the fact that the solar surface is strongly magnetized and thereby structured, proposals have been made that envision a significant lower-temperature contribution to the energy balance from (ion) heat flux across an arcade of different temperature loops. In this paper, we describe a self-consistent 2D MHD simulation, which includes the full thermal effects of parallel stability and anisotropic conduction, of a nested-loop model of the thermal and magnetic structure of the transition region. We then demonstrate that the predicted DEM agrees with observations in the conceptually elusive T less than 10 exp 5 K regime.

Van Hoven, G.↗
NASA NTRS1994

Influence of plasma sheet thickening toward the tail flanks on magnetotail stability and dynamics

On the basis of resistive MHD simulations we investigate the stability and dynamic evolution of the magnetotail for configurations with various degrees of plasma sheet thickening from midnight toward the tail flanks, associated with a similar increase of the B(sub z) component. This increase is varied by factors between 2 and 3.5, while the magnitude of B(sub z) at midnight is left unchanged. The increase of B(sub z) has a strong effect on the spatial extend of resistive tearing and reconnection across the tail and on the unstable evolution. For stronger thickening the cross-tail extent of the near-Earth neutral line, formed by reconnection, and correspondingly the amount of reconnected magnetic flux get reduced and the dipolarization effects earthward of the neutral line are more concentrated near midnight. An increase of B(sub z) from midnight toward the flanks by more than a factor of about 3 can possibly even suppress the resistive tearing instability. This indicates the possibility of a stability transition due to a reduction of B(sub z) in the flank regions in response to changes in the solar wind, even when the magnitude of B(sub z) near midnight and the fluctuation level in that region are unchanged (provided that a sufficient fluctuation level exists).

Birn, Joachim↗
NASA NTRS1994

A mechanism to produce a dawn-dusk component of plasma flow during magnetic reconnection in the magnetotail

Magnetic reconnection between antiparallel field lines in the magnetotail is generally thought to produce plasma acceleration in the earthward-tailward direction. However, measurements of the plasma velocity in the magnetotail during substorm activity sometimes reveal a dawn-dusk component of plasma flow. In this paper, we show that a dawn-dusk component of plasma acceleration may be produced during reconnection if the neutral line is not perpendicular to the magnetic field. In this case, Magnetohydrodynamic (MHD) simulations show that reconnection between antiparallel field lines will initially produce plasma acceleration that is nearly parallel to the neutral line because the magnetic tension force is not opposed by a pressure gradient force in this direction. As the magnetic field topology evolves to a steady state, the plasma flow direction rotates until it is nearly parallel to the plane that initially contained the antiparallel magnetic field lines before reconnection (hereafter referred to as the initial field plane). However, the time required to reach a steady state (typically several hundred seconds in the magnetotail region) may be greater than the time during which the reconnection process is active. Consequently, bursts of plasma flow with a dawn-dusk component may occur in the magnetotail. The initial acceleration along the neutral line depends on the angle theta (sub B) between the neutral line and the initial field plane, with the largest burst of plasma flow along the neutral line occuring when theta (s ub B) = 45 degs.

Hawkins, J. G.↗
NASA NTRS1994

MHD modeling of magnetotail instability for localized resistivity

We present results of a three-dimensional magnetohydrodynamic (MHD) simulation of magnetotail evolution initiated by a sudden occurrence or increase of spatially localized resistivity as the major expected concequence of some localized microinstability. Because of the absence of a quantitative model, possible variations of resistivity levels with current density, or the reduction thereof, are not incorporated in the present investigation. The emphasis of the study is on an investigation of the changes to the overall evolution brought about by this localization, in particular, on the disruption and diversion of the cross-tail current and the nonlinear evolution of the magnetotail instability. The immediate consequences of the occurrence of the localized resistance and the resulting electric field are a reduction and diversion of the electric current around the region of high resistivity, associated with an increase of B(sub z) ('dipolarization') at the earthward edge and a decrease of B(sub z) at the tailward edge of this region. These effects, however, are localized and do not involve a reduction of the total cross-tail current and hence do not lead to the global development of a 'substorm current wedge,' which includes not only the reduction of the cross-tail current but also the buildup of a global field-aligned current system of 'regional 1' type (toward the Earth on the dawnside and away on the duskside of the tail). Such signatures develop at a later time, as consequences of a three-dimensional tearing instability, which is triggered by the occurrence of the resistivity. These features are found in combination with plasmoid formation and ejection, quite similar to results of earlier simulations with uniform resistivity. Differences are found in the timescale of the evolution, which tends to be shorter for localized resistivity, and in the propagation of the dipolarization effects in the equatorial plane. Whereas for uniform resistivity the temporal increase in B(sub z) tends to propagate tailward, apparently due to a pileup effect in the near-Earth region, an earthward propagation is found for the localized resistivity. This propagation results from an earthward convection of the gradient in B(sub z), which is set up by the reconnection process further tailward.

Hesse, Michael↗
NASA NTRS1994

Tearing mode instability in a multiple current sheet system

The tearing mode and magnetic reconnection are studied for multiple current sheet systems by two-dimensional magnetohydrodynamic (MHD) simulations. Both the linear and nonlinear evolution of this process are anaylsed for laminar perturbations. The results illustrate the existence of a linear regime with a symmetric and antisymmetric mode and agree with previous analytic results (Otto and Birk, 1992). The nonlinear evolution shows a number of interesting new features and may explain some properties in corresponding studies of turbulent reconnection. For wavelengths larger than twice the current sheet separation the evolution of antisymmetric modes leads to an entire reconfiguration of the magnetic field and converts a major portion of the magnetic energy into kinetic energy. Antisymmetric modes with smaller wavelengths and symmetric modes are found to saturate. The influence of the value of the resistivity on the reconnection rate decreases in the nonlinear evolution, and the ratio of current sheet separation to wavelength seems to be of major importance. A comparion of the dynamics of periodic current sheets with the evolution of only two current sheets indicates that some of the results for the periodic system also apply to the evolution of only two interacting current sheets. The results are discussed with respect to observations of large-scale plasma and magnetic field reconfigurations in the magnetosheath and near the Earth's bow shock.

Yan, M.↗
NASA NTRS1994

The population of the magnetosphere by solar winds ions when the interplanetary magnetic field is northward

We have examined some possible entry mechanisms of solar wind ions into the magnetosphere by calculating the trajectories of thousands of non-interacting ions in the magnetic and electric fields from a three dimensional global magnetohydrodynamic (MHD) simulation of the magnetosphere and the magnetosheath, under northward interplanetary magnetic field (IMF) conditions. Particles, launched in the solar wind, entered the magnetosphere and formed the low latitude boundary layer (LLBL), plasma sheet and a region of trapped particles near the Earth. The densities and temperatures we obtained in these regions were realistic, with the exception of trapped particle densities. The dominant entry mechanism was convection into the magnetosphere on reconnecting field lines.

Richard, Robert L.↗
NASA NTRS1995

Magnetic reconnection in the presence of sheared flow and density asymmetry: Applications to the Earth's magnetopause

Classical models of magnetic reconnection consist of a small diffusion region bounded by two symmetric slow shocks, across which the plasma is accelerated. Asymmetries often present in space plasmas are sheared plasma flow and dissimilar plasma densities on the two sides of current sheets. In this paper, we investigate magnetic reconnection in the presence of a shear flow and an asymmetric density across the current sheet using two-dimensional magnetohydrodynamic (MHD) simulations. The results demonstrate that magnetic reconnection can occur only for a plasma flow velocity (in the frame of the X line) which is below the Alfven speed in each inflow region. This limits the velocity of the X line to a certain range for a given flow shear and provides an upper limit to the total velocity shear at which reconnection ceases to operate. Depending on the direction of the flow in the adjacent inflow region, the effects from the sheared flow and from the density asymmetry will compete with or enhance each other in respect to the magnitude and location of the currents which bound the outflow regions. The results are applied to the dayside and flank regions of the magnetosphere. For the dayside region where the magnetosheath flow is slow, the magnetic field transition region is thin and the accelerated flow is earthward of the sharp current layer (magnetopause). At the flanks tailward of the X line, shear flow and density asymmetry effects compete making the magnetic field transition layer broad with the high-speed flow contained within the transition region which explains corresponding observations. At the flanks sunward of the X line, shear flow and density asymmetry effects enhance each other and lead to a strong current sheet on the magnetosheath side of the accelerated flow. The total volume affected by magnetic reconnection is much larger than the steady state region. A large bulge region precedes the steady state region. Qualitatively, the bulge and the steady state region have similar signatures and both can explain observations. We provide criteria in order to distinguish between the bulge and the steady state region in observations.

La Belle-Hamer, A. L.↗
NASA NTRS1995

Flow properties along field lines in a 3-D tilted-dipole geometry

A 3-D MHD simulation of a global, tilted-dipole solar wind flow pattern is analyzed to determine flow properties along individual magnetic field lines. In the model, flow conditions near the Sun are chosen to provide a reasonable match to the interplanetary configuration prevailing during the recent south polar passage by Ulysses, i.e., a streamer belt inclined approximately 30 deg to the solar equator and speeds ranging from 325-800 km/s. Field lines all across the stream pattern are traced from 1 to 10 AU by following the motion of marker particles embedded in the flow. It is found that those field lines threading the core of the interaction region are subject to significant latitudinal and relative longitudinal displacement over this range of heliocentric distance. Thus, observations taken at a fixed latitude in the inner solar system sample, over the course of a solar rotation, field lines which connect to a range of latitudes in the outer heliosphere. Maps of the field line displacements are presented to help visualize these connections. In addition, it is found that depending upon the location relative to the CIR structure, the radial evolution of fluid elements frozen to different field lines can deviate considerably from that of the canonical solar wind. That is, for selected subsets of field lines, large speed changes (not just at shocks) can be experienced; the density variation can be far from 1/r(exp 2), and the magnetic field intensity need not decay monotonically with distance.

Pizzo, V. J.↗
NASA NTRS1995

CME stimulated by eruptive prominence

A model of CME arising due to drift motion in the corona in the presence of an eruptive prominence is presented. Magnetic field configuration is in accordance with a model of inverse polarity. In a region where a magnetic pressure of the filament magnetic field is greater than a gas pressure, plasma motion can be assumed as a drift motion. Its characteristic is such that the further one gets from the filament, the higher is drift velocity of a plasma. That sort of motion leads to a rarefaction of plasma and formation of a cavity around the filament. But a current strength in eruptive prominences estimated from observations is such that the region b is less than 1 has a limited size. Near the boundary b = 1 plasma deceleration is occurred and as a result of it coronal density is increasing. Plasma condensation near the surface b = 1 leads to formation of a dense envelope which can be collated with an outer loop of CME. Two dimensional numerical MHD simulation displays a process of cavity and loop formation. If a current is large enough, two compact regions of compressed dense matter arise at both sides of the rising filament and two narrow jets are developed. This scenario, perhaps, corresponds to CMEs in which a top of the loop is faint or it is absent at all.

Filippov, B. P.↗
NASA NTRS1995

Coronal mass ejections and the evolution of the large-scale corona

Coronal mass ejections (CMEs) are dynamic, large-scale events in the solar corona that expel plasma and magnetic fields into the solar wind. The structure of the large-scale corona influences CME occurrence; for example, CMEs are most frequently observed to be disruptions of coronal or helmet streamer configurations on the Sun. We have investigated the evolution of the large-scale corona in the presence of differential rotation, using time-dependent magnetohydrodynamic (MHD) simulations in two and three dimensions. An equilibrium coronal configuration is first computed by superimposing a Parker solar wind flow on an initial potential field, the equations are integrated in time until a steady-state is reached. When differential rotation is introduced, we find that the configuration does not reach a steady state; the closed field regions (helmet streamers) disrupt recurrently, ejecting plasmoids into the solar wind. Our results suggest that differential rotation may be one of the mechanisms by which mass ejections are initiated. We will describe the evolution of the polarization brightness (pB) for the resulting configurations, as well as the propagation of the ejected plasmoids in the solar wind, and we will compare the effects of differential rotation with other possible initiation mechanisms, such as the emergence of new magnetic flux.

Linker, J. A.↗
NASA NTRS1995

Scientific visualization tools for the ISTP project: Mission planning, data analysis and model interpretation

Visualization tools are being developed to meet the challenges of mission planning and data analysis presented by the International Solar-Terrestrial Physics (ISTP) program. ISTP encompasses a large number of spacecraft, multiple ground-based observatories, and several theoretical investigations, with the goal of understanding the global behavior of the solar wind/magnetosphere/ionosphere system. The tools include three-dimensional displays of key boundaries in geospace along with spacecraft trajectories, which can be animated and synchronized to universal time. Magnetic field models and MHD simulation results can be invoked to reveal the magnetic topology or to identify magnetic conjunctions between spacecraft and/or ground-based facilities. Simultaneous displays of satellite trajectories, spacecraft-borne observations, and model predictions are available to facilitate data processing and interpretation efforts. The current status of these tools is described, and their implementation at the ISTP Science Planning and Operations Facility and distribution to the entire ISTP community are discussed.

Peredo, M.↗
NASA NTRS1997

Studies of Plasma Flow Past Jupiters Satellite Io

We have investigated the interaction of Io, Jupiter's innermost Galilean satellite, with the Io plasma torus, and the interaction of Ganymede with the corotating Jovian plasma. With the successful insertion of the Galileo spacecraft into orbit around Jupiter, many new observations have been made of the Jovian magnetosphere. Some of the most exciting results thus far have been in regards to Jupiter's satellites, Io and Ganymede. In both cases the large perturbations to the background (Jovian) magnetic field have been consistent with the satellites' possession of an intrinsic magnetic field. The gravity measurements implying a differentiated core at both Io and Ganymede makes internal generation of a magnetic field by dynamo action in these satellites plausible, and, in the case of Ganymede, the identification of an intrinsic field is apparently unambiguous. For Io the situation is less clear, and further analysis is necessary to answer this important question. During the past year, we have used time-dependent three-dimensional magnetohydrodynamic (MHD) simulations to study these plasma-moon interactions. The results from these simulations have been used directly in the analysis of the Galileo magnetometer data. Our primary emphasis has been on the Io interaction, but we recently presented results on the Ganymede interaction as well. In this progress summary we describe our efforts on these problems to date.

Linker, Jon A.↗
NASA NTRS1996

Latitudinal Dependence of the Radial IMF Component: Coronal Imprint

Measurements by Ulysses have confirmed that there is no significant gradient with respect to heliomagnetic latitude in the radial component, B(sub r,) of the interplanetary magnetic field. In the corona, the plasma, beta is much less than 1, except directly above streamers, so longitudinal and latitudinal gradients in field strength will relax due to the transverse magnetic pressure gradient force as the solar wind carries magnetic flux away from the Sun. This happens quickly enough so that the field is essentially uniform by 5 - 10 solar radius, apparently remaining so as it is carried to beyond 1 AU. Here, we illustrate the coronal relaxation with a qualitative physical argument and by reference to a detailed Magneto HydroDynamics (MHD) simulation.

Suess, S. T.↗
NASA NTRS1997

Cross-field Current Instability for Substorm Expansions

The funding provided by the above-referenced NASA grant has enabled us: (1) to investigate the quasi-linear evolution of the IWI [Lui et al., 1993] and that of the generalized MTSI/IWI [Yoon and Lui, 1993], (2) to carry out the linear analysis of the LHDI to elucidate the difference between it and the MTSI/PM instability [Yoon et al., 1994], (3) to conduct some preliminary nonlocal analyses of the MTSI [Lui et al., 1995] and the IWI [Yoon and Lui, 1996] modes, (4) to study low-frequency shear-driven instability and its nonlinear evolution, which might compete with the CCI [Yoon et al., 1996], and (5) to study the evolution of current sheet during late substorm growth phase by means of 2-D Hall-MHD simulation in order to obtain a better understanding of the current sheet equilibrium crucial for CCI theory [Yoon and Lui, 1997].

Lui, Anthony↗
NASA NTRS1998

Theoretical Investigation of the High-Altitude Cusp Region using Observations from Interball and ISTP Spacecraft

A fundamental goal of magnetospheric physics is to understand the transport of plasma through the solar wind-magnetosphere-ionosphere system. To attain such an understanding, we must determine the sources of the plasma, the trajectories of the particles through the magnetospheric electric and magnetic fields to the point of observation, and the acceleration processes they undergo enroute. This study employed plasma distributions observed in the near-Earth plasma sheet by Interball and Geotail spacecraft together with theoretical techniques to investigate the ion sources and the transport of plasma. We used ion trajectory calculations in magnetic and electric fields from a global Magnetohydrodynamics (MHD) simulation to investigate the transport and to identify common ion sources for ions observed in the near-Earth magnetotail by the Interball and Geotail spacecraft. Our first step was to examine a number of distribution functions and identify distinct boundaries in both configuration and phase space that are indicative of different plasma sources and transport mechanisms. We examined events from October 26, 1995, November 29-30, 1996, and December 22, 1996. During the first event Interball and Geotail were separated by approximately 10 R(sub E) in z, and during the second event the spacecraft were separated by approximately 4(sub RE). Both of these events had a strong IMF By component pointing toward the dawnside. On October 26, 1995, the IMF B(sub Z) component was northward, and on November 1-9-30, 1996, the IMF B sub Z) component was near 0. During the first event, Geotail was located near the equator on the dawn flank, while Interball was for the most part in the lobe region. The distribution function from the Coral instrument on Interball showed less structure and resembled a drifting Maxwellian. The observed distribution on Geotail, on the other hand, included a great number of structures at both low and high energies. During the third event (December 22, 1996) both spacecraft were in the plasma sheet and were separated bY approximately 20 R(sub E) in the y direction. During this event the IMF was southward.

Ashour-Abdalla, Maha↗