SEARCH · Search NASA

Results for “Numerical magnetohydrodynamics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS2009

Three-Dimensional Numerical Modeling of Magnetohydrodynamic Augmented Propulsion Experiment

Over the past several years, NASA Marshall Space Flight Center has engaged in the design and development of an experimental research facility to investigate the use of diagonalized crossed-field magnetohydrodynamic (MHD) accelerators as a possible thrust augmentation device for thermal propulsion systems. In support of this effort, a three-dimensional numerical MHD model has been developed for the purpose of analyzing and optimizing accelerator performance and to aid in understanding critical underlying physical processes and nonideal effects. This Technical Memorandum fully summarizes model development efforts and presents the results of pretest performance optimization analyses. These results indicate that the MHD accelerator should utilize a 45deg diagonalization angle with the applied current evenly distributed over the first five inlet electrode pairs. When powered at 100 A, this configuration is expected to yield a 50% global efficiency with an 80% increase in axial velocity and a 50% increase in centerline total pressure.

Turner, M. W.↗

NASA NTRS1991

Lattice Boltzmann model for simulation of magnetohydrodynamics

A numerical method, based on a discrete Boltzmann equation, is presented for solving the equations of magnetohydrodynamics (MHD). The algorithm provides advantages similar to the cellular automaton method in that it is local and easily adapted to parallel computing environments. Because of much lower noise levels and less stringent requirements on lattice size, the method appears to be more competitive with traditional solution methods. Examples show that the model accurately reproduces both linear and nonlinear MHD phenomena.

Chen, Shiyi↗

NASA NTRS1986

Turbulent relaxation to a force-free field-reversed state

The evolution of nonequilibrium initial conditions of an incompressible magnetohydrodynamic Z pinch is described by a three-dimensional, pseudospectral numerical code. Magnetohydrodynamic turbulence develops in the resistive, nonviscous magnetofluid, resulting in the selective decay of the energy relative to the magnetic helicity, at Lundquist numbers of only a few hundred. An interior force-free region grows with time and achieves spontaneous reversal of the toroidal magnetic field at the wall, without the necessity of an external electric field.

Dahlburg, J. P.↗

NASA NTRS1994

Numerical simulation of torus-driven plasma transport in the Jovian magnetosphere

The Rice convection model has been modified for application to the transport of Io-generated plasma through the Jovian magnetosphere. The new code, called the RCM-J, has been used for several ideal-magnetohydrodynamic (MHD) numerical simulations to study how interchange instability causes an initially assumed torus configuration to break up. In simulations that start from a realistic torus configuration but include no energetic particles, the torus disintegrates too quickly (approximately 50 hours). By adding an impounding distribution of energetic particles to suppress the interchange instability, resonable lifetimes were obtained. For cases in which impoundment is insufficient to produce ideal-MHD stability, the torus breaks up predominantly into long fingers, unless the initial condition strongly favors some other geometrical form. If the initial torus has more mass on one side of the planet than the other, fingers form predominatly on the heavy side (which we associate with the active sector). Coriolis force bends the fingers to lag corotation. The simulation results are consistent with the idea that the fingers are formed with a longitudinal thickness that is roughly equal to the latitudinal distance over which the invariant density declines at the outer edges of the initial torus. Our calculations give an average longitudinal distance between plasma fingers of about 15 deg which corresponds to 20 to 30 minutes of rotation of the torus. We point to some Voyager and Ulysses data that are consistent with this scale of torus longitudinal irregularity.

Yang, Y. S.↗

NASA NTRS1982

Numerical experimentation on spherically symmetric one-dimensional magnetohydrodynamic /MHD/ wave propagation

Radial propagation of one-dimensional magnetohydrodynamic (MHD) waves are analyzed numerically on the basis of the Implicit-Continuous-Fluid-Eulerian (ICE) scheme. Accuracy of the numerical method and other properties are tested through the study of MHD wave propagation. The three different modes of MHD waves (i.e., fast-, slow- and Alfven (transverse) mode) are generated by applying physically consistent boundary perturbations derived from MHD compatibility relations. It is shown that the resulting flow following these waves depend upon the relative configurations of the initial magnetic field and boundary perturbations.

Han, S. M.↗

NASA NTRS1993

Broken symmetry in ideal magnetohydrodynamic turbulence

A numerical study of the long-time evolution of a number of cases of inviscid, isotropic, incompressible, three-dimensional fluid, and magneto-fluid turbulence has been completed. The results confirm that ideal magnetohydrodynamic turbulence is non-ergodic if there is no external magnetic field present. This is due essentially to a canonical symmetry being broken in an arbitrary dynamical representation. The broken symmetry manifests itself as a coherent structure, i.e., a non-zero time-averaged part of the turbulent magnetic field. The coherent structure is observed, in one case, to contain about eighteen percent of the total energy.

Shebalin, John V.↗

NASA NTRS1986

Formation of active region and quiescent prominence magnetic field configurations

To investigate the formation of prominences, researchers studied chromospheric mass injection into an overlying coronal dipole magnetic field using a 2-D ideal magnetohydrodynamic (MHD) numerical model. Researchers propose that active region prominences are formed by chromospheric plasmas injected directly into the overlying coronal magnetic field and that quiescent prominences are formed by plasmas evaporated at the interface between spicules and corona. Hence, for the simulation of an active region prominence magnetic field we inject the mass from one side, but use a symmetric mass injection to form a quiescent prominence field configuration. Researchers try to find optimum conditions for the formation of Kippenhahn-Schuluter(K-S)type field configuration for stable support of the injection plasmas. They find that the formation of K-S type field configuration by mass injection requires a delicate balance between injection velocity, density, and overlying magnetic fields. These results may explain why a prominence does not form on every neutral line.

An, C.-H.↗

NASA NTRS1985

Shock drift acceleration

This is a review of the fundamental physics of the interactions of charged particles treated individually while they interact with fast mode magnetohydrodynamic shocks. Numerical simulation and analytical theory are used to develop predictions of the expected characteristics of this process strong upstream anisotropies directed along the magnetic field and downstream anisotropies tending to be peaked more perpendicular to the field; relatively more of the enhancement of higher-energy particles occurring upstream; sensitive dependence on shock normal to magnetic field angle of the efficiency of energization. Observations which display all of the above characteristics are reviewed. Also discussed is the relationship of shock drift acceleration to the models for stochastic transport of charged particles in the vicinity of shocks. Extensions of this work in both the observational and theoretical approaches are discussed.

Armstrong, Thomas P.↗

NASA NTRS1987

Comparison of MHD simulation for the February 1986 events with interplanetary observations by the spacecraft Sakigake

During the period 3-10 February 1986 a series of major solar flares occurred on the Sun and several intense geomagnetic storms took place on the Earth. To examine the causality between the solar activity and the geomagnetic activity in this period, a magnetohydrodynamic (MHD) numerical simulation was performed using a 2 1/2 -D numerical code. In that period of February 1986, the Japanese spacecraft Sakigake was at 0.84 AU, 57 deg west of the Earth. Besides the in-situ measurements of the interplanetary plasma, Sakigake also provided Doppler scintillation observations. Comparisons between the results of the MHD simulation and the measurements made by the spacecraft Sakigake are presented.

Smith, Z.↗

NASA NTRS1994

Global, quasi-steady dynamics of the distant solar wind 2: Deformation of the heliospheric current sheet

A three-dimensional (3-D) magnetohydrodynamic (MHD) numerical model is used to trace global deformations of the heliospheric current sheet (HCS) caused by large-scale dynamical interactions associated with corotating solar wind flows. Configurations incorporating the tilted-dipole geometry are investigated out past 30 AU for a variety of dipole tilt angles, alpha. Inclusion of the full, 3-D interplanetary dynamics allows north-south displacements and the east-west warping of the HCS by advective corotational effects to be accurately assessed for the first time. It is found that large-scale spatial correlations between velocity and density imposed at the coronal source (i.e., the geometric arrangement whereby regions of high-speed low-density material lying adjacent to areas of slow, dense flow interact obliquely under the influence of solar rotation) result in a distinctive pattern of deformation of the HCS. For an alpha = 30 deg tilted-dipole example, it is shown that typical zonal variations in radial velocity lead to significant folding of the HCS within about 5 AU of the Sun. By 10 AU, additional sharp bends appear near the latitudinal extremes of the HCS surface, where it is overtaken by shock fronts driven by 3-D corotating interaction regions (CIRs). Moreover, the model suggests that inside about 20 AU, major plasma structures are systematically organized, about the HCS, such that the greatest concentrations of material and magnetic field (the centroids of the 3-D CIR structures) are coincident with the folded crests of the HCS (near heliographic latitudes lambda = +/- alpha). Thus, in these circumstances many of the more interesting dynamical features inconveniently lie well away from the heliographic equator. At larger heliocentric distances, where neighboring CIRs begin to interact strongly, the warping of the HCS abates dramatically and the association between folds in the HCS and major field and density concentrations is weakened and ultimately breaks down. On the basis of the model calculations, quasi-steady deformations and related phenomena should figure prominently in the interplanetary medium, except during those periods when the magnetic dipole and solar spin axes are nearly aligned (alpha is approximately less than 10 deg). Even then, they may be of significance if there are any substantial local warps or kinks in the streamer belt.

Pizzo, V. J.↗

NASA NTRS1994

3-D simulation of high-latitude interaction regions: Comparison with Ulysses results

A three-dimensional (3-D) magnetohydrodynamic (MHD) numerical model is used to simulate the global evolution of a steady, tilted-dipole solar wind flow configuration similar to that prevalent in interplanetary space in 1993. Systematic latitudinal changes in the structure of a corotating interaction region (CIR) near 5 AU is shown to agree well with recent Ulysses observations. The abrupt disappearance of forward shocks and continued persistence of reverse shocks poleward of the latitude where Ulysses crossed the southern edge of the coronal streamer belt is explained as a natural consequence of the 3-D flow geometry.

Pizzo, V. J.↗

NASA NTRS1998

Multidimensional MHD Model Studies of the Ionospheres of Venus and Mars

Continuing efforts have been made towards an increased understanding of the solar wind interaction and ionospheric processes at Venus and Mars. This work centered on a systematic development of a new generation of three dimensional magnetohydrodynamic (MHD) numerical code, which models the interaction processes of the solar wind with non-magnetic planets, such as Venus and Mars. We have also worked on a number of different, more specific and discrete studies, as various opportunities arose. We have developed new numerical codes for magnetospheric and cometary studies. As a first step in this process we built an axisymmetric model in which the solar wind interacts with a hard, perfectly conducting sphere. Even that model provided, in certain respects, significant improvements over previous ones.

Nagy, Andrew↗

NASA NTRS2000

Two-Species, 3D, MHD Simulation of Europa's Interaction with Jupiter's Magnetosphere

The interaction of Europa with the Jovian a magnetosphere has been studied by using a two species in ideal magnetohydrodynamic (MHD) numerical model. This model considers the upstream plasma in the Jovian magnetosphere and the molecular oxygen ions in the ionosphere of Europa, separately. We present results a from simulation studies, which take into account impact ionization, recombination, and the effect of a possible induced dipole magnetic field of Europa. The total mass loading of the magnetospheric flow and the ionization frequency used in the model are consistent with the estimates of Europa's ionosphere and atmosphere. The multi-species MHD equations are solved by using a finite volume, high-order, Godunov-type method on an adoptively refined unstructured grid, which allows a detailed modeling of the region near Europa's surface, while still resolving both the upstream region and the satellite's wake. We have paid special attention to the wake of Europa, in order to be able to make comparisons with the Galileo's E4 flyby observations, as well as other model calculations. The calculated escape flux of a O2+ down the tail was found to be about 5.6 x 10(exp 25) s(sup -1).

Liu, Yifan↗

NASA NTRS2016

A Model for Straight and Helical Solar Jets: II. Parametric Study of the Plasma Beta

Context. Jets are dynamic, impulsive, well-collimated plasma events that develop at many different scales and in different layers of the solar atmosphere. Aims. Jets are believed to be induced by magnetic reconnection, a process central to many astrophysical phenomena. Within the solar atmosphere, jet-like events develop in many different environments, e.g. in the vicinity of active regions as well as in coronal holes, and at various scales, from small photospheric spicules to large coronal jets. In all these events, signatures of helical structure and/or twisting/rotating motions are regularly observed. The present study aims to establish that a single model can generally reproduce the observed properties of these jet-like events. Methods. In this study, using our state-of-the-art numerical solver ARMS, we present a parametric study of a numerical tridimensional magnetohydrodynamic (MHD) model of solar jet-like events. Within the MHD paradigm, we study the impact of varying the atmospheric plasma beta on the generation and properties of solar-like jets. Results. The parametric study validates our model of jets for plasma beta ranging from 10(sup 3) to 1, typical of the different layers and magnetic environments of the solar atmosphere. Our model of jets can robustly explain the generation of helical solar jet-like events at various beta less than or equal to 1. We show that the plasma beta modifies the morphology of the helical jet, explaining the different observed shapes of jets at different scales and in different layers of the solar atmosphere. Conclusions. Our results allow us to understand the energisation, triggering, and driving processes of jet-like events. Our model allows us to make predictions of the impulsiveness and energetics of jets as determined by the surrounding environment, as well as the morphological properties of the resulting jets.

Magnetic Field↗

NASA NTRS2016

Do We Know the Actual Magnetopause Position for Typical Solar Wind Conditions?

We compare predicted magnetopause positions at the subsolar point and four reference points in the terminator plane obtained from several empirical and numerical MHD (magnetohydrodynamics) models. Empirical models using various sets of magnetopause crossings and making different assumptions about the magnetopause shape predict significantly different magnetopause positions (with a scatter greater than 1 Earth radius (R (sub E)) even at the subsolar point. Axisymmetric magnetopause models cannot reproduce the cusp indentations or the changes related to the dipole tilt effect, and most of them predict the magnetopause closer to the Earth than non axisymmetric models for typical solar wind conditions and zero tilt angle. Predictions of two global non axisymmetric models do not match each other, and the models need additional verification. MHD models often predict the magnetopause closer to the Earth than the non axisymmetric empirical models, but the predictions of MHD simulations may need corrections for the ring current effect and decreases of the solar wind pressure that occur in the foreshock. Comparing MHD models in which the ring current magnetic field is taken into account with the empirical Lin et al. model, we find that the differences in the reference point positions predicted by these models are relatively small for B (sub z) equals 0 (note: B (sub z) is when the Earth's magnetic field points north versus Sun's magnetic field pointing south). Therefore, we assume that these predictions indicate the actual magnetopause position, but future investigations are still needed.

Samsonov, A. A.↗

NASA NTRS1981

Magnetohydrodynamic interaction of high-speed streams

Numerical solutions of a magnetohydrodynamic model are carried out to describe the nonlinear interaction of corotating high-speed streams near the solar equatorial plane. Two problems are studied. The first problem is to simulate the evolution of an idealized high-speed stream. Numerical solutions are obtained to represent the variations of flow velocity, magnetic field, plasma density, temperature, and conduction heat flux in the interaction region. They demonstrate that the dynamical interaction and heat conduction process are responsible for the thermal structure of a high-speed stream. The second problem deals with the formation of corotating shock waves near the leading edge of a broad stream resulting from the merging of characteristic curves. Corotating shocks do not necessarily occur in pairs; a reverse shock can be formed without a forward shock nearby.

Whang, Y. C.↗

NASA NTRS1972

Numerical calculation of the lunar wake in a magnetohydrodynamic model.

Following previous interpretation, a model is developed which has the solar wind flowing directly into the lunar surface and the magnetic field passing unimpeded through the moon. The collapse of the plasma void behind the moon is calculated numerically using the equations of magnetohydrodynamics. The model is three-dimensional, allowing complete calculations to be made without special restrictions on the magnetic field direction. A cursory quantitative comparison with data shows agreement to be satisfactory. The calculations indicate that the rate of collapse of the void does not depend strongly on the solar wind-magnetic field angle. The penumbral decrease depends very strongly on both the solar wind-magnetic field angle and the angle of passage through the wake.

Beers, B. L.↗

NASA NTRS1994

Oblique magnetohydrodynamic cosmic-ray-modified shocks: Two-fluid numerical simulations

We present the first results of time-dependent two-fluid cosmic-ray (CR) modified magnetohydrodynamic (MHD) shock simulations. The calculations were carried out with a new numerical code for one-dimensional ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks, including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they apppear to have converged to dynamical steady states. We also present results from an initial time-dependent simulation which extends the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescales than parallel shocks.

Frank, Adam↗