Search NASA⌕ Search

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.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Numerical Modeling and Testing of an Inductively-Driven and High-Energy Pulsed Plasma Thrusters

Pulsed Plasma Thrusters (PPTs) are advanced electric space propulsion devices that are characterized by simplicity and robustness. They suffer, however, from low thrust efficiencies. This summer, two approaches to improve the thrust efficiency of PPTs will be investigated through both numerical modeling and experimental testing. The first approach, an inductively-driven PPT, uses a double-ignition circuit to fire two PPTs in succession. This effectively changes the PPTs configuration from an LRC circuit to an LR circuit. The LR circuit is expected to provide better impedance matching and improving the efficiency of the energy transfer to the plasma. An added benefit of the LR circuit is an exponential decay of the current, whereas a traditional PPT s under damped LRC circuit experiences the characteristic "ringing" of its current. The exponential decay may provide improved lifetime and sustained electromagnetic acceleration. The second approach, a high-energy PPT, is a traditional PPT with a variable size capacitor bank. This PPT will be simulated and tested at energy levels between 100 and 450 joules in order to investigate the relationship between efficiency and energy level. Arbitrary Coordinate Hydromagnetic (MACH2) code is used. The MACH2 code, designed by the Center for Plasma Theory and Computation at the Air Force Research Laboratory, has been used to gain insight into a variety of plasma problems, including electric plasma thrusters. The goals for this summer include numerical predictions of performance for both the inductively-driven PPT and high-energy PFT, experimental validation of the numerical models, and numerical optimization of the designs. These goals will be met through numerical and experimental investigation of the PPTs current waveforms, mass loss (or ablation), and impulse bit characteristics.

Parma, Brian↗
NASA NTRS1982

Potentials on large spacecraft in LEO

A system-oriented computer code is used to predict surface charging due to voltages generated within a satellite operating in the typical dense plasma environment of LEO. The use of this code is demonstrated by predicting the expansion of electric fields onto a kapton surface from a pinhole over a biased conductor in a LEO environment. The results are compared to a more-exact solution and experimental data.

Mandell, M. J.↗
NASA NTRS1980

Strongly turbulent stabilization of electron beam-plasma interactions

The stabilization of electron beam interactions due to strongly turbulent nonlinearities is studied analytically and numerically for a wide range of plasma parameters. A fluid mode coupling code is described in which the effects of electron and ion Landau damping and linear growth due to the energetic electron beam are included in a phenomenological manner. Stabilization of the instability is found to occur when the amplitudes of the unstable modes exceed the threshold of the oscillating two-stream instability. The coordinate space structure of the turbulent spectrum which results clearly shows that soliton-like structures are formed by this process. Phenomenological models of both the initial stabilization and the asymptotic states are developed. Scaling laws between the beam-plasma growth rate and the fluctuations in the fields and plasma density are found in both cases, and shown to be in good agreement with the results of the simulation.

Freund, H. P.↗
NASA NTRS1992

Numerical simulations of the electrodynamic interactions between the Tethered-Satellite-System and space plasma

The first Tethered-Satellite-System (TSS-1), scheduled for a flight in late 1992, is expected to provide relevant information related to the concept of generating an emf in a 20-km-long (or longer) conducting wire. This paper presents numerical simulations of the electrodynamic interactions between the TSS system and space plasma, using a 2D and 3D models of the system. The 2D case code simulates the motion of a long cylinder past a plasma, which is composed of electrons and H(+) ions. The system is solved by allowing the plasma to flow past the cylinder with an imposed magnetic field. The more complex 3D case is considered to study the dynamics in great detail. Results of 2D simulation show that the interaction of a satellite with plasma flowing perpendicularly to the magnetic field results in an enhancement in the current collection.

Vashi, Bharat I.↗
NASA NTRS1977

NASCAP user's manual

The NASCAP (NASA Charging Analyzer Program) code simulates the charging process for a complex object in either tenuous plasma or ground test environment. Detailed specifications needed to run the code are presented. The object definition section, OBJDEF, allows the test object to be easily defined in the cubic mesh. The test object is composed of conducting sections which may be wholly or partially covered with thin dielectric coatings. The potential section, POTENT, obtains the electrostatic potential in the space surrounding the object. It uses the conjugate gradient method to solve the finite element formulation of Poisson's equation. The CHARGE section of NASCAP treats charge redistribution among the surface cells of the object as well as charging through radiation bombardment. NASCAP has facilities for extensive graphical output, including several types of object display plots, potential contour plots, space charge density contour plots, current density plots, and particle trajectory plots.

Mandell, M. J.↗
NASA NTRS1982

Theory and observation of a dynamically evolving negative ion plasma

The study described here examines the full range of negative-ion chemistry in the upper ionosphere by using current reaction rate data to investigate the many chemical paths SF6-type injections might take in an F region environment. Special attention is given to the conditions required to create heavy negative ions that persist long enough to affect the dynamical properties of the F region. The ambipolar diffusion characteristics of a three-component plasma are described, and estimates are given of the incoherent scatter spectra obtained from such a plasma. Model calculations using a first-order chemical code are defined and tested in order to investigate the actual types of negative-ion plasma capable of being created under nighttime conditions.

Mendillo, M.↗
NASA NTRSDate not supplied

Development of a Comprehensive Physics-Based Model for Study of NASA Gateway Lunar Dust Contamination

NASA is committed to landing the first woman and first person of color on the Moon by 2025 to begin a sustaining presence. A major component of this mission is NASA’s orbiting lunar outpost, Gateway, that will be subjected to the harsh environment of space during its 7-day orbit around the Moon. One aspect of this environment that is not well quantified is microscopic lunar regolith particles, or simply, lunar dust. As Apollo 17 Mission Commander Eugene Cernan said, “I think dust is probably one of our greatest inhibitors to a nominal operation on the moon. I think we can overcome other physiological or physical or mechanical problems except dust” [1]. These dust particles, most of which are smaller than the width of two human hairs and of highly irregular shape, are composed mainly of dielectric materials. Such characteristics can allow the dust particles to collect electric charge from the surrounding environment, resulting in electrostatic interactions between the dust, electromagnetic fields, and electrically charged surfaces [2-3]. Due to the nature of these fields, dust can collect on and contaminate charged surfaces. Lunar dust introduced into the Gateway environment by a lunar ascent vehicle element returning from the surface of the Moon presents risks to Gateway hardware such as radiators, solar arrays, antennae, and docking mechanisms. To quantify this risk and inform NASA, International Partner, and Commercial Provider stakeholders, Booz Allen, in partnership with NASA, is developing a toolset to model the interaction of charged lunar dust particles with the cis-lunar and deep-space environments. This includes a comprehensive physics model of the space plasma and solar radiation environment and electromagnetic interaction of spacecraft and lunar dust. Spacecraft charging and plasma environment are solved using the open-source code from The French Aerospace Lab (ONERA) and the European Space Agency (ESA), known as Spacecraft Plasma Interaction Software (SPIS) [4]. The physics of particle charging and motion is handled by a time-dependent implementation of Orbital-Motion-Limited (OML) theory which accounts for plasma flows and positive potential grains [5-7], using Siemens STAR-CCM+ as the framework. Model validation includes lab experiments by NASA experts and academia, as well as future on-orbit dust detecting payloads on the exterior of Gateway and Commercial Lunar Payload Services (CLPS) missions to the lunar surface [7-9]. The results and analyses will inform Gateway Program system owners at risk for lunar dust contamination.

Complex Plasma↗
NASA NTRS1982

Numerical simulation of the interaction of the plasma sheet with the lobes of the Earth's magnetotail

Codes involving one and two spatial dimensions and three velocity dimensions were used to model the Earth's magnetotail. It was shown that the magnetotail can become inflated as a consequence of low energy plasma convection toward the neutral plane. The computer study exhibits a conversion of both magnetic field energy and of energy supplied by the convection electric field into particle energy. The numerical simulations suggest that much of the magnetotail substorm morphology may be a simple consequence of an increase, followed by a decrease, in the convection electric field, without the requirement of any magnetospheric size scale plasma instability or other disruptive processes. It is also concluded that the presence of the convection electric field and a continuing replenishment of low energy particles in the magnetotail are both necessary for maintenance of the magnetotail.

Swift, D. W.↗
NASA NTRS2004

Validation of NASCAP-2K Spacecraft-Environment Interactions Calculations

The recently released Nascap-2k, version 2.0, three-dimensional computer code models interactions between spacecraft surfaces and low-earth-orbit, geosynchronous, auroral, and interplanetary plasma environments. It replaces the earlier three-dimensional spacecraft interactions codes NASCAP/GEO, NASCAP/LEO, POLAR, and DynaPAC. Nascap-2k has improved numeric techniques, a modern user interface, and a simple, interactive satellite surface definition module (Object ToolKit). We establish the accuracy of Nascap-2k both by comparing computed currents and potentials with analytic results and by comparing Nascap-2k results with published calculations using the earlier codes. Nascap-2k predicts Langmuir-Blodgett or Parker-Murphy current collection for a nearly spherical (100 surfaces) satellite in a short Debye length plasma depending on the absence or presence of a magnetic field. A low fidelity (in geometry and time) Nascap-2k geosynchronous charging calculation gives the same results as the corresponding low fidelity NASCAP/GEO calculation. A high fidelity calculation (using the Nascap-2k improved geometry and time stepping capabilities) gives higher potentials, which are more consistent with typical observations. Nascap-2k predicts the same current as a function of applied potential as was observed and calculated by NASCAP/LEO for the SPEAR I rocket with a bipolar sheath. A Nascap-2k DMSP charging calculation gives results similar to those obtained using POLAR and consistent with observation.

Davis, V. A.↗
NASA NTRS1984

Jovian decametric arcs - An estimate of the required wave normal angles from three-dimensional ray tracing.

Three-dimensional ray tracing is applied to an analysis of variable radio wave normal angle effects on the associated decametric (DAM) arc structures in the Jovian magnetosphere. The analysis is bed on 1-40 MHz radio signature recorded during Voyage 1 and 2 passages. The frequencies considered are above the R-X cut-off, and several ratios of the emission frequency/source frequency. The ray tracing code is based on a cold plasma formula and integration of the Hasselgrove (1955) equations. It is assumed that the emission is in the R-X mode, the source lies at the foot of an Io flux tube, and the emission cone is hollow. Attention is focused on data for two intense, vertex-late, high curvature DAM arc. A possible source for the arcs is found to be doppler-shifted gyroemission from a beam of electrons with an energy of 10 keV. A value of 1.1 is set as the limit of the doppler shift of the DAM emissions.

Menietti, J. D.↗
NASA NTRS2001

Fluid Aspects of Solar Wind Disturbances Driven by Coronal Mass Ejections

Transient disturbances in the solar wind initiated by coronal eruptions have been modeled for many years, beginning with the self-similar analytical models of Parker and Simon and Axford. The first numerical computer code (one-dimensional, gas dynamic) to study disturbance propagation in the solar wind was developed in the late 1960s, and a variety of other codes ranging from simple one-dimensional gas dynamic codes through three-dimensional gas dynamic and magnetohydrodynamic codes have been developed in subsequent years. For the most part, these codes have been applied to the problem of disturbances driven by fast CMEs propagating into a structureless solar wind. Pizzo provided an excellent summary of the level of understanding achieved from such simulation studies through about 1984, and other reviews have subsequently become available. More recently, some attention has been focused on disturbances generated by slow CMEs, on disturbances driven by CMEs having high internal pressures, and disturbance propagation effects associated with a structured ambient solar wind. Our purpose here is to provide a brief tutorial on fluid aspects of solar wind disturbances derived from numerical gas dynamic simulations. For the most part we illustrate disturbance evolution by propagating idealized perturbations, mimicking different types of CMEs, into a structureless solar wind using a simple one-dimensional, adiabatic (except at shocks), gas dynamic code. The simulations begin outside the critical point where the solar wind becomes supersonic and thus do not address questions of how the CMEs themselves are initiated. Limited to one dimension (the radial direction), the simulation code predicts too strong an interaction between newly ejected solar material and the ambient wind because it neglects azimuthal and meridional motions of the plasma that help relieve pressure stresses. Moreover, the code ignores magnetic forces and thus also underestimates the speed with which pressure disturbances propagate in the wind.

Gosling, J. T.↗
NASA NTRS2013

Numerical Simulations of the Partially-Ionized Gas in a 100-A LaB6 Hollow Cathode

Numerical simulations of a hollow cathode with a LaB6 emitter operating at 100 A have been performed for the first time using the 2-D Orificed Cathode (OrCa2D) code. Results for a variety of plasma properties are presented and compared with laboratory measurements. The large size of the device permits peak electron number densities in the cathode interior that are lower than those established in the NSTAR hollow cathode, which operates with a 7.3x lower discharge current and 3.2x lower mass flow rate. Also, despite the higher discharge current in the LaB6 cathode, the maximum electron current density is lower, by 4.2x, than that in the NSTAR cathode due to the larger orifice size. Simulations and direct measurements show that at 12 sccm of xenon flow the peak emitter temperature is in the range of 1594-1630 C. It is also found that the conditions for the excitement of current-driven streaming instabilities and ion-acoustic turbulence (IAT) are satisfied in this cathode, similarly to what was found in the past in its smaller counterparts like the NSTAR cathode. Based on numerical simulations, it has long been argued that these instabilities may be responsible for the anomalously large ion energies that have been measured in these discharges as well as for the enhancement of the plasma resistivity. Confirmation of the presence of IAT in this cathode is presented for the first time in a companion paper.

robotic exploration↗
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 NTRS1991

Gyro-phase effects near the storm-time boundary of energetic plasma

The nonadiabatic acceleration of plasma sheet ions during the expansion phase of substorms is examined by means of single-particle codes. It is shown that, in the near-earth plasma sheet, the gyration phase at substorm onset controls the net ion energization, as the gyro-period is locally comparable to the field variation time scale. This can yield a particularly significant decrease of the particle magnetic moment. It is accordingly argued that, via adiabatic invariant violation, the dipolarization of magnetospheric field lines can induce a 'de-trapping' of inner plasma sheet populations and, hence, give rise to short-lived precipitation enhancements over the auroral zone. As this effect depends upon cyclotron frequency, it is expected to occur at distinct latitudes for different ion species. This contrasts with nonadiabatic ion behavior further out into the geotail, which is characterized by intense perpendicular heating and a likely collective trapping.

Delcourt, D. C.↗
NASA NTRS1989

The soft x ray telescope for Solar-A

The Solar-A satellite being prepared by the Institute for Sapce and Astronautical Sciences (ISAS) in Japan is dedicated to high energy observations of solar flares. The Soft X Ray Telescope (SXT) is being prepared to provide filtered images in the 2 to 60 A interval. The flight model is now undergoing tests in the 1000 foot tunnel at MSFC. Launch will be in September 1991. Earlier resolution and efficiency tests on the grazing incidence mirror have established its performance in soft x rays. The one-piece, two mirror grazing incidence telescope is supported in a strain free mount separated from the focal plane assembly by a carbon-epoxy metering tube whose windings and filler are chosen to minimize thermal and hygroscopic effects. The CCD detector images both the x ray and the concentric visible light aspect telescope. Optical filters provide images at 4308 and 4700 A. The SXT will be capable of producing over 8000 of the smallest partial frame images per day, or fewer but larger images, up to 1024 x 1024 pixel images. Image sequence with two or more of the five x ray analysis filters, with automatic exposure compensation to optimize the charge collection by the CCD detector, will be used to provide plasma diagnostics. Calculations using a differential emission measure code were used to optimize filter selection over the range of emission measure variations and to avoid redundancy, but the filters were chosen primarily to give ratios that are monotonic in plasma temperature.

Brown, W. A.↗
NASA NTRS1994

Three-dimensional numerical simulation of current collection by a probe in a magnetized plasma

A three-dimensional numerical model for current collection in a magnetized plasma is reported. The model is based on an electrostatic particle-in-cell code. The model yields self-consistent sheath structure including distributions of plasma and the electric potential around the body and the collection of electrons. The analytical theory of current collection by a body in a magnetized plasma yields an upper bound for the collected current determined by the conservation of energy and canonical angular momentum. The theory shows that the collected charged particles come from a cylindrical volume aligned with the magnetic shadow of the body; the maximum radius r(sub o) of this volume is determined by the body size, body potential, and the ambient magnetic field. This theory does not deal with the sheath structure around the body. The condition for the actual current to approach the upper-bound value has been a matter of debate. Our simulations reveal when and why the collected current becomes equal to its upper-bound value. Sheath size in the radial direction perpendicular to the axial ambient magnetic field is determined by the current-limiting radius r(sub o). Our simulation yields time-average current in good agreement with its upper bound. This feature of the current collection is explained as follows: Once electrons enter the sheath, some of them are freely accelerated perpendicular to the magnetic field because they are demagnetized by the large gradients in the perpendicular electric fields. Simulations show a large perpendicular acceleration, producing perpendicular energy as large as that determined by the potential on the body, especially in the region where perpendicular electric fields (E perpendicular) are the strongest. An analysis shows that the demagnetization of electrons occurs above a threshold potential on the body. This threshold condition follows from the breakdown of the adiabaticity of the electron dynamics inside the sheath.

Singh, Nagendra↗
NASA NTRS1992

A more accurate nonequilibrium air radiation code - NEQAIR second generation

Two experiments, one an equilibrium flow in a plasma torch at Stanford, the other a nonequilibrium flow in a SDIO/IST Bow-Shock-Ultra-Violet missile flight, have provided the basis for modifying, enhancing, and testing the well-known radiation code, NEQAIR. The original code, herein termed NEQAIR1, lacked computational efficiency, accurate data for some species and the flexibility to handle a variety of species. The modified code, herein termed NEQAIR2, incorporates recent findings in the spectroscopic and radiation models. It can handle any number of species and radiative bands in a gas whose thermodynamic state can be described by up to four temperatures. It provides a new capability of computing very fine spectra in a reasonable CPU time, while including transport phenomena along the line of sight and the characteristics of instruments that were used in the measurements. Such a new tool should allow more accurate testing and diagnosis of the different physical models used in numerical simulations of radiating, low density, high energy flows.

Moreau, Stephane↗