Introduction to MagnetoHydroDynamics (MHD) Antoine Cerfon, Courant Institute, New York University
|
|
- Hope Jennings
- 5 years ago
- Views:
Transcription
1 Introduction to MagnetoHydroDynamics (MHD) Antoine Cerfon, Courant Institute, New York University SULI Introductory Course in Plasma Physics, June 6, 2016
2 PART I: DESCRIBING A FUSION PLASMA
3 METHOD I: SELF-CONSISTENT PARTICLE PUSHING Intuitive idea: solve for the motion of all particles iteratively, combining Newton s law with Maxwell equations At each time step i, solve m d2 x (i) k dt ( E ( ) = q k E(x k ) (i 1) + dx(i) k B (i 1) (x dt k ) ) i N = c 2 B (i 1) µ 0 c 2 dx i k q k dt δ(x xi k ) t ( ) B i = E i t k=1 k = 1,..., N Even with best available numerical solvers, problem still not tractable even with the most powerful computers when N as in magnetic fusion grade plasmas
4 METHOD II: COARSE-GRAIN AVERAGE IN PHASE SPACE (From G. Lapenta s: u /weather/pic.pdf ) For hot and diffuse systems with large # of particles, following every single particle is a waste of time and resources Replace the discrete particles with smooth distribution function f (x, v, t) defined so that f (x, v, t)dxdv is the # of particles in the infinitesimal six-dimensional phase-space volume dxdv.
5 DISTRIBUTION FUNCTION AND VLASOV EQUATION Macroscopic (fluid) quantities are velocity moments of f n(x, t) = f (x, v, t)dv Density nv(x, t) = vf (x, v, t)dv Mean flow P(x, t) = m (v V) (v V) fdv Pressure tensor Conservation of f along the phase-space trajectories of the particles determines the time evolution of f : df dt = f t + dx dv f + dt dt vf = 0 dx dt = v dv dt = q (E + v B) m f t + v f + q m (E + v B) vf = 0 This is the Vlasov equation
6 THE BOLTZMANN EQUATION In fusion plasmas, we separate, leading to the Boltzmann equation: f t + v f + q ( ) f m (E + v B) vf = t c This equation to be combined with Maxwell s equations: E = B t B = µ 0 J + 1 c 2 E t Nonlinear, integro-differential, 6-dimensional PDE Describes phenomena on widely varying length ( m) and time ( s) scales Still not a piece of cake, and never solved as such for fusion plasmas
7 MOMENT APPROACH f t + v f + q ( ) f m (E + v B) vf = t c Taking the integrals dv, mvdv and mv 2 /2dv of this equation, we obtain the exact fluid equations: n s t + (n sv s ) = 0 Continuity ( ) Vs mn + V s V s = q s n s (E + V s B) P s + R s Momen t ( ) d 3 dt 2 p s p s V s + π s : V s + q s = 0 (Energy) with P s = p s I + π s. Closure problem: for each moment, we introduce a new unknown End up with too many unknowns Need to make approximations to close the moment hierarchy
8 KINETIC MODELS VS FLUID MODELS For some fusion applications/plasma regimes (heating and current drive, transport), kinetic treatment cannot be avoided Simplify and reduce dimensionality of the Vlasov equation with approximations: Strong magnetization : Gyrokinetic equation Small gyroradius compared to relevant length scales : Drift kinetic equation Vanishing gyroradius : Kinetic MHD In contrast, fluid models are based on approximate expressions for higher order moments (off-diagonal entries in pressure tensor, heat flux) in terms of lower order quantities(density, velocity, diagonal entries in pressure tensor) We will now focus on the relevant regime and the approximations made to derive a widely used fluid model: the ideal MHD model
9 PART II: THE IDEAL MHD MODEL
10 The maximum p is limited by the stability properties Job of MHD The maximum τ E is determined by the confinement properties Job of kinetic models LAWSON CRITERION AND MHD Condition for ignition: pτ E 8 bar.s T min 15keV
11 PHILOSOPHY The purpose of ideal MHD is to study the macroscopic behavior of the plasma Use ideal MHD to design machines that avoid large scale instabilities Regime of interest Typical length scale: the minor radius of the device a 1m Wave number k of waves and instabilitities considered: k 1/a Typical velocities: Ion thermal velocity speed v T 500km/s Typical time scale: τ MHD a/v T 2µs Frequency ω MHD of associated waves/instabilities ω MHD 500kHz
12 EXAMPLE: VERTICAL INSTABILITY Figure from F. Hofmann et al., Nuclear Fusion (1997)
13 IDEAL MHD - MAXWELL S EQUATIONS a λ D, the distance over which charge separation can take place in a plasma On the MHD length scale, the plasma is neutral : n i = n e ω MHD /k c and v Ti v Te c so we can neglect the displacement current in Maxwell s equations: n i = n e B = 0 E = B t B = µ 0 J
14 IDEAL MHD - MOMENTUM EQUATION a λ D and a r Le (electron Larmor radius) ω MHD ω pe, ω MHD ω ce The ideal MHD model assumes that on the time and length scales of interest, the electrons have an infinitely fast response time to changes in the plasma Mathematically, this can be done by taking the limit m e 0 Adding the ion and electron momentum equation, we then get ρ dv dt J B + p = (π i + π e ) where ρ = m i n and V is the ion fluid velocity If the condition v Ti τ ii /a 1 is satisfied in the plasma ρ dv dt = J B p (Ideal MHD momentum equation)
15 IDEAL MHD - ELECTRONS In the limit m e 0, the electron momentum equation can be written as E + V B = 1 en (J B p e π e + R e ) This is called the generalized Ohm s law Different MHD models (resitive MHD, Hall MHD) keep different terms in this equation If r Li /a 1, v Ti τ ii /a 1, and (m e /m i ) 1/2 (r Li /a) 2 (a/v Ti τ ii ) 1, the momentum equation becomes the ideal Ohm s law E + V B = 0 The ideal MHD plasma behaves like a perfectly conducting fluid
16 ENERGY EQUATION Define the total plasma pressure p = p i + p e Add electron and ion energy equations Under the conditions r Li /a 1 and v Ti τ ii /a 1, this simplifies as ( ) d p dt ρ 5/3 = 0 Equation reminiscent of pv γ = Cst: the ideal MHD plasma behaves like a monoatomic ideal gas undergoing a reversible adiabatic process
17 IDEAL MHD - SUMMARY Valid under the conditions ( ) 1/2 (vi mi τ ) ii 1 a m e ρ t + (ρv) = 0 ρ dv dt = J B p ( ) d p dt ρ 5/3 = 0 E + V B = 0 E = B t B = µ 0 J B = 0 r Li a 1 ( rli a ) 2 ( me m i ) 1/2 a v Ti τ ii 1
18 VALIDITY OF THE IDEAL MHD MODEL (I) Are the conditions for the validity of ideal MHD ( ) 1/2 (vi mi τ ) ii r ( Li 1 a a 1 rli a m e mutually compatible? Define x = (m i /m e ) 1/2 (v Ti τ ii /a), y = r Li /a. ) 2 ( me m i ) 1/2 a v Ti τ ii 1 x 1 (High collisionality) y 1 (Small ion Larmor radius) y 2 /x 1 (Small resistivity) There exists a regime for which ideal MHD is justified (Figure from Ideal MHD by J.P. Freidberg, CUP, 2014) Is that the regime of magnetic confinement fusion?
19 VALIDITY OF THE IDEAL MHD MODEL (II) Express three conditions in terms of usual physical parameters: n, T, a For tokamak-like pressures and a = 1m, we find: The regime of validity of ideal MHD does NOT coincide with the fusion plasma regime (Figure from Ideal MHD by J.P. Freidberg, CUP, 2014) The collisionality of fusion plasmas is too low for the ideal MHD model to be valid. Is that a problem?
20 VALIDITY OF THE IDEAL MHD MODEL (III) It turns out that ideal MHD often does a very good job at predicting stability limits for macroscopic instabilities This is not due to luck but to subtle physical reasons One can show that collisionless kinetic models for macroscopic instabilities are more optimistic than ideal MHD This is because ideal MHD is accurate for dynamics perpendicular to the fields lines Designs based on ideal MHD calculations are conservative designs
21 FROZEN IN LAW E + V B = 0: in the frame moving with the plasma, the electric field is zero The plasma behaves like a perfect conductor The magnetic field lines are frozen into the plasma motion
22 MAGNETIC RECONNECTION Image from Principles of Magnetohydrodynamics With Applications to Laboratory and Astrophysical Plasmas by J.P. Goedbloed and S. Poedts, Cambridge University Press (2004) Magnetic reconnection: a key phenomenon in astrophysical, space, and fusion plasmas Cannot happen according to ideal MHD Need to add additional terms in Ohm s law to allow reconnection: resistivity, off-diagonal pressure tensor terms, electron inertia,... E + V B = ηj + 1 ne J B 1 ne P e + m e J ne 2 t Associated instabilities take place on longer time scales than τ MHD
23
24 PART III: MHD EQUILIBRIUM
25 EQUILIBRIUM STATE By equilibrium, we mean steady-state: / t = 0 Often, for simplicity and/or physical reasons, we focus on static equilibria: V = 0 A more condensed form is B = 0 B = µ 0 J J B = p B = 0 ( B) B = µ 0 p Note that the density profile does not appear
26 1D EQUILIBRIA (I) θ-pinch Z-pinch Combine the two to get... (Figure from Ideal MHD by J.P. Freidberg, CUP, 2014)
27 1D EQUILIBRIA (II) Screw pinch (Figure from Ideal MHD by J.P. Freidberg, CUP, 2014) Equilibrium quantities only depend on r Plug into B = 0, ( B) B = µ 0 p to find: ( d p + B2 θ + ) B2 z + B2 θ dr 2µ 0 µ 0 r = 0 Balance between plasma pressure, magnetic pressure, and magnetic tension Two free functions define equilibrium: e.g. B z and p, or B θ and B z
28 GENERAL EQUILIBRIA Equilibrium relation J B = p (J B = p) B B p = 0 Magnetic field is tangent to surfaces of constant pressure (J B = p) J = 0 J p = 0 Current density is tangent to surfaces of constant pressure
29 (Figure from Plasma Physics and Fusion Energy by J.P. Freidberg, CUP, 2008)
30 The regions of constant pressure are nested toroidal surfaces Magnetic fields and currents lie on these nested surfaces
31 GRAD-SHAFRANOV EQUATION Surfaces of constant pressure coincide with surfaces of constant magnetic flux Ψ These are given by the Grad-Shafranov equation R ( ) 1 Ψ + 2 Ψ R R R Z 2 = µ 0R 2 dp dψ F df dψ Second-order, nonlinear, elliptic PDE. Derived independently by H. Grad 1 and V.D. Shafranov 2. The free functions p and F determine the nature of the equilibrium In general, the GSE has to be solved numerically 1 Proceedings of the Second United Nations Conference on the Peaceful Uses of Atomic Energy, Vol. 31, p Sov. Phys. JETP 6, 545 (1958)
32 EXAMPLES (I) Z/R R/R 0 Grad-Shafranov equilibrium for JET tokamak
33 EXAMPLES (II) 5 Y X Grad-Shafranov equilibrium for Field Reversed Configuration
34 3D EQUILIBRIA / φ 0
35 3D EQUILIBRIA Equilibrium equations B = 0, ( B) B = µ 0 p still hold Existence of nested toroidal surfaces not guaranteed anymore Tokamak Stellarator (3D) Z/R R/R 0 Computing 3D equilibria fast and accurately still a challenge Several existing codes, based on different assumptions/approximations and used to design and study stellarators: VMEC, PIES, SPEC, HINT, NSTAB
36 PART IV: MHD STABILITY
37 WHAT DO WE MEAN BY MHD STABILITY? That the plasma is initially in equilibrium does not mean it is going to remain there The plasma is constantly subject to perturbations, small and large The purpose of stability studies is to find out how the plasma will react to these perturbations Will it try to return to the initial steady-state? Will it find a new acceptable steady-state? Will it collapse?
38 A MECHANICAL ANALOG Figure from J.P. Freidberg, Ideal MHD, Cambridge University Press (2014)
39 SOLVING FULL NONLINEAR MHD EQUATIONS Here is an idea to study stability of a magnetically confined plasma: Choose a satisfying plasma equilibrium Perturb it Solve the full MHD equations with a computer Analyze results Such an approach provides knowledge of the entire plasma dynamics There exist several numerical codes that can do that, for various MHD models (not only ideal): M3D, M3D-C1, NIMROD Computationally intensive Get more information than one needs?
40 LINEAR STABILITY (I) Ideal MHD dynamics can be so fast and detrimental that one may often require linear stability for the equilibrium This can simplify the mathematical analysis tremendously Start with an MHD equilibrium: B 0 = 0, ( B 0 ) B 0 = µ 0 p 0 Take full ideal MHD equations, and write Q = Q 0 (r) + Q 1 (r, t) for each physical quantity, where Q 1 is considered very small compared to Q 0 Drop all the terms that are quadratic or higher orders in the quantities Q 1 (linearization)
41 LINEAR STABILITY (II) ρ 1 t + (ρ 0v 1 ) = 0 v 1 ρ 0 = J 1 B 0 + J 0 B 1 p 1 t p 1 t + v 1 p + γp v 1 = 0 B 1 = (v 1 B 0 ) t B 1 = 0 µ 0 J 1 = B 1 By design, the system is now linear in the unknown quantities ρ 1, v 1, J 1, B 1, p 1 Much easier to solve in a computer or by hand!
42 ILLUSTRATION: WAVES IN IDEAL MHD (I) Consider the stability of an infinite, homogeneous plasma: Write v 1 = ξ(r)/ t B = B 0 ez J = 0 p = p 0 ρ = ρ 0 v = 0 Expand ξ(r) as ξ(r) = ˆξe ik r Plug into linear system, and solve by hand Dynamics is anisotropic because of the magnetic field: k + k e z
43 ILLUSTRATION: WAVES IN IDEAL MHD (II) Writing the expression for each component, we get the system ω 2 k 2 v2 A ω 2 k 2 v 2 A k2 v2 S k k v 2 S k k v 2 S ω 2 k 2 v2 S ˆξ x ˆξ y = ˆξ z Two key velocities appear: v A = B 2 0 µ 0 ρ 0 v S = γ p 0 ρ 0 v A is called the Alfvén velocity, in honor of Hannes Alfvén, the Swedish scientist who first described MHD waves v s is the adiabatic sound speed
44 ILLUSTRATION: WAVES IN IDEAL MHD (III) ω 2 k 2 v2 A ω 2 k 2 v 2 A k2 v2 S k k v 2 S k k v 2 S ω 2 k 2 v2 S ˆξ x ˆξ y ˆξ z = For nontrivial solutions, determinant of the matrix should be 0 This leads to the following three possibilities for ω 2 : ω 2 = k 2 v2 A, ω2 = k2 ( k v 2 2 A + v 2 S) 1 ± 2 v A v2 S k 2 ( v 2 S + v 2 2 A) One can see that ω 2 0 The infinite homogeneous magnetized plasma is always MHD stable Some of these modes become unstable in magnetic fusion configurations, because of gradients and field line curvature
45 SHEAR ALFVÉN WAVE Branch ω 2 = k 2 v2 A (Figure from Ideal MHD by J.P. Freidberg, CUP, 2014) Transverse wave Balance between plasma inertia and field line tension Incompressible often the most unstable MHD mode in fusion devices
46 FAST MAGNETOSONIC WAVE Fast magnetosonic wave given by ω 2 = k2 ( k v 2 2 A + v 2 S) v A v2 S k 2 ( v 2 S + v 2 2 A) Simplifies in the limit v 2 S v2 A : it becomes the compressional Alfvén wave, ω 2 = k 2 v 2 A (Figure from Ideal MHD by J.P. Freidberg, CUP, 2014) Oscillation between plasma kinetic energy and magnetic compressional energy
47 SLOW MAGNETOSONIC WAVE Slow magnetosonic wave given by ω 2 = k2 ( k v 2 2 A + v 2 S) 1 2 v A v2 S k 2 ( v 2 S + v 2 2 A) Physics simplifies in the limit v 2 S v2 A : it is then called the sound wave, with dispersion relation ω 2 = k 2 v 2 S (Figure from Ideal MHD by J.P. Freidberg, CUP, 2014) Plasma motion parallel to field lines, compressible Oscillation between plasma kinetic energy and plasma internal energy (plasma pressure)
48 COMMON MHD INSTABILITIES (I) Interchange instability (Figure from Plasma Physics and Fusion Energy by J.P. Freidberg, CUP, 2008)
49 (Figure from Plasma Physics and Fusion Energy by J.P. Freidberg, CUP, 2008) COMMON MHD INSTABILITIES (II) Ballooning instability
50 COMMON MHD INSTABILITIES (III) Kink instability (Figure from Plasma Physics and Fusion Energy by J.P. Freidberg, CUP, 2008)
51 BALLOONING MODES IN TOKAMAKS (Right image from by J.W. Connor et al.)
52 KINK BALLOONING MODES IN TOKAMAKS Reconstruction of experimentally observed kink ballooning mode (Figure from U.S. Burning Plasma Organization enews February 28, 2013 (Issue 69) by S.A. Sabbagh et al.)
53 KEY ROLE OF MHD IN REACTOR DESIGN Maximum achievable pressure set by ballooning modes Maximum achievable pressure set by kink modes Maximum elongation set by vertical instability
54 POST SCRIPTUM: THE COURANT INSTITUTE OF MATHEMATICAL SCIENCES (CIMS) AT NYU
55 CIMS IN MANHATTAN
56 Abel prize in 2005, 2007, 2009 and members of the National Academy of Sciences 5 members of the National Academy of Engineering
57 Specialization in applied math, scientific computing, mathematical analysis Particular emphasis on Partial Differential Equations PhD programs in Mathematics, Atmosphere and Ocean Science, Computational Biology Masters of Science in Mathematics, Masters of Science in Scientific Computing, Masters of Science in Data Science, Masters of Science in Math Finance 60 faculty 100 PhD students
58
59 MFD DIVISION AT CIMS Founded by Harold Grad in faculty, 3 post-docs, 2 graduate students Work on MHD, wave propagation, kinetic models Analytic pen and paper work Development of new numerical solvers Collaboration with colleagues specialized in scientific computing, computational fluid dynamics, stochastic calculus, etc. Funding currently available for PhD students Feel free to contact me if you have any questions
Simple examples of MHD equilibria
Department of Physics Seminar. grade: Nuclear engineering Simple examples of MHD equilibria Author: Ingrid Vavtar Mentor: prof. ddr. Tomaž Gyergyek Ljubljana, 017 Summary: In this seminar paper I will
More 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 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 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 informationTransition From Single Fluid To Pure Electron MHD Regime Of Tearing Instability
Transition From Single Fluid To Pure Electron MHD Regime Of Tearing Instability V.V.Mirnov, C.C.Hegna, S.C.Prager APS DPP Meeting, October 27-31, 2003, Albuquerque NM Abstract In the most general case,
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 informationPlasma instabilities. Dr Ben Dudson, University of York 1 / 37
Plasma instabilities Dr Ben Dudson, University of York 1 / 37 Previously... Plasma configurations and equilibrium Linear machines, and Stellarators Ideal MHD and the Grad-Shafranov equation Collisional
More informationDispersive Media, Lecture 7 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2017-02-14 Dispersive Media, Lecture 7 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasmas as a coupled system Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas
More informationIdeal Magnetohydrodynamics (MHD)
Ideal Magnetohydrodynamics (MHD) Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 1, 2016 These lecture notes are largely based on Lectures in Magnetohydrodynamics
More informationHeating and current drive: Radio Frequency
Heating and current drive: Radio Frequency Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 13 th February 2012 Dr Ben Dudson Magnetic Confinement Fusion (1 of 26)
More 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 informationFundamentals of Magnetohydrodynamics (MHD)
Fundamentals of Magnetohydrodynamics (MHD) Thomas Neukirch School of Mathematics and Statistics University of St. Andrews STFC Advanced School U Dundee 2014 p.1/46 Motivation Solar Corona in EUV Want to
More informationMAGNETOHYDRODYNAMIC EQUILIBRIUM AND STABILITY OF PLASMA
University of Ljubljana Faculty of Mathematics and Physics Seminar 1 b -1st year, II. cycle MAGNETOHYDRODYNAMIC EQUILIBRIUM AND STABILITY OF PLASMA Author: Lino alamon Advisor: prof. dr. Tomaº Gyergyek
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 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 informationIdeal MHD Equilibria
CapSel Equil - 01 Ideal MHD Equilibria keppens@rijnh.nl steady state ( t = 0) smoothly varying solutions to MHD equations solutions without discontinuities conservative or non-conservative formulation
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 informationThe Field-Reversed Configuration (FRC) is a high-beta compact toroidal in which the external field is reversed on axis by azimuthal plasma The FRC is
and Stability of Field-Reversed Equilibrium with Toroidal Field Configurations Atomics General Box 85608, San Diego, California 92186-5608 P.O. APS Annual APS Meeting of the Division of Plasma Physics
More informationMagnetohydrodynamic waves in a plasma
Department of Physics Seminar 1b Magnetohydrodynamic waves in a plasma Author: Janez Kokalj Advisor: prof. dr. Tomaž Gyergyek Petelinje, April 2016 Abstract Plasma can sustain different wave phenomena.
More informationTwo Fluid Dynamo and Edge-Resonant m=0 Tearing Instability in Reversed Field Pinch
1 Two Fluid Dynamo and Edge-Resonant m= Tearing Instability in Reversed Field Pinch V.V. Mirnov 1), C.C.Hegna 1), S.C. Prager 1), C.R.Sovinec 1), and H.Tian 1) 1) The University of Wisconsin-Madison, Madison,
More informationTokamak Fusion Basics and the MHD Equations
MHD Simulations for Fusion Applications Lecture 1 Tokamak Fusion Basics and the MHD Equations Stephen C. Jardin Princeton Plasma Physics Laboratory CEMRACS 1 Marseille, France July 19, 21 1 Fusion Powers
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 informationEquilibrium and transport in Tokamaks
Equilibrium and transport in Tokamaks Jacques Blum Laboratoire J.-A. Dieudonné, Université de Nice Sophia-Antipolis Parc Valrose 06108 Nice Cedex 02, France jblum@unice.fr 08 septembre 2008 Jacques Blum
More informationMagnetohydrodynamic Waves
Magnetohydrodynamic Waves Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 17, 2016 These slides are largely based off of 4.5 and 4.8 of The Physics of
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 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 informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationSW103: Lecture 2. Magnetohydrodynamics and MHD models
SW103: Lecture 2 Magnetohydrodynamics and MHD models Scale sizes in the Solar Terrestrial System: or why we use MagnetoHydroDynamics Sun-Earth distance = 1 Astronomical Unit (AU) 200 R Sun 20,000 R E 1
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 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 informationKinetic, Fluid & MHD Theories
Lecture 2 Kinetic, Fluid & MHD Theories The Vlasov equations are introduced as a starting point for both kinetic theory and fluid theory in a plasma. The equations of fluid theory are derived by taking
More 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 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 information(a) (b) (c) (d) (e) (f) r (minor radius) time. time. Soft X-ray. T_e contours (ECE) r (minor radius) time time
Studies of Spherical Tori, Stellarators and Anisotropic Pressure with M3D 1 L.E. Sugiyama 1), W. Park 2), H.R. Strauss 3), S.R. Hudson 2), D. Stutman 4), X-Z. Tang 2) 1) Massachusetts Institute of Technology,
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationCalculation of alpha particle redistribution in sawteeth using experimentally reconstructed displacement eigenfunctions
Calculation of alpha particle redistribution in sawteeth using experimentally reconstructed displacement eigenfunctions R. Farengo, H. E. Ferrari,2, M.-C. Firpo 3, P. L. Garcia-Martinez 2,3, A. F. Lifschitz
More informationSensors Plasma Diagnostics
Sensors Plasma Diagnostics Ken Gentle Physics Department Kenneth Gentle RLM 12.330 k.gentle@mail.utexas.edu NRL Formulary MIT Formulary www.psfc.mit.edu/library1/catalog/ reports/2010/11rr/11rr013/11rr013_full.pdf
More information20. Alfven waves. ([3], p ; [1], p ; Chen, Sec.4.18, p ) We have considered two types of waves in plasma:
Phys780: Plasma Physics Lecture 20. Alfven Waves. 1 20. Alfven waves ([3], p.233-239; [1], p.202-237; Chen, Sec.4.18, p.136-144) We have considered two types of waves in plasma: 1. electrostatic Langmuir
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 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 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 informationLecture # 3. Introduction to Kink Modes the Kruskal- Shafranov Limit.
Lecture # 3. Introduction to Kink Modes the Kruskal- Shafranov Limit. Steve Cowley UCLA. This lecture is meant to introduce the simplest ideas about kink modes. It would take many lectures to develop the
More informationPLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT. Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford
1 PLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford 2 Overview Why do we need plasmas? For fusion, among other things
More informationTurbulence and Transport The Secrets of Magnetic Confinement
Turbulence and Transport The Secrets of Magnetic Confinement Presented by Martin Greenwald MIT Plasma Science & Fusion Center IAP January 2005 FUSION REACTIONS POWER THE STARS AND PRODUCE THE ELEMENTS
More informationIntroduction to Magnetohydrodynamics (MHD)
Introduction to Magnetohydrodynamics (MHD) Tony Arber University of Warwick 4th SOLARNET Summer School on Solar MHD and Reconnection Aim Derivation of MHD equations from conservation laws Quasi-neutrality
More informationExponential Growth and Filamentary Structure of Nonlinear Ballooning Instability 1
Exponential Growth and Filamentary Structure of Nonlinear Ballooning Instability 1 Ping Zhu in collaboration with C. C. Hegna and C. R. Sovinec University of Wisconsin-Madison Sherwood Conference Denver,
More informationPLASMA ASTROPHYSICS. ElisaBete M. de Gouveia Dal Pino IAG-USP. NOTES: (references therein)
PLASMA ASTROPHYSICS ElisaBete M. de Gouveia Dal Pino IAG-USP NOTES:http://www.astro.iag.usp.br/~dalpino (references therein) ICTP-SAIFR, October 7-18, 2013 Contents What is plasma? Why plasmas in astrophysics?
More informationReduced MHD. Nick Murphy. Harvard-Smithsonian Center for Astrophysics. Astronomy 253: Plasma Astrophysics. February 19, 2014
Reduced MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 19, 2014 These lecture notes are largely based on Lectures in Magnetohydrodynamics by Dalton
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 informationAnalysis of Jeans Instability of Partially-Ionized. Molecular Cloud under Influence of Radiative. Effect and Electron Inertia
Adv. Studies Theor. Phys., Vol. 5, 2011, no. 16, 755-764 Analysis of Jeans Instability of Partially-Ionized Molecular Cloud under Influence of Radiative Effect and Electron Inertia B. K. Dangarh Department
More informationIssues in Neoclassical Tearing Mode Theory
Issues in Neoclassical Tearing Mode Theory Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Austin, TX Tearing Mode Stability in Tokamaks According to standard (single-fluid)
More informationProgressing Performance Tokamak Core Physics. Marco Wischmeier Max-Planck-Institut für Plasmaphysik Garching marco.wischmeier at ipp.mpg.
Progressing Performance Tokamak Core Physics Marco Wischmeier Max-Planck-Institut für Plasmaphysik 85748 Garching marco.wischmeier at ipp.mpg.de Joint ICTP-IAEA College on Advanced Plasma Physics, Triest,
More informationThe RFP: Plasma Confinement with a Reversed Twist
The RFP: Plasma Confinement with a Reversed Twist JOHN SARFF Department of Physics University of Wisconsin-Madison Invited Tutorial 1997 Meeting APS DPP Pittsburgh Nov. 19, 1997 A tutorial on the Reversed
More informationModels for Global Plasma Dynamics
Models for Global Plasma Dynamics F.L. Waelbroeck Institute for Fusion Studies, The University of Texas at Austin International ITER Summer School June 2010 Outline 1 Models for Long-Wavelength Plasma
More informationMHD Linear Stability Analysis Using a Full Wave Code
US-Japan JIFT Workshop on Progress of Extended MHD Models NIFS, Toki,Japan 2007/03/27 MHD Linear Stability Analysis Using a Full Wave Code T. Akutsu and A. Fukuyama Department of Nuclear Engineering, Kyoto
More informationPHYS 643 Week 4: Compressible fluids Sound waves and shocks
PHYS 643 Week 4: Compressible fluids Sound waves and shocks Sound waves Compressions in a gas propagate as sound waves. The simplest case to consider is a gas at uniform density and at rest. Small perturbations
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 informationMHD RELATED TO 2-FLUID THEORY, KINETIC THEORY AND MAGANETIC RECONNECTION
MHD RELATED TO 2-FLUID THEORY, KINETIC THEORY AND MAGANETIC RECONNECTION Marty Goldman University of Colorado Spring 2017 Physics 5150 Issues 2 How is MHD related to 2-fluid theory Level of MHD depends
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 informationBasic plasma physics
Basic plasma physics SPAT PG Lectures Jonathan Eastwood 10-14 October 2016 Aims Provide new PhD students in SPAT and the SPC section with an overview of the most important principles in space plasma physics,
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 informationIntroduction to Fusion Physics
Introduction to Fusion Physics Hartmut Zohm Max-Planck-Institut für Plasmaphysik 85748 Garching DPG Advanced Physics School The Physics of ITER Bad Honnef, 22.09.2014 Energy from nuclear fusion Reduction
More informationIssues of Perpendicular Conductivity and Electric Fields in Fusion Devices
Issues of Perpendicular Conductivity and Electric Fields in Fusion Devices Michael Tendler, Alfven Laboratory, Royal Institute of Technology, Stockholm, Sweden Plasma Turbulence Turbulence can be regarded
More informationFast Secondary Reconnection and the Sawtooth Crash
Fast Secondary Reconnection and the Sawtooth Crash Maurizio Ottaviani 1, Daniele Del Sarto 2 1 CEA-IRFM, Saint-Paul-lez-Durance (France) 2 Université de Lorraine, Institut Jean Lamour UMR-CNRS 7198, Nancy
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 informationSawteeth in Tokamaks and their relation to other Two-Fluid Reconnection Phenomena
Sawteeth in Tokamaks and their relation to other Two-Fluid Reconnection Phenomena S. C. Jardin 1, N. Ferraro 2, J. Chen 1, et al 1 Princeton Plasma Physics Laboratory 2 General Atomics Supported by the
More informationInnovative Concepts Workshop Austin, Texas February 13-15, 2006
Don Spong Oak Ridge National Laboratory Acknowledgements: Jeff Harris, Hideo Sugama, Shin Nishimura, Andrew Ware, Steve Hirshman, Wayne Houlberg, Jim Lyon Innovative Concepts Workshop Austin, Texas February
More informationGeneralized Solovev equilibrium with sheared flow of arbitrary direction and stability consideration
Generalized Solovev equilibrium with sheared flow of arbitrary direction and stability consideration D.A. Kaltsas and G.N. Throumoulopoulos Department of Physics, University of Ioannina, GR 451 10 Ioannina,
More informationSelf-consistent modeling of ITER with BALDUR integrated predictive modeling code
Self-consistent modeling of ITER with BALDUR integrated predictive modeling code Thawatchai Onjun Sirindhorn International Institute of Technology, Thammasat University, Klong Luang, Pathumthani, 12121,
More informationMHD Modes of Solar Plasma Structures
PX420 Solar MHD 2013-2014 MHD Modes of Solar Plasma Structures Centre for Fusion, Space & Astrophysics Wave and oscillatory processes in the solar corona: Possible relevance to coronal heating and solar
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 informationINTERACTION OF DRIFT WAVE TURBULENCE AND MAGNETIC ISLANDS
INTERACTION OF DRIFT WAVE TURBULENCE AND MAGNETIC ISLANDS A. Ishizawa and N. Nakajima National Institute for Fusion Science F. L. Waelbroeck, R. Fitzpatrick, W. Horton Institute for Fusion Studies, University
More informationComputational modeling of fully ionized magnetized plasmas using the fluid approximation a
PHYSICS OF PLASMAS 13, 058103 2006 Computational modeling of fully ionized magnetized plasmas using the fluid approximation a D. D. Schnack b Center for Energy Space Science, Science Applications International
More informationPerformance limits. Ben Dudson. 24 th February Department of Physics, University of York, Heslington, York YO10 5DD, UK
Performance limits Ben Dudson Department of Physics, University of York, Heslington, York YO10 5DD, UK 24 th February 2014 Ben Dudson Magnetic Confinement Fusion (1 of 24) Previously... In the last few
More informationRecapitulation: Questions on Chaps. 1 and 2 #A
Recapitulation: Questions on Chaps. 1 and 2 #A Chapter 1. Introduction What is the importance of plasma physics? How are plasmas confined in the laboratory and in nature? Why are plasmas important in astrophysics?
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 informationMagnetic Reconnection Controlled by Multi-Hierarchy Physics in an Open System
Magnetic Reconnection Controlled by Multi-Hierarchy Physics in an Open System Ritoku HORIUCHI 1,2), Shunsuke USAMI 1), Hiroaki OHTANI 1,2) and Toseo MORITAKA 3) 1) National Institute for Fusion Science,
More informationSimulation of alpha particle current drive and heating in spherical tokamaks
Simulation of alpha particle current drive and heating in spherical tokamaks R. Farengo 1, M. Zarco 1, H. E. Ferrari 1, 1 Centro Atómico Bariloche and Instituto Balseiro, Argentina. Consejo Nacional de
More informationExponential Growth of Nonlinear Ballooning Instability. Abstract
Exponential Growth of Nonlinear Ballooning Instability P. Zhu, C. C. Hegna, and C. R. Sovinec Center for Plasma Theory and Computation University of Wisconsin-Madison Madison, WI 53706, USA Abstract Recent
More informationHigh-m Multiple Tearing Modes in Tokamaks: MHD Turbulence Generation, Interaction with the Internal Kink and Sheared Flows
TH/P3-3 High-m Multiple Tearing Modes in Tokamaks: MHD Turbulence Generation, Interaction with the Internal Kink and Sheared Flows A. Bierwage 1), S. Benkadda 2), M. Wakatani 1), S. Hamaguchi 3), Q. Yu
More informationLinear and non-linear evolution of the gyroresonance instability in Cosmic Rays
Linear and non-linear evolution of the gyroresonance instability in Cosmic Rays DESY Summer Student Programme, 2016 Olga Lebiga Taras Shevchenko National University of Kyiv, Ukraine Supervisors Reinaldo
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 informationTheory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus. K. C. Shaing
Theory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus K. C. Shaing Plasma and Space Science Center, and ISAPS, National Cheng Kung University, Tainan, Taiwan 70101, Republic
More informationModeling of ELM Dynamics for ITER
Modeling of ELM Dynamics for ITER A.Y. PANKIN 1, G. BATEMAN 1, D.P. BRENNAN 2, A.H. KRITZ 1, S. KRUGER 3, P.B. SNYDER 4 and the NIMROD team 1 Lehigh University, 16 Memorial Drive East, Bethlehem, PA 18015
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 informationLinear stability of MHD configurations
Linear stability of MHD configurations Rony Keppens Centre for mathematical Plasma Astrophysics KU Leuven Rony Keppens (KU Leuven) Linear MHD stability CHARM@ROB 2017 1 / 18 Ideal MHD configurations Interested
More informationSpace Physics. ELEC-E4520 (5 cr) Teacher: Esa Kallio Assistant: Markku Alho and Riku Järvinen. Aalto University School of Electrical Engineering
Space Physics ELEC-E4520 (5 cr) Teacher: Esa Kallio Assistant: Markku Alho and Riku Järvinen Aalto University School of Electrical Engineering The 6 th week: topics Last week: Examples of waves MHD: Examples
More informationIntegrated Simulation of ELM Energy Loss Determined by Pedestal MHD and SOL Transport
1 Integrated Simulation of ELM Energy Loss Determined by Pedestal MHD and SOL Transport N. Hayashi, T. Takizuka, T. Ozeki, N. Aiba, N. Oyama Japan Atomic Energy Agency, Naka, Ibaraki-ken, 311-0193 Japan
More informationCollisionless nonideal ballooning modes
PHYSICS OF PLASMAS VOLUME 6, NUMBER 1 JANUARY 1999 Collisionless nonideal ballooning modes Robert G. Kleva and Parvez N. Guzdar Institute for Plasma Research, University of Maryland, College Park, Maryland
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 informationOn existence of resistive magnetohydrodynamic equilibria
arxiv:physics/0503077v1 [physics.plasm-ph] 9 Mar 2005 On existence of resistive magnetohydrodynamic equilibria H. Tasso, G. N. Throumoulopoulos Max-Planck-Institut für Plasmaphysik Euratom Association
More informationEnergetic Particle Physics in Tokamak Burning Plasmas
Energetic Particle Physics in Tokamak Burning Plasmas presented by C. Z. (Frank) Cheng in collaboration with N. N. Gorelenkov, G. J. Kramer, R. Nazikian, E. Fredrickson, Princeton Plasma Physics Laboratory
More informationOverview of FRC-related modeling (July 2014-present)
Overview of FRC-related modeling (July 2014-present) Artan Qerushi AFRL-UCLA Basic Research Collaboration Workshop January 20th, 2015 AFTC PA Release# 15009, 16 Jan 2015 Artan Qerushi (AFRL) FRC modeling
More informationProf. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit
Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit Rough breakdown of MHD shocks Jump conditions: flux in = flux out mass flux: ρv n magnetic flux: B n Normal momentum flux: ρv n
More informationMHD equilibrium calculations for stellarators. Antoine Cerfon, MIT PSFC with F. Parra, J. Freidberg (MIT NSE)
MHD equilibrium calculations for stellarators Antoine Cerfon, MIT PSFC with F. Parra, J. Freidberg (MIT NSE) March 20, 2012 MAGNETIC FIELD LINE HAMILTONIAN Consider a general toroidal coordinate system
More informationMHD. Jeff Freidberg MIT
MHD Jeff Freidberg MIT 1 What is MHD MHD stands for magnetohydrodynamics MHD is a simple, self-consistent fluid description of a fusion plasma Its main application involves the macroscopic equilibrium
More informationA Study of 3-Dimensional Plasma Configurations using the Two-Fluid Plasma Model
A Study of 3-Dimensional Plasma Configurations using the Two-Fluid Plasma Model B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program University of Washington IEEE International Conference
More informationHeat Transport in a Stochastic Magnetic Field. John Sarff Physics Dept, UW-Madison
Heat Transport in a Stochastic Magnetic Field John Sarff Physics Dept, UW-Madison CMPD & CMSO Winter School UCLA Jan 5-10, 2009 Magnetic perturbations can destroy the nested-surface topology desired for
More informationHighlights from (3D) Modeling of Tokamak Disruptions
Highlights from (3D) Modeling of Tokamak Disruptions Presented by V.A. Izzo With major contributions from S.E. Kruger, H.R. Strauss, R. Paccagnella, MHD Control Workshop 2010 Madison, WI ..onset of rapidly
More information