Engineering Papers⌕ Search

Engineering topics

Mikellides, Ioannis G.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Self-Consistent Implementation of a Zero-Equation Transport Model Into a Predictive Model for a Hall Effect Thruster

The performance of an axisymmetric multi-fluid Hall thruster code that incorporates a self-consistent, data-driven closure model for the anomalous electron transport is investigated. Five different operating conditions of the H9 magnetically shielded Hall thruster are simulated with the Jet Propulsion Laboratory’s Hall2De. In order to capture the inherent uncertainty associated with the closure model, O(100) simulations are run for each condition, each of using a coefficient set sampled randomly from a probability distribution. The results of these simulations provide probabilistic predictions of thruster performance quantities including thrust, and discharge current, as well as several component efficiencies and centerline plasma properties. The model is found to yield converged solutions at all conditions, with large 10 kHzrange oscillations and performance trends with voltage and flow rate similar to experiment. The model under-predicts the thrust by 15-25% and over-predicts the discharge current by 20% on average compared to experiments at the same discharge voltage and mass flow rate. This performance discrepancy is due to lower beam utilization, mass utilization, and divergence efficiency than experiment, resulting from high Hall parameters in the acceleration region, which lead to a protracted ion acceleration region. The physical processes underlying this result are discussed in the context of future data-driven modeling efforts.

Jorns, Benjamin A.↗
NASA NTRS2020

Hall2De Simulations of a Magnetic Nozzle

The two-dimensional, fluid-based Hall thruster code Hall2De is adapted to simulate a low temperature propulsive magnetic nozzle. Electron cyclotron resonance (ECR) thrusters developed at the Office National d’Etudes et de Recherches Aérospatiales (ONERA) in France and the University of Michigan are investigated. For each device, the simulated domain is defined from the nozzle throat to a location 10 thruster radii downstream. A numerical mesh is defined based on the measured magnetic field, and the governing fluid equations for singlycharge ions, electrons, and neutrals are solved. The boundary conditions for the throat of the magnetic nozzle are based on measured values of electron temperature and ion velocity. Multiple boundary conditions for the electron temperature in the downstream regions of the domain are considered. The simulated ion velocity and potential agrees well with experiment, while the simulated divergence efficiencies and thrusts are 25% lower and 45-55% higher than experiment. When a Neumann boundary condition for electron temperature is applied to the far-plume boundaries, the nozzle was seen to be purely isothermal. This suggests that the thermal conductivity along field lines in real devices is much lower than classical simulations predict, and anomalous resistivity may be needed to reproduce the temperature distribution observed in experiment. The larger beam divergence observed in our simulation may also be related to the lack of anomalous resistivity.

Jorns, Benjamin A.↗
NASA NTRS2020

Electric Propulsion for the Psyche Mission: Development Activities and Status

NASA’s Psyche mission will launch in 2022 and begin a 3.6-year cruise to the metallic asteroid Psyche, where it will examine this unique body. The baseline spacecraft design is a hybrid of JPL’s deep-space heritage subsystems with commercial partner Maxar’s electric propulsion, power, and structure subsystems. All primary propulsion will be done with SPT-140 thrusters, which will be the first use of Hall thrusters for a NASA mission. The electric propulsion subsystem and its implementation for the Psyche mission are described here. Major testing activities have included the successful completion of subsystem integrated testing with the design modifications required for Psyche, and a series of low-power thrust repeatability tests that were performed in support of navigation analyses. Thruster performance models have been further validated with new SPT-140 flight data, and new analyses of thruster swirl torque have been performed that result in much higher values than previously estimated. Analysis of recent Maxar flight data has also provided a new understanding of in-flight propellant usage uncertainties. Subsystem integration and test activities are now underway and the status and plans are discussed.

Johnson, Ian↗
NASA NTRS2019

HERMeS Thruster Magnetic Field Topology Optimization Study: Performance, Stability, and Wear Results

NASA's Hall Effect Rocket with Magnetic Shielding (HERMeS) 12.5kW Technology Demonstration Unit-1 (TDU-1) has been the subject of extensive technology maturation in preparation for flight system development. The TDU-1 thruster implements a magnetically shielded field topology and has demonstrated the elimination of the discharge channel erosion as a life limiting mechanism. Extensive wear testing the TDU Hall thrusters has identified the thruster front pole covers as the next life limiting component. This effort aims to explore and investigate alternate magnetic field topologies to assess whether reductions in the front pole cover erosion can be attained while still maintaining very low erosion rates on the discharge channel walls. Four candidate magnetic field topologies that reduce the effectiveness of the shielding along the discharge channel walls with the intent to also reduce the erosion rates along the front pole covers were designed. Three of the four candidate magnetic field topologies (B1, B2, and B4) have been manufactured and were subjected to an extensive test campaign that includes laser induced fluorescence (LIF), performance, stability, wear, plume, and thermal characterization. In Phase I, LIF measurements along the discharge chamber centerline found that upstream retraction of the thruster's peak magnetic field does result in an upstream shift of the acceleration zone but the magnitude of the shift does not correspond one-to-one to the shift in the location of the peak radial magnetic field magnitude. Phase II test segment results found at a normalized thruster magnetic field setting of 1, the thruster performance was very similar for all configurations. Discharge current waveforms ripple data indicated that configurations Bo, B1, and B2 has similar oscillatory profiles with the B2 configuration transitioning to higher oscillatory mode at 400V instead of the 450V observed for configurations B0 and B1. Configuration B4 ripple data indicates that the thruster was operating in a very oscillatory mode above 350V. Inner front pole cover erosion rates for configuration B1 were approximately 65% relative to B0 and the erosion rates for configuration B2 were 40% relative to B0. Analysis of the performance, stability, and wear test results of configurations B0, B1, B2, and B4 indicates that configuration B2 may present an alternative option for B0 because it mostly maintained the same performance and stability as configuration B0 but with 40% of the inner front pole cover erosion rate.

Kamhawi, Hani↗
NASA NTRS2019

Numerical investigations of background pressure effects and channel erosion in the SPT-140 Hall thruster for the Psyche mission

Four SPT-140 Hall thrusters are the main components of the electric propulsion system for the Psyche mission. We employ the 2-D (r-z) code Hall2De to conduct numerical simulations of the SPT-140 with the goal of estimating erosion of the acceleration channel rings for the Psyche mission profile. The simulations are informed by direct measurements of the plasma conditions in the acceleration channel that were obtained using laser-induced fluorescence. We validate our simulation results with additional plasma measurements, wear test erosion rates, and performance data as a function of background pressure. We then make use of a series of numerical simulations to compute the progressive erosion of the channel walls, using as input parameters a profile for the Psyche mission that includes operation at 4.5, 2.5 and 1.0 kW with duration of all segments increased by a 50% margin. We find that the channel walls are not completely eroded at the end of the mission (plus 50% margin). We also observe a progressive decrease in thrust as the channel walls erode that is consistent with the wear test and flight data, with a rate of ~4.5 mN/kh for the first 2,500 h and constant thrust after that. The reasons for these thrust trends are explained.

Lenguito, Giovanni↗
NASA NTRS2019

Mechanism behind the dependence of thrust on facility backpressure and implications on the operation of the SPT-140 onboard the Psyche mission

The Psyche mission will employ Stationary Plasma Thrusters (SPT-140) for primary propulsion in space and will be the first to use Hall thrusters beyond lunar orbit. One of the longest standing risks in the implementation of this EP technology on missions like Psyche is our limited understanding of the effects of ground test facilities on the performance and wear of the thruster.

Lenguito, Giovanni↗
NASA NTRS2019

Plasma simulations for the assessment of pole erosion in the Magnetically Shielded Miniature Hall Thruster (MaSMi)

Numerical simulations have been conducted to assess the erosion rates at the pole covers of the magnetically shielded miniature (MaSMi) Hall thruster. MaSMi is part of the Ascendant Sub-kW Transcelestial Electric Propulsion System (ASTRAEUS), whose design requirements demand a xenon throughput of 100 kg, thus enabling deep space missions using SmallSats. The results of our simulations show that the proposed thickness of the graphite pole covers will not be completely eroded at end of life. We find that the erosion rates are larger near the corners of the covers, which are directly exposed to the plasma in the acceleration channel.

Chaplin, Vernon H.↗
NASA NTRS2018

Investigations of Pole Erosion Mechanisms in the 12.5 kW HERMeS Hall Thruster Using Numerical Simulations and Ion Velocity Measurements

When the Hall Effect Rocket with Magnetic Shielding (HERMeS) was tested at the NASA GRC’s VF-5 the erosion rate on the pole covers was found to be at least a factor of two higher at 300 V and 20.8 A than at other operating conditions in the 400 V to 600 V range. Simulations using our hydrodynamics code Hall2De do not predict increased erosion rates at 300 V but accurately compute erosion rates similar to the measurements at the other operating conditions. We investigate the source of the discrepancy between measurements and simulations at 300 V using a combination of numerical simulations and experimental measurements of the ion velocity fields in the acceleration region of the thruster. By examining previous simulations that predicted the measured erosion rates for other thrusters and operating conditions, we determine that sputtering of the pole surfaces by high energy ions is the most likely mechanism behind the erosion rates at 300 V. High energy ions sputter the pole surfaces when the acceleration region of the thruster is downstream of the pole surface plane. When the latter occurs, the curvature of the plasma potential contours at the edges of the channel accelerates a small fraction of the high energy ions radially. We find that neither our Hall2De simulations nor the experimental measurements produce ions of sufficiently high energy and flux to the inner pole cover to explain the measured erosion there. Only at the outer edge of the cover is where we find the ions needed to yield simulation results that are comparable to the erosion measurements. We argue that one possible mechanism for the higher erosion is high energy ions that graze the channel corner of the pole and become trapped in a sheath that develops downstream. Within the sheath, we find that the electric field is large enough to turn the ions towards the surface of the pole. A simplified calculation shows that the erosion rates produced by sputtering of trapped ions are similar to the measurements. We also propose that the presence of local oscillations and high-energy cathode ions may be sources of the enhanced erosion We conclude this article by proposing a combination of analysis, simulation, and experimental measurements that can be used to address the validity of these hypotheses.

Mikellides, Ioannis G.↗
NASA NTRS2018

Hollow Cathode Thermal Modelling and Self-Consistent Plasma Solution: Work Function Evaluation for a LaB6 Cathode.

Self-heating hollow cathodes are central components in modern electric thrusters. The plasma discharge inside these devices heats the internal components, thus maintaining the temperatures required for electron emission. Precise knowledge of the physical phenomena governing hollow cathode operation is key to predict their lifetime, specifically, their thermionic emission characteristics. A simulation platform has been built to couple plasma and thermal models of the self-heating hollow cathode to produce a self-consistent solution. A self-consistent solution has been found for a LaB6 hollow cathode operating at 25A and 13 sccm where the work function is assumed to be spatially uniform along the emitter with a value which is allowed to vary as the coupled model iterates to a self-consistent solution. The emitter temperature from the converged solution does not agree with experimental temperature measurements, however. The results of a sensitivity analysis suggest that none of the tolerances in the measurement are responsible for the discrepancy. We hypothesize that either the work function needs to be a function of position along the emitting surfaces or the heat fluxes have been overestimated in the plasma solver.

Meiron, Daniel I.↗
NASA NTRS2018

Investigations of Spot-to-plume Mode Transition in a Hollow Cathode Discharge Using 2-D Axisymmetric Plasma Simulations

After decades of laboratory investigations that provided invaluable measurements and insight, the physics behind the transition from spot to plume modes in hollow cathodes remains one of the longest standing theoretical problems in electric propulsion. This has prohibited the development of ab initio models that allow for the prediction of the transition across different cathodes and operating conditions. Since the beginning of its development over a decade ago, simulations with the 2-D axisymmetric Orificed Cathode (OrCa2D) code have helped elucidate a wide range of processes in hollow cathode discharges. However, the code has never been used to investigate the onset of plume mode. We present results from the first OrCa2D simulations of a 25-A LaB6 cathode for a range of flow rates (5-20 sccm) in which transition from spot to plume modes is known to occur. The cathode in this study was one of the two technologies considered for the 12.5 kW Hall Effect Rocket with Magnetic Shielding (HERMeS) and operates nominally at 21 A and 14.8 sccm. The simulations capture the characteristic rise of the peak-to-peak amplitude in the keeper voltage oscillations and underscore the significance of the plume neutral gas in the transition. The plasma inside the cathode is found to be relatively quiescent throughout the transition, in agreement with previous experimental observations. The computed keeper voltage fluctuations at low flow rates (<8 sccm) are found to be driven by oscillations of the same frequency in the plasma plume with the following main characteristics: (1) they are of low frequency (<10 kHz), and associated with small longitudinal motion in the direction of the applied magnetic field, (2) they occur in a region of the plume where the neutral gas provided by the cathode has been fully depleted, and (3) they have a (small) wave velocity of about 100 m/s, which is at least ~10 smaller than the drift, thermal and acoustic speeds of the ions. At 8 sccm, when the transition to the large-amplitude oscillations begins, the ionization frequency in the neutral-depleted plume region ranges ~2-100 kHz. The frequency of the oscillations in the plasma (and keeper voltage) is found to be equal to the ionization frequency (~5 kHz) at the center of this region. The findings suggest that the transition to plume mode is driven by ionization processes in the near-plume of the cathode, in line with previous conjectures that were based solely on laboratory observations.

Polk, James E.↗
NASA NTRS2018

Metallic Wall Hall Thrusters

A Hall thruster apparatus having walls constructed from a conductive material, such as graphite, and having magnetic shielding of the walls from the ionized plasma has been demonstrated to operate with nearly the same efficiency as a conventional nonmagnetically shielded design using insulators as wall components. The new design is believed to provide the potential of higher power and uniform operation over the operating life of a thruster device.

Goebel, Dan Michael↗
NASA NTRS2016

Metallic Wall Hall Thrusters

A Hall thruster apparatus having walls constructed from a conductive material, such as graphite, and having magnetic shielding of the walls from the ionized plasma has been demonstrated to operate with nearly the same efficiency as a conventional non-magnetically shielded design using insulators as wall components. The new design is believed to provide the potential of higher power and uniform operation over the operating life of a thruster device.

Goebel, Dan Michael↗
NASA NTRS2016

Facility Effect Characterization Test of NASA's HERMeS Hall Thruster

A test to characterize the effect of varying background pressure on NASA's 12.5-kW Hall Effect Rocket with Magnetic Shielding had being completed. This thruster is the baseline propulsion system for the Solar Electric Propulsion Technology Demonstration Mission (SEP TDM). Potential differences in thruster performance and oscillation characteristics when in ground facilities versus on-orbit are considered a primary risk for the propulsion system of the Asteroid Redirect Robotic Mission, which is a candidate for SEP TDM. The first primary objective of this test was to demonstrate that the tools being developed to predict the zero-background-pressure behavior of the thruster can provide self-consistent results. The second primary objective of this test was to provide data for refining a physics-based model of the thruster plume that will be used in spacecraft interaction studies. Diagnostics deployed included a thrust stand, Faraday probe, Langmuir probe, retarding potential analyzer, Wien filter spectrometer, and high-speed camera. From the data, a physics-based plume model was refined. Comparisons of empirical data to modeling results are shown.

Electric Propulsion↗
NASA NTRS2016

Hall2De Simulations with a First-Principles Electron Transport Model Based on the Electron Cyclotron Drift Instability

Several years of work have combined plasma measurements in a laboratory Hall thruster and r-z numerical simulations with Hall2De to isolate the spatial variation of the anomalous collision frequency needed in Ohm’s law to produce the observed thruster behavior. This numerical solution is used here to test the validity of a first-principles model of the anomalous transport in these devices before such model is implemented self-consistently in r-z fluid codes like Hall2De. The first-principles model employs quasi-linear theory and is based on the hypothesis that the Electron Cyclotron Drift Instability (ECDI) excites ion acoustic turbulence that, in turn, enhances the effective collision frequency in these devices. We find that a model of the ECDI that assumes Maxwellian velocity distributions for electrons and singly-charged, main-beam, cold ions (Ti=0.07 eV) is insufficient to explain the expected variation of the anomalous collision frequency both in the interior and exterior of the acceleration channel. When warm ions (~0.5-3 eV) are accounted for, the ECDI model in the channel interior appears more promising but fails by orders of magnitude in the near plume region due to the much higher Landau damping of the ion acoustic waves there. This implies that either (a) some process allows the ECDI instability to remain uninhibited by classical Landau damping or, (b) that a different instability (or instabilities) altogether, also insusceptible to Landau damping, is/are active in this region. A previous hypothesis, that convection of wave energy generated by the ECDI in the channel plays a significant role in the near plume, is not supported by the results of the simulations.

Mikellides, Ioannis G.↗