PPTX

Download Report

Transcript PPTX

Instability and Disappearance of Debye Sheaths Due to Secondary Electron Emission (SEE) 10/5/12 LTP Teleseminar by Mike Campanell, Grad Student, PPPL

Φ (x) (V)

wall plasma interior wall wall plasma interior

Sheath

x (m) Figure: Animation of the abrupt change of the electrostatic potential function during sheath instability.

No Sheath?

Figure: Potential profiles in a plasma slab bounded by floating walls with moderate emission (black plot) and very strong emission (

gray plot

Overview:

• General treatments of SEE effects on PSI [1,2] usually assume two things: (1) A familiar Debye sheath exists with surface potential below the plasma potential (2) The sheath is static in time • In this teleseminar, we will present simulations and basic theory showing that both assumptions can easily break down in practice.

Outline of Presentation:

• Quick background on SEE effects • Study of the inverse sheath effect (violation of assumption 1).

• Study of sheath instabilities (violation of assumption 2).

[1] G.D. Hobbs and J.A. Wesson, Plasma Phys. 9, 85 (1967). [2] L. A. Schwager, Phys. Fluids B 5, 631 (1993).

Background and Motivation:

• SEE is critical in many plasma physics applications. • The larger the SEE coefficient γ , the more plasma electrons must hit the wall to balance the ion flux: Γ e = Γ e,in – Γ e,out = Γ e,in (1 γ ) = Γ ion • Time-independent solutions of sheaths with fixed γ implicitly assume stability!

• For hot plasmas and/or highly emitting materials, γ may exceed unity [1]. In this case, a classical sheath, see image (a), cannot maintain zero current.

• Conventional SEE theories usually assume a “space charge limited” sheath arises . A dip forms in the potential profile Φ (x) near the wall to repel the “extra” secondaries back to the wall, see image (b).

• But this is not the only solution.

A Φ (x) in image (c) is possible. Ions are not drawn to the wall and more importantly, plasma electrons are unconfined, so the wall flux is thermal (extremely large)!

[1] C.A. Ordonez and R.E. Peterkin Jr., J. Appl. Physics. 79, 2270 (1996).

Fig: Qualitative differences between (a) the classical sheath, (b) SCL sheath and (c) the new inverse sheath.

Description of the Simulation Model:

• We use the EDIPIC code, written by D. Sydorenko [1], to simulate a planar xenon plasma bounded by walls made of boron-nitride.

• The plasma is heated by an ExB background field, as in a Hall thruster.

• When the ExB drift energy is low, classical sheaths form at the two walls.

• When the ExB drift energy is large, the electrons reaching the walls will have large energy parallel to the wall. So this could cause γ > 1.

We Simulate Two Cases:

Simulation A – Moderate drift energy case E z = 200V/cm, B x ,= 100G, Neutral Density n n = 10 12 cm -3 ,Plasma Density n p = 10 11 cm -3 Simulation B – High drift energy case E z = 250V/cm (larger V

Drift

), otherwise similar to Sim. A [1] D. Sydorenko, Ph.D. thesis, U. of Saskatchewan, 2006.

Results:

• The plasma in Sim. A has γ < 1 and features classical sheaths at the walls.

• The presheath/sheath structures in Sim. A form a potential well of large magnitude Φ = 21V to confine electrons and draw ions to the wall.

• In Sim. B, γ > 1. The sheaths vanish .

Φ(x) is small throughout the simulation domain. Just fluctuations appear.

• It is found that Φ(x) always decreases from the wall outward in Sim. B. The surface charge is positive, the net space charge near the interface is negative.

• Conclusion: Sim. B is an “inverse sheath”. Its amplitude is Φ -1 = -1V.

Fig: Comparison of Sim. A and Sim. B.

(b) Potential profiles Φ(x) .

(d) Charge density profiles in Sim. A.

(e) Charge density profiles in Sim. B

Further Analysis:

The reason for the disappearance of the Debye sheath is as follows:

• The wall flux in the simulations consists of the following components: Γ ion : ion flux Γ CE : collision-ejected electrons (CEE’s) Γ b : secondary “beam” electrons from the other wall • γ b γ b and γ CE denote the partial SEE yield of beam electrons and CEE’s. is determined by the ExB drift energy gained by secondaries crossing the plasma.

• It can be shown γ = γ CE /(1 + γ CE γ b ) • γ b = 0.94 in Sim. A. But

γ b

≈ 1.2 in Sim. B. If γ b > 1, classical sheaths cannot exist. • The “inverse sheath” forms to suppress some emission and maintain current balance.

• Γ ion is negligible because ions are not drawn to the wall.

• Note in the plot that γ net = 1 at all times even as the dominant electron flux component Γ b varies. Γ o , the cold secondaries pulled back to the wall by the inverse sheath, self-adjusts.

Time evolution of the flux components in Sim. B.

Conclusion:

• When the SEE yield of hot plasma electrons hitting the wall exceeds unity, a sheath is not needed to maintain current balance. • The “inverse sheath” (c) is fundamentally different from both classical sheaths (a) and SCL sheaths (b). The inverse sheath has a negative plasma potential . • Losses are far larger in the inverse sheath regime ~10-20 times that of Sim. A.

compared to classical and SCL regimes! It is found that the particle flux, energy flux and axial transport in Sim. B is Possible Applications: • The inverse sheath effect may be connected to the well known phenomenon of temperature saturation [1] in HT devices. Loss of sheaths should lead to enhanced wall losses and near-wall conductivity. • Inverse sheath is also possible for surfaces with thermionic or photoemission.

• In applications where ion sputtering is more damaging than electron heat flux, the inverse sheath regime would be beneficial.

[1] Y. Raitses et. al., IEEE Trans. on Plasma Sci. 39, 4, (2011).

Appendix: Why does the inverse sheath not appear in familiar SEE theories?

General SEE theories usually assume ions are drawn to the wall. Poisson’s eq. is then solved with charge densities in the sheath written in terms of a potential Φ(x) that is assumed below the plasma potential.

[1] G.D. Hobbs and J.A. Wesson, Plasma Phys. 9, 85 (1967).

Sheath Instability in Weakly Collisional Plasmas:

Sheath instabilities can (a) trigger abrupt changes in the plasma, (b) drive oscillations and (c) dramatically increase energy loss and cross-field transport. [1]

(V)

wall plasma interior Fig: Illustration of the sheath rapidly collapsing over 5ns.

Notice the abrupt drop in the sheath amplitude and the change in structure of the Φ(x) near-wall x (m) [1] M. Campanell, A. Khrabrov and I. Kaganovich, “General Cause of Sheath Instability Identified for Low Collisionality Plasmas in Devices with Secondary Electron Emission,” PRL 108, 235001 (2012).

Why Does Instability Occur?:

• Consider a planar plasma bounded by walls with SEE in a regime with classical sheaths.

• The arrangement of particles in energy space looks as in the phase plot to the right.

• Bulk plasma electrons (red) W x < e Φ are confined.

with Plot of particle energy components parallel to the walls (W y + W z ) and perpendicular to the walls (W x ). • The loss cone (W x > e Φ ) is strongly depleted in low collisionality. However, the secondary electrons form “beams” that pass between the walls. blue

Mechanism of Instability:

A decrease of the sheath potential Φ allows “weakly confined electrons” (WCE’s) to reach the walls. If the SEE coefficient of the WCE’s exceeds unity, the decrease of Φ causes reduction of the electron surface charge, which further lowers Φ, allowing more WCE’s to reach the wall etc… This leads to a runaway collapse of the sheath.

Verification of the Theory

• The electron flux to the wall consists of 3 components: (1) Collision ejected electrons (Γ CE ), scattered into the loss cone by collisions with neutrals.

(2) SEE “beam” electrons from other wall (Γ b ) (3) Weakly Confined Electrons nudged into the loss cone by fluctuations (Γ WC ) • By tracking each component and their SEE yield separately, it is seen that the system stability is dependent on the energies of the WCE’s.

• It was verified over many simulations that γ WC = 1 is the critical point of instability. Fig. Temporal evolution of Neutral Density n n Φ and partial SEE coefficients in a simulation with E z , = 140V/cm, B x ,= 100G, = 10 12 cm -3 ,Plasma Density n p = 10 11 cm -3 , turbulent heating rate ν turb = 5.6 MHz, H = 2.5cm

Consequence of Instability

• Instability triggers a rapid drop of Φ and a rapid increase of Γ WC (see Fig. below). • The instability quenches when the surface charge becomes positive and emitted secondaries start to get pulled back to the wall. • The system then may restore to a new quasisteady state, as is discussed in more detail in our paper. The new state generally has enhanced losses .

• Note that there are 3 different types of instability effects in the data plot below.

• All of these effects are caused by γ WC reaching or exceeding unity.

E z , = 200V/cm, B x ,= 100G, n n = 10 12 cm -3 , n p = 10 11 cm -3 ν turb = 4.2 MHz, H = 2.5cm

WCE Theory Reveals Cause of Past Observed Instabilities:

Beam Phase Instability 1 : • A one-time abrupt transition from one quasisteady state to another (point 2 in earlier Fig) • Φ drops drastically and the beam flight time switches to a new value.

Relaxation Sheath Oscillations 2 : • Quasiperiodic instabilities with steady behavior in between. (interval 3 of earlier fig) • During instability, there is a rapid decrease of Φ, an increase of Γ e,in transition to the SCL state.

and a [1] D. Sydorenko et al., Phys. Plasmas. 15, 053506, (2008). [2] D. Sydorenko et al., PRL 103, 14, (2009).

For More Information on Inverse Sheath:

If you would like a more complete discussion of the inverse sheath phenomenon, please see: M. D. Campanell, A. Khrabrov and I. D. Kaganovich, “Absence of Debye Sheaths Due To Secondary Electron Emission,” PRL 108, 2

001 (2012).

(there was a typo in the original email announcement of the talk )

For More Information on Instabilities:

For more detail on the instability phenomena, please see our other paper: M. D. Campanell, A. V. Khrabrov and I.D. Kaganovich, “General Cause of

Sheath Instability Identified for Low Collisionality Plasmas in Devices with

Secondary Electron Emission,” PRL 108, 2