Magnetic Control of Perturbed Plasma Equilibria

Size: px
Start display at page:

Download "Magnetic Control of Perturbed Plasma Equilibria"

Transcription

1 Magnetic Control of Perturbed Plasma Equilibria Nikolaus Rath February 17th, 2012 N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

2 The HBT-EP Tokamak S. Angelini, J. Bialek, P. Byrne, B. DeBono, P. Hughes, J. Levesque, B. Li, M. Mauel, G. Navratil, Q. Peng, D. Rhodes, D. Shiraki, C. Stoafer N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

3 Recent upgrade allows high-resolution magnetic measurement and control Inside the HBT-EP vacuum chamber N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

4 Plasma dynamics are complicated f ( α t + v α f α e E + 1 c v B) f α p = 0 B = 4π j c + 1 E c t, E = 1 c E = 4πρ, B = 0 ρ = e (f i f e ) d 3 p j = e (f i f e ) v d 3 p p/m α v α = 1 +p 2 /m 2 αc 2 B t N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

5 MHD makes things easier With a few not-quite-appropriate approximations: ρ d u dt = ( B) B p d B dt = ( u B) d dt dρ = (ρ u) dt ) (1) = 0 (2) ( p ρ γ But even better, in a Tokamak: u B µ0 ρ meters (3) microseconds d u dt 0 (4) N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

6 MHD makes things easier With a few not-quite-appropriate approximations: ρ d u dt = ( B) B p d B dt = ( u B) d dt dρ = (ρ u) dt ) (1) = 0 (2) ( p ρ γ But even better, in a Tokamak: u B µ0 ρ meters (3) microseconds d u dt 0 (4) N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

7 MHD makes things easier With a few not-quite-appropriate approximations: ρ d u dt = ( B) B p d B dt = ( u B) d dt dρ = (ρ u) dt ) (1) = 0 (2) ( p ρ γ But even better, in a Tokamak: u B µ0 ρ meters (3) microseconds d u dt 0 (4) N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

8 Tokamak plasmas transition through MHD equilibria Plasma reaches force-balance within microseconds Longer lived plasmas must thus be in equilibrium Evolution is caused by all our approximations The plasma always has minimum MHD energy, but the properties of the minimum-energy-state are changing slowly N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

9 Tokamak equilibria are tractable This talk is about: ( B) B = p (5) N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

10 In Tokamaks, the plasma state is axisymmetric with a 3d perturbation Solving ( B) B = p in axisymmetry is easy so we assume someone did it. Solving ( B) B = p in 3d is hard, so we don t want to do it. So will solve ( B) B = p with just small 3d perturbation - Tokamaks are great! A small 3d displacement ζ must minimize (Boozer, 1999 & 2003) dw( ζ) = ζ F( ζ)dv (6) plasma F is known from the axisymmetric solution, minimization is just number crunching. Just need boundary conditions on plasma surface N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

11 In Tokamaks, the plasma state is axisymmetric with a 3d perturbation Solving ( B) B = p in axisymmetry is easy so we assume someone did it. Solving ( B) B = p in 3d is hard, so we don t want to do it. So will solve ( B) B = p with just small 3d perturbation - Tokamaks are great! A small 3d displacement ζ must minimize (Boozer, 1999 & 2003) dw( ζ) = ζ F( ζ)dv (6) plasma F is known from the axisymmetric solution, minimization is just number crunching. Just need boundary conditions on plasma surface N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

12 In Tokamaks, the plasma state is axisymmetric with a 3d perturbation Solving ( B) B = p in axisymmetry is easy so we assume someone did it. Solving ( B) B = p in 3d is hard, so we don t want to do it. So will solve ( B) B = p with just small 3d perturbation - Tokamaks are great! A small 3d displacement ζ must minimize (Boozer, 1999 & 2003) dw( ζ) = ζ F( ζ)dv (6) plasma F is known from the axisymmetric solution, minimization is just number crunching. Just need boundary conditions on plasma surface N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

13 In Tokamaks, the plasma state is axisymmetric with a 3d perturbation Solving ( B) B = p in axisymmetry is easy so we assume someone did it. Solving ( B) B = p in 3d is hard, so we don t want to do it. So will solve ( B) B = p with just small 3d perturbation - Tokamaks are great! A small 3d displacement ζ must minimize (Boozer, 1999 & 2003) dw( ζ) = ζ F( ζ)dv (6) plasma F is known from the axisymmetric solution, minimization is just number crunching. Just need boundary conditions on plasma surface N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

14 The permeability matrix defines the plasma response to external fields The plasma shape follows magnetic field lines The non-axisymmetric field on the plasma surface is the boundary conditions for the plasma displacement. The field has components from both the plasma and the environment, write it as b( ζ) = b plasma +b external (7) For easier math, write surface functions as vectors of expansion coefficients, e.g. b(θ,φ) = i Φ i f i (θ,φ) = Φ f(φ,θ) (8) Define permeability matrix P by: P. Φ external = Φ (9) N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

15 The permeability matrix defines the plasma response to external fields The plasma shape follows magnetic field lines The non-axisymmetric field on the plasma surface is the boundary conditions for the plasma displacement. The field has components from both the plasma and the environment, write it as b( ζ) = b plasma +b external (7) For easier math, write surface functions as vectors of expansion coefficients, e.g. b(θ,φ) = i Φ i f i (θ,φ) = Φ f(φ,θ) (8) Define permeability matrix P by: P. Φ external = Φ (9) N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

16 Plasma response example External field N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

17 Plasma response example Total field N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

18 Half-Time Summary We want to control the plasma In Tokamaks, the plasma is mostly in MHD equilibrium In Tokamaks, the equilibrium is also approximately axisymmetric The axisymmetric part is constant, the 3d perturbation changes in time The axisymmetric part is easily measured and calculated, we assume it is known The 3d perturbation is uniquely defined by magnetic field due to external currents on the plasma surface To control Tokamak plasmas, we therefore have to measure, track and control external currents N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

19 Half-Time Summary We want to control the plasma In Tokamaks, the plasma is mostly in MHD equilibrium In Tokamaks, the equilibrium is also approximately axisymmetric The axisymmetric part is constant, the 3d perturbation changes in time The axisymmetric part is easily measured and calculated, we assume it is known The 3d perturbation is uniquely defined by magnetic field due to external currents on the plasma surface To control Tokamak plasmas, we therefore have to measure, track and control external currents N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

20 The system model consists of plasma, wall, control coils, and sensors Control Coils: I c Wall: I w (finite element model) Plasma: I p Sensors: Φ s External Currents I x = ( I c, I x ) N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

21 Maxwell s equations allow to represent all system components as vectors Maxwell: the normal magnetic field on a closed surface uniquely defines the external field in the enclosed volume Plasma sees all external currents as Φ x = M px. I x External circuits see plasma as Φ p = (P 1). Φ x. We can express all interactions as matrices. The complete system obeys R [ ] x. I x + L x +M xp.l 1 p.(p 1).M px. d I x dt = V x (10) M and L matrices encapsulate all geometric information. P depends on the axisymmetric equilibrium and contains information about the plasma. N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

22 Maxwell s equations allow to represent all system components as vectors Maxwell: the normal magnetic field on a closed surface uniquely defines the external field in the enclosed volume Plasma sees all external currents as Φ x = M px. I x External circuits see plasma as Φ p = (P 1). Φ x. We can express all interactions as matrices. The complete system obeys R [ ] x. I x + L x +M xp.l 1 p.(p 1).M px. d I x dt = V x (10) M and L matrices encapsulate all geometric information. P depends on the axisymmetric equilibrium and contains information about the plasma. N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

23 Maxwell s equations allow to represent all system components as vectors Maxwell: the normal magnetic field on a closed surface uniquely defines the external field in the enclosed volume Plasma sees all external currents as Φ x = M px. I x External circuits see plasma as Φ p = (P 1). Φ x. We can express all interactions as matrices. The complete system obeys R [ ] x. I x + L x +M xp.l 1 p.(p 1).M px. d I x dt = V x (10) M and L matrices encapsulate all geometric information. P depends on the axisymmetric equilibrium and contains information about the plasma. N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

24 Maxwell s equations allow to represent all system components as vectors Maxwell: the normal magnetic field on a closed surface uniquely defines the external field in the enclosed volume Plasma sees all external currents as Φ x = M px. I x External circuits see plasma as Φ p = (P 1). Φ x. We can express all interactions as matrices. The complete system obeys R [ ] x. I x + L x +M xp.l 1 p.(p 1).M px. d I x dt = V x (10) M and L matrices encapsulate all geometric information. P depends on the axisymmetric equilibrium and contains information about the plasma. N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

25 Controlling the plasma To control the solution of we need to control a system of the form ( B) B = p (11) Required steps: d I x dt = A. I x +B. V x (12) N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

26 Controlling the plasma To control the solution of we need to control a system of the form Required steps: ( B) B = p (11) d I x dt = A. I x +B. V x (12) 1 Subtract axisymmetric fields from measurements 2 Use measurements to estimate system state 3 Compare current plasma shape with target shape 4 Find control coil voltages that bring shape closer to target N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

27 Controlling the plasma To control the solution of we need to control a system of the form Required steps: ( B) B = p (11) d I x dt = A. I x +B. V x (12) 1 Subtract axisymmetric fields from measurements 2 Use measurements to estimate system state 3 Compare current plasma shape with target shape 4 Find control coil voltages that bring shape closer to target (not discussed in this talk) N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

28 Equilibrium fields are subtracted using continuous polynomial prediction Measurement errors in the (large) axisymmetric contribution easily exceed non-axisymmetric fluctuations. Solution: Continuously fit last n measurements of every sensor to quadratic polynomial, subtract predicted amplitude. N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

29 SVD of measurement matrix determines relative visibility of states Sensors measure magnetic flux: Φ s = M sx. I x +M sp. I p [ = M sx +M sp.l 1 p.(p 1).M px ]. I x =: C. (13) I x Singular value decomposition of C gives 1 v 1 U.S.V = u 1 u 2. s s 2. v 2 (14).... v i is the basis of directly measurable states s i is the relative visibility of a state We could drop invisible states and invert, but there is a better solution N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

30 The system model is used to estimate the state Control system continuously calculates: Error in calculated I x obeys d I x dt = A. I x +B. V x +K. ( Φ s C. ) I x (15) d e dt = (AT C T.K T ). e (16) Could pick K to obtain specific eigenvalues but which exactly? Better choice: let R be measurement uncertainties and S mode visibilities, and choose K to minimize S. e(t) 2 + R.K T. e 2 dt (17) 0 More visible states will converge faster, more reliable sensors will be used preferentially N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

31 Experimental results: wall model and shot reproducibility needs to be improved HBT-EP plasmas are short-lived, so equilibrium needs to be computed offline Variations in major radius change equilibrium considerably between shots, preventing pre-computation Plan: address by extending control to plasma major radius Vacuum testing show reasonably agreement, better wall model may give further improvements. Work in progress to add vacuum vessel. N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

32 Summary 1 The evolution of Tokamak plasmas can be approximated as a transition through a sequence of MHD equilibria. 2 These equilibria are approximately axisymmetric and consist of a static, axisymmetric part and a changing, 3-dimensional perturbation. 3 Transitions caused by changing external fields can be modeled as a linear, time invariant system. This system can be controlled with standard techniques from control theory, using magnetic sensors and control coils. 4 Experimental tests show reasonable agreement for vacuum model. Plasma tests were impaired by low shot reproducibility, which will be addressed by control of plasma major radius. N. Rath (Columbia University) Magnetic Control of Perturbed Plasma Equilibria February 17th, / 19

Effects of Noise in Time Dependent RWM Feedback Simulations

Effects of Noise in Time Dependent RWM Feedback Simulations Effects of Noise in Time Dependent RWM Feedback Simulations O. Katsuro-Hopkins, J. Bialek, G. Navratil (Department of Applied Physics and Applied Mathematics, Columbia University, New York, NY USA) Building

More information

The Effects of Noise and Time Delay on RWM Feedback System Performance

The Effects of Noise and Time Delay on RWM Feedback System Performance The Effects of Noise and Time Delay on RWM Feedback System Performance O. Katsuro-Hopkins, J. Bialek, G. Navratil (Department of Applied Physics and Applied Mathematics, Columbia University, New York,

More information

RESISTIVE WALL MODE STABILIZATION RESEARCH ON DIII D STATUS AND RECENT RESULTS

RESISTIVE WALL MODE STABILIZATION RESEARCH ON DIII D STATUS AND RECENT RESULTS RESISTIVE WALL MODE STABILIZATION RESEARCH ON STATUS AND RECENT RESULTS by A.M. Garofalo1 in collaboration with J. Bialek,1 M.S. Chance,2 M.S. Chu,3 T.H. Jensen,3 L.C. Johnson,2 R.J. La Haye,3 G.A. Navratil,1

More information

Dynamical plasma response of resistive wall modes to changing external magnetic perturbations

Dynamical plasma response of resistive wall modes to changing external magnetic perturbations Dynamical plasma response of resistive wall modes to changing external magnetic perturbations M. Shilov, C. Cates, R. James, A. Klein, O. Katsuro-Hopkins, Y. Liu, M. E. Mauel, D. A. Maurer, G. A. Navratil,

More information

Dynamical plasma response of resistive wall modes to changing external magnetic perturbations a

Dynamical plasma response of resistive wall modes to changing external magnetic perturbations a PHYSICS OF PLASMAS VOLUME 11, NUMBER 5 MAY 2004 Dynamical plasma response of resistive wall modes to changing external magnetic perturbations a M. Shilov, b) C. Cates, R. James, A. Klein, O. Katsuro-Hopkins,

More information

35. RESISTIVE INSTABILITIES: CLOSING REMARKS

35. RESISTIVE INSTABILITIES: CLOSING REMARKS 35. RESISTIVE INSTABILITIES: CLOSING REMARKS In Section 34 we presented a heuristic discussion of the tearing mode. The tearing mode is centered about x = 0, where F x that B 0 ( ) = k! B = 0. Note that

More information

Plasma Response Control Using Advanced Feedback Techniques

Plasma Response Control Using Advanced Feedback Techniques Plasma Response Control Using Advanced Feedback Techniques by M. Clement 1 with J. M. Hanson 1, J. Bialek 1 and G. A. Navratil 1 1 Columbia University Presented at 59 th Annual APS Meeting Division of

More information

RWM FEEDBACK STABILIZATION IN DIII D: EXPERIMENT-THEORY COMPARISONS AND IMPLICATIONS FOR ITER

RWM FEEDBACK STABILIZATION IN DIII D: EXPERIMENT-THEORY COMPARISONS AND IMPLICATIONS FOR ITER GA A24759 RWM FEEDBACK STABILIZATION IN DIII D: EXPERIMENT-THEORY COMPARISONS AND IMPLICATIONS FOR ITER by A.M. GAROFALO, J. BIALEK, M.S. CHANCE, M.S. CHU, D.H. EDGELL, G.L. JACKSON, T.H. JENSEN, R.J.

More information

Derivation of dynamo current drive in a closed current volume and stable current sustainment in the HIT SI experiment

Derivation of dynamo current drive in a closed current volume and stable current sustainment in the HIT SI experiment Derivation of dynamo current drive and stable current sustainment in the HIT SI experiment 1 Derivation of dynamo current drive in a closed current volume and stable current sustainment in the HIT SI experiment

More information

0 Magnetically Confined Plasma

0 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 information

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

Lecture # 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 information

RWM Control in FIRE and ITER

RWM Control in FIRE and ITER RWM Control in FIRE and ITER Gerald A. Navratil with Jim Bialek, Allen Boozer & Oksana Katsuro-Hopkins MHD Mode Control Workshop University of Texas-Austin 3-5 November, 2003 OUTLINE REVIEW OF VALEN MODEL

More information

Plasma Current, Position and Shape Control June 2, 2010

Plasma Current, Position and Shape Control June 2, 2010 , Position and Shape Control June 2, 2010 June 2, 2010 - ITER International Summer School 2010 Control 1 in collaboration with CREATE and EFDA-JET PPCC contributors 1 CREATE, Università di Napoli Federico

More information

Resistive Wall Mode Control in DIII-D

Resistive Wall Mode Control in DIII-D Resistive Wall Mode Control in DIII-D by Andrea M. Garofalo 1 for G.L. Jackson 2, R.J. La Haye 2, M. Okabayashi 3, H. Reimerdes 1, E.J. Strait 2, R.J. Groebner 2, Y. In 4, M.J. Lanctot 1, G.A. Navratil

More information

MHD Linear Stability Analysis Using a Full Wave Code

MHD 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 information

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

More information

A MIMO architecture for integrated control of plasma shape and flux expansion for the EAST tokamak

A MIMO architecture for integrated control of plasma shape and flux expansion for the EAST tokamak A MIMO architecture for integrated control of plasma shape and flux expansion for the EAST tokamak R. Albanese 1 R. Ambrosino 2 G. Calabrò 3 A. Castaldo 1 F. Crisanti 3 G. De Tommasi 1 L. Liu 4 Z. P. Luo

More information

Model based optimization and estimation of the field map during the breakdown phase in the ITER tokamak

Model based optimization and estimation of the field map during the breakdown phase in the ITER tokamak Model based optimization and estimation of the field map during the breakdown phase in the ITER tokamak Roberto Ambrosino 1 Gianmaria De Tommasi 2 Massimiliano Mattei 3 Alfredo Pironti 2 1 CREATE, Università

More information

Reconstruction of Pressure Profile Evolution during Levitated Dipole Experiments

Reconstruction of Pressure Profile Evolution during Levitated Dipole Experiments Reconstruction of Pressure Profile Evolution during Levitated Dipole Experiments M. E. Mauel, A. Boxer, J. Ellsworth, D. Garnier, J. Kesner ICC Conference: Reno, Nevada (June 24, 28) 1 Abstract Magnetic

More information

DIII D. by M. Okabayashi. Presented at 20th IAEA Fusion Energy Conference Vilamoura, Portugal November 1st - 6th, 2004.

DIII D. by M. Okabayashi. Presented at 20th IAEA Fusion Energy Conference Vilamoura, Portugal November 1st - 6th, 2004. Control of the Resistive Wall Mode with Internal Coils in the Tokamak (EX/3-1Ra) Active Measurement of Resistive Wall Mode Stability in Rotating High Beta Plasmas (EX/3-1Rb) by M. Okabayashi Presented

More information

Ideal MHD Equilibria

Ideal 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 information

Modeling of active control of external magnetohydrodynamic instabilities*

Modeling of active control of external magnetohydrodynamic instabilities* PHYSICS OF PLASMAS VOLUME 8, NUMBER 5 MAY 2001 Modeling of active control of external magnetohydrodynamic instabilities* James Bialek, Allen H. Boozer, M. E. Mauel, and G. A. Navratil Department of Applied

More information

AN UPDATE ON DIVERTOR HEAT LOAD ANALYSIS

AN UPDATE ON DIVERTOR HEAT LOAD ANALYSIS AN UPDATE ON DIVERTOR HEAT LOAD ANALYSIS T.K. Mau, A. Grossman, A.R. Raffray UC-San Diego H. McGuinness RPI ARIES-CS Project Meeting June 4-5, 5 University of Wisconsin, Madison OUTLINE Divertor design

More information

Reconstruction of the Pressure Profile of LDX High Beta Plasma

Reconstruction of the Pressure Profile of LDX High Beta Plasma Reconstruction of the Pressure Profile of LDX High Beta Plasma M. Mauel, I. Karim, A. Boxer, J. Ellsworth, D. Garnier, A. Hansen, J. Kesner, and E. Ortiz Dept. of Applied Physics and Applied Mathematics,

More information

Design of next step tokamak: Consistent analysis of plasma flux consumption and poloidal field system

Design of next step tokamak: Consistent analysis of plasma flux consumption and poloidal field system Design of next step tokamak: Consistent analysis of plasma flux consumption and poloidal field system J.M. Ané 1, V. Grandgirard, F. Albajar 1, J.Johner 1 1Euratom-CEA Association, Cadarache, France Euratom-EPFL

More information

KSTAR Equilibrium Operating Space and Projected Stabilization at High Normalized Beta

KSTAR Equilibrium Operating Space and Projected Stabilization at High Normalized Beta 1 THS/P2-05 KSTAR Equilibrium Operating Space and Projected Stabilization at High Normalized Beta Y.S. Park 1), S.A. Sabbagh 1), J.W. Berkery 1), J.M. Bialek 1), Y.M. Jeon 2), S.H. Hahn 2), N. Eidietis

More information

Going with the flow: A study of Lagrangian derivatives

Going with the flow: A study of Lagrangian derivatives 1 Going with the flow: A study of Lagrangian derivatives Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc/ 12 February

More information

Tokamak Divertor System Concept and the Design for ITER. Chris Stoafer April 14, 2011

Tokamak Divertor System Concept and the Design for ITER. Chris Stoafer April 14, 2011 Tokamak Divertor System Concept and the Design for ITER Chris Stoafer April 14, 2011 Presentation Overview Divertor concept and purpose Divertor physics General design considerations Overview of ITER divertor

More information

where the last equality follows from the divergence theorem. Now since we did this for an arbitrary volume τ, it must hold locally:

where the last equality follows from the divergence theorem. Now since we did this for an arbitrary volume τ, it must hold locally: 8 Electrodynamics Read: Boas Ch. 6, particularly sec. 10 and 11. 8.1 Maxwell equations Some of you may have seen Maxwell s equations on t-shirts or encountered them briefly in electromagnetism courses.

More information

Effect of an error field on the stability of the resistive wall mode

Effect of an error field on the stability of the resistive wall mode PHYSICS OF PLASMAS 14, 022505 2007 Effect of an error field on the stability of the resistive wall mode Richard Fitzpatrick Institute for Fusion Studies, Department of Physics, University of Texas at Austin,

More information

Formation 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 ) 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

The Virial Theorem, MHD Equilibria, and Force-Free Fields

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

More information

- Effect of Stochastic Field and Resonant Magnetic Perturbation on Global MHD Fluctuation -

- Effect of Stochastic Field and Resonant Magnetic Perturbation on Global MHD Fluctuation - 15TH WORKSHOP ON MHD STABILITY CONTROL: "US-Japan Workshop on 3D Magnetic Field Effects in MHD Control" U. Wisconsin, Madison, Nov 15-17, 17, 2010 LHD experiments relevant to Tokamak MHD control - Effect

More information

Plasma models for the design of the ITER PCS

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

More information

Requirements for Active Resistive Wall Mode (RWM) Feedback Control

Requirements for Active Resistive Wall Mode (RWM) Feedback Control Requirements for Active Resistive Wall Mode (RWM) Feedback Control Yongkyoon In 1 In collaboration with M.S. Chu 2, G.L. Jackson 2, J.S. Kim 1, R.J. La Haye 2, Y.Q. Liu 3, L. Marrelli 4, M. Okabayashi

More information

Three Dimensional Effects in Tokamaks How Tokamaks Can Benefit From Stellarator Research

Three Dimensional Effects in Tokamaks How Tokamaks Can Benefit From Stellarator Research 1 TH/P9-10 Three Dimensional Effects in Tokamaks How Tokamaks Can Benefit From Stellarator Research S. Günter, M. Garcia-Munoz, K. Lackner, Ph. Lauber, P. Merkel, M. Sempf, E. Strumberger, D. Tekle and

More information

Computations with Discontinuous Basis Functions

Computations with Discontinuous Basis Functions Computations with Discontinuous Basis Functions Carl Sovinec University of Wisconsin-Madison NIMROD Team Meeting November 12, 2011 Salt Lake City, Utah Motivation The objective of this work is to make

More information

Resistive Wall Mode Observation and Control in ITER-Relevant Plasmas

Resistive Wall Mode Observation and Control in ITER-Relevant Plasmas Resistive Wall Mode Observation and Control in ITER-Relevant Plasmas J. P. Levesque April 12, 2011 1 Outline Basic Resistive Wall Mode (RWM) model RWM stability, neglecting kinetic effects Sufficient for

More information

Before you begin read these instructions carefully.

Before you begin read these instructions carefully. MATHEMATICAL TRIPOS Part IB Thursday, 9 June, 2011 9:00 am to 12:00 pm PAPER 3 Before you begin read these instructions carefully. Each question in Section II carries twice the number of marks of each

More information

Transmission Lines and E. M. Waves Prof. R. K. Shevgaonkar Department of Electrical Engineering Indian Institute of Technology, Bombay

Transmission Lines and E. M. Waves Prof. R. K. Shevgaonkar Department of Electrical Engineering Indian Institute of Technology, Bombay Transmission Lines and E. M. Waves Prof. R. K. Shevgaonkar Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture 18 Basic Laws of Electromagnetics We saw in the earlier lecture

More information

Model based estimation of the Eddy currents for the ITER tokamak *

Model based estimation of the Eddy currents for the ITER tokamak * Model based estimation of the Eddy currents for the ITER tokamak * Alfredo Pironti 1 Consorzio CREATE - University of Naples Federico II Via Claudio 21, 80125 Napoli, Italy E-mail: pironti@unina.it Roberto

More information

Plasma Current, Position and Shape Control Hands-on Session June 2, 2010

Plasma Current, Position and Shape Control Hands-on Session June 2, 2010 , Position and Shape Control Hands-on Session June 2, 2010 June 2, 2010 - ITER International Summer School 2010 Shape 1 1 CREATE, Università di Napoli Federico II 1 for the JET tokamak PF Currents for

More information

Imposed Dynamo Current Drive

Imposed Dynamo Current Drive Imposed Dynamo Current Drive by T. R. Jarboe, C. Akcay, C. J. Hansen, A. C. Hossack, G. J. Marklin, K. Morgan, B. A. Nelson, D. A. Sutherland, B. S. Victor University of Washington, Seattle, WA 98195,

More information

Bunno, M.; Nakamura, Y.; Suzuki, Y. Matsunaga, G.; Tani, K. Citation Plasma Science and Technology (

Bunno, M.; Nakamura, Y.; Suzuki, Y. Matsunaga, G.; Tani, K. Citation Plasma Science and Technology ( Title The Finite Beta Effects on the Toro Tokamak Plasma Author(s) Bunno, M.; Nakamura, Y.; Suzuki, Y. Matsunaga, G.; Tani, K. Citation Plasma Science and Technology (2013 Issue Date 2013-02 URL http://hdl.handle.net/2433/173038

More information

Simple examples of MHD equilibria

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 information

GERALD A. NAVRATIL. Thomas Alva Edison Professor of Applied Physics Department of Applied Physics and Applied Mathematics Columbia University

GERALD A. NAVRATIL. Thomas Alva Edison Professor of Applied Physics Department of Applied Physics and Applied Mathematics Columbia University 1 GERALD A. NAVRATIL Thomas Alva Edison Professor of Applied Physics Department of Applied Physics and Applied Mathematics Columbia University 209 S. W. Mudd Building office: (212) 854-4496 Columbia University

More information

Disruption mitigation in ITER

Disruption mitigation in ITER Disruption mitigation in ITER S.Putvinski, R.Pitts, M.Sugihara, L.Zakharov Page 1 Outline Introduction Massive gas injection for mitigation of thermal loads Forces on VV and in vessel components Suppression

More information

Review problems for MA 54, Fall 2004.

Review problems for MA 54, Fall 2004. Review problems for MA 54, Fall 2004. Below are the review problems for the final. They are mostly homework problems, or very similar. If you are comfortable doing these problems, you should be fine on

More information

Effects of stellarator transform on sawtooth oscillations in CTH. Jeffrey Herfindal

Effects of stellarator transform on sawtooth oscillations in CTH. Jeffrey Herfindal Effects of stellarator transform on sawtooth oscillations in CTH Jeffrey Herfindal D.A. Ennis, J.D. Hanson, G.J. Hartwell, E.C. Howell, C.A. Johnson, S.F. Knowlton, X. Ma, D.A. Maurer, M.D. Pandya, N.A.

More information

Review of similarity transformation and Singular Value Decomposition

Review of similarity transformation and Singular Value Decomposition Review of similarity transformation and Singular Value Decomposition Nasser M Abbasi Applied Mathematics Department, California State University, Fullerton July 8 7 page compiled on June 9, 5 at 9:5pm

More information

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

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

More information

ELECTROMAGNETIC LIQUID METAL WALL PHENOMENA

ELECTROMAGNETIC LIQUID METAL WALL PHENOMENA ELECTROMAGNETIC LIQUID METAL WALL PHENOMENA BY BOB WOOLLEY 15-19 FEBRUARY 1999 APEX-6 MEETING LIQUID WALLS A sufficiently thick, flowing, liquid first wall and tritium breeding blanket which almost completely

More information

From tokamaks to stellarators: understanding the role of 3D shaping

From tokamaks to stellarators: understanding the role of 3D shaping Under consideration for publication in J. Plasma Phys. From tokamaks to stellarators: understanding the role of 3D shaping Samuel A. Lazerson and John C. Schmitt 2 Princeton Plasma Physics Laboratory,

More information

Disruption dynamics in NSTX. long-pulse discharges. Presented by J.E. Menard, PPPL. for the NSTX Research Team

Disruption dynamics in NSTX. long-pulse discharges. Presented by J.E. Menard, PPPL. for the NSTX Research Team Disruption dynamics in NSTX long-pulse discharges Presented by J.E. Menard, PPPL for the NSTX Research Team Workshop on Active Control of MHD Stability: Extension of Performance Monday, November 18, 2002

More information

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

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

More information

Toroidal confinement of non-neutral plasma. Martin Droba

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

More information

MHD 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) 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 information

The Status of the Design and Construction of the Columbia Non-neutral Torus

The Status of the Design and Construction of the Columbia Non-neutral Torus The Status of the Design and Construction of the Columbia Non-neutral Torus J. P. Kremer,T.S.Pedersen,N.Pomphrey,W.Reiersen and F. Dahlgren Dept. of Applied Physics and Applied Mathematics, Columbia University,

More information

Current Drive Experiments in the HIT-II Spherical Tokamak

Current Drive Experiments in the HIT-II Spherical Tokamak Current Drive Experiments in the HIT-II Spherical Tokamak T. R. Jarboe, P. Gu, V. A. Izzo, P. E. Jewell, K. J. McCollam, B. A. Nelson, R. Raman, A. J. Redd, P. E. Sieck, and R. J. Smith, Aerospace & Energetics

More information

Behavior of Compact Toroid Injected into the External Magnetic Field

Behavior of Compact Toroid Injected into the External Magnetic Field Behavior of Compact Toroid Injected into the External Magnetic Field M. Nagata 1), N. Fukumoto 1), H. Ogawa 2), T. Ogawa 2), K. Uehara 2), H. Niimi 3), T. Shibata 2), Y. Suzuki 4), Y. Miura 2), N. Kayukawa

More information

Shear Alfvén Wave Continuous Spectrum in the Presence of a Magnetic Island

Shear Alfvén Wave Continuous Spectrum in the Presence of a Magnetic Island 1 TH/P3-5 Shear Alfvén Wave Continuous Spectrum in the Presence of a Magnetic Island A. Biancalani 1), L. Chen 2) 3), F. Pegoraro 1), F. Zonca 4), S. V. Annibaldi 5), A. Botrugno 4), P. Buratti 4) and

More information

Linear Least-Squares Data Fitting

Linear Least-Squares Data Fitting CHAPTER 6 Linear Least-Squares Data Fitting 61 Introduction Recall that in chapter 3 we were discussing linear systems of equations, written in shorthand in the form Ax = b In chapter 3, we just considered

More information

Nonlinear MHD Stability and Dynamical Accessibility

Nonlinear MHD Stability and Dynamical Accessibility Nonlinear MHD Stability and Dynamical Accessibility Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University Philip J. Morrison Department of Physics and Institute

More information

STABILIZATION OF THE RESISTIVE WALL MODE IN DIII D BY PLASMA ROTATION AND MAGNETIC FEEDBACK

STABILIZATION OF THE RESISTIVE WALL MODE IN DIII D BY PLASMA ROTATION AND MAGNETIC FEEDBACK GA A24014 STABILIZATION OF THE RESISTIVE WALL MODE IN DIII D BY PLASMA ROTATION AND MAGNETIC FEEDBACK by M. Okabayashi, J. Bialek, M.S. Chance, M.S. Chu, E.D. Fredrickson, A.M. Garofalo, R. Hatcher, T.H.

More information

W15D1: Poynting Vector and Energy Flow. Today s Readings: Course Notes: Sections 13.6,

W15D1: Poynting Vector and Energy Flow. Today s Readings: Course Notes: Sections 13.6, W15D1: Poynting Vector and Energy Flow Today s Readings: Course Notes: Sections 13.6, 13.12.3-13.12.4 1 Announcements Final Math Review Week 15 Tues from 9-11 pm in 32-082 Final Exam Monday Morning May

More information

Issues in Neoclassical Tearing Mode Theory

Issues 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 information

Experiment 1 1. Charge- to- Mass Ratio of the Electron Physics 2150 Experiment No. 1 University of Colorado

Experiment 1 1. Charge- to- Mass Ratio of the Electron Physics 2150 Experiment No. 1 University of Colorado Experiment 1 1 Introduction Charge- to- Mass Ratio of the Electron Physics 2150 Experiment No. 1 University of Colorado Both the charge and the mass of the electron are fundamental constants of considerable

More information

On existence of resistive magnetohydrodynamic equilibria

On 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 information

Current-driven instabilities

Current-driven instabilities Current-driven instabilities Ben Dudson Department of Physics, University of York, Heslington, York YO10 5DD, UK 21 st February 2014 Ben Dudson Magnetic Confinement Fusion (1 of 23) Previously In the last

More information

Performance limits. Ben Dudson. 24 th February Department of Physics, University of York, Heslington, York YO10 5DD, UK

Performance 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 information

Control of Sawtooth Oscillation Dynamics using Externally Applied Stellarator Transform. Jeffrey Herfindal

Control of Sawtooth Oscillation Dynamics using Externally Applied Stellarator Transform. Jeffrey Herfindal Control of Sawtooth Oscillation Dynamics using Externally Applied Stellarator Transform Jeffrey Herfindal D.A. Ennis, J.D. Hanson, G.J. Hartwell, S.F. Knowlton, X. Ma, D.A. Maurer, M.D. Pandya, N.A. Roberds,

More information

Effect of non-axisymmetric magnetic perturbations on divertor heat and particle flux profiles

Effect of non-axisymmetric magnetic perturbations on divertor heat and particle flux profiles 1 EXD/P3-01 Effect of non-axisymmetric magnetic perturbations on divertor heat and particle flux profiles J-W. Ahn 1, J.M. Canik 1, R. Maingi 1, T.K. Gray 1, J.D. Lore 1, A.G. McLean 1, J.-K. Park 2, A.L.

More information

Effects of q(r) on the Alpha Particle Ripple Loss in TFTR

Effects of q(r) on the Alpha Particle Ripple Loss in TFTR Effects of q(r) on the Alpha Particle Ripple Loss in TFTR S.J. Zweben, D.S. Darrow, S. Batha, R.V. Budny, M. Diesso, H. W. Herrmann, J. Giarrusso, M. Redi, H. Takahashi, S. von Goeler, R.B. White, and

More information

The Levitated Dipole Experiment: Towards Fusion Without Tritium

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

The compact dipole configuration for plasma confinement

The compact dipole configuration for plasma confinement The compact dipole configuration for plasma confinement D. A. Baver Lodestar Research Corporation, Boulder, Colorado, 80301 March, 2011 Submitted to Journal of Fusion Energy LRC-11-141 Lodestar Research

More information

N H I. 3.2 l When a conducting coil is placed in a magnetic field, the magnetic flux is

N H I. 3.2 l When a conducting coil is placed in a magnetic field, the magnetic flux is Experiment No : EM 8 Experiment Name: Inductance of a Solenoid Objective: Investigation of the inductance of different solenoids and their dependence on certain parameters of solenoids Theoretical Information

More information

Simulation Study of Interaction between Energetic Ions and Alfvén Eigenmodes in LHD

Simulation Study of Interaction between Energetic Ions and Alfvén Eigenmodes in LHD 1 Simulation Study of Interaction between Energetic Ions and Alfvén Eigenmodes in LHD Y. Todo 1), N. Nakajima 1), M. Osakabe 1), S. Yamamoto 2), D. A. Spong 3) 1) National Institute for Fusion Science,

More information

Modeling Prey and Predator Populations

Modeling Prey and Predator Populations Modeling Prey and Predator Populations Alison Pool and Lydia Silva December 15, 2006 Abstract In this document, we will explore the modeling of two populations based on their relation to one another. Specifically

More information

Before you begin read these instructions carefully:

Before you begin read these instructions carefully: NATURAL SCIENCES TRIPOS Part IB & II (General) Friday, 31 May, 2013 9:00 am to 12:00 pm MATHEMATICS (2) Before you begin read these instructions carefully: You may submit answers to no more than six questions.

More information

Experiments with a Supported Dipole

Experiments 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 information

Robust control of resistive wall modes using pseudospectra

Robust control of resistive wall modes using pseudospectra Robust control of resistive wall modes using pseudospectra M. Sempf, P. Merkel, E. Strumberger, C. Tichmann, and S. Günter Max-Planck-Institut für Plasmaphysik, EURATOM Association, Garching, Germany GOTiT

More information

Exp. #2-4 : Measurement of Characteristics of Magnetic Fields by Using Single Coils and a Computer Interface

Exp. #2-4 : Measurement of Characteristics of Magnetic Fields by Using Single Coils and a Computer Interface PAGE 1/17 Exp. #2-4 : Measurement of Characteristics of Magnetic Fields by Using Single Coils and a Computer Interface Student ID Major Name Team No. Experiment Lecturer Student's Mentioned Items Experiment

More information

NIMROD simulations of dynamo experiments in cylindrical and spherical geometries. Dalton Schnack, Ivan Khalzov, Fatima Ebrahimi, Cary Forest,

NIMROD simulations of dynamo experiments in cylindrical and spherical geometries. Dalton Schnack, Ivan Khalzov, Fatima Ebrahimi, Cary Forest, NIMROD simulations of dynamo experiments in cylindrical and spherical geometries Dalton Schnack, Ivan Khalzov, Fatima Ebrahimi, Cary Forest, 1 Introduction Two experiments, Madison Plasma Couette Experiment

More information

DIAGNOSTICS FOR ADVANCED TOKAMAK RESEARCH

DIAGNOSTICS FOR ADVANCED TOKAMAK RESEARCH DIAGNOSTICS FOR ADVANCED TOKAMAK RESEARCH by K.H. Burrell Presented at High Temperature Plasma Diagnostics 2 Conference Tucson, Arizona June 19 22, 2 134 /KHB/wj ROLE OF DIAGNOSTICS IN ADVANCED TOKAMAK

More information

Electrode and Limiter Biasing Experiments on the Tokamak ISTTOK

Electrode and Limiter Biasing Experiments on the Tokamak ISTTOK Electrode and Limiter Biasing Experiments on the Tokamak ISTTOK C. Silva, H. Figueiredo, J.A.C. Cabral,, I. Nedzelsky, C.A.F. Varandas Associação Euratom/IST, Centro de Fusão Nuclear, Instituto Superior

More information

Modelling the Electromechanical Interactions in a Null-Flux EDS Maglev System

Modelling the Electromechanical Interactions in a Null-Flux EDS Maglev System Modelling the Electromechanical Interactions in a Null-Flux EDS Maglev System Jeroen de Boeij,2, Maarten Steinbuch 2 and Hector Gutiérrez Department of Mechanical & Aerospace Engineering Florida Institute

More information

GA A27910 THE SINGLE-DOMINANT MODE PICTURE OF NON-AXISYMMETRIC FIELD SENSITIVIITY AND ITS IMPLICATIONS FOR ITER GEOMETRIC TOLERANCES

GA A27910 THE SINGLE-DOMINANT MODE PICTURE OF NON-AXISYMMETRIC FIELD SENSITIVIITY AND ITS IMPLICATIONS FOR ITER GEOMETRIC TOLERANCES GA A279 THE SINGLE-DOMINANT MODE PICTURE OF NON-AXISYMMETRIC FIELD SENSITIVIITY AND ITS IMPLICATIONS FOR ITER by C. PAZ-SOLDAN, K.H. BURRELL, R.J. BUTTERY, J.S. DEGRASSIE, N.M. FERRARO, A.M. GAROFALO,

More information

Observation of Neo-Classical Ion Pinch in the Electric Tokamak*

Observation of Neo-Classical Ion Pinch in the Electric Tokamak* 1 EX/P6-29 Observation of Neo-Classical Ion Pinch in the Electric Tokamak* R. J. Taylor, T. A. Carter, J.-L. Gauvreau, P.-A. Gourdain, A. Grossman, D. J. LaFonteese, D. C. Pace, L. W. Schmitz, A. E. White,

More information

Role of the Electron Temperature in the Current Decay during Disruption in JT-60U )

Role of the Electron Temperature in the Current Decay during Disruption in JT-60U ) Role of the Electron Temperature in the Current Decay during Disruption in JT-60U ) Yoshihide SHIBATA, Akihiko ISAYAMA, Go MATSUNAGA, Yasunori KAWANO, Seiji MIYAMOTO 1), Victor LUKASH 2), Rustam KHAYRUTDINOV

More information

NIMEQ: MHD Equilibrium Solver for NIMROD

NIMEQ: MHD Equilibrium Solver for NIMROD NIMEQ: MHD Equilibrium Solver for NIMOD E.C.Howell, C..Sovinec University of Wisconsin-Madison 5 th Annual Meeting of Division of Plasma Physics Dallas, Texas, Nov. 17-Nov. 1,8 1 Abstract A Grad-Shafranov

More information

Princeton Plasma Physics Laboratory. Multi-mode analysis of RWM feedback with the NMA Code

Princeton Plasma Physics Laboratory. Multi-mode analysis of RWM feedback with the NMA Code Princeton Plasma Physics Laboratory Multi-mode analysis of RWM feedback with the NMA Code M. S. Chance, M.Okabayashi, M. S. Chu 12 th Workshop on MHD Stability Control: Improved MHD Control Configurations

More information

Physics of fusion power. Lecture 13 : Diffusion equation / transport

Physics of fusion power. Lecture 13 : Diffusion equation / transport Physics of fusion power Lecture 13 : Diffusion equation / transport Many body problem The plasma has some 10 22 particles. No description is possible that allows for the determination of position and velocity

More information

Radiative & Magnetohydrodynamic Shocks

Radiative & Magnetohydrodynamic Shocks Chapter 4 Radiative & Magnetohydrodynamic Shocks I have been dealing, so far, with non-radiative shocks. Since, as we have seen, a shock raises the density and temperature of the gas, it is quite likely,

More information

1 of 11 4/27/2013 8:48 PM

1 of 11 4/27/2013 8:48 PM A Dash of Maxwell's: A Maxwell's Equations Primer - Part 2: Why Things Radiate Written by Glen Dash, Ampyx LLC In Chapter I, I introduced Maxwell s Equations for the static case, that is, where charges

More information

Chapter 5 Summary 5.1 Introduction and Definitions

Chapter 5 Summary 5.1 Introduction and Definitions Chapter 5 Summary 5.1 Introduction and Definitions Definition of Magnetic Flux Density B To find the magnetic flux density B at x, place a small magnetic dipole µ at x and measure the torque on it: N =

More information

ITER operation. Ben Dudson. 14 th March Department of Physics, University of York, Heslington, York YO10 5DD, UK

ITER operation. Ben Dudson. 14 th March Department of Physics, University of York, Heslington, York YO10 5DD, UK ITER operation Ben Dudson Department of Physics, University of York, Heslington, York YO10 5DD, UK 14 th March 2014 Ben Dudson Magnetic Confinement Fusion (1 of 18) ITER Some key statistics for ITER are:

More information

Tidal Effects on Earth s Surface

Tidal Effects on Earth s Surface Tidal Effects on Earth s Surface Tom Murphy February, 1 This treatment follows the conventions of F. D. Stacey s Physics of the Earth, and is largely an elaboration on this work. 1 Tidal Potential The

More information

Experiment 2-6. Magnetic Field Induced by Electric Field

Experiment 2-6. Magnetic Field Induced by Electric Field Experiment 2-6. Magnetic Field Induced by Electric Field - Biot-Savart law and Ampere s Law - Purpose of Experiment We introduce concept called charge to describe electrical phenomenon. The simplest electrical

More information

Physics 4. Magnetic Induction. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Physics 4. Magnetic Induction. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB Physics 4 Magnetic Induction Before we can talk about induction we need to understand magnetic flux. You can think of flux as the number of field lines passing through an area. Here is the formula: flux

More information

Mathematical Properties of Stiffness Matrices

Mathematical Properties of Stiffness Matrices Mathematical Properties of Stiffness Matrices CEE 4L. Matrix Structural Analysis Department of Civil and Environmental Engineering Duke University Henri P. Gavin Fall, 0 These notes describe some of the

More information