Solution of time-dependent Boltzmann equation for electrons in non-thermal plasma
|
|
- Oswald Boone
- 5 years ago
- Views:
Transcription
1 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, Czech Republic Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
2 Description of plasma Kinetic description: most fundamental way providing f(r, v, t) achieved by solving the Boltzmann equation self-consistently determined the electromagnetic fields via Maxwell s equations Fluid description: description of the plasma based on macroscopic quantities fluid equations obtained by taking velocity moments of the BE Hybrid description: a combination of fluid and kinetic models treating some components of the system as a fluid and others kinetically Gyrokinetic Description: appropriate to systems with a strong background magnetic field the kinetic equations are averaged over the fast circular motion of the gyroradius Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
3 What will I talk about? weakly ionized plasma (small fraction of the atoms are ionized) moderate pressure ( 1 Torr) non-thermal or cold plasma (T e >> T ion = T gas ) non-magnetized plasma (B=0) solution of linear Boltzmann equation for electrons Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
4 What will I talk about? weakly ionized plasma (small fraction of the atoms are ionized) moderate pressure ( 1 Torr) non-thermal or cold plasma (T e >> T ion = T gas ) non-magnetized plasma (B=0) solution of linear Boltzmann equation for electrons Are you still interested? Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
5 The Boltzmann equation The Boltzmann equation is an equation for the time evolution of the distribution function f(r, v, t) in phase space, where r and v are position and velocity f t + v rf + a v f = f t We will discussed the Boltzmann equation with a collision term for two-body collisions between particles of different types that are assumed to be uncorrelated prior to the collision: the linear Boltzmann equation. c Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
6 The Boltzmann equation The Boltzmann equation is an equation for the time evolution of the distribution function f(r, v, t) in phase space, where r and v are position and velocity f t + v rf + a v f = f t We will discussed the Boltzmann equation with a collision term for two-body collisions between particles of different types that are assumed to be uncorrelated prior to the collision: the linear Boltzmann equation. How we can solve it? c Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
7 How to solve the Boltzmann equation? There are generally two ways: Monte Carlo methods (TPMC, PIC/MCC) A way how to solve it without solving it Actual solution of The Equation difficult but in simple cases and with approximations it is feasible Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
8 Monte Carlo methods a way how to obtain results with minimum of approximations MC: Solves a particle charge transport by the numerical simulation of large particle ensemble, collision processes are introduced by generating appropriately distributed random numbers. Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
9 Monte Carlo methods a way how to obtain results with minimum of approximations MC: Solves a particle charge transport by the numerical simulation of large particle ensemble, collision processes are introduced by generating appropriately distributed random numbers. How does it work? Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
10 Consider n-dimensional integral I = f(x 1,..., x n ) d n x, in n-cube of unit volume. It is the mean value I = f. MC gives us an approximation of exact value of the mean value with statistic estimate based on a numerical experiment. The mean value is estimated on an average over N trials I N i=1 f(η 1,...,η n ), which is exact for N N Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
11 Consider n-dimensional integral I = f(x 1,..., x n ) d n x, in n-cube of unit volume. It is the mean value I = f. MC gives us an approximation of exact value of the mean value with statistic estimate based on a numerical experiment. The mean value is estimated on an average over N trials I N i=1 f(η 1,...,η n ), which is exact for N N But why random numbers? Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
12 Consider n-dimensional integral I = f(x 1,..., x n ) d n x, in n-cube of unit volume. It is the mean value I = f. MC gives us an approximation of exact value of the mean value with statistic estimate based on a numerical experiment. The mean value is estimated on an average over N trials I N i=1 f(η 1,...,η n ), which is exact for N N Because a finite distribution of regular points gives a poor resolution whent the dimension n is high. Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
13 Propagation method Problem to be solved: determination of distribution function f of charged particles based on the initial condition f 0. Initial condition is a sum of mathematical particles (δ-functions) which are propagated in time and space to represent a solution at later times Consider N simulated superparticles, these represent N real real particles w = N real N weight of superparticles Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
14 Propagation method Problem to be solved: determination of distribution function f of charged particles based on the initial condition f 0. Initial condition is a sum of mathematical particles (δ-functions) which are propagated in time and space to represent a solution at later times Consider N simulated superparticles, these represent N real real particles w = N real N weight of superparticles How to determine the density and the energy distribution function? Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
15 Estimation of particle density in a given cell of the mathematical mesh n = N cellw x y z Determination of energy distribution function by introducing a mash on the kinetic energy axis Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
16 How to propagate particles? performing a numerical solution of equations of motion under the effect of Lorentz force and random force which represents collision events between charged particles and neutral particles background dv t dr dt = v = q m E + F coll m collision events are instantaneous, successions of kicks Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
17 How to propagate particles? performing a numerical solution of equations of motion under the effect of Lorentz force and random force which represents collision events between charged particles and neutral particles background dv t dr dt = v = q m E + F coll m collision events are instantaneous, successions of kicks How are collision times, t c, distributed? Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
18 How to propagate particles? performing a numerical solution of equations of motion under the effect of Lorentz force and random force which represents collision events between charged particles and neutral particles background dv t dr dt = v = q m E + F coll m collision events are instantaneous, successions of kicks How are collision times, t c, distributed? The answer is simple for constant collision frequency! Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
19 Distribution of collision times for constant collision frequency α P(t c ) = α exp( αt c ) free flight time can be generated easily t c = lnη α. Implementation of propagation: particle collisionless equation of motion between (t, t + t c ) is solved dr = v dt dv = q t m E Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
20 Distribution of collision times for constant collision frequency α P(t c ) = α exp( αt c ) free flight time can be generated easily t c = lnη α. Implementation of propagation: particle collisionless equation of motion between (t, t + t c ) is solved dr = v dt dv = q t m E And... Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
21 Distribution of collision times for constant collision frequency α P(t c ) = α exp( αt c ) free flight time can be generated easily t c = lnη α. Implementation of propagation: particle collisionless equation of motion between (t, t + t c ) is solved dr = v dt dv = q t m E the collision is introduced by laws of physics, based on appropriate random distribution of impact parameters Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
22 But... Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
23 Collision frequency is not a constant! Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
24 Null collision technique A new collision frequency is introduced α max = max α(v, r). r,v w α(v, r) = p n(r)vσ p (v) α(v, r) is a sum of all contributions from different collision processes p A new null collision process is added in order to equate α(v, r) and α max Non-physical collisions are eliminated by considering null collision with probability 1 α(v, r)/α max Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
25 Null collision technique selection of an event with a probability p: Uniform η (0, 1] : (p > η? accepted : rejected) after the collision a new direction of velocity and its magnitude is determined Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
26 Null collision technique Isotropic scattering can be assumed in many cases, momentum transfer cross section is used insted of DCS. Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
27 Null collision technique Isotropic scattering can be assumed in many cases, momentum transfer cross section is used insted of DCS. Is the Monte Carlo technique merely a solution on the bases of physical analogy? Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
28 Formal view of Monte Carlo technique Consider linear Boltzmann equation in the form df dt = α maxf + α max Af A is a linear combination of unit operator and a Boltzmann collision operator The integral form: t f(t) = exp( α max t) + α max exp( α max τ)a(τ) f(τ) dτ Monte carlo method can be derived from numerical solution of this equation. MC proceure is applied to calculate this integral. 0 Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
29 Solving the BE by solving it Approximations: spatialy homogeneous electrons collide with neutral gas particles collisions are isotropic ionization is taken to be conservative gas particles at rest The Boltzmann equation solved by the expansion into the Legendre polynomials Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
30 Expansion of the BE into the LPs E z-axis BE has a form = f(v, t) f(v, θ, t) (rotational symmetry) f t + ee m f = f v z t. c Distribution function expanded into Legendre polynomials: f(v, θ, t) = f n (v, t) P n (cosθ). n=0 n=0 [ fn t + ee m + ( n f n 1 + n + 1 f n+1 2n 1 v 2n + 3 v (n + 2)(n + 1) f n+1 2n + 3 v n(n 1) 2n 1 f n 1 v + )] P n (cos θ) = f t. c Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
31 Expansion of BE into LPs Resulting hierarchy of PDEs f n(v, t) = t n» ee(t) 2n 1 m n + 1» ee(t) 2n + 3 m v n 1 v 1 v n+2 v fn 1 (v, t) v n 1 «`v n+2 f n+1 (v, t) + + vn(t)`q n(v) Q 0 (v) f n(v, t) + + m 1 h i v 4 N(t)`Q M v 2 1n (v) Q n(v) f n(v, t) + v + X» (v k A ) 2 N(t)Q v n(v k A)f k n(va, k t) k vn(t)q0(v)f k n(v, t), n 0, f 1 0, where... Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
32 ... where Q n(v) = Q 1n (v) = Q k n(v) = Z 2π 0 Z 2π 0 Z 2π 0 dφ dφ dφ Z π 0 Z π 0 Z π 0 σ el (v, θ)p n(cos θ) sin θ dθ, σ el (v, θ)p 1 (cos θ)p n(cos θ) sin θ dθ, σ k (v, θ)p n(cos θ) sin θ dθ, where σ el (v) σ k (v) are differential collision cross sections for elastic and k-th inelastic collision process and «va k = v 2 + 2Ak, A k > 0 m A k is excitation energy of k-th state. For numerical solution kinetic energy U used instead of velocity v: U = (v/γ) 2, γ = (2/m) 1/2. Isotropic scattering (introducing mt) Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
33 Different approximations of the expansion of BE Two-term approximations: conventional two-term approx. f 1 / t 0, f 1 (U, t) E(t) f 0 (U, t)/ U you cannot choose f 1 (U, 0)! strict time-dependent two-term approx. f 2 (U, t) 0 Multi-term approximation: Expansion truncated after arbitrary (> 2) number of terms Effective field approximation (2-term): E(t) = E 0 exp jωt Time variation f 0 a f 1 slow enough compared to ω. Generally: ω/2π GHz or higher Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
34 Solution procedure Hierarchy of PDEs solved as a initial boundary value problem f 0 U = 0, f n (U = 0, t) = 0, n 1 U=0 f n (U U, t) = 0, n 0, U sufficiently large energy, by the method of lines BE discretized in energy space on uniform mesh U (0, U ) system of ODEs received ODEs integrated by Runge Kutta method of order 4(5) with step-size control Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
35 Case study: Reid ramp model Reid ramp model: Classical benchmark problem for solvers of BE, (Reid I. D., 1979, Aust. J. Phys. 32 p 231) Hypothetical gas mass: 4 amu cross section for elastic collisions σ el = m 2 cross section for inelastic collisions { 0 ε 0.2 ev σ in = 10(ε 0.2) m 2 ε > 0.2 ev All collisions assumed to be isotropic Thermal motion of gas particles neglected Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
36 Reid ramp model: results Time independent electric field Initial distribution of electrons: Maxwellian 300 K. Steady state for E/N = 24 Td (1 Td = V m 1 ) BE MC Reid Ness Drift velocity (10 4 m s 1 ) Mean energy (ev) Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
37 v d [10 4 m/s] Reid ramp model: E/N = 24 Td, 1 Torr BE MC ε [ev] time [s] Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
38 Rare gases and the effect of negative mobility f n f 0 final f 0 initial Argon Q el Q(ε) [10-20 m 2 ] v d (10 4 m/s) v d ε ε (ev) Argon: energy [ev] time (s) initial distribution: Gaussian, center 5.5 ev, width 1 ev gas density m 3 (pressure 1 Torr, temperature 273 K) E = 2 kv m 1, E/N = 58 Td. Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
39 Rare gases and the effect of negative mobility ε 1/2 f 0 /n e [ev -1 ], ε f 1 /n e [ev -1/2 ] Xenon: s f 0 f 1 > 0 f 1 < ε [ev] coll. cross section [a.u.] v d (10 4 m/s) ε 2.0 v d time (s) ε (ev) initial distribution: Gaussian, center 15.0 ev, width 1 ev gas density m 3 (pressure 1 Torr, temperature 273 K) E = 2 kv m 1, E/N = 58 Td Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
40 Conclusions For a solution of time-dependent Boltzmann equation the Monte Carlo Method: is easy to implement, rather straightforward method, even for complex boundary conditions, may be rather time consuming Solution of the Equation: is a challenge, it is a fun, it gives good results (multi-term), for many problems it is impossible Z. Bonaventura (UFE, F.Sci. MUNI) Simulation of Boltzmann equation Ondřejov, / 41
Solution of Time-dependent Boltzmann Equation
WDS'5 Proceedings of Contributed Papers, Part III, 613 619, 5. ISBN 8-8673-59- MATFYZPRESS Solution of Time-dependent Boltzmann Equation Z. Bonaventura, D. Trunec, and D. Nečas Masaryk University, Faculty
More informationThe Role of Secondary Electrons in Low Pressure RF Glow Discharge
WDS'05 Proceedings of Contributed Papers, Part II, 306 312, 2005. ISBN 80-86732-59-2 MATFYZPRESS The Role of Secondary Electrons in Low Pressure RF Glow Discharge O. Brzobohatý and D. Trunec Department
More informationA theoretical study of the energy distribution function of the negative hydrogen ion H - in typical
Non equilibrium velocity distributions of H - ions in H 2 plasmas and photodetachment measurements P.Diomede 1,*, S.Longo 1,2 and M.Capitelli 1,2 1 Dipartimento di Chimica dell'università di Bari, Via
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 informationScattering in Cold- Cathode Discharges
Simulating Electron Scattering in Cold- Cathode Discharges Alexander Khrabrov, Igor Kaganovich*, Vladimir I. Demidov**, George Petrov*** *Princeton Plasma Physics Laboratory ** Wright-Patterson Air Force
More informationModelling of low-temperature plasmas: kinetic and transport mechanisms. L.L. Alves
Modelling of low-temperature plasmas: kinetic and transport mechanisms L.L. Alves llalves@tecnico.ulisboa.pt Instituto de Plasmas e Fusão Nuclear Instituto Superior Técnico, Universidade de Lisboa Lisboa,
More informationPHYSICS OF HOT DENSE PLASMAS
Chapter 6 PHYSICS OF HOT DENSE PLASMAS 10 26 10 24 Solar Center Electron density (e/cm 3 ) 10 22 10 20 10 18 10 16 10 14 10 12 High pressure arcs Chromosphere Discharge plasmas Solar interior Nd (nω) laserproduced
More informationPlasmas as fluids. S.M.Lea. January 2007
Plasmas as fluids S.M.Lea January 2007 So far we have considered a plasma as a set of non intereacting particles, each following its own path in the electric and magnetic fields. Now we want to consider
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 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 informationMulti-fluid Simulation Models for Inductively Coupled Plasma Sources
Multi-fluid Simulation Models for Inductively Coupled Plasma Sources Madhusudhan Kundrapu, Seth A. Veitzer, Peter H. Stoltz, Kristian R.C. Beckwith Tech-X Corporation, Boulder, CO, USA and Jonathan Smith
More informationParallel Ensemble Monte Carlo for Device Simulation
Workshop on High Performance Computing Activities in Singapore Dr Zhou Xing School of Electrical and Electronic Engineering Nanyang Technological University September 29, 1995 Outline Electronic transport
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 informationANGULAR DEPENDENCE OF ELECTRON VELOCITY DISTRIBUTIONS IN LOW-PRESSURE INDUCTIVELY COUPLED PLASMAS 1
ANGULAR DEPENDENCE OF ELECTRON VELOCITY DISTRIBUTIONS IN LOW-PRESSURE INDUCTIVELY COUPLED PLASMAS 1 Alex V. Vasenkov 2, and Mark J. Kushner Department of Electrical and Computer Engineering Urbana, IL
More information26. Non-linear effects in plasma
Phys780: Plasma Physics Lecture 26. Non-linear effects. Collisionless shocks.. 1 26. Non-linear effects in plasma Collisionless shocks ([1], p.405-421, [6], p.237-245, 249-254; [4], p.429-440) Collisionless
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationKinetic theory of gases
Kinetic theory of gases Toan T. Nguyen Penn State University http://toannguyen.org http://blog.toannguyen.org Graduate Student seminar, PSU Jan 19th, 2017 Fall 2017, I teach a graduate topics course: same
More informationMacroscopic plasma description
Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationTheory of Gas Discharge
Boris M. Smirnov Theory of Gas Discharge Plasma l Springer Contents 1 Introduction 1 Part I Processes in Gas Discharge Plasma 2 Properties of Gas Discharge Plasma 13 2.1 Equilibria and Distributions of
More informationKinetic Solvers with Adaptive Mesh in Phase Space for Low- Temperature Plasmas
Kinetic Solvers with Adaptive Mesh in Phase Space for Low- Temperature Plasmas Vladimir Kolobov, a,b,1 Robert Arslanbekov a and Dmitry Levko a a CFD Research Corporation, Huntsville, AL 35806, USA b The
More informationCharacteristics and classification of plasmas
Characteristics and classification of plasmas PlasTEP trainings course and Summer school 2011 Warsaw/Szczecin Indrek Jõgi, University of Tartu Partfinanced by the European Union (European Regional Development
More informationKINETIC DESCRIPTION OF MAGNETIZED TECHNOLOGICAL PLASMAS
KINETIC DESCRIPTION OF MAGNETIZED TECHNOLOGICAL PLASMAS Ralf Peter Brinkmann, Dennis Krüger Fakultät für Elektrotechnik und Informationstechnik Lehrstuhl für Theoretische Elektrotechnik Magnetized low
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 informationSlowing down the neutrons
Slowing down the neutrons Clearly, an obvious way to make a reactor work, and to make use of this characteristic of the 3 U(n,f) cross-section, is to slow down the fast, fission neutrons. This can be accomplished,
More information14. Energy transport.
Phys780: Plasma Physics Lecture 14. Energy transport. 1 14. Energy transport. Chapman-Enskog theory. ([8], p.51-75) We derive macroscopic properties of plasma by calculating moments of the kinetic equation
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 informationSupplementary Information
1 Supplementary Information 3 Supplementary Figures 4 5 6 7 8 9 10 11 Supplementary Figure 1. Absorbing material placed between two dielectric media The incident electromagnetic wave propagates in stratified
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 informationAdvances in the kinetic simulation of microwave absorption in an ECR thruster
Advances in the kinetic simulation of microwave absorption in an ECR thruster IEPC-2017-361 Presented at the 35th International Electric Propulsion Conference Georgia Institute of Technology Atlanta, Georgia
More informationCHAPTER 9. Microscopic Approach: from Boltzmann to Navier-Stokes. In the previous chapter we derived the closed Boltzmann equation:
CHAPTER 9 Microscopic Approach: from Boltzmann to Navier-Stokes In the previous chapter we derived the closed Boltzmann equation: df dt = f +{f,h} = I[f] where I[f] is the collision integral, and we have
More informationOn the use of classical transport analysis to determine cross-sections for low-energy e H 2 vibrational excitation
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS B: ATOMIC, MOLECULAR AND OPTICAL PHYSICS J. Phys. B: At. Mol. Opt. Phys. 35 (2002) 605 626 PII: S0953-4075(02)24907-X On the use of classical transport
More informationMonte Carlo radiation transport codes
Monte Carlo radiation transport codes How do they work? Michel Maire (Lapp/Annecy) 23/05/2007 introduction to Monte Carlo radiation transport codes 1 Decay in flight (1) An unstable particle have a time
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 informationAccurate representation of velocity space using truncated Hermite expansions.
Accurate representation of velocity space using truncated Hermite expansions. Joseph Parker Oxford Centre for Collaborative Applied Mathematics Mathematical Institute, University of Oxford Wolfgang Pauli
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 Dec The Maxwell Velocity Distribution
Physics 301 7-Dec-2005 29-1 The Maxwell Velocity Distribution The beginning of chapter 14 covers some things we ve already discussed. Way back in lecture 6, we calculated the pressure for an ideal gas
More informationPlasma collisions and conductivity
e ion conductivity Plasma collisions and conductivity Collisions in weakly and fully ionized plasmas Electric conductivity in non-magnetized and magnetized plasmas Collision frequencies In weakly ionized
More informationNARROW GAP ELECTRONEGATIVE CAPACITIVE DISCHARGES AND STOCHASTIC HEATING
NARRW GAP ELECTRNEGATIVE CAPACITIVE DISCHARGES AND STCHASTIC HEATING M.A. Lieberman, E. Kawamura, and A.J. Lichtenberg Department of Electrical Engineering and Computer Sciences University of California
More informationIntroduction Statistical Thermodynamics. Monday, January 6, 14
Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can
More informationElectron Transport Coefficients in a Helium Xenon Mixture
ISSN 1068-3356, Bulletin of the Lebedev Physics Institute, 2014, Vol. 41, No. 10, pp. 285 291. c Allerton Press, Inc., 2014. Original Russian Text c S.A. Mayorov, 2014, published in Kratkie Soobshcheniya
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 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 informationPhysics Fall Mechanics, Thermodynamics, Waves, Fluids. Lecture 20: Rotational Motion. Slide 20-1
Physics 1501 Fall 2008 Mechanics, Thermodynamics, Waves, Fluids Lecture 20: Rotational Motion Slide 20-1 Recap: center of mass, linear momentum A composite system behaves as though its mass is concentrated
More informationDESIGN CONSIDERATIONS FOR A LABORATORY DUSTY PLASMA WITH MAGNETIZED DUST PARTICLES
1 DESIGN CONSIDERATIONS FOR A LABORATORY DUSTY PLASMA WITH MAGNETIZED DUST PARTICLES Bob Merlino Department of Physics and Astronomy The University of Iowa Iowa City, IA 52242 October 15, 2009 The design
More informationAnalysis of recombination and relaxation of non-equilibrium air plasma generated by short time energetic electron and photon beams
22 nd International Symposium on Plasma Chemistry July 5-10, 2015; Antwerp, Belgium Analysis of recombination and relaxation of non-equilibrium air plasma generated by short time energetic electron and
More informationMonte Carlo methods for kinetic equations
Monte Carlo methods for kinetic equations Lecture 4: Hybrid methods and variance reduction Lorenzo Pareschi Department of Mathematics & CMCS University of Ferrara Italy http://utenti.unife.it/lorenzo.pareschi/
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 informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More informationExtremely far from equilibrium: the multiscale dynamics of streamer discharges
Extremely far from equilibrium: the multiscale dynamics of streamer discharges Ute Ebert 1,2 1 Centrum Wiskunde & Informatica Amsterdam 2 Eindhoven University of Technology http://www.cwi.nl/~ebert www.cwi.nl/~ebert
More informationChapter V: Interactions of neutrons with matter
Chapter V: Interactions of neutrons with matter 1 Content of the chapter Introduction Interaction processes Interaction cross sections Moderation and neutrons path For more details see «Physique des Réacteurs
More informationRandom Averaging. Eli Ben-Naim Los Alamos National Laboratory. Paul Krapivsky (Boston University) John Machta (University of Massachusetts)
Random Averaging Eli Ben-Naim Los Alamos National Laboratory Paul Krapivsky (Boston University) John Machta (University of Massachusetts) Talk, papers available from: http://cnls.lanl.gov/~ebn Plan I.
More informationEXAMINATION QUESTION PAPER
Faculty of Science and Technology EXAMINATION QUESTION PAPER Exam in: Fys-2009 Introduction to Plasma Physics Date: 20161202 Time: 09.00-13.00 Place: Åsgårdvegen 9 Approved aids: Karl Rottmann: Matematisk
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 informationPLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT. Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford
1 PLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford 2 Overview Why do we need plasmas? For fusion, among other things
More informationDispersive Media, Lecture 7 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2017-02-14 Dispersive Media, Lecture 7 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasmas as a coupled system Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas
More informationModeling and Simulation of Plasma Based Applications in the Microwave and RF Frequency Range
Modeling and Simulation of Plasma Based Applications in the Microwave and RF Frequency Range Dr.-Ing. Frank H. Scharf CST of America What is a plasma? What is a plasma? Often referred to as The fourth
More informationUltra-Cold Plasma: Ion Motion
Ultra-Cold Plasma: Ion Motion F. Robicheaux Physics Department, Auburn University Collaborator: James D. Hanson This work supported by the DOE. Discussion w/ experimentalists: Rolston, Roberts, Killian,
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 informationThe Physics of Fluids and Plasmas
The Physics of Fluids and Plasmas An Introduction for Astrophysicists ARNAB RAI CHOUDHURI CAMBRIDGE UNIVERSITY PRESS Preface Acknowledgements xiii xvii Introduction 1 1. 3 1.1 Fluids and plasmas in the
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 informationDatabase for Imaging and Transport with Energetic Neutral Atoms (ENAs)
Database for Imaging and Transport with Energetic Neutral Atoms (ENAs) Nicholas Lewkow Harvard-Smithsonian Center for Astrophysics University of Connecticut April 26, 213 Collaborators Vasili Kharchenko
More informationTwo-stage Rydberg charge exchange in a strong magnetic field
Two-stage Rydberg charge exchange in a strong magnetic field M. L. Wall, C. S. Norton, and F. Robicheaux Department of Physics, Auburn University, Auburn, Alabama 36849-5311, USA Received 21 June 2005;
More informationTitle. Author(s)Sugawara, Hirotake. CitationPlasma Sources Science and Technology, 26(4): Issue Date Doc URL. Rights.
Title A relaxation-accelerated propagator method for calcu transport parameters in gas under dc electric fields Author(s)Sugawara, Hirotake CitationPlasma Sources Science and Technology, 26(4): 044002
More informationModeling of noise by Legendre polynomial expansion of the Boltzmann equation
Modeling of noise by Legendre polynomial expansion of the Boltzmann equation C. Jungemann Institute of Electromagnetic Theory RWTH Aachen University 1 Outline Introduction Theory Noise Legendre polynomial
More informationMONTE CARLO SIMULATION OF RADIATION TRAPPING IN ELECTRODELESS LAMPS: A STUDY OF COLLISIONAL BROADENERS*
MONTE CARLO SIMULATION OF RADIATION TRAPPING IN ELECTRODELESS LAMPS: A STUDY OF COLLISIONAL BROADENERS* Kapil Rajaraman** and Mark J. Kushner*** **Department of Physics ***Department of Electrical and
More informationSCALING OF PLASMA SOURCES FOR O 2 ( 1 ) GENERATION FOR CHEMICAL OXYGEN-IODINE LASERS
SCALING OF PLASMA SOURCES FOR O 2 ( 1 ) GENERATION FOR CHEMICAL OXYGEN-IODINE LASERS D. Shane Stafford and Mark J. Kushner Department of Electrical and Computer Engineering Urbana, IL 61801 http://uigelz.ece.uiuc.edu
More information2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasma physics Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas Transverse waves
More informationEffect of Applied Electric Field and Pressure on the Electron Avalanche Growth
Effect of Applied Electric Field and Pressure on the Electron Avalanche Growth L. ZEGHICHI (), L. MOKHNACHE (2), and M. DJEBABRA (3) () Department of Physics, Ouargla University, P.O Box.5, OUARGLA 3,
More informationCharged particle motion in external fields
Chapter 2 Charged particle motion in external fields A (fully ionized) plasma contains a very large number of particles. In general, their motion can only be studied statistically, taking appropriate averages.
More informationPIC-MCC/Fluid Hybrid Model for Low Pressure Capacitively Coupled O 2 Plasma
PIC-MCC/Fluid Hybrid Model for Low Pressure Capacitively Coupled O 2 Plasma Kallol Bera a, Shahid Rauf a and Ken Collins a a Applied Materials, Inc. 974 E. Arques Ave., M/S 81517, Sunnyvale, CA 9485, USA
More informationCHAPTER 4. Basics of Fluid Dynamics
CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
More informationVOLUMETRIC CONTROL OF ANISOTROPIC ELECTRON DISTRIBUTION FUNCTION IN PLASMAS WITH LANGMUIR OSCILLATIONS
VOLUMETRIC CONTROL OF ANISOTROPIC ELECTRON DISTRIBUTION FUNCTION IN PLASMAS WITH LANGMUIR OSCILLATIONS A.S. Mustafaev St. Petersburg State Mining University, Russia V.I. Demidov West Virginia University,
More informationANTENNAS. Vector and Scalar Potentials. Maxwell's Equations. E = jωb. H = J + jωd. D = ρ (M3) B = 0 (M4) D = εe
ANTENNAS Vector and Scalar Potentials Maxwell's Equations E = jωb H = J + jωd D = ρ B = (M) (M) (M3) (M4) D = εe B= µh For a linear, homogeneous, isotropic medium µ and ε are contant. Since B =, there
More informationBeams and magnetized plasmas
Beams and magnetized plasmas 1 Jean-Pierre BOEUF LAboratoire PLAsma et Conversion d Energie LAPLACE/ CNRS, Université Paul SABATIER, TOULOUSE Beams and magnetized plasmas 2 Outline Ion acceleration and
More informationA comparative study of no-time-counter and majorant collision frequency numerical schemes in DSMC
Purdue University Purdue e-pubs School of Aeronautics and Astronautics Faculty Publications School of Aeronautics and Astronautics 2012 A comparative study of no-time-counter and majorant collision frequency
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 informationLong range interaction --- between ions/electrons and ions/electrons; Coulomb 1/r Intermediate range interaction --- between
Collisional Processes Long range interaction --- between ions/electrons and ions/electrons; Coulomb 1/r Intermediate range interaction --- between ions/electrons and neutral atoms/molecules; Induced dipole
More informationGeo-nanomaterials ION DISTRIBUTION FUNCTION IN THEIR OWN GAS PLASMA
DOI.8454/PMI.6.6.864 UDC 533.9.8 ION DISTRIBUTION FUNCTION IN THEIR OWN GAS PLASMA А.S. MUSTAFAEV, V.S. SOUKHOMLINOV Saint-Petersburg Mining Universit Saint-Petersburg, Russia Saint-Petersburg State Universit
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 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 information8.333: Statistical Mechanics I Fall 2007 Test 2 Review Problems
8.333: Statistical Mechanics I Fall 007 Test Review Problems The second in-class test will take place on Wednesday 10/4/07 from :30 to 4:00 pm. There will be a recitation with test review on Monday 10//07.
More informationSimilarities and differences:
How does the system reach equilibrium? I./9 Chemical equilibrium I./ Equilibrium electrochemistry III./ Molecules in motion physical processes, non-reactive systems III./5-7 Reaction rate, mechanism, molecular
More informationElectron Temperature Modification in Gas Discharge Plasma
Electron Temperature Modification in Gas Discharge Plasma Valery Godyak University of Michigan and RF Plasma Consulting egodyak@comcast.net Workshop: Control of Distribution Functions in Low Temperature
More informationKinetic, Fluid & MHD Theories
Lecture 2 Kinetic, Fluid & MHD Theories The Vlasov equations are introduced as a starting point for both kinetic theory and fluid theory in a plasma. The equations of fluid theory are derived by taking
More informationThe Computational Simulation of the Positive Ion Propagation to Uneven Substrates
WDS' Proceedings of Contributed Papers, Part II, 5 9,. ISBN 978-8-778-85-9 MATFYZPRESS The Computational Simulation of the Positive Ion Propagation to Uneven Substrates V. Hrubý and R. Hrach Charles University,
More informationSingle Particle Motion
Single Particle Motion C ontents Uniform E and B E = - guiding centers Definition of guiding center E gravitation Non Uniform B 'grad B' drift, B B Curvature drift Grad -B drift, B B invariance of µ. Magnetic
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 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 informationAMSC 663 Project Proposal: Upgrade to the GSP Gyrokinetic Code
AMSC 663 Project Proposal: Upgrade to the GSP Gyrokinetic Code George Wilkie (gwilkie@umd.edu) Supervisor: William Dorland (bdorland@umd.edu) October 11, 2011 Abstract Simulations of turbulent plasma in
More informationRigid Body Motion in a Special Lorentz Gas
Rigid Body Motion in a Special Lorentz Gas Kai Koike 1) Graduate School of Science and Technology, Keio University 2) RIKEN Center for Advanced Intelligence Project BU-Keio Workshop 2018 @Boston University,
More informationSemiclassical Electron Transport
Semiclassical Electron Transport Branislav K. Niolić Department of Physics and Astronomy, University of Delaware, U.S.A. PHYS 64: Introduction to Solid State Physics http://www.physics.udel.edu/~bniolic/teaching/phys64/phys64.html
More informationarxiv: v1 [physics.plasm-ph] 10 Nov 2014
arxiv:1411.2464v1 [physics.plasm-ph] 10 Nov 2014 Effects of fast atoms and energy-dependent secondary electron emission yields in PIC/MCC simulations of capacitively coupled plasmas A. Derzsi 1, I. Korolov
More informationFokker-Planck collision operator
DRAFT 1 Fokker-Planck collision operator Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3NP, UK (This version is of 16 April 18) 1. Introduction In these
More informationModelling and numerical methods for the diffusion of impurities in a gas
INERNAIONAL JOURNAL FOR NUMERICAL MEHODS IN FLUIDS Int. J. Numer. Meth. Fluids 6; : 6 [Version: /9/8 v.] Modelling and numerical methods for the diffusion of impurities in a gas E. Ferrari, L. Pareschi
More informationElectric and Magnetic Forces in Lagrangian and Hamiltonian Formalism
Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Benjamin Hornberger 1/26/1 Phy 55, Classical Electrodynamics, Prof. Goldhaber Lecture notes from Oct. 26, 21 Lecture held by Prof. Weisberger
More informationLecture 6 Gas Kinetic Theory and Boltzmann Equation
GIAN Course on Rarefied & Microscale Gases and Viscoelastic Fluids: a Unified Framework Lecture 6 Gas Kinetic Theory and Boltzmann Equation Feb. 23 rd ~ March 2 nd, 2017 R. S. Myong Gyeongsang National
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 information