SEARCH · Search NASA

Results for “Plasma codes”

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 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1972

Monte Carlo Simulation of Nonlinear Radiation Induced Plasmas

A Monte Carlo simulation model for radiation induced plasmas with nonlinear properties due to recombination was, employing a piecewise linearized predict-correct iterative technique. Several important variance reduction techniques were developed and incorporated into the model, including an antithetic variates technique. This approach is especially efficient for plasma systems with inhomogeneous media, multidimensions, and irregular boundaries. The Monte Carlo code developed has been applied to the determination of the electron energy distribution function and related parameters for a noble gas plasma created by alpha-particle irradiation. The characteristics of the radiation induced plasma involved are given.

Wang, B. S.↗

NASA NTRS2001

Numerical Model of the Plasma Sheath Generated by the Plasma Source Instrument Aboard the Polar Satellite

The plasma sheath generated by the operation of the Plasma Source Instrument (PSI) aboard the Polar satellite is studied by using a three-dimensional particle-in-cell (PIC) code. When the satellite passes through the region of low-density plasma, the satellite charges to positive potentials as high as 40-50 V, owing to the photoelectron emission. In such a case, ambient core ions cannot accurately be measured or detected. The goal of the onboard PSI is to reduce the floating potential of the satellite to a sufficiently low value so that the ions in the polar wind become detectable. When the PSI is operated, ion-rich xenon plasma is ejected from the satellite, such that the floating potential of the satellite is reduced and is maintained at approximately 2 V. Accordingly, in our three-dimensional PIC simulation we considered that the potential of the satellite is 2 V as a fixed bias. Considering the relatively high density of the xenon plasma in the sheath (10-10(exp 3)/cc), the ambient plasma of low density (<1/cc) is neglected. In the simulations the electric fields and plasma dynamics are calculated self-consistently. We found that an 'apple'-shape positive potential sheath forms surrounding the satellite. In the region near the PSI emission a high positive potential hill develops. Near the Thermal Ion Dynamics Experiment detector away from the PSI, the potentials are sufficiently low for the ambient polar wind ions to reach it. In the simulations it takes only about a couple of tens of electron gyroperiods for the sheath to reach a quasi steady state. This time is approximately the time taken by the heavy Xe(+) ions to expand up to about one average Larmor radius of electrons from the satellite surface. After this time the expansion of the sheath in directions transverse to the ambient magnetic field slows down because the electrons are magnetized. Using the quasi steady sheath, we performed trajectory calculations to characterize the detector response to a highly supersonic polar wind flow. The detected ions' velocity distribution shows significant deviations from a shifted Maxwellian in the ambient polar wind population. The deviations are caused by the effects of electric fields on the ions' motion as they traverse the sheath.

Singh, N.↗

NASA NTRS1993

Kinetic simulation of plasma flows in the inner magnetosphere

A one-dimensional hybrid particle code is used to study the interactions between upflowing thermal ions from conjugate ionospheres. The simulation model allows for multiple species, convection of plasmaspheric flux tubes, and Coulomb self-collisions which conserve momentum and energy locally. The model incorporates a variable-flux boundary condition where the flux, at the boundaries, approaches zero as the plasmasphere fills and equilibrium conditions are reached. The effects of two important processes on plasmaspheric refilling have been considered. The first includes convection of the plasmaspheric flux tube. The second is the interaction of ionospheric thermal plasma and particle injection from an external source. Particle injection seems to play an important role in the evolution of the total particle distribution on the early timescales (t less than 1 hour); however, for late timescales (t larger than 8 days) the thermal plasma from the ionosphere dominates the particle distribution.

Miller, Ronald H.↗

NASA NTRS1975

Computer simulation of plasma and N-body problems

The following FORTRAN language computer codes are presented: (1) efficient two- and three-dimensional central force potential solvers; (2) a three-dimensional simulator of an isolated galaxy which incorporates the potential solver; (3) a two-dimensional particle-in-cell simulator of the Jeans instability in an infinite self-gravitating compressible gas; and (4) a two-dimensional particle-in-cell simulator of a rotating self-gravitating compressible gaseous system of which rectangular coordinate and superior polar coordinate versions were written.

Harries, W. L.↗

NASA NTRS1981

Additional extensions to the NASCAP computer code, volume 3

The ION computer code is designed to calculate charge exchange ion densities, electric potentials, plasma temperatures, and current densities external to a neutralized ion engine in R-Z geometry. The present version assumes the beam ion current and density to be known and specified, and the neutralizing electrons to originate from a hot-wire ring surrounding the beam orifice. The plasma is treated as being resistive, with an electron relaxation time comparable to the plasma frequency. Together with the thermal and electrical boundary conditions described below and other straightforward engine parameters, these assumptions suffice to determine the required quantities. The ION code, written in ASCII FORTRAN for UNIVAC 1100 series computers, is designed to be run interactively, although it can also be run in batch mode. The input is free-format, and the output is mainly graphical, using the machine-independent graphics developed for the NASCAP code. The executive routine calls the code's major subroutines in user-specified order, and the code allows great latitude for restart and parameter change.

Mandell, M. J.↗

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 NTRS2012

General Relativistic Simulations of Magnetized Plasmas Around Merging Supermassive Black Holes

Coalescing supermassive black hole binaries are produced by the mergers of galaxies and are the most powerful sources of gravitational waves accessible to space-based gravitational observatories. Some such mergers may occur in the presence of matter and magnetic fields and hence generate an electromagnetic counterpart. In this paper we present the first general relativistic simulations of magnetized plasma around merging supermassive black holes using the general relativistic magnetohydrodynamic code Whisky. By considering different magnetic field strengths, going from non-magnetically dominated to magnetically dominated regimes, we explore how magnetic fields affect the dynamics of the plasma and the possible emission of electromagnetic signals. In particular we observe, total amplification of the magnetic field of approx 2 orders of magnitude which is driven by the accretion onto the binary and that leads to stronger electromagnetic signals than in the force-free regime where such amplifications are not possible.

Giacomazzo, Bruno↗

NASA NTRS1997

Nonlinear Waves and Particle Heating in the Auroral Ionosphere

In the two year period, the researchers have developed innovative 2-D simulation codes for modeling the interaction between electron beams and plasma waves and particles in the auroral ionosphere. These new simulations include kinetic wave particles as well as wave-wave interactions near and above approximately 500 km. The principal numerical model consists of partial-differential equations which evolve large amplitude magnetized 2-D Langmuir waves self-consistently with the electron distribution function. Progress was also made in developing 2-D Particle in Cell (PIC) codes.

Goldman, Martin↗

NASA NTRS1986

Simulating AOTV heating environment in an arc jet

The application of arc-jet facilities to the simulation of the combined radiative and convective heating environments calculated for AOTV return missions is presented. The AOTV mission requirement and design configuration determine the magnitude of the radiative and convective heating rates. This paper considers methods for obtaining radiation from the high temperature plasma in arc-jet flows. A computer code that predicts arc-heater performance was used to determine the effect of gas mixtures on arc-heater performance, with the intent of increasing the gas temperature and resultant radiative emission. Numerical radiation solutions were compared with direct measurements of the radiative intensity from the gas cap. Convective heating measurements were made as well. The arc-jet was operated at AOTV enthalpy and pressure levels for these measurements to assess the simulation capabilities for AOTV material development/testing. A parametric study is presented that defines the capability of the Interaction Heating Facility arc jet at Ames to simulate both the radiative and the convective heating for current AOTV configurations.

Balboni, J.↗

NASA NTRS1987

Electron heating by ion acoustic turbulence in simulated low Mach number shocks

Explicit and fully electromagnetic particle-in-cell simulations of perpendicular, collisionless, and nominally subcritical shocks are performed in one and two spatial dimensions using the code wave. Shock parameters are chosen to maximixe the growth rates of the current-driven ion acoustic instability in the shock. Electron heating by ion acoustic turbulence is observed at the shocks, at rates in agreement with second-order Vlasov theory predictions. However, the amount of resistive electron heating is small and ion reflection provides the major source of dissipation. Strictly resistive shocks do not exist for the parameters suitable for explicit particle codes running on today's supercomputers, because the plasma convects through these shocks so quickly that current-driven instabilities have little time to be amplified and to heat the electrons resistively. This effect is primarily a result of the relatively small values of omega(pe)/omega(ce) that can be analyzed.

Tokar, Robert L.↗

NASA NTRS1994

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

The nature of turbulent cascades arising from the parametric instabilities of a monochromatic field-aligned large-amplitude circularly polarized Alfven wave is investigated via direct numerical simulation for the case of low plasma Beta and no wave dispersion. The magnetohydrodynamic 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. Anisotropic turbulent cascades, similar to those found in early incompressible two-dimensional simulations, occur after nonlinear saturation of the parallel propagating decay instability. The turbulent spectrum can be divided into three regimes: the lowest wave numbers are dominated by lower sideband remnants of the parametric process, intermediate wave numbers display nearly incompressible dynamics, and the highest wave numbers are dominated by acoustic turbulence.

Ghosh, S.↗

NASA NTRS2004

Computational Ion Optics Design Evaluations

Ion optics computational models are invaluable tools in the design of ion optics systems. In this study a new computational model developed by an outside vendor for use at the NASA Glenn Research Center (GRC) is presented. This computational model is a gun code that has been modified to model the plasma sheaths both upstream and downstream of the ion optics. The model handles multiple species (e.g. singly and doubly-charged ions) and includes a charge-exchange model to support erosion estimations. The model uses commercially developed solid design and meshing software to allow high flexibility in ion optics geometric configurations. The results from this computational model are applied to the NEXT project to investigate the effects of crossover impingement erosion seen during the 2000-hour wear test.

Malone, Shane P.↗

NASA NTRS2004

Computational Modeling of NEXT 2000-Hour Wear Test Results

Ion optics computational models are invaluable tools for the design of ion optics systems. In this study, a new computational model developed by an outside vendor for NASA Glenn Research Center (GRC) is presented. This model is a gun code which has been modified to model the plasma sheaths both upstream and downstream of the ion optics. The model handles multiple species (e.g. singly and doubly-charged ions) and includes a charge-exchange model for erosion estimates. The model uses commercially available solid design and meshing software, allowing high flexibility in ion optics geometric configurations. This computational model is compared to experimental results from the NASA Evolutionary Xenon Thruster (NEXT) 2000-hour wear test, including over-focusing along the edge apertures, pit-and-groove erosion due to charge exchange, and beamlet distortion at the edge of the hole pattern.

Malone, Shane P.↗

NASA NTRS2021

Fe xvii 2p–3s Line Ratio Diagnostic of Shock Formation Radius in O Stars

The 2p–3s lines of Fe xvii in the X-ray spectrum of the O-type star ζ Puppis exhibit an anomalous (3G + M2)/(3F) line ratio of ∼1.4, in comparison with ∼2.4 for almost all other collisionally excited astrophysical spectra. Based on the work of Mauche et al., we conjectured that the strong UV field of ζ Puppis produces the observed ratio by depopulation of metastable 3s excited states, and that the ratio can potentially be used as an independent diagnostic of plasma formation radius. We used the Flexible Atomic Code collisional-radiative model to model the effect of UV photoexcitation from O stars on the Fe xvii lines. We compared our model calculations to archival spectra of coronal and hot stars from the Chandra HETGS and XMM-Newton RGS to benchmark our calculations for various electron densities and UV field intensities. Our calculations show that UV photoexcitation does not produce a sufficiently large dynamic range in the 3F / (3F + 3G + M2) fraction to explain the difference in the observed ratio between coronal stars and ζ Pup. Thus, this effect likely cannot explain the observed line ratio of ζ Pup, and its origin is still unexplained.

Gabriel J. Grell↗

NASA NTRS2021

Fe XVII 2p–3s Line Ratio Diagnostic of Shock Formation Radius in O Stars

Gabriel Grell↗

NASA NTRS2004

Simulations of Solar Wind Plasma Flow Around a Simple Solar Sail

In recent years, a number of solar sail missions of various designs and sizes have been proposed (e.g., Geostorm). Of importance to these missions is the interaction between the ambient solar wind plasma environment and the sail. Assuming a typical 1 AU solar wind environment of 400 km/s velocity, 3.5 cu cm density, ion temperature of approx.10 eV, electron temperature of 40 eV, and an ambient magnetic field strength of 10(exp -4) G, a first order estimate of the plasma interaction with square solar sails on the order of the sizes being considered for a Geostorm mission (50 m x 50 m and 75 m x 75 m corresponding to approx.2 and approx.3 times the Debye length in the plasma) is carried out. First, a crude current balance for the sail surface immersed in the plasma environment and in sunlight was used to estimate the surface potential of the model sails. This gave surface potentials of approx.10 V positive relative to the solar wind plasma. A 3-D, Electrostatic Particle-in-Cell (PIC) code was then used to simulate the solar wind flowing around the solar sail. It is assumed in the code that the solar wind protons can be treated as particles while the electrons follow a Boltzmann distribution. Next, the electric field and particle trajectories are solved self-consistently to give the proton flow field, the electrostatic field around the sail, and the plasma density in 3-D. The model sail was found to be surrounded by a plasma sheath within which the potential is positive compared to the ambient plasma and followed by a separate plasma wake which is negative relative to the plasma. This structure departs dramatically from a negatively charged plate such as might be found in the Earth s ionosphere on the night side where both the plate and its negative wake are contiguous. The implications of these findings are discussed as they apply to the proposed Geostorm solar sail mission.

Garrett, Henry B.↗

NASA NTRS1994

A simulation study of radial expansion of an electron beam injected into an ionospheric plasma

Injections of nonrelativistic electron beams from a finite equipotential conductor into an ionospheric plasma have been simulated using a two-dimensional electrostatic particle code. The purpose of the study is to survey the simulation parameters for understanding the dependence of beam radius on physical variables. The conductor is charged to a high potential when the background plasma density is less than the beam density. Beam electrons attracted by the charged conductor are decelerated to zero velocity near the stagnation point, which is at a few Debye lengths from the conductor. The simulations suggest that the beam electrons at the stagnation point receive a large transverse kick and the beam expands radially thereafter. The buildup of beam electrons at the stagnation point produces a large electrostatic force responsible for the transverse kick. However, for the weak charging cases where the background plasma density is larger than the beam density, the radial expansion mechanism is different; the beam plasma instability is found to be responsible for the radial expansion. The simulations show that the electron beam radius for high spacecraft charging cases is of the order of the beam gyroradius, defined as the beam velocity divided by the gyrofrequency. In the weak charging cases, the beam radius is only a fraction of the beam gyroradius. The parameter survey indicates that the beam radius increases with beam density and decreases with magnetic field and beam velocity. The beam radius normalized by the beam gyroradius is found to scale according to the ratio of the beam electron Debye length to the ambient electron Debye length. The parameter dependence deduced would be useful for interpreting the beam radius and beam density of electron beam injection experiments conducted from rockets and the space shuttle.

Koga, J.↗

NASA NTRS1996

3-D Simulations of NSTAR Ion Thruster Plasma Interactions

Described is a Particle-in-Cell with Monte Carlo Collision code developed to perform detailed three-dimensional ion thruster simulations. To capture the full kinetic behavior of ion thruster plumes, both the electrons and ions are treated as test particles. Simulation results are given of the NSTAR ion thruster under ground test and in space conditions. Numerical results are compared.

NSTAR Ion Thruster Plasma Ionization Particles↗