Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thrusters”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A theoretical thrust density limit for Hall thrusters

Abstract Hall Thrusters typically operate at thrust densities on the order of 10 N/m $$^2$$ 2 , which appear to be orders of magnitude below the thrust density limits suggested in previous literature. These limits have been considered here and each component of thrust density is analyzed to demonstrate the relative contribution to the total thrust density. Dependencies of the thrust density limits upon the thruster geometry, electron mobility, and the applied magnetic field are revealed and compared with experimental measurements of thrust density. This analysis reveals that with conventional applied magnetic field strengths, Hall thruster thrust density appears to be on the order of 1000 N/m $$^2$$ 2 . It is shown that this limit can be further increased through higher applied magnetic fields, applied voltage, and suppression of anomalous electron transport. This suggests Hall thrusters can be made much more compact and operated at higher power densities, given improvements to the thermal management and materials.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Effect of plasma initialization on 3D PIC simulation of Hall thruster azimuthal instability

Abstract The lack of understanding of the azimuthal instability and the resulting electron anomalous transport limits further improvement of Hall thrusters. Compared to theoretical and experimental approaches, the numerical particle-in-cell (PIC) simulation is a suitable and powerful tool, which has been widely applied to investigate the azimuthal instability, and great progress has been made in the past decades. However, PIC simulations are intrinsically computationally expensive, and it is realized that the Hall thruster azimuthal instability has a three dimensional nature. Therefore, massive 3D PIC simulation must be carried out to completely reveal the mechanism of the instability. In this paper, the effect of plasma initialization on 3D PIC simulation of Hall thruster azimuthal instability is studied as a starting point. It is found that by initializing with ion density and velocity fitting functions to the steady-state simulation results, a faster convergence can be obtained and the computational time can be reduced by about 1.5 times. Typical fitting functions of ion density, drifting velocity, and temperature are given, and the influence of different initialization profiles is presented.

Physics↗
DOE OSTI2023

Multiple frequency electron cyclotron resonance thruster

An electron cyclotron resonance (ECR) thruster includes a magnetic field source configured to generate a magnetic field, a thruster body that defines a chamber, the thruster body being disposed relative to the magnetic field source such that the magnetic field is present in the chamber and such that a magnetic nozzle is established, an antenna configured to propagate radio frequency (RF) power within the chamber, and a waveform generator coupled to the antenna to generate an RF waveform for the RF power. The waveform generator is configured such that the RF waveform includes multiple frequencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

GPU-accelerated kinetic Hall thruster simulations in WarpX

Abstract Two-dimensional (axial-azimuthal) simulations of a Hall thruster are performed using the open-source particle-in-cell code WarpX. The simulation conditions are chosen to match those of the axial-azimuthal benchmark first reported by Charoy et al. in 2019. A range of numerical and solver parameters is investigated in order to find those which yield the best performance. It is found that WarpX completes the benchmark case in 3.8 days on an Nvidia V100 GPU, and in as low as 1.5 days on a more recent Nvidia H100 GPU. Of the numerical parameters investigated, it is determined that the field-solve tolerance and particle resampling thresholds have the largest effect on the simulation wall time and that particle resampling may artificially widen electron velocity distribution functions, leading to unphysical heating. A semi-implicit scheme for the electrostatic field solve is tested and is found to produce results consistent to within 10% of the benchmark in less than twelve hours. The scaling properties of the electrostatic solver to multiple GPUs are also assessed on a uniform plasma test problem. The results of this work are discussed in the context of advancements in GPU hardware and the suitability of kinetic Hall thruster simulations for engineering applications.

Marks, Thomas A.↗
DOE OSTI2026

Impact of reverse flow induced by a sawtooth anode on the performance of an argon Hall thruster

This study reports the performance comparison of the RAIJIN-66 TAL-type Hall thruster with argon propellant, using two hollow anode designs, one with a simple straight shape and another one with a sawtooth shape. Simulations of the rarefied gas in the sawtooth hollow anode region suggest a successful reversal of the direction of neutral particle flow, with an average increase in the neutral particle density by 15.5% when compared to the straight-shaped hollow anode. This contributes to the improvement in performance across a range of discharge voltages that was measured experimentally. When the thruster is operated with a discharge voltage of 150 V and a flow rate of 70 SCCM, the propellant utilization efficiency with the straight anode is 14%, while with the sawtooth anode it is 25.1%. The anode efficiency reaches in the vicinity of 15% with the sawtooth anode in high voltage conditions, exceeding the efficiencies achieved with the straight anode in the same operating conditions. The optimal magnetic field condition for argon operation and the tuning parameters of the sawtooth anode design for RAIJIN-66 are also discussed.

Satpathy, Dibyesh [Univ. of Tokyo (Japan)] (ORCID:↗
DOE OSTI2022

Fluid and hybrid simulations of the ionization instabilities in Hall thruster

Low-frequency axial oscillations in the range of 5–50 kHz stand out as a pervasive feature observed in many types of Hall thrusters. While it is widely recognized that the ionization effects play the central role in this mode, as manifested via the large-scale oscillations of neutral and plasma density, the exact mechanism(s) of the instabilities remain unclear. To gain further insight into the physics of the breathing mode and evaluate the role of kinetic effects, a one-dimensional time-dependent full nonlinear low-frequency model describing neutral atoms, ions, and electrons is developed in full fluid formulation and compared to the hybrid model in which the ions and neutrals are kinetic. Both models are quasi-neutral and share the same electron fluid equations that include the electron diffusion, mobility across the magnetic field, and the electron energy evolution. The ionization models are also similar in both approaches. Further, the predictions of fluid and hybrid simulations are compared for different test cases. Two main regimes are identified in both models: one with pure low-frequency behavior and the other one, where the low-frequency oscillations coexist with high-frequency oscillations in the range of 100–200 kHz, with the characteristic time scale of the ion channel fly-by time, 100–200 kHz. The other test case demonstrates the effect of a finite temperature of injected neutral atoms, which has a substantial suppression effect on the oscillation amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Study of the Design and Characteristics of a Modified Pulsed Plasma Thruster with Graphite and Tungsten Trigger Electrodes

The paper presents experimental results for a modified pulsed plasma thruster (PPT) with solid propellant, using a coaxial anode–cathode design. Graphite from pencil leads served as propellant, and a tungsten trigger electrode was tested to reduce carbonization effects. Experiments were performed in a vacuum chamber at 0.001 Pa, employing diagnostics such as discharge current/voltage recording, power measurement, ballistic pendulum, time-of-flight (TOF) method, and a Faraday cup. Current and voltage waveforms matched an oscillatory RLC circuit with variable plasma channel resistance. Key discharge parameters were measured, including current pulse duration/amplitude and plasma channel formation/decay dynamics. Impulse bit values, obtained with a ballistic pendulum, reached up to 8.5 μN·s. Increasing trigger capacitor capacitance reduced thrust due to unstable “pre-plasma” formation and partial pre-discharge energy loss. Using TOF and Faraday cup diagnostics, plasma front velocity, ion current amplitude, current density, and ion concentration were determined. Tungsten electrodes produced lower charged particle concentrations than graphite but offered better adhesion resistance, minimal carbonization, and stable long-term performance. The findings support optimizing trigger electrode materials and PPT operating modes to extend lifetime and stabilize thrust output.

Faraday cup↗
DOE OSTI2024

Segmented wall-less hall thruster

Disclosed is a miniaturized plasma propulsion device with minimized surface area of the thruster walls exposed to the plasma and, as a result, reduced plasma-surface interactions including a set of segmented electrodes to facilitate the following improvements compared to relevant existing technologies: 1) control of the plasma flow including focusing of the plasma plume 2) increase of the thrust 2) reduction of inefficiencies associated with the electron cross field current, and 3) mitigation of low frequency oscillations. The electrodes affect all these actions when a DC or modulated voltage is applied to one or all of them with the same or different amplitudes, with the same or different frequencies or phases which are all optimized to realize the best performance through changes in the acceleration and/or ionization regions. In addition, the applied voltage to the main electrodes may also be modulated.

Simmonds, Jacob↗
DOE OSTI2023

Anomalous cross-field transport in a Hall thruster inferred from direct measurement of instability growth rates

The contribution of the electron drift instability to anomalous electron transport is experimentally assessed in a Hall effect discharge. The transport is represented by an anomalous collision frequency, which is related through quasilinear theory to the energy and growth rate of the instability. The wave energy is measured directly with ion saturation probes, while estimates of the growth rate are employed based on both linearized theory and direct measurement. The latter measurement is performed with a bispectral analysis method. The wave-driven collision frequency is compared to measurements of the actual collision frequency inferred from a method based on laser- induced fluorescence. It is found that estimates for transport using linearized theory for the growth differ by over an order of magnitude from the actual anomalous collision frequency in the plasma. The wave-driven anomalous collision frequency with measured growth, however, is shown to agree with the electron collision frequency in magnitude and capture aspects of the trends in spatial variation. This result demonstrates experimentally that wave-driven effects ultimately can explain the observed cross-field transport in these devices. As a result, the implications of this finding are discussed in the context of the key lengthscales that drive the transport as well as the implications identifying reduced fidelity models that could be used to predict anomalous collision frequency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Systems and methods for satellite movement

A satellite includes a plurality of thrusters disposed about the satellite, each of the plurality of thrusters having a minimum thruster firing time, and a control circuit connected to the plurality of thrusters. The control circuit is configured to identify violations of the minimum thruster firing time in a non-compliant thruster firing pattern selected to achieve a specified movement, generate a plurality of compliant thruster firing patterns by replacing each of the violations of the non-compliant thruster firing pattern by zero and a minimum time in different combinations, select a compliant thruster firing pattern from the plurality of compliant thruster firing patterns to produce a satellite movement that is within a predetermined range of the specified movement, and cause the plurality of thrusters to fire according to the compliant thruster firing pattern.

Cheng, Haifei↗
DOE OSTI2026

Theoretical assessment of the transition between electron emission mechanisms in a collisional crossed-field gap

Understanding the transition to the limiting current under vacuum and collisional conditions for thermal-field emission in a crossed-field gap is critical for characterizing the operating regime for numerous applications, such as high-power microwaves and Hall thrusters. By deriving the exact solution for the current density as a function of cathode temperature (and, by extension, initial velocity), cathode electric field, applied voltage, collision frequency, and magnetic field, we determine the transition from the canonical current source (thermal-field emission) to the vacuum and collisional limiting currents under different operating conditions. Here, we further show that the exact solution recovers the asymptotic solutions for the canonical and limiting currents for vacuum and collisional crossed-field gaps under appropriate limits. The exact solution demonstrates that the critical current decreases with both increasing magnetic field strength and increasing collisionality. Nexus phase space plots constructed by matching the individual asymptotic solutions provide the conditions for the transitions between mechanisms, elucidating when a single mechanistic equation would suffice and when theories combining multiple mechanisms are necessary.

Breen, Lorin I. [Purdue Univ., West Lafayette, IN ↗
DOE OSTI2022

Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling

In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E × B) flow relative to the unmagnetized ions can cause the so-called electron cyclotron drift instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in, for example, collisionless shock ramps in space, and in E × B discharge devices such as Hall thrusters. A prominent feature of ECDI is its capability to induce an electron flow parallel to the background >E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport is important due to its role in particle thermalization at space shocks, and in causing plasma flows towards the walls of E × B devices, leading to unfavorable erosion and performance degradation, etc. The development of ECDI and anomalous transport is often considered requiring a fully kinetic treatment. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a collisionless two-fluid framework without any adjustable collision parameter. By treating both electron and ion species on an equal footing, the free energy due to the inter-species velocity shear allows the growth of an anomalous electron flow parallel to the background E field. We will first present linear analyses of the instability in the two-fluid five- and ten-moment models, and compare them against the fully-kinetic theory. At low temperatures, the two-fluid models predict the fastest-growing mode in good agreement with the kinetic result. Also, by including more (> = 10) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on ion-to-electron mass ratio, plasma temperature, and background B field strength is also thoroughly explored. We then carry out direct numerical simulations of the cross-field setup using the five-moment model. The development of the instability, as well as the anomalous transport, is confirmed and in excellent agreement with theoretical predictions. The force balance properties are also studied using the five-moment simulation data. Here this work casts new insights into the nature of ECDI and the associated anomalous transport and demonstrates the potential of the two-fluid moment model in efficient modeling of E × B plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Alternative Power and Propulsion for Space Systems

What are the main sources of power and propulsion for space systems and what are their advantages, disadvantages, how do they work, and what are the politics behind them? The following topics will be analyzed solar panels, batteries, radioisotope thermoelectric generator (RTG), advanced Stirling radioscopic generator (ASRG), reactors, nuclear thermal rockets (NTP), thermionic experiment with conversion (TOPAZ), electrothermal plasma thruster (VASIMR), ion propulsion, and finally nuclear pulse propulsion. The production of Pu-238 is extremely controversial but a vital source of fuel for many of these space systems. Soral and Batteries are the most common and the least controversial forms of power used in space systems today. However, is ineffective for deep space exploration as they only function for up to four years without needing to be recharged. This is where nuclear systems come into play. Nuclear systems such as an RTG will “live” up to 88 years. For example, Voyager 1 is powered by an RTG It has been traveling since 1977 and has left out solar system in 2012 and has yet to run out of power. However, this type of technology (nuclear) has many hurtles it needs to jump through before advancements in can go much further.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Anomalous Retrograde Drifts in Obstructed Magnetron Microdischarges: a Consequence of a Field Reversal in the Anode Sheath?

The goal of this project was to carry out experiments at the Plasma Collaborative Research Facility (Princeton Plasma Physics Laboratory) aimed at confirming or refuting the hypothesis that an electric field reversal occurs in the near-anode region of a strongly obstructed planar magnetron microdischarges. Magnetron discharges are magnetized plasmas with E x B field topologies, and variants of such discharges are widely used in materials processing. The topologies are also similar to that of Hall thrusters used in space propulsion. Understanding the structure of these discharges will aid in developing reliable simulations that are intended to facilitate the further development of such technologies, as the field reversals, presumed to be due to strong axial-driven electron transport, provides further insight into transport properties. The possible existence of such a field reversal is hinted at because of the retrograde-propagating (i.e., negative E x B) spokes or gradient-driven drift instabilities seen in the high speed videos and segmented anode probing of the plasma. The presence of this field reversal should be reflected in the ion dynamics, particularly, in an ion drift towards the anode. The proposed experiments to be conducted at the PCRF are to characterize this anode-directed ion drift using state-of-the-art laser-induced fluorescence velocimetry. In particular, we seek to see if the near-anode ion velocity reverses direction when the external circuit is varied to alter the electron current demand imposed at the anode. Such an ion velocity reversal would confirm the mechanism, which is that the highly obstructed (small gapped) discharge produces strong axial gradients in plasma density, driving more electrons towards the anode then that demanded by the external circuit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Conference Support for 2024 MagNetUS Meeting

The 4th annual MagNetUS meeting was held from April 14-18, 2024 in Los Angeles, California at the UCLA Lake Arrowhead Conference Center. The goal of the meeting was to continue the successful template of a mix of invited talks, contributed talks, posters, and discussion sessions established by the first three MagNetUS meetings, in order to 1) communicate physics advances relevant for basic magnetized plasma research which can be performed at small and mid-size experimental facilities; 2) allow collaborator facilities to connect and interact with existing users and potential new users; and 3) grow the community of potential users of collaborator facilities, as well as scientists who focus on theory/modeling and spacecraft observations. The meeting was attended by 77 participants (71 in person) from 23 unique institutions and included 14 invited talks, 19 contributed talks, 22 posters, and six discussion sessions. Topics covered by the meeting included fusion research, space and astrophysical systems (including laboratory analogs), magnetic reconnection, plasma thrusters, turbulence, waves, dusty plasmas, laser plasmas, plasma vortices, plasma sources, and software ecosystems. Discussion sessions included discussions on science communication, with NSF and DOE program managers, on the MagNetUS application process for new users, on experimental facility capabilities, and as a tie-in to the Solar Wind Machine meeting.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Optimal Frequency-domain Analysis for Spacecraft Time Series: Introducing the Missing-data Multitaper Power Spectrum Estimator

While the Lomb–Scargle periodogram is foundational to astronomy, it has a significant shortcoming: the variance in the estimated power spectrum does not decrease as more data are acquired. Statisticians have a 60 yr history of developing variance-suppressing power spectrum estimators, but most are not used in astronomy because they are formulated for time series with uniform observing cadence and without seasonal or daily gaps. Here we demonstrate how to apply the missing-data multitaper power spectrum estimator to spacecraft data with uniform time intervals between observations but missing data during thruster fires or momentum dumps. The F-test for harmonic components may be applied to multitaper power spectrum estimates to identify statistically significant oscillations that would not rise above a white noise–based false alarm probability. Multitapering improves the dynamic range of the power spectrum estimate and suppresses spectral window artifacts. We show that the multitaper–F-test combination applied to Kepler observations of KIC 6102338 detects differential rotation without requiring iterative sinusoid fitting and subtraction. Significant signals reside at harmonics of both fundamental rotation frequencies and suggest an antisolar rotation profile. Next we use the missing-data multitaper power spectrum estimator to identify the oscillation modes responsible for the complex "scallop-shell" shape of the K2 light curve of EPIC 203354381. We argue that multitaper power spectrum estimators should be used for all time series with regular observing cadence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Vehicle for installing anchors in an underwater substrate

A method of installing one or more anchors in an underwater substrate in a body of water including installing an anchor into the underwater substrate by rotating an anchor installation vehicle about a central axis Y to drive the anchor coupled to the anchor installation vehicle into the underwater substrate. The anchor installation vehicle includes a vehicle frame having a top end and bottom end, a plurality of arms extending outward from the vehicle frame, one or more rotational thrusters disposed at distal ends of the respective arms, and an anchor system that holds the anchor extending from the bottom end of the vehicle frame with the anchor aligned with a central axis Y.

42 ENGINEERING↗