Cylindrical Hall Thruster (cylHallThruster.sdf)

Keywords:

electric propulsion, Hall thruster channel, erosion models.

Problem description

Electric Hall thrusters are used for in-space propulsion and satellite station-keeping needs. The discharge plasma inside the Hall thruster channel is produced by the ionization of electrons with a neutral propellant gas such as xenon. The electrons are emitted from the neutralizer cathode placed at the exit of the Hall thruster (cathode end). The neutral gas is fed into the channel from the anode end of the Hall thruster channel. The electrons are confined inside the Hall thruster channel by the radial magnetic field applied through the solenoidal magnetic fields. Plasma xenon ions are accelerated out of the channel at high velocity, which produces the thrust necessary for space propulsion. Recently these thrusters are being designed to support long life time, high-power and high-thrust operations. The channel wall erosion occurring inside of the Hall thruster is one of the main limitations to these design needs. It becomes important to understand the plasma discharge processes occurring inside the Hall thruster channel and predict the lifetime of the Hall thruster based on the calculations of sputtered material from the Hall thruster channel.

This example demonstrates the xenon discharge plasma processes of a Stationary Plasma Thruster (SPT-100) channel. The outer cylinder has a radius of 5 cm and the inner cylinder has a radius of 3.5 cm. The radial channel gap is 1.5 cm. The channel length is 2.5 cm and our simulation domain is extended up to 1 cm outside the channel region (to simulate both channel and plume plasma). The left wall is biased at 300 V anode potential. Both inner and outer cylinders are maintained as dielectric cylinders (with hexagonal boron-nitride dielectric material coating) with a dielectric permittivity ratio of 4.6. The exit boundary is taken at 0 V. An electron source is placed at the channel exit for considering the electron emission from the neutralizer cathode. The cathode emission current is taken as 4.5 A. The xenon neutral gas is considered as a static fluid background. The neutral gas density is taken to be linearly decreasing trend with a maximum gas density is taken at the anode end of the channel. The simulation is initiated from a uniform plasma with both electrons and xenon ions. We enable a self-similar scaling system laws for the simulation of Hall thruster channel described by figure 1 in [Mahalingham2011]. This is based on earlier work by Taccogna [Taccogna2005] [Taccogna2004]. We have taken a scale factor of 1/50, i.e., the thruster dimensions are scaled by 1/50. This scaling is followed so that the kinetic Hall thruster channel plasma simulations can be performed in a reasonable run time. The scaling affects the physical dimensions and external potentials.

This example does demonstrate ionization of neutral fluids in cylindrical geometry.

This simulation can be performed with a VSimPD license.

Opening the Simulation

The Cylindrical Hall Thruster example is accessed from within VSimComposer by the following actions:

• Select the New → From Example… menu item in the File menu.

• In the resulting Examples window expand the VSim for Plasma Discharges option.

• Expand the Spacecraft option.

• Select Cylindrical Hall Thruster and press the Choose button.

• In the resulting dialog, create a New Folder if desired, and press the Save button to create a copy of this example.

All of the properties and values that create the simulation are now available in the Setup Window as shown in Fig. 558. You can expand the tree elements and navigate through the various properties, making any changes you desire. The right pane shows a 3D view of the geometry, if any, as well as the grid, if actively shown. To show or hide the grid, expand the Grid element and select or deselect the box next to Grid.

Fig. 558 Setup Window for the Cylindrical Hall Thruster Channel example.

Simulation Properties

This example contains many user defined Constants which help simplify the setup and make it easy to modify. These include constants such as:

• B0: The amplitude of the background magnetic field

• anodeV: the anode voltage

• innerRad and outerRad: inner and outer cylinder radius

• xeMaxDensity: the maximum density of the background Xe fluid

There are also several SpaceTimeFunctions that are used to define spatially and/or temporally varying inputs to other properties. These include:

• By: the magnetic field profile

• initialGasDensity: the profile for the background gas density

The self-consistent electric field is solved from Poisson’s equation by the electrostatic solver in a cylindrical coordinate system. The simulation is performed in axisymmetric 2-D fashion. The plasma is represented by macro-particles which are moved using the Boris pusher in cylindrical coordinate system. Various types of elastic and inelastic collisions of the particles are also taken into account.

Running the simulation

After performing the above actions, continue as follows:

• Proceed to the Run Window by pressing the Run button in the left column of buttons.

• To run the file, click on the Run button in the upper left corner of the Logs and Output Files pane. You will see the output of the run in this pane. The run has completed when you see the output, “Engine completed successfully.” A snapshot of the simulation run completion is shown in Fig. 559.

Fig. 559 The Run Window at the end of execution.

On an 8 core I7-6000 CPU, it takes about 40 minutes to run for 6,000 time steps. To reach steady state, about 100,000 time steps are required.

Analyzing the Results

If it is desired to calculate the density of the electrons or ions the analysis script computePtclNumDensity.py must be used.

• First click on the Analyze Tab.

• From the Available Anaylzers list, choose computePtclNumDensity.py. Then click Open.

• This script accepts the simulationName (Name of the input file) and speciesName to be calculated (species of particles).

• To calculate the density of the electrons, set the simulationName to “cylHallThruster” and the speciesName to “electrons”.

• Click on the Analyze button at the top right of the Analysis Results pane.

• A snapshot of the simulation run completion is shown in Fig. 560.

Fig. 560 The Analyze Window at the end of execution.

The resulting data will be visualizable as electronsDensity under the Scalar Data menu in the Visualize Tab. The density of other particles such as heavyIons, or Xeplus can also be calculated if those species names are used in place of electrons.

Visualizing the Results

To visualize the results, continue as follows:

• Proceed to the Visualize Window by pressing the Visualize button in the left column of buttons.

There are many different fields, particles, and histories that can be visualized in this example. The horizontal axis represents Z direction and the vertical axis represents R direction.

To view the electric potential, switch to the Data View to Data Overview in the Controls pane. Expand Scalar Data and choose Phi. Move the Dump slider to the right most position. Fig. 561 shows the visualization seen for the electric potential of the cylindrical Hall thruster channel and in the exit region.

Fig. 561 Visualization of Cylindrical Hall thruster channel electric potential results.

In the Hall thruster channel plasma, the electrons injected from the right end (i.e., exit of the channel) are accelerated towards the anode biased wall at the left end. To plot the particles, expand Particle Data and select “electrons” and “XePlus” check boxes The figure below, Fig. 562, shows the distribution of positively charged xenon ions (red dots) and electrons (blue dots).

Fig. 562 Visualization of Cylindrical Hall thruster channel particle phase-space distribution results.

The static radial magnetic field distribution considered for the SPT-100 Hall thruster channel set up is shown in Fig. 563. To reproduce this plot, expand Scalar Data then nodalB and select nodalB_r. The magnetic field is strong near the inner cylinder and has a Gaussian bell-shaped field distribution both inside and at the exit of the channel.

Fig. 563 Visualization of radial magnetic field in the Cylindrical Hall thruster channel.

The background xenon neutral fluid density distribution (plottable as XeNeutralFluid under Scalar Data) used in the simulation set up is shown in Fig. 564. The maximum neutral fluid density is taken at the left end of the channel near the anode wall. A linearly varying neutral fluid density is assumed.

Fig. 564 Visualization of xenon neutral fluid density in the Cylindrical Hall thruster channel.

Further Experiments

This example can be modified to test different design parameters such as varying anode voltages, varying background neutral gas densities and varying electron emission currents. This will allow users to study high-to-low power and high-to-low throttle levels.

Also the background fluid type can be changed to investigate other fluid species’ in this simulation set up.

Taccogna2005

Taccogna, F., Longo, S., Capitelli, M., and Schneider, R., “Plasma Flow in a Hall Thruster,” Physics of Plasmas, 112, 2005.

Taccogna2004

Taccogna, F., Longo, S., Capitelli, M., and Schneider, R., “Stationary plasma thruster simulation,” Comp. Phys. Comm., 165, pp. 160-170, 2004.