Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma 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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Magnetic field reconnection in a collisionless plasma

A reasonably consistent model of steady-state magnetic-field-line reconnection in a collisionless plasma is constructed by incorporating ion-acoustic anomalous resistance into the hydromagnetic flow in the vicinity of the x-type neutral line. The Petschek-Vasyliunas (1975) reconnection theory is applied, and properties of the ion-acoustic instability are reviewed for the case of comparable ion and electron temperatures. Nonlinear saturation of the instability is examined, the saturation wave intensity is determined as a function of electron drift speed and electron/ion temperature ratio, and the computed wave intensities are used to estimate the steady electric field in the neutral region. Ion-acoustic anomalous resistance is shown to limit the electron drift speed to slightly above the marginally stable value. A model for the resistive-diffusion region is constructed which incorporates the properties of ion-acoustic anomalous resistance, and an approximate solution for the external flow region is matched to the resistive-region solution. It is found that the two solutions are sensibly matched only for a restricted range of upstream plasma parameters. Limitations and possible extensions of the model are discussed.

Coroniti, F. V.↗
NASA NTRS1972

Nonlinear interaction of resonant plasma oscillations.

Using the water-bag model of an electron plasma the nonlinear interaction of two plasma wave pulses (each of sufficiently narrow wave-number spectrum) is discussed. It is shown that for the special case of a resonant (zero group velocity) pulse interacting with a nonresonant pulse, the interaction produces an effective slowing down of the group velocity of the nonresonant pulse. An application of this theory to the explanation of the apparent resonances at fractions of the plasma frequency observed in the ionosphere is discussed.

Lewak, G. J.↗
NASA NTRS1998

Flow Characteristics of UHV Nitrogen Arcjet

Flow characteristics in the nozzle of a nitrogen arcjet have been simulated by a model. Electrons, ions, and neutral atoms and molecules are represented as separate species. Dissociation, ionization, and recombination are treated as separate reactions. Thermal non-equilibrium is represented by a two-temperature model. Energy input to plasma from the arc is modeled as a source determined by local current and electrical conductivity. Momentum and energy losses and recombinations at the nozzle wall are included as source/sink determined by fluxes to the wall. The results show essentially frozen gas-phase reactions and thermal nonequilibrium due to Joule heating of electrons and in the expansion part of the nozzle. Wall interactions have strong effects on the results, indicating that they play important roles in the flow due to very high surface area relative to the volume of the plasma. For example, plasma velocity decreases from 6000 m per s without wall interactions to 3000 m per s with modest amount of wall interactions. These results also show reasonable agreement with optical emission measurements, which confirms that the arcjet plasma is far from LTE. The spectra suggest nitrogen dissociation levels of 0.3% - 9%, depending on nitrogen flow rate and arc plasma current, which also determine the relative amounts of excited atomic and molecular nitrogen. Langmuir probe studies of the source show that electron and ion fluxes increase with increasing power, and that the ion energy distribution shifts to lower energies. Typical ion fluxes were on the order of 4E-9 A/cm**2 with a maximum ion kinetic energy of 3.5eV. The median electron energy was 1eV, with a maximum of less than 4eV.

Chang, C. H.↗
NASA NTRS1995

Cooling of solar flares plasmas. 1: Theoretical considerations

Theoretical models of the cooling of flare plasma are reexamined. By assuming that the cooling occurs in two separate phase where conduction and radiation, respectively, dominate, a simple analytic formula for the cooling time of a flare plasma is derived. Unlike earlier order-of-magnitude scalings, this result accounts for the effect of the evolution of the loop plasma parameters on the cooling time. When the conductive cooling leads to an 'evaporation' of chromospheric material, the cooling time scales L(exp 5/6)/p(exp 1/6), where the coronal phase (defined as the time maximum temperature). When the conductive cooling is static, the cooling time scales as L(exp 3/4)n(exp 1/4). In deriving these results, use was made of an important scaling law (T proportional to n(exp 2)) during the radiative cooling phase that was forst noted in one-dimensional hydrodynamic numerical simulations (Serio et al. 1991; Jakimiec et al. 1992). Our own simulations show that this result is restricted to approximately the radiative loss function of Rosner, Tucker, & Vaiana (1978). for different radiative loss functions, other scaling result, with T and n scaling almost linearly when the radiative loss falls off as T(exp -2). It is shown that these scaling laws are part of a class of analytic solutions developed by Antiocos (1980).

Cargill, Peter J.↗
NASA NTRS2007

Numerical Simulations of Flow Separation Control in Low-Pressure Turbines using Plasma Actuators

A recently introduced phenomenological model to simulate flow control applications using plasma actuators has been further developed and improved in order to expand its use to complicated actuator geometries. The new modeling approach eliminates the requirement of an empirical charge density distribution shape by using the embedded electrode as a source for the charge density. The resulting model is validated against a flat plate experiment with quiescent environment. The modeling approach incorporates the effect of the plasma actuators on the external flow into Navier Stokes computations as a body force vector which is obtained as a product of the net charge density and the electric field. The model solves the Maxwell equation to obtain the electric field due to the applied AC voltage at the electrodes and an additional equation for the charge density distribution representing the plasma density. The new modeling approach solves the charge density equation in the computational domain assuming the embedded electrode as a source therefore automatically generating a charge density distribution on the surface exposed to the flow similar to that observed in the experiments without explicitly specifying an empirical distribution. The model is validated against a flat plate experiment with quiescent environment.

Suzen, Y. B.↗
NASA NTRS2022

Hybrid Kinetic Model of the Interaction Between the Dense Plasma Clouds and Magnetospheric Plasma on Large Time and Spatial Scales, and Comparison With MMS Observations

We present a new simulation results of the cloud dynamics in the ambient magnetospheric plasma on the large time and spatial scales. It was assumed that these impulsive structures observed by the MMS spacecraft originally were created because of the reconnection at the magnetopause. Our new 3-D hybrid kinetic modeling on the large time and spatial scales captures several of these processes: an excitation of the electromagnetic waves (whistler and shear-Alfvén waves) and plasma instabilities (mirror and flute); a formation of shock waves, and collapsing diamagnetic cavity; particle acceleration. A strong overshoot in plasma density profile was observed in the modeling and MMS observation at the interface between the cloud and magnetospheric plasma. The cloud expansion into ambient magnetospheric plasma causes the flute waves connected with excitation of the Rayleigh-Taylor instability observed at the overshoot in plasma density profile across the external magnetic field. The modeling demonstrates a formation of the whistler waves at the initial stage which propagate in the external magnetic field direction. At the later stage, a formation of shear-Alfvén waves was observed

A. S. Lipatov↗
NASA NTRS1994

Evaluation Of Mathematical Models Of VPPA Welding

Report describes study of three mathematical models of variable-polarity plasma arc welding VPPAW. Models based on physics of welding process considered at three different levels of abstraction.

Mccutcheon, Kimble D.↗
NASA NTRS1994

A model of Triton's role in Neptune's magnetosphere

Escape of neutral hydrogen (H) and nitrogen (N) from Triton's maintains a large neutral cloud, called the Triton torus, in Neptune's magnetosphere. We have developed the first detailed Monte Carlo simulation model of the Triton torus that includes the collisionality, the complex geometry, the injection of two neutral species from Triton (H and N), and the combined effects of photoionization, electron impact ionization, and charge exchange. Ionization in Neptune's plasma sheet was modeled using Voyager plasma observations. Collisions cause both the H and N neutral clouds to become more radially extended, both toward Neptune and out beyond the magnetopause, as well as more extended in latitude, when compared with collisionless models. Moreover, collisions of H with the much more massive N greatly enhance the collisional ejection of H from the system and into Neptune's atmosphere. This effect decreases the probability of H ionization within the magnetosphere relative to that for N, and furthermore causes model results for two-species injection from Triton to differ significantly from those for H injection alone. For a hydrogen escape rate from Triton of 5 x 10(exp 25)/s, as given by photo-chemical models of Triton's upper atmosphere, a nitrogen escape rate of 5 x 10(exp 24)/s gives proton and N(+) sources of 5.6 x 10(exp 24)/s and 3.3 x 10(exp 24)/s, respectively, whose ratio is close to the observed ratio of protons to heavies. A nitrogen escape rate of 2 x 10(exp 25)/s, yields an N(+) source more than twice that of protons, inconsistent with the Voyager data.

Decker, R. B.↗
NASA NTRS1994

Centrifugal acceleration of the polar wind

The effect of parallel ion acceleration associated with convection was first applied to energization of test particle polar ions by Cladis (1986). However, this effect is typically neglected in 'self-consistent' models of polar plasma outflow, apart from the fluid simulation by Swift (1990). Here we include approximations for this acceleration, which we broadly characterize as centrifugal in nature, in our time-dependent, semikinetic model of polar plasma outflow and describe the effects on the bulk parameter profiles and distribution functions of H+ and O+. For meridional convection across the pole the approximate parallel force along a polar magnetic field line may be written as F(sub cent, pole) = 1.5m(E(sub i))/B(sub i))squared (r(squared)/r(sup 3)(sub i)) where m is ion mass, r is geometric distance; and E(sub i), B(sub i) and r(sub i) refer to the electric and magnetic field magnitudes and geocentric distance at the ionosphere, respectively. For purely longitudinal convection along a constant L shell the parallel force is F(cent. long) = F(sub cent, pole)(1 - (r/(r(sub i)L))(sup 3/2)/(1 - 3r/(4 r(sub i)L))(sup 5/2). For high latitudes the difference between these two cases is relatively unimportant below approximately 5 R(sub E). We find that the steady state O+ bulk velocities and parallel temperatures strongly increase and decrease, respectively, with convection strength. In particular, the bulk velocities increase from near 0 km/s at 4000 km altitude to approximately 10 km/s at 5 R(sub E) geocentric distance for 50-mV/m ionospheric convection electric field. However, the centrifugal effect on the steady O+ density profiles depends on the exobase ion and electron temperatures: for low-base temperatures (T(sub i) = T(sub e) = 3000 K) the O+ density at high altitudes increases greatly with convection, while for higher base temperatures (T(sub i) = 5000 K, T(sub e) = 9000 K), the high-altitude O+ density decreases somewhat as convection is enhanced. The centrifugal force further has a pronounced effect on the escaping O+ flux, especially for cool exobase conditions; as referenced to the 4000-km altitude, the steady state O+ flux increases from 10(exp 5) ions/sq cm/s when the ionospheric convection field E(sub i) = 0 mV/m to approximately 10(exp 7) ions/sq cm/s when E(sub i) = 100 mV/m. The centrifugal effect also decreases the time scale for approach to steady-state. For example, in the plasma expansion for T(sub i) = T(sub e) = 3000 K, the O+ density at 7 R(sub E) reaches only 10(exp -7) of it final value approximately 1.5 hours after expansion onset for E(sub i) = 0. For meridional convection driven by E(sub i) = 50 mV/m, the density at the same time after initial injection is 30-50% of its asymptotic level. The centrifugal acceleration described here is a possible explanation for the large (up to approximately 10 km/s or more) o+ outflow velocities observed in the midlatitude polar magnetosphere with the Dynamics Explorer 1 and Akebono spacecraft.

Horwitz, J. L.↗
NASA NTRS2011

Computational Validation of a Two-Dimensional Semi-Empirical Model for Inductive Coupling in a Conical Pulsed Inductive Plasma Thruster

A two-dimensional semi-empirical model of pulsed inductive thrust efficiency is developed to predict the effect of such a geometry on thrust efficiency. The model includes electromagnetic and gas-dynamic forces but excludes energy conversion from radial motion to axial motion, with the intention of characterizing thrust efficiency loss mechanisms that result from a conical versus a at inductive coil geometry. The range of conical pulsed inductive thruster geometries to which this model can be applied is explored with the use of finite element analysis. A semi-empirical relation for inductance as a function of current sheet radial and axial position is the limiting feature of the model, restricting the applicability as a function of half cone angle to a range from ten degrees to about 60 degrees. The model is nondimensionalized, yielding a set of dimensionless performance scaling parameters. Results of the model indicate that radial current sheet motion changes the axial dynamic impedance parameter at which thrust efficiency is maximized. This shift indicates that when radial current sheet motion is permitted in the model longer characteristic circuit timescales are more efficient, which can be attributed to a lower current sheet axial velocity as the plasma more rapidly decouples from the coil through radial motion. Thrust efficiency is shown to increase monotonically for decreasing values of the radial dynamic impedance parameter. This trend indicates that to maximize the radial decoupling timescale should be long compared to the characteristic circuit timescale.

Hallock, Ashley K.↗
NASA NTRS1996

Kinetic Space Weather: Toward a Global Hybrid Model of the Polar Ionosphere-Lower Magnetosphere Plasma Transport

During the indicated period of performance, we had a number of publications concerned with kinetic polar ionosphere-lower magnetosphere plasma transport. For the IUGG 1991-4 Quadrennial Report, we reviewed aspects of U.S. accomplishments concerned with polar plasma transport, among other issues. In another review, we examined the computer simulations of multiple-scale processes in space plasmas, including polar plasma outflow and transport. We also examined specifically multiscale processes in ionospheric outflows. We developed a Generalized Semi-Kinetic(GSK) model for the topside-lower magnetosphere which explored the synergistic action of wave heating and electric potentials in the formation of auroral Ion conics, in particular the "pressure cooker" mechanism. We extended the GSK model all the way down to 120 km and applied this code to illustrate the response of the ionosphere- magnetosphere to soft-electron precipitation and convection-driven frictional ion heating, respectively. Later, the convection-driven heating work was extended to a paper for the Journal of Geophysical Research. In addition to the above full published papers, we also presented the first developments of the coupled fluid-semikinetic model for polar plasma transport during this period. The results from a steady-state treatment were presented, with the second presentation being concerned with the effects of photo-electrons on the polar wind, and the first garnering an outstanding student paper award from the American Geophysical Union. We presented the first results from a time-dependent version of this coupled fluid-semikinetic model.

Horwitz, James L.↗
NASA NTRS1977

Methodology of modeling and measuring computer architectures for plasma simulations

A brief introduction to plasma simulation using computers and the difficulties on currently available computers is given. Through the use of an analyzing and measuring methodology - SARA, the control flow and data flow of a particle simulation model REM2-1/2D are exemplified. After recursive refinements the total execution time may be greatly shortened and a fully parallel data flow can be obtained. From this data flow, a matched computer architecture or organization could be configured to achieve the computation bound of an application problem. A sequential type simulation model, an array/pipeline type simulation model, and a fully parallel simulation model of a code REM2-1/2D are proposed and analyzed. This methodology can be applied to other application problems which have implicitly parallel nature.

Wang, L. P. T.↗
NASA NTRS1999

Semiempirical Two-Dimensional Magnetohydrodynamic Model of the Solar Corona and Interplanetary Medium

We have developed a two-dimensional semiempirical MHD model of the solar corona and solar wind. The model uses empirically derived electron density profiles from white-light coronagraph data measured during the Skylub period and an empirically derived model of the magnetic field which is fitted to observed streamer topologies, which also come from the white-light coronagraph data The electron density model comes from that developed by Guhathakurta and coworkers. The electron density model is extended into interplanetary space by using electron densities derived from the Ulysses plasma instrument. The model also requires an estimate of the solar wind velocity as a function of heliographic latitude and radial component of the magnetic field at 1 AU, both of which can be provided by the Ulysses spacecraft. The model makes estimates as a function of radial distance and latitude of various fluid parameters of the plasma such as flow velocity V, effective temperature T(sub eff), and effective heat flux q(sub eff), which are derived from the equations of conservation of mass, momentum, and energy, respectively. The term effective indicates that wave contributions could be present. The model naturally provides the spiral pattern of the magnetic field far from the Sun and an estimate of the large-scale surface magnetic field at the Sun, which we estimate to be approx. 12 - 15 G. The magnetic field model shows that the large-scale surface magnetic field is dominated by an octupole term. The model is a steady state calculation which makes the assumption of azimuthal symmetry and solves the various conservation equations in the rotating frame of the Sun. The conservation equations are integrated along the magnetic field direction in the rotating frame of the Sun, thus providing a nearly self-consistent calculation of the fluid parameters. The model makes a minimum number of assumptions about the physics of the solar corona and solar wind and should provide a very accurate empirical description of the solar corona and solar wind Once estimates of mass density rho, flow velocity V, effective temperature T(sub eff), effective heat flux q(sub eff), and magnetic field B are computed from the model and waves are assumed unimportant, all other plasma parameters such as Mach number, Alfven speed, gyrofrequency, etc. can be derived as a function of radial distance and latitude from the Sun. The model can be used as a planning tool for such missions as Slar Probe and provide an empirical framework for theoretical models of the solar corona and solar wind The model will be used to construct a semiempirical MHD description of the steady state solar corona and solar wind using the SOHO Large Angle Spectrometric Coronagraph (LASCO) polarized brightness white-light coronagraph data, SOHO Extreme Ultraviolet Imaging Telescope data, and Ulysses plasma data.

Sittler, Edward C., Jr.↗
NASA NTRS2004

Electron Density and Two-Channel Neutron Emission Measurements in Steady-State Spherical Inertial-Electrostatically Confined Plasmas, with Review of the 1-D Kinetic Model

Electron density measurements have been made in steady-state plasmas in a spherical inertial electrostatic confinement (IEC) discharge using microwave interferometry. Plasma cores interior to two cathodes, having diameters of 15 and 23 cm, respectively, were probed over a transverse range of 10 cm with a spatial resolution of about 1.4 cm for buffer gas pressures from 0.2 to 6 Pa in argon and deuterium. The transverse profiles are generally flat, in some cases with eccentric symmetric minima, and give mean densities of from approx. = 0.4 to 7x 10(exp 10)/cu cm, the density generally increasing with the neutral gas pressure. Numerical solutions of the 1-D Poisson equation for EC plasmas are reviewed and energy distribution functions are identified which give flat transverse profiles. These functions are used with the plasma approximation to obtain solutions which also give densities consistent with the measurements, and a double potential well solution is obtained which has minima qualitatively similar to those observed. Explicit consideration is given to the compatibility of the solutions interior and exterior to the cathode, and to grid transparency. Deuterium fusion neutron emission rates were also measured and found to be isotropic, to within the measurement error, over two simultaneous directions. Anisotropy was observed in residual emissions during operation with non-fusing hydrogen-1. The deuterium rates are consistent with predictions from the model.

Dobson, Chris C.↗
NASA NTRS2001

Soluble TNF-alpha receptor 1 and IL-6 plasma levels in humans subjected to the sleep deprivation model of spaceflight

BACKGROUND: The extent to which sleep loss may predispose astronauts to a state of altered immunity during extended space travel prompts evaluation with ground-based models. OBJECTIVE: We sought to measure plasma levels of selected cytokines and their receptors, including the putative sleep-regulation proteins soluble TNF-alpha receptor (sTNF-alpha R) I and IL-6, in human subjects undergoing 2 types of sleep deprivation during environmental confinement with performance demands. METHODS: Healthy adult men (n = 42) were randomized to schedules that varied in severity of sleep loss: 4 days (88 hours) of partial sleep deprivation (PSD) involving two 2-hour naps per day or 4 days of total sleep deprivation (TSD). Plasma samples were obtained every 6 hours across 5 days and analyzed by using enzyme-linked immunoassays for sTNF-alpha RI, sTNF-alpha RII, IL-6, soluble IL-2 receptor, IL-10, and TNF-alpha. RESULTS: Interactions between the effects of time and sleep deprivation level were detected for sTNF-alpha RI and IL-6 but not for sTNF-alpha RII, soluble IL-2 receptor, IL-10, and TNF-alpha. Relative to the PSD condition, subjects in the TSD condition had elevated plasma levels of sTNF-alpha RI on day 2 (P =.04), day 3 (P =.01), and across days 2 to 4 of sleep loss (P =.01) and elevated levels of IL-6 on day 4 (P =.04). CONCLUSIONS: Total sleep loss produced significant increases in plasma levels of sTNF-alpha RI and IL-6, messengers that connect the nervous, endocrine, and immune systems. These changes appeared to reflect elevations of the homeostatic drive for sleep because they occurred in TSD but not PSD, suggesting that naps may serve as the basis for a countermeasures approach to prolonged spaceflight.

Non-NASA Center↗
NASA NTRS1991

Computer simulation of astrophysical plasmas

The role of sophisticated numerical models and simulations in the field of plasma astrophysics is discussed. The need for an iteration between microphysics and macrophysics in order for astrophysical plasma physics to produce quantitative results that can be related to astronomical data is stressed. A discussion on computational requirements for simulations of astrophysical plasmas contrasts microscopic plasma simulations with macroscopic system models. An overview of particle-in-cell simulations (PICS) is given and two examples of PICS of astrophysical plasma are discussed including particle acceleration by collisionless shocks in relativistic plasmas and magnetic field reconnection in astrophysical plasmas.

Max, Claire E.↗
NASA NTRS2018

Effects of Multiscale Phase-Mixing and Interior Conductance in the Lunar-like Pickup Ion Plasma Wake. First Results from 3-D Hybrid Kinetic Modeling

The study of multiscale pickup ion phase-mixing in the lunar plasma wake with a hybrid model is the main subject of our investigation in this paper. Photoionization and charge exchange of protons with the lunar exosphere are the ionization processes included in our model. The computational model includes the self-consistent dynamics of the light (hydrogen ions or hydrogen molecule ion or helium ion), and heavy (sodium ion ) pickup ions. The electrons are considered as a fluid. The lunar interior is considered as a weakly conducting body. In this paper we considered for the first time the cumulative effect of heavy neutrals in the lunar exosphere (e.g., Aluminum, Argon), an effect which was simulated with one species of sodium ion but with a tenfold increase in total production rates. We find that various species produce various types of plasma tail in the lunar plasma wake. Specifically, sodium ion and helium ion pickup ions form a cycloid-like tail, whereas the hydrogen ion or hydrogen molecule ion pickup ions form a tail with a high density core and saw-like periodic structures in the flank region. The length of these structures varies from 1:5 R(sub M) to 3:3 R(sub M) depending on the value of gyro radius for hydrogen ion or hydrogen molecule ion pickup ions. The light pickup ions produce more symmetrical jump in the density and magnetic field at the Mach cone which is mainly controlled by the conductivity of the interior, an effect previously unappreciated. Although other pickup ion species had little effect on the nature of the interaction of the Moon with the solar wind, the global structure of the lunar tail in these simulations appeared quite different when the hydrogen molecule ion production rate was high.

Lipatov, A.S.↗