MHD Linear Stability Analysis Using a Full Wave Code
|
|
- Brenda Burns
- 5 years ago
- Views:
Transcription
1 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 University Present affiliation: SGI Japan, Ltd
2 Outlines Introduction TASK/WA code Full Wave Analysis Derivation of resistive MHD dielectric tensor Benchmark tests Resistive wall mode including the effect of ferromagnetism Dispersion relation Comparisons with TASK/WA code Effect of plasma rotation Summary
3 Introduction Various kinds of extended MHD models Ideal MHD Resistive MHD Hall MHD Finite-electron-mass MHD Two-fluid MHD Kinetic models Uniform kinetic model: ω ω,ω d,ω B No FLR, Reduced FLR, Differential OP, Integral OP, Spectral Gyrokinetic model: ω ω B Hamiltonian model: three (adiabatic) constants of motion
4 Full Wave Analysis Electric field E with a complex frequency ω Maxwell s equation and the dielectric tensor 1 μ(r) [ E(r)] ω2 ɛ 0 dr ɛ (r, r ) E(r ) = i ω j ext vr) Type of analyses Antenna excitation: ω: real (boundary-value problem) Analysis of wave heating and current drive Loading impedance analysis (weakly stable eigenmode) Spontaneous excitation: ω: complex (eigen-value problem) Analysis of linear stability
5 Previous Analysis TASK/WM Fourier mode expansion in poloidal and toroidal directions Maxwell s equation and the dielectric tensor 1 μ 0 [ Emn (ρ) ] ω 2 ɛ 0 m n ɛ m m,n n (ρ) E m n (ρ) = i ω j ext(r) (Uniform kinetic + gyrokinetic) dielectric tensor Plasma dispersion function for parallel interaction Electronfinite mass Nofinite Larmor radius effects Finite difference method in radial direction Alfvén eigenmode analysis Excitation by energetic particle No coupling with drift waves
6 Present Analysis TASK/WA Fourier mode expansion in poloidal and toroidal directions Maxwell s equation and the dielectric tensor 1 [ Emn (ρ) ] ω 2 ɛ μ(ρ) 0 ɛ m m,n n (ρ) E m n (ρ) = i ω j ext(r) m n Resistive MHD dielectric tensor Finite element method in radial direction Higher numerical accuracy Comparison with previous Resistive MHD analyses Analysis of Resistive Wall Mode (RWM)
7 TASK/WA: Coordinates Flux coordinate: (s, Θ,ϕ) = (x 1, x 2, x 3 ) ϕ: Toroidal angle s = ψ p : Radial coordinate Θ: Poloidal angle (Boozer coordinate) k R = n B ϕ B + mr B Θ r B : constant k = n R B ϕ B + m r B Θ B : not constant k = 0 on a rational surface
8 TASK/WA: Finite Element Method Weak form of Maxwell s equation: (w: weighting function) ( w) 1 μ ( E) dv + ω2 ɛ 0 w ɛ E dv + i ω w j ext dv Interpolation function In plasma: second-order Lagrange polynomial In vacuum: linear function Boundary conditions On magnetic axis: E 2 se θ 0 for s 0 Finiteness of i ωb = E On plasma surface: Surface term due to a derivative of discontinuous parameters. On perfectly conducting surface: E 2 = 0, E 3 = 0
9 TASK/WA: Eigenmode analysis Excitation, j ext, proportional to the electron density Find complex ω which maximize the volume-integrated electric field amplitude Inverse of volume-integrate wave amplitude, phase and amplitude of the determinant of coefficient matrix as a function of the growth rate ω
10 Resistive MHD Dielectric Tensor Resistive MHD equation Linearization: X = X 0 + X 1 Laplace Transform: complex ω E = B t dp dt +Γp u = 0 ρ m du dt + p = j B E + u B = η j Flux coordinate: (s, Θ,ϕ)=(x 1, x 2, x 3 ) Co-variant and contra-variant expression
11 Matrix Expression (I) Definitions ū 0 = (0, u 2 0, u3 0 ), ū 1 = (u 1 1, u2 1, u3 1 ) B 0 = (0, B 2 0, B3 0 ), Ē 1 = (E 1,1, E 1,2, E 1,3 ) g = {g ij }, = (d/dx 1, d/dx 2, d/dx 3 ) u i 0 u j 1 Λα ij ū 0 Λ ū 1, ū 1 ū 0 + 2ū 0 Λ ū 1 U ū 1 j 1 B 0 B axis R j 1 1 ( ) k e j ω 0 k e j 0 I K ( j 0 ), k = ( id/dx 1, m, n): wave number ω = ω ū 0 k
12 Matrix expression (II) Adiabatic equation of state i ω p 1 = dp 0 dx 1u1 1 +Γp 0 u 1 Equation of motion ρ 0 ( i ω + ) U ū 1 + g p 1 = JB axis g R j 1 + g K ( j 0 ) Ē 1 Ohm s law g Ē 1 + JB axis g R ū1 + g K (ū 0 ) Ē 1 = η j 1 Finally we obtain the linear relation between ū 1 and Ē 1 [ηρ 0 ( i ω + U ) (JB axis ) 2 g R g R ] = ū 1 + η i ω g ( dp0 dx 1u1 1 +Γp 0 u 1 [ JB axis g R g + JBaxis g R g K (ū0 ) + η g K ( j 0 ) ) ] Ē 1
13 Dielectric Tensor Induced current: J 1 = σ E 1 Electrical conductivity tensor σ = σ 0 + d σ x dx 1 + d σ 1 R dx d d σ 1 R dx 1 1 dx 1 + d d σ 2 L dx 1 22 R dx d d d σ 2 L dx 1 22 R dx 1 1 dx 1 where σ 0 and σ x are 3 3 matrixes, σ 1 and σ 2 are 3 2 matrixes, R 0 and R 1 are 2 3 matrixes, L 22 is a 2 2 matrix. Dielectric tensor ɛ = I + i ωɛ 0 s σ s
14 Circular Waveguide Analytic solution k 2 z + ( j ml )2 a 2 = ω2 c MHz, 128.2MHz, 204.3MHz for k z = 1/3m 1, a = 2mCm = 1 Inverse of volume-integrated electric field amplitude 1.40E+01 log 10 (1/E 2 ) 1.00E E E E E E E E E+02 Re(ω)
15 Spectrum in a Cylindrical Plasma Square of frequency normalized by the Alfvén time τ A = R 0 /v A, Ω 2 = (ωτ A ) 2 From high frequency: fast wave, shear Alfvén wave, slow wave Good agreement with Fig. 5 in M. Chance, et. al, Nucl. Fusion 17 (1977) 65 Parameters Uniform density, B z, and q Poloidal mode number m = 2 nq = 1.9 Ω
16 External Kink Mode in a Cylindrical Plasma Flat current profile Analytic solution [ γ 2 q 2 aτ 2 A = 2(m nq a) 1 m nq ] a 1 (a/b) 2m a = 1m, b = 1.2aCL = 2πR 0 = 40π mcm = 2Cn = 1 Growth rate as a function of nq Analytic solution TASK/WA 0.20 γ 2 τa nq a
17 m = 1 Internal Kink Mode and Tearing Mode Growth rate as a function of plasma resistivity 10 1 TASK/WA Ideal internal kink τ H -2/3τ R -1/3 γτa η = Ω m η plasma (Ω m) E 2 / g 22 E 2 / g E E+03 η = Ω m 3.0E E E E s(ψ) 3.0E E E E s(ψ)
18 Spectrum in a Tokamak Plasma Spectrum of Incompressible plasma in a Solovev equilibrium cf. Fig. 2 in M. Chance, et. al, J. Comput. Phys 28 (1978) 1 ɛ = a/r 0 = 1/20, κ = 1, q 0 = , n = 10, m = 0, 1, TASK/WA Chance et al Ω
19 plasma vacuum I wall vacuum II Dispersion Relation of RWM Cylindrical plasma: (radius: a) Axial length: 2πR 0 (periodic) Poloidal mode number: m Toroidal mode number: n, k = n/r 0,(k 2 r 2 1) Surrounded by a resistive and ferromagnetic wall: (radius: b) r a b b + d 0 Wall thickness: d Wall permeability: μ W,(ˆμ μ W /μ 0 ) Wall resistivity: η W
20 Dispersion Relation of RWM (I) Response of the vacuum regions and the wall [ ] λ 2 (b + d) m2 ˆμ 2 (1 e 2dλ ) λmˆμ(1 + e 2dλ ) d b b ( a 2m α = [ ] λ 2 (b + d) + m2 ˆμ 2 ( (1 e 2dλ ) + λmˆμ(1 + e 2dλ ) 2 + d ) b) b b where λ 2 = ˆμγτ W bd Response of the plasma α = m2 b 2 + k2 2m (nq a + m)/(nq a m) +Δ a (γ)(1 + ˆγ 2 ) + ˆγ 2 m where the normalized growth rate is ˆγ = γτ A q a /(m nq a )
21 Dispersion Relation of RWM (II) Value of Δ a (γ) When the q profile is flat: Δ a (γ) = m 1 When the rational surface q(r mn ) = m/n is close to the plasma surface, r mn a: Δ a (γ) = Δ a (0)(1 + ˆγ 2 ) 1 1 Δ a (0) ( 1 1 (y2 + ˆγ 2 ) 1 dy ), Δ a (0) where ξ is a solution of the Euler equation f = rf2 k 2 0 d dr ( f dξ dr ), g = 2 k2 k0 2 (μ 0 p) + gξ = 0 ( k 2 0 r 2 ) 1 rf k2 k 2 0 r2 rk 4 0 ( ) rξ ξ ˆFF a F = kb z + m r B θ, k 2 0 = k 2 + m2 r 2, ˆF = kb z m r B θ
22 Numerical Results for RWM Comparison between the dispersion relation and the numerical results by TASK/WA Plasma resistivity η plasma is non-zero in TASK/WA Ideally conducting wall at r = c in TASK/WA Er /10 5 (AU) Typical Electric Field Profile (Re/Im) for b/a = 1.23, c/a = 2.0, m = 2, n = 1 E r E θ E θ in the wall r/a Eθ/10 4 (AU) r/a Eθ (AU) 1.24 r/a 1.25
23 Dependence on the Wall Radius b/a γτw ratio q profile: Flat q profile: Wesson type q(r) = q a = 1.35 q 0 = 0.8, q a = 1.8 for various values of d/a q(r) = q a r 2 /a dispersion relation (μˆ =1) (μˆ =10) TASK/WA (μˆ =1) (μˆ =10) (DR) / (WA code) (μˆ=1) (μˆ =10) b/a γτw (1 r 2 /a 2 ) q a/q 0 dispersion relation (μˆ =1) (μˆ =10) TASK/WA (μˆ =1) (μˆ =10) There is a soft threshold of the wall radius b/a. The dispersion relation correctly predicts γ. b/a
24 Destabilization by µ Wall > 1 (b/a = 1.23) γτw ratio q profile: Flat q profile: Wesson type q(r) = q a = 1.35 q 0 = 0.8, q a = 1.8 for various values of d/a wall thickness d/a = dispersion relation d/a=0.02 d/a=0.01 d/a=0.005 TASK/WA d/a=0.02 d/a=0.01 d/a=0.005 (DR)/(WA code) d/a=0.02 d/a=0.01 d/a= μˆ γτw dispersion relation TASK/WA γ increases with the increase of μ Wall. For flat q profiles, the dispersion relation correctly predicts γ. μˆ
25 Stabilization by Toroidal Rotation Plasma rotation stabilizes RWM when plasma dissipation exists It was confirmed that the kinetic dissipation is not important when the plasma is incompressible near the plasma edge. q profile: Wesson type (q 0 = 0.8Cq a = 1.8) Toroidal rotation dependence of the RWM growth rate for various values of η edge γτw η edge (Ω m) = Ωτ A
26 Comparison of Dissipation Mechanisms Resistivity only TASK/WA Resistivity and perp. viscosity R. Fitzpatrick and A. Aydemir Nucl. Fusion (1996) γτw Ωτ A r w =1.090 r w =1.200 r w =1.275 Perpendicular viscosity reduces the growth rate for faster rotation and longer wall distance.
27 Summary We have developed a full wave code TASK/WA in order to study the applicability of the full wave analysis to the linear stability study of global modes in a tokamak. Resistive MHD dielectric tensor was used in the present analysis and the results are compared with previous MHD studies. Good agreement was shown for cylindrical plasmas and for incompressible modes in tokamak plasmas. We analyzed the ferromagnetic effect on the resistive wall mode in a cylindrical plasma using a resistive MHD dielectric tensor. The growth rate of RWMs calculated from the dispersion relation agrees well with the results of the TASK/WA code. The critical plasma rotation which stabilizes RWMs increases almost linearly with the increase of ˆμ.
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 information2D full wave analysis of wave structure by TASK/WF2
17th NEXT (Numerical Experiment of Tokamak) Meeting, University of Tokyo, Kashiwa, Japan, March 15-16, 2012 2D full wave analysis of wave structure by TASK/WF2 Y. Maruyama, A. Fukuyama Department of Nuclear
More informationThe Linear Theory of Tearing Modes in periodic, cyindrical plasmas. Cary Forest University of Wisconsin
The Linear Theory of Tearing Modes in periodic, cyindrical plasmas Cary Forest University of Wisconsin 1 Resistive MHD E + v B = ηj (no energy principle) Role of resistivity No frozen flux, B can tear
More informationIntegrated Full Wave Analysis of RF Heating and Current Drive in Toroidal Plasmas
Integrated Full Wave Analysis of RF Heating and Current Drive in Toroidal Plasmas IAEA Fusion Energy Conference Chengdu, China 2006/10/20 A. Fukuyama, S. Murakami, A. Sonoda, M. Honda// Department of Nuclear
More informationHeating and current drive: Radio Frequency
Heating and current drive: Radio Frequency Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 13 th February 2012 Dr Ben Dudson Magnetic Confinement Fusion (1 of 26)
More 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 informationModeling of ELM Dynamics for ITER
Modeling of ELM Dynamics for ITER A.Y. PANKIN 1, G. BATEMAN 1, D.P. BRENNAN 2, A.H. KRITZ 1, S. KRUGER 3, P.B. SNYDER 4 and the NIMROD team 1 Lehigh University, 16 Memorial Drive East, Bethlehem, PA 18015
More informationControl of linear modes in cylindrical resistive MHD with a resistive wall, plasma rotation, and complex gain
Control of linear modes in cylindrical resistive MHD with a resistive wall, plasma rotation, and complex gain Dylan Brennan 1 and John Finn 2 contributions from Andrew Cole 3 1 Princeton University / PPPL
More informationActive and Fast Particle Driven Alfvén Eigenmodes in Alcator C-Mod
Active and Fast Particle Driven Alfvén Eigenmodes in Alcator C-Mod JUST DID IT. J A Snipes, N Basse, C Boswell, E Edlund, A Fasoli #, N N Gorelenkov, R S Granetz, L Lin, Y Lin, R Parker, M Porkolab, J
More informationA simple model of the resistive wall mode in tokamaks
A simple model of the resistive wall mode in tokamaks Richard Fitzpatrick Institute for Fusion Studies, Department of Physics, University of Texas at Austin, Austin TX 78712 (February 18, 2003) A simple
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 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 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 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 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 informationHigh-m Multiple Tearing Modes in Tokamaks: MHD Turbulence Generation, Interaction with the Internal Kink and Sheared Flows
TH/P3-3 High-m Multiple Tearing Modes in Tokamaks: MHD Turbulence Generation, Interaction with the Internal Kink and Sheared Flows A. Bierwage 1), S. Benkadda 2), M. Wakatani 1), S. Hamaguchi 3), Q. Yu
More informationFundamentals of Magnetic Island Theory in Tokamaks
Fundamentals of Magnetic Island Theory in Tokamaks Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Austin, TX, USA Talk available at http://farside.ph.utexas.edu/talks/talks.html
More informationGyrokinetic Transport Driven by Energetic Particle Modes
Gyrokinetic Transport Driven by Energetic Particle Modes by Eric Bass (General Atomics) Collaborators: Ron Waltz, Ming Chu GSEP Workshop General Atomics August 10, 2009 Outline I. Background Alfvén (TAE/EPM)
More informationSMR/ Summer College on Plasma Physics. 30 July - 24 August, Introduction to Magnetic Island Theory.
SMR/1856-1 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 Introduction to Magnetic Island Theory. R. Fitzpatrick Inst. for Fusion Studies University of Texas at Austin USA Introduction
More informationModelling of the penetration process of externally applied helical magnetic perturbation of the DED on the TEXTOR tokamak
INSTITUTE OF PHYSICS PUBLISHING Plasma Phys. Control. Fusion 8 (6) 69 8 PLASMA PHYSICS AND CONTROLLED FUSION doi:.88/7-/8// Modelling of the penetration process of externally applied helical magnetic perturbation
More informationL Aquila, Maggio 2002
Nonlinear saturation of Shear Alfvén Modes and energetic ion transports in Tokamak equilibria with hollow-q profiles G. Vlad, S. Briguglio, F. Zonca, G. Fogaccia Associazione Euratom-ENEA sulla Fusione,
More informationDispersive Media, Lecture 7 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2017-02-14 Dispersive Media, Lecture 7 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasmas as a coupled system Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas
More informationDOPPLER RESONANCE EFFECT ON ROTATIONAL DRIVE BY ION CYCLOTRON MINORITY HEATING
DOPPLER RESONANCE EFFECT ON ROTATIONAL DRIVE BY ION CYCLOTRON MINORITY HEATING V.S. Chan, S.C. Chiu, Y.A. Omelchenko General Atomics, San Diego, CA, U.S.A. 43rd Annual APS Division of Plasma Physics Meeting
More informationThe RFP: Plasma Confinement with a Reversed Twist
The RFP: Plasma Confinement with a Reversed Twist JOHN SARFF Department of Physics University of Wisconsin-Madison Invited Tutorial 1997 Meeting APS DPP Pittsburgh Nov. 19, 1997 A tutorial on the Reversed
More informationKinetic damping in gyro-kinetic simulation and the role in multi-scale turbulence
2013 US-Japan JIFT workshop on New Aspects of Plasmas Kinetic Simulation NIFS, November 22-23, 2013 Kinetic damping in gyro-kinetic simulation and the role in multi-scale turbulence cf. Revisit for Landau
More informationINTERACTION OF AN EXTERNAL ROTATING MAGNETIC FIELD WITH THE PLASMA TEARING MODE SURROUNDED BY A RESISTIVE WALL
INTERACTION OF AN EXTERNAL ROTATING MAGNETIC FIELD WITH THE PLASMA TEARING MODE SURROUNDED BY A RESISTIVE WALL S.C. GUO* and M.S. CHU GENERAL ATOMICS *Permanent Address: Consorzio RFX, Padova, Italy **The
More informationNonperturbative Effects of Energetic Ions on Alfvén Eigenmodes
1 TH/3-1Ra Nonperturbative Effects of Energetic Ions on Alfvén Eigenmodes Y. Todo 1), N. Nakajima 1), K. Shinohara 2), M. Takechi 2), M. Ishikawa 2), S. Yamamoto 3) 1) National Institute for Fusion Science,
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 informationMHD WAVES AND GLOBAL ALFVÉN EIGENMODES
MHD WVES ND GLOBL LFVÉN EIGENMODES S.E. Sharapov Euratom/CCFE Fusion ssociation, Culham Science Centre, bingdon, Oxfordshire OX14 3DB, UK S.E.Sharapov, Lecture 3, ustralian National University, Canberra,
More information2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasma physics Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas Transverse waves
More informationINTERACTION OF DRIFT WAVE TURBULENCE AND MAGNETIC ISLANDS
INTERACTION OF DRIFT WAVE TURBULENCE AND MAGNETIC ISLANDS A. Ishizawa and N. Nakajima National Institute for Fusion Science F. L. Waelbroeck, R. Fitzpatrick, W. Horton Institute for Fusion Studies, University
More informationMHD. Jeff Freidberg MIT
MHD Jeff Freidberg MIT 1 What is MHD MHD stands for magnetohydrodynamics MHD is a simple, self-consistent fluid description of a fusion plasma Its main application involves the macroscopic equilibrium
More informationTheoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions
Theoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions Tom Elsden 1 Andrew Wright 1 1 Dept Maths & Stats, University of St Andrews DAMTP Seminar - 8th May 2017 Outline Introduction Coordinates
More informationStability of the resistive wall mode in HBT-EP plasmas
Stability of the resistive wall mode in HBT-EP plasmas Richard Fitzpatrick Institute for Fusion Studies Department of Physics University of Texas at Austin Austin, TX 78712 James Bialek Department of Applied
More informationEquilibrium and transport in Tokamaks
Equilibrium and transport in Tokamaks Jacques Blum Laboratoire J.-A. Dieudonné, Université de Nice Sophia-Antipolis Parc Valrose 06108 Nice Cedex 02, France jblum@unice.fr 08 septembre 2008 Jacques Blum
More informationHybrid Kinetic-MHD simulations with NIMROD
simulations with NIMROD 1 Yasushi Todo 2, Dylan P. Brennan 3, Kwang-Il You 4, Jae-Chun Seol 4 and the NIMROD Team 1 University of Washington, Seattle 2 NIFS, Toki-Japan 3 University of Tulsa 4 NFRI, Daejeon-Korea
More informationNonlinear magnetohydrodynamic effects on Alfvén eigenmode evolution and zonal flow
Home Search Collections Journals About Contact us My IOPscience Nonlinear magnetohydrodynamic effects on Alfvén eigenmode evolution and zonal flow generation This article has been downloaded from IOPscience.
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 informationGyrokinetic Simulations of Tearing Instability
Gyrokinetic Simulations of Tearing Instability July 6, 2009 R. NUMATA A,, W. Dorland A, N. F. Loureiro B, B. N. Rogers C, A. A. Schekochihin D, T. Tatsuno A rnumata@umd.edu A) Center for Multiscale Plasma
More informationSupported by. Role of plasma edge in global stability and control*
NSTX Supported by College W&M Colorado Sch Mines Columbia U CompX General Atomics INL Johns Hopkins U LANL LLNL Lodestar MIT Nova Photonics New York U Old Dominion U ORNL PPPL PSI Princeton U Purdue U
More informationGlobal particle-in-cell simulations of Alfvénic modes
Global particle-in-cell simulations of Alfvénic modes A. Mishchenko, R. Hatzky and A. Könies Max-Planck-Institut für Plasmaphysik, EURATOM-Association, D-749 Greifswald, Germany Rechenzentrum der Max-Planck-Gesellschaft
More information20. Alfven waves. ([3], p ; [1], p ; Chen, Sec.4.18, p ) We have considered two types of waves in plasma:
Phys780: Plasma Physics Lecture 20. Alfven Waves. 1 20. Alfven waves ([3], p.233-239; [1], p.202-237; Chen, Sec.4.18, p.136-144) We have considered two types of waves in plasma: 1. electrostatic Langmuir
More informationA Study of 3-Dimensional Plasma Configurations using the Two-Fluid Plasma Model
A Study of 3-Dimensional Plasma Configurations using the Two-Fluid Plasma Model B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program University of Washington IEEE International Conference
More informationFast Ion Confinement in the MST Reversed Field Pinch
Fast Ion Connement in the MST Reversed Field Pinch Gennady Fiksel B. Hudson, D.J. Den Hartog, R.M. Magee, R. O'Connell, S.C. Prager MST Team - University of Wisconsin - Madison Center for Magnetic Self-Organization
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 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 informationA Hybrid Inductive Scenario for a Pulsed- Burn RFP Reactor with Quasi-Steady Current. John Sarff
A Hybrid Inductive Scenario for a Pulsed- Burn RFP Reactor with Quasi-Steady Current John Sarff 12th IEA RFP Workshop Kyoto Institute of Technology, Kyoto, Japan Mar 26-28, 2007 The RFP fusion development
More informationTwo Fluid Dynamo and Edge-Resonant m=0 Tearing Instability in Reversed Field Pinch
1 Two Fluid Dynamo and Edge-Resonant m= Tearing Instability in Reversed Field Pinch V.V. Mirnov 1), C.C.Hegna 1), S.C. Prager 1), C.R.Sovinec 1), and H.Tian 1) 1) The University of Wisconsin-Madison, Madison,
More informationEnergetic Particle Physics in Tokamak Burning Plasmas
Energetic Particle Physics in Tokamak Burning Plasmas presented by C. Z. (Frank) Cheng in collaboration with N. N. Gorelenkov, G. J. Kramer, R. Nazikian, E. Fredrickson, Princeton Plasma Physics Laboratory
More informationTransition From Single Fluid To Pure Electron MHD Regime Of Tearing Instability
Transition From Single Fluid To Pure Electron MHD Regime Of Tearing Instability V.V.Mirnov, C.C.Hegna, S.C.Prager APS DPP Meeting, October 27-31, 2003, Albuquerque NM Abstract In the most general case,
More informationFlow and dynamo measurements in the HIST double pulsing CHI experiment
Innovative Confinement Concepts (ICC) & US-Japan Compact Torus (CT) Plasma Workshop August 16-19, 211, Seattle, Washington HIST Flow and dynamo measurements in the HIST double pulsing CHI experiment M.
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 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 informationHeat Transport in a Stochastic Magnetic Field. John Sarff Physics Dept, UW-Madison
Heat Transport in a Stochastic Magnetic Field John Sarff Physics Dept, UW-Madison CMPD & CMSO Winter School UCLA Jan 5-10, 2009 Magnetic perturbations can destroy the nested-surface topology desired for
More informationInteraction of scrape-off layer currents with magnetohydrodynamical instabilities in tokamak plasmas
Interaction of scrape-off layer currents with magnetohydrodynamical instabilities in tokamak plasmas Richard Fitzpatrick Institute for Fusion Studies Department of Physics University of Texas at Austin
More informationSOLAR MHD Lecture 2 Plan
SOLAR MHD Lecture Plan Magnetostatic Equilibrium ü Structure of Magnetic Flux Tubes ü Force-free fields Waves in a homogenous magnetized medium ü Linearized wave equation ü Alfvén wave ü Magnetoacoustic
More informationStabilization of sawteeth in tokamaks with toroidal flows
PHYSICS OF PLASMAS VOLUME 9, NUMBER 7 JULY 2002 Stabilization of sawteeth in tokamaks with toroidal flows Robert G. Kleva and Parvez N. Guzdar Institute for Plasma Research, University of Maryland, College
More informationMHD Simulation of High Wavenumber Ballooning-like Modes in LHD
1 TH/P9-16 MHD Simulation of High Wavenumber Ballooning-like Modes in LHD H. Miura and N. Nakajima National Institute for Fusion Science, 322-6 Oroshi, Toki, Gifu 509-5292, JAPAN e-mail contact of main
More informationNonlinear MHD effects on TAE evolution and TAE bursts
Nonlinear MHD effects on TAE evolution and TAE bursts Y. Todo (NIFS) collaborating with H. L. Berk and B. N. Breizman (IFS, Univ. Texas) GSEP 3rd Annual Meeting (remote participation / Aug. 9-10, 2010)
More informationDT Fusion Ignition of LHD-Type Helical Reactor by Joule Heating Associated with Magnetic Axis Shift )
DT Fusion Ignition of LHD-Type Helical Reactor by Joule Heating Associated with Magnetic Axis Shift ) Tsuguhiro WATANABE National Institute for Fusion Science, 322-6 Oroshi-cho, Toki 509-5292, Japan (Received
More informationHybrid Kinetic-MHD simulations with NIMROD
in NIMROD simulations with NIMROD Charlson C. Kim 1 Dylan P. Brennan 2 Yasushi Todo 3 and the NIMROD Team 1 University of Washington, Seattle 2 University of Tulsa 3 NIFS, Toki-Japan December 2&3, 2011
More informationExperimental Study of the Stability of Alfvén Eigenmodes on JET
IAEA FEC, Paper EX/P-7 Experimental Study of the Stability of Alfvén Eigenmodes on JET D.Testa,, A.Fasoli,, G.Fu 4, A.Jaun 3, D.Borba, P.de Vries 6, and JET-EFDA contributors [] Plasma Science and Fusion
More informationMicrotearing Simulations in the Madison Symmetric Torus
Microtearing Simulations in the Madison Symmetric Torus D. Carmody, P.W. Terry, M.J. Pueschel - University of Wisconsin - Madison dcarmody@wisc.edu APS DPP 22 Overview PPCD discharges in MST have lower
More informationEffect of the Toroidal Asymmetry on the Structure of TAE Modes in Stellarators
Effect of the Toroidal Asymmetry on the Structure of TAE Modes in Stellarators Yu.V. Yakovenko 1), Ya.I. Kolesnichenko 1), V.V. Lutsenko 1), A. Weller 2), A. Werner 2) 1) Institute for Nuclear Research,
More information22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 19
. tability of the straight tokamak.65, MHD Theory of Fusion ystems Prof. Freidberg Lecture 9. ressure driven modes (uydams Criterion). internal modes 3. external modes. Tokamak Ordering Bθ ar B μ q or
More informationResponse of a Resistive and Rotating Tokamak to External Magnetic Perturbations Below the Alfven Frequency
Response of a Resistive and Rotating Tokamak to External Magnetic Perturbations Below the Alfven Freuency by M.S. Chu In collaboration with L.L. Lao, M.J. Schaffer, T.E. Evans E.J. Strait (General Atomics)
More informationarxiv: v1 [physics.plasm-ph] 11 Mar 2016
1 Effect of magnetic perturbations on the 3D MHD self-organization of shaped tokamak plasmas arxiv:1603.03572v1 [physics.plasm-ph] 11 Mar 2016 D. Bonfiglio 1, S. Cappello 1, M. Veranda 1, L. Chacón 2 and
More informationMeasuring from electron temperature fluctuations in the Tokamak Fusion Test Reactor
PHYSICS OF PLASMAS VOLUME 5, NUMBER FEBRUARY 1998 Measuring from electron temperature fluctuations in the Tokamak Fusion Test Reactor C. Ren, a) J. D. Callen, T. A. Gianakon, and C. C. Hegna University
More informationWaves in plasma. Denis Gialis
Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.
More informationControl of resistive wall modes in a cylindrical tokamak with plasma rotation and complex gain
Control of resistive wall modes in a cylindrical tokamak with plasma rotation and complex gain arxiv:146.5245v1 [physics.plasm-ph] 2 Jun 214 D. P. Brennan and J. M. Finn June 23, 214 Department of Astrophysical
More informationPLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT. Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford
1 PLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford 2 Overview Why do we need plasmas? For fusion, among other things
More informationW.A. HOULBERG Oak Ridge National Lab., Oak Ridge, TN USA. M.C. ZARNSTORFF Princeton Plasma Plasma Physics Lab., Princeton, NJ USA
INTRINSICALLY STEADY STATE TOKAMAKS K.C. SHAING, A.Y. AYDEMIR, R.D. HAZELTINE Institute for Fusion Studies, The University of Texas at Austin, Austin TX 78712 USA W.A. HOULBERG Oak Ridge National Lab.,
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 informationA Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1
A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1 B. Srinivasan 2, U. Shumlak Aerospace and Energetics Research Program University of Washington,
More informationENERGETIC PARTICLES AND BURNING PLASMA PHYSICS
ENERGETIC PARTICLES AND BURNING PLASMA PHYSICS Reported by J. Van Dam Institute for Fusion Studies The University of Texas at Austin US-Japan JIFT Workshop on Theory-Based Modeling and Integrated Simulation
More informationLow Beta MHD Equilibrium Including a Static Magnetic Island for Reduced MHD Equations in a Straight Heliotron Configuration
Low Beta MHD Equilibrium Including a Static Magnetic Island for Reduced MHD Equations in a Straight Heliotron Configuration Kinya SAITO 1,a), Katsuji ICHIGUCHI 1,2) and Ryuichi ISHIZAKI 1,2) 1) The Graduate
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationSurface currents associated with external kink modes in tokamak
Surface currents associated with external kink modes in tokamak plasmas during a major disruption C. S. Ng 1,2,a), and A. Bhattacharjee 1 1 Princeton Plasma Physics Laboratory, Princeton University, Princeton,
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 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 informationObservation of Co- and Counter Rotation Produced by Lower Hybrid Waves in Alcator C-Mod*
Observation of Co- and Counter Rotation Produced by Lower Hybrid Waves in Alcator C-Mod* R. R. Parker, Y. Podpaly, J. Lee, M. L. Reinke, J. E. Rice, P.T. Bonoli, O. Meneghini, S. Shiraiwa, G. M. Wallace,
More informationNonlinear hydrid simulations of precessional Fishbone instability
Nonlinear hydrid simulations of precessional Fishbone instability M. Faganello 1, M. Idouakass 1, H. L. Berk 2, X. Garbet 3, S. Benkadda 1 1 Aix-Marseille University, CNRS, PIIM UMR 7345, Marseille, France
More informationResponse of a Resistive and Rotating Tokamak to External Magnetic Perturbations Below the Afvén Frequency
1 THS/P5-04 Response of a Resistive and Rotating Tokamak to External Magnetic Perturbations Below the Afvén Frequency M.S. Chu, 1 L.L. Lao, 1, M.J. Schaffer 1 T.E. Evans, 1 E.J. Strait, 1 Y.Q. Liu 2, M.J.
More informationarxiv: v1 [physics.plasm-ph] 24 Nov 2017
arxiv:1711.09043v1 [physics.plasm-ph] 24 Nov 2017 Evaluation of ideal MHD mode stability of CFETR baseline scenario Debabrata Banerjee CAS Key Laboratory of Geospace Environment and Department of Modern
More informationEffect of magnetic reconnection on CT penetration into magnetized plasmas
Earth Planets Space, 53, 547 55, 200 Effect of magnetic reconnection on CT penetration into magnetized plasmas Yoshio Suzuki, Takaya Hayashi 2, and Yasuaki Kishimoto Naka Fusion esearch Establishment,
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 informationDamping and drive of low frequency modes in tokamak plasmas
Damping and drive of low frequency modes in tokamak plasmas Ph. Lauber, S. Günter Max-Planck-Institut für Plasmaphysik, EURATOM-Association, Garching, Germany E-mail: Philipp.Lauber@ipp.mpg.de Abstract.
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 informationCurrent Profile Control by ac Helicity Injection
Current Profile Control by ac Helicity Injection Fatima Ebrahimi and S. C. Prager University of Wisconsin- Madison APS 2003 Motivations Helicity injection is a method to drive current in plasmas in which
More informationM. T. Beidler 1, J. D. Callen 1, C. C. Hegna 1, C. R. Sovinec 1, N. M. Ferraro 2
P1.015 Nonlinear Modeling Benchmarks of Forced Magnetic Reconnection with NIMROD and M3D-C1 M. T. Beidler 1, J. D. Callen 1, C. C. Hegna 1, C. R. Sovinec 1, N. M. Ferraro 2 1 Department of Engineering
More informationCalculation of alpha particle redistribution in sawteeth using experimentally reconstructed displacement eigenfunctions
Calculation of alpha particle redistribution in sawteeth using experimentally reconstructed displacement eigenfunctions R. Farengo, H. E. Ferrari,2, M.-C. Firpo 3, P. L. Garcia-Martinez 2,3, A. F. Lifschitz
More 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 informationMomentum transport from magnetic reconnection in laboratory an. plasmas. Fatima Ebrahimi
Momentum transport from magnetic reconnection in laboratory and astrophysical plasmas Space Science Center - University of New Hampshire collaborators : V. Mirnov, S. Prager, D. Schnack, C. Sovinec Center
More informationWater-bag reduced gyrokinetic model for the study of the spatial structure of ITG instability linear modes
Water-bag reduced gyrokinetic model for the study of the spatial structure of ITG instability linear modes PhD work directed by N. Besse with invaluable input from P. Bertrand, E. Gravier, P. Morel,R.
More informationAPPENDIX Z. USEFUL FORMULAS 1. Appendix Z. Useful Formulas. DRAFT 13:41 June 30, 2006 c J.D Callen, Fundamentals of Plasma Physics
APPENDIX Z. USEFUL FORMULAS 1 Appendix Z Useful Formulas APPENDIX Z. USEFUL FORMULAS 2 Key Vector Relations A B = B A, A B = B A, A A = 0, A B C) = A B) C A B C) = B A C) C A B), bac-cab rule A B) C D)
More informationPlasma heating in stellarators at the fundamental ion cyclotron frequency
PHYSICS OF PLASMAS VOLUME 7, NUMBER FEBRUARY 000 Plasma heating in stellarators at the fundamental ion cyclotron frequency V. A. Svidzinski and D. G. Swanson Department of Physics, Auburn University, Auburn,
More informationThe Field-Reversed Configuration (FRC) is a high-beta compact toroidal in which the external field is reversed on axis by azimuthal plasma The FRC is
and Stability of Field-Reversed Equilibrium with Toroidal Field Configurations Atomics General Box 85608, San Diego, California 92186-5608 P.O. APS Annual APS Meeting of the Division of Plasma Physics
More informationExponential Growth and Filamentary Structure of Nonlinear Ballooning Instability 1
Exponential Growth and Filamentary Structure of Nonlinear Ballooning Instability 1 Ping Zhu in collaboration with C. C. Hegna and C. R. Sovinec University of Wisconsin-Madison Sherwood Conference Denver,
More informationarxiv:physics/ v1 [physics.plasm-ph] 14 Nov 2005
arxiv:physics/0511124v1 [physics.plasm-ph] 14 Nov 2005 Early nonlinear regime of MHD internal modes: the resistive case M.C. Firpo Laboratoire de Physique et Technologie des Plasmas (C.N.R.S. UMR 7648),
More informationlim = F F = F x x + F y y + F z
Physics 361 Summary of Results from Lecture Physics 361 Derivatives of Scalar and Vector Fields The gradient of a scalar field f( r) is given by g = f. coordinates f g = ê x x + ê f y y + ê f z z Expressed
More information