Hamiltonian and Non-Hamiltonian Reductions of Charged Particle Dynamics: Diffusion and Self-Organization
|
|
- Morris Hutchinson
- 5 years ago
- Views:
Transcription
1 NNP 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 The University of Tokyo
2 INTRODUCTION z θ r Magnetospheres are ubiquitous in the universe. Dipole magnetic fields confine charged particles. Spontaneous confinement (self-organization) is achieved by inward diffusion. 2/38
3 MAGNETOSPHERIC CONFINEMENT OF CHARGED PARTICLES Magnetospheric confinement is a promising source of applications. Advanced fusion reactor based on D- 3 H reaction [1]. Trap for matter and antimatter plasmas [2]. Study of astrophysical systems [3]. The creation of the density gradient seemingly violates the entropy principle. OUTLINE 1. Diffusion in a dipole field: correct expression of diffusion operator in magnetic coordinates. 2. Diffusion in non-integrable magnetic field B B 0. [1] A. Hasegawa et al., A D- 3 H Fusion Reactor Based on a Dipole Magnetic Field, Nuc. Fus. 30, 11 (1990). [2] T. S. Pedersen and A. H. Boozer, Confinement of Nonneutral Plasmas on Magnetic Surfaces, Phys. Rev. Lett. 88, 20 (2002). [3] A. M. Persoon et al., The Plasma Density Distribution in the Inner Region of Saturn s Magnetosphere, J. Geo. Res.: Space Physics 118, (2013). 3/38
4 THE RT-1 MACHINE (I) RT-1 produces a laboratory magnetosphere. A similar device is LDX at MIT. [4] Z. Yoshida et al., Self-organized confinement by magnetic dipole: recent results from RT-1 and theoretical modeling, Plasma Phys. Control. Fus (2012). [5] Z. Yoshida et al., Magnetospheric Vortex Formation: Self-Organized Confinement of Charged Particles, Phys. Rev. Lett. 104, (2010). 4/31
5 THE RT-1 MACHINE (II) 5/31
6 EXPERIMENTAL OBSERVATIONS [6] H. Saitoh et al., Observation of a New High-β and High-Density State of a Magnetospheric Plasma in RT-1, Phys. Plasmas 21, (2014). [7] Y. Kawazura et al., Observation of Particle Acceleration in Laboratory Magnetosphere, Phys. Plasmas 22, (2015). 6/31
7 THE POINT DIPOLE MAGNETIC FIELD B = ψ θ B l ψ FLUX FUNCTION ψ r, z = r 2 r 2 + z 2 3/2 LENGTH OF A FIELD LINE [7] l r, z = 1 2ψ ቈ 1 log ቆ 3 1 rψ 2/3 3 θ [8] N. Sato and Z. Yoshida, A Stochastic Model of Inward Diffusion in Magnetospheric Plasmas, J. Phys. A: Math. Theor (2015). 7/31
8 INWARD DIFFUSION Fickian Diffusion Inward Diffusion Flux-tube ψ It is thought that particle number per flux-tube volume tends to be homogenized [9,10]. [9] A. C. Boxer et al., Turbulent Inward Pinch of Plasma Confined by a Levitated Dipole Magnet, Nature Physics, Vol. 6, (2010). [10] A. Hasegawa, Motion of a Charged Particle and Plasma Equilibrium in a Dipole Magnetic Field, Phys. Scr. T (2005). 8/31
9 THE MECHANISM OF INWARD DIFFUSION In a dipole field the motion of a charged particle is constrained by 3 adiabatic invariants: ψ J μ STRENGTH OF CONSTRAINT Each invariant is preserved on a different time scale: ω d ω b ω c In a dipole plasma the key normal transport mechanism is diffusion by E B drift: v = E B. 9/31
10 DIFFUSION ON A FOLIATED PHASE SPACE Conservation of adiabatic invariants foliates phase space: particle motion is constrained on leaves [11,12]. Electromagnetic fluctuations break the weakest topological constraints. The result is a random process among leaves. A diffusion operator must be formulated on the corresponding invariant measure [8,13]. Entropy should be maximized on the symplectic leaf [14]. [8] N. Sato and Z. Yoshida, A Stochastic Model of Inward Diffusion in Magnetospheric Plasmas, J. Phys. A: Math. Theor (2015). [11] N. Sato, N. Kasaoka, and Z. Yoshida, Thermal Equilibrium of Non-Neutral Plasma in Dipole Magnetic Field, Phys. Plasmas 22 (2015), [12] Z. Yoshida and S. M. Mahajan, Self-Organization in Foliated Phase Space: Construction of a Scale Hierarchy by, Prog. Theor. Exp. Phys (2014), 073J01 [13] N. Sato, Z. Yoshida, and Y. Kawazura, Self Organization and Heating by Inward Diffusion in Magnetospheric Plasmas, Plasma Fus. Res (2016). 10/31 [14] N. Sato and Z. Yoshida, Up-Hill Diffusion, Creation of Density Gradients: Entropy Measure for Systems with Topological Constraints, Phys. Rev. E (2016).
11 INWARD DIFFUSION: A NON-CANONICAL HAMILTONIAN SYSTEM The constraint causing non-canonicality is the conservation of μ: the pair θ c, μ can be reduced: dx dy dz dp x dp y dp z = m 2 dl dv dθ dψ dθ c dμ. Dynamical variables are reduced to four: l, v, θ, ψ. The reduced equations are the guiding center equations: ቊ ሶ v = m 1 μb + eφ l + v v E B k, v = v + v E B + v B + v k. Transforming to the new coordinates: ሶ l = v qφ θ, v ሶ = m 1 μb + eφ l + v q l φ θ, θ ሶ = ψ + q l e 1 μb + φ e 1 mv 2 q l, ψ ሶ = φ θ. 11/31
12 CONSTRUCTION OF CANONICAL COORDINATES The reduced system should be Hamiltonian because the reduced variable μ is a dynamical constant and the phase θ c can be neglected. We look for the nondegenerate closed two-form ω = dλ such that: i X ω = dh 0, with H 0 = μb + m v eφ the Hamiltonian function. Setting q = l ψ : ω = mdv dl + dψ d eθ + qv. Setting: η = eθ + qv, the canonical pairs are l, mv and η, ψ. They describe the motion parallel and perpendicular to the magnetic field respectively. 12/31
13 INVARIANT MEASURE AND APPLICATION OF THE ERGODIC HYPOTHESIS Observe that: Thus, the flow X = dl dv dθ dψ = e 1 dl dv dη dψ. l, ሶ vሶ, θ, ሶ ψሶ is measure preserving in virtue of Liouville s theorem [15]. To formulate a diffusion operator on the symplectic leaf μ, consider two hypothesis: 1. Perturbations with vanishing ensemble average E = Ergodic Hypothesis [16] on the symplectic leaf: E = dφ = m e D 1/2 1/2 1/2 Γ dl + D Γ dθ + D θ Γθ dψ, with dw = Γdt a Wiener process. The E B drift velocity causing diffusion is then written as: 1/2 dx = 1 δφ rb θ dt = D rb dw. [15] J. W. Gibbs, Elementary Principles in Statistical Mechanics, chap. 1 (Scribner s Sons, New York, 1902), pp [16] G. D. Birkhoff, Proof of the Ergodic Theorem, Proc. N. A. S. 17, pp (1931). 13/31
14 FOKKER-PLANCK EQUATION ON THE SYMPLECTIC LEAF (I) dl = v + C l α D ψ + q l q + q ψ + q l log rb dt + qd 1/2 dw dv = μ B m l + γv C v dt + D 1/2 dw D 1/2 v q l dw f t = l + v v + C l α D ψ + q l q + q ψ + q l log rb f + m2 2e 2 D D dθ = μ e ψ + q l B + C θ m e v 2 q l μ B m l + γv C v ψ D f θ μ e ψ + q l B + C θ m e v 2 q l 1 2 α ψ + q l log rb + C ψ f D 1/2 dt D θ dwθ m e qd 1/2 dw dψ = D 1 2 α ψ + q l log rb + C ψ dt + D 1/2 dw 2 θ 2 q2 f m e D 2 v θ qf + D 2 l ψ qf αd l 2 2 v v q 2 2 l f D v ψ v q l f+αd v v q l + ψ q l v q 2 l f 2 f l 2 q2 f D 2 f 2 v D 2 f ψ D 2 f θ θ 2 2 ψ + q l q f D l v v qq l f f dl dψ dθ dv dμ dθ c = fb dx dy dz dv dμ dθ c 14/31
15 CHANGE OF COORDINATES TO THE LABORATORY FRAME Diffusion is controlled by D (strength of fluctuations), α (time-properties of fluctuations), and q = l ψ (non-orthogonality of tangent vectors). Approximating D = D, α = 1/2, q~0, θ = 0, we have: f t = v f l + μ B m l f v D Note: assuming a diffusion operator arbitrarily, such as 2 2 D x f, gives wrong results. The creation of the density gradient is caused by the inhomogeneous Jacobian B linking symplectic leaf to laboratory frame: dl dψ dθ dv = B dx dy dz dv. The density ρ x, y, z in the laboratory frame x, y, z is: x 2 f ψ 2 ρ = B u = B න f dv dμ dθ c, where u is the density on the symplectic leaf. If u becomes flat due to inward diffusion, ρ must become inhomogeneous. 15/31
16 CREATION OF DENSITY GRADIENT Symplectic Leaf Density u Laboratory Frame Density ρ = B u BEFORE z AFTER r 16/31
17 EMERGENCE OF TEMPERATURE ANISOTROPY Laboratory Frame Density ρ Temperature Anisotropy T /T 17/31
18 ENTROPY (I) The system satisfies Boltzmann s H-theorem on the invariant measure [14]: dv IM = dl dψ dθ dv dμ dθ c = B dx dy dz dv dμ dθ c = B dv. Therefore, the proper entropy measure is: Σ = න Ω f log f dv IM. The entropy measure in the laboratory coordinates reads: ሚS = න Ω fb log fb dv = Σ log B. Recalling that t f = div fz, the rate of change in Σ and ሚS are: dσ dt = div Z Φ Ω = σ Φ Ω, d ሚS dt = div Z Φ Ω d dt log B = σ Φ Ω d dt log B. Here σ is the entropy production rate. [14] N. Sato and Z. Yoshida, Up-Hill Diffusion, Creation of Density Gradients: Entropy Measure for Systems with Topological Constraints, Phys. Rev. E (2016). 18/31
19 ENTROPY (II) Pictures from [13]. (a) Time evolution of entropy measures Σ and ሚS during formation of a radiation belt. (b) Time evolution of entropy production σ. This result can be extended to all Poisson operators I (and to any operator endowed with an invariant measure for any choice of H 0 ): the proper entropy measure Σ is defined on dv IM, and the solution f to the FPE satisfies: lim I d log f + βdh 0 = 0 a. e. t Darboux s theorem [14] ensures that a Poisson operator has an integrable kernel I dc i = 0. In regions U M where rank I is constant: lim f = exp βh 0 γ i C i a. e. t [14] N. Sato and Z. Yoshida, Up-Hill Diffusion, Creation of Density Gradients: Entropy Measure for Systems with Topological Constraints, Phys. Rev. E (2016). [17] M. de León and M. R. Rodrigues, Methods of Differential Geometry in Analytical Mechanics, Elsevier, p (1989). 19/31
20 TOPOLOGICAL CONSTRAINTS AND INTEGRABILITY A dipole magnetic field is integrable, i.e. it does not have helicity: B = ζ B B = 0. What happens to diffusion in an arbitrarily complex magnetic field B B 0? [20] Microscopic systems arise in Canonical Hamiltonian form. μ leaf [18] ρ x, y Macroscopic systems result from the reduction of negligible degrees of freedom that impart topological constraints. In non-canonical Hamiltonian systems, selforganization is driven by integrable constraints. The reduction process may destroy the Hamiltonian form: the standard formulation of statistical mechanics does not apply due to absence of phase space. The simplest non-hamiltonian behavior arises in 3 dimensions: V = w H 0. [18] Z. Yoshida, Self-Organization by Topological Constraints: Hierarchy of Foliated Phase Space, Adv. Phys. X 1, pp (2016). 20/31
21 6D phase space: ω 2 = dp i dx i = J c 1 E B DRIFT IS NOT AN HAMILTONIAN SYSTEM Motion of a charged particle in an e-m field m dv dt Charged particle Hamiltonian: H 0 = 1 p qa 2 + qφ 2m equations of motion: i V ω 2 = dh with V = p, ሶ xሶ = q V B + E Canonical Hamiltonian form Reduction m 0 (non canonical Hamiltonian form) 3D phase space: ω 2 = da i dx i = da = B Non inertial Hamiltonian: H 0 = qφ Equations of motion: i V da = dφ with V B + E = 0 and X = x, ሶ y, ሶ zሶ w w = Reduction V 0 V = w H 0 = B B 2 E. B B B 4 0 violation of Jacobi identity. 21/31
22 Consider a Beltrami field w = μ w. TOPOLOGICAL CONSTRAINTS IN 3D Since V = w H 0, the motion is constrained to be orthogonal to w: w V = 0. Euler s equation for the motion of a rigid body with angular momentum x, moment of inertia I = I x, I y, I z and Hamiltonian H 0 = 1 2 This system has a non-integrable constraint since the Jacobi identity is not satisfied: x 2 V = x H 0. I x + y2 I y + z2 I z : This is a non-canonical Hamiltonian system with integrable constraint x V = 0 which can be integrated to the Casimir invariant C = x 2 /2: ሶ C = C V = x V = 0. w w = μw /31
23 CONSTRAINED ORBITS Euler Rigid Body (Hamiltonian) w = x Beltrami (Non-Hamiltonian) w = cos z sin y, sin z, cos y X = w H = w H 0 + Γ 23/31
24 E B DIFFUSION AND CLASSIFICATION OF CONSERVATIVE DYNAMICS Assume charge neutrality H = H 0 + Γ = e φ + Γ = Γ. Then: V = w Γ, ρ t = 1 w ρw. 2 The form of ρ x depends on the geometrical properties of w = B/B 2. Magnetic field with helicity B B 0. Magnetic field without helicity B B = 0. NON-BELTRAMI w = w ρ + b ρ + bρ = 0 BELTRAMI b = 0 ρ = constant POISSON (Hamiltonian) w w = 0 w = λ C ρ λ = f C NULL w = 0 ρ = ρ NON-BELTRAMINESS b = b = w w 24/31
25 NUMERICAL TEST OF 3D DIFFUSION (I) The theory is tested by comparing the numerical integration of the stochastic equation of motion with the analytical solution of the Fokker-Planck equation. (a) (b) (c) Sample size: particles per simulation. (a) Initial uniform particle distribution in 3D cube with periodic boundaries. (b) Noise Γ i as a function of time t. (c) 250 sample paths in x, y plane for w = 1 + cos 2 x z cos x cos y. 25/31
26 w B/B 2 NUMERICAL TEST OF 3D DIFFUSION (II) w 1 x, y B x, y ρ x, y Poisson Operator (Constant Magnetic Field) w = B = z. Analytical Solution ρ = constant. Poisson Operator (Gradient Magnetic Field) w = z cos x cos y. Analytical Solution ρ = constant. 26/31
27 w w = B/B 2 B/B 2 NUMERICAL TEST OF 3D DIFFUSION (III) w 1 x, y B x, y B x, y ρ x, y Poisson Operator (Gradient Magnetic Field) w = λ C = 1 + cos 2 x z cos x cos y. Analytical Solution ρ = λ 1 f C. Beltrami Operator (Beltrami Magnetic Field) w = 2B = cos z + sin z cos z sin z 0 Analytical Solution ρ = constant. 27/31
28 w w = B/B 2 B/B 2 NUMERICAL TEST OF 3D DIFFUSION (IV) w 1 x, y B x, y B x, y ρ x, y Non-Beltrami Operator (Non-Integrable Magnetic Field) b 0 w = 1, sin x + cos y, cos x. Non-Beltrami Operator (Non-Integrable Magnetic Field) (free boundaries) b 0 w = 1 + y sin y cos y 1,0, 2 y sin y cos y 2 sin x sin x 2. 28/31
29 NUMERICAL TEST OF 3D DIFFUSION (V) In the most general case of a non-integrable operator w w 0 with b 0 the density profile is determined by field strength w and non-beltraminess b. To elucidate the role of b, we consider an operator of uniform strength: w 1 = constant. w w = B/B 2 B/B 2 b = b b x, y = w w ρ x, y Non-Beltrami Operator (Non-Integrable Magnetic Field) b 0 w = cos y, cos x, sin y 1 + cos 2 x w = constant 29/31
30 CONCLUSION Dipole fields can efficiently trap charged particles by the mechanism of inward diffusion. The inward diffusion occurs in virtue of the topological constraint imposed by μ. The creation of the density gradient is consistent with the second law of thermodynamics. The density profile resulting from diffusion by a white noise electric field in an arbitrary magnetic field B is determined by its geometric properties. POISSON (Hamiltonian) BELTRAMI w w = 0 w = λ C b = 0 ρ = ρ λ, C ρ = constant In 3D, three classes of w = B/B 2 can be identified: Poisson, Beltrami, and Non-Beltrami. NON- BELTRAMI b 0 ρ = ρ w, b The construction of a diffusion operator requires careful treatment, especially when the system is not endowed with a canonical Hamiltonian structure. The statistical theory developed here can be applied to general conservative systems. Diffusion by E B drift is of special mathematical interest because it involves a non-elliptic partial differential equation. 30/31
31 End
What is macro? self-organized structures in physical systems
What is macro? self-organized structures in physical systems Z. Yoshida U. Tokyo 2016.5.24 Z. Yoshida (U. Tokyo) Self-organization in macro systems 2016/05/24 1 / 15 Outline Self-organization without blueprints
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 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 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 informationarxiv: v1 [math-ph] 23 Aug 2018
arxiv:1808.07706v1 [math-ph] 23 Aug 2018 Beltrami Operators N. Sato 1 1 Research Institute for Mathematical Sciences, Kyoto University, Kyoto 606-8502, Japan Email: sato@kurims.kyoto-u.ac.jp August 24,
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 informationConfinement of toroidal non-neutral plasma in Proto-RT
Workshop on Physics with Ultra Slow Antiproton Beams, RIKEN, March 15, 2005 Confinement of toroidal non-neutral plasma in Proto-RT H. Saitoh, Z. Yoshida, and S. Watanabe Graduate School of Frontier Sciences,
More informationConfinement of toroidal non-neutral plasma in Proto-RT
Workshop on Physics with Ultra Slow Antiproton Beams, RIKEN, March 15, 2005 Confinement of toroidal non-neutral plasma in Proto-RT H. Saitoh, Z. Yoshida, and S. Watanabe Graduate School of Frontier Sciences,
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 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 informationHamiltonian Dynamics
Hamiltonian Dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS Feb. 10, 2009 Joris Vankerschaver (CDS) Hamiltonian Dynamics Feb. 10, 2009 1 / 31 Outline 1. Introductory concepts; 2. Poisson brackets;
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 informationFormation of High-b ECH Plasma and Inward Particle Diffusion in RT-1
J Fusion Energ (2010) 29:553 557 DOI 10.1007/s10894-010-9327-6 ORIGINAL RESEARCH Formation of High-b ECH Plasma and Inward Particle Diffusion in RT-1 H. Saitoh Z. Yoshida J. Morikawa Y. Yano T. Mizushima
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamics in High Energy Accelerators Part 2: Basic tools and concepts Nonlinear Single-Particle Dynamics in High Energy Accelerators This course consists of eight lectures: 1.
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 informationSelf-Organization of Plasmas with Flows
Self-Organization of Plasmas with Flows ICNSP 2003/ 9/10 Graduate School of Frontier Sciences,, National Institute for Fusion Science R. NUMATA, Z. YOSHIDA, T. HAYASHI ICNSP 2003/ 9/10 p.1/14 Abstract
More informationGEOMETRIC QUANTIZATION
GEOMETRIC QUANTIZATION 1. The basic idea The setting of the Hamiltonian version of classical (Newtonian) mechanics is the phase space (position and momentum), which is a symplectic manifold. The typical
More informationCaltech Ph106 Fall 2001
Caltech h106 Fall 2001 ath for physicists: differential forms Disclaimer: this is a first draft, so a few signs might be off. 1 Basic properties Differential forms come up in various parts of theoretical
More informationCANONICAL EQUATIONS. Application to the study of the equilibrium of flexible filaments and brachistochrone curves. By A.
Équations canoniques. Application a la recherche de l équilibre des fils flexibles et des courbes brachystochrones, Mem. Acad. Sci de Toulouse (8) 7 (885), 545-570. CANONICAL EQUATIONS Application to the
More informationList of Comprehensive Exams Topics
List of Comprehensive Exams Topics Mechanics 1. Basic Mechanics Newton s laws and conservation laws, the virial theorem 2. The Lagrangian and Hamiltonian Formalism The Lagrange formalism and the principle
More informationProspects for Driven Particle Convection Tests in LDX. Abstract
Prospects for Driven Particle Convection Tests in LDX M.E. Mauel, A.C. Boxer, J.L. Ellsworth, D.T. Garnier, J. Kesner Columbia University and Plasma Science and Fusion Center, MIT 49th Meeting of the APS
More informationComplete integrability of geodesic motion in Sasaki-Einstein toric spaces
Complete integrability of geodesic motion in Sasaki-Einstein toric spaces Mihai Visinescu Department of Theoretical Physics National Institute for Physics and Nuclear Engineering Horia Hulubei Bucharest,
More informationControl of chaos in Hamiltonian systems
Control of chaos in Hamiltonian systems G. Ciraolo, C. Chandre, R. Lima, M. Vittot Centre de Physique Théorique CNRS, Marseille M. Pettini Osservatorio Astrofisico di Arcetri, Università di Firenze Ph.
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 informationWhat we ve learned so far about the Stability of Plasma Confined by a Laboratory Dipole Magnet
What we ve learned so far about the Stability of Plasma Confined by a Laboratory Dipole Magnet M. E. Mauel and the CTX and LDX Experimental Teams Annual Meeting of the Division of Plasma Physics Philadelphia,
More informationThe Virial Theorem, MHD Equilibria, and Force-Free Fields
The Virial Theorem, MHD Equilibria, and Force-Free Fields Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 10 12, 2014 These lecture notes are largely
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 information(a) Sections 3.7 through (b) Sections 8.1 through 8.3. and Sections 8.5 through 8.6. Problem Set 4 due Monday, 10/14/02
Physics 601 Dr. Dragt Fall 2002 Reading Assignment #4: 1. Dragt (a) Sections 1.5 and 1.6 of Chapter 1 (Introductory Concepts). (b) Notes VI, Hamilton s Equations of Motion (to be found right after the
More informationMAGNETOHYDRODYNAMICS
Chapter 6 MAGNETOHYDRODYNAMICS 6.1 Introduction Magnetohydrodynamics is a branch of plasma physics dealing with dc or low frequency effects in fully ionized magnetized plasma. In this chapter we will study
More informationanalysis for transport equations and applications
Multi-scale analysis for transport equations and applications Mihaï BOSTAN, Aurélie FINOT University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Numerical methods for kinetic equations Strasbourg
More informationThe linear equations can be classificed into the following cases, from easier to more difficult: 1. Linear: u y. u x
Week 02 : Method of C haracteristics From now on we will stu one by one classical techniques of obtaining solution formulas for PDEs. The first one is the method of characteristics, which is particularly
More informationLecture 16 March 29, 2010
Lecture 16 March 29, 2010 We know Maxwell s equations the Lorentz force. Why more theory? Newton = = Hamiltonian = Quantum Mechanics Elegance! Beauty! Gauge Fields = Non-Abelian Gauge Theory = Stard Model
More informationMassachusetts Institute of Technology Department of Physics. Final Examination December 17, 2004
Massachusetts Institute of Technology Department of Physics Course: 8.09 Classical Mechanics Term: Fall 004 Final Examination December 17, 004 Instructions Do not start until you are told to do so. Solve
More informationSolutions of M3-4A16 Assessed Problems # 3 [#1] Exercises in exterior calculus operations
D. D. Holm Solutions to M3-4A16 Assessed Problems # 3 15 Dec 2010 1 Solutions of M3-4A16 Assessed Problems # 3 [#1] Exercises in exterior calculus operations Vector notation for differential basis elements:
More informationGeometric Gyrokinetic Theory and its Applications to Large-Scale Simulations of Magnetized Plasmas
Geometric Gyrokinetic Theory and its Applications to Large-Scale Simulations of Magnetized Plasmas Hong Qin Princeton Plasma Physics Laboratory, Princeton University CEA-EDF-INRIA School -- Numerical models
More informationThe Levitated Dipole Experiment: Towards Fusion Without Tritium
The Levitated Dipole Experiment: Towards Fusion Without Tritium Jay Kesner MIT M.S. Davis, J.E. Ellsworth, D.T. Garnier, M.E. Mauel, P.C. Michael, P.P. Woskov MCP I3.110 Presented at the EPS Meeting, Dublin,
More informationHelicity fluctuation, generation of linking number and effect on resistivity
Helicity fluctuation, generation of linking number and effect on resistivity F. Spineanu 1), M. Vlad 1) 1) Association EURATOM-MEC Romania NILPRP MG-36, Magurele, Bucharest, Romania spineanu@ifin.nipne.ro
More informationRadiation belt particle dynamics
Radiation belt particle dynamics Prepared by Kevin Graf Stanford University, Stanford, CA IHY Workshop on Advancing VLF through the Global AWESOME Network Basic Motion Motion of charged particle q in presence
More information(a) Write down the total Hamiltonian of this system, including the spin degree of freedom of the electron, but neglecting spin-orbit interactions.
1. Quantum Mechanics (Spring 2007) Consider a hydrogen atom in a weak uniform magnetic field B = Bê z. (a) Write down the total Hamiltonian of this system, including the spin degree of freedom of the electron,
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 informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
More informationThe Equipartition Theorem
Chapter 8 The Equipartition Theorem Topics Equipartition and kinetic energy. The one-dimensional harmonic oscillator. Degrees of freedom and the equipartition theorem. Rotating particles in thermal equilibrium.
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 informationSingle Particle Motion
Single Particle Motion Overview Electromagnetic fields, Lorentz-force, gyration and guiding center, drifts, adiabatic invariants. Pre-requisites: Energy density of the particle population smaller than
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 informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 1 Defining statistical mechanics: Statistical Mechanics provies the connection between microscopic motion of individual atoms of matter and macroscopically
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 informationIntroduction. Statistical physics: microscopic foundation of thermodynamics degrees of freedom 2 3 state variables!
Introduction Thermodynamics: phenomenological description of equilibrium bulk properties of matter in terms of only a few state variables and thermodynamical laws. Statistical physics: microscopic foundation
More informationSingle particle motion and trapped particles
Single particle motion and trapped particles Gyromotion of ions and electrons Drifts in electric fields Inhomogeneous magnetic fields Magnetic and general drift motions Trapped magnetospheric particles
More informationPhysics 607 Exam 2. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2
Physics 607 Exam Please be well-organized, and show all significant steps clearly in all problems. You are graded on your work, so please do not just write down answers with no explanation! Do all your
More informationConservation Laws in Ideal MHD
Conservation Laws in Ideal MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 3, 2016 These lecture notes are largely based on Plasma Physics for Astrophysics
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 1 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationTHE LONGITUDINAL KAM-COCYCLE OF A MAGNETIC FLOW
THE LONGITUDINAL KAM-COCYCLE OF A MAGNETIC FLOW GABRIEL P. PATERNAIN Abstract. Let M be a closed oriented surface of negative Gaussian curvature and let Ω be a non-exact 2-form. Let λ be a small positive
More informationStochastic equations for thermodynamics
J. Chem. Soc., Faraday Trans. 93 (1997) 1751-1753 [arxiv 1503.09171] Stochastic equations for thermodynamics Roumen Tsekov Department of Physical Chemistry, University of Sofia, 1164 Sofia, ulgaria The
More informationMULTIVARIABLE INTEGRATION
MULTIVARIABLE INTEGRATION (SPHERICAL POLAR COORDINATES) Question 1 a) Determine with the aid of a diagram an expression for the volume element in r, θ, ϕ. spherical polar coordinates, ( ) [You may not
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationFundamentals of Magnetic Island Theory in Tokamaks
Fundamentals of Magnetic Island Theory in Tokamaks Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Austin, TX, USA Talk available at http://farside.ph.utexas.edu/talks/talks.html
More informationExperiments with a Supported Dipole
Experiments with a Supported Dipole Reporting Measurements of the Interchange Instability Excited by Electron Pressure and Centrifugal Force Introduction Ben Levitt and Dmitry Maslovsky Collisionless Terrella
More informationClassical Mechanics and Statistical/Thermodynamics. January 2015
Classical Mechanics and Statistical/Thermodynamics January 2015 1 Handy Integrals: Possibly Useful Information 0 x n e αx dx = n! α n+1 π α 0 0 e αx2 dx = 1 2 x e αx2 dx = 1 2α 0 x 2 e αx2 dx = 1 4 π α
More informationEva Miranda. UPC-Barcelona and BGSMath. XXV International Fall Workshop on Geometry and Physics Madrid
b-symplectic manifolds: going to infinity and coming back Eva Miranda UPC-Barcelona and BGSMath XXV International Fall Workshop on Geometry and Physics Madrid Eva Miranda (UPC) b-symplectic manifolds Semptember,
More informationcos 6 λ m sin 2 λ m Mirror Point latitude Equatorial Pitch Angle Figure 5.1: Mirror point latitude as function of equatorial pitch angle.
Chapter 5 The Inner Magnetosphere 5.1 Trapped Particles The motion of trapped particles in the inner magnetosphere is a combination of gyro motion, bounce motion, and gradient and curvature drifts. In
More informationFormation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )
Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University
More information[variable] = units (or dimension) of variable.
Dimensional Analysis Zoe Wyatt wyatt.zoe@gmail.com with help from Emanuel Malek Understanding units usually makes physics much easier to understand. It also gives a good method of checking if an answer
More informationStatus of A Positron-Electron Experiment (APEX) towards the formation of pair plasmas
Status of A Positron-Electron Experiment (APEX) towards the formation of pair plasmas H. Saitoh 1,5, J. Stanja 1, T. Sunn Pedersen 1,3, U. Hergenhahn 1, E. V. Stenson 1, H. Niemann 1,3, N. Paschkowski
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 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 informationTowards Multi-field Inflation with a Random Potential
Towards Multi-field Inflation with a Random Potential Jiajun Xu LEPP, Cornell Univeristy Based on H. Tye, JX, Y. Zhang, arxiv:0812.1944 and work in progress 1 Outline Motivation from string theory A scenario
More informationPHYS 352 Homework 2 Solutions
PHYS 352 Homework 2 Solutions Aaron Mowitz (, 2, and 3) and Nachi Stern (4 and 5) Problem The purpose of doing a Legendre transform is to change a function of one or more variables into a function of variables
More informationIV. Classical Statistical Mechanics
IV. Classical Statistical Mechanics IV.A General Definitions Statistical Mechanics is a probabilistic approach to equilibrium macroscopic properties of large numbers of degrees of freedom. As discussed
More informationChapter 2 Radiation of an Accelerated Charge
Chapter 2 Radiation of an Accelerated Charge Whatever the energy source and whatever the object, (but with the notable exception of neutrino emission that we will not consider further, and that of gravitational
More informationMotion of a charged particle in an EM field
Department of physics Seminar - 4 th year Motion of a charged particle in an EM field Author: Lojze Gačnik Mentor: Assoc. Prof. Tomaž Gyergyek Ljubljana, November 2011 Abstract In this paper I will present
More informationMultiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015
Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction
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 informationStability Subject to Dynamical Accessibility
Stability Subject to Dynamical Accessibility P. J. Morrison Department of Physics and Institute for Fusion Studies The University of Texas at Austin morrison@physics.utexas.edu http://www.ph.utexas.edu/
More informationDynamics of charged particles in spatially chaotic magnetic fields
PSFC/JA-1-38 Dynamics of charged particles in spatially chaotic magnetic fields Abhay K. Ram and Brahmananda Dasgupta a October 21 Plasma Science and Fusion Center, Massachusetts Institute of Technology
More information### dv where ρ is density, R is distance from rotation axis, and dv is
Comments This is one of my favorite problem assignments in our Solid Earth Geophysics class, typically taken by junior and senior concentrators and by first-year graduate students. I encourage students
More informationFinal Review Prof. WAN, Xin
General Physics I Final Review Prof. WAN, Xin xinwan@zju.edu.cn http://zimp.zju.edu.cn/~xinwan/ About the Final Exam Total 6 questions. 40% mechanics, 30% wave and relativity, 30% thermal physics. Pick
More informationAn introduction to General Relativity and the positive mass theorem
An introduction to General Relativity and the positive mass theorem National Center for Theoretical Sciences, Mathematics Division March 2 nd, 2007 Wen-ling Huang Department of Mathematics University of
More informationDynamic and Thermodynamic Stability of Black Holes and Black Branes
Dynamic and Thermodynamic Stability of Black Holes and Black Branes Robert M. Wald with Stefan Hollands arxiv:1201.0463 Commun. Math. Phys. 321, 629 (2013) (see also K. Prabhu and R.M. Wald, Commun. Math.
More informationRadiative Processes in Astrophysics
Radiative Processes in Astrophysics 6. Relativistic Covariance & Kinematics Eline Tolstoy http://www.astro.rug.nl/~etolstoy/astroa07/ Practise, practise, practise... mid-term, 31st may, 9.15-11am As we
More informationMicroscopic electrodynamics. Trond Saue (LCPQ, Toulouse) Microscopic electrodynamics Virginia Tech / 46
Microscopic electrodynamics Trond Saue (LCPQ, Toulouse) Microscopic electrodynamics Virginia Tech 2015 1 / 46 Maxwell s equations for electric field E and magnetic field B in terms of sources ρ and j The
More informationSolution. For one question the mean grade is ḡ 1 = 10p = 8 and the standard deviation is 1 = g
Exam 23 -. To see how much of the exam grade spread is pure luck, consider the exam with ten identically difficult problems of ten points each. The exam is multiple-choice that is the correct choice of
More informationLow Emittance Machines
CERN Accelerator School Advanced Accelerator Physics Course Trondheim, Norway, August 2013 Low Emittance Machines Part 1: Beam Dynamics with Synchrotron Radiation Andy Wolski The Cockcroft Institute, and
More informationBACKGROUND IN SYMPLECTIC GEOMETRY
BACKGROUND IN SYMPLECTIC GEOMETRY NILAY KUMAR Today I want to introduce some of the symplectic structure underlying classical mechanics. The key idea is actually quite old and in its various formulations
More information[#1] R 3 bracket for the spherical pendulum
.. Holm Tuesday 11 January 2011 Solutions to MSc Enhanced Coursework for MA16 1 M3/4A16 MSc Enhanced Coursework arryl Holm Solutions Tuesday 11 January 2011 [#1] R 3 bracket for the spherical pendulum
More informationFundamental Concepts of Particle Accelerators III : High-Energy Beam Dynamics (2) Koji TAKATA KEK. Accelerator Course, Sokendai. Second Term, JFY2012
.... Fundamental Concepts of Particle Accelerators III : High-Energy Beam Dynamics (2) Koji TAKATA KEK koji.takata@kek.jp http://research.kek.jp/people/takata/home.html Accelerator Course, Sokendai Second
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 informationCurves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,
Physics 6010, Fall 2010 Hamiltonian Formalism: Hamilton s equations. Conservation laws. Reduction. Poisson Brackets. Relevant Sections in Text: 8.1 8.3, 9.5 The Hamiltonian Formalism We now return to formal
More informationApplications to simulations: Monte-Carlo
Applications to simulations: Monte-Carlo A.C. Maggs ESPCI, Paris June 2016 Summary Algorithms Faster/simpler codes Thermodynamics of Electric fields Partition function of electric field Thermal Casimir/Lifshitz
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 informationDivergence Theorem December 2013
Divergence Theorem 17.3 11 December 2013 Fundamental Theorem, Four Ways. b F (x) dx = F (b) F (a) a [a, b] F (x) on boundary of If C path from P to Q, ( φ) ds = φ(q) φ(p) C φ on boundary of C Green s Theorem:
More informationJP Sta,onary Density Profiles in Alcator C Mod
JP8.00072 Sta,onary Density Profiles in Alcator C Mod 1 In the absence of an internal particle source, plasma turbulence will impose an intrinsic relationship between an inwards pinch and an outwards diffusion
More informationThe Linear Theory of Tearing Modes in periodic, cyindrical plasmas. Cary Forest University of Wisconsin
The Linear Theory of Tearing Modes in periodic, cyindrical plasmas Cary Forest University of Wisconsin 1 Resistive MHD E + v B = ηj (no energy principle) Role of resistivity No frozen flux, B can tear
More informationProblem 1, Lorentz transformations of electric and magnetic
Problem 1, Lorentz transformations of electric and magnetic fields We have that where, F µν = F µ ν = L µ µ Lν ν F µν, 0 B 3 B 2 ie 1 B 3 0 B 1 ie 2 B 2 B 1 0 ie 3 ie 2 ie 2 ie 3 0. Note that we use the
More informationTopics for the Qualifying Examination
Topics for the Qualifying Examination Quantum Mechanics I and II 1. Quantum kinematics and dynamics 1.1 Postulates of Quantum Mechanics. 1.2 Configuration space vs. Hilbert space, wave function vs. state
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 informationELECTROHYDRODYNAMICS IN DUSTY AND DIRTY PLASMAS
ELECTROHYDRODYNAMICS IN DUSTY AND DIRTY PLASMAS Gravito-Electrodynamics and EHD by HIROSHI KIKUCHI Institute for Environmental Electromagnetics, Tokyo, Japan KLUWER ACADEMIC PUBLISHERS DORDRECHT / BOSTON
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 informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More information3.320 Lecture 18 (4/12/05)
3.320 Lecture 18 (4/12/05) Monte Carlo Simulation II and free energies Figure by MIT OCW. General Statistical Mechanics References D. Chandler, Introduction to Modern Statistical Mechanics D.A. McQuarrie,
More information