# Lecture # 3. Introduction to Kink Modes the Kruskal- Shafranov Limit.

Save this PDF as:

Size: px
Start display at page:

Download "Lecture # 3. Introduction to Kink Modes the Kruskal- Shafranov Limit."

## Transcription

1 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 theory in all its detail, but hopefully the key ideas can be communicated in 1 hour. We will consider a very simple zero pressure cylindrical equilibrium with nearly constant current in the z direction. When the current is ramped up the plasma becomes unstable to a helical kinking of the plasma column. The dynamics of this instability is our topic. These notes provide considerable detail a lot of this detail is impossible to convey in the lecture but I hope you will read the notes at your leisure. A good exposition of kink modes is given in J.P. Freiberg Ideal Magneto-Hydrodynamics Plenum (1987). and in H. Goedbloed and S. Poedts Principles of Magneto-hydrodynamics Cambridge (2004). 1 Equilibrium We consider an equilibrium with constant current in the z direction J z = J 0 =constant for r < a. The basic expansion parameter is inverse aspect ratio, ɛ = a/r 0 1. The magnetic field is given by: B 0 = (B 0 + B 2 (r))ẑ + B p (r) with B 2 ɛ 2 B 0 and B p ɛb 0 (1) Where B 0 is a constant the dominant toroidal field in a tokamak and r is the cylindrical radius. This ordering of the equilibrium is sometimes called the tokamak ordering since it is roughly true for many tokamak experiments. From B 0 = µ 0 J 0 we obtain (for r < a): B p = µ 0J 0 2 (ẑ r) = B 0 (ẑ r) = B p (ẑ r). (2) qr 0 r where q is constant (see below) and using the force balance equation we obtain to O(ɛ 2 ), The field lines obey the equation: J 0 B 0 = 0 (3) B 2 = B 0 (a 2 r 2 ). (4) q 2 R0 2 dθ dz = B p rb 0 = 1 qr 0 or equivalently θ = z qr 0 + θ 0. (5)

2 2πR 0 a b Figure 1: Cylindrical Screw-pinch The plasma is contained in the region r < a and the wall at r = b is considered perfectly conducting. The region between r = a and r = b is the vacuum region. In the model considered here the (black) field lines in the plasma all have the same rotational transform. The blue line is the vacuum field at the wall. The two ends of the pinch are identified i.e. it is periodic in z. Thus all the plasma field lines have the same twist every time we go from one end of the cylinder to the other the plasma field lines go around 1/q times in the θ poloidal direction. In tokamaks q is called the safety factor. In the Fig. (1) the q is 0.4. Clearly the greater the current the smaller the q an upper limit on the current is then a lower limit on q. In the vacuum region the current is, of course, zero and the magnetic field is: B 0 = B 0 ẑ + B 0a 2 The vacuum field has a q profile q v (r) = qr 2 /a 2. qr 2 R 0 (ẑ r). (6)

3 2 Plasma Perturbation To investigate stability we imagine pushing the plasma a small distance away from its equilibrium state. Does the plasma move further away from or return to its equilibrium state? We model the plasma dynamics with Ideal Magneto-hydrodynamics (MHD). The magnetic field obeys: the frozen-in law which when linearized yields: B t = (v B) (7) δb = (ξ B 0 ) = B 0 ξ ξ B 0 B 0 ξ. (8) where ξ is the displacement of the plasma ( ξ/ t = v). The perturbed J B force accelerates the displacement. Thus the momentum equation becomes: ρ 0 2 ξ t 2 = J 0 δb + ( δb) B 0 = F(ξ). (9) The plasma density ρ 0 is taken to be constant and the pressure is set to zero. If we were doing this formally we would now expand the components of this equation and solve the eigenvalue problem for the growth rate. However we can get to the answer quicker by using our intuition to figure out a good first approximation to the unstable perturbation. We shall look at helical perturbations that move the column. Simply compressing the strong z (toroidal) field makes a stable compressional Alfvén wave. Thus we take the displacement to be dominantly a rigid helical shift of the plasma column, i.e. ξ 0 = ξ 0 [cos ( z )ˆx + sin ( z ] )ŷ = ξ 0 ˆξ0. (10) R 0 R 0 Each poloidal (x, y) plane of the plasma is shifted without distortiona constant distance ξ 0. The radial displacement is ξ r0 = ξ 0 ˆr = ξ 0 cos (θ z R 0 ) (11) This is the usual way to write the cylindrical displacement. We often expand as helical modes with displacements proportional to exp (imθ in z R 0 ) here (clearly) we are examining m = n = 1 modes. The field and current perturbations from the rigid displacement are from Eq. (14), using Eqs. (2) and (4): δb 0 = B [ 0 (1 q)(ξ qr 0 ẑ) + 2ξ ] 0 r ẑ 0 qr 0 δj 0 = B [ 0 (2 q + q 2 )(ξ µ 0 q 2 R0 2 0 ẑ) ] (12)

4 2πR 0 ξ Figure 2: Helical m=1 n=1 Kink displacement Rigid displacement of each poloidal plane of the plasma by displacement of the form given by Eq. (10). The first term in he field perturbation comes from bending the z (toroidal) field. The second term is from moving the varying part toroidal field (B 2 ). Thus the perturbed force from this displacement is (using Eq. (9)) is given by F 0 = B2 0 µ 0 (1 + q)ξ qr0 2 0 this is a stabilizing force since it opposes the motion! This displacement is not quite right. Outside the plasma the vacuum field is perturbed and produces a pressure on the outside of the plasma this pressure will vary in θ in the same way as the perturbation. We must match the magnetic pressure just inside the plasma with the magnetic pressure just outside the plasma by slightly compressing the plasma. Thus we need to add to our plasma displacement a small compressive component whose magnitude will be determined by the matching of pressures at the boundary. For this compressive perturbation to produce a constant force on the plasma and match the pressure from the vacuum at the boundary it must be quadratic in r and point in the direction of ξ 0. Thus we set: ξ 2 = α r2 ξ 2qR (13) We shall determine the constant α so that the magnetic pressures balance at the boundary. α is O(1) so that the compressional displacement is small, ξ 2 O(ɛ 2 )ξ 0. With this small compressional displacement the perturbed field and current become:

5 δb = B [ 0 (1 q)(ξ qr 0 ẑ) + (2 αq)ξ ] 0 r ẑ 0 qr 0 δj = B [ 0 (2 qα q + q 2 )(ξ µ 0 q 2 R0 2 0 ẑ) ] (14) We keep the perturbed z magnetic field to O(ɛ 2 ) and the perturbed poloidal field to O(ɛ) since they contribute to the magnetic pressure and the forces at the same order. The compression contributes through compressing the z (toroidal) field. After a little algebra the momentum equation becomes: 2 ξ ρ 0 0 = F(ξ) = B2 0 (1 + q α)ξ t 2 µ 0 qr (15) 3 Vacuum Perturbation In the vacuum δb = 0 so we set δb = χ. From δb = 0 we obtain 2 χ = 0. Taking χ = ˆχ(r) sin (θ z R 0 ) we find: 1 r ( d dr r dˆχ dr ) ˆχ r 2 ˆχ R 2 0 = 0 (16) We can drop the third term on the left hand side of this equation to lowest order and obtain: χ = (C 1 r + C 2 r + O(ɛ2 )) sin (θ z R 0 ) (17) Where C 1 and C 2 are constants to be determined. The perturbed vacuum field is: δb v = (C 1 C 2 r 2 ) sin (θ z R 0 )ˆr + (C 1 r + C 2 r ) cos (θ z R 0 )[ 1 r ẑ ˆr 1 R 0 ẑ] (18) At the conducting wall (r = b) we can have no radial field thus C 2 = b 2 C 1.

6 n a ξ r ξ Figure 3: Vacuum-Plasma Interface. The shift of the boundary by the displacement to the order we keep the displacement is just ξ 0. The boundary is given by the line r = a + ξ r and the unit vector normal to the boundary is n. 4 Matching at the Vacuum-Plasma Interface The normal to the perturbed flux surfaces (labeled by their original radius r 0 ) is in the direction of the r 0 = (r ξ r ). Then to the order we need: n = ˆr ξ r (19) To join the plasma and vacuum solutions we must satisfy two conditions: 1) no field through the vacuum-plasma interface and, 2) magnetic pressure continuity across the vacuum-plasma interface. No Vacuum Field Through the Vacuum-Plasma Interface. (B 0 + δb v ) (ˆr ξ r ) = 0 B 0 ξ r = δb v ˆr (20) Substituting into this expression we determine C 1 in terms of ξ 0 : C 1 = ξ 0B 0 (1 1) q R 0 (1 b2 a 2 ) (21) Magnetic Pressure Continuity Across the Vacuum-Plasma Interface. Force balance across the plasma boundary becomes the continuity of magnetic pressure.

7 Thus on the surface r = r 0 + ξ r we have: B 2 plasma = B 2 vacuum (22) Using Eq. (14) and a little algebra we find: B 2 plasma = B a2 B 2 0 q 2 R 2 0 and from Eqs. (18) and (21) we obtain: B 2 vacuum = B a2 B 2 0 q 2 R 2 0 2B2 0ξ r qr B2 0ξ r (1 α) (23) qr0 2 1 q q(1 1 q b2 ( )2 ( b2 + 1) a 2 (24) 1) a 2 Equating these expressions we obtain α = [1 + 1 q(1 1 q q )2 ( b 2 a 2 +1) Eq. (15) we obtain: ( b2 a 2 1)] and substituting into 2 ξ ρ 0 0 = F(ξ) = 2B2 ( b2 1) 0 a 2 q t 2 µ 0 qr0 2 ( b2 1) (1 q 1)ξ 0. (25) a 2 Clearly we have a growing mode when 1 > q > a 2 /b 2 with growth rate, γ = 2B2 ( b2 1) 0 a 2 q µ 0 ρ 0 qr0 2 ( b2 1) (1 1). (26) q a 2 When we set the conducting wall to go to infinity we get the Kruskal-Shafranov limit for stability: q > 1 Stable. (27) This is a limit on the current, specifically we must have µ 0 J 0 < 2B 0 /R 0 to be stable. Although the stability criterion is very simple all the regions and forces play a role in this instability. Specifically note that the field line bending forces in the plasma the 1 + q part of the right hand side of Eq. (15) are stabilizing. The vacuum pressure is destabilizing since where ξ r > 0 the magnetic pressure on the boundary has decreased.

8 This arises because the vacuum field decreases with distance from the axis (as 1/r). The decrease in pressure on the boundary is balanced in the plasma by a reduced z field caused by an expansion of the plasma (where the radial displacement is positive ξ r > 0) this corresponds to (α > 0). The magnetic pressure forces from the expansion/compression inside the plasma are then destabilizing. When q < 1 these pressure forces exceed the field line bending forces and the plasma is unstable. Thus this is an external kink releasing energy in the vacuum field indeed the vacuum field lines straighten. When the wall is close there is a narrow region of instability which corresponds to having the q = 1 radius in the vacuum region i.e. q v (r 1 ) = 1 = qr 2 1/a 2. When the wall is present the flux trapped between the plasma and the wall cannot escape and makes the vacuum field pressure less destabilizing. 5 Other Kink Modes In Tokamaks kink modes are driven by both current and pressure (Professor Brennan will describe the situation in some detail) and the plasma is not cylindrical. Although the m = n = 1 is often the most visible mode number the poloidal (m) harmonics are coupled and none of the modes are a pure poloidal harmonic. The external kinks that are observed can be predominantly higher harmonics. Often what one sees is a mode with m/n q edge where q edge is the value of q at the edge of the plasma. It is instructive to draw the kinked boundary for the m = 2 mode. Figure 4: Helical m=2 n=1 Kink displacement Eliptical displacement of each flux surface in the plasma by displacement of the form given by ξ r cos (2θ z R 0 ).

9 6 Problem Vertical and Horizontal Stability. Tokamaks have to be stable to moving up/down or side to side. This puts constraints on the coil arrangement next to the plasma. Lets look at the stability of current filaments as a model of this situation. Consider 5 current filaments (current flowing in the z direction): Plasma Filament initially at x = y = 0 representing the plasma carrying current I plasma 2 Horizontal Coils at x = ±1, y = 0 carrying current I 0 2 Vertical Coils at y = ±1, x = 0 carrying current I 0 Let us imagine that the plasma and coil currents are fixed i.e. they do not change during any motion of the plasma. 1. Draw the configuration. 2. Calculate the force on the undisplaced plasma? 3. Calculate the force on the plasma when it is displaced by a distance δ horizontally. 4. Calculate the force on the plasma when it is displaced by a distance δ vertically. 5. In which direction is it unstable? What is the growth rate if the plasma has a mass M? 6. How could you stabilize this instability? Discuss!

### The 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

### The 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

### Magnetohydrodynamic 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

### Introduction 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

### The 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

### Physics 505 Fall 2005 Homework Assignment #7 Solutions

Physics 505 Fall 005 Homework Assignment #7 Solutions Textbook problems: Ch. 4: 4.10 Ch. 5: 5.3, 5.6, 5.7 4.10 Two concentric conducting spheres of inner and outer radii a and b, respectively, carry charges

### Lecture IV: Magnetized plasma stability

Outline Lecture IV: Magnetized plasma stability Wojciech Fundamenski UKAEA Fusion, Culham Science Centre, Abingdon and Plasma Physics Group, Imperial College, London 22nd May 2009 Outline Part I: 1 Hydrodynamic

### 22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 19

. tability of the straight tokamak.65, MHD Theory of Fusion ystems Prof. Freidberg Lecture 9. ressure driven modes (uydams Criterion). internal modes 3. external modes. Tokamak Ordering Bθ ar B μ q or

### Bifurcated states of a rotating tokamak plasma in the presence of a static error-field

Bifurcated states of a rotating tokamak plasma in the presence of a static error-field Citation: Physics of Plasmas (1994-present) 5, 3325 (1998); doi: 10.1063/1.873000 View online: http://dx.doi.org/10.1063/1.873000

### arxiv: v2 [physics.plasm-ph] 15 Nov 2016

arxiv:1608.04939v2 [physics.plasm-ph] 15 Nov 2016 Spectrum of multi-region-relaxed magnetohydrodynamic modes in topologically toroidal geometry R. L. Dewar, L. H. Tuen, M. J. Hole Centre for Plasmas and

### Prof. 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

### 6.1. Linearized Wave Equations in a Uniform Isotropic MHD Plasma. = 0 into Ohm s law yields E 0

Chapter 6. Linear Waves in the MHD Plasma 85 Chapter 6. Linear Waves in the MHD Plasma Topics or concepts to learn in Chapter 6:. Linearize the MHD equations. The eigen-mode solutions of the MHD waves

### PHYS 432 Physics of Fluids: Instabilities

PHYS 432 Physics of Fluids: Instabilities 1. Internal gravity waves Background state being perturbed: A stratified fluid in hydrostatic balance. It can be constant density like the ocean or compressible

### Momentum transport from magnetic reconnection in laboratory an. plasmas. Fatima Ebrahimi

Momentum transport from magnetic reconnection in laboratory and astrophysical plasmas Space Science Center - University of New Hampshire collaborators : V. Mirnov, S. Prager, D. Schnack, C. Sovinec Center

### Stellarators. Dr Ben Dudson. 6 th February Department of Physics, University of York Heslington, York YO10 5DD, UK

Stellarators Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 6 th February 2014 Dr Ben Dudson Magnetic Confinement Fusion (1 of 23) Previously... Toroidal devices

### Modelling of the penetration process of externally applied helical magnetic perturbation of the DED on the TEXTOR tokamak

INSTITUTE OF PHYSICS PUBLISHING Plasma Phys. Control. Fusion 8 (6) 69 8 PLASMA PHYSICS AND CONTROLLED FUSION doi:.88/7-/8// Modelling of the penetration process of externally applied helical magnetic perturbation

### Toroidal confinement of non-neutral plasma. Martin Droba

Toroidal confinement of non-neutral plasma Martin Droba Contents Experiments with toroidal non-neutral plasma Magnetic surfaces CNT and IAP-high current ring Conclusion 2. Experiments with toroidal non-neutral

### SMR/ 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

### arxiv: v1 [physics.plasm-ph] 11 Mar 2016

1 Effect of magnetic perturbations on the 3D MHD self-organization of shaped tokamak plasmas arxiv:1603.03572v1 [physics.plasm-ph] 11 Mar 2016 D. Bonfiglio 1, S. Cappello 1, M. Veranda 1, L. Chacón 2 and

### Active and Fast Particle Driven Alfvén Eigenmodes in Alcator C-Mod

Active and Fast Particle Driven Alfvén Eigenmodes in Alcator C-Mod JUST DID IT. J A Snipes, N Basse, C Boswell, E Edlund, A Fasoli #, N N Gorelenkov, R S Granetz, L Lin, Y Lin, R Parker, M Porkolab, J

### Simulations of Sawteeth in CTH. Nicholas Roberds August 15, 2015

Simulations of Sawteeth in CTH Nicholas Roberds August 15, 2015 Outline Problem Description Simulations of a small tokamak Simulations of CTH 2 Sawtoothing Sawtoothing is a phenomenon that is seen in all

### Control of resistive wall modes in a cylindrical tokamak with plasma rotation and complex gain

Control of resistive wall modes in a cylindrical tokamak with plasma rotation and complex gain arxiv:146.5245v1 [physics.plasm-ph] 2 Jun 214 D. P. Brennan and J. M. Finn June 23, 214 Department of Astrophysical

### SOLAR MHD Lecture 2 Plan

SOLAR MHD Lecture Plan Magnetostatic Equilibrium ü Structure of Magnetic Flux Tubes ü Force-free fields Waves in a homogenous magnetized medium ü Linearized wave equation ü Alfvén wave ü Magnetoacoustic

### Amplification of magnetic fields in core collapse

Amplification of magnetic fields in core collapse Miguel Àngel Aloy Torás, Pablo Cerdá-Durán, Thomas Janka, Ewald Müller, Martin Obergaulinger, Tomasz Rembiasz Universitat de València; Max-Planck-Institut

### Heat 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

### Plasma models for the design of the ITER PCS

Plasma models for the design of the ITER PCS G. De Tommasi 1,2 on behalf of the CREATE team 1 Consorzio CREATE, Naples, Italy 2 Department of Electrical Engineering and Information Technology, University

### INITIAL EVALUATION OF COMPUTATIONAL TOOLS FOR STABILITY OF COMPACT STELLARATOR REACTOR DESIGNS

INITIAL EVALUATION OF COMPUTATIONAL TOOLS FOR STABILITY OF COMPACT STELLARATOR REACTOR DESIGNS A.D. Turnbull and L.L. Lao General Atomics (with contributions from W.A. Cooper and R.G. Storer) Presentation

### Measuring from electron temperature fluctuations in the Tokamak Fusion Test Reactor

PHYSICS OF PLASMAS VOLUME 5, NUMBER FEBRUARY 1998 Measuring from electron temperature fluctuations in the Tokamak Fusion Test Reactor C. Ren, a) J. D. Callen, T. A. Gianakon, and C. C. Hegna University

### ELM control with RMP: plasma response models and the role of edge peeling response

ELM control with RMP: plasma response models and the role of edge peeling response Yueqiang Liu 1,2,3,*, C.J. Ham 1, A. Kirk 1, Li Li 4,5,6, A. Loarte 7, D.A. Ryan 8,1, Youwen Sun 9, W. Suttrop 10, Xu

### 7 Curvilinear coordinates

7 Curvilinear coordinates Read: Boas sec. 5.4, 0.8, 0.9. 7. Review of spherical and cylindrical coords. First I ll review spherical and cylindrical coordinate systems so you can have them in mind when

### Tokamak 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

### Figure 11.1: A fluid jet extruded where we define the dimensionless groups

11. Fluid Jets 11.1 The shape of a falling fluid jet Consider a circular orifice of a radius a ejecting a flux Q of fluid density ρ and kinematic viscosity ν (see Fig. 11.1). The resulting jet accelerates

### Toroidal 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

### GA A26887 ADVANCES TOWARD QH-MODE VIABILITY FOR ELM-FREE OPERATION IN ITER

GA A26887 ADVANCES TOWARD QH-MODE VIABILITY FOR ELM-FREE OPERATION IN ITER by A.M. GAROFALO, K.H. BURRELL, M.J. LANCTOT, H. REIMERDES, W.M. SOLOMON and L. SCHMITZ OCTOBER 2010 DISCLAIMER This report was

### STABILIZATION OF m=2/n=1 TEARING MODES BY ELECTRON CYCLOTRON CURRENT DRIVE IN THE DIII D TOKAMAK

GA A24738 STABILIZATION OF m=2/n=1 TEARING MODES BY ELECTRON CYCLOTRON CURRENT DRIVE IN THE DIII D TOKAMAK by T.C. LUCE, C.C. PETTY, D.A. HUMPHREYS, R.J. LA HAYE, and R. PRATER JULY 24 DISCLAIMER This

### Generalized 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,

### arxiv:physics/ v1 [physics.plasm-ph] 14 Nov 2005

arxiv:physics/0511124v1 [physics.plasm-ph] 14 Nov 2005 Early nonlinear regime of MHD internal modes: the resistive case M.C. Firpo Laboratoire de Physique et Technologie des Plasmas (C.N.R.S. UMR 7648),

### Extended Lumped Parameter Model of Resistive Wall Mode and The Effective Self-Inductance

Extended Lumped Parameter Model of Resistive Wall Mode and The Effective Self-Inductance M.Okabayashi, M. Chance, M. Chu* and R. Hatcher A. Garofalo**, R. La Haye*, H. Remeirdes**, T. Scoville*, and T.

### Homework Assignment 4 Solution Set

Homework Assignment 4 Solution Set PHYCS 442 7 February, 24 Problem (Griffiths 2.37 If the plates are sufficiently large the field near them does not depend on d. The field between the plates is zero (the

### INGENIERÍA EN NANOTECNOLOGÍA

ETAPA DISCIPLINARIA TAREAS 385 TEORÍA ELECTROMAGNÉTICA Prof. E. Efren García G. Ensenada, B.C. México 206 Tarea. Two uniform line charges of ρ l = 4 nc/m each are parallel to the z axis at x = 0, y = ±4

### 3 Hydrostatic Equilibrium

3 Hydrostatic Equilibrium Reading: Shu, ch 5, ch 8 31 Timescales and Quasi-Hydrostatic Equilibrium Consider a gas obeying the Euler equations: Dρ Dt = ρ u, D u Dt = g 1 ρ P, Dɛ Dt = P ρ u + Γ Λ ρ Suppose

### SW103: 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

### generalfusion Characterization of General Fusion's Plasma Devices 2015 Nimrod Summer Workshop

Characterization of General Fusion's Plasma Devices 2015 Nimrod Summer Workshop Aaron Froese, Charlson Kim, Meritt Reynolds, Sandra Barsky, Victoria Suponitsky, Stephen Howard, Russ Ivanov, Peter O'Shea,

### The Madison Dynamo Experiment: magnetic instabilities driven by sheared flow in a sphere. Cary Forest Department of Physics University of Wisconsin

The Madison Dynamo Experiment: magnetic instabilities driven by sheared flow in a sphere Cary Forest Department of Physics University of Wisconsin February 28, 2001 Planets, stars and perhaps the galaxy

### Physics 221. Exam III Spring f S While the cylinder is rolling up, the frictional force is and the cylinder is rotating

Physics 1. Exam III Spring 003 The situation below refers to the next three questions: A solid cylinder of radius R and mass M with initial velocity v 0 rolls without slipping up the inclined plane. N

### xy 2 e 2z dx dy dz = 8 3 (1 e 4 ) = 2.62 mc. 12 x2 y 3 e 2z 2 m 2 m 2 m Figure P4.1: Cube of Problem 4.1.

Problem 4.1 A cube m on a side is located in the first octant in a Cartesian coordinate system, with one of its corners at the origin. Find the total charge contained in the cube if the charge density

### Plasma Stability in Tokamaks and Stellarators

Plasma Stability in Tokamaks and Stellarators Gerald A. Navratil GCEP Fusion Energy Workshop Princeton, NJ 1- May 006 ACKNOWLEDGEMENTS Borrowed VGs from many colleagues: J. Bialek, A. Garofalo,R. Goldston,

### BRIEF COMMUNICATION. Near-magnetic-axis Geometry of a Closely Quasi-Isodynamic Stellarator. Greifswald, Wendelsteinstr. 1, Greifswald, Germany

BRIEF COMMUNICATION Near-magnetic-axis Geometry of a Closely Quasi-Isodynamic Stellarator M.I. Mikhailov a, J. Nührenberg b, R. Zille b a Russian Research Centre Kurchatov Institute, Moscow,Russia b Max-Planck-Institut

### Figure 5.16 Compound pendulum: (a) At rest in equilibrium, (b) General position with coordinate θ, Freebody

Lecture 27. THE COMPOUND PENDULUM Figure 5.16 Compound pendulum: (a) At rest in equilibrium, (b) General position with coordinate θ, Freebody diagram The term compound is used to distinguish the present

### RESISTIVE BALLOONING MODES AND THE SECOND REGION OF STABILITY

Plasma Physics and Controlled Fusion, Vol. 29, No. 6, pp. 719 to 121, 1987 Printed in Great Britain 0741-3335/87\$3.00+.OO 1OP Publishing Ltd. and Pergamon Journals Ltd. RESISTIVE BALLOONING MODES AND THE

### Concept Question: Normal Force

Concept Question: Normal Force Consider a person standing in an elevator that is accelerating upward. The upward normal force N exerted by the elevator floor on the person is 1. larger than 2. identical

### Observation of tearing mode deceleration and locking due to eddy currents induced in a conducting shell

PHYSICS OF PLASMAS VOLUME 11, NUMBER 5 MAY 2004 Observation of tearing mode deceleration and locking due to eddy currents induced in a conducting shell B. E. Chapman Department of Physics, University of

### Selected Topics in Plasma Astrophysics

Selected Topics in Plasma Astrophysics Range of Astrophysical Plasmas and Relevant Techniques Stellar Winds (Lecture I) Thermal, Radiation, and Magneto-Rotational Driven Winds Connections to Other Areas

### Direct 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

### Linjin Zheng, Infernal Modes at Tokamak H- mode Pedestal A Physics Interpreta;on for Edge Harmonic Oscilla;on (EHO)

International Sherwood Fusion Theory Conference, Austin, May 2-4, 2011 Infernal Modes at Tokamak H- mode Pedestal A Physics Interpreta;on for Edge Harmonic Oscilla;on (EHO) Linjin Zheng, M. T. Kotschenreuther,

### Chapter IX: Nuclear fusion

Chapter IX: Nuclear fusion 1 Summary 1. General remarks 2. Basic processes 3. Characteristics of fusion 4. Solar fusion 5. Controlled fusion 2 General remarks (1) Maximum of binding energy per nucleon

### Müeller Matrix for Imperfect, Rotated Mirrors

From: Steve Scott, Howard Yuh To: MSE Enthusiasts Re: MSE Memo #0c: Proper Treatment of Effect of Mirrors Date: February 1, 004 (in progress) Memo 0c is an extension of memo 0b: it does a calculation for

### The Magnetorotational Instability

The Magnetorotational Instability Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics March 10, 2014 These slides are based off of Balbus & Hawley (1991), Hawley

### Neoclassical Tearing Modes

Neoclassical Tearing Modes O. Sauter 1, H. Zohm 2 1 CRPP-EPFL, Lausanne, Switzerland 2 Max-Planck-Institut für Plasmaphysik, Garching, Germany Physics of ITER DPG Advanced Physics School 22-26 Sept, 2014,

### Physics 20 Homework 2 SIMS 2016

Physics 20 Homework 2 SIMS 2016 Due: Saturday, August 20 th 1. In class, we ignored air resistance in our discussion of projectile motion. Now, let s derive the relevant equation of motion in the case

### Curvilinear Coordinates

University of Alabama Department of Physics and Astronomy PH 106-4 / LeClair Fall 2008 Curvilinear Coordinates Note that we use the convention that the cartesian unit vectors are ˆx, ŷ, and ẑ, rather than

### 2/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

### Edwin Soeryadjaya Problem Theoretical 3: Mirage

The refractive index of the air varies with temperature. Cold air is denser than warm air and has therefore a greater refractive index. Thus a temperature gradient in the atmosphere is always associated

### Saturated ideal modes in advanced tokamak regimes in MAST

Saturated ideal modes in advanced tokamak regimes in MAST IT Chapman 1, M-D Hua 1,2, SD Pinches 1, RJ Akers 1, AR Field 1, JP Graves 3, RJ Hastie 1, CA Michael 1 and the MAST Team 1 EURATOM/CCFE Fusion

### Massachusetts Institute of Technology 22.68J/2.64J Superconducting Magnets. February 27, Lecture #4 Magnetic Forces and Stresses

Massachusetts Institute of Technology.68J/.64J Superconducting Magnets February 7, 003 Lecture #4 Magnetic Forces and Stresses 1 Forces For a solenoid, energy stored in the magnetic field acts equivalent

### VIII. Phase Transformations. Lecture 38: Nucleation and Spinodal Decomposition

VIII. Phase Transformations Lecture 38: Nucleation and Spinodal Decomposition MIT Student In this lecture we will study the onset of phase transformation for phases that differ only in their equilibrium

### Vertical Displacement Events in Shaped Tokamaks. Abstract

Vertical Displacement Events in Shaped Tokamaks A. Y. Aydemir Institute for Fusion Studies The University of Texas at Austin Austin, Texas 78712 USA Abstract Computational studies of vertical displacement

### Internal Transport Barrier Triggering by Rational Magnetic Flux Surfaces in Tokamaks

EFDA JET CP(0)07/09 E. Joffrin, C.D. Challis, G.D. Conway, X. Garbet, A. Gude, S. Guenther, N. C. Hawkes, T.C. Hender, D. Howell, G.T.A. Huysmans, E. Lazarro, P. Maget, M. Marachek, A.G. Peeters, S.D.

### Ph.D. QUALIFYING EXAMINATION PART A. Tuesday, January 3, 2012, 1:00 5:00 P.M.

PhD QUALIFYING EXAMINATION PART A Tuesday, January 3, 212, 1: 5: PM Work each problem on a separate sheet(s) of paper and put your identifying number on each page Do not use your name Each problem has

### D. To the right (Total 1 mark)

1. An electron passes the north pole of a bar magnet as shown below. What is the direction of the magnetic force on the electron? A. Into the page B. Out of the page C. To the left D. To the right 2. A

### Magnetohydrodynamic 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.

### Summary: Curvilinear Coordinates

Physics 2460 Electricity and Magnetism I, Fall 2007, Lecture 10 1 Summary: Curvilinear Coordinates 1. Summary of Integral Theorems 2. Generalized Coordinates 3. Cartesian Coordinates: Surfaces of Constant

### Single 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

### Chapter 4. MHD Equilibrium and Stability. 4.1 Basic Two-Dimensional Equilibrium Equations and Properties. Resistive Diffusion

Chapter 4 MHD Equilibrium and Stability Resistive Diffusion Before discussing equilibrium properties let us first consider effects of electric resistivity. Using the resistive form of Ohm s law with constant

### CONTRIBUTIONS FROM THE NEUTRAL BEAMLINE IRON TO PLASMA NON-AXISYMMETRIC FIELDS

GA A24838 CONTRIBUTIONS FROM THE NEUTRAL BEAMLINE IRON TO PLASMA NON-AXISYMMETRIC FIELDS by J.L. LUXON MARCH 2005 QTYUIOP DISCLAIMER This report was prepared as an account of work sponsored by an agency

### Quiz 4 (Discussion Session) Phys 1302W.400 Spring 2018

Quiz 4 (Discussion ession) Phys 1302W.400 pring 2018 This group quiz consists of one problem that, together with the individual problems on Friday, will determine your grade for quiz 4. For the group problem,

### Electromagnetism Answers to Problem Set 8 Spring Jackson Prob. 4.1: Multipole expansion for various charge distributions

Electromagnetism 76 Answers to Problem Set 8 Spring 6. Jackson Prob. 4.: Multipole expansion for various charge distributions (a) In the first case, we have 4 charges in the xy plane at distance a from

### Thermonuclear Fusion with the Sheared Flow Stabilized Z-Pinch. F. Winterberg and L. F. Wanex

Thermonuclear Fusion with the Sheared Flow Stabilized Z-Pinch F. Winterberg and L. F. Wanex University of Nevada Reno, Department of Physics / MS 37, Reno, Nevada, 89557, USA ABSTRACT Two basic approaches

### Impact of Energetic-Ion-Driven Global Modes on Toroidal Plasma Confinements

Impact of Energetic-Ion-Driven Global Modes on Toroidal Plasma Confinements Kazuo TOI CHS & LHD Experimental Group National Institute for Fusion Science Toki 59-5292, Japan Special contributions from:

### Accelerator Physics Homework #7 P470 (Problems: 1-4)

Accelerator Physics Homework #7 P470 (Problems: -4) This exercise derives the linear transfer matrix for a skew quadrupole, where the magnetic field is B z = B 0 a z, B x = B 0 a x, B s = 0; with B 0 a

### ME-B41 Lab 1: Hydrostatics. Experimental Procedures

ME-B41 Lab 1: Hydrostatics In this lab you will do four brief experiments related to the following topics: manometry, buoyancy, forces on submerged planes, and hydraulics (a hydraulic jack). Each experiment

### PHYS 110B - HW #5 Fall 2005, Solutions by David Pace Equations referenced equations are from Griffiths Problem statements are paraphrased

PHYS 0B - HW #5 Fall 005, Solutions by David Pace Equations referenced equations are from Griffiths Problem statements are paraphrased [.] Imagine a prism made of lucite (n.5) whose cross-section is a

### Energy build- up & release erup3ons. MHD instabili3es Lecture 2 Jan. 23, 2017

Energy build- up & release erup3ons MHD instabili3es Lecture 2 Jan. 23, 2017 Our mission: explain what has happened here Our tool: magnetohydrodynamics (MHD) MHD describes Plasma as a single fluid (combining

### Candidacy Exam Department of Physics February 6, 2010 Part I

Candidacy Exam Department of Physics February 6, 2010 Part I Instructions: ˆ The following problems are intended to probe your understanding of basic physical principles. When answering each question,

### (b) Show that electric field and magnetic field (in Heavyside Lorentz) have units (force)/area or energy/volume.

Problem. Units (a) The Gaussian unit system (or cgs) is often used in Quantum mechanics. Look it up on the internet, and write down the Coulomb and Biot Savart laws in this system of units. (i) What is

### School and Conference on Analytical and Computational Astrophysics November, Angular momentum transport in accretion disks

2292-13 School and Conference on Analytical and Computational Astrophysics 14-25 November, 2011 Angular momentum transport in accretion disks Gianluigi Bodo Osservatorio Astronomico, Torino Italy Angular

### Solutions to Merav Opher (2010) Problems

Solutions to Merav Opher 00 Problems. The normal of the shock is Since from the plot you can obtain all the three components of Bu and Bd, the normal can be easily found. The shock speed is: With the shock

### Global Gyrokinetic Simulation of Tokamak Transport

Global Gyrokinetic Simulation of Tokamak Transport G. Furnish, W. Horton, Y. Kishimoto, M. LeBrun, and T. Tajima Institute for Fusion Studies, The University of Texas at Austin Austin, Texas 78712 and

### Chapter 28 Sources of Magnetic Field

Chapter 28 Sources of Magnetic Field In this chapter we investigate the sources of magnetic of magnetic field, in particular, the magnetic field produced by moving charges (i.e., currents). Ampere s Law

### 1 Introduction. PACS: Kq, y, Kj, a, Cv, Fv DOI: / /14/11/07

Plasma Science and Technology, Vol.14, No.11, Nov. 2012 Three-Dimensional Simulation of Plasma Deformation During Contact Opening in a Circuit Breaker, Including the Analysis of Kink Instability and Sausage

### UNIVERSITY OF MISSOURI-COLUMBIA PHYSICS DEPARTMENT. PART I Qualifying Examination. August 20, 2013, 5:00 p.m. to 8:00 p.m.

UNIVERSITY OF MISSOURI-COLUMBIA PHYSICS DEPARTMENT PART I Qualifying Examination August 20, 2013, 5:00 p.m. to 8:00 p.m. Instructions: The only material you are allowed in the examination room is a writing

### Physics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I

Physics 342 Lecture 23 Radial Separation Lecture 23 Physics 342 Quantum Mechanics I Friday, March 26th, 2010 We begin our spherical solutions with the simplest possible case zero potential. Aside from

### GA A24166 SUPER-INTENSE QUASI-NEUTRAL PROTON BEAMS INTERACTING WITH PLASMA: A NUMERICAL INVESTIGATION

GA A24166 SUPER-INTENSE QUASI-NEUTRAL PROTON BEAMS INTERACTING WITH PLASMA: A NUMERICAL INVESTIGATION by H. RUHL, T.E. COWAN, and R.B. STEPHENS OCTOBER 2 DISCLAIMER This report was prepared as an account

### Analyses on the Ionization Instability of Non-Equilibrium Seeded Plasma in an MHD Generator

Plasma Science and Technology, Vol.18, No.6, Jun. 2016 Analyses on the Ionization Instability of Non-Equilibrium Seeded Plasma in an MHD Generator Le Chi KIEN Ho Chi Minh City University of Technology

### Magnetic Deflection of Ionized Target Ions

Magnetic Deflection of Ionized Target Ions D. V. Rose, A. E. Robson, J. D. Sethian, D. R. Welch, and R. E. Clark March 3, 005 HAPL Meeting, NRL Solid wall, magnetic deflection 1. Cusp magnetic field imposed

### Rotation and Neoclassical Ripple Transport in ITER

Rotation and Neoclassical Ripple Transport in ITER Elizabeth J. Paul 1 Matt Landreman 1 Francesca Poli 2 Don Spong 3 Håkan Smith 4 William Dorland 1 1 University of Maryland 2 Princeton Plasma Physics

### Design of the PROTO-SPHERA experiment and of its first step (MULTI-PINCH)

INSTITUTE OF PHYSICS PUBLISHING and INTERNATIONAL ATOMIC ENERGY AGENCY NUCLEAR FUSION Nucl. Fusion 46 (2006) S613 S624 Design of the PROTO-SPHERA experiment and of its first step (MULTI-PINCH) doi:10.1088/0029-5515/46/8/s07

### UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 110A. Homework #7. Benjamin Stahl. March 3, 2015

UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS A Homework #7 Benjamin Stahl March 3, 5 GRIFFITHS, 5.34 It will be shown that the magnetic field of a dipole can written in the following