Methods of plasma description Particle motion in external fields

Similar documents
Single Particle Motion

Motion of Charged Particles in Fields

Lectures on basic plasma physics: Hamiltonian mechanics of charged particle motion

Single particle motion

Radiation belt particle dynamics

Charged particle motion in external fields

Single Particle Motion

CHARGED PARTICLE MOTION IN CONSTANT AND UNIFORM ELECTROMAGNETIC FIELDS

Mechanics Departmental Exam Last updated November 2013

Waves in plasma. Denis Gialis

10. The dimensional formula for c) 6% d) 7%

Symmetries 2 - Rotations in Space

Physics: Dr. F. Wilhelm E:\Excel files\230 lecture\230 tests\final 230 F07 Practice.doc page 1 of 9

Magnetospheric Physics - Homework, 02/21/ Local expansion of B: Consider the matrix B = (1)

Chapter 2 Single-Particle Motions

Motion of a charged particle in an EM field

Newton s Laws of Motion, Energy and Oscillations

SOLUTIONS FOR THEORETICAL COMPETITION Theoretical Question 1 (10 points) 1A (3.5 points)

Neoclassical transport

5. Charge Particle Motion in Fields

Exam 2 Solutions. Applying the junction rule: i 1 Applying the loop rule to the left loop (LL), right loop (RL), and the full loop (FL) gives:

Oscillatory Motion. Solutions of Selected Problems

Basics of electromagnetic response of materials

Fundamentals of Applied Electromagnetics. Chapter 2 - Vector Analysis

= qe cos(kz ωt). ω sin(αt) This is further averaged over the distribution of initial ( ) α + ω. = 2m k P g(α) sin(αt) g(α) = g(0) + αg (0) + α2 2 g +

Louisiana State University Physics 2102, Exam 3 April 2nd, 2009.

Part III. Interaction with Single Electrons - Plane Wave Orbits

Single particle motion and trapped particles

Plasmas as fluids. S.M.Lea. January 2007

Lecture 04. Curl and Divergence

Single Particle Motion in a Magnetized Plasma

Oscillations. Tacoma Narrow Bridge: Example of Torsional Oscillation

PY3107 Experimental Physics II

13.1 Ion Acoustic Soliton and Shock Wave

cos 6 λ m sin 2 λ m Mirror Point latitude Equatorial Pitch Angle Figure 5.1: Mirror point latitude as function of equatorial pitch angle.

Foundation Course Unit. Fields and Waves. SOLUTIONS. Examiners: Dr. A. G. Polnarev Prof. F. Di Lodovico

Translational vs Rotational. m x. Connection Δ = = = = = = Δ = = = = = = Δ =Δ = = = = = 2 / 1/2. Work

TURBULENT TRANSPORT THEORY

FOUNDATION STUDIES EXAMINATIONS June PHYSICS Semester One February Main

Accelerator School Transverse Beam Dynamics-2. V. S. Pandit

Solutions to PS 6 Physics 401a

Distance travelled time taken and if the particle is a distance s(t) along the x-axis, then its instantaneous speed is:

The Physics of Field-Aligned Currents

y=1/4 x x=4y y=x 3 x=y 1/3 Example: 3.1 (1/2, 1/8) (1/2, 1/8) Find the area in the positive quadrant bounded by y = 1 x and y = x3

Chapter 15 - Oscillations

Physics GRE: Electromagnetism. G. J. Loges 1. University of Rochester Dept. of Physics & Astronomy. xkcd.com/567/

Practice Test 3. Multiple Choice Identify the choice that best completes the statement or answers the question.

Effects of an Inhomogeneous Magnetic Field (E =0)

(a) Use washer cross sections: a washer has

Summer College on Plasma Physics August Introduction to Nonlinear Gyrokinetic Theory

FOUNDATION STUDIES EXAMINATIONS September 2009

SAFE HANDS & IIT-ian's PACE EDT-12 (JEE) SOLUTIONS

Pan Pearl River Delta Physics Olympiad 2005

Classical Mechanics Comprehensive Exam Solution

PHYSICS PART II SECTION- I. Straight objective Type

Gyrokinetic simulations of magnetic fusion plasmas

THE UNIVERSITY OF SYDNEY FACULTY OF SCIENCE INTERMEDIATE PHYSICS PHYS 2912 PHYSICS 2B (ADVANCED) ALL QUESTIONS HAVE THE VALUE SHOWN

E = K + U. p mv. p i = p f. F dt = p. J t 1. a r = v2. F c = m v2. s = rθ. a t = rα. r 2 dm i. m i r 2 i. I ring = MR 2.

PHY 5246: Theoretical Dynamics, Fall Assignment # 9, Solutions. y CM (θ = 0) = 2 ρ m

Harmonic Oscillator - Model Systems

More accurate ion flow coefficients in the Pfirsch Schluter regime

2 Particle dynamics. The particle is called free body if it is no in interaction with any other body or field.

FYS 3120: Classical Mechanics and Electrodynamics

Observations on the ponderomotive force

Physics 2101 S c e t c i cti n o 3 n 3 March 31st Announcements: Quiz today about Ch. 14 Class Website:

Solutions Assignment 2

Problem Set 4 Solutions

Laser-plasma particle accelerators. Zulfikar Najmudin

1300 (W/m 2 ) (V/cm) = 275 (V/m) (A/cm) = (A/m). E = 990 (V/m), H = 2.63 (A/m).

ENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 3 Solutions [Multiple Integration; Lines of Force]

Lecture 2: Plasma particles with E and B fields

Basic Mathematics and Units

FOUNDATION STUDIES EXAMINATIONS June PHYSICS Semester One February Main

Large scale flows and coherent structure phenomena in flute turbulence

Integrals along a curve in space. (Sect. 16.1)

r r 1 r r 1 2 = q 1 p = qd and it points from the negative charge to the positive charge.

Final Exam. June 10, 2008, 1:00pm

The distance of the object from the equilibrium position is m.

Apr 14, Calculus with Algebra and Trigonometry II Lecture 20More physics Aprapplications

is conserved, calculating E both at θ = 0 and θ = π/2 we find that this happens for a value ω = ω given by: 2g

Classical Mechanics. FIG. 1. Figure for (a), (b) and (c). FIG. 2. Figure for (d) and (e).

arxiv: v2 [math-ph] 14 Apr 2008

Equations involving an unknown function and its derivatives. Physical laws encoded in differential equations INTRODUCTION TO DIFFERENTIAL EQUATIONS

Newton's Laws You should be able to state these laws using both words and equations.

The curious properties of spin

Department of Physics

Chapter 14 Periodic Motion

The Calculus of Vec- tors

EXERCISES Chapter 15: Multiple Integrals. Evaluating Integrals in Cylindrical Coordinates

Figure 1 A) 2.3 V B) +2.3 V C) +3.6 V D) 1.1 V E) +1.1 V Q2. The current in the 12- Ω resistor shown in the circuit of Figure 2 is:

NST1A: Mathematics II (Course A) End of Course Summary, Lent 2011

Waves in plasmas. S.M.Lea

APPENDIX Z. USEFUL FORMULAS 1. Appendix Z. Useful Formulas. DRAFT 13:41 June 30, 2006 c J.D Callen, Fundamentals of Plasma Physics

Equipotential surfaces

Cluster fusion in a high magnetic field

Lecture XXVI. Morris Swartz Dept. of Physics and Astronomy Johns Hopkins University November 5, 2003

Indian National Physics Olympiad 2017

High Energy Astrophysics

. (70.1) r r. / r. Substituting, we have the following equation for f:

Reminders: Show your work! As appropriate, include references on your submitted version. Write legibly!

Transcription:

Methods of plasma description (suggested reading D.R. Nicholson, chap., Chen chap., 8.4) Charged particle motion in eternal electromagnetic (elmg) fields Charged particle motion in self-consistent elmg fields Kinetic equations Fluid (hdrodnamic) equations Particle motion in eternal fields A) Homogeneous fields dv a) E = m qv dt = mv = qv mv = qv v q = m v v - fluids - 1 fluid - magnetohdrodnamics = zˆ m v = q = m v ω c z q m cclotron frequenc PČ 1

v = v ep( ± iω t+ δ ) v = v ep( iω t), c, c m iωc v = v = ± iv e q t v i e ω i ct = ωc v i c =± e ω ω c t Larmor radius r L v = = ω c mv q µ when kt v = vt = m 1 r L = ( mkt ) q 1/ 1 q q mv mv µ = r Id l = IS= Sn= π n= n T π m q (,,z) gration center DIAMAGNETIC MEDIUM PČ

b) E dv m qe ( v ) dt = + d v = q E ±ω v dt m c E = ( E,, E ) dv dt z d dt v z = q E m = ωcv v = ω v z c q E =± ω E ± m ω = ω + v c cv c v d E E v v + = ωc + dt c v = v i( t ) e ω + δ E c v =± iv e i( ω t + δ) PČ 3

v gc - station. E+ v = E gc= gration center v gc = v E drift in E field here q >, v E = ω C r L v E < ω C r L v E > ω C r L potential energ - 1 ma, 3 min when v(t=) =, then v E = ω c r L = v (ccloid), otherwise trochoid, often case v v T > v E general force e.g. gravit force 1 F mg v f = v g = q q gravit drift different direction for electrons and ions ( ) g j= nm+ m gravit current PČ 4

) Inhomogeneous mv mv F ˆ cf = r = R R R k curvature drift k k (centrifugal force) v R 1 = = q q R Fcf mv Rk k div = rot = curved field cannot be constant a) grad- drift z = z( ) = + linear approimation of field during the particle motion on the Larmor circle F = qv z( ) = qv cos ct ± rl cos ct ( ω ) ( ω ) PČ 5

= + ( r ) +... z = + ( z ) +... cos ωct = 1 mv 1 F =± qv rl = v =± v rl 1 curvature drift and grad- drift often complement each other Rk m Rk 1 v v v v + R = + R qr k b) Magnetic mirrors clindrical smmetr div = k 1 z ( rr ) + = r r z => in the neighborhood of ais r z r z r = z z r PČ 6

v r θ rl mv Fz = qv r = qv = µ = magnetic moment m µ? µ = J S S = πr = π q Invariance of µ µ v L 1 mv µ = (s trajector along the line of force) F ( µ ) = µ J q qωc q = = = T π πm dv d 1 v d m = µ m = µ v = µ dt s dt s dt d 1 1 d 1 d d dµ m m m ( ) dt dt dt dt dt adiabatic invariant v + v = v + µ = µ + µ = = PČ 7

v v = = sin v v θ Where is particle with v from area reflected? 1 1 mv = mv v = v = v sin 1 θ m = = m Rm where R m is mirror ratio, it defines loss cone for θ, < θm particle is not trapped Adiabatic invariant quantit that is conserved during slow spatial and temporal variations of the sstem Classical mechanics during periodic motion, action J = pdq is conserved. PČ 8

Gration motion p= mv; q = π ω c π mv mπ J = mvd = mv sin ( ω ct) dt = = µ µ = const. ω q When is the adiabatic invariant µ not conserved? a) cclotron heating ω ω c, E, oscillates ω << ωc does not hold µ const. b) magnetic pumping varies sinusoidall in time, the invariance of µ is broken during particle collisions with moving magnetic mirror If compression (field increase) occurs during = collision, then partl v v in epansion v is unchanged c) magnetic cusp in the middle = ω c = µ const. c PČ 9

Second adiabatic invariant a, b turnabout points J b = v ds longitudinal invariant a Third adiabatic invariant v,v R, J R k 3 v d - drift in direction of angle ϕ = dl - 3. adiabatic invariant PČ 1

C) Inhomogeneous E E = ˆ cos k E = z ˆ dv ( m = q E( ) + v ) = ± rlcos ωct dt E v = = ωc v ωc cos k( ± rlcos ωct) 1 cosk 1 k r 4 L E 1 1 E ve = 1 kr 1 L = + rl 4 4 Polarization drift (temporall varing E) E E t = zˆ E() t = Et ˆ PČ 11

mv = qe + qv v= v + v ˆ ˆ E+ vp assumption v p = const. m v +v = qet ˆ+ qv qv ˆ+ qv ˆ ( ) E E p mv = qv cclotron rotation mv v ˆ E = q p v p = polarization drift = q ˆ v ˆ set qs E v E E drift Et E E me 1 m v ˆ ˆ E= = v v E = p = = E q q m d Mi ne de ρm de v = E J = ne(v v ) = m + = q dt Z dt dt p p e pi pe e polarization current PČ 1

PONDEROMOTIVE FORCE = low frequenc force acting on charged particles in inhomogeneous high frequenc electromagnetic field Oscillation energ of charged particles in high frequenc field is given b the particle position so it is some kind of potential energ U and there is force F = U that epels charged particles from the area of strong field. Ponderomotive force acts on an dielectric if its dielectric constant depends on densit (electrostriction)!! Ponderomotive force consists of parts force caused b particle oscillation in inhomogeneous electric field and force due to magnetic field action on oscillating particle First, we derive it force caused b particle oscillation in inhomogeneous electric field E of frequenc ω : E = E( r)cosωt mr r = r = qe = qe ( r )cosωt r = r + r 1 PČ 13

We linearize field variations in oscillation area r1 and equations of motion are mr ( + r1) = q E + ( r1 ) E cosωt qe q ( ) q mr = qe cosωt r = cosωt r r cosωt E ( E ) E mω = m = m ω q F ( ) E = E E So low frequenc force is mω 1 1 1 High frequenc magnetic field acts on oscillating particle 1 q = ( r)sinωt = rot E v = Esinωt ω m ω The force is then given b the following epression q 1 q F = qv = E sin ωt E rot E mω = mω PČ 14

The total ponderomotive force is then given b the sum of the forces 1 q ( ) 1 q 1 q Fp = FE + F = E E + E rot E = E = E mω 4 mω 4 mω We have used a real amplitude during the derivation, however, the field ma be phase shifted in general, thus, the general epression includes absolute value of the field amplitude. Low frequenc force acting on the particle thus reads, as follows F p q = 4mω E F p = W 1 1 qe 1 q Wosc = mv = m t = E m ω 4 mω osc cos ω force equal to -gradient of potential energ There eists also high frequenc force of frequenc ω. For field containing frequencies there forces with sum and difference of ω. PČ 15