Nonlinear hydrid simulations of precessional Fishbone instability
|
|
- Easter Grant
- 5 years ago
- Views:
Transcription
1 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 2 Institute for Fusion Studies, University of Texas, USA 3 CEA, IRFM, 13115, St Paul lès Durance, France
2 Outline Introduction The Precessional Fishbone instability Tokamak geometry and particle motion A reduced Fishbone model Linear theory Analytic results Numerical results Numerical code Linear benchmark Nonlinear results Conclusions
3 Context: the Fishbone instability First observation on the PDX tokamak, with near perpendicular neutral beam injection active Bursts of electromagnetic instability, associated with losses of fast particles Mode frequency consistent with trapped fast ion precession frequency Dominant m = n = 1, kink-shaped, mode structure Observation of frequency chirping Figure : [McGuire et al., PRL, 1983]
4 Trapped particle trajectories (passing are confined too!) ω c ω bounce > ω precession possible decoupling of the three dynamics: Full kinetic, Gyrokinetic or Bounce-averaged descriptions For elegant derivation, using Lie algebra, see [Littlejohn, PS, 1982] (x, v) (X gc, v, (J c, ζ c)) ω ωc (X gc, v ; J c) (X gc, v ; J c) (α, β, (J, ζ ); J c) ω ω b (α, β; J, J c) ((J p, ζ p); J, J c))
5 The small inverse aspect ratio Tokamak ɛ = a/r 0 1 ( B = B T ˆϕ ) + ψ ϕ R B T = B 0 0 R B0 1 r R 0 cos θ & ψ = RA ϕ R 0A ϕ(r, θ, ϕ, t) with ψ eq = ψ eq(r) (concentric circular cross sections at equilibrium) Safety factor: q(r) = rb θ,eq R 0 B T where B θ,eq = R 1 0 r ψ eq = R 1 0 ψ eq i.e. magnetic winding number Note: B θ B T & q = O(1)
6 Equilibrium trajectories: bounce motion Let s take ( r ) 1 ρ c J c µ = 1 2 mv 2 gyro/b T cst E = 1 2 mv mv mv 2 + µb T (r, θ) }{{}}{{} potential well, r cst E kin along field line Circular cross section A classical pendulum :-)) ( ) v 2 = (r 2 + R0 2 q 2 (r)) θ 2 R0 2 q 2 (r) θ 2 & B T = B 0 1 r R 0 cos θ E mr = 1 0 2q2 2 θ 2 µb0r mr 3 0 q2 }{{} = g/l Trapping parameter: k 2 = e+1, e = E µb 0 R 0 2 µb 0 r particles) J = 1 2π Bounce period: T b = 2πω 1 b cos θ + µb0 mr 2 0 q2 }{{} =cst (k 2 [0, 1] for trapped mv ds = 8 π R 1/2 0 q ] mµb 0 [(k 2 1)K(k) + E(k) mr = 2π J / E = 4R 0q 0 µb 0 K(k)
7 Banana width & precessional motion P ϕ = mrv ϕ + e ψ is an exact invariant (ϕ is cyclic at the equilibrium) c v ϕ v changes along the trajectory ψ (i.e. r) must change too P ϕ = mr 0v θ=0 + e c ψeq(r) + e c ψ eqδr θ=0 = P ϕ turning point = e c ψeq(r) Banana orbit width: δr θ=0 = qr 0 ω c r v θ=0 Two possible reasons for toroidal precession: 1) Magnetic shear: forward/backward motions are not along the same line ϕ = dϕ = qdθ }{{} =0 + q δrdθ q (r)δr θ=0 ω D,1 q (r) mω cr J ω b 2) The magnitude of v is not the same during forward/backward motion δe k, θ=0 = µb 0δr θ=0 ω D,2 δv θ=0 R 0 = µb0q mω cr 0r It is possible to obtain ω D rigorously ( ω D = ω D,1 + ω D,2 ) using ω D = E J ω b J p E = J ; J p = P ϕdϕ = P ϕ = e J p c ψeq(r) ω D = c e 1 ψ eq J r
8 Deeply trapped particles Let s take v 0, thus: J = 0 and stays zero (ω ω b ), as well as v µ and J = 0 are parameters (ϕ, Pϕ ) are the natural canonical variables: 2D phase space :-) P ϕ e c ψ, i.e. P ϕ r & Pϕ r ( The Hamiltonian H = 1 P e A ) 2 2m c + µb T + eφ reduces to H(ϕ, P ϕ) = µb T (P ϕ) + eφ(ϕ, P ϕ, t) }{{}}{{} equilibrium mode The coupling is done only via the electric potential φ ϕ = H P ϕ = µb T P ϕ + eφ = µ r B T (r) P ϕ P ϕ r }{{} ω D = µb 0 q mωc R 0 r P ϕ = H φ(r, ϕ, t) = e ϕ ϕ +e r φ(r, ϕ, t) P ϕ r
9 A reduced Fishbone model: Fast particle response Take into account only deeply trapped particles (v = 0) with a single value for the magnetic moment, µ = µ Reduction of the phase space from 6D to 2D (ϕ, P ϕ ) All fast particles are contained well inside the q = 1 surface ( core region ) The dominant mode has a kink-like shape (mode numbers m = n = 1, electric potential φ/r cst well inside of the q = 1 surface) f t + [H, f ] = 0 ; H(ϕ, Pϕ) = µ B T (P ϕ) + eφ(ϕ, P ϕ) µ B T P α = ω D (P α) = µb T q(p ϕ) ω C mrr(p ϕ) ; φ = r(p φ0(t) ϕ) e iϕ r 0 Note: particles interact with a single mode. All the other ones vanish far quicker in the core region where particles are.
10 A reduced Fishbone model: Bulk plasma response Fluid description for the bulk of the plasma, neglecting the thermal pressure effects and density variations Reduced-MHD description. φ t ψ t + {φ, ψ} = 0 + {φ, φ} {ψ, ψ} = ρ = α normalization [( ˆϕ κ) P,h ] Toroidal effects are retained only for the contribution given by fast particles (κ is the toroidal curvature and P,h the fast particle pressure) Cylindrical geometry As before we set φ = r φ 0(t) r 0 e iϕ iθ well inside the q = 1 surface On the contrary strong variations are allowed across the q = 1 surface and finally φ 0 for r a One mode evolution: only Fast particle (kinetic) nonlinearities are retained, MHD nonlinearities are neglected
11 Linear theory: Analytic results I Let s take f = F eq + δf, with δf F eq, and t iω. The mode equation reads ω 2 ( φ r ) + v 2 A,T R 2 0 ( 1 1 ) 2 ( ) φ = iω 1 r d r r 2 ρ(δf ) q(r) r r 3 0 where δf = er φ 0 df eq r 0 dr dr dp ϕ ω D (r) ω Note as a spatial gradient corresponds to the usual velocity gradient. Here a decreasing density is equivalent to a bump on the tail. Finally a general dispersion relation is obtained : ( ) 1 y 2 q(y) 2 dfeq dy i = K dy (1) q(y) 0 yω ω D (y=1) where K takes into account the energetic content for the fast particles, the MHD and geometric parameters and y = r/r (r being so that q(r ) = 1).
12 Linear theory: analytic results II Let s take: F eq = n 0 2 δ(µ µ )δv (1 erf(β(y y 0)) analytic values the threshold condition and mode frequency at the threshold K 0 = 1 ( ; ω πβy 3 0 = ω D (y 0) 1 1 ) 0 β 2 y0 2 where y 0 is the position of the highest radial gradient in the distribution function Close to the threshold K = K 0 + δk ; ω = ω 0 + δω + iγ we obtain a growth rate and a correction for the real frequency: γ δk π ω 0 ; δω 1 π γ K 0 2βy 0 2 βy 0
13 Nonlinear numerical code Based on the same domain decomposition: 1) The core region I Nonlinear kinetic description: Semi-lagrangian code assuming φ = r φr bound I Linearized MHD response, including the fast particle pressure: i.e. an evolution equation for r φr bound plus the frozen-in equation for ψ bound 2) Thin annular region around q = 1 surface I No fast particle here I Cylindrical slab description I Semi-spectral MHD code (only one mode at present time) φ & ψ bound as B.C. r bound I Uses r I Provides ψ & φ, in particular φ bound (the Hamiltonian)
14 γ Linear benchmark With numerical simulations, the linear results are recovered Good agreement with frequency value Good agreement with growth rate and mode shape 2.5 x x φ K/K x Figure : Growth rate as a function of K/K 0. Stars are the numerical values, the green line is the analytic prediction. Figure : φ profile in the annular layer. In blue the real part, in red the imaginary part.
15 Nonlinear results: Mode saturation level The first local maximum of the amplitude is proportional to γ 2 i.e when the phase-space island width φ reaches the resonance width γ [Zonca et al., NJP, 2015] Amplitude oscillations are far larger compared to the usual Bump on Tail case, with [Berk et al., PLA, 1997] or without dissipation [O Neil, PF, 1965] x 10 5 Figure : Evolution of the kinetic energy of the mode versus time, in logscale
16 Frequency chirping and particle ejection Chirping is observed during the saturated phase (case studied here: K/K 0 = 1.2, giving γ/ω = 2.8%) Phase space structure motion matches frequency change Asymmetric system Higher amplitude for the down chirping mode, i.e. particle ejection Outgoing particles continue to interact with the mode but actually are not trapped into the mode well Figure : On the left : Evolution of the distribution function averaged over angle. On the right : spectrogram of the mode.
17 e e e x e e Structures in phase-space Dynamics close to the first maximum and minimum of the mode amplitude Only partial folding of f inside the phase-space island Strong stretching in ϕ & P ϕ directions Figure : Distribution function at different stages of the saturation.
18 Island contraction and slippage
19 Contraction and slippage: The Fishbone peculiarity The dispersion relation of the Fishbone is really different from the usual Bump on Tail one: BoT: marginally stable plasma wave + particle driver leading to instability i.e. RD(ω, k) 0 & a small ID gives the instability Thus when the energetic driver drops, the mode response is not dramatic. Fishbone exists only because fast particles are there, indeed ω R = ω D,h It is a genuine energetic particle mode [Zonca et al., NJP, 2015] The mode response to the energetic driver variation can be strong: 2 t 2 ( ) φ + ωa(r) 2 r ( ) φ = i r t 1 r d r r 2 ρ(f ) = RHS r 3 0 In the core ω 2 A,bound tt ω 2 R thus the mode is slave, compared to the particle driver: ( ) φ r bound RHS ωa 2 bound giving φ bound = φ 0 RHS Amplitude and phase simply follow: the mode is slave to the energetic driver.
20 Conclusions The Precessional Fishbone instability can be described by an hybrid Fluid-Hamiltonian model where the phase-space coordinates are quite unusual. This reduced model is able to catch the qualitative dynamics of the mode: the frequency downchirping of the mode and the gradual ejection of fast particles. It has the great advantage, over more complete models, of permitting an easier analysis of the structure dynamics in phase-space. Future perspectives: 1) Allowing MHD nonlinearities to develops around the q = 1 surface. 2) Looking at a more general approach for energetic particle modes.
Gyrokinetic 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 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 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 informationModelling of Frequency Sweeping with the HAGIS code
Modelling of Frequency Sweeping with the HAGIS code S.D.Pinches 1 H.L.Berk 2, S.E.Sharapov 3, M.Gryaznavich 3 1 Max-Planck-Institut für Plasmaphysik, EURATOM Assoziation, Garching, Germany 2 Institute
More informationEffects of drag and diffusion on nonlinear behavior of EP-driven instabilities.
IAEA-TM EP 2011 / 09 / 07 Effects of drag and diffusion on nonlinear behavior of EP-driven instabilities. Maxime Lesur Y. Idomura, X. Garbet, P. Diamond, Y. Todo, K. Shinohara, F. Zonca, S. Pinches, M.
More informationInteraction between EGAMs and turbulence in full-f gyrokinetic simulations
Interaction between EGAMs and turbulence in full-f gyrokinetic simulations David Zarzoso 1 X Garbet 1, Y Sarazin 1, V Grandgirard 1, J Abiteboul 1, A Strugarek 1,2, G Dif-Pradalier 1, R Dumont 1, G Latu
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 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 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 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 informationGlobal gyrokinetic particle simulations with kinetic electrons
IOP PUBLISHING Plasma Phys. Control. Fusion 49 (2007) B163 B172 PLASMA PHYSICS AND CONTROLLED FUSION doi:10.1088/0741-3335/49/12b/s15 Global gyrokinetic particle simulations with kinetic electrons Z Lin,
More informationNonlinear Consequences of Weakly Driven Energetic Particle Instabilities
2008 International Sherwood Fusion Theory Conference March 30 - April 2, 2008, Boulder, Colorado Nonlinear Consequences of Weakly Driven Energetic Particle Instabilities Boris Breizman Institute for Fusion
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 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 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 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 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 informationMHD instabilities and fast particles
ENEA F. Zonca 1 MHD instabilities and fast particles Fulvio Zonca Associazione Euratom-ENEA sulla Fusione, C.R. Frascati, C.P. 65-44 - Frascati, Italy. July 11.th, 25 : Turbulence overshoot and resonant
More informationSingle Particle Motion
Single Particle Motion C ontents Uniform E and B E = - guiding centers Definition of guiding center E gravitation Non Uniform B 'grad B' drift, B B Curvature drift Grad -B drift, B B invariance of µ. Magnetic
More informationEnergetic particle modes: from bump on tail to tokamak plasmas
Energetic particle modes: from bump on tail to tokamak plasmas M. K. Lilley 1 B. N. Breizman 2, S. E. Sharapov 3, S. D. Pinches 3 1 Physics Department, Imperial College London, London, SW7 2AZ, UK 2 IFS,
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 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 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 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 informationLectures on basic plasma physics: Hamiltonian mechanics of charged particle motion
Lectures on basic plasma physics: Hamiltonian mechanics of charged particle motion Department of applied physics, Aalto University March 8, 2016 Hamiltonian versus Newtonian mechanics Newtonian mechanics:
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 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 informationA theory for localized low-frequency ideal MHD modes in axisymmetric toroidal systems is generalized to take into account both toroidal and poloidal
MHD spectra pre-history (selected results I MHD spectra pre-history (selected results II Abstract A theory for localized low-frequency ideal MHD modes in axisymmetric toroidal systems is generalized to
More informationProgress and Plans on Physics and Validation
Progress and Plans on Physics and Validation T.S. Hahm Princeton Plasma Physics Laboratory Princeton, New Jersey Momentum Transport Studies: Turbulence and Neoclassical Physics Role of Trapped Electrons
More informationPFC/JA-93-5 ANOMALOUS ELECTRON STREAMING DUE TO WAVES IN TOKAMAK PLASMAS. April S. D. Schultz, A. Bers, and A. K. Ram
PFC/JA-93-5 ANOMALOUS ELECTRON STREAMING DUE TO WAVES IN TOKAMAK PLASMAS S. D. Schultz, A. Bers, and A. K. Ram April 1993 Plasma Fusion Center Massachusetts Institute of Technology Cambridge, Massachusetts
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 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 informationPredictions of fusion α-particle transport due to Alfvén eigenmodes in ITER
Predictions of fusion α-particle transport due to Alfvén eigenmodes in ITER M. Fitzgerald, S.E. Sharapov, P. Rodrigues 2, A. Polevoi 3, D. Borba 2 2 Instituto de Plasmas e Fusão Nuclear, Instituto Superior
More informationGyrokinetic simulations with GYSELA: Main current issues in physics & numerics
Gyrokinetic simulations with GYSELA: Main current issues in physics & numerics Y. Sarazin, Y. Asahi 2, N. Bouzat, G. Dif-Pradalier, P. Donnel, C. Ehrlacher, C. Emeriau 3, X. Garbet, Ph. Ghendrih, V. Grandgirard,
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 informationIon plateau transport near the tokamak magnetic axis. Abstract
Ion plateau transport near the tokamak magnetic axis K.C. Shaing and R.D. Hazeltine Institute for Fusion Studies, The University of Texas at Austin, Austin, Texas 78712 (December 2, 1997) Abstract Conventional
More informationSpectroscopic determination of the internal amplitude of frequency sweeping TAE
INSTITUTE OF PHYSICS PUBLISHING Plasma Phys. Control. Fusion 46 (2004) S47 S57 PLASMA PHYSICS AND CONTROLLED FUSION PII: S0741-3335(04)72680-9 Spectroscopic determination of the internal amplitude of frequency
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 informationBursty Transport in Tokamaks with Internal Transport Barriers
Bursty Transport in Tokamaks with Internal Transport Barriers S. Benkadda 1), O. Agullo 1), P. Beyer 1), N. Bian 1), P. H. Diamond 3), C. Figarella 1), X. Garbet 2), P. Ghendrih 2), V. Grandgirard 1),
More informationNonlinear magnetohydrodynamical effects in precessional fishbone oscillations
Nonlinear magnetohydrodynamical effects in precessional fishbone oscillations A. Ödblom, B. N. Breizman, S. E. Sharapov, T. C. Hender, and V. P. Pastukhov Citation: Phys. Plasmas 9, 155 (22); doi: 1.163/1.1421373
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 informationb c. d a e Λ h α
MHD Phenomena and Transport of Energetic Ions in Spherical Tori Ya.I. Kolesnichenko ), V.V. Lutsenko ), V.S. Marchenko ), R.B. White ), Yu.V. Yakovenko ) ) Institute for Nuclear Research, National Academy
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 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 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 informationGeometric Gyrokinetic Theory and its Applications to Large-Scale Simulations of Magnetized Plasmas
Geometric Gyrokinetic Theory and its Applications to Large-Scale Simulations of Magnetized Plasmas Hong Qin Princeton Plasma Physics Laboratory, Princeton University CEA-EDF-INRIA School -- Numerical models
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 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 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 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 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 informationOverview of Gyrokinetic Theory & Properties of ITG/TEM Instabilities
Overview of Gyrokinetic Theory & Properties of ITG/TEM Instabilities G. W. Hammett Princeton Plasma Physics Lab (PPPL) http://w3.pppl.gov/~hammett AST559: Plasma & Fluid Turbulence Dec. 5, 2011 (based
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 informationPhysical Dynamics (SPA5304) Lecture Plan 2018
Physical Dynamics (SPA5304) Lecture Plan 2018 The numbers on the left margin are approximate lecture numbers. Items in gray are not covered this year 1 Advanced Review of Newtonian Mechanics 1.1 One Particle
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 informationMotion of Charged Particles in Fields
Chapter Motion of Charged Particles in Fields Plasmas are complicated because motions of electrons and ions are determined by the electric and magnetic fields but also change the fields by the currents
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 informationControlling chaotic transport in Hamiltonian systems
Controlling chaotic transport in Hamiltonian systems Guido Ciraolo Facoltà di Ingegneria, Università di Firenze via S. Marta, I-50129 Firenze, Italy Cristel Chandre, Ricardo Lima, Michel Vittot CPT-CNRS,
More informationParticle Pinch Model of Passing/Trapped High-Z Impurity with Centrifugal Force Effect )
Particle Pinch Model of Passing/Trapped High-Z Impurity with Centrifugal Force Effect ) Yusuke SHIMIZU, Takaaki FUJITA, Atsushi OKAMOTO, Hideki ARIMOTO, Nobuhiko HAYASHI 1), Kazuo HOSHINO 2), Tomohide
More informationModels for Global Plasma Dynamics
Models for Global Plasma Dynamics F.L. Waelbroeck Institute for Fusion Studies, The University of Texas at Austin International ITER Summer School June 2010 Outline 1 Models for Long-Wavelength Plasma
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 informationRF cavities (Lecture 25)
RF cavities (Lecture 25 February 2, 2016 319/441 Lecture outline A good conductor has a property to guide and trap electromagnetic field in a confined region. In this lecture we will consider an example
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 informationLecture 5. Frequency sweeping, phase locking and supra-thermal particle transport:
IFTS Intensive Course on Advanced Plasma Physics-Spring 212 Lecture 5 1 Lecture 5 Frequency sweeping, phase locking and supra-thermal particle transport: heuristic models for understanding nonlinear behaviors
More informationSummer College on Plasma Physics August Introduction to Nonlinear Gyrokinetic Theory
2052-24 Summer College on Plasma Physics 10-28 August 2009 Introduction to Nonlinear Gyrokinetic Theory T.S. Hahm Princeton Plasma Physics Laboratory Princeton University USA Introduction to Nonlinear
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 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 informationPresentation by Herb Berk University of Texas at Austin Institute for Fusion Studies in Vienna, Austria Sept. 1-4, 2015
Review of Theory Papers at 14 th IAEA technical meeting on Engertic Particles in Magnetic Confinement systems Presentation by Herb Berk University of Texas at Austin Institute for Fusion Studies in Vienna,
More informationMAGNETOHYDRODYNAMICS
Chapter 6 MAGNETOHYDRODYNAMICS 6.1 Introduction Magnetohydrodynamics is a branch of plasma physics dealing with dc or low frequency effects in fully ionized magnetized plasma. In this chapter we will study
More informationHamiltonian formulation of reduced Vlasov-Maxwell equations
Hamiltonian formulation of reduced Vlasov-Maxell equations Cristel CHANDRE Centre de Physique Théorique CNRS, Marseille, France Contact: chandre@cpt.univ-mrs.fr importance of stability vs instability in
More informationcos 6 λ m sin 2 λ m Mirror Point latitude Equatorial Pitch Angle Figure 5.1: Mirror point latitude as function of equatorial pitch angle.
Chapter 5 The Inner Magnetosphere 5.1 Trapped Particles The motion of trapped particles in the inner magnetosphere is a combination of gyro motion, bounce motion, and gradient and curvature drifts. In
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 informationWhat place for mathematicians in plasma physics
What place for mathematicians in plasma physics Eric Sonnendrücker IRMA Université Louis Pasteur, Strasbourg projet CALVI INRIA Nancy Grand Est 15-19 September 2008 Eric Sonnendrücker (U. Strasbourg) Math
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 informationLarge scale flows and coherent structure phenomena in flute turbulence
Large scale flows and coherent structure phenomena in flute turbulence I. Sandberg 1, Zh. N. Andrushcheno, V. P. Pavleno 1 National Technical University of Athens, Association Euratom Hellenic Republic,
More informationSingle particle motion and trapped particles
Single particle motion and trapped particles Gyromotion of ions and electrons Drifts in electric fields Inhomogeneous magnetic fields Magnetic and general drift motions Trapped magnetospheric particles
More informationAdvances in stellarator gyrokinetics
Advances in stellarator gyrokinetics Per Helander and T. Bird, F. Jenko, R. Kleiber, G.G. Plunk, J.H.E. Proll, J. Riemann, P. Xanthopoulos 1 Background Wendelstein 7-X will start experiments in 2015 optimised
More informationSpectroscopic Determination of the Internal Amplitude of Frequency Sweeping TAE
EFDA JET PR(3)58 S.D. Pinches, H.L. Berk, M.P. Gryaznevich, S.E. Sharapov and JET EFDA contributors Spectroscopic Determination of the Internal Amplitude of Frequency Sweeping TAE . Spectroscopic Determination
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 informationPhysical Dynamics (PHY-304)
Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.
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 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 informationInternal transport barriers in some Hamiltonian systems modeling the magnetic lines dynamics in tokamak
Journal of Physics: Conference Series Internal transport barriers in some Hamiltonian systems modeling the magnetic lines dynamics in tokamak To cite this article: D Constantinescu et al 2005 J. Phys.:
More informationIntroduction to Plasma Physics
Introduction to Plasma Physics Hartmut Zohm Max-Planck-Institut für Plasmaphysik 85748 Garching DPG Advanced Physics School The Physics of ITER Bad Honnef, 22.09.2014 A simplistic view on a Fusion Power
More informationStability of a plasma confined in a dipole field
PHYSICS OF PLASMAS VOLUME 5, NUMBER 10 OCTOBER 1998 Stability of a plasma confined in a dipole field Plasma Fusion Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 Received
More informationNonlinear Alfvén Wave Physics in Fusion Plasmas
ASIPP 40 th Anniversary Nonlinear Alfvén Wave Physics in Fusion Plasmas 1 Nonlinear Alfvén Wave Physics in Fusion Plasmas Institute for Fusion Theory and Simulation, Zhejiang University, Hangzhou 310027,
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 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 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 informationNeoclassical transport
Neoclassical transport Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 28 th January 2013 Dr Ben Dudson Magnetic Confinement Fusion (1 of 19) Last time Toroidal devices
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationModeling of Transport Barrier Based on Drift Alfvén Ballooning Mode Transport Model
9th IAEA TM on H-mode Physics and Transport Barriers Catamaran Resort Hotel, San Diego 3-9-5 Modeling of Transport Barrier Based on Drift Alfvén Ballooning Mode Transport Model A. Fukuyama, M. Uchida and
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 informationEdge Rotational Shear Requirements for the Edge Harmonic Oscillation in DIII D Quiescent H mode Plasmas
Edge Rotational Shear Requirements for the Edge Harmonic Oscillation in DIII D Quiescent H mode Plasmas by T.M. Wilks 1 with A. Garofalo 2, K.H. Burrell 2, Xi. Chen 2, P.H. Diamond 3, Z.B. Guo 3, X. Xu
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 informationCharged particle motion in external fields
Chapter 2 Charged particle motion in external fields A (fully ionized) plasma contains a very large number of particles. In general, their motion can only be studied statistically, taking appropriate averages.
More informationTheory and Simulation of Neoclassical Transport Processes, with Local Trapping
Theory and Simulation of Neoclassical Transport Processes, with Local Trapping Daniel H. E. Dubin Department of Physics, University of California at San Diego, La Jolla, CA USA 92093-0319 Abstract. Neoclassical
More informationMicroturbulence in optimised stellarators
Q Josefine H. E. Proll, Benjamin J. Faber, Per Helander, Samuel A. Lazerson, Harry Mynick, and Pavlos Xanthopoulos Many thanks to: T. M. Bird, J. W. Connor, T. Go rler, W. Guttenfelder, G.W. Hammett, F.
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 informationChapter 5 MAGNETIZED PLASMAS. 5.1 Introduction. 5.2 Diamagnetic current
Chapter 5 MAGNETIZED PLASMAS 5.1 Introduction We are now in a position to study the behaviour of plasma in a magnetic field. In the first instance we will re-examine particle diffusion and mobility with
More information