Engineering Papers⌕ Search

Engineering topics

Mikellides, Ioannis G.

Publications and source records attributed to Mikellides, Ioannis G..

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Hybrid-PIC Computer Simulation of the Plasma and Erosion Processes in Hall Thrusters

HPHall software simulates and tracks the time-dependent evolution of the plasma and erosion processes in the discharge chamber and near-field plume of Hall thrusters. HPHall is an axisymmetric solver that employs a hybrid fluid/particle-in-cell (Hybrid-PIC) numerical approach. HPHall, originally developed by MIT in 1998, was upgraded to HPHall-2 by the Polytechnic University of Madrid in 2006. The Jet Propulsion Laboratory has continued the development of HPHall-2 through upgrades to the physical models employed in the code, and the addition of entirely new ones. Primary among these are the inclusion of a three-region electron mobility model that more accurately depicts the cross-field electron transport, and the development of an erosion sub-model that allows for the tracking of the erosion of the discharge chamber wall. The code is being developed to provide NASA science missions with a predictive tool of Hall thruster performance and lifetime that can be used to validate Hall thrusters for missions.

Hofer, Richard R.↗
NASA NTRS2010

Magnetic Shielding of the Acceleration Channel Walls in a Long-Life Hall Thruster

In a Qualification Life Test (QLT) of the BPT-4000 Hall thruster that recently accumulated greater than 10,000 h it was found that the erosion of the acceleration channel practically stopped after approximately 5,600 h. Numerical simulations of this thruster using a 2-D axisymmetric, magnetic field-aligned-mesh (MFAM) plasma solver reveal that the process that led to this significant reduction of the erosion was multifaceted. It is found that when the channel receded from its early-in-life geometry to its steady-state configuration several changes in the near-wall plasma and sheath were induced by the magnetic field that, collectively, constituted an effective shielding of the walls from any significant ion bombardment. Because all such changes in the behavior of the ionized gas near the eroding surfaces were caused by the topology of the magnetic field there, we term this process "magnetic shielding."

life modeling↗
NASA NTRS2010

Barium Depletion in the NSTAR Discharge Cathode After 30,000 Hours of Operation

Dispenser hollow cathodes rely on a consumable supply of barium released by impregnant materials in the pores of a tungsten matrix to maintain a low work function surface. Examinations of cathode inserts from long duration ion engine tests show deposits of tungsten at the downstream end that appear to block the flow of barium from the interior. In addition, a numerical model of barium transport in the insert plasma indicates that the barium partial pressure in the insert may exceed the equilibrium vapor pressure of the dominant barium-producing reaction, and it was postulated previously that this would suppress barium loss in the upstream part of the insert. New measurements of the depth of barium depletion from a cathode insert operated for 30,352 hours reveal that barium loss is confined to a narrow region near the downstream end, confirming this hypothesis.

plasma transport↗
NASA NTRS2009

Hall-Effect Thruster Simulations with 2-D Electron Transport and Hydrodynamic Ions

A computational approach that has been used extensively in the last two decades for Hall thruster simulations is to solve a diffusion equation and energy conservation law for the electrons in a direction that is perpendicular to the magnetic field, and use discrete-particle methods for the heavy species. This "hybrid" approach has allowed for the capture of bulk plasma phenomena inside these thrusters within reasonable computational times. Regions of the thruster with complex magnetic field arrangements (such as those near eroded walls and magnets) and/or reduced Hall parameter (such as those near the anode and the cathode plume) challenge the validity of the quasi-one-dimensional assumption for the electrons. This paper reports on the development of a computer code that solves numerically the 2-D axisymmetric vector form of Ohm's law, with no assumptions regarding the rate of electron transport in the parallel and perpendicular directions. The numerical challenges related to the large disparity of the transport coefficients in the two directions are met by solving the equations in a computational mesh that is aligned with the magnetic field. The fully-2D approach allows for a large physical domain that extends more than five times the thruster channel length in the axial direction, and encompasses the cathode boundary. Ions are treated as an isothermal, cold (relative to the electrons) fluid, accounting for charge-exchange and multiple-ionization collisions in the momentum equations. A first series of simulations of two Hall thrusters, namely the BPT-4000 and a 6-kW laboratory thruster, quantifies the significance of ion diffusion in the anode region and the importance of the extended physical domain on studies related to the impact of the transport coefficients on the electron flow field.

plasma↗
NASA NTRS2009

Neutralizer Hollow Cathode Simulations and Comparisons with Ground Test Data

The fidelity of electric propulsion physics-based models depends largely on the validity of their predictions over a range of operating conditions and geometries. In general, increased complexity of the physics requires more extensive comparisons with laboratory data to identify the region(s) that lie outside the validity of the model assumptions and to quantify the uncertainties within its range of application. This paper presents numerical simulations of neutralizer hollow cathodes at various operating conditions and orifice sizes. The simulations were performed using a two-dimensional axisymmetric model that solves numerically a relatively extensive system of conservation laws for the partially ionized gas in these devices. A summary of the comparisons between simulation results and Langmuir probe measurements is provided. The model has also been employed to provide insight into recent ground test observations of the neutralizer cathode in NEXT. It is found that a likely cause of the observed keeper voltage drop is cathode orifice erosion. However, due to the small magnitude of this change, is approx. 0.5 V (less than 5% of the beginning-of-life value) over 10 khrs, and in light of the large uncertainties of the cathode material sputtering yield at low ion energies, other causes cannot be excluded. Preliminary simulations to understand transition to plume mode suggest that in the range of 3-5 sccm the existing 2-D model reproduces fairly well the rise of the keeper voltage in the NEXT neutralizer as observed in the laboratory. At lower flow rates the simulation produces oscillations in the keeper current and voltage that require prohibitively small time-steps to resolve with the existing algorithms.

electric propulsion↗
NASA NTRS2009

Neutralizer Hollow Cathode Simulations and Validation with Experiments

The fidelity of electric propulsion physics-based models depends largely on the validity of their predictions over a range of operating conditions and geometries. In general, increased complexity of the physics requires more extensive comparisons with laboratory data to identify the region(s) that lie outside the validity of the model assumptions and to quantify the uncertainties within its range of application. This paper presents numerical simulations of neutralizer hollow cathodes at various operating conditions and orifice sizes. The simulations were performed using a two-dimensional axisymmetric model that solves numerically a relatively extensive system of conservation laws for the partially-ionized gas in these devices. The results for the plasma are compared directly with Langmuir probe measurements. The computed keeper voltages are also compared with the observed values. Wherever model inputs and/or specific physics of the cathode discharge are uncertain, additional sensitivity calculations have been performed to quantify the uncertainties.

partially ionized gases↗
NASA NTRS2009

Barium Depletion in Hollow Cathode Emitters

The effect of tungsten erosion, transport and redeposition on the operation of dispenser hollow cathodes was investigated in detailed examinations of the discharge cathode inserts from an 8200 hour and a 30,352 hour ion engine wear test. Erosion and subsequent re-deposition of tungsten in the electron emission zone at the downstream end of the insert reduces the porosity of the tungsten matrix, preventing the ow of barium from the interior. This inhibits the interfacial reactions of the barium-calcium-aluminate impregnant with the tungsten in the pores. A numerical model of barium transport in the internal xenon discharge plasma shows that the barium required to reduce the work function in the emission zone can be supplied from upstream through the gas phase. Barium that flows out of the pores of the tungsten insert is rapidly ionized in the xenon discharge and pushed back to the emitter surface by the electric field and drag from the xenon ion flow. This barium ion flux is sufficient to maintain a barium surface coverage at the downstream end greater than 0.6, even if local barium production at that point is inhibited by tungsten deposits. The model also shows that the neutral barium pressure exceeds the equilibrium vapor pressure of the impregnant decomposition reaction over much of the insert length, so the reactions are suppressed. Only a small region upstream of the zone blocked by tungsten deposits is active and supplies the required barium. These results indicate that hollow cathode failure models based on barium depletion rates in vacuum dispenser cathodes are very conservative.

plasma discharges↗
NASA NTRS2008

Tungsten and Barium Transport in the Internal Plasma of Hollow Cathodes

The effect of tungsten erosion, transport and redeposition on the operation of dispenser hollow cathodes was investigated in detailed examinations of the discharge cathode inserts from an 8200 hour and a 30,352 hour ion engine wear test. Erosion and subsequent re-deposition of tungsten in the electron emission zone at the downstream end of the insert reduces the porosity of the tungsten matrix, preventing the flow of barium from the interior. This inhibits the interfacial reactions of the barium-calcium-aluminate impregnant with the tungsten in the pores. A numerical model of barium transport in the internal xenon discharge plasma shows that the barium required to reduce the work function in the emission zone can be supplied from upstream through the gas phase. Barium that flows out of the pores of the tungsten insert is rapidly ionized in the xenon discharge and pushedback to the emitter surface by the electric field and drag from the xenon ion flow. Thisbarium ion flux is sufficient to maintain a barium surface coverage at the downstream endgreater than 0.6, even if local barium production at that point is inhibited by tungsten deposits. The model also shows that the neutral barium pressure exceeds the equilibrium vapor pressure of the impregnant decomposition reaction over much of the insert length,so the reactions are suppressed. Only a small region upstream of the zone blocked by tungsten deposits is active and supplies the required barium. These results indicate that hollowcathode failure models based on barium depletion rates in vacuum dispenser cathodes are very conservative.

Hollow cathodes↗
NASA NTRS2008

Assessments of Hollow Cathode Wear in the Xenon Ion Propulsion System (XIPs(c)) by Numerical Analyses and Wear Tests

The standard approach presently followed by NASA to qualify electric propulsion for the required mission throughput has been based largely on life tests, which can be costly and time consuming. Revised electric propulsion lifequalification approaches are being formulated that combine analytical and/or computational methods with (shorter-duration) wear tests. As a model case, a wear test is being performed at JPL to assess the lifetime of the discharge hollow cathode in the Xenon Ion Propulsion System (XIPS(c)), a 25-cm ion engine developed by L-3 Communications Electron Technologies, Inc. for commercial applications. Wear and plasma data accumulated throughout this life-assessment program are being used to validate the existing 2-D hollow cathode code OrCa2D. We find that the OrCa2D steady-state solution predicts very well the time-averaged plasma data and the keeper voltage after 5500 hrs of operation in high-power mode. When the wave motion that occurs naturally in these devices is accounted for, based on an estimate of the maximum wave amplitude, the molybdenum-keeper erosion profile observed in the XIPS(c) discharge cathode is also reproduced within a factor of two of the observation. When the same model is applied to predict the erosion of a tantalum keeper we find that erosion is reduced by more than two orders of magnitude compared to the molybdenum keeper due the significantly lower sputtering yield of tantalum. A tantalum keeper would therefore allow keeper lifetimes that greatly exceed the present requirements for deep-space robotic missions considered by NASA. Moreover, such large reduction of the erosion renders the largest uncertainties in the models, which are associated with the wave amplitude estimates and the electron transport model, negligible.

plasma simulations↗
NASA NTRS2008

Qualification of Commercial XIPS(R) Ion Thrusters for NASA Deep Space Missions

Electric propulsion systems based on commercial ion and Hall thrusters have the potential for significantly reducing the cost and schedule-risk of Ion Propulsion Systems (IPS) for deep space missions. The large fleet of geosynchronous communication satellites that use solar electric propulsion (SEP), which will approach 40 satellites by year-end, demonstrates the significant level of technical maturity and spaceflight heritage achieved by the commercial IPS systems. A program to delta-qualify XIPS(R) ion thrusters for deep space missions is underway at JPL. This program includes modeling of the thruster grid and cathode life, environmental testing of a 25-centimeter electromagnetic (EM) thruster over DAWN-like vibe and temperature profiles, and wear testing of the thruster cathodes to demonstrate the life and benchmark the model results. This paper will present the delta-qualification status of the XIPS thruster and discuss the life and reliability with respect to known failure mechanisms.

electric propulsion↗
NASA NTRS2008

Wear Mechanisms in Electron Sources for Ion Propulsion, 2: Discharge Hollow Cathode

The wear of the keeper electrode in discharge hollow cathodes is a major impediment to the implementation of ion propulsion onboard long-duration space science missions. The development of a predictive theoretical model for hollow cathode keeper life has long been sought, but its realization has been hindered by the complexities associated with the physics of the partially ionized gas and the associated erosion mechanisms in these devices. Thus, although several wear mechanisms have been hypothesized, a quantitative explanation of life test erosion profiles has remained incomplete. A two-dimensional model of the partially ionized gas in a discharge cathode has been developed and applied to understand the mechanisms that drove the erosion of the keeper in two long-duration life tests of a 30-cm ion thruster. An extensive set of comparisons between predictions by the numerical simulations and measurements of the plasma properties and of the erosion patterns is presented. It is found that the near-plume plasma oscillations, predicted by theory and observed by experiment, effectively enhance the resistivity of the plasma as well as the energy of ions striking the keeper.

hollow cathode↗
NASA NTRS2008

Wear Mechanisms in Electron Sources for Ion Propulsion, 1: Neutralizer Hollow Cathode

Upon the completion of two long-duration life tests of a 30-cm ion engine, the orifice channel of the neutralizer hollow cathode was eroded away to as much as twice its original diameter. Whereas the neutralizer cathode orifice opened significantly, no noticeable erosion of the discharge cathode orifice was observed. Noquantitative explanation of these erosion trends has been established since the completion of the two life tests. A two-dimensional model of the partially ionized gas inside these devices has been developed and applied to the neutralizer hollow cathode. The numerical simulations show that the main mechanism responsible for the channel erosion is sputtering by Xe+. These ions are accelerated by the sheath along the channel and bombard the surface with kinetic energy/charge of about 17 V at the beginning of cathode life. The density of the ions inside the neutralizer orifice is computed to be as high as 2.1 x 10(sup 22) m(sup -3). Because of the 3.5-times larger diameter of the discharge cathode orifice, the ion density inside the orifice is more than 40 times lower and the sheath drop 7 V lower compared with the values in the neutralizer. At these conditions, Xe+ can cause no significant sputtering of the surface.

ionized gas↗
NASA NTRS2007

Potential Fluctuations and Energetic Ion Production in Hollow Cathode Discharges

Ions with energies significantly in excess of the applied discharge voltage have been reported for many years in hollow cathode discharges. Models of dc potential hills downstream of the cathode and instabilities in postulated double layers in the cathode orifice have been proposed to explain this, but have not been substantiated. Measurements of the dc and rf plasma density and potential profiles near the exit of hollow cathodes by miniature fast-scanning probes suggests that turbulent ion acoustic fluctuations and ionization instabilities in the cathode plume significantly increase the energy of the ions that flow from this region. Increases in the discharge current and/or decreases in the cathode gas flow enhance the amplitude of the fluctuations and increase the number and energy of the energetic ions, which increases the erosion rate of the cathode electrodes. The transition from the quiescent 'spot mode' to the noisy 'plume mode' characteristic of these discharges is found to be a gradual transition of increasing fluctuation amplitudes.

keeper↗
NASA NTRS2006

Towards the Identification of the Keeper Erosion Cause(s): Numerical Simulations of the Plasma and Neutral Gas Using the Global Cathode Model OrCa2D-II

Numerical simulations with the time-dependent Orificed Cathode (OrCa2D-II) computer code show that classical enhancements of the plasma resistivity can not account for the elevated electron temperatures and steep plasma potential gradients measured in the plume of a 25-27.5 A discharge hollow cathode. The cathode, which employs a 0.11-in diameter orifice, was operated at 5.5 sccm without an applied magnetic field using two different anode geometries. It is found that anomalous resistivity based on electron-driven instabilities improves the comparison between theory and experiment. It is also estimated that other effects such as the Hall-effect from the self-induced magnetic field, not presently included in OrCa2D-II, may contribute to the constriction of the current density streamlines thus explaining the higher plasma densities observed along the centerline.

plasma modeling↗
NASA NTRS2005

Ion thruster life models

In this paper we review the most recent 2-D and 3-D plasma simulation models for ion thruster hollow cathodes, discharge chambers, and ion optics.

grid life↗
NASA NTRS2004

Model of a Hollow Cathode Insert Plasma

A 2-D axisymmetric fluid model of the plasma in the insert region of a hollow cathode is presented. The level of sophistication included in the model is motivated in part by the need to determine quantitatively plasma fluxes to the emitter surface. The ultimate goal is to assess whether plasma effects can degrade the life of impregnated inserts beyond those documented throughout the 30-50 year history of vacuum cathode technologies. Results from simulations of a 1.2-cm diameter cathode operating at a discharge current of 25 A, and a gas flow rate of 5 sccm, suggest that approximately 10 A of electron current, and 3.5 A of ion current return to the emitter surface. The total emitted electron current computed by the model is about 35 A. Comparisons with plasma measurements suggest that anomalous heating of the plasma due to two-stream instabilities is possible near the orifice region. Solution to the heavy species energy equation, with classical transport and no viscous effects, predicts heavy species temperatures as high as 2640 K.

hollow cathodes↗