Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 217 records · Page 12

High voltage spheres in an unmagnetized plasma - Fluid and PIC simulations

The basic physics involving the interaction of a high voltage sphere with a LEO-Plasma was investigated via computer experiments. Two approaches were used: the fluid description and the more general particle pushing method. Both techniques displayed qualitatively similar features. Thus, the initial time-dependent response of the plasma predicted by the fluid model, including the initial current surge and sheath formation, the formation of ion and electron density shells about the sphere with a double layer potential structure, and the subsequent propagation of the density shells away from the sphere. The PIC results also displayed an oscillatory character and slightly different time constants for the growth and decay mechanisms involved.

Thiemann, H.↗

High amplitude waves in the expanding solar wind plasma

We simulated the 1-D nonlinear time-evolution of high-amplitude Alfven, slow and fast magnetoacoustic waves in the solar wind propagating outward at different angles to the mean magnetic (spiral) field, using the expanding box model. The simulation results for Alfven waves and fast magnetoacustic waves fit the observational constraints in the solar wind best, showing decreasing trends for energies and other rms-quantities due to expansion and the appearance of inward propagating waves as minor species in the wind. Inward propagating waves are generated by reflection of Alfven waves propagating at large angles to the magnetic field or they coincide with the occurrence of compressible fluctuations. In our simulations, fast and slow magnetoacoustic waves seem to have a level in the density-fluctuations which is too high when we compare with the observations. Furthermore, the evolution of energies for slow magnetoacoustic waves differs strongly from the evolution of fluctuation energies in situ.

Schmidt, J. M.↗

Plasma Processing of Lunar Regolith Simulant for Diverse Applications

Versatile manufacturing technologies for extracting resources from the moon are needed to support future space missions. Of particular interest is the production of gases and metals from lunar resources for life support, propulsion, and in-space fabrication. Deposits made from lunar regolith could yield highly emissive coatings and near-net shaped parts for replacement or repair of critical components. Equally important is development of high fidelity lunar simulants for ground based validation of potential lunar surface operations. Described herein is an innovative plasma processing technique for insitu production of gases, metals, coatings, and deposits from lunar regolith, and synthesis of high fidelity lunar simulant from NASA issued lunar simulant JSC-1. Initial plasma reduction trials of JSC-1 lunar simulant have indicated production of metallic iron and magnesium. Evolution of carbon monoxide has been detected subsequent to reduction of the simulant using the plasma process. Plasma processing of the simulant has also resulted in glassy phases resembling the volcanic glass and agglutinates found in lunar regolith. Complete and partial glassy phase deposits have been obtained by varying the plasma process variables. Experimental techniques, product characterization, and process gas analysis will be discussed.

Schofield, Elizabeth C.↗

Sheath effects observed on a 10 meter high voltage panel in simulated low earth orbit plasma

A large (1m x 10m) flat surface of conductive material was biased to high voltage (+ or - 3000 V) to simulate the behavior of a large solar array in low earth orbit. The model array was operated in a plasma environment of 1,000 to 1,000,000/cu cm, with sufficient free space around it for the resulting plasma sheaths to develop unimpeded for 5-10 meters into the surrounding plasma. Measurements of the resulting sheath thickness were obtained. The observed thickness varied approximately as V to the 3/4 power and N to the 1/2 power. This effect appears to limit total current leakage from the test array until sheath dimensions exceed about 1 meter. Total leakage current was also measured with the array.

Mccox, J. E.↗

The Plasma Wake Downstream of Lunar Topographic Obstacles: Preliminary Results from 2D Particle Simulations

Anticipating the plasma and electrical environments in permanently shadowed regions (PSRs) of the moon is critical in understanding local processes of space weathering, surface charging, surface chemistry, volatile production and trapping, exo-ion sputtering, and charged dust transport. In the present study, we have employed the open-source XOOPIC code [I] to investigate the effects of solar wind conditions and plasma-surface interactions on the electrical environment in PSRs through fully two-dimensional pattic1e-in-cell simulations. By direct analogy with current understanding of the global lunar wake (e.g., references) deep, near-terminator, shadowed craters are expected to produce plasma "mini-wakes" just leeward of the crater wall. The present results (e.g., Figure I) are in agreement with previous claims that hot electrons rush into the crater void ahead of the heavier ions, fanning a negative cloud of charge. Charge separation along the initial plasma-vacuum interface gives rise to an ambipolar electric field that subsequently accelerates ions into the void. However, the situation is complicated by the presence of the dynamic lunar surface, which develops an electric potential in response to local plasma currents (e.g., Figure Ia). In some regimes, wake structure is clearly affected by the presence of the charged crater floor as it seeks to achieve current balance (i.e. zero net current to the surface).

Zimmerman, Michael I.↗

Simulation of radial expansion of an electron beam injected into a background plasma

A 2-D electrostatic particle code was used to study the beam radial expansion of a nonrelativistic electron beam injected from an isolated equipotential conductor into a background plasma. The simulations indicate that the beam radius is generally proportional to the beam electron gyroradius when the conductor is charged to a large potential. The simulations also suggest that the charge buildup at the beam stagnation point causes the beam radial expansion. From a survey of the simulation results, it is found that the ratio of the beam radius to the beam electron gyroradius increases with the square root of beam density and decreases inversely with beam injection velocity. This dependence is explained in terms of the ratio of the beam electron Debye length to the ambient electron Debye length. These results are most applicable to the SEPAC electron beam injection experiments from Spacelab 1, where high charging potential was observed.

Koga, J.↗

Simulating Sources of Superstorm Plasmas

We evaluated the contributions to magnetospheric pressure (ring current) of the solar wind, polar wind, auroral wind, and plasmaspheric wind, with the surprising result that the main phase pressure is dominated by plasmaspheric protons. We used global simulation fields from the LFM single fluid ideal MHD model. We embedded the Comprehensive Ring Current Model within it, driven by the LFM transpolar potential, and supplied with plasmas at its boundary including solar wind protons, polar wind protons, auroral wind O+, and plasmaspheric protons. We included auroral outflows and acceleration driven by the LFM ionospheric boundary condition, including parallel ion acceleration driven by upward currents. Our plasmasphere model runs within the CRCM and is driven by it. Ionospheric sources were treated using our Global Ion Kinetics code based on full equations of motion. This treatment neglects inertial loading and pressure exerted by the ionospheric plasmas, and will be superceded by multifluid simulations that include those effects. However, these simulations provide new insights into the respective role of ionospheric sources in storm-time magnetospheric dynamics.

Fok, Mei-Ching↗

Integration Assessment of Visiting Vehicle Induced Electrical Charging of the International Space Station Structure

The International Space Station (ISS) undergoes electrical charging in low Earth orbit (LEO) due to positively biased, exposed conductors on solar arrays that collect electrical charges from the space plasma. Exposed solar array conductors predominately collect negatively charged electrons and thus drive the metal ISS structure electrical ground to a negative floating potential (FP) relative to plasma. This FP is variable in location and time as a result of local ionospheric conditions. ISS motion through Earth s magnetic field creates an addition inductive voltage up to 20 positive and negative volts across ISS structure depending on its attitude and location in orbit. ISS Visiting Vehicles (VVs), such as the planned Orion crew exploration vehicle, contribute to the ISS plasma charging processes. Upon physical contact with ISS, the current collection properties of VVs combine with ISS. This is an ISS integration concern as FP must be controlled to minimize arcing of ISS surfaces and ensure proper management of extra vehicular activity crewman shock hazards. This report is an assessment of ISS induced charging from docked Orion vehicles employing negatively grounded, 130 volt class, UltraFlex (ATK Space Systems) solar arrays. To assess plasma electron current collection characteristics, Orion solar cell test coupons were constructed and subjected to plasma chamber current collection measurements. During these tests, coupon solar cells were biased between 0 and 120 V while immersed in a simulated LEO plasma. Tests were performed using several different simulated LEO plasma densities and temperatures. These data and associated theoretical scaling of plasma properties, were combined in a numerical model which was integrated into the Boeing Plasma Interaction Model. It was found that the solar array design for Orion will not affect the ISS FP by more than about 2 V during worst case charging conditions. This assessment also motivated a trade study to determine acceptable plasma electron current levels that can be collected by a single or combined fleet of ISS-docked VVs.

Kramer, Leonard↗

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↗

Magnetogasdynamic compression of a coaxial plasma accelerator flow for micrometeoroid simulation.

A new configuration of a coaxial plasma accelerator with self-energized magnetic compressor coil attached is described. It is shown that the circuit may be treated theoretically by analyzing an equivalent circuit mesh. The results obtained from the theoretical analysis are shown to compare favorably with the results measured experimentally. Using this accelerator configuration, glass beads of 125-micron diameter were accelerated to velocities as high as 11 km/sec, while 700-micron-diam glass beads were accelerated to velocities as high as 5 km/sec. The velocities are within the hypervelocity regime of meteroids.

Igenbergs, E. B.↗

Magnetogasdynamic compression of a coaxial plasma accelerator flow for micrometeoroid simulation

A new configuration of a coaxial plasma accelerator with self-energized magnetic compressor coil attached is described. It is shown that the circuit may be treated theoretically by analyzing an equivalent circuit mesh. The results obtained from the theoretical analysis compare favorably with the results measured experimentally. Using this accelerator configuration, glass beads of 125 micron diameter were accelerated to velocities as high as 11 kilometers per second, while 700 micron diameter glass beads were accelerated to velocities as high as 5 kilometers per second. The velocities are within the hypervelocity regime of meteoroids.

Igenbergs, E. B.↗

Titan's organic chemistry: Results of simulation experiments

Recent low pressure continuous low plasma discharge simulations of the auroral electron driven organic chemistry in Titan's mesosphere are reviewed. These simulations yielded results in good accord with Voyager observations of gas phase organic species. Optical constants of the brownish solid tholins produced in similar experiments are in good accord with Voyager observations of the Titan haze. Titan tholins are rich in prebiotic organic constituents; the Huygens entry probe may shed light on some of the processes that led to the origin of life on Earth.

Sagan, Carl↗

Numerical simulations of a nonequilibrium argon plasma in a shock-tube experiment

A code developed for the numerical modeling of nonequilibrium radiative plasmas is applied to the simulation of the propagation of strong ionizing shock waves in argon gas. The simulations attempt to reproduce a series of shock-tube experiments which will be used to validate the numerical models and procedures. The ability to perform unsteady simulations makes it possible to observe some fluctuations in the shock propagation, coupled to the kinetic processes. A coupling mechanism by pressure waves, reminiscent of oscillation mechanisms observed in detonation waves, is described. The effect of upper atomic levels is also briefly discussed.

Cambier, Jean-Luc↗

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.↗

Numerical simulation of a conducting disk partially covered by an insulator interacting with a plasma

A numerical simulation has been made of the 'pinhole effect' which produces the characteristic 'S-shaped' current-voltage curve. A disk-shaped conducting probe immersed in a plasma is modeled using a particle-in-cell (PIC) code. A probe partially covered by a very thin insulating layer is considered, as well as a probe mounted on an insulating disk. The simulation uses a cylindrical particle mover and allows for a variable number of particles in the system. The simulation space grid uses three different mesh sizes, the coarsest being away from the probe and the finest near the probe, in order to accurately calculate the trajectories of the simulation particles contributing the current to the probe and to the surface charge density on the dielectric. The calculation of the electrostatic potential is done self-consistently using successive over-relaxation (SOR). Backscattering and secondary electron emission are included for the case of positive probe voltage.

Chaky, R. C.↗

3-D Simulations of NSTAR Ion Thruster Plasma Environment

Full three-dimensional numerical models based on partical-in-cell with Monte Carlo collisions (PIC-MCC) simulations are developed to study ion thruster plasma interactions. 3-D simulation results for the NSTAR ion thruster plasma environment are presented.

NSTAR↗

Ground-Based Simulation of Low-Earth Orbit Plasma Conditions: Plasma Generation and Characterization

A 16-cm diameter plasma source operated on argon is described that is capable of producing a plasma environment that closely simulates the low Earth orbit (LEO) conditions experienced by satellites in the altitude range between 300 to 500 km. The plasma source uses a transverse-field magnetic filter, and has been successful in producing low electron temperature plasmas that contain streaming ion populations. Both of these characteristics are important because the plasma in LEO is relatively cold (e.g., Te approx. 0.1 eV) and the ram energy of the ions due to the motion of the satellite relative to the LEO plasma is high (e.g., 7,800 m/s which corresponds to approx. 5 eV for O+ ions). Plasma source operational conditions of flow rate and discharge power are presented that allow the electron temperature to be adjusted over a range from 0.14 to 0.4 eV. The expanding plasma flow field downstream of the source contains both low-energy, charge-exchange ions and streaming ions with energies that are adjustable over a range from 4 eV to 6 eV. At low flow rates and low facility pressures, the streaming ion component of the ion population comprises over 90% of the total plasma density. In the work described herein, a large area retarding potential analyzer was used to measure both electron and ion energy distribution functions in the low density, expanding plasma produced downstream of the plasma source. The benefits of using this type of plasma diagnostic tool in easily perturbed, low-density plasma are identified, and techniques are also discussed that can be used to perform real-time measurements of electron temperature. Finally, recommendations are made that may enable lower electron temperatures to be produced while simultaneously decreasing the plasma source flow rate below 1 to 2 sccm.

Williams, John D.↗