Magnetic Control of Perturbed Plasma Equilibria
|
|
- Darcy Heath
- 5 years ago
- Views:
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 O. Katsuro-Hopkins, J. Bialek, G. Navratil (Department of Applied Physics and Applied Mathematics, Columbia University, New York, NY USA) Building
More informationThe 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 informationRESISTIVE 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 informationDynamical 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 informationDynamical 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 information35. 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 informationPlasma 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 informationRWM 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 informationDerivation 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 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 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 informationRWM 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 informationPlasma 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 informationResistive 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 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 informationExtended 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 informationA 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 informationModel 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 informationReconstruction 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 informationDIII 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 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 informationModeling 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 informationAN 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 informationReconstruction 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 informationDesign 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 informationKSTAR 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 informationGoing 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 informationTokamak 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 informationwhere 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 informationEffect 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 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 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 information- 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 informationPlasma 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 informationRequirements 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 informationThree 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 informationComputations 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 informationResistive 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 informationBefore 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 informationTransmission 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 informationModel 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 informationPlasma 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 informationImposed 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 informationBunno, 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 informationSimple 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 informationGERALD 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 informationDisruption 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 informationReview 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 informationEffects 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 informationReview 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 informationINITIAL 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 informationELECTROMAGNETIC 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 informationFrom 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 informationDisruption 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 informationStellarators. 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 informationToroidal 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 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 informationThe 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 informationCurrent 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 informationBehavior 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 informationShear 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 informationLinear 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 informationNonlinear 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 informationSTABILIZATION 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 informationW15D1: 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 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 informationExperiment 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 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 informationCurrent-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 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 informationControl 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 informationEffect 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 informationEffects 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 informationThe Levitated Dipole Experiment: Towards Fusion Without Tritium
The Levitated Dipole Experiment: Towards Fusion Without Tritium Jay Kesner MIT M.S. Davis, J.E. Ellsworth, D.T. Garnier, M.E. Mauel, P.C. Michael, P.P. Woskov MCP I3.110 Presented at the EPS Meeting, Dublin,
More informationThe 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 informationN 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 informationSimulation 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 informationModeling 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 informationBefore 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 informationExperiments with a Supported Dipole
Experiments with a Supported Dipole Reporting Measurements of the Interchange Instability Excited by Electron Pressure and Centrifugal Force Introduction Ben Levitt and Dmitry Maslovsky Collisionless Terrella
More informationRobust 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 informationExp. #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 informationNIMROD 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 informationDIAGNOSTICS 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 informationElectrode 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 informationModelling 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 informationGA 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 informationObservation 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 informationRole 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 informationNIMEQ: 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 informationPrinceton 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 informationPhysics 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 informationRadiative & 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 information1 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 informationChapter 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 informationITER 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 informationTidal 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 informationExperiment 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 informationPhysics 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 informationMathematical 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