Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma fluid 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Extended Magnetohydrodynamics with Embedded Particle-in-Cell Simulation of Ganymede's Magnetosphere

We have recently developed a new modeling capability to embed the implicit particle-in-cell (PIC) model iPIC3D into the Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme magnetohydrodynamic (MHD) model. The MHD with embedded PIC domains (MHO-EPIC) algorithm Is a two-way coupled kinetic-fluid model. As one of the very first applications of the MHD-EPIC algorithm, we simulate the Interaction between Jupiter's magnetospherlc plasma and Ganymede's magnetosphere. We compare the MHO-EPIC simulations with pure Hall MHD simulations and compare both model results with Galileo observations to assess the Importance of kinetic effects In controlling the configuration and dynamics of Ganymede's magnetosphere. We find that the Hall MHD and MHO-EPIC solutions are qualitatively similar, but there are significant quantitative differences. In particular. the density and pressure inside the magnetosphere show different distributions. For our baseline grid resolution the PIC solution is more dynamic than the Hall MHD simulation and it compares significantly better with the Galileo magnetic measurements than the Hall MHD solution. The power spectra of the observed and simulated magnetic field fluctuations agree extremely well for the MHD-EPIC model. The MHO-EPIC simulation also produced a few flux transfer events (FTEs) that have magnetic signatures very similar to an observed event. The simulation shows that the FTEs often exhibit complex 3-0 structures with their orientations changing substantially between the equatorial plane and the Galileo trajectory, which explains the magnetic signatures observed during the magnetopause crossings. The computational cost of the MHO-EPIC simulation was only about 4 times more than that of the Hall MHD simulation.

Toth, Gabor↗
NASA NTRS1982

Collisionless electrostatic interchange instabilities

The linear Vlasov dispersion equation for electrostatic plasma instabilities driven by gravity and weak density gradients perpendicular to a uniform magnetic field is derived and solved numerically. Two interchange instabilities emerge: the well-known fluid mode at long wavelengths and a kinetic model at wavelengths short compared with the ion gyroradius. The properties of both instabilities are studied, as well as the effects of gravity on the universal and lower-hybrid density drift instabilities. The results show that the kinetic interchange generally has a larger growth rate than the fluid interchange instability, indicating that, whenever the latter is present in a collisionless plasma, the former may also be found.

Gray, S. P.↗
NASA NTRS2011

Background and Pickup Ion Velocity Distribution Dynamics in Titan's Plasma Environment: 3D Hybrid Simulation and Comparison with CAPS T9 Observations

In this report we discuss the ion velocity distribution dynamics from the 3D hybrid simulation. In our model the background, pickup, and ionospheric ions are considered as a particles, whereas the electrons are described as a fluid. Inhomogeneous photoionization, electron-impact ionization and charge exchange are included in our model. We also take into account the collisions between the ions and neutrals. The current simulation shows that mass loading by pickup ions H(+); H2(+), CH4(+) and N2(+) is stronger than in the previous simulations when O+ ions are introduced into the background plasma. In our hybrid simulations we use Chamberlain profiles for the atmospheric components. We also include a simple ionosphere model with average mass M = 28 amu ions that were generated inside the ionosphere. The moon is considered as a weakly conducting body. Special attention will be paid to comparing the simulated pickup ion velocity distribution with CAPS T9 observations. Our simulation shows an asymmetry of the ion density distribution and the magnetic field, including the formation of the Alfve n wing-like structures. The simulation also shows that the ring-like velocity distribution for pickup ions relaxes to a Maxwellian core and a shell-like halo.

Lipatov, A. S.↗
NASA NTRS1967

A study of seeded rotating arcs.

Neutral carrier gas and slightly ionized plasma binary mixture model for nonhomogeneous seeding effects on rotating arc performance

ROTATING FLUID↗
NASA NTRS1994

A time-dependent gyro-kinetic model of thermal ion upflows in the high-latitude F region

Ample evidence supports the significance of the high-latitude ionospheric contribution to magnetospheric plasma. Assuming flux conservation along a flux tube, the upward field-aligned ion flows observed in the magnetosphere require high-latitude ionospheric field-aligned ion upflows of the order of 10(exp 8) to 10(exp 9)/sq cm/s. Since radar and satellite observations of high-latitude F region flows at times exceed this flux requirement by an order of magnitude, the thermal ionospheric upflows are not simply the ionospheric response to a magnetospheric flux requirement. Several ionospheric ion upflow mechanisms have been proposed, but simulations based on fluid theory do not reproduce all the observed features of ionospheric ion upflows. Certain asymmetries in the statistical morphology of high-latitude F region ion upflows suggest that the ion upflows may be generated by ion-neutral frictional heating. We developed a single-component (O(+)), time-dependent gyro-kinetic model of the high-latitude F region response to frictional heating in which the neutral exobase is a discontinuous boundary between fully collisional and collisionless plasmas. The concept of a discontinuous neutreal exobase and the assumption of a constant and uniform polarization electric field reduce the ion velocity distribution function, from which we can compute the ion density, parallel velocity, parallel and perpendicular temperature, and parallel flux. Using our model, we simulated the response of a convecting flux tube between 500 km and 2500 km to various frictional heating inputs; the results were both qualitatively and quantitatively different from fluid model results, which may indicate an inadequacy of the fluid theory approach. The gyro-kinetic frictional heating model responses to the various simulations were qualitatively similar: (1) initial perturbations of all the modeled parameters propagated rapidly up the flux tube, (2) transient values of the ion parallel velocity, temperature, and flux exceeded 3 km/s, 2 x 10(exp 4) K, and 10(exp 9)/sq cm/s, respectively, (3) a second transient regime developed wherein the parallel temperature drops to very low values (a few hundred Kelvins), and (4) well after heating ceased, large parallel temperatures and large downward parallel velocities and fluxes developed as the flux tube slowly returned to diffusive equilibrium. The ion velocity distributions during the simulation are often non-Maxwellian and are sometimes composed of two distinct ion populations.

Loranc, M.↗
NASA NTRS1997

Predictive Modeling in Plasma Reactor and Process Design

Research continues toward the improvement and increased understanding of high-density plasma tools. Such reactor systems are lauded for their independent control of ion flux and energy enabling high etch rates with low ion damage and for their improved ion velocity anisotropy resulting from thin collisionless sheaths and low neutral pressures. Still, with the transition to 300 mm processing, achieving etch uniformity and high etch rates concurrently may be a formidable task for such large diameter wafers for which computational modeling can play an important role in successful reactor and process design. The inductively coupled plasma (ICP) reactor is the focus of the present investigation. The present work attempts to understand the fundamental physical phenomena of such systems through computational modeling. Simulations will be presented using both computational fluid dynamics (CFD) techniques and the direct simulation Monte Carlo (DSMC) method for argon and chlorine discharges. ICP reactors generally operate at pressures on the order of 1 to 10 mTorr. At such low pressures, rarefaction can be significant to the degree that the constitutive relations used in typical CFD techniques become invalid and a particle simulation must be employed. This work will assess the extent to which CFD can be applied and evaluate the degree to which accuracy is lost in prediction of the phenomenon of interest; i.e., etch rate. If the CFD approach is found reasonably accurate and bench-marked with DSMC and experimental results, it has the potential to serve as a design tool due to the rapid time relative to DSMC. The continuum CFD simulation solves the governing equations for plasma flow using a finite difference technique with an implicit Gauss-Seidel Line Relaxation method for time marching toward a converged solution. The equation set consists of mass conservation for each species, separate energy equations for the electrons and heavy species, and momentum equations for the gas. The sheath is modeled by imposing the Bohm velocity to the ions near the walls. The DSMC method simulates each constituent of the gas as a separate species which would be analogous in CFD to employing separate species mass, momentum, and energy equations. All particles including electrons are moved and allowed to collide with one another with the stipulation that the electrons remain tied to the ions consistent with the concept of ambipolar diffusion. The velocities of the electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. These benefits come at a price in terms of computational time and memory. The DSMC and CFD are made as consistent as possible by using similar chemistry and power deposition models. Although the comparison of CFD and DSMC is interesting, the main goal of this work is the increased understanding of high-density plasma flowfields that can then direct improvements in both techniques. This work is unique in the level of the physical models employed in both the DSMC and CFD for high-density plasma reactor applications. For example, the electrons are simulated in the present DSMC work which has not been done before for low temperature plasma processing problems. In the CFD approach, for the first time, the charged particle transport (discharge physics) has been self-consistently coupled to the gas flow and heat transfer.

Hash, D. B.↗
NASA NTRS1992

Current-driven plasma acceleration versus current-driven energy dissipation. III - Anomalous transport

In the present paper the linear stability description and weak turbulence theory are used to develop a second order description of wave-particle transport and anomalous dissipation. The goal is to arrive at anomalous transport coefficients that can be readily included in fluid flow codes. In particular, expressions are derived for the heating rates of ions and electrons by the unstable waves and for the electron-wave momentum exchange rate that controls the anomalous resistivity effect. Comparative calculations were undertaken assuming four different saturation models: ion trapping, electron trapping, ion resonance broadening, and thermodynamic bound. A foremost finding is the importance of the role of electron Hall parameter in scaling the level of anomalous dissipation for the parameter range of the MPD thruster plasma. Polynomial expressions of the relevant transport coefficients cast solely in terms of macroscopic parameters are also obtained for inclusion in plasma fluid codes for the self-consistent numerical simulation of real thruster flows including microturbulent effects.

Choueiri, Edgar Y.↗
NASA NTRS1990

Kinetic effects on Alfven wave nonlinearity. II - The modified nonlinear wave equation

A previously developed Vlasov theory is used here to study the role of resonant particle and other kinetic effects on Alfven wave nonlinearity. A hybrid fluid-Vlasov equation approach is used to obtain a modified version of the derivative nonlinear Schroedinger equation. The differences between a scalar model for the plasma pressure and a tensor model are discussed. The susceptibilty of the modified nonlinear wave equation to modulational instability is studied. The modulational instability normally associated with the derivative nonlinear Schroedinger equation will, under most circumstances, be restricted to left circularly polarized waves. The nonlocal term in the modified nonlinear wave equation engenders a new modulational instability that is independent of beta and the sense of circular polarization. This new instability may explain the occurrence of wave packet steepening for all values of the plasma beta in the vicinity of the earth's bow shock.

Spangler, Steven R.↗
NASA NTRS1985

In situ measurements of the plasma bulk velocity near the Io flux tube

The flow around the Io flux tube was studied by analyzing the eleven spectra taken by the Voyager 1 Plasma Science (PLS) experiment in its vicinity. The bulk plasma parameters were determined using a procedure that uses the full response function of the instrument and the data in all four PLS sensors. The mass density of the plasma in the vicinity of Io is found to be 22,500 + or - 2,500 amu/cu cm and its electron density is found to be 1500 + or - 200/cu cm. The Alfven speed was determined using three independent methods; the values obtained are consistent and taken together yield V sub A = 300 + or - 50 km/sec, corresponding to an Alfven Mach number of 0.19 + or - 0.02. For the flow pattern, good agreement was found with the model of Neubauer (1980), and it was concluded that the plasma flows around the flux tube with a pattern similar to the flow of an incompressible fluid around a long cylinder obstacle of radius 1.26 + or - 0.1 R sub Io.

Barnett, A.↗
NASA NTRS2013

Radial Profiles of Plasma Electron Characteristics in a Low-Power Arcjet

Since 1960, the Arc Jet Complex at NASA Ames Research Center has been a source of long-duration, high-enthalpy flow for materials testing with application to the thermal protection of aerospace vehicle components. From their inception the facilities have played an integral role supporting many of NASA's space flight programs and numerous DoD projects. In recent years advancements in computational fluid dynamics (CFO) have made the resultant models a valuable tool for assessing and predicting performance, however, the inherent limitation of models to compensate for a dissociated, transitionally high temperature nonequilibrated flowfield have made further measurements necessary. The use of electrostatic probe diagnostics within similarly harsh plasma environments in previous studies have been met with much success. In this study, the use of a single Langmuir probe was implemented to characterize the plasma parameters of interest as they vary radially within a large volume of the plume. Classical Langmuir probe theory was applied to achieve first order estimates of the heavy particle temperature, the ratio T(sub i)/T(sub e), and the ionization fraction. As expected, both the electron temperature and electron density measurements show a dependence on radial distance from the plume centerline, with electron density profiles showing the largest dependence. This paper aims to validate and strengthen modeling work conducted at the NASA Ames Research Center by measuring the critical plasma parameters (electron temperature, electron density, and plasma potential) within the arc plume of a subscale arc jet. These parameters are intended to give physical insight into the flow characteristics while providing the necessary boundary conditions to validate full scale simulations.

Codron, Douglas↗
NASA NTRS1974

The electromagnetic interchange mode in a partially ionized collisional plasma

A collisional electromagnetic dispersion relation is derived from two-fluid theory for the interchange mode coupled to the Alfven, acoustic, drift and entropy modes in a partially ionized plasma. The fundamental electromagnetic nature of the interchange model is noted; coupling to the intermediate Alfven mode is strongly stabilizing for finite k sub z. Both ion viscous and ion-neutral stabilization are included, and it was found that collisions destroy the ion finite Larmor radius cutoff at short perpendicular wavelengths.

Hudson, M. K.↗
NASA NTRS1990

A mathematical model of the structure and evolution of small-scale discrete auroral arcs

A three-dimensional fluid model for the structure and evolution of small-scale discrete auroral arcs originating from Alfven waves is developed and used to study the nonlinear macroscopic plasma dynamics of these auroral arcs. The results of simulations show that stationary auroral arcs can be unstable to a collisionless tearing mode which may be responsible for the observed transverse structuring in the form of folds and curls. At late times, the plasma becomes turbulent having transverse electric field power spectra that tend toward a universal k exp -5/3 spectral form.

Seyler, Charles E.↗
NASA NTRS1990

A new kinetic model for time-dependent polar plasma outflow - Initial results

A new time-dependent kinetic plasma outflow model has been developed, which uses a kinetic description of the parallel motions of the ion guiding centers, while assuming the electrons are a massless neutralizing fluid. The ions, O(+) and H(+) are followed as individual particles which respond to the gravitational, magnetic mirror and ambipolar electric forces as they move in one dimension along a magnetic flux tube. Results are presented for a case where the electron temperature in the flux tube is raised from a value near the ion temperature (3000 K) to a value of 10,000 K.

Wilson, G. R.↗
NASA NTRS1995

Early time evolution of a chemically produced electron depletion

The early time evolution of an ionospheric electron depletion produced by a radially expanding electron attachment chemical release is studied with a two-dimensional simulation model. The model includes electron attachment chemistry, incorporates fluid electrons, particle ions and neutrals, and considers the evolution in a plane perpendicular to the geomagnetic field for a low beta plasma. Timescales considered are of the order of or less than the cyclotron period of the negative ions that result as a by-product of the electron attacment reaction. This corresponds to time periods of tenths of seconds during recent experiemts. Simulation results show that a highly sheared azimuthal electron flow velocity develops in the radially expanding depletion boundary. This sheared electron flow velocity and the steep density gradients in the boundary give rise to small-scale irregulatities in the form of electron density cavities and spikes. The nonlinear evolution of these irregularities results in trapping and ultimately turbulent heating of the negative ions.

Scales, W. A.↗
NASA NTRS2017

A First-Principles Model Based on Saturation of the Electron Cyclotron Drift Instability for Electron Transport in Hydrodynamics Simulations of Hall Thruster Plasmas

We present a computational model based on the hypothesis that the higher-than-expected electron mobility across magnetic field lines observed in Hall thrusters is due to the growth and eventual saturation of the electron cyclotron drift instability. The key improvement of this model with respect to previous work is that it decouples the saturation of the instability from the correlation that exists between mobility and magnitude of the wave perturbations in linear theory, thereby enabling first-principles simulations. The model has been incorporated in the 2-D (r-z) multi-fluids code Hall2De and simulations of the H6 thruster at 300 V, 20 A are numerically stable and achieve a steady-state solution at a computational cost that is not significantly higher than previous Hall2De simulations. The computed location of the acceleration region is found to be within 10% of the length of the acceleration channel compared to that inferred by experiments. The simulations also capture well the plasma gradients along the channel centerline of this thruster. To further establish the validity of this model, we plan to conduct simulations of other thrusters and/or operating conditions.

Lopez Ortega, Alejandro↗
NASA NTRS2009

Efficacy of Compression Garments to Simulate Fluid Shifts during Lunar Bedrest

The effectiveness of JOBST(Registered TradeMark) compression stockings for fluid redistribution was examined in the lunar bed rest analog based on Digital Astronaut computer model predictions. NASA's future goals include missions to the moon which will require extended exposure to the lunar gravity environment. To prepare for these missions, physiological adaptations in various systems must be resolved. A bed rest at a 10 degree head-up tilt was used to simulate lunar gravity. Plasma volume (PV) loss may cause some of the cardiovascular adaptations which occur during space flight so PV is measured in the bed rest analog to evaluate fluid loss. PV varies from individual to individual so PV index (PVI) is used to determine the magnitude and time course of fluid shift and cardiovascular adaptation to 1/6 g. The Digital Astronaut, a computer simulation tool, predicts a 6% PV loss during an extended simulated lunar mission for a male with a body surface area of 1.95 m(exp 2). Simple geometry calculations suggest that 10 degree head-up tilt is most useful for simulated measurement of deconditioning in bone and muscle, however, 2 degree head-up tilt may best imitate cardiovascular fluid shifts. In order to reconcile these different models, compression stocking must be used in the 10 degree paradigm to better approximate expected cardiovascular changes.

Brinley, Alaina↗
NASA NTRS2002

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗
NASA NTRS2000

Laser Lightcraft Performance

The purpose of this study is to establish the technical ground for modeling the physics of laser powered pulse detonation phenomenon. The principle of the laser power propulsion is that when high-powered laser is focused at a small area near the surface of a thruster, the intense energy causes the electrical breakdown of the working fluid (e.g. air) and forming high speed plasma (known as the inverse Bremsstrahlung, IB, effect). The intense heat and high pressure created in the plasma consequently causes the surrounding to heat up and expand until the thrust producing shock waves are formed. This complex process of gas ionization, increase in radiation absorption and the forming of plasma and shock waves will be investigated in the development of the present numerical model. In the first phase of this study, laser light focusing, radiation absorption and shock wave propagation over the entire pulsed cycle are modeled. The model geometry and test conditions of known benchmark experiments such as those in Myrabo's experiment will be employed in the numerical model validation simulations. The calculated performance data will be compared to the test data.

Chen, Yen-Sen↗