Hybrid Kinetic-MHD simulations with NIMROD
|
|
- Logan Walton
- 5 years ago
- Views:
Transcription
1 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
2 Outline in NIMROD 1 in NIMROD intro to NIMROD hybrid kinetic-mhd the equations of motion and δf 2 3 Components of δf 3 3 regions of stability map
3 Outline in NIMROD intro to NIMROD hybrid kinetic-mhd the equations of motion and δf 1 in NIMROD intro to NIMROD hybrid kinetic-mhd the equations of motion and δf 2 3 Components of δf 3 3 regions of stability map
4 in NIMROD NIMROD C.R. Sovinec, JCP, 195, 2004 intro to NIMROD hybrid kinetic-mhd the equations of motion and δf parallel, 3-D, initial value extended MHD code 2D high order finite elements + Fourier in symmetric direction linear and nonlinear simulations semi-implicit and implicit time advance operators model fusion relevant geometry and physical parameters magnetic Lundquist number - S = τ η /τ A large thermal anisotropy - χ χ 10 8 active developer and user base continually expanding capabilities large and growing Verification&Validation
5 in NIMROD intro to NIMROD hybrid kinetic-mhd the equations of motion and δf Example NIMROD grid - D3D with conformal vacuum
6 in NIMROD intro to NIMROD hybrid kinetic-mhd the equations of motion and δf NIMROD s Extended MHD Equations B t = E+κ divb ( B) J = 1 µ 0 B E = V B+ηJ [ + me e (J B p e) n ee 2 m e + J ] t + (JV+VJ) n α + (nv) α = D n α ( t ) V ρ t +V V = J B p + ρν V Π p h ( ) n α Tα +V α T α = p α V α Γ 1 t q α+q α Π α : V α resistive MHD Hall and 2-fluid Braginski and beyond closures energetic particles
7 in NIMROD Equations intro to NIMROD hybrid kinetic-mhd the equations of motion and δf momentum equation modified by hot particle pressure tensor: ( ) U ρ t +U U = J B p b p h b,h denote bulk plasma and hot particles ρ,u for entire plasma, both bulk and hot particle quasi-neutrality n e = n i +n h steady state equation J 0 B 0 = p 0 = p b0 + p h0 p b0 is scaled to accomodate hot particles scalar pressure requires isotropic velocity distribution alternative J h current coupling possible
8 in NIMROD intro to NIMROD hybrid kinetic-mhd the equations of motion and δf Schematic of Hybrid δf PIC-MHD model advance particles and δf 1 deposit δp(η) = space z n+1 i δf n+1 i = z n i +ż(z i ) t = δf n i + δf(z i ) t N δf i m(v i V h ) 2 S(η η i ) on FE logical i=1 advance NIMROD hybrid kinetic-mhd with modified momentum equation ρ s δu t = J s δb+δj B s δp b δp h 1 S. E. Parker and W. W. Lee, A fully nonlinear characteristic method for gyro-kinetic simulation, Physics of Fluids B, 5, 1993
9 in NIMROD Summary of PIC Capabilities intro to NIMROD hybrid kinetic-mhd the equations of motion and δf tracers, linear, (nonlinear) two equations of motion drift kinetic (v,µ), Lorentz force ( v) multiple spatial profiles - loading in x proportional to MHD profile, uniform, peaked gaussian multiple distribution functions - loading in v slowing down distribution, Maxwellian, monoenergetic multispecies capability initial implementation as self-consistent tracers e.g. - Maxwellian+slowing down, drift+lorentz
10 in NIMROD Drift Kinetic Equation of Motion intro to NIMROD hybrid kinetic-mhd the equations of motion and δf follows gyrocenter in limit of zero Larmour radius [ ] 1 2 reduces 6D to 4D +1 x(t),v (t),µ = mv2 B drift kinetic equations of motion ẋ = v ˆb+v D +v E B v D = m eb 4 ( v 2 + v2 2 v E B = E B B 2 m v = ˆb (µ B ee) ) )(B B2 + µ 0mv 2 2 eb 2 J
11 in NIMROD Drift Kinetic δp h PIC Moment intro to NIMROD hybrid kinetic-mhd the equations of motion and δf for drift kinetic equations, assume CGL-like δp 0 0 δp h = 0 δp δp evaluate pressure moment at x δp(x) = m v V h v V h δf(x,v)d 3 v δf is perturbed phase space density, m mass of particle, and V h is COM velocity of particles implemented as δp = δp I+δ pbb, δ p = δp δp
12 in NIMROD δf PIC reduces 1/ N noise S. E. Parker and W. W. Lee, PFB 5 (1993) intro to NIMROD hybrid kinetic-mhd the equations of motion and δf δf PIC reduces the discrete 1/ N particle noise Vlasov Equation f(z) t +ż f(z) z = 0 split phase space f = f eq (z)+δf(z,t) δf evolves along the characteristics ż (equations of motion) δf = δz f eq z using ż = z eq + δz and ż eq feq z = 0 moments of distribution function are A n = v n δf dv
13 in NIMROD intro to NIMROD hybrid kinetic-mhd the equations of motion and δf Slowing Down Distribution for Hot Particles slowing down distribution function f eq = P 0exp( P ζ ψ 0 ) ε 3/2 +ε 3/2 c P ζ = gρ ψ canonical toroidal momentum, ε energy, ψ p poloidal flux, ψ 0 gradient scale length, ε c critical energy { [( ) ] δf mg = f eq eψ 0 B 3 v 2 + v2 δb B µ 0 v 2 J δe } + δv ψ p + 3 ψ 0 2 eε 1/2 ε 3/2 +ε 3/2 c v D δe δv = δe B B 2 +v δb B v D = m ( ) eb 3 v 2 + v2 (B B)+ µ 0mv 2 2 eb 2 J
14 Outline in NIMROD 3 Components of δf 1 in NIMROD intro to NIMROD hybrid kinetic-mhd the equations of motion and δf 2 3 Components of δf 3 3 regions of stability map
15 in NIMROD 3 Components of δf Benchmark of Drift Kinetic (1,1) Kink With M3D C. C. Kim, PoP (2008) circular, monotonic q, q 0 =.6, q a = 2.5, β 0 =.08, R/a = 2.76, E hmax = 10KeV, dt=1e-7, τ A = 1.e6
16 in NIMROD 3 Components of δf movie1
17 in NIMROD 3 Components of δf Evolution of n = 1 δp hot δp hot t = 0 t.2ω h t.4ω h δ p hot tr = 0 t.2ω h t.4ω h
18 in NIMROD 3 Components of δf Orthogonal view of n = 1, δf(v,v ) t = 0 t.1ω h t.2ω h t.25ω h t.35ω h t.4ω h
19 in NIMROD 3 Components of δf Summary of Observations of Moments of δf n=1 energetic particles cause real frequency response phase lag between particles and fluid δ p concentrated on outboard side trapped particles concentration of δf n=1 (v,v ) activity to trapped cone? what is origin and role of the asymmetry? what is role of trapped vs. passing particles? what is the structure in the trapped cone
20 in NIMROD 3 Components of δf 3 Components of δf 1 f eq δf = [( ) ] mg v 2 eψ 0B 3 + v2 δb B µ 0v J δe (1) 2 + δv ψp (2)+ 3 ψ 0 2 eε 1/2 ε 3/2 +ε 3/2 c v D δe(3) δv = δe B +v B 2 δb B ( ) m v D = v 2 eb 3 + v2 (B B)+ µ0mv2 J 2 eb 2 δf has 3 components 1 gρ - kinetic term 2 δv E B ψ - radial particle flux 3 v D δe - energy exchange examine δf moments of components, e.g. (v D δe)δfdv convolution is axisymmetric, i.e. n = 0 and static
21 in NIMROD 3 Components of δf 0D Comparison of (1, 1) benchmark normalized amplitude (wrt δv E B ψ) %(passing/trapped) contribution β frac gρ δv E B ψ v D δe (-13/113) (69/31) (45/55) (ε/100) (58/42) (27/73) (0/100) (52/48) (18/82) δv E B ψ comparable passing trapped (destabilizing) gρ &v D δe small, negative, and trapped (stabilizing)
22 in NIMROD 3 Components of δf 3 Components in Velocity Space (β frac =.50)
23 Outline in NIMROD 3 regions of stability map 1 in NIMROD intro to NIMROD hybrid kinetic-mhd the equations of motion and δf 2 3 Components of δf 3 3 regions of stability map
24 in NIMROD 3 regions of stability map β N DIII-D experiments indicate sensitivity to increases in δβ N La Haye - Physics of Plasmas, 17 (2010), Nuclear Fusion, 51 (2011) n = 1 MHD Stability Map in (q min,β N ) Ideal MHD Stability Boundary Reconstruction Flattop γτ A MHD Only (n=1) S=τ R /τ A ~ Pr=100 Experimental Trajectory q min contours of MHD growth rate (NIMROD) - γτ A as q min Ideal MHD Stability Boundary in red (PEST) - stable to right Experimental Trajectory primarily in stable region of MHD stability map β N trigger 2/1 in experiments
25 in NIMROD n = 1 DIII-D linear simulations 3 regions of stability map 6 a) x b) c) P 0 e -2ψ/ψ 0 10 q P (pa) (m) 0 q=1.5 (m) q= ψ (Wb) R (m) 1 2 R (m) understand the stability properties of D3D hybrid discharge hybrid discharge sensitive to small changes in β N explore range of q = [.9,1.4] β N = [2,4]
26 in NIMROD q scan at fixed q-profile increase in q accompanied by decrease in β N 3 regions of stability map (γ,ω)τ A q min γ and ω a) NIMROD Pressure 1.2 x Experiment γτ A MHD Only γτ A with particles ωτ A with particles β N 4 3 β N energetic particles stabilizing in ideal region map divides into 3 regions ideal MHD unstable high frequency low frequency abrupt transition at q 1.25 m = 1 m = 2
27 in NIMROD 3 regions of stability map β N Energetic particles may help explain sensitivity to δβ N Kinetic-MHD Stability Map in (q min,β N ) Ideal MHD Stability Boundary Flattop Reconstruction γτ A Energetic Particle Driven MHD Modes (n=1) q min 0.01 Experimental Trajectory contours of γτ A have dependence on β N Ideal MHD Stability Boundary in red EP s drive instabilities in MHD stable region
28 a) in NIMROD 3 regions of stability map Spatial Eigenmodes show distinct topologies Br Vr P P per P ani b) c) R B r R Br R V r R Vr R R R P Pper Pani P R R Pper R Pani starts as dominantly core m = 1 at higher q more m = 2 abrupt shift to m = 2 dominant at q > 1.25 even some m = 3 d) R R Br Vr R P R Pper R Pani R R R R R
29 in NIMROD 3 regions of stability map Velocity space eigenmodes have more distinct features v (m/s) v (m/s) v (m/s) x x x10 δf total δf total δf total a) (m/s) x v δf r<r(q min ) δf r<r(q min ) δf r<r(q min ) (m/s) v x δf r>r(q min ) v 1 (m/s) 2 x (m/s) 2 x (m/s) 2 x10 6 b) (m/s) x v v x10 6 c) (m/s) x v v v (m/s) v x δf r>r(q min ) v x v 1 2 x10 6 δf r>r(q min ) (m/s) v x v x10 6 v x10 6 v x10 6 δf(v,v ) = r 2 r 1 δf(r,v)dr velocity eigenmodes have distinct structures r qmin fixed for q-scan divide into inner and outer region q min =[1.04,1.21,1.36] transition from inner to outer dominance for q > 1.25 transition from co-direction to counter-direction dominance
30 in NIMROD q = D diagnostic vs. time 3 regions of stability map phase volume integrated vs time, total(green), passing(blue), trapped(red) comparable amplitudes from δv E B ψ and v D δe passing particles have larger contribution!
31 in NIMROD 3 regions of stability map Velocity space plots detail 0D diagnostic significant v resonance
32 in NIMROD q = D diagnostic vs. time 3 regions of stability map δv E B ψ > v D δe, 5 larger passing particles have larger contribution!
33 in NIMROD 3 regions of stability map Velocity space plots detail 0D diagnostic δv E B ψ shows (v,v ) localized activity
34 Conclusion in NIMROD 3 regions of stability map new phase space diagnostics have been presented identifying interesting phenomenological features eigenmode in (v,v ) often more distinctive than spatial eigenmode distinct resonances features in trapped cone co-/counter- passing dominance both trapped and passing particles is important passing particles are crucial and often a dominant contribution from phenomenology, move to quantification, then theory
35 in NIMROD 3 regions of stability map movie2
Hybrid 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 informationHybrid Kinetic-MHD simulations Status and Updates
in NIMROD simulations Status and Updates Charlson C. Kim 1,2 Yasushi Todo 2 and the NIMROD Team 1. University of Washington, Seattle 2. National Institute for Fusion Science NIMROD Team Meeting Austin,
More informationHybrid Kinetic-MHD simulations with NIMROD
simulations with NIMROD Charlson C. Kim and the NIMROD Team Plasma Science and Innovation Center University of Washington CSCADS Workshop July 18, 2011 NIMROD C.R. Sovinec, JCP, 195, 2004 parallel, 3-D,
More informationThe δf PIC-MHD model in NIMROD
The δf PIC-MHD model in NIMROD Charlson C. Kim PSI Center U. Washington NIMROD Team Meeting Madison, WI July 29, 2009 Outline hybrid kinetic-mhd model PIC and δf and the δf PIC-MHD model the details drift
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 informationMHD-particle simulations and collective alpha-particle transport: analysis of ITER scenarios and perspectives for integrated modelling
MHD-particle simulations and collective alpha-particle transport: analysis of ITER scenarios and perspectives for integrated modelling G. Vlad, S. Briguglio, G. Fogaccia, F. Zonca Associazione Euratom-ENEA
More informationEffects of Alpha Particle Transport Driven by Alfvénic Instabilities on Proposed Burning Plasma Scenarios on ITER
Effects of Alpha Particle Transport Driven by Alfvénic Instabilities on Proposed Burning Plasma Scenarios on ITER G. Vlad, S. Briguglio, G. Fogaccia, F. Zonca Associazione Euratom-ENEA sulla Fusione, C.R.
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 informationEffects of fast ion phase space modifications by instabilities on fast ion modeling
Effects of fast ion phase space modifications by instabilities on fast ion modeling M. Podestà, M. Gorelenkova, E. Fredrickson, N. Gorelenkov, R. White PPPL, Princeton USA Acknowledgements: NSTX-U and
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 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 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 informationEnergetic-Ion-Driven MHD Instab. & Transport: Simulation Methods, V&V and Predictions
Energetic-Ion-Driven MHD Instab. & Transport: Simulation Methods, V&V and Predictions 7th APTWG Intl. Conference 5-8 June 2017 Nagoya Univ., Nagoya, Japan Andreas Bierwage, Yasushi Todo 14.1MeV 10 kev
More informationNonlinear Simulation of Energetic Particle Modes in JT-60U
TH/P6-7 Nonlinear Simulation of Energetic Particle Modes in JT-6U A.Bierwage,N.Aiba 2, K.Shinohara 2, Y.Todo 3,W.Deng 4,M.Ishikawa 2,G.Matsunaga 2 and M. Yagi Japan Atomic Energy Agency (JAEA), Rokkasho,
More informationFinite-Orbit-Width Effect and the Radial Electric Field in Neoclassical Transport Phenomena
1 TH/P2-18 Finite-Orbit-Width Effect and the Radial Electric Field in Neoclassical Transport Phenomena S. Satake 1), M. Okamoto 1), N. Nakajima 1), H. Sugama 1), M. Yokoyama 1), and C. D. Beidler 2) 1)
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 informationBounce-averaged gyrokinetic simulations of trapped electron turbulence in elongated tokamak plasmas
Bounce-averaged gyrokinetic simulations of trapped electron turbulence in elongated tokamak plasmas Lei Qi a, Jaemin Kwon a, T. S. Hahm a,b and Sumin Yi a a National Fusion Research Institute (NFRI), Daejeon,
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 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 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 informationAnalysis and modelling of MHD instabilities in DIII-D plasmas for the ITER mission
Analysis and modelling of MHD instabilities in DIII-D plasmas for the ITER mission by F. Turco 1 with J.M. Hanson 1, A.D. Turnbull 2, G.A. Navratil 1, C. Paz-Soldan 2, F. Carpanese 3, C.C. Petty 2, T.C.
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 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 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 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 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 informationUnderstanding Edge Harmonic Oscillation Physics Using NIMROD
Understanding Edge Harmonic Oscillation Physics Using NIMROD J. King With contributions from S. Kruger & A. Pankin (Tech-X); K. Burrell, A. Garofalo, R. Groebner & P. Snyder (General Atomics) Work supported
More informationGyrokinetic simulations of magnetic fusion plasmas
Gyrokinetic simulations of magnetic fusion plasmas Tutorial 2 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr
More 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 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 informationHighlights from (3D) Modeling of Tokamak Disruptions
Highlights from (3D) Modeling of Tokamak Disruptions Presented by V.A. Izzo With major contributions from S.E. Kruger, H.R. Strauss, R. Paccagnella, MHD Control Workshop 2010 Madison, WI ..onset of rapidly
More informationComparison of Kinetic and Extended MHD Models for the Ion Temperature Gradient Instability in Slab Geometry
Comparison of Kinetic and Extended MHD Models for the Ion Temperature Gradient Instability in Slab Geometry D. D. Schnack University of Wisconsin Madison Jianhua Cheng, S. E. Parker University of Colorado
More informationMAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT
MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT ABSTRACT A. G. Tarditi and J. V. Shebalin Advanced Space Propulsion Laboratory NASA Johnson Space Center Houston, TX
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 informationIntroduction to Fusion Physics
Introduction to Fusion Physics Hartmut Zohm Max-Planck-Institut für Plasmaphysik 85748 Garching DPG Advanced Physics School The Physics of ITER Bad Honnef, 22.09.2014 Energy from nuclear fusion Reduction
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationReduced MHD. Nick Murphy. Harvard-Smithsonian Center for Astrophysics. Astronomy 253: Plasma Astrophysics. February 19, 2014
Reduced MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 19, 2014 These lecture notes are largely based on Lectures in Magnetohydrodynamics by Dalton
More informationSawteeth in Tokamaks and their relation to other Two-Fluid Reconnection Phenomena
Sawteeth in Tokamaks and their relation to other Two-Fluid Reconnection Phenomena S. C. Jardin 1, N. Ferraro 2, J. Chen 1, et al 1 Princeton Plasma Physics Laboratory 2 General Atomics Supported by the
More informationNIMROD FROM THE CUSTOMER S PERSPECTIVE MING CHU. General Atomics. Nimrod Project Review Meeting July 21 22, 1997
NIMROD FROM THE CUSTOMER S PERSPECTIVE MING CHU General Atomics Nimrod Project Review Meeting July 21 22, 1997 Work supported by the U.S. Department of Energy under Grant DE-FG03-95ER54309 and Contract
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 informationGA A THERMAL ION ORBIT LOSS AND RADIAL ELECTRIC FIELD IN DIII-D by J.S. degrassie, J.A. BOEDO, B.A. GRIERSON, and R.J.
GA A27822 THERMAL ION ORBIT LOSS AND RADIAL ELECTRIC FIELD IN DIII-D by J.S. degrassie, J.A. BOEDO, B.A. GRIERSON, and R.J. GROEBNER JUNE 2014 DISCLAIMER This report was prepared as an account of work
More informationTheory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus. K. C. Shaing
Theory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus K. C. Shaing Plasma and Space Science Center, and ISAPS, National Cheng Kung University, Tainan, Taiwan 70101, Republic
More informationInnovative Concepts Workshop Austin, Texas February 13-15, 2006
Don Spong Oak Ridge National Laboratory Acknowledgements: Jeff Harris, Hideo Sugama, Shin Nishimura, Andrew Ware, Steve Hirshman, Wayne Houlberg, Jim Lyon Innovative Concepts Workshop Austin, Texas February
More informationProgress in the Plasma Science and Innovation Center
1 THP/2-02 Progress in the Plasma Science and Innovation Center U. Shumlak 1, C. Akcay 1, A. H. Glasser 1, C. J. Hansen 1, E. D. Held 2, T. R. Jarboe 1, J.-Y. Ji 2, C. Kim 1, W. Lowrie 1, V. S. Lukin 4,
More informationElectromagnetic Turbulence Simulations with Kinetic Electrons from the the Summit Framework
19th IAEA Fusion Energy Conference Tuesday, October 15, 2002 Paper: TH/P1-13 Electromagnetic Turbulence Simulations with Kinetic Electrons from the the Summit Framework Scott Parker and Yang Chen University
More informationIntegrated Heat Transport Simulation of High Ion Temperature Plasma of LHD
1 TH/P6-38 Integrated Heat Transport Simulation of High Ion Temperature Plasma of LHD S. Murakami 1, H. Yamaguchi 1, A. Sakai 1, K. Nagaoka 2, H. Takahashi 2, H. Nakano 2, M. Osakabe 2, K. Ida 2, M. Yoshinuma
More informationAMSC 663 Project Proposal: Upgrade to the GSP Gyrokinetic Code
AMSC 663 Project Proposal: Upgrade to the GSP Gyrokinetic Code George Wilkie (gwilkie@umd.edu) Supervisor: William Dorland (bdorland@umd.edu) October 11, 2011 Abstract Simulations of turbulent plasma in
More informationOverview of edge modeling efforts for advanced divertor configurations in NSTX-U with magnetic perturbation fields
Overview of edge modeling efforts for advanced divertor configurations in NSTX-U with magnetic perturbation fields H. Frerichs, O. Schmitz, I. Waters, G. P. Canal, T. E. Evans, Y. Feng and V. Soukhanovskii
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 informationDirect drive by cyclotron heating can explain spontaneous rotation in tokamaks
Direct drive by cyclotron heating can explain spontaneous rotation in tokamaks J. W. Van Dam and L.-J. Zheng Institute for Fusion Studies University of Texas at Austin 12th US-EU Transport Task Force Annual
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 informationOverview of Tokamak Rotation and Momentum Transport Phenomenology and Motivations
Overview of Tokamak Rotation and Momentum Transport Phenomenology and Motivations Lecture by: P.H. Diamond Notes by: C.J. Lee March 19, 2014 Abstract Toroidal rotation is a key part of the design of ITER
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 informationControl of Neo-classical tearing mode (NTM) in advanced scenarios
FIRST CHENGDU THEORY FESTIVAL Control of Neo-classical tearing mode (NTM) in advanced scenarios Zheng-Xiong Wang Dalian University of Technology (DLUT) Dalian, China Chengdu, China, 28 Aug, 2018 Outline
More informationHybrid Simulations: Numerical Details and Current Applications
Hybrid Simulations: Numerical Details and Current Applications Dietmar Krauss-Varban and numerous collaborators Space Sciences Laboratory, UC Berkeley, USA Boulder, 07/25/2008 Content 1. Heliospheric/Space
More informationLinear stability of MHD configurations
Linear stability of MHD configurations Rony Keppens Centre for mathematical Plasma Astrophysics KU Leuven Rony Keppens (KU Leuven) Linear MHD stability CHARM@ROB 2017 1 / 18 Ideal MHD configurations Interested
More informationTURBULENT TRANSPORT THEORY
ASDEX Upgrade Max-Planck-Institut für Plasmaphysik TURBULENT TRANSPORT THEORY C. Angioni GYRO, J. Candy and R.E. Waltz, GA The problem of Transport Transport is the physics subject which studies the physical
More 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 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 informationCoarse-graining the electron distribution in turbulence simulations of tokamak plasmas
Coarse-graining the electron distribution in turbulence simulations of tokamak plasmas Yang Chen and Scott E. Parker University of Colorado at Boulder Gregory Rewoldt Princeton Plasma Physics Laboratory
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 informationEFFECT OF ION CYCLOTRON HEATING ON FAST ION TRANSPORT AND PLASMA ROTATION IN TOKAMAKS
EFFECT OF ION CYCLOTRON HEATING ON FAST ION TRANSPORT AND PLASMA ROTATION IN TOKAMAKS by V.S. Chan, S.C. Chiu, and Y.A. Omelchenko Presented at the American Physical Society Division of Plasma Physics
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 informationOVERVIEW OF THE ALCATOR C-MOD PROGRAM. IAEA-FEC November, 2004 Alcator Team Presented by Martin Greenwald MIT Plasma Science & Fusion Center
OVERVIEW OF THE ALCATOR C-MOD PROGRAM IAEA-FEC November, 2004 Alcator Team Presented by Martin Greenwald MIT Plasma Science & Fusion Center OUTLINE C-Mod is compact, high field, high density, high power
More informationSupported by. Validation of a new fast ion transport model for TRANSP. M. Podestà - PPPL
NSTX-U Supported by Validation of a new fast ion transport model for TRANSP Coll of Wm & Mary Columbia U CompX General Atomics FIU INL Johns Hopkins U LANL LLNL Lodestar MIT Lehigh U Nova Photonics Old
More informationAlpha Particle Transport Induced by Alfvénic Instabilities in Proposed Burning Plasma Scenarios
Alpha Particle Transport Induced by Alfvénic Instabilities in Proposed Burning Plasma Scenarios G. Vlad, S. Briguglio, G. Fogaccia and F. Zonca Associazione Euratom-ENEA sulla Fusione, C.R. Frascati C.P.
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 informationRWM Control Code Maturity
RWM Control Code Maturity Yueqiang Liu EURATOM/CCFE Fusion Association Culham Science Centre Abingdon, Oxon OX14 3DB, UK Work partly funded by UK EPSRC and EURATOM. The views and opinions expressed do
More informationProperties of freely decaying and driven turbulence of fusion plasmas using gyrokinetic particle simulation
J. Plasma Fusion Res. SERIES, Vol. 9 () Properties of freely decaying and driven turbulence of fusion plasmas using gyrokinetic particle simulation R. Ganesh Institute for Plasma Research, Bhat Village,
More informationGTC Simulation of Turbulence and Transport in Tokamak Plasmas
GTC Simulation of Turbulence and Transport in Tokamak Plasmas Z. Lin University it of California, i Irvine, CA 92697, USA and GPS-TTBP Team Supported by SciDAC GPS-TTBP, GSEP & CPES Motivation First-principles
More information(a) (b) (c) (d) (e) (f) r (minor radius) time. time. Soft X-ray. T_e contours (ECE) r (minor radius) time time
Studies of Spherical Tori, Stellarators and Anisotropic Pressure with M3D 1 L.E. Sugiyama 1), W. Park 2), H.R. Strauss 3), S.R. Hudson 2), D. Stutman 4), X-Z. Tang 2) 1) Massachusetts Institute of Technology,
More informationVerification and validation of integrated simulation of energetic particles in fusion plasmas I: linear simulations
Page of 0 AUTHOR SUBMITTED MANUSCRIPT - NF- 0 0 0 Verification and validation of integrated simulation of energetic particles in fusion plasmas I: linear simulations S. Taimourzadeh, E. M. Bass, Y. Chen,
More informationA.G. PEETERS UNIVERSITY OF BAYREUTH
IN MEMORIAM GRIGORY PEREVERZEV A.G. PEETERS UNIVERSITY OF BAYREUTH ESF Workshop (Garching 2013) Research areas Grigory Pereverzev. Current drive in magnetized plasmas Transport (ASTRA transport code) Wave
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 informationThe Levitated Dipole Experiment: Experiment and Theory
The Levitated Dipole Experiment: Experiment and Theory Jay Kesner, R. Bergmann, A. Boxer, J. Ellsworth, P. Woskov, MIT D.T. Garnier, M.E. Mauel Columbia University Columbia University Poster CP6.00083
More informationRotation and Neoclassical Ripple Transport in ITER
Rotation and Neoclassical Ripple Transport in ITER Elizabeth J. Paul 1 Matt Landreman 1 Francesca Poli 2 Don Spong 3 Håkan Smith 4 William Dorland 1 1 University of Maryland 2 Princeton Plasma Physics
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 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 informationDIII D. by F. Turco 1. New York, January 23 rd, 2015
Modelling and Experimenting with ITER: the MHD Challenge by F. Turco 1 with J.M. Hanson 1, A.D. Turnbull 2, G.A. Navratil 1, F. Carpanese 3, C. Paz-Soldan 2, C.C. Petty 2, T.C. Luce 2, W.M. Solomon 4,
More informationAdvanced Tokamak Research in JT-60U and JT-60SA
I-07 Advanced Tokamak Research in and JT-60SA A. Isayama for the JT-60 team 18th International Toki Conference (ITC18) December 9-12, 2008 Ceratopia Toki, Toki Gifu JAPAN Contents Advanced tokamak development
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More 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 informationIntrinsic rotation due to non- Maxwellian equilibria in tokamak plasmas. Jungpyo (J.P.) Lee (Part 1) Michael Barnes (Part 2) Felix I.
Intrinsic rotation due to non- Maxwellian equilibria in tokamak plasmas Jungpyo (J.P.) Lee (Part 1) Michael Barnes (Part 2) Felix I. Parra MIT Plasma Science & Fusion Center. 1 Outlines Introduction to
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 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 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 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 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 informationSimulation study on the nonlinear EMIC waves
SH21B-2210 Simulation study on the nonlinear EMIC waves Kicheol Rha 1*, Chang-Mo Ryu 1 and Peter H Yoon 2 * lancelot@postech.ac.kr 1 Department of Physics, Pohang University of Science and Technology,
More informationELMs and Constraints on the H-Mode Pedestal:
ELMs and Constraints on the H-Mode Pedestal: A Model Based on Peeling-Ballooning Modes P.B. Snyder, 1 H.R. Wilson, 2 J.R. Ferron, 1 L.L. Lao, 1 A.W. Leonard, 1 D. Mossessian, 3 M. Murakami, 4 T.H. Osborne,
More informationPHYSICS OF HOT DENSE PLASMAS
Chapter 6 PHYSICS OF HOT DENSE PLASMAS 10 26 10 24 Solar Center Electron density (e/cm 3 ) 10 22 10 20 10 18 10 16 10 14 10 12 High pressure arcs Chromosphere Discharge plasmas Solar interior Nd (nω) laserproduced
More informationTurbulence in Tokamak Plasmas
ASDEX Upgrade Turbulence in Tokamak Plasmas basic properties and typical results B. Scott Max Planck Institut für Plasmaphysik Euratom Association D-85748 Garching, Germany Uni Innsbruck, Nov 2011 Basics
More informationExtended MHD simulation of Rayleigh-Taylor/Kelvin-Helmholtz instability
Extended MHD simulation of Rayleigh-Taylor/Kelvin-Helmholtz instability R. Goto (a), T. Hatori (a), H. Miura (a, b), A. Ito (a, b) and M. Sato (b) (a) The Graduate University for Advanced Studies (Sokendai)
More informationThe Physics of Collisionless Accretion Flows. Eliot Quataert (UC Berkeley)
The Physics of Collisionless Accretion Flows Eliot Quataert (UC Berkeley) Accretion Disks: Physical Picture Simple Consequences of Mass, Momentum, & Energy Conservation Matter Inspirals on Approximately
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 informationG. Rewoldt, W.X. Wang, M. Bell, S. Kaye, W. Solomon, R. Nazikian, and W.M. Tang Princeton Plasma Physics Lab 1
Studies with GTC-Neo, including: Recent Applications of GTC-Neo for: (1)Studies of Toroidal Angular Momentum and Ion Heat Transport, and (2)Implications for CHERS Temperature Measurements G. Rewoldt, W.X.
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 informationGlobal gyrokinetic modeling of geodesic acoustic modes and shear Alfvén instabilities in ASDEX Upgrade.
1 EX/P1-18 Global gyrokinetic modeling of geodesic acoustic modes and shear Alfvén instabilities in ASDEX Upgrade. A. Biancalani 1, A. Bottino 1, S. Briguglio 2, G.D. Conway 1, C. Di Troia 2, R. Kleiber
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 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 information