1 / 55

Complex phenomena in magnetized plasmas with an electron emission

Complex phenomena in magnetized plasmas with an electron emission. Yevgeny Raitses Princeton Plasma Physics Laboratory. Michigan Institute for Plasma Science and Engineering. Ann Arbor, December 5, 2012. Plasma Science & Technology Research

bella
Download Presentation

Complex phenomena in magnetized plasmas with an electron emission

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Complex phenomena in magnetized plasmas with an electron emission Yevgeny Raitses Princeton Plasma Physics Laboratory Michigan Institute for Plasma Science and Engineering Ann Arbor, December 5, 2012

  2. Plasma Science & Technology Research at Princeton Plasma Physics Laboratory (PPPL) Heavy ion beam MRI

  3. Outline • EB plasma devices: • Configurations • Electron rotating effects • Maximizing electric field applied in plasma • Anomalous electron cross-field transport: • Secondary electron emission effects • Turbulent fluctuations and coherent structures • Suppression of anomalous electron transport • Summary and concluding remarks

  4. PPPL DC-RF EB discharge of Penning-type DC E×B fields applied in a 20 cm × 50 cm st. steel chamber with ceramic side walls Plasma cathode: 2 MHz, 50-200 W Ferromagnetic ICP Operating parameters: Bkg. pressure: 0.1-1 mtorr RF-power: 50-60 W DC voltage/current: 0-100 V/0-3 A Magnetic field: up to 500 Gauss Anode Coils Coils Magnetically shielded RF-plasma cathode B Insulator E Axis

  5. Plasma in E ×B region: weakly collisional, non-equilibrium, with magnetized electrons and non-magnetized ions Neutral density ~ 1013 cm3 Plasma density ~ 0.5-3  1011 cm-3 Electron temperature ~ 3-5 eV Magnetic field: 5-500 Gauss • ea/L ~ 1-2 • ei/L ~ 10 • ee/L ~ 20-50 • ia/L~ 0.5-3 • Energy relaxation length in inelastic range  > * • */L ~ 2 Electron cross-field displacement during time loss (inelastic or wall collisions) X ~ 2RLe (scat /2loss )0.5 For B = 35 Gauss ce/coll ~ 150-200 4/14

  6. Examples of E ×B devices Large Plasma Device (LaPD) at UCLA 20-meter long, 1 meter diameter Penning Gauge Sputtering magnetron discharge

  7. Hall Thruster (HT) – fuel effective plasma propulsion device for space applications e e e<< L << i E =-ve B Diam ~ 1 -100 cm B ~ 100 Gauss Working gases: Xe, Kr Pressure ~ 10-1 mtorr Vd ~ 0.2 – 1 kV Power ~ 0.1- 50 kW Thrust ~ 10-3 - 1N Isp ~ 1000-3000 sec Efficiency ~ 6-70% • Unlike ion thruster, HT is not space-charge limited • Thrust density is limited by B2/2

  8. Parameters of HT plasma ea/h ~ 20 – 200 ei/h ~ 4103 ia/h ~ 10-100 Energy relaxation length in the inelastic range */h ~ 30 - 300 Neutral density ~ 1012-1013 cm3 Plasma density ~ 1011-1012 cm-3 Highly ionized flow: ion/n~ 80% Electron temperature ~ 20-60 eV Ion temperature ~ 1 eV Ion kinetic energy ~ 102-103 eV Collisionless, non-equilibrium plasma with magnetized electrons and non-magnetized ions

  9. Comparison of different E B plasma devices

  10. Cylindrical Hall thruster (CHT) – EB plasma in diverging magnetic field • Similar to conventional HTs, the CHT operation is based on closed electron EB drift. • Fundamentally differences from conventional HTs: Electrons are confined in the magneto-electrostatic trap. • Ions are accelerated in a large volume-to-surface channel Related concepts DCF by MIT and HEMP by Thales, CHT by Osaka, etc. Raitses and Fisch, Phys. Plasmas 8, 2579 (2001)

  11. Unusual focusing of the plasma flow in diverging magnetic field of CHT LIF measurements of ion velocity Ion current in plume Spektor et al., Phys. Plasmas 17, 093502 (2010) Raitses at al., Appl. Phys. Lett. 90, 221502 (2007)

  12. Plasma with azimuthal symmetric magnetic field and E×B rotating electrons is common in industrial and laboratory plasmas: non-neutral plasmas, solar physics, magnetic mirrors, magnetic fusion devices, plasma centrifuges and, most recently, plasma thrusters

  13. Rotating electron effects Isorotation For magnetized electrons and non-magnetized ions, common assumption is that magnetic surfaces are equipotential surfaces leads to a force field that is perpendicular to the magnetic surfaces, a good assumption for non-rotating cold magnetized plasma

  14. Ion focusing due to rotating electron effects Pressure gradient Centrifugal force effect on electrons Non-magnetized ions are not affected by the magnetic field, but the addition of the field Es results in focusing deflection of the original electric field En Ion focusing should benefit from supersonic electrons Fisch et al., PPCF, 53, 124038 (2011)

  15. Challenging requirements for the generation of supersonically rotating electrons in a steady state - Strong electric field and low magnetic field to get high EB speed - Colder plasma Common approach: Control of E-field with biased electrodes

  16. HT is capable to generate supersonically rotating electrons

  17. Electric field and thruster performance are affected by anomalous electron cross-field transport • Thruster efficiency With all other parameters held constant, HTs efficiency reduces with increasing electron current across the magnetic field • Classical collisional mechanism can not explain the discharge current measured for Hall thrusters: e-a~ 106 s-1 < eff~ 107 s-1 • Enhanced cross-field conductivity in HTs usually attributed to • SEE induced near-wall conductivity • Anomalous (Bohm-type) diffusion induced by high frequency azimuthal plasma oscillations • A new route for electron transport across magnetic field - low frequency rotating spoke oscillations

  18. Effect of the channel wall material on the discharge characteristic Carbon segments drastically change V-I characteristics • Boron nitride - high SEE • - Carbon velvet - zero SEE Raitses, Staack, Dunaevsky, Fisch, Phys. Plasmas (2006)

  19. Wall material affects the maximum electron temperature in the thruster PPPL Hall thruster setup Electron temperature from emissive probe measurements Raitses , Staack, Smirnov, Fisch Phys. Plasmas ,2005

  20. Pz26 + Pz26 -  Boron Nitride Teflon SEE from dielectrics reaches 1 at lower energies (< 50 eV) of primary electrons than for metals PPPL SEE setup Note: for boron nitride, if primary electrons are Maxwellian (Te) 1 at Te = 18.3 eV Dunaevsky, Raitses, Fisch, Phys. Plasmas (2003)

  21. SEE can significantly enhance electron flux from plasma to the wall When When SEE turns sheath to space-charge limited regime [Hobbs and Wesson, 1967] Fluid Approach w(x) scs Te i e see

  22. SEE effect on plasma electrons: comparing experiment with predictions Fluid theory Temax 18.3 eV According to fluid theories, the maximum electron temperature should not be above 18.3 eV (for BN and Xenon) Large quantitative disagreement with fluid theory!

  23. EVDF in HT is strongly anisotropic with beams of SEE electrons Hall thruster plasma, 2D-EVDF Isotropic Maxwellian plasma, 2D-EVDF Loss cones and beams Sydorenko et al, Kaganovich et al., Phys. Plasmas (2005, 2006, 2007 2009), Ahedo, Phys. Plasmas (2005)

  24. (x) i i 1p 1p 2 1b 1b 2 Total emission coefficient: Electron fluxes have several components, including counter-streaming SEE beams from opposite walls 1- primary 2- secondary SEE coefficients: p 2p / 1p - SEE due to plasma electrons b 2b / 1b - SEE due to beam electrons  1b / 2 - Penetration of the SEE beams Note, p can be > cr if eff < cr

  25. Conditions for the existence of self-sustained counter-streaming SEE electron beams 1) Weak two-stream and plasma beam instabilities PIC simulations predict: • EVDF is decreasing f (vx) • Beam penetration is high, 0.9 unstable f(vx) vx stable f(vx) vx Sydorenko et al., Phys. Plasmas 2007

  26. Conditions for self-sustained counter-streaming SEE electron beams (Cont’d) 2) Sufficiently strong electric field • SEE electrons gain additional energy during the flight between the channel walls due to EB motion • This energy must be high enough to induce strong SEE on opposite wall • The maximum additional energy is scaled as • For typical HT conditions: E = 100-200 V/cm, B ~ 100 Gauss • Bmax~ 30-60 eV enough for strong SEE from any ceramic material

  27. Near-wall conductivity SEE-induced cross-field current Wall collisionality - exchange of primary magnetized electrons by non-magnetized SEE electrons during the flight timeH/ubx The displacement , , gives average velocity B and current E Kaganovich, Raitses, Sydorenko. Smolyakov, Phys. Plasmas (2007)

  28. Two profiles for two regimes of SEE-induced electron cross-field current Predicted profiles of the cross-field current density: Classical sheath with SEE E = 200 V/cm Inverse sheath at a very strong SEE  > 1, E = 250 V/cm

  29. Disappearance of near-wall sheath at a very strong SEE  > 1 Qualitative differences between the potential profile, relative to the wall, of a classical sheath (a), SCL sheath (b) and the new inverse sheath (c). Note that plasma electrons are still confined by the SCL sheath, but not confined by the inverse sheath. Results of particle-in-cell simulations of Hall thruster discharge: a comparison of results with classical (Sim. A), E = 200 V/cm, and inverse sheath (Sim. B) E = 250 V/cm. M. Campanell et al., Phys. Rev. Lett. 108, 255001 (2012)

  30. When plasma is bounded with non-emitting and zero-recycling (100% absorbing) walls Engineered materials to mitigate plasma-surface interaction effects, e.g. carbon velvet material Carbon fibers bonded to carbon substrate • Low SEE because: • Carbon has low SEE • SEE electrons are trapped • in inter-fiber micro cavities • Low back flux of contamination: • Ion grazing incidence • Redep. is trapped in velvet texture Raitses, Staack, Dunaevsky, Fisch, Phys. Plasmas (2006)

  31. Without SEE, the magnetized plasma can withstand much stronger electric field No-SEE High-SEE 2.5 cm Probe path -4.6 cm 0 With No-SEE walls, the electric field at high voltages,  1 kV/cm, approaches a fundamental limit for a quasineutral plasma: E ~ Te/D (Te ~ 100 eV, ne ~ 1011 cm-3)

  32. Without SEE, the cross-field mobility reduces to almost classical collisional level Experimental cross-field mobility estimated using measured data and 1-D Ohm’s law at the placement of Emax Possibly EB shear effect?* For No-SEE, the shearing frequency, d(Ez/Br)/dz, reaches 5-8 nsec-1 at 600 V Such a large shear may affect the dynamics of all instabilities, which were previously predicted for Hall thrusters at moderate voltages *Fernandez, Cappellli, et al., Phys. Plasmas 15, 2008

  33. HT is capable to generate supersonically rotating electrons

  34. How azimuthal oscillations can cause cross-field transport? • In principle, HT discharge is azimuthallyy symmetric • If there are azimuthal oscillations of ne and  and they are correlated so that their time average over one period is nonzero, a wave-based azimuthal force appears: • For • Therefore, the F×B drift of that wave-based force could be responsible for collisionless cross-field transport

  35. Hall Thruster Oscillations Oscillations in Hall thruster plasma

  36. Imaging of HT operation PPPL Hall Thruster Experiment (HTX) Xenon operation of 12 cm diameter 2 kW PPPL Hall thruster Phantom camera V7.3 • Records 400,000 fps • Unfiltered emission • ~ 7.5 m away

  37. High speed imaging of HT operation 12 cm diameter PPPL HT 300 V, 20 sccm Xenon 100 Gauss 700 W Steady-state operation

  38. Rotating spoke Azimuthal non-uniformity of visible light emission and plasma density rotating in EB direction (~ 10 kHz) observed using fast cameras and electrostatic probes for different types of HTs Low voltage operation (< 200 V), probes Janes and Lowder, Phys. of Fluids 9 (1966) Morozov, et al, Sov. Phys. Tech. Phys. 5 (1973) Meezan, Hargus, Cappelli, Phys. Rev E 63 (2001) Modern HTs, > 200 V, fast imaging and probes Parker, Raitses, Fisch, Appl. Phys. Lett., 97 (2010) McDonald and Gallimore, IEEE TPS, 11 (2011) Ellison, Raitses, Fisch, Phys. Plasmas 19 (2012) Griswold et al Phys. Plasmas 19, (2012) Theory and simulations of low frequency azimuthal oscillations Escobar and Ahedo, IEPC 2011 Matyash , Schneider et al., IEPC 2011 Vesselovszorov, IEPC 2011 Spoke is always ~ 10 times slower than local EB speed !

  39. A possible mechanism of cross-field transport through the spoke E0z×B + + - - + - + - + - + - Br Eθ E0z • Possible transport mechanism through the spoke: • Initial density perturbation • Only electrons undergo azimuthal drift motion • Eθ generated across the perturbation • Eθ×B drift across the magnetic field, towards the anode Eθ×B • Correlated density and azimuthal electric field fluctuations would explain enhanced electron transport

  40. Cross-field transport through coherent plasma structures in magnetically controlled plasmas Non-diffusive transport - particles are not moving by a random walk (drift wave fluctuations), but rather form coherent structures (or blobs) that convect towards the walls Evolution of turbulent structures at the edge of the NSTX tokamak Serfanni et al, PPCF 49, 2007, Photo: Courtesy of S. Zweben UCLA LAPD MISTRAL, Aix-Marseille Univ. (EB linear device) HIPIMS RU, Bohum Jaeger, Pierre, Rebont, Phys. Plasmas 16 (2000) Carter, Phys. Plasmas 13 (2006)

  41. Cylindrical Hall thruster (CHT) • Mirror-cusp magnetic field topology • Similar to conventional HTs, the operation involves closed EB electron drift • Electrons are confined in the hybrid magneto-electrostatic trap • Ions are accelerated in a large • volume-to-surface area channel • (potentially lower erosion) 100 W 2.6 cm CHT Raitses and Fisch, Phys. Plasmas 8, (2001) Cathode

  42. Rotating spoke in CHT • Direction: ExB • Frequency: 15-35 kHz • Velocity: 1.2-2.8 km/s • E/B: 10-30 km/sec • E/B frequency 100-500 kHz • Size: 1.0-1.6 cm Cusp: Enhanced Radial Field

  43. Does spoke conduct current ? • Rotating spoke can not be observed in the discharge current traces • Segmented anode (4 isolated segments) allows to see the rotating spoke • Synchronized measurements with the fast camera reveal spoke-induced current More than 50% of the discharge current is conducted via the spoke Similar results were obtained for cylindrical and annular Hall thrusters Ellison, Raitses, Fisch, Phys. Plasmas 19 (2012)

  44. Slow movable probe Stationary probe arrays Insight of spoke with probes Plasma density oscillations by planar tungsten probes Plasma potential oscillations by floating emissive probe Inside the channel: probe tips are flush with channel wall Outside the channel: probe tips are at radial position of channel wall 3 azimuthal probes, 90 degrees apart, per axial location 2 azimuthal probes, 30 degrees apart, on a movable positioner outside the channel back middle front 13 mm 23 mm

  45. Spoke is everywhere along the channel, but the coherent rotation is only near the anode Density fluctuations: S(kθ,ω) kθ>0 corresponds to E×B direction Anode region Channel middle Cathode region • An azimuthal mode does exist in all three regions. • Mode is strongest in the back, although also “noisy” and extends over a large frequency range Parker, Raitses, Fisch, Appl. Phys. Lett., 97 (2010)

  46. - Potential and density fluctuations - Cross-field current estimation • The density oscillates in-phase with the spoke current • The potential is ~45 out of phase • The azimuthal electric field • The current to the anode: • where d =E/B • The drift current is ~¼ the discharge current, explaining a large fraction of the electron cross-field current to the anode Ellison, Raitses, Fisch, Phys. Plasmas 19 (2012)

  47. Do we know how to explain the spoke instability in Hall thrusters? • A linear stability analysis of the ionization region in HT • An extension of Morozov’s linear analysis for collisionless instability • Spoke appears when the ionization and E-field make it possible to have positive gradients of plasma density and ion velocity Escobar and Ahedo, IEPC 2011 • 3-D Full PIC with MC collisions relate the spoke to neutral depletion Matyash , Schneider et al., IEPC 2011

  48. Potential explanation was given 20 years ago Simon-Hoh instability (SHI) for Penning discharge Conditions for SHI neo 0 neo Er0 > 0 F. C. Hoh, Phys. Fluids 6, 1963 • Modified Simon-Hoh instability (MSHI)- electrostatic instability in a plasma with magnetized electrons and unmagnetized ions due to finite ion Larmor radius effect on azimuthal velocity difference between electrons and ions • Y. Sakawa, C. Josh, P. K. Kaw, F. F. Chen, V. K. Jain, Phys. Fluids B 5, 1993

  49. Can MSHI be excited in CHT plasma ? From the dispersion relation for MSHI, the instability is excited when Azimuthal ion velocity at the location of instability Y. Sakawa et al Phys. Fluids B 5, 1993 From probe measurements of plasma properties and spoke in near- anode region of the Xenon CHT thruster: Br 900 Gauss, Ez 10-20 V/cm, k  1 cm-1, b  30 Not far from our observations

More Related