Open access peer-reviewed chapter

Waves and Instabilities in E × B Dusty Plasma

By Sukhmander Singh

Submitted: April 24th 2019Reviewed: November 6th 2019Published: December 12th 2019

DOI: 10.5772/intechopen.90397

Downloaded: 61

Abstract

Hall thrusters are common examples of E × B configuration, where electron trajectory gets trapped along the external magnetic field lines. This significantly increases the residence time of electrons in the plasma discharge channel. Hall thrusters are potential candidates for spacecraft station keeping, rephrasing and orbit topping applications because of its high thrust resolutions and efficiency. The goal of this chapter is to explain the working principle of Hall thrusters and to characterize the resistive instability in hot dusty plasma. The studies of these instabilities are useful to design efficient Hall thrusters and to understand the solar dusty plasma. The large amplitude of these oscillations has an adverse effect on the power processing unit of the devices. This reduces the efficiency and specific impulse and shortens the operating life of the Hall thruster. The theory of linearization of fluid equation for small oscillation has been given. The chapter also discusses the origin of plasma oscillation in a plasma discharge mechanics.

Keywords

  • plasma oscillations
  • instabilities
  • Hall thrusters
  • resistive plasma
  • growth rate

1. Introduction to E × B devices

There are some devices where plasma is weakly magnetized because of the larger Larmor radius of the massive ions against the length of the discharge channel. Electrons get trapped along the magnetic field lines in the channel. Hall thrusters and magnetrons are the common examples of E × B devices. The E × B configuration is used to confine electrons, increasing the electron residence time and allowing ionization and plasma sustainment. Hall thrusters have emerged as an integral part of propulsion technology. Unlike chemicals and electric rockets, in a Hall thruster, the propulsion thrust is achieved by gas which has high atomic number and low ionization potential. For this Xenon is mostly used. In the discharge channel, Xenon is ionized and then accelerated by electrostatic forces. Hall thrusters are versatile electric propulsion devices, where thrust efficiencies can exceed 50% and specific impulses are typically between 20 min and 1 h. The specific impulse has the dimension of time and is a measure for the effective lifetime of the thruster, when lifting its own propellant from the earth’s surface. The specific impulse is defined as Isp=Tṁpg, where ṁpis the mass flow rate; T is the thrust, which is the total force undergone by the Hall Thruster in relation to the acceleration of the ions; and gis the acceleration due to gravity. If the propellant is fully ionized, the specific impulse is equal to the mean axial component of the final ion velocity divided by g. It is interesting that these thrusters adjust their thrust and impulse by varying the acceleration voltage and the flow rate of the propellant. The fine tuning of the thrust correction can be done for the compensation of atmospheric drag for low-flying satellites using its high thrust resolution. The necessary thrust for such applications ranges from micro-Newtons to some Newtons with electric input powers of some 10 to some 10,000 Watts [1, 2, 3, 4].

Figure 1 shows the internal components of a Hall thruster which is generally made of an axis-symmetric cylindrical discharge chamber. A cathode is fixed outside to produce electrons to neutralize the outer surface of the device to overcome the space-charging problems. A high atomic weight number and low ionization potential gases are preferred propellant (Xenon, Argon) for Hall thrusters to get more thrust. The propellant enters from the left side of the channel via anode and gets ionized through the hollow cathode of the device. The electric field of strength ∼1000 V/m gets generated inside the discharge channel along the axial direction of the device. By using magnets around the annular channel and along the thruster centreline, a radial magnetic field of moderate strength (∼150–200 G) is created, which is strong enough for the electrons to get magnetized, i.e. they are able to gyrate within the discharge channel, but the ions remain unaffected due to their Larmor radius much larger than the dimension of the thruster [3].

Figure 1.

Schematic diagram of a typical Hall plasma thruster [3].

We used a Cartesian coordinate system to understand the different forces on the particles inside the channel and let us suppose, the X-axis represents, the axis of the thruster. Generally, the applied electric and magnetic fields are in axial (along X axis) and radial (along Z axis) directions, respectively, of the device. Therefore because of the perpendicular electric and magnetic fields, the Lorentz forces act on the electrons along the Y axis (E×Bazimuthal direction). Since electrons have smaller Larmor radius than the length of the channel, therefore electrons rotate along the azimuthal direction under the influences of the Lorentz force [1, 2, 3, 4].

2. Review on plasma instabilities in Hall thrusters

It is well known that plasma pressure drives the instabilities in plasma. Therefore the confined plasma is prone to non-equilibrium thermodynamic state. Therefore it must be important to know the consequences of these instabilities. It has been established that the amplitude and frequency of the oscillations in the Hall thrusters depend on mass flow rate, discharge voltage, geometry, magnetic field profile and cathode operation mode. On the other hand, the plasma in a Hall thruster does not stay uniform, and an inhomogeneous plasma immersed in the external electric and magnetic fields is not in the thermodynamically equilibrium state; this deviation in general is a source of plasma instabilities.

In Hall thrusters, from low frequencies (few Hertz) to high frequencies (few GHz), oscillation spectra have been observed on theoretically as well as experimentally based studies. The oscillations in the range of 10–20 kHz are called as discharge oscillations, and oscillations in the range of 5–25 kHz are said to be ionization-driven oscillations. The drift instabilities and density gradient plasma are responsible to produce oscillations in the range of 20–60 kHz in a Hall thruster. The oscillations in the range of 70–500 kHz are also called transient time oscillations and are the order of ion residence time in the channel of the device. The oscillations associated with azimuthal waves are represented by high-frequency (0.5–5 MHz) oscillations [5]. Litvak and Fisch [6] have developed an analytical model for electrostatic and electromagnetic resistive instabilities in a Hall plasma for azimuthal disturbances. Singh and Malik investigated resistive instabilities for axial and azimuthal disturbances in a Hall thrusters [7, 8]. Fernandez et al. [9] did simulations for the growth of resistive instability. Litvak and Fisch [10] have analysed gradient-driven Rayleigh-type instabilities in a Hall thruster using two fluid hydrodynamic equations. Ducrocq et al. [11] have investigated high-frequency electron drift instability in the cross-field configuration of a Hall thruster. Barral and Ahedo [12] have developed a low-frequency model of breathing oscillations in Hall discharges, where they observed that unstable modes are strongly nonlinear and are characterized by frequencies obeying a scaling law different from that of linear modes. Chesta et al. [13] have developed a theoretical model to obtain the growth rate and frequencies of axial and azimuthally propagating plasma disturbances.

3. Studies of fine particles in plasma

The presence of heavy fine particles with a size of 1–50 microns and mass of orders 10−10 to 10−15 kg in a classical plasma acts as a external component in plasma. If the density of the dust particles is less than the plasma density, the system is called dusty plasma. These fine particles acquire some charges from the electrons to get charged. The magnitude of charge on dust grain is not constant. It depends on the type of dust grain, the surface properties of dust grain, the dust dynamics, the temperature, the density of plasma and the wave motion in the medium. The presence of fine particles in a plasma makes it more complex and these particles alter the dynamics of the plasma species which generate new propagating modes by exhibit their own dynamics. The dusty plasmas have an exciting property which has attracted researchers over the world in this area [14, 15, 16, 17, 18, 19, 20, 21, 22]. The presence of charged dust grains modifies the ion-acoustic waves, lower hybrid waves, ion-acoustic and introduces dust acoustic waves and dust ion acoustic waves [22]. Verma et al. have studied the electrostatic oscillation in the presence of grain charge perturbation in a dusty plasma [23]. They studied the property of electrostatic oscillation and instability phenomena taking into account the temporal evolution of the grain charge in an unmagnetized dusty plasma. Cui and Goree have studied the effect of fluctuations of the charge on a dust grain in plasma [24]. Sharma and Sugawa studied the effect of ion beam on dust charge fluctuations [25]. It is observed that growth rate of the instability increases with the relative density of negatively charged dust. If dust particle charge is eZd. The quasineutrality condition is given for dust particles by ne0=ni0+Zdnd0.

4. Electron plasma discharge oscillation

When electrons are displaced from the equilibrium position of the charged particles relative to the uniform background of the ions in plasma, an electric field is developed in such a direction that it tries to pull the electrons back to its equilibrium position to restore the neutrality. Because of the inertia effect, the electrons overshoot the equilibrium position, and now the electric field is developed in the opposite direction which again tries to pull back the electrons to their position of equilibrium. The massive ions are supposed to be fixed in the background and are not capable to respond the oscillating field generated by the oscillation of electrons. If n0is the number of electrons per unit volume in infinite sheet plasma of thickness x, let the electron be displaced (as shown in Figure 2) to the right from their equilibrium position which results to generate surface charge density σ=en0xon the left side of the sheet and equal and opposite surface charge density on the right side of the sheet. The generated electric field E=en0xε0tries to pull the electrons back, and thus oscillation takes place in plasma. From Newton’s second law, we write med2xdt2=eE=e2n0xε0. Thus the solution of the above second-order differential equation is given by xt=Acosωpet, where ωpe=e2n0meε0is the plasma frequency or character frequency at which imbalance charges oscillate. Using the values of various parameters, plasma frequency f(Hz) for electron is in the order 9n0(Hz), when the electron density is taken per m3. For example, for plasma having an electron density of 1018/m3, we have f9 GHz. Generally the plasma frequency lies in the microwave region. The above equation also shows that group velocity of these oscillations is zero, and hence there is no propagation of information. These waves are called stationary waves.

Figure 2.

Schematic of electron displacement in plasma sheet of thickness x.

5. Plasma oscillation when the motion of ions is also taken into account

Let x1and x2are the displacement from the equilibrium positions of electrons and ions, respectively, then the equation of motions is

med2x1dt2=e2n0x1x2ε0E1
mid2x2dt2=e2n0x1x2ε0E2

By combining the above equations, we obtain

d2x1x2dt2=e2n0x1x2ε01me+1miE3

or

d2Xdt2=e2n0Xmeε01+memiE4

where X=x1x2.

Therefore the frequency of oscillation is

ω2pie=e2n0meε01+memi=ω2pe1+memiE5

5.1 Frequency of oscillation for pair plasma

The pair plasma comprise of particles with opposite charge but equal mass, which gives plasma frequency ωie=2ωpe.

6. Concept of plasma resistivity

The equation of motion for electron in unmagnetized cold plasma can be given by the equation mυet+υeυe=eEmvυe. In a steady state, the collisional friction balances the electric acceleration, and the above equation results to E=mvυee. Since electrons move with respect to the ions, they carry the current density J=neυe. Substituting the current density into electric field equation yields E=ηJ. Here we have defined the plasma resistivity η=mvne2. Therefore plasma resistivity depends on the collision frequency between the neutral particles with plasma species.

7. Plasma model and basic equations

A Hall thruster with two-component plasma consisting of ions and electrons is considered in which only the electrons are magnetized and the ions are not. For the case of simplicity, the presence of dust particles has been ignored; otherwise the mathematical expression would become cumbersome. In order to realize the exact behaviour and the consequences of finite temperature on the thruster efficiency, it is of much importance to investigate the plasma disturbances in Hall thrusters by including the finite temperatures of the plasma species.

As discussed in Section 1, the electrons experience force along the azimuthal direction, and ions are accelerated along the exit side of the device to produce thrust by the external electric field. We use the common symbols to write the continuity and equation of motion for the ions and electrons under the thermal effects of ion and electron pressure gradient forces. The collision momentum transfer frequency (v) between the electrons and neutral atoms is also taken into account to see the resistive effects in the plasma:

nit+υini=0E6
υit+υiυi=eEMpiMniE7
net+υene=0E8
υet+υeυe=emE+υe×BvυepemneE9

8. Linearization of fluid equations

We consider the perturbed densities for ions and electrons as ni1and ne1velocities as υi1and υe1indicated by subscript 1 along with their unperturbed values as υ0and u0in the X and Y directions respectively. The unperturbed part represents the state of the plasma in the absence of oscillations and is indicated by a subscript 0. The unperturbed density is taken as n0, the electric field (magnetic field) as E0(B0) and the perturbed value of the electric field (magnetic field) as E1(B1). To linearize all the equations, let us write ni=n0+ni1, υi=υi1+υ0, B=B1+B0and E=E1+E0. In view of small variations of both the density and magnetic field along the channel, the plasma inhomogeneities are neglected. The perturbations of the ion and electron densities are taken small enough ni1ne1<<n0so that the collisional effect due to the velocity perturbations dominate over the one due to the density perturbation. Since υ0and u0are constant, the terms υ0n0, n0υ0and n1υ0are equal to zero. Further the terms υ1n1and n1υ1are neglected as they are quadratic in perturbation. The linearized form of the above equations thus reads

ni1t+υ0ni1x+n0υi1=0E10
υi1t+υ0υi1x=eE1MpiMn0E11
ne1t+u0ne1y+n0υe1=0E12
υe1t+u0υe1y=emE1+υe1×B0+u0×B1vυe1pemn0E13

The initial drifts υ0and u0of the ions and electrons in the channel are related to the electric and magnetic fields according to υ0υ0x=eE0Mand u0=E0B0obtained from the unperturbed part of Eqs. (7) and (9). The electron pressure in Eqs. (9) and (13) is given by pe=YeneTetogether with Teas the electron temperature, which we consider to be constant, and Yeas the ratio of specific heats.

Normal mode analysis: We seek the sinusoidal solution of the above equations; therefore the perturbed quantities are taken as f1expiωtikr. Then the time derivative /tcan be replaced by and the gradient by ik. Here f1ni1,ne1, υi1, υe1, E1and B1together with ωas the frequency of oscillations and kas the propagation vector.

9. Dispersion equation and growth rate of electrostatic oscillations

Since, we are only interested in electrostatic oscillations, and therefore in the meanwhile, the perturbed magnetic field can be ignored in Eq. (13). By using Fourier analysis in Eq. (10) and Eq. (13), the perturbed ion and electron densities are given as follows:

ni1=n0ωkxυ0kxυi1x+kyυi1yE14
ne1=n0ωkyu0kxυe1x+kyυe1yE15

Using Eq. (11) into Eq. (14) gives

ni1=ek2n0φMωkxυ021k2VthI2ωkxυ021E16

where we used k2=kx2+ky2.

The expression for the electron density ne1contains the velocity components υe1xand υe1y, which are derived in terms of the potential φunder the assumption Ω>>ω, kyu0and vin view of the oscillations observed in Hall thrusters [13, 22, 23]. Letting ωkyu0ivω̂, eB0mΩ, YeTemVth, n0e2mε0ωeand n0e2Mε0ωi.

we use Eq. (13) to write the velocity components

iωkyu0ivυe1x=emφxΩυe1y+ikxne1VthE2n0E17
iωkyu0ivυe1y=emφy+Ωυe1x+ikyne1VthE2n0E18

Further simplification gives

υe1x=1mΩ2iekyΩφimω̂kxVthE2ne1n0+kxω̂imΩkyVthE2ne1n0+iω̂2mΩ3ekyφmkyVthE2ne1n0E19
υe1y=1mΩ2imkxVthE2Ωne1n0+kyω̂ikxΩmω̂kyVthE2ne1n0+iω̂2mΩ3mkxVthE2ne1n0kxE20

With the above velocity components, the perturbed electron density ne1can be expressed in terms of perturbed potential as follows:

ne1=en0ω̂k2φmΩ2ωkyu0+mω̂k2VthE2E21

Finally, we use the expressions for the perturbed ion density ni1and electron density ne1in Poisson’s equation ε02φ=ene1ni1in order to obtain

k2φ=ωe2ω̂k2φΩ2ωkyu0+ω̂k2VthE2ωi2k2φωkxυ02k2VthI2E22

Since the perturbed potential is φ0, we have from Eq. (21)

ωe2ω̂Ω2ωkyu0+ω̂k2VthE2+ωkxυ02k2VthI2ωi2ωkxυ02k2VthI2=0E23

This is the dispersion relation that governs the electrostatic waves in the Hall thruster’s channel.

9.1 The limiting case

For smaller oscillations, that is, ω<<kyu0, the above relation yields

ωkxυ02=kyu0+ivk2VthI2ωe2+k2VthE2ωi2+k2VthI2+Ω2kyu0ωi2+k2VthI2kyu0+ivωe2+k2VthE2+Ω2kyu0E24

Now, using the conditions Ω<ωe>ωiand VthI<<VthEin the above equation and letting ω1=Ω2ωi2+ωe2k2VthI2Ω2+ωe2,

ωkxυ02ω121+ivωi2k2VthE2kyu0Ω2ωi2+ωe2k2VthI21+ivωe2kyu0Ω2+ωe2E25

Since the last terms in the second brackets of the numerator and denominator in the right-hand side of Eq. (25) are small, we obtain the following

ωkxυ0±ω11+ivk2VthE2ωi22kyu0Ω2ωi2+ωe2k2VthI21ivωe22kyu0Ω2+ωe2E26

10. Instability analysis

The roots ωjof the above polynomial may be real and/or complex number. For a real root, sinusoidal behaviour gives expiωjtikrshowing an oscillatory behaviour. For complex roots which always occur in complex conjugate, we write ωj=ωrjiγj, where ωrjand γjare real numbers and the sinusoidal behaviour becomes etγjexpiωrjtikrwhich is showing an exponentially growing wave for a positive value of γjor an exponentially damped wave for a negative value of γj. Since the complex roots always appear in a conjugate pair, i.e. we have both negative and positive values of γjsimultaneously. Hence, for the complex roots, one of the waves is always unstable.

Finally, the growth rate γof the resistive instability is calculated from Eq. (26) as below:

γvω12kyu0ωe2Ω2+ωe2ωi2k2VthE2Ω2ωi2+ωe2k2VthI2E27

The corresponding real frequency is obtained as

ωrkxυ0±ω11+v2k2VthE2ωi2ωe24ky2u02Ω2+ωe2Ω2ωi2+ωe2k2VthI2E28

From Eq. (27), it is obvious that growth rate is directly proportional to the collisional (dissipative effects) frequency of the electrons which depends on various plasma parameters.

The results obtained in Eq. (27) matches with Litvak and Fisch [6] when the thermal effects become ignorable (i.e. Ti=Te=0). Under this situation, the growth rate takes the form

γv2kyu0Ω2ωi2Ω2+ωe2ωe2Ω2+ωe2.E29

In terms of lower hybrid frequency ωLH=Ω2ωi2Ω2+ωe2and in the limit Ω<<ωe, the growth rate can be written as γ±ωLHv2kyu0.

The above relation matches with Eq. (21) of [6].

Since it is not possible to find an analytical solution of the above equation, we look for the numerical solution along with typical values of B0, n0, Te, u0,ky, vand υ0. In Hall plasma thrusters, these parameters can have the values as thruster channel diameter = 4–10 cm, B0=100200 G, n0=5×10171018/m3, Te=1015 eV, u0106 m/s, v106/s and υ0=2×1045×104 m/s [6, 7, 8, 9, 10]. With regard to the value of ky, we constraint ky=m/r(where ris the radius of thruster channel) together with m=1for the azimuthal mode propagation [6, 7]. Accordingly we set ky=20/m; however, for higher mode (m > 1) or larger value of ky, the wavelength would be much smaller than the azimuthal dimension of the channel. The variations of the resistive growth rate are seen in the figures by solving Eq. (23) numerically. Some propagating modes and instabilities with smaller growth are also observed during the analysis of Eq. (23).

Figure 3 confirms that the growth rate of the instability gets enhanced with the increase of collision frequency of the electrons due to the resistive coupling of the oscillations to the electron azimuthal drift. The growth rate also increases with the increase of the electron temperature (Figure 4), and it also increases with higher electron density of plasma. Therefore it can be concluded that the collisional effect is responsible to unstable the plasma system. The numerical value of the growth under the collision frequency is observed in the order of ∼105/s.

Figure 3.

Variation of growth rate γ with collision frequency for different values of azimuthal wave number in a plasma having Xe ions (M = 131 amu), when Te=10 eV, n0=1018/m3, u0=106 m/s and B=0.02 T.

Figure 4.

Variation of growth rate γ with electron temperature for different values of plasma density in a plasma having Xe ions (M = 131 amu), when v=106/s, ky=20/m, u0=106 m/s and B=0.02 T.

11. Conclusions

In conclusions, we can say that the waves propagating in azimuthal and axial direction in a Hall thruster channel become unstable due to the resistive coupling to the electrons’ E×Bflow in the presence of their collisions. By controlling the various parameters, the growth rate and the propagating frequency of the oscillation may be controlled to optimize the performance and lifetime of the device.

Acknowledgments

The University Grants Commission (UGC), New Delhi, India, is thankfully acknowledged for providing the startup grant (No. F. 30-356/2017/BSR).

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Sukhmander Singh (December 12th 2019). Waves and Instabilities in E × B Dusty Plasma, Thermophysical Properties of Complex Materials, Aamir Shahzad, IntechOpen, DOI: 10.5772/intechopen.90397. Available from:

chapter statistics

61total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Structural Disorder as Control of Transport Properties in Metallic Alloys

By Eric R. Kaiser and Yong W. Kim

Related Book

First chapter

Introductory Chapter: A Novel Approach to Compute Thermal Conductivity of Complex System

By Aamir Shahzad, Syed Irfan Haider, Maogang He and Yan Feng

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us