Sheath stability conditions for PIC simulation
|
|
- Holly Carroll
- 5 years ago
- Views:
Transcription
1 Plasma_conference (2011/11/24,Kanazawa,Japan) PIC Sheath stability conditions for PIC simulation Osaka Prefecture. University H. Matsuura, and N. Inagaki,
2 Background Langmuir probe is the most old but widely-used tool for plasma monitoring. Many improvement work to monitor many parameters are still conducting all over the world. Double probe, triple probe --- density, Te Plasma absorption probe, Frequency shift probe --- ne E missive probe --- Vs Mach probe, directional probe --- flow Ion sensitive probe, tunnel probe --- Ti Combined force-mach probe --- Pressure, Ti Thermal probe --- Heat flux, Ti, negative ion Understanding of sheath is key issue.
3 Schematic view of kinetic/pic model Plasma side Inflection point (ni=ne) Wall side concave(ni<ne) convex(ni>ne) i e 1D in space, stationary. no collision, particle source(sink) potential is zero at plasma boundary and negatively biased externally at wall boundary ref. R.D.Smirnov; Dr.Thesis,
4 Solution of Boltzmann equation 1D steady collisionless Boltzmann equation v x f j (x, v x ) x + q j m j dφ dx f j (x, v x ) v x = 0, can be transferred so that ɛ jx = 1 2 m jv 2 x + q j Φ is used as independent valuable instead of v x. v x (x, ɛ jx ) f j(x, ɛ jx ) x = 0, So energy distribution functions of particle j is the same as those at plasma source boundary.(ex. Maxwellian)
5 Boundary condition Plasma particles( electrons, positive ions, negative ions ) are supplied from boundary. In order to simulate the sheath, choose particle current intensity( or density ) at source boundary so that Debye length is much smaller than geometry size. choose current( or density) ratio of positive and negative particles so that potential profile has at least one inflection point( that is perfect neutral point n+ = n-) in the geometry. At wall boundaries, ions and electrons are perfectly absorbed.(no reflection/no secondary electrons) Initial velocity distribution is Maxwellian.
6 Velocity distribution function Ion f(v x ) = Electron f(v x ) = n is mi 2πT i exp( 1 T i ( 1 2 m iv 2 x + eφ(x))) 0 (v x > 2e(Φ (v x < s Φ(x)) m i ) 2e(Φ s Φ(x)) m i ) 2e(Φ(x) Φ 0 (v x < w ) n es me 2πT e exp( 1 T e ( 1 2 m ev 2 x eφ(x))) (v x > m e ) 2e(Φ(x) Φ w ) m e ) Number density n j = 0 v x v 2 x + 2Z j e(φ s Φ) m j f js(v x )dv x
7 Plasma parameter Ion density n i = n is exp( e e (Φ s Φ))erfc( (Φ s Φ)) T i T i Electron density n e = n es exp( e e (Φ s Φ)){1 + erf( (Φ Φ w ))} T e T e Particle flux Γ i = n is T i 2πm i, Γ e = n es T e 2πm e exp( e T e (Φ s Φ w )) Heat flux Q i = {e(φ s Φ) + 2T i }Γ i Q e = 2T e Γ e
8 Inflection point of potential From normalized Poisson equation d2 φ ds 2 = ρ(s) = (n i n e )/n 0, potential profile is convex in the sheath(ρ(s) > 0). Though, ρ(s) may be equal to be 0 in presheath of real plasma, we set ρ(s) to be slightly positive ( potential is concave ) at plasma source boundary so as that inflection point of potential ( and flat region around it ) exists middle of the simulation geometry. If n is = n es and T i = T e, Normalized potential depth φ = (Φ s Φ)/T e 1.7 at the inflection point. (In the fluid model, φ = 0.5 at presheath end.) And this value does not depend upon mass ration m i /m e and slightly depends upon T i /T e.
9 potential profile calculation Semi-implicit Runge-Kutta method( Kaps-Rentrop-Shapine formula ) is used to solve stiff Poisson equation. T.Watabe et al., "Numerical software by fortran 77",(Maruzen, Tokyo, 1990)[in Japanese]. phiw=10.0d0 eφ w T e denie=1.0d0 ( n i n e ) s tie=1.0d0 ( T i T e ) s xst = 0.0d0 (x/λ D ) s y(1) = 0.0d0 eφ s T e y(2) = 0.7d0 ( eλ D Te dφ dx ) s The effect of "y(2)"(boundary E-field) and "tie"(ion temperature) is mainly studied.
10 Sheath potential from kinetic model Distance is normalized with Debye length.(ti=te) The legends are the normalized electric field(y(2))at source boundary.
11 Condition for sheath establishment System length/debye length System length from kinetic model No sheath E field/(t e /Debye length) at boundary Absolute value of electric field at source boundary must exceed the critical value. The long flat potential is obtain for lower electric field.
12 Berkeley code(xoopic) 2-Dimensional PIC code developed and distributed by PTSG group( Prof. C. K. Birdsall ) J.P.Verboncoueur et al.;comp.phys.comm.,87(1995)
13 parameter limits for sheath simulation For 1D sheath simulation, spatial mesh size x λ D. time step size t 2π/ω pe. Courant condition v e t < x. Sheath stability( derived from generalized Bohm condition ) L x /λ D < 50.
14 1D Simulation model coordinate Cartesian geometry size : L x 4.0, 6.0, 8.0, 10.0[cm] : L y 10.0[cm] mesh size x 150 y 50 time step t [s] particle weight : g Source current : I 0.30, 0.40, 0.55, 0.70[A] Ion He + Boundary condition : Top & bottom Perfect absorbing dielectric : Left Plasma source( and Exit ) with Φ = 0[V] : right conductor with Φ = 100[V]
15 Sheath stability condition 20 z-range=0.04 z-range=0.06 z-range=0.08 z-range=0.10 z-range= current=0.30[a] current=0.40[a] current=0.55[a] current=0.70[a] current=0.70[a] Potential -40 Potential z[m] If system size/debye length becomes too large, the simulation becomes unstable due to plasma oscillation. z[m] If plasma density increase and Debye length becomes small, similar oscillation occurs.
16 Generalized Bohm condition From Poisson equation ɛ d2 Φ dx 2 = e(n e Z i n i ) we obtain ɛ Φ 2 (E2 x Ex0) 2 = e (n e (Φ ) Z i n i (Φ ))dφ Φ 0 As for electrons, Boltzmann relation is assumed here. (Inclusion of the effect of truncation of electron velocity distribution is straightforward, but rather complex.) n e = n es exp( e(φ s Φ) T es ) e Φ Φ 0 n e (Φ )dφ = n es (exp( e(φ s Φ) T es ) exp( e(φ s Φ 0 ) T es ))
17 Generalized Bohm condition(2) For general ion velocity distribution, similar derivation is possible. Z i e Φ Φ 0 n i (Φ )dφ = 0 ( m i u 2 )((1 + 2Z ie(φ s Φ) m i u 2 ) 1 2 (1 + 2Z ie(φ s Φ 0 ) m i u 2 ) 1 2 )fis (u)du By setting x 0 = x s and expanding around here, we obtain, ɛ 2 (E2 x E2 xs ) = (n ese 2 2T es Z2 i n is 2m i < 1 u 2 >)(Φ s Φ) 2 + So generalized Bohm condition is 1 Z it es m i < 1 u 2 > Unfortunately for some kind of distributions( ex. half Maxwellian), f is (u) 0 for u = 0. For such f is (u), < 1 u >= 1 2 n is 0 u f 2 is (u)du will diverge for Φ Φ s and generalized Bohm condition is not satisfied at plasma boundary(x = x s ). 1
18 half-maxwellian ion distribution For some finite values of E xs, there exists an inflection point (x = x n ) and perfect charge neutrality n e = Z i n i is satisfied. Hereafter we use normalized potential depth φ = e(φ s Φ)/T e and normalized ion temperature at plasma boundary τ = T is /T es. ɛ φ 2 (E2 x E2 x0 ) = T es (n es exp( φ ) Z i n is exp( φ φ φ 0 τ )erfc( τ ))dφ = T es n es (exp( φ) exp( φ 0 )) +T es Z i n is τ(exp( φ φ τ )erfc( τ ) + 2 φ π τ exp( φ 0 τ )erfc( φ0 τ ) 2 φ0 π τ ) Then by setting φ 0 = φ n and using n es exp( φ n ) = Z i n is exp( φ n τ )erfc( we obtain ɛ 2 (E2 x E2 xn ) = T esz i n is (( τ ) exp(φ n τ )erfc( φn τ ) 1 τ π φ n τ ) τ φ n )(φ φ n ) 2 +
19 half-maxwellian ion distribution(2) So Bohm condition for half Maxwell ion distribution 10 φn π(τ + 1) τ exp(φ n τ )erfc( φn τ ) 1 In right figure, Green line is the left hand side of above equation plotted as the function of normalized potential depth at neutral point φ n.(τ = 1 and Z i = 1 case ) Another line (Blue ) is Z i n i /n e as as the function of φ. This graph shows φ n, at which blue line cross unity, satisfies above relation. However this may be trivial. There is no reason that so-called sheath boundary where bohm condition is satisfied coincide with an inflection point (x = x n ) (Phi_s-Phi)/Te
20 half-maxwellian ion distribution(3) Instead of expanding with φ φ n, we set φ = φ s. Since E xn = 0 and φ s = 0, we obtain ɛ φn 2 E2 xs = e (n e (Φ ) Z i n i (Φ ))dφ φ s = T es n es (1 exp( φ n ) +T es Z i n is τ(1 exp( φ n τ )erfc( φn By normalized with T es /eλ D, λ 2 D = e2 n es /m e T es, τ ) 2 π φn τ ) E 2 xs (T es /eλ D ) 2 = 2(1 exp( φ n ) + Z in is τ n es (1 exp( φ n τ )erfc( φn τ ) 2 π φn τ ) This value is also plotted in above figure(red line) as the function of φ n. For τ = 1 and Z i = 1 case, φ n 1.7 and E xs 1.3. There values is well agree with kinetic model.
21 2D Simulation model coordinate Cylinder geometry size : L z 4.5, 4.8, 6.0[cm] : L r 1.0[cm] mesh size x 150 y 50 time step t [s] particle weight : g Source current : left 0.55[A] :top [A] Ion He + Boundary condition : Top Plasma source( and Exit ) : bottom exit(not axis) or dielectric : Left Plasma source( and Exit ) with Φ = 0[V] : right conductor or dielectric
22 2D Sheath stability 2D simulation must obey stability condition simultaneously both in r-direction and in z-direction. lz=20[mm] is somehow too long, and potential oscillation and electron energy burst are observed.
23 2D Sheath stability (2) When gap between probe head and left boundary is reduced to lz=5[mm], instability is suppressed. This length corresponds to about 30 times Debye length.
24 Conclusion A kinetic sheath model with Boltzmann - Poisson equations and XOOPIC simulation are compared. Spatial profile of sheath potential, is determined by characteristic electric field strength, which is assigned in artificial boundary condition. Checking validity of modeling ( or simulation )with only potential profile shape is meaningless. In order to keep stable sheath, the necessary value of potential drop is well agree to the estimation with generalized Bohm condition. Simple Bohm condition or fluid model gives us wrong value.
25 Conclusions(cont.) In our simple geometry, the distance between plasma boundary and solid boundary can not exceed over 50 times of Debye length. 2D simulation must obey this condition simultaneously both in r-direction and in z-direction. We acknowledges Prof. Y.Tomita and Prof. D.Tskhakaya for their advic on sheath kinetic model, and Dr.K.Asano for his advice on XOOPIC. We also thank Plasma Theory and Simulation Group in UC-Berkeley for their excellent PIC codes.
26 Appendix
27 Bug of XOOPIC Current version may still contain this bug.
28 Input file correction Input files for old versions must be used carefully.
29 Source file correction physics/vmaxwelln.cpp
30 "simple" Bohm condition If ion velocity distribution is beam like ( f is (u) = n is δ(u V s ), V s > 0 ) and charge neutrality at x = x s is assumed (n es = Z i n is ), n i = n is (1 + 2Z ie(φ s Φ) m i Vs 2 ) 1 2 Z i e Φ n i (Φ )dφ = n is m i Vs 2 ((1+ 2Z ie(φ s Φ) Φ 0 m i V 2 s ) 1 2 (1+ 2Z i e(φ s Φ 0 ) m i V 2 s and, by setting x 0 = x s and expanding around here, we obtain, ɛ 2 (E2 x E2 xs ) = n ese 2 (1 Z it es 2T es m i Vs 2 )(Φ s Φ) 2 + As in stable sheath E 2 x E 2 xs 0, conventional Bohm condition is obtained. ) 1 2 ) V 2 s Z it es m i
31 cutted Maxwellian ion distribution mi f is (u) = Cn is H(u u c ) exp( m i u 2 ), u c > 0 2πT is 2T is Normalization Constant C is determined from n is = 0 f is (u)du m i u C = 2(erfc( 2 c )) 1, erfc(x) = 2 exp( t 2 )dt 2πT is π ion density x n i = Cn is 2 erfc( x c (Φ)) exp( Z ie(φ s Φ) T is ) x c (Φ) = m iu 2 c 2πT is + Z ie(φ s Φ) T is By setting y = Z ie(φ s Φ) T is and y c = m iu 2 c 2πT is and expanding around y s = 0, Z i e Φ n i (Φ )dφ = ( T is ) Cn is Φ s 2 = T is n is { y y s erfc( y c + y ) exp(y )dy y (1 1 πyc exp(y c )erfc( y c ) )y2 + }
32 cutted Maxwellian ion distribution(2) So Bohm condition for cutted ion distribution 10 ɛ 2 (E2 x E2 xs ) = n ese 2 (Φ s Φ) 2 + Z2 i n ise 2 1 (1 2T es 2T is πyc exp(y c )erfc( y c ) )(Φ s Φ) 2 + And, in order to get stable sheath potential, T is Z i T es πyc exp(y c )erfc( y c ) In right figure, right hand side of above equation is plotted as the function of y c = m iu 2 c 2πT is. If T is = T es and Z i = 1, y c > 0.2 will be necessary. And, if T is decreases, y c will increase. (1) mu^2/2ti
33 example of VDF 0.6 vdist.dat using 1:2 vdist.dat using 1:3 0.6 vdist.dat using 1:2 vdist.dat using 1: vdist.dat using 1:2 vdist.dat using 1:3 vdist.dat using 1: Red line: ion, Blue line: green line:electr e(φ s Φ(x)) T e = 0.0, 1.7, 8.0 n is /n es = 1.0, T i /T e = 1.0, m i /m e =
Kinetic Theory of the Presheath and the Bohm Criterion
UW-CPTC 10-3 Kinetic Theory of the Presheath and the Bohm Criterion S D Baalrud and C C Hegna Department of Engineering Physics, University of Wisconsin-Madison, 1500 Engineering Drive, Madison, WI 53706,
More informationOne dimensional hybrid Maxwell-Boltzmann model of shearth evolution
Technical collection One dimensional hybrid Maxwell-Boltzmann model of shearth evolution 27 - Conferences publications P. Sarrailh L. Garrigues G. J. M. Hagelaar J. P. Boeuf G. Sandolache S. Rowe B. Jusselin
More informationPotential Formation in a Plasma Diode Containing Two-electron Temperature Plasma Comparison of Analytical and Numerical Solutions and PIC Simulations
International onference Nuclear Energy for New Europe 006 ortorož, Slovenia, September 18-1, 006 http://www.djs.si/port006 otential Formation in a lasma Diode ontaining Two-electron Temperature lasma omparison
More informationMinimization formulation of a bi-kinetic sheath
Minimization formulation of a bi-kinetic sheath Mehdi Badsi (PHD), B. Després and Martin Campos-Pinto Laboratoire Jacques-Louis Lions UPMC-Paris VI, CNRS UMR 7598, Paris, France Sponsors: ANR Chrome, FRFCM
More informationComputational Methods in Plasma Physics
Computational Methods in Plasma Physics Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Purpose of Talk Describe use of numerical methods to solve simple problem in plasma
More informationA Kinetic Theory of Planar Plasma Sheaths Surrounding Electron Emitting Surfaces
A Kinetic Theory of Planar Plasma Sheaths Surrounding Electron Emitting Surfaces J. P. Sheehan1, I. Kaganovich2, E. Barnat3, B. Weatherford3, H. Wang2, 4 1 2 D. Sydorenko, N. Hershkowitz, and Y. Raitses
More information13.1 Ion Acoustic Soliton and Shock Wave
13 Nonlinear Waves In linear theory, the wave amplitude is assumed to be sufficiently small to ignore contributions of terms of second order and higher (ie, nonlinear terms) in wave amplitude In such a
More informationSpin Stability of Aysmmetrically Charged Plasma Dust. I.H. Hutchinson. February 2004
PSFC/JA-04-3 Spin Stability of Aysmmetrically Charged Plasma Dust I.H. Hutchinson February 2004 Plasma Science and Fusion Center Massachusetts Institute of Technology Cambridge, MA 02139 USA This work
More informationPIC Algorithm with Multiple Poisson Equation Solves During One Time Step
Journal of Physics: Conference Series PAPER OPEN ACCESS PIC Algorithm with Multiple Poisson Equation Solves During One Time Step To cite this article: Junxue Ren et al 2015 J. Phys.: Conf. Ser. 640 012033
More informationSimulation of a two-dimensional sheath over a flat wall with an insulatorõconductor interface exposed to a high density plasma
JOURNAL OF APPLIED PHYSICS VOLUME 94, NUMBER 5 1 SEPTEMBER 2003 Simulation of a two-dimensional sheath over a flat wall with an insulatorõconductor interface exposed to a high density plasma Doosik Kim
More informationA comparison of emissive probe techniques for electric potential measurements in a Hall thruster plasma
A comparison of emissive probe techniques for electric potential measurements in a Hall thruster plasma J. P. Sheehan*, Y. Raitses**, N. Hershkowitz*, I. Kaganovich**, and N. J. Fisch** *University of
More informationLow Temperature Plasma Technology Laboratory
Low Temperature Plasma Technology Laboratory Equilibrium theory for plasma discharges of finite length Francis F. Chen and Davide Curreli LTP-6 June, Electrical Engineering Department Los Angeles, California
More informationAn introduction to Langmuir probe diagnostics of plasmas
An introduction to Langmuir probe diagnostics of plasmas Luis Conde Departamento de Física Aplicada E.T.S. Ingenieros Aeronáuticos 284 Madrid, Spain. May 28, 211 Abstract In this short review is introduced
More informationEnergy-Conserving Numerical Simulations of Electron Holes in Two-Species Plasmas
Energy-Conserving Numerical Simulations of Electron Holes in Two-Species Plasmas Yingda Cheng Andrew J. Christlieb Xinghui Zhong March 18, 2014 Abstract In this paper, we apply our recently developed energy-conserving
More informationDevelopment of a Hall Thruster Fully Kinetic Simulation Model Using Artificial Electron Mass
Development of a Hall Thruster Fully Kinetic Simulation Model Using Artificial Electron Mass IEPC-013-178 Presented at the 33rd International Electric Propulsion Conference, The George Washington University
More informationA note on the plasma sheath and the Bohm Criterion
A note on the plasma sheath and the Bohm Criterion G.D. Severn Dept. of Physics, University of San Diego, San Diego CA 92110 (Dated: April 6, 2006) PACS numbers: 52.27.Aj, 52.27.Cm The word sheath in connection
More informationAnalytic Expression For the Electric Potential in the Plasma Sheath
884 IEEE TRANSACTIONS ON PLASMA SCIENCE. VOL. 17. NO. 6. DECEMBER 1989 Analytic Expression For the Electric Potential in the Plasma Sheath TERRENCE E. SHERIDAN, JR. AND JOHN A GOREE Abstract-An expression
More informationAccurately Determining the Plasma Potential Using Emissive Probes
Accurately Determining the Plasma Potential Using Emissive Probes IEPC-2013-313 Presented at the 33 rd International Electric Propulsion Conference, The George Washington University, Washington, D.C.,
More informationPlasma-Wall Interaction: Sheath and Pre-sheath
Plasa-Wall Interaction: Sheath and Pre-sheath Under ost conditions, a very thin negative sheath appears in the vicinity of walls, due to accuulation of electrons on the wall. This is in turn necessitated
More informationHigh Order Semi-Lagrangian WENO scheme for Vlasov Equations
High Order WENO scheme for Equations Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Andrew Christlieb Supported by AFOSR. Computational Mathematics Seminar, UC Boulder
More informationBoundary Conditions for the Child Langmuir Sheath Model
IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 28, NO. 6, DECEMBER 2000 2207 Boundary Conditions for the Child Langmuir Sheath Model Mikhail S. Benilov Abstract A collision-free space-charge sheath formed by
More informationThe Invalidity of a Mach Probe Model
The Invalidity of a Mach Probe Model I. H. Hutchinson Plasma Science and Fusion Center Massachusetts Institute of Technology Cambridge, MA, USA Abstract Despite its recent application to the interpretation
More informationParticle-In-Cell Simulations of a Current-Free Double Layer
Particle-In-Cell Simulations of a Current-Free Double Layer S. D. Baalrud 1, T. Lafleur, C. Charles and R. W. Boswell American Physical Society Division of Plasma Physics Meeting November 10, 2010 1Present
More informationFigure 1.1: Ionization and Recombination
Chapter 1 Introduction 1.1 What is a Plasma? 1.1.1 An ionized gas A plasma is a gas in which an important fraction of the atoms is ionized, so that the electrons and ions are separately free. When does
More informationPositive ion flux from a low-pressure electronegative discharge
Plasma Sources Sci. Technol. 8 (1999) 457 462. Printed in the UK PII: S0963-0252(99)03994-8 Positive ion flux from a low-pressure electronegative discharge T E Sheridan, P Chabert and R W Boswell Space
More informationModeling Electron Characteristics in an Ion Thruster Plume: Fully Kinetic PIC vs. Hybrid PIC
Modeling Electron Characteristics in an Ion Thruster Plume: Fully Kinetic PIC vs. Hybrid PIC IEPC-2017-301 Presented at the 35th International Electric Propulsion Conference Georgia Institute of Technology,
More informationSheaths: More complicated than you think a
PHYSICS OF PLASMAS 12, 055502 2005 Sheaths: More complicated than you think a Noah Hershkowitz b University of Wisconsin-Madison, Madison, Wisconsin 53706 Received 7 December 2004; accepted 7 February
More informationThe ideal Maxwellian plasma
The ideal Maxwellian plasma Dr. L. Conde Departamento de Física Aplicada. E.T.S. Ingenieros Aeronáuticos Universidad Politécnica de Madrid Plasmas are,... The plasma state of matter may be defined as a
More information1358 IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 27, NO. 5, OCTOBER 1999
1358 IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 27, NO. 5, OCTOBER 1999 Sheath Thickness Evaluation for Collisionless or Weakly Collisional Bounded Plasmas Shiang-Bau Wang and Amy E. Wendt Abstract The
More informationFundamentals of Plasma Physics Transport in weakly ionized plasmas
Fundamentals of Plasma Physics Transport in weakly ionized plasmas APPLAuSE Instituto Superior Técnico Instituto de Plasmas e Fusão Nuclear Luís L Alves (based on Vasco Guerra s original slides) 1 As perguntas
More informationKinetic Simulations of Ion Beam Neutralization
Kinetic Simulations of Ion Beam Neutralization O. Chang and J. Wang Astronautical Engineering Department University of Southern California Los Angeles, CA 90089-1192, USA Abstract. Full particle PIC simulations
More informationNumerical Simulation of Faraday Probe Measurements in a Multi-component Non-equilibrium Plasma
Numerical Simulation of Faraday Probe Measurements in a Multi-component Non-equilibrium Plasma IEPC-005-85 Presented at the 9 th International Electric Propulsion Conference, Princeton University, Jeremiah
More informationSummer College on Plasma Physics. 30 July - 24 August, The particle-in-cell simulation method: Concept and limitations
1856-30 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 The particle-in-cell M. E. Dieckmann Institut fuer Theoretische Physik IV, Ruhr-Universitaet, Bochum, Germany The particle-in-cell
More informationSensors Plasma Diagnostics
Sensors Plasma Diagnostics Ken Gentle Physics Department Kenneth Gentle RLM 12.330 k.gentle@mail.utexas.edu NRL Formulary MIT Formulary www.psfc.mit.edu/library1/catalog/ reports/2010/11rr/11rr013/11rr013_full.pdf
More informationNegative plasma potential relative to electron-emitting surfaces
PHYSICAL REVIEW E 88, 033103 (2013) Negative plasma potential relative to electron-emitting surfaces M. D. Campanell * Princeton Plasma Physics Laboratory, Princeton University, Princeton, New Jersey 08543,
More informationLandau Damping Simulation Models
Landau Damping Simulation Models Hua-sheng XIE (u) huashengxie@gmail.com) Department of Physics, Institute for Fusion Theory and Simulation, Zhejiang University, Hangzhou 310027, P.R.China Oct. 9, 2013
More informationLow Temperature Plasma Technology Laboratory
Low Temperature Plasma Technology Laboratory CENTRAL PEAKING OF MAGNETIZED GAS DISCHARGES Francis F. Chen and Davide Curreli LTP-1210 Oct. 2012 Electrical Engineering Department Los Angeles, California
More informationPhysics and Modelling of a Negative Ion Source Prototype for the ITER Neutral Beam Injection
1 ITR/P1-37 Physics and Modelling of a Negative Ion Source Prototype for the ITER Neutral Beam Injection J.P. Boeuf a, G. Fubiani a, G. Hagelaar a, N. Kohen a, L. Pitchford a, P. Sarrailh a, and A. Simonin
More informationKinetic simulation of the stationary HEMP thruster including the near field plume region
Kinetic simulation of the stationary HEMP thruster including the near field plume region IEPC-2009-110 Presented at the 31st International Electric Propulsion Conference, University of Michigan Ann Arbor,
More informationSpacecraft Environment Interaction Engineering
Spacecraft Environment Interaction Engineering Spacecraft Charging Analysis Mengu Cho Laboratory of Spacecraft Environment Interaction Engineering Kyushu Institute of Technology cho@ele.kyutech.ac.jp http://laseine.ele.kyutech.ac.jp
More informationKinetic Theory of Instability-Enhanced Collisions and Its Application to Langmuir s Paradox and the Multi-Species Bohm Criterion
Kinetic Theory of Instability-Enhanced Collisions and Its Application to Langmuir s Paradox and the Multi-Species Bohm Criterion Scott D. Baalrud in collaboration with Chris C. Hegna and James D. Callen
More informationQuasi-neutral limit for Euler-Poisson system in the presence of plasma sheaths
in the presence of plasma sheaths Department of Mathematical Sciences Ulsan National Institute of Science and Technology (UNIST) joint work with Masahiro Suzuki (Nagoya) and Chang-Yeol Jung (Ulsan) The
More informationFundamentals of Plasma Physics
Fundamentals of Plasma Physics Definition of Plasma: A gas with an ionized fraction (n i + + e ). Depending on density, E and B fields, there can be many regimes. Collisions and the Mean Free Path (mfp)
More informationSolution of time-dependent Boltzmann equation for electrons in non-thermal plasma
Solution of time-dependent Boltzmann equation for electrons in non-thermal plasma Z. Bonaventura, D. Trunec Department of Physical Electronics Faculty of Science Masaryk University Kotlářská 2, Brno, CZ-61137,
More informationAPPENDIX Z. USEFUL FORMULAS 1. Appendix Z. Useful Formulas. DRAFT 13:41 June 30, 2006 c J.D Callen, Fundamentals of Plasma Physics
APPENDIX Z. USEFUL FORMULAS 1 Appendix Z Useful Formulas APPENDIX Z. USEFUL FORMULAS 2 Key Vector Relations A B = B A, A B = B A, A A = 0, A B C) = A B) C A B C) = B A C) C A B), bac-cab rule A B) C D)
More informationPhysique des plasmas radiofréquence Pascal Chabert
Physique des plasmas radiofréquence Pascal Chabert LPP, Ecole Polytechnique pascal.chabert@lpp.polytechnique.fr Planning trois cours : Lundi 30 Janvier: Rappels de physique des plasmas froids Lundi 6 Février:
More informationModeling neutral-plasma interactions in scrape-off layer (SOLT) simulations*
Modeling neutral-plasma interactions in scrape-off layer (SOLT) simulations* D. A. Russell and J. R. Myra Research Corporation Boulder CO USA Presented at the US Transport Task Force Workshop Williamsburg
More informationHeating and current drive: Radio Frequency
Heating and current drive: Radio Frequency Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 13 th February 2012 Dr Ben Dudson Magnetic Confinement Fusion (1 of 26)
More informationPRINCIPLES OF PLASMA DISCHARGES AND MATERIALS PROCESSING
PRINCIPLES OF PLASMA DISCHARGES AND MATERIALS PROCESSING Second Edition MICHAEL A. LIEBERMAN ALLAN J, LICHTENBERG WILEY- INTERSCIENCE A JOHN WILEY & SONS, INC PUBLICATION CONTENTS PREFACE xrrii PREFACE
More informationFluid Neutral Momentum Transport Reference Problem D. P. Stotler, PPPL S. I. Krasheninnikov, UCSD
Fluid Neutral Momentum Transport Reference Problem D. P. Stotler, PPPL S. I. Krasheninnikov, UCSD 1 Summary Type of problem: kinetic or fluid neutral transport Physics or algorithm stressed: thermal force
More informationLecture 2. Introduction to plasma physics. Dr. Ashutosh Sharma
Preparation of the concerned sectors for educational and R&D activities related to the Hungarian ELI project Ion acceleration in plasmas Lecture 2. Introduction to plasma physics Dr. Ashutosh Sharma Zoltán
More informationWaves in plasma. Denis Gialis
Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.
More informationKinetic Simulation of Effects of Secondary Electron Emission on Electron Temperature in Hall Thrusters
Kinetic Simulation of Effects of Secondary Electron Emission on Electron Temperature in Hall Thrusters IEPC-25-78 Presented at the 29 th International Electric Propulsion Conference, Princeton University
More informationMagnetized ion collection by oblique surfaces including self-consistent drifts: Mach-probes of arbitrary shape.
1 Magnetized ion collection by oblique surfaces including self-consistent drifts: Mach-probes of arbitrary shape I H Hutchinson Plasma Science and Fusion Center and and Engineering Department MIT APS DPP
More informationCharacteristics of Positive Ions in the Sheath Region of Magnetized Collisional Electronegative Discharges
Plasma Science and Technology, Vol.6, No.6, Jun. 204 Characteristics of Positive Ions in the Sheath Region of Magnetized Collisional Electronegative Discharges M. M. HATAMI, A. R. NIKNAM 2 Physics Department
More informationTurbulent Transport due to Kinetic Ballooning Modes in High-Beta Toroidal Plasmas
1 TH/P-3 Turbulent Transport due to Kinetic allooning Modes in High-eta Toroidal Plasmas A. Ishizawa 1, S. Maeyama, T.-H. Watanabe 1, H. Sugama 1 and N. Nakajima 1 1 National Institute for Fusion Science,
More informationEXPERIMENTAL STUDIES ON ELECTRO- POSITIVE/NEGATIVE PLASMAS PRODUCED BY RF DISCHARGES WITH EXTERNAL ANTENNA IN LOW-PRESSURE GASES
EXPERIMENTAL STUDIES ON ELECTRO- POSITIVE/NEGATIVE PLASMAS PRODUCED BY RF DISCHARGES WITH EXTERNAL ANTENNA IN LOW-PRESSURE GASES SEPTEMBER 2004 DEPARTMENT OF ENERGY AND MATERIALS SCIENCE GRADUATE SCHOOL
More informationApplication of a Laser Induced Fluorescence Model to the Numerical Simulation of Detonation Waves in Hydrogen-Oxygen-Diluent Mixtures
Supplemental material for paper published in the International J of Hydrogen Energy, Vol. 30, 6044-6060, 2014. http://dx.doi.org/10.1016/j.ijhydene.2014.01.182 Application of a Laser Induced Fluorescence
More informationA SHORT COURSE ON THE PRINCIPLES OF PLASMA DISCHARGES AND MATERIALS PROCESSING
A SHORT COURSE ON THE PRINCIPLES OF DISCHARGES AND MATERIALS PROCESSING Michael A. Lieberman Department of Electrical Engineering and Computer Sciences, CA 94720 LiebermanShortCourse15 TABLE OF CONTENTS
More informationModelling of plasma tank and related langmuir probe calibration MATEO-VELEZ J.-C, ROUSSEL J.-F., SARRAIL D, BOULAY F., INGUIMBERT V.
Modelling of plasma tank and related langmuir probe calibration MATEO-VELEZ J.-C, OUSSEL J.-F., SAAIL D, BOULAY F., INGUIMBET V. PAYAN D. ONEA CNES Objectives Initial: Validation of SPIS modelling (LEO
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationPlasma-Wall Interaction Controlled by Secondary Electron Emission
Plasma-Wall Interaction Controlled by Secondary Electron Emission IEPC-0-/ISTS-0-b- Presented at Joint Conference of 0th International Symposium on Space Technology and Science, th International Electric
More informationLecture Note 1. 99% of the matter in the universe is in the plasma state. Solid -> liquid -> Gas -> Plasma (The fourth state of matter)
Lecture Note 1 1.1 Plasma 99% of the matter in the universe is in the plasma state. Solid -> liquid -> Gas -> Plasma (The fourth state of matter) Recall: Concept of Temperature A gas in thermal equilibrium
More informationEP2Plus: a hybrid plasma. plume/spacecraft. interaction code. F. Cichocki, M. Merino, E. Ahedo
EP2Plus: a hybrid plasma plume/spacecraft interaction code F. Cichocki, M. Merino, E. Ahedo 24 th SPINE meeting ESTEC, Noordwijk, October 23 rd, 2017 Contents Introduction to EP2PLUS Overall structure
More informationWhat happens to ions at the plasma boundary in multiple-ion species plasmas?
What happens to ions at the plasma boundary in multiple-ion species plasmas? How diode lasers help illumine the problem of sheath formation Greg Severn Plasma Theory and Simulation Group Seminar, University
More informationParallel transport and profile of boundary plasma with a low recycling wall
1 TH/P4-16 Parallel transport and profile of boundary plasma with a low recycling wall Xian-Zhu Tang 1 and Zehua Guo 1 1 Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545, U.S.A.
More informationModeling of planar plasma diode
International Letters of Chemistry, Physics and Astronomy Online: 2013-05-03 ISSN: 2299-3843, Vol. 13, pp 220-242 doi:10.18052/www.scipress.com/ilcpa.13.220 2013 SciPress Ltd., Switzerland Modeling of
More informationTwo-dimensional Particle-In-Cell model of the extraction region of the PEGASES ion-ion plasma source
Two-dimensional Particle-In-Cell model of the extraction region of the PEGASES ion-ion plasma source IEPC-2013-249 Presented at the 33rdInternational Electric Propulsion Conference, The George Washington
More informationSimple examples of MHD equilibria
Department of Physics Seminar. grade: Nuclear engineering Simple examples of MHD equilibria Author: Ingrid Vavtar Mentor: prof. ddr. Tomaž Gyergyek Ljubljana, 017 Summary: In this seminar paper I will
More informationFeature-level Compensation & Control
Feature-level Compensation & Control 2 Plasma Eray Aydil, UCSB, Mike Lieberman, UCB and David Graves UCB Workshop November 19, 2003 Berkeley, CA 3 Feature Profile Evolution Simulation Eray S. Aydil University
More informationKinetic theory of ions in the magnetic presheath
Kinetic theory of ions in the magnetic presheath Alessandro Geraldini 1,2, Felix I. Parra 1,2, Fulvio Militello 2 1. Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, Oxford
More informationEquilibrium model for two low-pressure electronegative plasmas connected by a double layer
PHYSICS OF PLASMAS 13, 093504 2006 Equilibrium model for two low-pressure electronegative plasmas connected by a double layer P. Chabert, a N. Plihon, C. S. Corr, and J.-L. Raimbault Laboratoire de Physique
More informationPlasma Astrophysics Chapter 1: Basic Concepts of Plasma. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University
Plasma Astrophysics Chapter 1: Basic Concepts of Plasma Yosuke Mizuno Institute of Astronomy National Tsing-Hua University What is a Plasma? A plasma is a quasi-neutral gas consisting of positive and negative
More informationChapter 1 Nature of Plasma
Chapter 1 Nature of Plasma Abstract Charge neutrality is one of fundamental property of plasma. Section 1.2 explains Debye length λ D in (1.2), a measure of shielding distance of electrostatic potential,
More informationTheory and simulation of electron sheaths and anode spots in low pressure laboratory plasmas
University of Iowa Iowa Research Online Theses and Dissertations Summer 217 Theory and simulation of electron sheaths and anode spots in low pressure laboratory plasmas Brett Stanford Scheiner University
More informationStructure of Velocity Distribution of Sheath-Accelerated Secondary Electrons in Asymmetric RF-DC Discharge
Structure of Velocity Distribution of Sheath-Accelerated Secondary Electrons in Asymmetric RF-DC Discharge Alexander V. Khrabrov 1, Igor D. Kaganovich 1, Peter L. G. Ventzek 2, Alok Ranjan 3, and Lee Chen
More informationA floating potential method for measuring ion density
A floating potential method for measuring ion density Francis F. Chen, John D. Evans, and Donald Arnush Electrical Engineering Department, University of California Los Angeles Los Angeles, California 995-1594
More informationPlasma Physics Prof. V. K. Tripathi Department of Physics Indian Institute of Technology, Delhi
Plasma Physics Prof. V. K. Tripathi Department of Physics Indian Institute of Technology, Delhi Lecture No. # 03 DC Conductivity and Negative Differential Conductivity Well friends, in this lecture, I
More informationNon-Equidistant Particle-In-Cell for Ion Thruster Plumes
Non-Equidistant Particle-In-Cell for Ion Thruster Plumes IEPC-213-67 Presented at the 33 rd International Electric Propulsion Conference, The George Washington University, Washington, D.C., USA October
More informationA particle-in-cell method with adaptive phase-space remapping for kinetic plasmas
A particle-in-cell method with adaptive phase-space remapping for kinetic plasmas Bei Wang 1 Greg Miller 2 Phil Colella 3 1 Princeton Institute of Computational Science and Engineering Princeton University
More informationExperimental Studies of Ion Beam Neutralization: Preliminary Results
Experimental Studies of Ion Beam Neutralization: Preliminary Results N. Ding, J. Polansky, R. Downey and J. Wang Department of Astronautical Engineering University of Southern California Los Angeles, CA
More informationIntroduction. Chapter Plasma: definitions
Chapter 1 Introduction 1.1 Plasma: definitions A plasma is a quasi-neutral gas of charged and neutral particles which exhibits collective behaviour. An equivalent, alternative definition: A plasma is a
More informationProgress on Quantitative Modeling of rf Sheaths
Progress on Quantitative Modeling of rf Sheaths D. A. D Ippolito, J. R. Myra, H. Kohno and J. C. Wright Lodestar Research Corporation, Boulder, Colorado, 80301 May, 2011 Prepared for the 19th Topical Conference
More informationxkcd.com It IS about physics. It ALL is.
xkcd.com It IS about physics. It ALL is. Introduction to Space Plasmas The Plasma State What is a plasma? Basic plasma properties: Qualitative & Quantitative Examples of plasmas Single particle motion
More informationPHYSICS Computational Plasma Physics
PHYSICS 78 - Computational Plasma Physics INSTRUCTOR Dr. Earl Scime (escime@wvu.edu) 93-34, ext. 1437 Office hours: MW :30 3:30 and whenever door is open Rm 18 & 05 Hodges Hall Class: MWF 1:30-:0 Rm 334
More informationKinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles
Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3. The 8th International Conference on Computational
More informationTheory of dust voids in plasmas
PHYSICAL REVIEW E VOLUME 59, NUMBER 6 JUNE 1999 Theory of dust voids in plasmas J. Goree,* G. E. Morfill, V. N. Tsytovich, and S. V. Vladimirov Max Planck Institut für Extraterrestrische Physik, Giessenbachstrasse,
More informationContents: 1) IEC and Helicon 2) What is HIIPER? 3) Analysis of Helicon 4) Coupling of the Helicon and the IEC 5) Conclusions 6) Acknowledgments
Contents: 1) IEC and Helicon 2) What is HIIPER? 3) Analysis of Helicon 4) Coupling of the Helicon and the IEC 5) Conclusions 6) Acknowledgments IEC:! IEC at UIUC modified into a space thruster.! IEC has
More informationSIMULATIONS OF ECR PROCESSING SYSTEMS SUSTAINED BY AZIMUTHAL MICROWAVE TE(0,n) MODES*
25th IEEE International Conference on Plasma Science Raleigh, North Carolina June 1-4, 1998 SIMULATIONS OF ECR PROCESSING SYSTEMS SUSTAINED BY AZIMUTHAL MICROWAVE TE(,n) MODES* Ron L. Kinder and Mark J.
More informationEquilibrium Properties of Matter and Radiation
Equilibrium Properties of Matter and Radiation Temperature What is it? A measure of internal energy in a system. Measure from (1) velocities of atoms/molecules () population of excited/ionized states (3)
More informationSimulation of a two-dimensional sheath over a flat insulator conductor interface on a radio-frequency biased electrode in a high-density plasma
JOURNAL OF APPLIED PHYSICS VOLUME 95, NUMBER 7 1 APRIL 2004 Simulation of a two-dimensional sheath over a flat insulator conductor interface on a radio-frequency biased electrode in a high-density plasma
More informationKinetic Simulations of Plasma Plume Potential in a Vacuum Chamber
Missouri University of Science and Technology Scholars' Mine Mechanical and Aerospace Engineering Faculty Research & Creative Works Mechanical and Aerospace Engineering 10-1-2013 Kinetic Simulations of
More informationPIC-MCC simulations for complex plasmas
GRADUATE SUMMER INSTITUTE "Complex Plasmas August 4, 008 PIC-MCC simulations for complex plasmas Irina Schweigert Institute of Theoretical and Applied Mechanics, SB RAS, Novosibirsk Outline GRADUATE SUMMER
More informationCollisional sheath dynamics in the intermediate rf frequency regime
Collisional sheath dynamics in the intermediate rf frequency regime N.Xiang, F.L.Waelbroeck Institute for Fusion Studies, University of Texas Austin,Texas 7871 email: nongx@mail.utexas.edu April 16, 4
More informationImpact of neutral atoms on plasma turbulence in the tokamak edge region
Impact of neutral atoms on plasma turbulence in the tokamak edge region C. Wersal P. Ricci, F.D. Halpern, R. Jorge, J. Morales, P. Paruta, F. Riva Theory of Fusion Plasmas Joint Varenna-Lausanne International
More informationarxiv: v1 [physics.plasm-ph] 16 Apr 2018
Consistent simulation of capacitive radio-frequency discharges and external matching networks Frederik Schmidt Institute of Theoretical Electrical Engineering, arxiv:1804.05638v1 [physics.plasm-ph] 16
More informationMHD RELATED TO 2-FLUID THEORY, KINETIC THEORY AND MAGANETIC RECONNECTION
MHD RELATED TO 2-FLUID THEORY, KINETIC THEORY AND MAGANETIC RECONNECTION Marty Goldman University of Colorado Spring 2017 Physics 5150 Issues 2 How is MHD related to 2-fluid theory Level of MHD depends
More informationNonlinear Diffusion in Magnetized Discharges. Francis F. Chen. Electrical Engineering Department
Nonlinear Diffusion in Magnetized Discharges Francis F. Chen Electrical Engineering Department PPG-1579 January, 1998 Revised April, 1998 Nonlinear Diffusion in Magnetized Discharges Francis F. Chen Electrical
More informationA Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1
A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1 B. Srinivasan 2, U. Shumlak Aerospace and Energetics Research Program University of Washington,
More information