Potential Formation in a Plasma Diode Containing Two-electron Temperature Plasma Comparison of Analytical and Numerical Solutions and PIC Simulations
|
|
- Emily Henderson
- 6 years ago
- Views:
Transcription
1 International onference Nuclear Energy for New Europe 006 ortorož, Slovenia, September 18-1, otential Formation in a lasma Diode ontaining Two-electron Temperature lasma omparison of Analytical and Numerical Solutions and I Simulations Tomaž Gyergyek 1,, Milan Čerček,3, Borut Jurčič-Zlobec 1 1) University of Ljubljana, Faculty of electrical engineering, Tržaška 5, SI-1000 Ljubljana, Slovenia ) Jožef Stefan Institute Jamova 39, SI-1000 Ljubljana, Slovenia 3) University of Maribor, Faculty of civil engineering, Smetanova 17, SI-000 Maribor, Slovenia Tomaz.gyergyek@fe.uni-lj.si, Milan.ercek@ijs.si ABSTRAT otential formation in front of a negatively biased electron emitting electrode (collector) is studied by a fully kinetic one dimensional model of a plasma diode that contains electrons with a two-temperature maxwellian velocity distribution function. Analytical and numerical results of the model are compared to the results of the I simulation. 1 INTRODUTION In the edge plasmas of fusion devices energetic electron populations are readily produced. The presence of energetic electrons has a remarkable effect on the potential formation in the plasma and consequently on particle losses to the wall. Also electron emission from the walls that limit the plasma is a common occurrence in fusion devices. In this work we present a fully kinetic model of plasma potential formation in a bounded plasma system containing energetic and emitted electrons and compare the results of the model with numerical solutions of the oisson equation and with the computer simulation. MODEL We consider a large planar electrode (collector) with its surface perpendicular to the x axis of the coordinate system. The collector is located at x = 0. This electrode absorbs all the particles that hit it. On the other hand it may also emit electrons. This electron emission may be thermal or secondary, triggered by the impact of incoming electrons and/or ions. The details of the emission mechanism are not essential for the model. An infinitely large planar plasma source also has its surface perpendicular to the x axis. The source is located at a certain distance x = L from the collector. The distance L is not essential for the model and may be taken arbitrarily. The source injects 3 groups of charged particles into the system: singly charged positive ions (index i), the cool electrons (index 1), the hot electrons (index ). The electrons that are emitted from the collector have index 3. The particles i, 1 and are injected from the source with half-maxwellian velocity distribution function with respective temperatures T i, T 1 and T. The emitted electrons also have a halfmaxwellian velocity distribution function at the collector with the temperature T 3. We assume 506.1
2 506. that T > T 1 and T 1 >> T 3. The ion temperature T i can in principle be arbitrary, but is usually taken smaller than or equal to T 1. The collector is biased to a certain (negative!) potential Φ. The potential an the electric field at the source are set to zero. This imposes the following boundary conditions: dφ Φ ( x= L) = 0, ( x= L) = 0, (1) dx for the oisson equation: d Φ e = 0 ( ni ( x) n1( x) n( x) n3( x) ). () dx ε 0 where e 0 is the elementary charge and n j (x) are particle densities. We are looking only for such solutions of the oisson equation (), where Φ(x) is a monotonically decreasing function in the direction towards the collector. An ion that has left the source with a negligibly small initial velocity has at the distance x < L from the collector the velocity: e0φ( x) vmi =. (3) m The distribution function for the ions can therefore be written in the following way: m i e0 Φ( x) mv i fi = ni exp exp H( v vmi ). (4) π kti kti kti Here n i is the ion density at the source. Velocity direction towards the collector is positive and away from the collector is negative. An electron that has left the collector source with a negligibly small initial velocity has at the distance x from the collector the velocity: v i ( Φ x Φ ) e ( ) 0 me =. (5) me So the distribution functions for the electrons are written in the following way: m e e0φ( x) mev f1 = n1 exp exp H( v vme), (6) π kt1 kt1 kt1 m e e0φ( x) mev f = n exp exp H( v vme), (7) π kt kt kt m e0 ( Φ( x) Φ ) e m 3 3 exp exp e v f = n H( vme v). (8) π kt3 kt3 kt3 Here n 1, n, and n 3 are the respective electron densities at the source and at the collector. We now introduce the following variables: me Ti T T3 e0φ( x) e0φ ( x= 0) e0φ μ =, τ =, Θ=, σ =, Ψ =, Ψ = =, mi T1 T1 T1 kt1 kt1 kt1 (9) ni n n3 kt1 v x ε 0kT α =, β =, ε =, v0 =, u =, z =, λ 1 D =. n n n m v λ ne e 0 D 1 0 With these variables the distribution functions are written in the following way: α Ψ u Fi ( u, Ψ ) = exp exp H( u μψ ), (10) πτμ τ μτ roceedings of the International onference Nuclear Energy for New Europe, 006
3 (, ) exp ( ) exp( ) ( ) F1 u Ψ = Ψ u H u+ Ψ Ψ, (11) π β Ψ u ( ) ( F u, Ψ = exp exp H u+ Ψ Ψ ), (1) π Θ Θ Θ ε Ψ Ψ u ( u ) ( F3, Ψ = exp exp H Ψ Ψ u. (13) πσ σ σ The n-th moment of the velocity distribution is obtained in the following way. The distribution function is multiplied by the selected power of the velocity u n and then integrated over the velocity space. In this way the moment of the distribution function is obtained as a function of the potential Ψ. The zero moments (n = 0) are particle densities and the first moments (n = 1) are fluxes. For the distribution functions (10) - (13) the first and zero moments are the following: α Ψ Ψ α μτ Ni( Ψ ) = exp Erfc, Ji( Ψ ) =, τ τ π 1 1 N1( Ψ ) = exp( Ψ )( 1+ Erf Ψ Ψ ), J1( Ψ ) = exp ( Ψ), π (14) β Ψ Ψ Ψ β Θ Ψ N( Ψ ) = exp 1+ Erf, J( Ψ ) = exp, Θ Θ π Θ ε Ψ Ψ Ψ Ψ ε σ N3 ( Ψ ) = exp Erfc σ σ, J3 ( Ψ ) =, π where ( ) x Erf x = exp ( t ) dt, Erfc( x) = exp ( t ) dt. (15) π 0 π x The total current density J t to the collector is then given by: α μτ ε σ 1 β Θ Ψ Jt = Ji+ J3 J1 J = + exp( Ψ ) exp. π π π π Θ (16) With the variables (9) and particle densities given by (14) the oisson equation () is written in the following form: d Ψ Ψ( z) Ψ( z) = α exp Erfc + dz τ τ ( z) ( z) ( z) ( z) +exp( ( z) ) 1 Erf ( ( z) Ψ Ψ Ψ Ψ + Ψ Ψ ) + β exp 1 + Erf + Θ Θ Ψ Ψ Ψ Ψ + εexp Erfc. σ σ At a certain distance z = z 0 from the collector the plasma is neutral and the potential at that point has the value Ψ(z = z 0 ) = Ψ. So the neutrality condition can be obtained directly from (17): ) (17) roceedings of the International onference Nuclear Energy for New Europe, 006
4 z= z d Ψ Ψ Ψ α exp Erfc = + dz τ τ Ψ Ψ Ψ +exp( Ψ ) 1 + Erf ( Ψ Ψ ) + β exp 1 + Erf + Θ Θ Ψ Ψ Ψ Ψ + εexp Erfc = 0. σ σ At z = z 0 the potential Ψ(z) obviously has an inflection point. Since: 1 d dψ dψ d Ψ =, dz dz dz dz the oisson equation (17) is multiplied by dψ/dz and integrated once over Ψ from Ψ = 0 (at the source) to Ψ = Ψ (at the inflection point). The following equation is obtained: Ψ Ψ Ψ ατ 1 exp Erfc + π τ τ τ ( ) ( )( ) +exp Ψ 1 Erf exp 1 Erf + Ψ Ψ Ψ Ψ Ψ Ψ + Ψ + π Ψ Ψ Ψ Ψ Ψ Ψ Ψ Ψ + βθ exp 1+ Erf exp 1+ Erf + Θ Θ π Θ Θ Θ Θ Ψ Ψ Ψ Ψ Ψ Ψ Ψ Ψ Ψ + εσ exp Erfc exp Erfc 0. π + σ σ = σ σ σ σ (19) This equation is called the zero electric field condition at the inflection point. If the electron emission from the collector increases, at a certain level of emission a negative space charge may start to accumulate at the collector surface. At certain critical electron emission the electric field that accelerates the negative electrons away from the collector becomes equal to zero. The zero electric field condition at the collector can be derived in a very similar way as the zero electric field condition at the inflection point (19). Again the oisson equation (17) is multiplied by dψ/dz and integrated once over Ψ, only now the integration limits go from Ψ = Ψ (at the inflection point) to Ψ = Ψ (at the collector). The following equation is obtained: Ψ Ψ Ψ Ψ Ψ Ψ ατ exp Erfc exp Erfc + τ τ τ τ π τ τ +exp( Ψ ) 1+ Ψ Ψ exp( Ψ ) 1 + Erf Ψ Ψ + π Ψ Ψ Ψ Ψ Ψ Ψ + βθ exp 1 + exp 1 + Erf Θ π Θ + Θ Θ Ψ Ψ Ψ Ψ Ψ Ψ + εσ 1 exp Erfc 0. π σ = σ σ roceedings of the International onference Nuclear Energy for New Europe, (18) (0)
5 506.5 In an experiment the parameters like μ, Θ, τ, σ and β are determined by the method of the plasma production. The potential of the collector Ψ can be determined by the external power supply connected to it. If the collector is floating, then the potential Ψ is determined by the floating condition of zero total current to the collector. This means that J t given by (16) is set to zero. The electron emission from the collector (parameter ε) can be considered as given by the experimental set-up, if the emission is thermal and temperature limited. If it is space charge limited, then the zero electric field condition at the collector (0) is fulfilled. To summarize: if the collector is floating and the emission is critical, equations (16) (with J t = 0), (18), (19) and (0) form a system of 4 equations for 4 unknown quantities: α, Ψ, Ψ and the critical value of ε for a given set of the parameters μ, Θ, τ, σ and β. If the collector potential Ψ is also given, we have to solve the system of 3 equations (18), (19) and (0) for 3 unknown quantities: α, Ψ and the critical value of ε. If also the emission ε is determined (e. g. by the collector temperature) and it is below the space charge limit then we only have to solve equations (18) and (19) for α and Ψ. The potential profile Ψ(z) can be found by solving the oisson equation (17) for a given set of the parameters α, Θ, τ, ε, σ and β. This can only be done numerically. The collector potential Ψ and the derivative dψ/dz (electric field at the collector) must in this case be given as the boundary conditions for the oisson equation (17). 3 RESULTS In this paper we are only interested in a floating collector with a given (temperature limited) electron emission. Our main objective is to compare numerical solutions of the oisson equation (17) with the analytical predictions of the model and with the results of the I simulations. We first select the following parameters: μ = 1/1836 (proton mass), τ = 0.1, β = 0.31, Θ = 0 and σ = 0.01 and then we solve the system of equations (16) (with J t = 0), (18), (19) and (0). In this case the system of equation (16), (18), (19) and (0) has 3 solutions which are the following: 1) Ψ = , Ψ = , α = , ε = ; ) Ψ = , Ψ = , α = , ε =.91603; 3) Ψ = , Ψ = , α = , ε = The solutions are most distinguished by the value of Ψ. We call the solution with the Ψ closest to zero the low solution. The solution with the most negative Ψ is called the high solution and solution with the intermediate Ψ is called the middle solution. The nature of these solutions is discussed elsewhere [1,]. Let us just mention that the low and the high solution correspond to the cases, when the plasma potential is determined by the cool or by the hot electron population, while the middle solution is not a physical one. We now focus on the comparison with numerical solutions. We proceed in the following way: the same parameters as before are selected: μ = 1/1836, τ = 0.1, β = 0.31, Θ = 0 and σ = The emission coefficient ε is now also selected and gradually increased, but always kept below the critical value, found in the above solutions. The system of equations (16), (18) and (19) is then solved for each set of parameters. The solutions are given in Table 1. In this way α, Ψ and Ψ are found. Then α is inserted into (17) together with other parameters and Ψ is used as the boundary condition. In order to solve the equation (17) we need an additional boundary condition, which is the value of the derivative dψ/dz Ψ. This boundary condition is found by guessing. Initial value for Ψ is usually taken close to 1. Then the equation (17) is solved numerically using the 4-th order Runge-Kutta method and the graph of the solution is plotted. From the graph it can be seen easily, whether the selected value for Ψ is to large or to small and it is corrected. When gradually more and more decimal places of Ψ are found the region where Ψ(z) has (almost) a constant value becomes longer and longer. roceedings of the International onference Nuclear Energy for New Europe, 006
6 506.6 Table 1: Solutions of the system of equations (16), (18) and (19) for μ = 1/1836, τ = 0.1, β = 0.31, Θ = 0 and σ = 0.01, while ε is gradually increased. ε(t) Ψ Ψ Ψ (t) Ψ (n) α(t) α(n) β(n) ε(n) Figure 1: Numerical solutions Ψ(z) of the oisson equation (17), the negative derivative (the electric field) -dψ/dz, the particle densities and ΔN for the following parameters: μ = 1/1836, τ = 0.1, β = 0.31, Θ = 0, σ = 0.01 and ε = are plotted versus z. In the left column the low solution is shown, in the middle column the middle solution is shown and in the right column the high solution is shown. roceedings of the International onference Nuclear Energy for New Europe, 006
7 506.7 In Fig. 1 we show a plot of Ψ(z) versus z of the numerical solutions that correspond to the low, middle and high solution of the system of equations (16), (18) and (19) for = 1/1836, τ = 0.1, β = 0.31, Θ = 0, σ = 0.01 and ε =. The corresponding values of Ψ, Ψ, Ψ and α are also given in the Table. In the Table 1 the values of Ψ and α that are obtained from the system of equations (16), (18) and (19) are labelled Ψ (t) and α(t). In separate columns we give the values of Ψ (n) and α(n). These are obtained from the numerical solutions of the equation (17). Note that for the low and the high solution the theoretical value Ψ (t) and the numerical value Ψ (n), which is simply read from the constant portion of the Ψ(z) versus z graph (see Fig. 1) are in very good agreement. For the middle solution the difference between Ψ (t) and Ψ (n) is considerable. In Fig. 1 we also show the electric field profile -dψ/dz, the density profiles of all 4 particle species found from (14) and ΔN = N i N e1 N e N e3. The potential profile Ψ(z) has a local maximum value for the low, the middle and the high solution. But only for the low solution the value of this maximum comes very close to zero and therefore (almost) fulfills the boundary condition (1). So for the low solution this point could be identified as the position of the plasma source. This position therefore depends on the precision of Ψ. Figure : omparison of I simulation (left) and numerical solution of (17) (right). For the low solution the particle densities at the source can be determined and so also the ratios α(n) and β(n) and ε(n) can be determined and compared to the corresponding theoretical values α(t) and β(t) and ε(t). For β the agreement is very good, for ε the difference is exactly a factor of and for α approximately a factor of. In Fig. we show a comparison between a result of a computer I simulation obtained by the XD1 code [3] and a numerical solution of (17). In the simulation τ and Θ at the roceedings of the International onference Nuclear Energy for New Europe, 006
8 506.8 source and σ at the collector can be selected independently. The other independently selected input parameters are influxes of the particle species from the source and from the collector. These influxes must be selected so, that the electric field at the source is zero. Then α, β and ε are obtained as results of the simulation, together with particle densities, collector floating potential and other parameters. In the simulation shown in Fig. we have selected kt 1 = 1 ev, kt = 0 ev, kt i = 0.1 ev, so that τ = 0.1, Θ = 0, ε = 0 and μ = 1/ Then we have found that the collector floating potential is V, which gives directly Ψ = From the simulation we also find β = With these parameters we solve the system of equations (16) and (18) to get α = and Ψ = The later result is in good agreement with simulation. Then we solve the equation (17) numerically for the parameters that τ = 0.1, Θ = 0, ε = 0, μ = 1/ , Ψ = and α = Qualitative and also quantitative agreement of the simulation and the numerical solution is very good. The ratio between the ion and the cool electron density at the source in the computer simulation is α sim =.97, while for the numerical solution it is α(n) =.9. 4 ONLUSIONS We have compared an analytical fully kinetic model of a bounded plasma system with numerical solutions of the corresponding oisson equation and the results of a I simulation. We found very good agreement in the predicted plasma potential. There is also a very good agreement between the numerical solution and the computer simulation in the predicted ion density at the source. AKNOWLEDGMENTS This work been carried out within the Association EURATOM-MHEST. The content of the publication is the sole responsibility of its authors and it does not necessarily represent the view of the ommission or its services. REFERENES [1] T. Gyergyek, M. Čerček, B. Jurčič-Zlobec, zechoslovak journal of physics, vol. 56, (006), suppl. B, pp. B733-B739J. [] T. Gyergyek, M. Čerček, zechoslovak Journal of hysics, vol. 54, (004), pp [3] J.. Verbocoeur, M. V. Alves, V. Vahedi and. K. Birdsall, J. omput. hys., vol. 104, (1993), pp roceedings of the International onference Nuclear Energy for New Europe, 006
Studying the Formation of the Pre-Sheath in an Oblique Magnetic Field using a Fluid Model and PIC Simulation
J. Plasma Fusion Res. SERIES, Vol. 8 (2009) Studying the Formation of the Pre-Sheath in an Oblique Magnetic Field using a Fluid Model and PIC Simulation Jernej KOVAI 1, Tomaž Gyergyek 2,1, Milan EREK 1,3
More information1) Jozef Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia
SI0100095 Nuclear Energy in Central Europe '98 Terme Catez, September 7 to 10, 1998 PLASMA RESPONSE TO A POSITIVE VOLTAGE STEP APPLIED TO AN ANODE IMMERSED IN A WEAKLY MAGNETIZED DISCHARGE PLASMA COLUMN
More informationSheath stability conditions for PIC simulation
Plasma_conference (2011/11/24,Kanazawa,Japan) PIC Sheath stability conditions for PIC simulation Osaka Prefecture. University H. Matsuura, and N. Inagaki, Background Langmuir probe is the most old but
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 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 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 informationMonte Carlo Collisions in Particle in Cell simulations
Monte Carlo Collisions in Particle in Cell simulations Konstantin Matyash, Ralf Schneider HGF-Junior research group COMAS : Study of effects on materials in contact with plasma, either with fusion or low-temperature
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 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 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 informationNeutral beam plasma heating
Seminar I b 1 st year, 2 nd cycle program Neutral beam plasma heating Author: Gabrijela Ikovic Advisor: prof.dr. Tomaž Gyergyek Ljubljana, May 2014 Abstract For plasma to be ignited, external heating is
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 informationToday in Physics 217: boundary conditions and electrostatic boundary-value problems
Today in Physics 17: boundary conditions and electrostatic boundary-value problems Boundary conditions in electrostatics Simple solution of Poisson s equation as a boundary-value problem: the space-charge
More informationFundamental Constants
Fundamental Constants Atomic Mass Unit u 1.660 540 2 10 10 27 kg 931.434 32 28 MeV c 2 Avogadro s number N A 6.022 136 7 36 10 23 (g mol) 1 Bohr magneton μ B 9.274 015 4(31) 10-24 J/T Bohr radius a 0 0.529
More informationFor the following statements, mark ( ) for true statement and (X) for wrong statement and correct it.
Benha University Faculty of Engineering Shoubra Electrical Engineering Department First Year communications. Answer all the following questions Illustrate your answers with sketches when necessary. The
More information6. Stellar spectra. excitation and ionization, Saha s equation stellar spectral classification Balmer jump, H -
6. Stellar spectra excitation and ionization, Saha s equation stellar spectral classification Balmer jump, H - 1 Occupation numbers: LTE case Absorption coefficient: = n i calculation of occupation numbers
More informationYou may not start to read the questions printed on the subsequent pages of this question paper until instructed that you may do so by the Invigilator.
NTURL SCIENCES TRIPOS Part I Saturday 9 June 2007 1.30 pm to 4.30 pm PHYSICS nswer the whole of Section and four questions from Sections B, C and D, with at least one question from each of these Sections.
More informationSemiconductor device structures are traditionally divided into homojunction devices
0. Introduction: Semiconductor device structures are traditionally divided into homojunction devices (devices consisting of only one type of semiconductor material) and heterojunction devices (consisting
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 informationDirect Energy Conversion: Thermionic Conversion
Direct Energy Conversion: Thermionic Conversion References: Direct Energy Conversion by Stanley W. Angrist, Allyn and Beacon, 1982. Direct Energy Conversion by Reiner Decher,Oxford University press, 1997.
More informationGeneral Physics (PHY 2140)
General Physics (PHY 2140) Lecture 4 Electrostatics Electric flux and Gauss s law Electrical energy potential difference and electric potential potential energy of charged conductors http://www.physics.wayne.edu/~apetrov/phy2140/
More informationDiffusion. Diffusion = the spontaneous intermingling of the particles of two or more substances as a result of random thermal motion
Diffusion Diffusion = the spontaneous intermingling of the particles of two or more substances as a result of random thermal motion Fick s First Law Γ ΔN AΔt Γ = flux ΔN = number of particles crossing
More informationSimulation study on the nonlinear EMIC waves
SH21B-2210 Simulation study on the nonlinear EMIC waves Kicheol Rha 1*, Chang-Mo Ryu 1 and Peter H Yoon 2 * lancelot@postech.ac.kr 1 Department of Physics, Pohang University of Science and Technology,
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 informationElectrodynamics of Radiation Processes
Electrodynamics of Radiation Processes 7. Emission from relativistic particles (contd) & Bremsstrahlung http://www.astro.rug.nl/~etolstoy/radproc/ Chapter 4: Rybicki&Lightman Sections 4.8, 4.9 Chapter
More informationThere is light at the end of the tunnel. -- proverb. The light at the end of the tunnel is just the light of an oncoming train. --R.
A vast time bubble has been projected into the future to the precise moment of the end of the universe. This is, of course, impossible. --D. Adams, The Hitchhiker s Guide to the Galaxy There is light at
More informationHigh Beta Discharges with Hydrogen Storage Electrode Biasing in the Tohoku University Heliac
J. Plasma Fusion Res. SERIES, Vol. 8 (2009) High Beta Discharges with Hydrogen Storage Electrode Biasing in the Tohoku University Heliac Hiroyasu UTOH, Kiyohiko NISHIMURA 1), Hajime UMETSU, Keiichi ISHII,
More information6. Stellar spectra. excitation and ionization, Saha s equation stellar spectral classification Balmer jump, H -
6. Stellar spectra excitation and ionization, Saha s equation stellar spectral classification Balmer jump, H - 1 Occupation numbers: LTE case Absorption coefficient: = n i calculation of occupation numbers
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 information3 Minority carrier profiles (the hyperbolic functions) Consider a
Microelectronic Devices and Circuits October 9, 013 - Homework #3 Due Nov 9, 013 1 Te pn junction Consider an abrupt Si pn + junction tat as 10 15 acceptors cm -3 on te p-side and 10 19 donors on te n-side.
More informationParticle Pinch Model of Passing/Trapped High-Z Impurity with Centrifugal Force Effect )
Particle Pinch Model of Passing/Trapped High-Z Impurity with Centrifugal Force Effect ) Yusuke SHIMIZU, Takaaki FUJITA, Atsushi OKAMOTO, Hideki ARIMOTO, Nobuhiko HAYASHI 1), Kazuo HOSHINO 2), Tomohide
More informationCALCULATION OF ISOTOPIC COMPOSITION DURING CONTINUOUS IRRADIATION AND SUBSEQUENT DECAY IN BIOLOGICAL SHIELD OF THE TRIGA MARK ΙΙ REACTOR
International Conference Nuclear Energy for New Europe 2002 Kranjska Gora, Slovenia, September 9-12, 2002 www.drustvo-js.si/gora2002 CALCULATION OF ISOTOPIC COMPOSITION DURING CONTINUOUS IRRADIATION AND
More informationDr. Kasra Etemadi September 21, 2011
Dr. Kasra Etemadi September, 0 - Velocity Distribution -Reaction Rate and Equilibrium (Saha Equation 3-E3 4- Boltzmann Distribution 5- Radiation (Planck s Function 6- E4 z r dxdydz y x Applets f( x r
More information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 7 (2/26/04) Neutron Elastic Scattering - Thermal Motion and Chemical Binding Effects
.54 Neutron Interactions and Applications (Spring 004) Chapter 7 (/6/04) Neutron Elastic Scattering - Thermal Motion and Chemical Binding Effects References -- J. R. Lamarsh, Introduction to Nuclear Reactor
More informationImprovement of a 1-D Fluid Probe Model for Mach Probe Measurements
Improvement of a -D Fluid Probe Model for Mach Probe Measurements P. Peleman, S. Jachmich, M. Van Schoor, G. Van Oost Ghent University Department of Applied Physics, B-9000 Ghent, Belgium Laboratory for
More informationParticle in one-dimensional box
Particle in the box Particle in one-dimensional box V(x) -a 0 a +~ An example of a situation in which only bound states exist in a quantum system. We consider the stationary states of a particle confined
More informationEEE4106Z Radiation Interactions & Detection
EEE4106Z Radiation Interactions & Detection 2. Radiation Detection Dr. Steve Peterson 5.14 RW James Department of Physics University of Cape Town steve.peterson@uct.ac.za May 06, 2015 EEE4106Z :: Radiation
More informationRelativistic Electron Heating in Focused Multimode Laser Fields with Stochastic Phase Purturbations
1 Relativistic Electron Heating in Focused Multimode Laser Fields with Stochastic Phase Purturbations Yu.A.Mikhailov, L.A.Nikitina, G.V.Sklizkov, A.N.Starodub, M.A.Zhurovich P.N.Lebedev Physical Institute,
More informationThermionic emission and Frank-Hertz oscillations
Thermionic emission and Frank-Hertz oscillations Matthew Krupcale, Evan Telford Department of Physics, Case Western Reserve University, Cleveland Ohio, 44106-7079 21 October 2012 Abstract Charges emitted
More informationarxiv: v1 [physics.plasm-ph] 23 May 2013
The characteristic shape of emission profiles of plasma spokes in HiPIMS: the role of secondary electrons A. Hecimovic 1, S. Gallian 2, R.-P. Brinkmann 2, M. Böke 1, J. Winter 1 1 Institute of Experimental
More informationPhysicsAndMathsTutor.com
PhysicsAndMathsTutor.com 1 1. Define electric field strength at a point in space....... [Total 1 marks] 2. Fig. 1 shows a square flat coil of insulated wire placed in a region of a uniform magnetic field
More informationRevision Guide for Chapter 14
Revision Guide for Chapter 14 Contents Revision Checklist Revision Notes Values of the energy kt...4 The Boltzmann factor...4 Thermal activation processes...5 Summary Diagrams Climbing a ladder by chance...7
More informationThomas Schwarz-Selinger Max-Planck-Institut für Plasmaphysik, Garching, Germany
Deuterium retention and isotope exchange studies in self-ion damaged tungsten exposed to neutral atoms Project: Hydrogen retention in self-damaged and Heirradiated tungsten and alloys for PFC Sabina Markelj,
More informationPreview from Notesale.co.uk Page 4 of 35
field 64 If a dielectric is inserted b/w the plates of a charged capacitor, its Remains Becomes infinite capacitance constant decreases increases 65 Selenium is an insulator in the dark but when exposed
More informationDopant Diffusion. (1) Predeposition dopant gas. (2) Drive-in Turn off dopant gas. dose control. Doped Si region
Dopant Diffusion (1) Predeposition dopant gas dose control SiO Si SiO Doped Si region () Drive-in Turn off dopant gas or seal surface with oxide profile control (junction depth; concentration) SiO SiO
More informationThe 5 basic equations of semiconductor device physics: We will in general be faced with finding 5 quantities:
6.012 - Electronic Devices and Circuits Solving the 5 basic equations - 2/12/08 Version The 5 basic equations of semiconductor device physics: We will in general be faced with finding 5 quantities: n(x,t),
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 informationProperties of Error Function erf(z) And Complementary Error Function erfc(z)
Properties of Error Function erf(z) And Complementary Error Function erfc(z) z erf (z) π e -y dy erfc (z) 1 - erf (z) erf () erf( ) 1 erf(- ) - 1 erf (z) d erf(z) dz π z for z
More informationThe Q Machine. 60 cm 198 cm Oven. Plasma. 6 cm 30 cm. 50 cm. Axial. Probe. PUMP End Plate Magnet Coil. Filament Cathode. Radial. Hot Plate.
1 The Q Machine 60 cm 198 cm Oven 50 cm Axial Probe Plasma 6 cm 30 cm PUMP End Plate Magnet Coil Radial Probe Hot Plate Filament Cathode 2 THE Q MACHINE 1. GENERAL CHARACTERISTICS OF A Q MACHINE A Q machine
More informationThermodynamics in Cosmology Nucleosynthesis
Thermodynamics in Cosmology Nucleosynthesis Thermodynamics Expansion Evolution of temperature Freeze out Nucleosynthesis Production of the light elements Potential barrier Primordial synthesis calculations
More informationQUESTION ANSWER. . e. Fourier number:
QUESTION 1. (0 pts) The Lumped Capacitance Method (a) List and describe the implications of the two major assumptions of the lumped capacitance method. (6 pts) (b) Define the Biot number by equations and
More informationTHE FRANCK-HERTZ EXPERIMENT
Rice University Physics 332 THE FRANCK-HERTZ EXPERIMENT I. INTRODUCTION... 2 II. THEORETICAL CONSIDERATIONS... 3 III. MEASUREMENT AND ANALYSIS... 5 Revised July 1991 I. Introduction By the early part of
More informationChapter 10: Steady Heat Conduction
Chapter 0: Steady Heat Conduction In thermodynamics, we considered the amount of heat transfer as a system undergoes a process from one equilibrium state to another hermodynamics gives no indication of
More information11 Quantum theory: introduction and principles
Part 2: Structure Quantum theory: introduction and principles Solutions to exercises E.b E.2b E.3b E.4b E.5b E.6b Discussion questions A successful theory of black-body radiation must be able to explain
More informationFundamental Stellar Parameters. Radiative Transfer. Stellar Atmospheres. Equations of Stellar Structure
Fundamental Stellar Parameters Radiative Transfer Stellar Atmospheres Equations of Stellar Structure Nuclear Reactions in Stellar Interiors Binding Energy Coulomb Barrier Penetration Hydrogen Burning Reactions
More informationEE 5344 Introduction to MEMS CHAPTER 5 Radiation Sensors
EE 5344 Introduction to MEMS CHAPTER 5 Radiation Sensors 5. Radiation Microsensors Radiation µ-sensors convert incident radiant signals into standard electrical out put signals. Radiant Signals Classification
More informationThermal Equilibrium in Nebulae 1. For an ionized nebula under steady conditions, heating and cooling processes that in
Thermal Equilibrium in Nebulae 1 For an ionized nebula under steady conditions, heating and cooling processes that in isolation would change the thermal energy content of the gas are in balance, such that
More informationMATHEMATICS OF DOPING PROFILES. C(x,t) t. = D 2 C(x,t) x 2. 4Dt dx '
EE43 MATHEMATICS OF DOPING PROFILES N. Cheung The diffusion equation with constant D : has the general solution: C(x,t) = C(x,t) = D 2 C(x,t) 4πDt F(x ' ) e -(x-x' ) 2 4Dt dx ' - where F(x') is the C(x,t)
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 informationenergy loss Ionization + excitation of atomic energy levels Mean energy loss rate de /dx proportional to (electric charge) 2 of incident particle
Lecture 4 Particle physics processes - particles are small, light, energetic à processes described by quantum mechanics and relativity à processes are probabilistic, i.e., we cannot know the outcome of
More informationMANIPAL INSTITUTE OF TECHNOLOGY
MANIPAL INSTITUT OF TCHNOLOGY MANIPAL UNIVRSITY, MANIPAL SCOND SMSTR B.Tech. ND-SMSTR XAMINATION - JULY 01 SUBJCT: NGINRING PHYSICS (PHY101/10) Time: Hrs. Max. Marks: 50 Note: Answer any FIV FULL questions.
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 informationGeneral Physics (PHY 2140)
General Physics (PHY 2140) Lecture 2 Electrostatics Electric flux and Gauss s law Electrical energy potential difference and electric potential potential energy of charged conductors http://www.physics.wayne.edu/~alan/
More informationMore Energetics of Alpha Decay The energy released in decay, Q, is determined by the difference in mass of the parent nucleus and the decay products, which include the daughter nucleus and the particle.
More informationTo remain at 0 K heat absorbed by the medium must be removed in the amount of. dq dx = q(l) q(0) dx. Q = [1 2E 3 (τ L )] σ(t T 4 2 ).
294 RADIATIVE HEAT TRANSFER 13.7 Two infinite, isothermal plates at temperatures T 1 and T 2 are separated by a cold, gray medium of optical thickness τ L = L (no scattering). (a) Calculate the radiative
More information12. MHD Approximation.
Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal
More informationChapter 13 Nuclear physics
OCR (A) specifications: 5.4.11i,j,k,l Chapter 13 Nuclear physics Worksheet Worked examples Practical: Simulation (applet) websites nuclear physics End-of-chapter test Marking scheme: Worksheet Marking
More informationStudies of the cathode sheath of a low pressure hollow cathode discharge
Journal of Physics: onference Series Studies of the cathode sheath of a low pressure hollow cathode discharge To cite this article: G Petraconi et al 0 J. Phys.: onf. Ser. 370 004 View the article online
More informationCONTRIBUTION FUNCTION OF MOLYBDENUM ATOM AND IONS IN ASTROPHYSICAL AND LABORATORY PLASMA AS A FUNCTION OF ELECTRON TEMPERATURE
COTRIBUTIO FUCTIO OF MOLYBDEUM ATOM AD IOS I ASTROPHYSICAL AD LABORATORY PLASMA AS A FUCTIO OF ELECTRO TEMPERATURE A.. Jadhav Department of Electronics, Yeshwant Mahavidyalaya, anded. Affiliated to Swami
More informationTheory of optically thin emission line spectroscopy
Theory of optically thin emission line spectroscopy 1 Important definitions In general the spectrum of a source consists of a continuum and several line components. Processes which give raise to the continuous
More informationPlasma Spectroscopy Inferences from Line Emission
Plasma Spectroscopy Inferences from Line Emission Ø From line λ, can determine element, ionization state, and energy levels involved Ø From line shape, can determine bulk and thermal velocity and often
More informationICPMS Doherty Lecture 1
ICPMS Doherty Lecture 1 Mass Spectrometry This material provides some background on how to measure isotope abundances by means of mass spectrometry. Mass spectrometers create and separate ionized atoms
More informationContents Motivation Particle In Cell Method Projects Plasma and Ion Beam Simulations
PIC Method for Numerical Simulation Ninad Joshi NNP Group 1 Contents Motivation Particle In Cell Method Projects Plasma and Ion Beam Simulations Motivation 3 Particle simulation Ion beams and Plasmas Accelerators
More informationControl of gas discharge plasma properties utilizing nonlocal electron energy distribution functions DOE-OFS Plasma Science Center
1 Control of gas discharge plasma properties utilizing nonlocal electron energy distribution functions DOE-OFS Vladimir Demidov Current cooperation with: Mark Koepke (WVU), Igor Kaganovich (PPPL), Yevgeny
More informationThermal Systems. What and How? Physical Mechanisms and Rate Equations Conservation of Energy Requirement Control Volume Surface Energy Balance
Introduction to Heat Transfer What and How? Physical Mechanisms and Rate Equations Conservation of Energy Requirement Control Volume Surface Energy Balance Thermal Resistance Thermal Capacitance Thermal
More informationHigh-Order Representation of Poincaré Maps
High-Order Representation of Poincaré Maps Johannes Grote, Martin Berz, and Kyoko Makino Department of Physics and Astronomy, Michigan State University, East Lansing, MI, USA. {grotejoh,berz,makino}@msu.edu
More information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 7 (2/26/04) Neutron Elastic Scattering - Thermal Motion and Chemical Binding Effects
.54 Neutron Interactions and Applications (Spring 004) Chapter 7 (/6/04) Neutron Elastic Scattering - Thermal Motion and Chemical Binding Effects References -- J. R. Lamarsh, Introduction to Nuclear Reactor
More informationA Method of Knock-on Tail Observation Accounting Temperature Fluctuation Using 6 Li+T/D+T Reaction in Deuterium Plasma
A Method of Knock-on Tail Observation Accounting Temperature Fluctuation Using 6 Li+T/D+T Reaction in Deuterium Plasma Yasuko KAWAMOTO and Hideaki MATSUURA Department of Applied Quantum Physics and Nuclear
More informationThe collision is elastic (KE is conserved)
004 HKAL Physics MC Suggested Solution Let F s : force acting on the object Q by the spring F t : force on the object Q by the thread W: weight of the object Q Vertical balance: F s cos = W.(i) Horizontal
More informationICF ignition and the Lawson criterion
ICF ignition and the Lawson criterion Riccardo Betti Fusion Science Center Laboratory for Laser Energetics, University of Rochester Seminar Massachusetts Institute of Technology, January 0, 010, Cambridge
More informationWhat are Lasers? Light Amplification by Stimulated Emission of Radiation LASER Light emitted at very narrow wavelength bands (monochromatic) Light
What are Lasers? What are Lasers? Light Amplification by Stimulated Emission of Radiation LASER Light emitted at very narrow wavelength bands (monochromatic) Light emitted in a directed beam Light is coherenent
More informationExamination Radiation Physics - 8N120, 2 November
Examination Radiation Physics - 8N0, November 0-4.00-7.00 Four general remarks: This exam consists of 6 assignments on a total of pages. There is a table on page listing the maximum number of that can
More informationP. Diomede, D. J. Economou and V. M. Donnelly Plasma Processing Laboratory, University of Houston
P. Diomede, D. J. Economou and V. M. Donnelly Plasma Processing Laboratory, University of Houston Acknowledgements: DoE Plasma Science Center, NSF Presented at the 57 th AVS Conference, Albuquerque, NM
More informationexp Compared to the values obtained in Example 2.1, we can see that the intrinsic carrier concentration in Ge at T = 300 K is 2.
.1 (a) k =8.617 10 5 ev/k n i (T = 300 K) = 1.66 10 15 (300 K) 3/ 66 ev exp (8.617 10 5 ev/k) (300 K) =.465 10 13 cm 3 n i (T = 600 K) = 1.66 10 15 (600 K) 3/ 66 ev exp (8.617 10 5 ev/k) (600 K) = 4.14
More information11 SEPTEMBER This document consists of printed pages.
S 11 SEPTEMBER 2017 6 Write your name, centre number, index number and class in the spaces at the top of this page and on all work you hand in. Write in dark blue or black pen on both sides of the paper.
More informationBremsstrahlung Radiation
Bremsstrahlung Radiation Wise (IR) An Example in Everyday Life X-Rays used in medicine (radiographics) are generated via Bremsstrahlung process. In a nutshell: Bremsstrahlung radiation is emitted when
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 informationLow Energy Nuclear Fusion Reactions in Solids
Kasagi, J., et al. Low Energy Nuclear Fusion Reactions in Solids. in 8th International Conference on Cold Fusion. 2000. Lerici (La Spezia), Italy: Italian Physical Society, Bologna, Italy. Low Energy Nuclear
More informationCathode thermal effects in one-dimensional space-charge-limited electron flow
Cathode thermal effects in one-dimensional space-charge-limited electron flow Stanley Humphries, Copyright 2012 Field Precision PO Box 13595, Albuquerque, NM 87192 U.S.A. Telephone: +1-505-220-3975 Fax:
More informationThermal Power Calibration of the TRIGA Mark II Reactor
ABSTRACT Thermal Power Calibration of the TRIGA Mark II Reactor Žiga Štancar Jožef Stefan Institute Jamova cesta 39 1000, Ljubljana, Slovenia ziga.stancar@gmail.com Luka Snoj Jožef Stefan Institute Jamova
More informationIntroduction. Introduction. Forces An. Forces An. Forces in Action. Forces in Action. Pressure and Pressure. Pressure and Pressure.
Forces An Introduction A force is a vector quantity. What does this mean? Forces An Introduction A vector quantity, such as force, has a direction as well as a magnitude. 1 1 Forces in Action The moment
More information16.55 Ionized Gases Problem Set #5
16.55 Ionized Gases Problem Set #5 Problem 1: A probe in a non-maxellian plasma The theory of cold probes, as explained in class, applies to plasmas with Maxwellian electron and ion distributions. However,
More informationNuclear Theory - Course 227 NUCLEAR STRUCTURE
Nuclear Theory - Course 227 NUCLEAR STRUCTURE The Nucleus, Nuclear Particles The atomic nucleus consists of Z protons and N neutrons, where Z and N are the atomic number and neutron number respectively.
More informationCyclotron Motion. We can also work-out the frequency of the cyclotron motion. f cyc =
Cyclotron Motion We can also work-out the frequency of the cyclotron motion f cyc = qb 2πm Neil Alberding (SFU Physics) Physics 121: Optics, Electricity & Magnetism Spring 2010 1 / 1 Cyclotron Motion We
More informationExplain how Planck resolved the ultraviolet catastrophe in blackbody radiation. Calculate energy of quanta using Planck s equation.
Objectives Explain how Planck resolved the ultraviolet catastrophe in blackbody radiation. Calculate energy of quanta using Planck s equation. Solve problems involving maximum kinetic energy, work function,
More informationComplex functions in the theory of 2D flow
Complex functions in the theory of D flow Martin Scholtz Institute of Theoretical Physics Charles University in Prague scholtz@utf.mff.cuni.cz Faculty of Transportation Sciences Czech Technical University
More informationMatti Laan Gas Discharge Laboratory University of Tartu ESTONIA
Matti Laan Gas Discharge Laboratory University of Tartu ESTONIA Outline 1. Ionisation 2. Plasma definition 3. Plasma properties 4. Plasma classification 5. Energy transfer in non-equilibrium plasma 6.
More informationExercise 1 Atomic line spectra 1/9
Exercise 1 Atomic line spectra 1/9 The energy-level scheme for the hypothetical one-electron element Juliettium is shown in the figure on the left. The potential energy is taken to be zero for an electron
More informationRutherford Backscattering Spectrometry
Rutherford Backscattering Spectrometry EMSE-515 Fall 2005 F. Ernst 1 Bohr s Model of an Atom existence of central core established by single collision, large-angle scattering of alpha particles ( 4 He
More informationANALYSIS OF PLASMA FACING MATERIALS IN CONTROLLED FUSION DEVICES. Marek Rubel
ANALYSIS OF PLASMA FACING MATERIALS IN CONTROLLED FUSION DEVICES Marek Rubel Alfvén Laboratory, Royal Institute of Technology, Association EURATOM VR, Stockholm, Sweden Acknowledgements Paul Coad and Guy
More information