Kinetic theory of ions in the magnetic presheath
|
|
- Jasper Jennings
- 5 years ago
- Views:
Transcription
1 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 OX1 3NP, UK 2. CCFE, Culham Science Centre, Abingdon, OX14 3DB, UK 9 th Plasma kinetic working group meeting (Vienna, 2016) 1
2 Introduction Boundary conditions for fluid codes used to simulate the SOL plasma Ones currently used* are obtained using fluid equations: aim is to obtain boundary conditions using a kinetic treatment Could be used with future drift kinetic codes of SOL Interesting problem from a purely theoretical point of view: generalizing gyrokinetics to strongly distorted orbits in the magnetic presheath geometry * J. Loizu, P. Ricci, F.D. Halpern and S. Jolliet, Phys. Plasmas 19, (2012). 2
3 Geometry flux surface ζ divertor target through the looking glass Figure (left) from: P.C. Stangeby, The plasma boundary of magnetic fusion devices (2000) 3
4 Boundary layers Width Estimate Collisional layer αλ mfp 100 mm Magnetic presheath ρ i 0.7 mm Debye sheath λ D 0.02 mm λ D << ρ i << αλ mfp Estimates using data from: F. Militello and W. Fundamenski, Plasma Phys. Control. Fusion, 53, (2011) N.B. wall = divertor target 4
5 Assump+ons Focus on ion gyroradius scale ρ i => derive magnetic presheath equations λ D << ρ i ~ αλ mfp => include collisions but keep quasineutral Angle α satisfies (1º=) 0.02 ~ (m e /m i ) ½ << α << 1 Magnetic field constant and presheath electrostatic Turbulent gradients parallel to the wall inherited from outside magnetic presheath neglect / y~1/l~δ/ρ i and / z~δ/l~δ 2 /ρ i Order / t~δ 2 % << presheath timescale δ% Maximal ordering α~δ<<1 α<<1 and δ=ρ i /l <<1 N.B. l~1cm => δ~0.1 5
6 Assumption: (not too) small angle Q: When does it take an ion and an electron the same time to reach the wall? Ion intersects wall after a single orbit: time~ρ i /v t,i Electrons travel along field at v t,e reaching wall after time d/v t,e Times equal when α ρ i /d= (m e /m i ) In order to assume negatively charged wall, require electrons to reach wall much faster than ions Corresponds to α>> (m e /m i ) 6
7 The zeroth order problem: α=δ=0 Equations of motion of single particle are v x = v y = v z =0 Constants of motion: Ze d (x) m i dx v x + v y Orbit position Perpendicular energy x = x + (1/*)v y U = ½v x 2 + ½v y 2 + Zeφ/m i Total energy U = U + ½v z 2 See also: R.H. Cohen and D.D. Ryutov, Phys. Plasmas 5, 808 (1998) 7
8 Gyrophase The zeroth order motion periodic when magnetic force large enough to make the ion turn Can write v x = σ x [2(U - χ(x))] ½ where σ x =±1 Both x and v x are periodic if particle is trapped around a minimum of the effective potential χ(x) = ½* 2 (x-x ) 2 + Zeφ(x)/m i Modified ion gyrofrequency % mod 2 mod =2 Gyrophase ϕ Z xt x b dx v x =2 Z xt x b p 2(U? dx ' = + mod t = x mod Z x (x)) 1 p x b 2(U? (x 0 )) dx0 8
9 Change of variables Can describe particle motion entirely using new set of variables (x, y, z, v x, v y, v z ) à (x, y, z, U, ϕ, U, σ ) x, U and U are constant and y, z are symmetry directions Define gyroaverage of a quantity as an average over ϕ while holding all other variables fixed h...i ' 1 2 I (...) d' = mod Z xt (...) p x b 2(U? (x)) dx 9
10 Single particle motion in system v x = v y = v z = with α~δ<<1 The exact equations of motion are m (x, y, z, neglect t)+ v 1 y cos Ze m (x, y, z, t) v 1 x cos v z sin Ze small m (x, y, z, t)+ v y sin α Changes in the orbit parameters: ExB drift parallel to wall Orbit drift normal to wall Perpendicular energy change Total energy change D U? E hẏi ' = 1 B ' 1 = v k B D UE h xi ' = v k 1 ' ' 10
11 First order constants: y * and 0 Frame exists in which zeroth order is periodic à => An adiabatic invariant * exists, with <d0/dt> = O(α 2 0) when α~δ<<1 Quantity y * proportional to z- canonical momentum derived by integrating dv z /dt = α*v y to obtain v z = α*(y y * ), with dy * /dt = O(α 2 y * ) Adiabatic invariant* y-star µ = 1 Z xt x b p 2(U? (x))dx y? = y h µi ' =0 1 v k hẏ? i ' =0 replaces U replaces y new change of variables: (x, y, z, U, ϕ, U, σ ) à (x, y *, z, *, ϕ, U, σ ) * R.H. Cohen and D.D. Ryutov, Phys. Plasmas 5, 808 (1998) D.L. Holland, B.D. Fried and G.J. Morales, Phys Fluids B, 5, 1723 (1993) 11
12 Gyrokine+c x h x + hẏi ' F <F> ϕ first order = C[F + D = hc[f ]i ' To solve it, require a form for the potential φ(x,y) (that we will determine) and boundary conditions: 1. F = 0 at yà± (corresponds to outside SOL) 2. F = 0 at U à 3. F = F (y, U, U) for <dx /dt> ϕ < 0 at x à 4. F = 0 for open orbits (ones that intersect wall) 12
13 Geometry flux surface ζ divertor target through the looking glass Figure (left) from: P.C. Stangeby, The plasma boundary of magnetic fusion devices (2000) 13
14 Gyrokine+c x +ẏ? F <F> ϕ first order @µ + '@F = C[F ] small small Boundary conditions are: F 0 = F 0 (y *, 0, U) for <dx /dt> ϕ < 0 at x à F 0 = 0 for open orbits 14
15 Quasineutrality The assumption α>>(m e /m i ) ½ => Boltzmann electrons n e =n e, exp[e(φ-φ )/T e ] λ D << ρ i => quasineutrality holds throughout magnetic presheath To obtain ion density, integrate F 0 in velocity space using (x, U, U) instead of (v x, v y, v z )* e ( (x, y) 1 (y)) n e1 (y)exp = Z X k=±1 Z 1 x m (x,y) Z d x T e M ( x,y) (x, x,y) 2 du p? 2(U? (x, x)) Z 1 F ( x, y, U?,U, k) p du U? 2(U U? ) Quasineutrality (above) and the gyrokinetic equation allow to solve for F and φ self consistently in the magnetic presheath * M.J. Gerver, S.E. Parker and K. Theilhaber, Phys. Fluids B, 2, 1069 (1990) 15
16 Collisionless magnetic presheath If ρ i << αλ mfp, gyrokinetic equation becomes (using x =0 The solution is ( therefore F 1 (y?,µ,u) for closed orbits, F = 0 for open orbits. An iteration scheme with quasineutrality would allow to determine the self-consistent φ(x,y) e ( (x, y) 1 (y)) n e1 (y)exp = Z X k=±1 Z 1 x m (x,y) Z d x T e M ( x,y) (x, x,y) will improve on this! 2 du p? 2(U? (x, x)) Z 1 U? F 1 (y?,µ,u)) p 2(U U? ) du See also: R.H. Cohen and D.D. Ryutov, Phys. Plasmas 5, 808 (1998). 16
17 Open orbits Orbit opening 17
18 Why are open orbits important close to the wall? At exit of quasineutral presheath, electric field diverges (breakdown of quasineutrality) but potential does not (φ ~ x)* Expect potential of this form à Electric force always overcomes magnetic force (=Zev y B) close to x=0 =>Effective potential always has maximum near wall à No closed orbits => need open orbit density * K.-U. Riemann, Journal of Physics D: Applied Physics 24, 493 (1991) 18
19 Quasineutrality without open orbits Closed orbit density n i,closed goes to zero at wall Electron density n e ~ n exp[e(φφ )/T e ] ~ n exp[e φ MPS /T e ] Gives φ MPS =- Density of ion open orbits required to obtain finite potential jump Zn i,open n e near wall n i,closed/e / n x / ρ i 19
20 χ(x) Open orbit density Conservation of distribution function F in phase space f open v x U? = M dv x dv y dv z = Density is Z 1 n i,open = 3 du F (U, U? = M, x) v xi v k 0 Z 1 0 d x Z x M x t U? U? = M dx i dxdu x M dx i (x i x M ) v x v xi F 1 ( x, U? = M,U) x U =χ M x i v i x M x t x x i Confirmed by calculation of velocity corrections near X-point Divergence at v x =0: calculate correction at X point v x = Z qv 2M +2( M (x)) with xm V 2M =2 2 s x v k p M x i 2( M x M y (s)) ds 20
21 Preliminary numerical results n i/e /n Green: electron density Blue: ion density Without With open orbits n i/e /n x/ρ i Analytically, n i,open ~ αn => φ MPS ~ (T/e) ln α Numerical results shown in plots obtained using: α=0.1 F ~v 2 exp(-m(v 2 +2*B)/2T i ) φ=1.05*ln(1/α)*exp(-0.19 x-0.49x) 21 x/ρ i
22 Conclusions Derived gyrokinetic equations of ions in magnetic presheath Assumed small magnetic field to wall angle α, small gradients parallel to the wall (δ<<1) and constant B field Assumed electron repelling wall => Boltzmann electrons Derived form of quasineutrality Proposed and currently applying iteration scheme that could solve for collisionless magnetic presheath Solution valid to lowest order in α and δ Open orbits are important close to the wall In progress: quantifying effects of open orbits near the wall + including open orbit density in numerical work Future work: study the purely collisional layer αλ mfp wide 22
23 BACKUP SLIDES 23
24 Assumption: turbulent gradients We include weak gradients parallel to the wall of the electrostatic potential φ and the ion and electron densities, due to the width of turbulent structures l=ρ i /δ with δ~α<<1 These gradients are in the y direction, across the magnetic field The direction (almost) parallel to the magnetic field is associated with even smaller gradients because turbulent structures are elongated along this direction: v z / z~α 2 * 0 Time derivatives are ordered very small: / t~α 2 * 0 neglect Figure: The elongation of a turbulent structure along the magnetic field causes the gradients in the z direction to be very small. The characteristic length scale in the y direction (out of the page) is l. 24
25 Assumption: (not too) small angle When does it take the same time for an ion and an electron to reach the wall? Suppose we are very close to the wall It takes an ion a single orbit to intersect the wall ~ρ i /v t,i Electrons have a much smaller Larmor radius, so they travel parallel to the field at v t,e and reach the wall after a time d/v t,e Times equal when α ρ i /d= (m e /m i ) 25 Figure: A schematic which shows that, when very close to the wall, ions intersect the wall during their gyromotion while electrons are tied much more closely to the field line and have to drift along B with speed ~v t,e much faster than the characteristic ion speed.. At small enough angle α, the ions reach the wall more quickly than thefaster-moving electrons, which travel almost parallel to the wall.
26 Assumption: (not too) small angle (2) In order to assume a negatively charged wall, we require electrons to reach it much faster than ions Corresponds to α>> (m e /m i ) Provided this is satisfied, almost all of the electrons are repelled by the wall We can therefore assume that electrons are in equilibrium and are Boltzmann distributed n e =n e, exp[e(φ-φ )/T e ] Figure: A schematic which shows that, when very close to the wall, ions intersect the wall during their gyromotion while electrons are tied much more closely to the field line and have to drift along B with speed ~v t,e much faster than the characteristic ion speed.. At small enough angle α, the ions reach the wall more quickly than thefaster-moving electrons, which travel almost parallel to the wall. 26
27 dv y /dt = - %v x => v y the orbit position Orbit parameters = %(x - x) where x is Perpendicular energy U and total energy U are conserved Orbit position Perpendicular energy x = x + (1/*)v y U = ½v x 2 + ½v y 2 + Zeφ/m i Total energy U = U + ½v z 2 27
28 Kine+c equa+on change variables small small small zeroth order ~%F expand F = F 0 + F 1 + O(α 2 F) to get => F 0 = <F 0 > ϕ 28
29 Changes in the orbit parameters ExB drift parallel to wall Orbit drift normal to wall Perpendicular energy change Total energy change hẏi ' = 1 B h xi ' = v k 1 B D U? E ' 1 = v k B D UE @y ' ' 29
30 Gyrokine+c equa+on expand F and recover F 0 = <F 0 > first order ~α%f and gyroaverage small small Boundary conditions are: F 0 = F 0 (y *, 0, U) for <dx /dt> ϕ < 0 at x à F 0 = 0 for open orbits 30
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
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 informationGyrokinetic simulations of magnetic fusion plasmas
Gyrokinetic simulations of magnetic fusion plasmas Tutorial 2 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr
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 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 informationIntrinsic rotation due to non- Maxwellian equilibria in tokamak plasmas. Jungpyo (J.P.) Lee (Part 1) Michael Barnes (Part 2) Felix I.
Intrinsic rotation due to non- Maxwellian equilibria in tokamak plasmas Jungpyo (J.P.) Lee (Part 1) Michael Barnes (Part 2) Felix I. Parra MIT Plasma Science & Fusion Center. 1 Outlines Introduction to
More informationMultiscale, multiphysics modeling of turbulent transport and heating in collisionless, magnetized plasmas
Multiscale, multiphysics modeling of turbulent transport and heating in collisionless, magnetized plasmas Michael Barnes Plasma Science & Fusion Center Massachusetts Institute of Technology Collaborators:
More informationTURBULENT TRANSPORT THEORY
ASDEX Upgrade Max-Planck-Institut für Plasmaphysik TURBULENT TRANSPORT THEORY C. Angioni GYRO, J. Candy and R.E. Waltz, GA The problem of Transport Transport is the physics subject which studies the physical
More informationNeoclassical transport
Neoclassical transport Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 28 th January 2013 Dr Ben Dudson Magnetic Confinement Fusion (1 of 19) Last time Toroidal devices
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 informationConfinement of toroidal non-neutral plasma
10th International Workshop on Non-neutral Plasmas 28 August 2012, Greifswald, Germany 1/20 Confinement of toroidal non-neutral plasma in magnetic dipole RT-1: Magnetospheric plasma experiment Visualized
More informationArbiTER studies of filamentary structures in the SOL of spherical tokamaks
ArbiTER studies of filamentary structures in the SOL of spherical tokamaks D. A. Baver, J. R. Myra, Research Corporation F. Scotti, Lawrence Livermore National Laboratory S. J. Zweben, Princeton Plasma
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 informationEdge Momentum Transport by Neutrals
1 TH/P3-18 Edge Momentum Transport by Neutrals J.T. Omotani 1, S.L. Newton 1,2, I. Pusztai 1 and T. Fülöp 1 1 Department of Physics, Chalmers University of Technology, 41296 Gothenburg, Sweden 2 CCFE,
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 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 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 informationBlob sizes and velocities in the Alcator C-Mod scrapeoff
P1-59 Blob sizes and velocities in the Alcator C-Mod scrapeoff layer R. Kube a,b,*, O. E. Garcia a,b, B. LaBombard b, J. L. Terry b, S. J. Zweben c a Department of Physics and Technology, University of
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 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 informationModelling of plasma edge turbulence with neutrals
Modelling of plasma edge turbulence with neutrals Ben Dudson 1 1 York Plasma Institute, Department of Physics, University of York, Heslington, York YO1 5DD, UK 7 th IAEA TM on Plasma Instabilities 4-6
More informationLimitations of gyrokinetics on transport time scales
PSFC/JA-08-16 Limitations of gyrokinetics on transport time scales Parra, F.I. and Catto, P.J February 2008 Plasma Science and Fusion Center Massachusetts Institute of Technology Cambridge MA 02139 USA
More informationGyrokinetic Field Theory without Lie Transforms
ASDEX Upgrade Gyrokinetic Field Theory without Lie Transforms for realisable tokamak codes B. Scott Max Planck Institut für Plasmaphysik Euratom Association D-85748 Garching, Germany Lectures in Gyrokinetic
More informationA kinetic neutral atom model for tokamak scrape-off layer tubulence simulations. Christoph Wersal, Paolo Ricci, Federico Halpern, Fabio Riva
A kinetic neutral atom model for tokamak scrape-off layer tubulence simulations Christoph Wersal, Paolo Ricci, Federico Halpern, Fabio Riva CRPP - EPFL SPS Annual Meeting 2014 02.07.2014 CRPP The tokamak
More informationTurbulence in Tokamak Plasmas
ASDEX Upgrade Turbulence in Tokamak Plasmas basic properties and typical results B. Scott Max Planck Institut für Plasmaphysik Euratom Association D-85748 Garching, Germany Uni Innsbruck, Nov 2011 Basics
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 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 informationLectures on basic plasma physics: Hamiltonian mechanics of charged particle motion
Lectures on basic plasma physics: Hamiltonian mechanics of charged particle motion Department of applied physics, Aalto University March 8, 2016 Hamiltonian versus Newtonian mechanics Newtonian mechanics:
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 informationTurbulent Transport of Toroidal Angular Momentum in Low Flow Gyrokinetics
PSFC/JA-09-27 Turbulent Transport of Toroidal Angular Momentum in Low Flow Gyrokinetics Felix I. Parra and P.J. Catto September 2009 Plasma Science and Fusion Center Massachusetts Institute of Technology
More informationPlasma Descriptions I: Kinetic, Two-Fluid
CHAPTER 5. PLASMA DESCRIPTIONS I: KINETIC, TWO-FLUID 1 Chapter 5 Plasma Descriptions I: Kinetic, Two-Fluid Descriptions of plasmas are obtained from extensions of the kinetic theory of gases and the hydrodynamics
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 informationTH/P4-9. T. Takizuka 1), K. Shimizu 1), N. Hayashi 1), M. Hosokawa 2), M. Yagi 3)
1 Two-dimensional Full Particle Simulation of the Flow Patterns in the Scrape-off-layer Plasma for Upper- and Lower- Null Point Divertor Configurations in Tokamaks T. Takizuka 1), K. Shimizu 1), N. Hayashi
More informationOn the locality of parallel transport of heat carrying electrons in the SOL
P1-068 On the locality of parallel transport of heat carrying electrons in the SOL A.V. Chankin* and D.P. Coster Max-Planck-Institut für Pasmaphysik, 85748 Garching, Germany Abstract A continuum Vlasov-Fokker-Planck
More informationA.G. PEETERS UNIVERSITY OF BAYREUTH
IN MEMORIAM GRIGORY PEREVERZEV A.G. PEETERS UNIVERSITY OF BAYREUTH ESF Workshop (Garching 2013) Research areas Grigory Pereverzev. Current drive in magnetized plasmas Transport (ASTRA transport code) Wave
More informationGas Dynamics: Basic Equations, Waves and Shocks
Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks Susanne Höfner Susanne.Hoefner@fysast.uu.se Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks
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 informationCollision Processes. n n The solution is 0 exp x/ mfp
Collision Processes Collisions mediate the transfer of energy and momentum between various species in a plasma, and as we shall see later, allow a treatment of highly ionized plasma as a single conducting
More informationDirect drive by cyclotron heating can explain spontaneous rotation in tokamaks
Direct drive by cyclotron heating can explain spontaneous rotation in tokamaks J. W. Van Dam and L.-J. Zheng Institute for Fusion Studies University of Texas at Austin 12th US-EU Transport Task Force Annual
More informationMHD WAVES AND GLOBAL ALFVÉN EIGENMODES
MHD WVES ND GLOBL LFVÉN EIGENMODES S.E. Sharapov Euratom/CCFE Fusion ssociation, Culham Science Centre, bingdon, Oxfordshire OX14 3DB, UK S.E.Sharapov, Lecture 3, ustralian National University, Canberra,
More informationVlasov simulations of plasma-wall interactions in a magnetized and weakly collisional plasma
PHYSICS OF PLASMAS 13, 083504 2006 Vlasov simulations of plasma-wall interactions in a magnetized and weakly collisional plasma S. Devaux a Laboratoire de Physique des Milieux Ionisés et Applications,
More informationGreg Hammett Imperial College, London & Princeton Plasma Physics Lab With major contributions from:
Greg Hammett Imperial College, London & Princeton Plasma Physics Lab With major contributions from: Steve Cowley (Imperial College) Bill Dorland (Imperial College) Eliot Quataert (Berkeley) LMS Durham
More information0 Magnetically Confined Plasma
0 Magnetically Confined Plasma 0.1 Particle Motion in Prescribed Fields The equation of motion for species s (= e, i) is written as d v ( s m s dt = q s E + vs B). The motion in a constant magnetic field
More informationToroidal confinement devices
Toroidal confinement devices Dr Ben Dudson Department of Physics, University of York, Heslington, York YO10 5DD, UK 24 th January 2014 Dr Ben Dudson Magnetic Confinement Fusion (1 of 20) Last time... Power
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 informationstable sheath Formation in Magnetized Plasma
J. Plasma Fusion Res. SERES, Vol.4 (2001) 578-582 stable sheath Formation in Magnetized Plasma T9MTA YuKihiro*, TAKAYAMA AriMiChi, TAKAMARU HiSANOTi ANd SATO TETSUYA National nstitute.for Fusion science,
More informationTurbulent Transport Analysis of JET H-mode and Hybrid Plasmas using QuaLiKiz, TGLF and GLF23
EFDA JET CP(1)/ B. Baiocchi, J. Garcia, M. Beurkens, C. Bourdelle, F. Crisanti, C. Giroud, J. Hobirk, F. Imbeaux, I. Nunes, EU-ITM ITER Scenario Modelling group and JET EFDA contributors Turbulent Transport
More informationPhysics 2210 Fall smartphysics Exam 3 Review smartphysics units /04/2015
Physics 22 Fall 25 smartphysics Exam 3 Review smartphysics units -3 /4/25 Review Problem The figure shown extends from x = to x = and is bounded on the left by the y-axis, on the bottom by the x-axis,
More informationHistory of PARASOL! T. Takizuka! Graduate School of Engineering, Osaka University!! PARASOL was developed at Japan Atomic Energy Agency!
History of PARASOL! T. Takizuka!!! Graduate School of Engineering, Osaka University!! PARASOL was developed at Japan Atomic Energy Agency! OSAKA UNIVERSITY! 20th NEXT Meeting, Kyoto Terrsa, Kyoto, 13-14
More informationTH/P6-14 Integrated particle simulation of neoclassical and turbulence physics in the tokamak pedestal/edge region using XGC a)
1 TH/P6-14 Integrated particle simulation of neoclassical and turbulence physics in the tokamak pedestal/edge region using XGC a) 1 Chang, C.S., 1 Ku, S., 2 Adams M., 3 D Azevedo, G., 4 Chen, Y., 5 Cummings,
More informationDiffusion equation, flux, diffusion coefficient scaling. Diffusion in fully ionized plasma vs. weakly ionized plasma. n => Coulomb collision frequency
Last Time Diffusion in plasma: the main reason why we need to control it (i.e. using magnetic field) Diffusion equation, flux, diffusion coefficient scaling o o t nx,t Dn D2 nx,t o D ~ L 2 T Diffusion
More informationFine-Scale Zonal Flow Suppression of Electron Temperature Gradient Turbulence
Fine-Scale Zonal Flow Suppression of Electron Temperature Gradient Turbulence S.E. Parker, J.J. Kohut, Y. Chen, Z. Lin, F.L. Hinton and W.W. Lee Center for Integrated Plasma Studies, University of Colorado,
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 informationPower Exhaust on JET: An Overview of Dedicated Experiments
Power Exhaust on JET: An Overview of Dedicated Experiments W.Fundamenski, P.Andrew, T.Eich 1, G.F.Matthews, R.A.Pitts 2, V.Riccardo, W.Sailer 3, S.Sipila 4 and JET EFDA contributors 5 Euratom/UKAEA Fusion
More informationResistive MHD, reconnection and resistive tearing modes
DRAFT 1 Resistive MHD, reconnection and resistive tearing modes Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3NP, UK (This version is of 6 May 18 1. Introduction
More informationStability of a plasma confined in a dipole field
PHYSICS OF PLASMAS VOLUME 5, NUMBER 10 OCTOBER 1998 Stability of a plasma confined in a dipole field Plasma Fusion Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 Received
More informationMotion of Charged Particles in Fields
Chapter Motion of Charged Particles in Fields Plasmas are complicated because motions of electrons and ions are determined by the electric and magnetic fields but also change the fields by the currents
More informationGyrokinetic simulations with GYSELA: Main current issues in physics & numerics
Gyrokinetic simulations with GYSELA: Main current issues in physics & numerics Y. Sarazin, Y. Asahi 2, N. Bouzat, G. Dif-Pradalier, P. Donnel, C. Ehrlacher, C. Emeriau 3, X. Garbet, Ph. Ghendrih, V. Grandgirard,
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 informationMicrotearing Simulations in the Madison Symmetric Torus
Microtearing Simulations in the Madison Symmetric Torus D. Carmody, P.W. Terry, M.J. Pueschel - University of Wisconsin - Madison dcarmody@wisc.edu APS DPP 22 Overview PPCD discharges in MST have lower
More informationVerification & Validation: application to the TORPEX basic plasma physics experiment
Verification & Validation: application to the TORPEX basic plasma physics experiment Paolo Ricci F. Avino, A. Bovet, A. Fasoli, I. Furno, S. Jolliet, F. Halpern, J. Loizu, A. Mosetto, F. Riva, C. Theiler,
More informationIntroduction to Plasma Physics
Introduction to Plasma Physics Hartmut Zohm Max-Planck-Institut für Plasmaphysik 85748 Garching DPG Advanced Physics School The Physics of ITER Bad Honnef, 22.09.2014 A simplistic view on a Fusion Power
More informationGyrokinetic Simulations of Tearing Instability
Gyrokinetic Simulations of Tearing Instability July 6, 2009 R. NUMATA A,, W. Dorland A, N. F. Loureiro B, B. N. Rogers C, A. A. Schekochihin D, T. Tatsuno A rnumata@umd.edu A) Center for Multiscale Plasma
More informationLectures on basic plasma physics: Introduction
Lectures on basic plasma physics: Introduction Department of applied physics, Aalto University Compiled: January 13, 2016 Definition of a plasma Layout 1 Definition of a plasma 2 Basic plasma parameters
More informationSingle particle motion
Single particle motion Plasma is a collection of a very large number of charged particles moving in, and giving rise to, electromagnetic fields. Before going to the statistical descriptions, let us learn
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 informationRelativistic magnetohydrodynamics. Abstract
Relativistic magnetohydrodynamics R. D. Hazeltine and S. M. Mahajan Institute for Fusion Studies, The University of Texas, Austin, Texas 78712 (October 19, 2000) Abstract The lowest-order description of
More informationSMR/ Summer College on Plasma Physics. 30 July - 24 August, Introduction to Magnetic Island Theory.
SMR/1856-1 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 Introduction to Magnetic Island Theory. R. Fitzpatrick Inst. for Fusion Studies University of Texas at Austin USA Introduction
More informationPlasma-neutrals transport modeling of the ORNL plasma-materials test stand target cell
Plasma-neutrals transport modeling of the ORNL plasma-materials test stand target cell J.M. Canik, L.W. Owen, Y.K.M. Peng, J. Rapp, R.H. Goulding Oak Ridge National Laboratory ORNL is developing a helicon-based
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 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 informationPlasma physics in noninertial frames
Plasma physics in noninertial frames PHYSICS OF PLASMAS 16, 092506 2009 A. Thyagaraja and K. G. McClements EURATOM/UKAEA Fusion Association, Culham Science Centre, Abingdon OX14 3DB, United Kingdom Received
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 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 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 Module No. # 01 Lecture No. # 22 Adiabatic Invariance of Magnetic Moment and Mirror Confinement Today, we
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 informationGyrokinetic neoclassical study of the bootstrap current in the tokamak edge pedestal with fully non-linear Coulomb collisions
Physics of Plasmas Gyrokinetic neoclassical study of the bootstrap current in the tokamak edge pedestal with fully non-linear Coulomb collisions Robert Hager 1, a) 1, b) and C.S. Chang Princeton Plasma
More informationNon-perturbative statistical theory of intermittency in ITG drift wave turbulence with zonal flows
Non-perturbative statistical theory of intermittency in ITG drift wave turbulence with zonal flows Johan Anderson 1 and Eun-jin Kim University of Sheffield Department of Applied Mathematics Hicks Building,
More informationSensitivity of Detachment Extent to Magnetic Configuration and External Parameters. September 2015
PSFC/JA-15-64 Sensitivity of Detachment Extent to Magnetic Configuration and External Parameters Bruce Lipschultz 1, Felix I. Parra 2, 3 and Ian H. Hutchinson 4 1 York Plasma Institute, University of York,
More informationarxiv: v1 [physics.plasm-ph] 2 Jun 2015
arxiv:1506.00863v1 [physics.plasm-ph] 2 Jun 2015 Turbulent momentum transport due to neoclassical flows Jungpyo Lee Plasma Science and Fusion Center, Massachusetts Institute of Technology, MA, USA E-mail:
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 informationOperational Phase Space of the Edge Plasma in Alcator C-Mod
Operational Phase Space of the Edge Plasma in B. LaBombard, T. Biewer, M. Greenwald, J.W. Hughes B. Lipschultz, N. Smick, J.L. Terry, Team Contributed talk RO.00008 Presented at the 47th Annual Meeting
More informationarxiv: v1 [physics.plasm-ph] 10 Feb 2017
Under consideration for publication in J. Plasma Phys. 1 arxiv:17.35v1 [physics.plasm-ph] 1 Feb 17 Gyrokinetic continuum simulation of turbulence in a straight open-field-line plasma E. L. Shi 1, G. W.
More informationCoupled radius-energy turbulent transport of alpha particles
Coupled radius-energy turbulent transport of alpha particles George Wilkie, Matt Landreman, Ian Abel, William Dorland 24 July 2015 Plasma kinetics working group WPI, Vienna Wilkie (Maryland) Coupled transport
More informationGyrokinetics an efficient framework for studying turbulence and reconnection in magnetized plasmas
Frank Jenko Gyrokinetics an efficient framework for studying turbulence and reconnection in magnetized plasmas Max-Planck-Institut für Plasmaphysik, Garching Workshop on Vlasov-Maxwell Kinetics WPI, Vienna,
More informationGyrokinetic theory for particle transport in fusion plasmas
Gyrokinetic theory for particle transport in fusion plasmas Matteo Valerio Falessi 1,2, Fulvio Zonca 3 1 INFN - Sezione di Roma Tre, Via della Vasca Navale, 84 (00146) Roma (Roma), Italy 2 Dipartimento
More informationION THERMAL CONDUCTIVITY IN TORSATRONS. R. E. Potok, P. A. Politzer, and L. M. Lidsky. April 1980 PFC/JA-80-10
ION THERMAL CONDUCTIVITY IN TORSATRONS R. E. Potok, P. A. Politzer, and L. M. Lidsky April 1980 PFC/JA-80-10 ION THERMAL CONDUCTIVITY IN TORSATRONS R.E. Potok, P.A. Politzer, and L.M. Lidsky Plasma Fusion
More informationMagnetic Deflection of Ionized Target Ions
Magnetic Deflection of Ionized Target Ions D. V. Rose, A. E. Robson, J. D. Sethian, D. R. Welch, and R. E. Clark March 3, 005 HAPL Meeting, NRL Solid wall, magnetic deflection 1. Cusp magnetic field imposed
More informationCHARGED PARTICLE MOTION IN CONSTANT AND UNIFORM ELECTROMAGNETIC FIELDS
CHARGED PARTICLE MOTION IN CONSTANT AND UNIFORM ELECTROMAGNETIC FIELDS In this and in the following two chapters we investigate the motion of charged particles in the presence of electric and magnetic
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 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 informationGuiding Center Orbit Studies in a Tokamak Edge Geometry Employing Boozer and Cartesian Coordinates
Contrib. Plasma Phys. 48, No. -3, 4 8 (8) / DOI./ctpp.839 Guiding Center Orbit Studies in a Tokamak Edge Geometry Employing Boozer and Cartesian Coordinates Y. Nishimura,Y.Xiao,and Z. Lin Department of
More informationParticle-in-cell simulations of electron transport from plasma turbulence: recent progress in gyrokinetic particle simulations of turbulent plasmas
Institute of Physics Publishing Journal of Physics: Conference Series 16 (25 16 24 doi:1.188/1742-6596/16/1/2 SciDAC 25 Particle-in-cell simulations of electron transport from plasma turbulence: recent
More informationModelling of JT-60U Detached Divertor Plasma using SONIC code
J. Plasma Fusion Res. SERIES, Vol. 9 (2010) Modelling of JT-60U Detached Divertor Plasma using SONIC code Kazuo HOSHINO, Katsuhiro SHIMIZU, Tomonori TAKIZUKA, Nobuyuki ASAKURA and Tomohide NAKANO Japan
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 informationGyrokinetic Theory and Dynamics of the Tokamak Edge
ASDEX Upgrade Gyrokinetic Theory and Dynamics of the Tokamak Edge B. Scott Max Planck Institut für Plasmaphysik D-85748 Garching, Germany PET-15, Sep 2015 these slides: basic processes in the dynamics
More informationGyrokinetic Transport Driven by Energetic Particle Modes
Gyrokinetic Transport Driven by Energetic Particle Modes by Eric Bass (General Atomics) Collaborators: Ron Waltz, Ming Chu GSEP Workshop General Atomics August 10, 2009 Outline I. Background Alfvén (TAE/EPM)
More informationHamiltonian and Non-Hamiltonian Reductions of Charged Particle Dynamics: Diffusion and Self-Organization
NNP2017 11 th July 2017 Lawrence University Hamiltonian and Non-Hamiltonian Reductions of Charged Particle Dynamics: Diffusion and Self-Organization N. Sato and Z. Yoshida Graduate School of Frontier Sciences
More informationarxiv: v1 [physics.plasm-ph] 3 Apr 2011
A comparison of Vlasov with drift kinetic and gyrokinetic theories arxiv:1104.0427v1 [physics.plasm-ph] 3 Apr 2011 H. Tasso 1, G. N. Throumoulopoulos 2 1 Max-Planck-Institut für Plasmaphysik, Euratom Association,
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 informationAdvances in stellarator gyrokinetics
Advances in stellarator gyrokinetics Per Helander and T. Bird, F. Jenko, R. Kleiber, G.G. Plunk, J.H.E. Proll, J. Riemann, P. Xanthopoulos 1 Background Wendelstein 7-X will start experiments in 2015 optimised
More information