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.

SourceResearch connectionFrom yearThrough year

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

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

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 NTRS1993

Theoretical studies of a molecular beam generator

The following is a proposed baseline model that is being developed for the simulation of hydrodynamic generator, which can be converted at a later date to magnetohydrodynamic MHD thruster by adding the necessary electric and magnetic fields. The following development will include the electric and magnetic terms, however, the initial computer program will not include these terms. The analysis that follows is a one species, single temperature model constructed over a domain that is defined by an enclosed region.

Heinbockel, John H.↗

NASA NTRS1994

Numerical Simulation of Cylindrical, Self-field MPD Thrusters with Multiple Propellants

A two-dimensional, two-temperature, single fluid MHD code was used to predict the performance of cylindrical, self-field magnetoplasmadynamic (MPD) thrusters operated with argon, lithium, and hydrogen propellants. A thruster stability equation was determined relating maximum stable J(sup 2)/m values to cylindrical thruster geometry and propellant species. The maximum value of J(sup 2)/m was found to scale as the inverse of the propellant molecular weight to the 0.57 power, in rough agreement with limited experimental data which scales as the inverse square root of the propellant molecular weight. A general equation which relates total thrust to electromagnetic thrust, propellant molecular weight, and J(sup 2)/m was determined using reported thrust values for argon and hydrogen and calculated thrust values for lithium. In addition to argon, lithium, and hydrogen, the equation accurately predicted thrust for ammonia at sufficiently high J(sup 2)/m values. A simple algorithm is suggested to aid in the preliminary design of cylindrical, self-field MPD thrusters. A brief example is presented to illustrate the use of the algorithm in the design of a low power MPD thruster.

Lapointe, Michael R.↗

NASA NTRS1994

Nonlinear evolution of a large-amplitude circularly polarized Alfven wave: High beta

The nonlinear dynamics following saturation of the parametric instabilities of a monochromatic field-aligned large-amplitude circularly polarized Alfven wave is investigated via direct numerical simulation in the case of high plasma beta and no wave dispersion. The magnetohydrodynamic (MHD) code permits nonlinear couplings in the parallel direction to the ambient magnetic field and one perpendicular direction. Compressibility is included in the form of a polytropic equation of state. Turbulent cascades develop after saturation of two coupled oblique three-wave parametric instabilities; one of which is an oblique filamentationlike instability reported earlier. Remnants of the parametric processes, as well as of the original Alfven pump wave, persist during late nonlinear times. Nearly incompressible MHD features such as spectral anisotropies appear as well.

Ghosh, S.↗

NASA NTRS1999

JANNAF 24th Airbreathing Propulsion Subcommittee and 36th Combustion Subcommittee Joint Meeting

Volume 1, the first of three volumes is a compilation of 16 unclassified/unlimited-technical papers presented at the Joint Army-Navy-NASA-Air Force (JANNAF) 24th Airbreathing Propulsion Subcommittee and 36th Combustion Subcommittee held jointly with the 181 Propulsion Systems Hazards Subcommittee. The meeting was held on 18-21 October 1999 at NASA Kennedy Space Center and The DoubleTree Oceanfront Hotel, Cocoa Beach, Florida. Topics covered include overviews of RBCC and PDE hypersonic technology, Hyper-X propulsion ground testing, development of JP-8 for hypersonic vehicle applications, numerical simulation of dual-mode SJ combustion, V&V of M&S computer codes, MHD SJ and Rocket Based Combined Cycle (RBCC) launch vehicle concepts, and Pulse Detonation Engine (PDE) propulsion technology development including fundamental investigations, modeling, aerodynamics, operation and performance.

Fry, Ronald S.↗

NASA NTRS1995

Magnetohydrodynamic turbulence in the solar wind

The fluctuations in magnetic field and plasma velocity in solar wind, which possess many features of fully developed magnetohydrodynamic (MHD) turbulence, are discussed. Direct spacecraft observations from 0.3 to over 20 AU, remote sensing radio scintillation observations, numerical simulations, and various models provide complementary methods that show that the fluctuations in the wind parameters undergo significant dynamical evolution independent of whatever turbulence might exist in the solar photosphere and corona. The Cluster mission, with high time resolution particle and field measurements and its variable separation strategies, should be able to provide data for answering many questions on MHD turbulence.

Goldstein, Melvyn L.↗

NASA NTRS1995

Coronal heating by the resonant absorption of Alfven waves: Wavenumber scaling laws.

The importance of global modes in coronal loop heating is well established. In the present work the scaling of the global-mode resonant heating rate with the perturbation wavenumbers is studied with the numerical solution of the linearized time-dependent MHD equations for a full compressible, low-beta, resistive plasma using an implicit integration scheme. The numerical simulations demonstrate that the dissipation on inhomogeneties in the background Alfven speed occurs in narrow resonant layer with the highest heating rate at the global-mode frequency. The global-mode heating rate H (sub r) was found to scale as H (Sub r) approximately k (sub y) (exp 1.03) when k (sub z) = 0.1, and as H (sub r) approximately k (sub y) (exp -1.93) when k (sub z) = 0.75, where k (sub y) and k (sub z) are the wavenumbers in the perpendicular and parallel to the magnetic field directions, respectively, while the dependence of H (sub r) on k (sub z) is more complex. The quality factor Q of the MHD resonance cavity scales as Q approximately k (sub y) (exp -1.8) for k (sub z) = 0.75 and as Q approximately k (sub y) (exp -1.46) for k (sub z) = 0.1. The numerically determined heating rate scaling, the global-mode fequency, and the quality factor are in good agreement with the analytical linear theory. The magnitude of the perturbed velocities was found to decrease with k (sub y). Assuming typical coronal loop parameters (B (sub 0) = 100-200 G, upsilon (sub A) = 2000-4000 km/s), the Alfven waves can supply the required heating to a low-Q loops.

Ofman, L.↗

NASA NTRS2014

MHD Modelling of Coronal Loops: Injection of High-Speed Chromospheric Flows

Context. Observations reveal a correspondence between chromospheric type II spicules and bright upward-moving fronts in the corona observed in the extreme-ultraviolet (EUV) band. However, theoretical considerations suggest that these flows are probably not the main source of heating in coronal magnetic loops. Aims. We investigate the propagation of high-speed chromospheric flows into coronal magnetic flux tubes and the possible production of emission in the EUV band. Methods. We simulated the propagation of a dense 104 K chromospheric jet upward along a coronal loop by means of a 2D cylindrical MHD model that includes gravity, radiative losses, thermal conduction, and magnetic induction. The jet propagates in a complete atmosphere including the chromosphere and a tenuous cool (approximately 0.8 MK) corona, linked through a steep transition region. In our reference model, the jet initial speed is 70 km per second, its initial density is 10(exp 11) per cubic centimeter, and the ambient uniform magnetic field is 10 G. We also explored other values of jet speed and density in 1D and different magnetic field values in 2D, as well as the jet propagation in a hotter (approximately 1.5 MK) background loop. Results. While the initial speed of the jet does not allow it to reach the loop apex, a hot shock-front develops ahead of it and travels to the other extreme of the loop. The shock front compresses the coronal plasma and heats it to about 10(exp 6) K. As a result, a bright moving front becomes visible in the 171 Angstrom channel of the SDO/AIA mission. This result generally applies to all the other explored cases, except for the propagation in the hotter loop. Conclusions. For a cool, low-density initial coronal loop, the post-shock plasma ahead of upward chromospheric flows might explain at least part of the observed correspondence between type II spicules and EUV emission excess.

Sun: chromosphere↗

NASA NTRS1998

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

Plumes and rays are magnetic field aligned density striations in coronal holes with different values of plasma beta. The overall plasma beta is very small in the low corona but exceeds unity beyond 15-20 solar radius. High speed solar wind reported beyond 0.3 AU is relatively smooth and uniform and known to originate from the much filamented coronal hole. Thus the obvious question is how to generate a smooth solar wind from seemingly filamentary structure. Hence one has to find a mechanism to substantiate this apparent observed (Ulysses) phenomenon. To do this we model plumes as jets (or wakes) of plasma emitted from the solar surface. The shear between a jet and its ambient is known to become unstable to the MHD Kelvin-Helmholtz ("KH") instability if the Alfven Mach number of the jet is greater than one and the uniform external magnetic field is small. Starting with a simple configuration we consider a jet of half thickness R, having uniform density and uniform internal magnetic field. The external medium has also a uniform density and uniform magnetic field. The jet is perturbed at the boundary with a linear amplitude and fixed frequency. We simulate the coronal jet using the 3D ZEUS code. The first results indicate the slab jet is unstable to the MHD KH instability at 5-10 solar radius for some angle of wave propagation. The propagating instability may smooth the filamented flow. It may also produce the entrained Alfvenic fluctuations observed by Ulysses in the high speed wind. We are at present determining the parameters which induce large growth rate. This may clarify the mystery behind the emergence of fast smooth solar wind from very filamentary structures in coronal holes. Also, using the dispersion relation already available for such a flow we obtain some general description of the instability criteria for the KH instability at a jet interface.

Parhi, Shyamsundar↗

NASA NTRS2017

Electron Drift Resonance in the MHD-Coupled Comprehensive Inner Magnetosphere-Ionosphere Model

Relativistic electrons in the outer radiation belt are highly dynamic and respond to interplanetary solar wind structures interacting with the Earth's magnetic field. A known mechanism dictating electron dynamics is the drift-resonant interaction with ultralow frequency (ULF) waves. The present work simulates the ring current and radiation belt electron populations in the bounce-averaged, kinetic Comprehensive Inner Magnetosphere-Ionosphere model coupled with the Block Adaptive Tree Solar Wind Roe-type Upwind Scheme global magnetospheric magnetohydrodynamic (MHD) code using an idealized ULF wave solar wind density driver. ULF waves generated with 10 min periods (at 1.67 mHz frequencies) in the MHD model are characterized and the corresponding energization of electrons and radial transport of electron phase space density is presented. The drift-resonant electron energy is determined in the simulation and is consistent with the electron resonance conditions in dipolar magnetic fields. The present results will be an important component of understanding inner magnetospheric dynamics and how these inner magnetospheric populations interact with ULF waves resulting from interplanetary solar wind structures.

Komar, C. M.↗

NASA NTRS1995

Coronal loop formation resulting from photospheric convection

We have demonstrated the dynamic formation of coronal magnetic loops in three dimensions as a result of horizontal vortex-like convection on the photosphere. Localized plasma motions twist bipolar magnetic field lines which are tied to the dense photosphere by high electrical conductivity. The twists propagate into the corona along the field and create a narrow quasi-toroidal region where the field lines interwind. At the same time, this tubeline region rises in altitude, expands in cross section, and distorts into a slight S shape before settling into an equilibrium state. The MHD stability of such line-tied magnetic loop structures is directly exhibited by this dynamic simulation.

Hoven, G. Van↗

NASA NTRS2024

Feasibility Assessment of Magnetohydrodynamic Control for Aerocapture at Neptune

The Planetary Science Decadal Survey released in 2022 posed a mission to one of the Ice Giants as the top priority for flagship missions for NASA. However, current technologies limit the amount of scientific payload available for future Uranian and Neptunian missions due to the need for fuel for orbit insertion maneuvers. Thus, to maximize the scientific potential of future missions, aerocapture has been heavily researched. While aerocapture simulations using only aerodynamic control have proven enabling for capturing around Ice Giants like Neptune, the deep atmospheric pass requires an aeroshell with robust thermal protection systems (TPS). Magnetohydrodynamically controlled (MHD) aerocapture serves as a potential improvement to the limitations of both fully propulsive orbit insertion and aerodynamically controlled aerocapture. Using NASA tools for modeling planetary exploration missions, both the aerodynamic-only and magnetohydrodynamic aerocapture methods were simulated and compared for identical missions to an Ice Giant, with Neptune chosen as the target planet. After applying a guidance algorithm for both methods, the results showed that magnetohydrodynamics has not only the control authority to successfully capture around Neptune, but also the unique advantage of a shallower atmospheric pass, decreasing the maximum heat load and the required TPS mass.

Aerocapture↗

NASA NTRS1998

Can Kelvin-Helmholtz Instabilities of Jet-Like Structures and Plumes Cause Solar Wind Fluctuations at 1 AU?

The long high latitude sampling of Ulysses provides the opportunity to study fine structures. At latitudes poleward of approx. -60 deg, the solar wind had fluctuations in velocity gradients which were attributed to "microstreams". The data also suggested fluctuations characterized by magnetic plus thermal pressure balance structures ("PBS"). At higher frequencies, MHD turbulence was observed and found to be less evolved than in the ecliptic, but essentially independent of heliographic latitude. It is argued here that microstreams, PBS, and MHD turbulence could all be the remnants of mixing due to shear instabilities associated with plumes and other filamentarystructures ("jets") in coronal holes. To show this we simulate a plume-like jet in the presence of an ambient magnetic field. We find the presence of the ambient field reduces the growth rate of the instability, but the shear between a jet and its ambient still becomes unstable to the MHD Kelvin-Helmholtz ("KH") instability when the shear speed is larger than the largest local magnetosonic speed - a condition probably satisfied for plumes.

Shyamsundar, Parhi↗

NASA NTRS1989

Monte Carlo simulation of particle acceleration at astrophysical shocks

A Monte Carlo code was developed for the simulation of particle acceleration at astrophysical shocks. The code is implemented in Turbo Pascal on a PC. It is modularized and structured in such a way that modification and maintenance are relatively painless. Monte Carlo simulations of particle acceleration at shocks follow the trajectories of individual particles as they scatter repeatedly across the shock front, gaining energy with each crossing. The particles are assumed to scatter from magnetohydrodynamic (MHD) turbulence on both sides of the shock. A scattering law is used which is related to the assumed form of the turbulence, and the particle and shock parameters. High energy cosmic ray spectra derived from Monte Carlo simulations have observed power law behavior just as the spectra derived from analytic calculations based on a diffusion equation. This high energy behavior is not sensitive to the scattering law used. In contrast with Monte Carlo calculations diffusive calculations rely on the initial injection of supra-thermal particles into the shock environment. Monte Carlo simulations are the only known way to describe the extraction of particles directly from the thermal pool. This was the triumph of the Monte Carlo approach. The question of acceleration efficiency is an important one in the shock acceleration game. The efficiency of shock waves efficient to account for the observed flux of high energy galactic cosmic rays was examined. The efficiency of the acceleration process depends on the thermal particle pick-up and hence the low energy scattering in detail. One of the goals is the self-consistent derivation of the accelerated particle spectra and the MHD turbulence spectra. Presumably the upstream turbulence, which scatters the particles so they can be accelerated, is excited by the streaming accelerated particles and the needed downstream turbulence is convected from the upstream region. The present code is to be modified to include a better description of particle scattering (pitch-angle instead of hard-sphere) and as iterative procedure for treating the self-excitation of the MHD turbulence.

Campbell, Roy K.↗

NASA NTRS1991

Relaxation processes in a low-order three-dimensional magnetohydrodynamics model

The time asymptotic behavior of a Galerkin model of 3D magnetohydrodynamics (MHD) has been interpreted using the selective decay and dynamic alignment relaxation theories. A large number of simulations has been performed that scan a parameter space defined by the rugged ideal invariants, including energy, cross helicity, and magnetic helicity. It is concluded that time asymptotic state can be interpreted as a relaxation to minimum energy. A simple decay model, based on absolute equilibrium theory, is found to predict a mapping of initial onto time asymptotic states, and to accurately describe the long time behavior of the runs when magnetic helicity is present. Attention is also given to two processes, operating on time scales shorter than selective decay and dynamic alignment, in which the ratio of kinetic to magnetic energy relaxes to values 0(1). The faster of the two processes takes states initially dominant in magnetic energy to a state of near-equipartition between kinetic and magnetic energy through power law growth of kinetic energy. The other process takes states initially dominant in kinetic energy to the near-equipartitioned state through exponential growth of magnetic energy.

Stribling, Troy↗

NASA NTRS1999

Can Kelvin-Helmholtz Instabilites of Jet-Like Structures and Plumes Cause Solar Wind Fluctuations at 1 AU?

The long high-latitude sampling of Ulysses provides the opportunity to study fine structures. At latitudes poleward of about -60 degrees the solar wind had fluctuations in velocity gradients which were attributed to "microstreams." The data also suggested fluctuations characterized by magnetic plus thermal pressure balance structures ('PBS'). At higher frequencies, MHD turbulence was observed and found to be less evolved than it is in the ecliptic but essentially independent of heliographic latitude. It is argued here that microstreams, PBS, and MHD turbulence could all be the remnants of mixing due to shear instabilities associated with plumes and other filamentary structures ("jets") in coronal holes. To show this, we simulate a plume-like jet in the presence of an ambient magnetic field. We find that the presence of the ambient field reduces the growth rate of the instability, but the shear between a jet and its ambient still becomes unstable to the MHD Kelvin-Helmholtz instability when the shear speed is larger than the largest local magnetosonic speed, a condition probably satisfied for plumes.

Parhi, S.↗

NASA NTRS1997

Laboratory Plasma Source as an MHD Model for Astrophysical Jets

The significance of the work described herein lies in the demonstration of Magnetized Coaxial Plasma Gun (MCG) devices like CPS-1 to produce energetic laboratory magneto-flows with embedded magnetic fields that can be used as a simulation tool to study flow interaction dynamic of jet flows, to demonstrate the magnetic acceleration and collimation of flows with primarily toroidal fields, and study cross field transport in turbulent accreting flows. Since plasma produced in MCG devices have magnetic topology and MHD flow regime similarity to stellar and extragalactic jets, we expect that careful investigation of these flows in the laboratory will reveal fundamental physical mechanisms influencing astrophysical flows. Discussion in the next section (sec.2) focuses on recent results describing collimation, leading flow surface interaction layers, and turbulent accretion. The primary objectives for a new three year effort would involve the development and deployment of novel electrostatic, magnetic, and visible plasma diagnostic techniques to measure plasma and flow parameters of the CPS-1 device in the flow chamber downstream of the plasma source to study, (1) mass ejection, morphology, and collimation and stability of energetic outflows, (2) the effects of external magnetization on collimation and stability, (3) the interaction of such flows with background neutral gas, the generation of visible emission in such interaction, and effect of neutral clouds on jet flow dynamics, and (4) the cross magnetic field transport of turbulent accreting flows. The applicability of existing laboratory plasma facilities to the study of stellar and extragalactic plasma should be exploited to elucidate underlying physical mechanisms that cannot be ascertained though astrophysical observation, and provide baseline to a wide variety of proposed models, MHD and otherwise. The work proposed herin represents a continued effort on a novel approach in relating laboratory experiments to astrophysical jet observation. There exists overwhelming similarity among these flows that has already produced some fascinating results and is expected to continue a high pay off in future flow similarity studies.

Mayo, Robert M.↗

NASA NTRS1986

Models for application of radiation boundary condition for MHD waves in collapse calculations

The problem of reflection of magnetohydrodynamic (MHD) waves at the boundary of a numerical grid has to be resolved in order to obtain reliable results for the end state of the (isothermal) collapse of a rotating, magnetic protostellar cloud. Since the goal of investigating magnetic braking in collapse simulations is to see if the transport of angular momentum via alfven waves is large enough to solve the angular momentum problem an approximation that artificially suppresses large amplitudes in the MHD waves can be self-defeating. For this reason, four alternate methods of handling reflected waves where no assumptions are made regarding the amplitudes of the waves were investigated. In order to study this problem (of reflection) without interference from other effects these methods were tried on two simpler cases. The four methods are discussed.

Vanajakshi, C. T.↗