The gyrokinetic turbulence code GENE - Numerics and applications
|
|
- Darleen Rice
- 5 years ago
- Views:
Transcription
1 Contributors: T. Dannert (1), F. Jenko (1),F. Merz (1), D. Told (1), X. Lapillonne (2), S. Brunner (2), and others T. Görler (1) The gyrokinetic turbulence code GENE - Numerics and applications (1) Max-Planck-Institut für Plasmaphysik, Garching (2) CRPP/EPFL, Lausanne 3 rd summer school of the Large Scale Initiative FUSION September 26 th -30 th, 2011 Paris
2 What to expect from this lecture? A brief motivation for and introduction to gyrokinetic theory A discussion of a massively parallelized implementation - the plasma microturbulence code GENE Examples of simulation results Summary & Outlook
3 Plasma turbulence modeling and gyrokinetic theory in brief
4 Plasma microturbulence Energy confinement time set by anomalous heat and particle transport one of the key physics problems Plasma density Plasma temperature Ion temperature (kev) Commonly attributed to plasma microturbulence ρ Microinstabilities driven by strong temperature and density background gradients Quasistationary state far from thermodynamic equilibrium
5 Plasma turbulence: Kinetic theory! Dilute and/or hot plasmas are almost collisionless. Thus, if kinetic effects (finite Larmor radius, Landau damping, magnetic trapping etc.) play a role, magnetohydrodynamics is not applicable, and one has to use a kinetic description! Vlasov-Maxwell equations Removing the fast gyromotion leads to a dramatic speed-up Charged rings as quasiparticles; gyrocenter coordinates; keep kinetic effects Brizard & Hahm, Rev. Mod. Phys. 79, 421 (2007)
6 The gyrokinetic ordering The gyrokinetic model is a Vlasov-Maxwell on which the GK ordering is imposed Slow time variation as compared to the gyromotion time scale Spatial equilibrium scale much larger than the Larmor radius Strong anisotropy, i.e. only perpendicular gradients of the fluctuating quantities can be large Small amplitude perturbations, i.e. energy of the perturbation much smaller than the thermal energy:
7 Theoretical framework: Gyrokinetic theory Gyrokinetic Vlasov equation Poisson equation with gyrocenter position Ampère s law parallel velocity and magnetic moment Nonlinear, 5-dimensional, partial integro-differential system of equations to be solved!
8 Major theoretical speedups G.Hammett relative to original Vlasov/pre-Maxwell system on a naïve grid, for ITER * = /a ~ 1/1000 Nonlinear gyrokinetic equations eliminate plasma frequency: pe / i ~ m i /m e x10 3 eliminate Debye length scale: ( i / De ) 3 ~ (m i /m e ) 3/2 x10 5 average over fast ion gyration: i / ~ 1/ * x10 3 Field-aligned coordinates adapt to elongated structure of turbulent eddies: / ~ 1/ * x10 3 Reduced simulation volume reduce toroidal mode numbers (i.e., 1/15 of toroidal direction) L r ~ a/6 ~ 160 r ~ 10 correlation lengths x15 x6 Total speedup x10 16 For comparison: Massively parallel computers ( ) x10 7
9 Status quo in gyrokinetic simulation over the last decade or so, GK has emerged as the standard approach to plasma turbulence a variety of nonlinear GK codes is being used and (further) developed these codes differ with respect to their numerics, physics, parallel scalability, and public availability GENE GYSELA ORB5 ELMFIRE... Europe US GYRO GTC GEM GS2... GKV G5D... Japan
10 The gyrokinetic Vlasov code GENE
11 The gyrokinetic code GENE (gene.rzg.mpg.de) GENE is a physically comprehensive Vlasov code: allows for kinetic electrons & electromagnetic fluctuations, collisions, and external ExB shear flows is coupled to various MHD codes and the transport code TRINITY can be used as initial value or eigenvalue solver supports local (flux-tube) and global (full-torus), gradient- and fluxdriven simulations GENE is well benchmarked and hyperscalable Strong scaling on BG/P 11
12 Concepts used within GENE for speed-up field-aligned (Clebsch-type) coordinate system to exploit the high anisotropy of plasma turbulence; parallel (z) derivatives can be taken small compared to perpendicular (x,y) ones(!) f-splitting: Apply same approach as in the derivation of the GKE and split the distribution function f = F 0 + f where F 0 : stationary background, here: local Maxwellian f: fluctuating part with f/f 0 << 1 Lowest-order nonlinearity kept next order can be switched on for testing (electrostatic version) 12
13 Numerical methods time scheme I Method of Lines: turn PDEs into ODEs by discretizing the spatial derivative first solve for the continuous time coordinate Time Solver: Full eigenvalue spectrum of the GK lin. operator Linear system: Iterative eigenvalue solver based on PETSc/SLEPc/Scalapack lib s solve for largest abs/re/im eigenvalues gain insights in linear stability/physics Explicit Runge-Kutta (ERK) schemes 13
14 Numerical methods time scheme II Time Solver: (Full) Nonlinear system: Linear stability regions for low-order RK-schemes Several ERK-methods, e.g. 4 th order Optimum linear time step can be precomputed using iterative EV-Solver Adaptive CFL time step adaption for nonlinearity 14
15 Phase space discretization GENE is a Eulerian code; thus solving the 5D (δf-splitted) distribution function on a fixed grid in (x, v, ) trapped/passing boundary radial direction x: toroidal direction y: parallel direction z: v -velocity space: -velocity space: equidistant grid (either configuration or Fourier space) equidistant grid in Fourier space equidistant grid points equidistant grid Gauss-Legendre or Gauss-Laguerre knots 15
16 Phase space discretization finite differences 4 th order finite differences are basic choice however, the more elaborate Arakawa-scheme is employed, if possible [A. Arakawa, JCP 135, 103 (1997), reprint] Exception: spectral methods are employed in the y direction - always x direction - depends on the type of operation (local/global) 16
17 Local vs. global GENE Local: in the radial direction Simulation domain small compared to machine size; thus, constant temperatures/densities and fixed gradients a Periodic boundary conditions; allows application of spectral methods Local sim. domain ρ*=ρ s /a << 1 Global sim. domain Global: adding nonlocal features in the radial direction Consider full temperature & density profiles; radially varying metric Dirichlet or v. Neumann boundary conditions Heat sources & sinks 17
18 Local vs. global GENE numerical point of view Local approach: Spectral methods derivatives: Gyroaverage & field operators can be given analytically: T(x) / T 0 Global approach: derivatives: finite difference scheme, typically 4 th order Use interpolation schemes for gyroaverage & field operators with gyromatrix
19 Gyro-averaging procedure in more detail discretize gyro-angle integration coordinate transform interpolation between grid-cells required (!) here: 1-dimensional problem (y remains in Fourier space) use finite-elements which allows easy extraction of gyroaveraged quantities at original grid points (~ Hermite polynomial interpol.) Gyromatrix is constructed at initialization only:
20 Boundary conditions - overview radial direction x: local: periodic global: Dirichlet v. Neumann T T r toroidal direction y: periodic r [Beer, PoP 95] parallel direction z: quasi-periodic v -velocity space: Dirichlet b.c. magnetic shear tilts the simulation box phase factor -velocity space: Not required (if collisions are neglected)
21 Sources/sink models in GENE Full domain: Radially dependent heat source/sink in whole domain (density and parallel momentum are unaffected) Similar to ORB5 [McMillan et al., PoP `08] Background temperature gradient is approx. maintained while avoiding direct influence on microturbulence (similar operator exists for density) Localized sources: Krook operator, boundaries, to damp fluctuations close to the Heat source model (cmp. GYSELA, [Grandgirard et al., PPCF 07]) with Example of a heat source profile: and floating boundary conditions 21
22 Geometry arbitrary flux-surface shapes can be considered as long as the metric is provided to compute, e.g., standard choice (implemented): circular concentric flux surfaces others can be read via interfaces to TRACER/GIST (field line tracer) and the equilibrium code CHEASE
23 Code parallelization and optimization MPI (message passing interface) parallelization available in all phase space directions + species; however, difficult to estimate interference, especially on large numbers of processors (>1k) automatic comparison at code initialization For better utilization of different architectures (different cache sizes!), strip mining techniques and message block adaption FFTW-like performance optimization during initialization
24 Applications
25 Why global? Cover a larger radial domain (instead of using several flux tubes) Check validity of local simulations: When do meso-/large scale events, i.e. avalanches or turbulence spreading, occur? Do they affect the transport scaling? machine-size events: Bohm scaling Gyroradius scale turbulence: Gyro-Bohm scaling Re-assess earlier results by [Z. Lin et al., PRL, 2002] and [Candy et al., PoP 2004] Allow for flux-driven simulations 25
26 Local limit test (ρ* 0) with L x /a =const. for 2 profile types local value Here: CBC-like parameters with kinetic electrons at k y ρ s ~ 0.3 The global results do approach the local limit
27 Local limit test (ρ* 0) with L x /a =const. for 2 profile types local value The profile shape (linear driving region) matters! Data points align with
28 Nonlinear investigation of finite size effects Candy et al., PoP 11 (2004) GENE-local using the circular, concentric equilibrium GENE-global ORB5 GENE-local using the s-α equilibrium McMillan et al., PRL 105, (2010) ORB5 (Lagrangian) and GENE (Eulerian) agree if the same geometry model is used long lasting controversy probably resolved Both, GENE and ORB5 converge towards the local limit Deviations (global/local) < 10% at ρ* < 1/300 28
29 Finite system size: Profile shape matters logarithmic temperature gradient Local (flux-tube) limit PRL 105, (2010) Both codes also show that it is the parameter which really matters this should be kept in mind when dealing, e.g. with Internal Transport Barriers Scaling cannot be explained with profile shearing (only weak dependence) Turbulence spreading, avalanches? 29
30 Can avalanches break gyro-bohm scaling? Global flux-driven simulations of ITG-ae turbulence with GENE (for ρ*=1/140) time flux-surface averaged radial heat flux associated ExB shearing rate (radial shear layers) radius radius avalanches are mesoscale ; radial extent ~20-40 ρ i their propagation speed is found to be ~ρ*v ti propagation direction is correlated with sign(ω E ) importance of low-frequency zonal flows & mean flows Same phenomenology as, e.g., McMillan (PoP 2009) 30
31 Heat flux for global, gradient-driven ITG-ae Radial extent and propagation velocities do not depend much on ρ* Δt = 100 R0/cs ρ*=1/140 ρ*=1/560 ρ*=1/280 Δt = 100 R0/cs ρ*=1/1000 (ITER-like) 31
32 Interesting insights from local simulations ITG turbulence (adiabatic electrons) ITG turbulence (kinetic electrons) Δt = 100 R 0 /c s Avalanches are not inherently nonlocal ETG turbulence (adiabatic ions) Avalanches are less pronounced in more complete model TEM turbulence Δt = 100 R 0 /c s Avalanches are not to be identified with streamers Avalanches tend to be absent
33 Conclusions
34 Summary & Outlook The nonlinear gyrokinetic plasma turbulence code GENE has been introduced and its features and recent extension to a nonlocal code have been discussed; Examples for applications of this code suite have been presented: Transition from nonlocal to local turbulence (ρ* 0) has been revisited cooperatively via Lagrangian & Eulerian codes; linear driving region important Heat flux avalanches seem to be mesoscale (ρ i -related) phenomena and are found not to break the gyro-bohm scaling Next applications: Take advantage of the comprehensive physics model (electromagnetic effects, collisions, etc.) and address specific experimental discharges, e.g. eitb in TCV consider non-axisymmetric stellarator equilibria keeping finite-size effects References: ( publications) [Görler et al., JCP 230, 7053 (2011)]
Gyrokinetics an efficient framework for studying turbulence and reconnection in magnetized plasmas
Frank Jenko Gyrokinetics an efficient framework for studying turbulence and reconnection in magnetized plasmas Max-Planck-Institut für Plasmaphysik, Garching Workshop on Vlasov-Maxwell Kinetics WPI, Vienna,
More informationTowards Multiscale Gyrokinetic Simulations of ITER-like Plasmas
Frank Jenko Max-Planck-Institut für Plasmaphysik, Garching Universität Ulm Towards Multiscale Gyrokinetic Simulations of ITER-like Plasmas 23 rd IAEA Fusion Energy Conference 11-16 October 2010, Daejeon,
More informationGyrokinetic Turbulence in Tokamaks and Stellarators
Gyrokinetic Turbulence in Tokamaks and Stellarators Frank Jenko IPP, Germany Acknowledgements: P. Xanthopoulos, F. Merz, T. Görler, M. Pueschel, D. Told; A. Boozer, G. Hammett, D. Mikkelsen, M. Zarnstorff,
More informationCo-existence and interference of multiple modes in plasma turbulence: Some recent GENE results
Co-existence and interference of multiple modes in plasma turbulence: Some recent GENE results Frank Jenko IPP Garching, Germany University of Ulm, Germany Acknowledgements: F. Merz, T. Görler, D. Told,
More informationInternational Workshop on the Frontiers of Modern Plasma Physics July On the Nature of Plasma Core Turbulence.
1953-43 International Workshop on the Frontiers of Modern Plasma Physics 14-25 July 2008 On the Nature of Plasma Core Turbulence. F. Jenko Max-Planck Institute fuer Plasmaphysik Garching bei Munchen Germany
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 informationFrank Jenko. Session on Plasma Physics: Petascale Computations of Turbulence in Laboratory and Astrophysical Plasmas
Frank Jenko www.ipp.mpg.de/~fsj Session on Plasma Physics: Petascale Computations of Turbulence in Laboratory and Astrophysical Plasmas Max-Planck-Institut für Plasmaphysik, Germany 3 rd EU-US HPC Summer
More informationGlobal Nonlinear Simulations of Ion and Electron Turbulence Usintg a Particle-In-Cell Approach
Global Nonlinear Simulations of Ion and Electron Turbulence Usintg a Particle-In-Cell Approach S. Jolliet 1), B. F. McMillan 1), T. M. Tran 1), X. Lapillonne 1), L. Villard 1), A. Bottino 2), P. Angelino
More informationEntropy evolution and dissipation in collisionless particle-in-cell gyrokinetic simulations
Max-Planck-Insititut für Plasmaphysik Entropy evolution and dissipation in collisionless particle-in-cell gyrokinetic simulations A. Bottino Objectives Develop a numerical tool able to reproduce and predict
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 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 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 informationarxiv: v1 [physics.plasm-ph] 3 Apr 2011
A comparison of Vlasov with drift kinetic and gyrokinetic theories arxiv:1104.0427v1 [physics.plasm-ph] 3 Apr 2011 H. Tasso 1, G. N. Throumoulopoulos 2 1 Max-Planck-Institut für Plasmaphysik, Euratom Association,
More informationMultiscale, multiphysics modeling of turbulent transport and heating in collisionless, magnetized plasmas
Multiscale, multiphysics modeling of turbulent transport and heating in collisionless, magnetized plasmas Michael Barnes Plasma Science & Fusion Center Massachusetts Institute of Technology Collaborators:
More informationHPC Simulations of Magnetic Fusion Plasmas
HPC Simulations of Magnetic Fusion Plasmas L. Villard S. Brunner, A. Casati, W.A. Cooper, J.P. Graves, M. Jucker, S. Khoshagdam, X. Lapillonne, B.F. McMillan, P. Popovich, O. Sauter, T.M. Tran, T. Vernay
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 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 informationParticle-in-cell simulations of electron transport from plasma turbulence: recent progress in gyrokinetic particle simulations of turbulent plasmas
Institute of Physics Publishing Journal of Physics: Conference Series 16 (25 16 24 doi:1.188/1742-6596/16/1/2 SciDAC 25 Particle-in-cell simulations of electron transport from plasma turbulence: recent
More informationGyrokinetic Simulations of Tokamak Microturbulence
Gyrokinetic Simulations of Tokamak Microturbulence W Dorland, Imperial College, London With key contributions from: S C Cowley F Jenko G W Hammett D Mikkelsen B N Rogers C Bourdelle W M Nevins D W Ross
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 including the centrifugal force in a strongly rotating tokamak plasma
Gyrokinetic simulations including the centrifugal force in a strongly rotating tokamak plasma F.J. Casson, A.G. Peeters, Y. Camenen, W.A. Hornsby, A.P. Snodin, D. Strintzi, G.Szepesi CCFE Turbsim, July
More informationResearch on Structure Formation of Plasmas Dominated by Multiple Hierarchical Dynamics
Chapter 4 Epoch Making Simulation Research on Structure Formation of Plasmas Dominated by Multiple Hierarchical Dynamics Group Representative Yasuaki Kishimoto Naka Fusion Research Establishment, Japan
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 informationGyrokinetic Turbulence Investigations Involving Ion and Electron Scales
Gyrokinetic Turbulence Investigations Involving Ion and Electron Scales T. Görler, F. Jenko, M.J. Pueschel, D. Told, and H. Lesch Abstract Plasma microinstabilities are one of the key physics problems
More informationSize Scaling and Nondiffusive Features of Electron Heat Transport in Multi-Scale Turbulence
Size Scaling and Nondiffusive Features of Electron Heat Transport in Multi-Scale Turbulence Z. Lin 1, Y. Xiao 1, W. J. Deng 1, I. Holod 1, C. Kamath, S. Klasky 3, Z. X. Wang 1, and H. S. Zhang 4,1 1 University
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 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 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 informationValidation studies on local gyrokinetic simulations of tokamak ITG-TEM driven turbulent transport
Validation studies on local gyrokinetic simulations of tokamak ITG-TEM driven turbulent transport *Motoki Nakata 1, Mitsuru Honda 1, Maiko Yoshida 1, Hajime Urano 1, Shinya Maeyama 1, Masanori Nunami 2,
More informationValidating Simulations of Multi-Scale Plasma Turbulence in ITER-Relevant, Alcator C-Mod Plasmas
Validating Simulations of Multi-Scale Plasma Turbulence in ITER-Relevant, Alcator C-Mod Plasmas Nathan Howard 1 with C. Holland 2, A.E. White 1, M. Greenwald 1, J. Candy 3, P. Rodriguez- Fernandez 1, and
More informationGyrokinetic Large Eddy Simulations
Gyrokinetic Large Eddy Simulations A. Bañón Navarro 1, P. Morel 1, M. Albrecht-Marc 1, D. Carati 1, F. Merz 2, T. Görler 2, and F. Jenko 2 1 Laboratoire de Physique Statistique et des Plasmas Université
More informationA Numerical Method for Parallel Particle Motions in Gyrokinetic Vlasov Simulations )
A Numerical Method for Parallel Particle Motions in Gyrokinetic Vlasov Simulations ) Shinya MAEYAMA, Akihiro ISHIZAWA 1), Tomo-Hiko WATANABE 1), Noriyoshi NAKAJIMA 1), Shunji TSUJI-IIO and Hiroaki TSUTSUI
More informationMechanisms of intrinsic toroidal rotation tested against ASDEX Upgrade observations
Mechanisms of intrinsic toroidal rotation tested against ASDEX Upgrade observations William A. Hornsby C. Angioni, E. Fable, P. Manas, R. McDermott, Z.X. Lu, S. Grosshauser 2, A. G. Peeters 2 and the ASDEX
More informationCoupled radius-energy turbulent transport of alpha particles
Coupled radius-energy turbulent transport of alpha particles George Wilkie, Matt Landreman, Ian Abel, William Dorland 24 July 2015 Plasma kinetics working group WPI, Vienna Wilkie (Maryland) Coupled transport
More informationDPG School The Physics of ITER Physikzentrum Bad Honnef, Energy Transport, Theory (and Experiment) Clemente Angioni
Max-Planck-Institut für Plasmaphysik DPG School The Physics of ITER Physikzentrum Bad Honnef, 23.09.2014 Energy Transport, Theory (and Experiment) Clemente Angioni Special acknowledgments for material
More informationBlock-Structured Grids in Full Velocity Space for Eulerian Gyrokinetic Simulations
Block-Structured Grids in Full Velocity Space for Eulerian Gyrokinetic Simulations D. Jarema a,, H.J. Bungartz a, T. Görler b, F. Jenko c,b, T. Neckel a, D. Told c,b a Institut für Informatik, Technische
More informationINTRODUCTION TO GYROKINETIC AND FLUID SIMULATIONS OF PLASMA TURBULENCE AND OPPORTUNITES FOR ADVANCED FUSION SIMULATIONS
INTRODUCTION TO GYROKINETIC AND FLUID SIMULATIONS OF PLASMA TURBULENCE AND OPPORTUNITES FOR ADVANCED FUSION SIMULATIONS G.W. Hammett, Princeton Plasma Physics Lab w3.pppl.gov/ hammett Fusion Simulation
More informationTRANSPORT PROGRAM C-MOD 5 YEAR REVIEW MAY, 2003 PRESENTED BY MARTIN GREENWALD MIT PLASMA SCIENCE & FUSION CENTER
TRANSPORT PROGRAM C-Mod C-MOD 5 YEAR REVIEW MAY, 2003 PRESENTED BY MARTIN GREENWALD MIT PLASMA SCIENCE & FUSION CENTER C-MOD - OPPORTUNITIES AND CHALLENGES Prediction and control are the ultimate goals
More informationBlock-Structured Grids for Eulerian Gyrokinetic Simulations
Block-Structured Grids for Eulerian Gyrokinetic Simulations D. Jarema a,, H.J. Bungartz a, T. Görler b, F. Jenko c,b, T. Neckel a, D. Told c,b a Institut für Informatik, Technische Universität München,
More informationPlasma instabilities. Dr Ben Dudson, University of York 1 / 37
Plasma instabilities Dr Ben Dudson, University of York 1 / 37 Previously... Plasma configurations and equilibrium Linear machines, and Stellarators Ideal MHD and the Grad-Shafranov equation Collisional
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 informationCharacterizing electron temperature gradient turbulence via numerical simulation
Characterizing electron temperature gradient turbulence via numerical simulation W. M. Nevins Lawrence Livermore National Laboratory, Livermore, California 94551 J. Candy General Atomics, San Diego, California
More informationDoppler Reflectometry Simulations for ASDEX Upgrade
Doppler Reflectometry Simulations for ASDEX Upgrade C. Lechte IGVP University of Stuttgart Pfaffenwaldring 31, 70569 Stuttgart Germany Phone +49 711 685 62306 Fax +49 711 685 63102 G. D. Conway, T. Görler,
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 informationAccurate representation of velocity space using truncated Hermite expansions.
Accurate representation of velocity space using truncated Hermite expansions. Joseph Parker Oxford Centre for Collaborative Applied Mathematics Mathematical Institute, University of Oxford Wolfgang Pauli
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 informationGyrokinetic Turbulence Simulations at High Plasma Beta
Gyrokinetic Turbulence Simulations at High Plasma Beta Moritz J. Pueschel Thanks to F. Jenko and M. Kammerer Ringberg Theory Meeting, Nov. 18, 2008 1 Motivation 2 3 The Beta Parameter Definition β β e
More informationZ. Lin University of California, Irvine, CA 92697, USA. Supported by SciDAC GPS-TTBP, GSEP & CPES
GTC Framework Development and Application Z. Lin University of California, Irvine, CA 92697, USA and dgpsttbp GPS-TTBP Team Supported by SciDAC GPS-TTBP, GSEP & CPES GPS-TTBP Workshop on GTC Framework
More informationNumKin, Strasbourg, October 17 th, 2016
F. Palermo 1 A.Biancalani 1, C.Angioni 1, F.Zonca 2, A.Bottino 1, B.Scott 1, G.D.Conway 1, E.Poli 1 1 Max Planck Institut für Plasmaphysik, Garching, Germany 2 ENEA C. R. Frascati - Via E. Fermi 45, CP
More informationTokamak Edge Turbulence background theory and computation
ASDEX Upgrade Tokamak Edge Turbulence background theory and computation B. Scott Max Planck Institut für Plasmaphysik Euratom Association D-85748 Garching, Germany Krakow, Sep 2006 Outline Basic Concepts
More informationGyrokinetic simulation of collisionless trapped-electron mode turbulence
PHYSICS OF PLASMAS 1, 07309 005 Gyrokinetic simulation of collisionless trapped-electron mode turbulence Tilman Dannert a and Frank Jenko Max-Planck Institut für Plasmaphysik, EURATOM Association, 85748
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 informationPerformance of the fusion code GYRO on three four generations of Crays. Mark Fahey University of Tennessee, Knoxville
Performance of the fusion code GYRO on three four generations of Crays Mark Fahey mfahey@utk.edu University of Tennessee, Knoxville Contents Introduction GYRO Overview Benchmark Problem Test Platforms
More informationSelf-consistent particle tracking in a simulation of the entropy mode in a Z pinch
Self-consistent particle tracking in a simulation of the entropy mode in a Z pinch K. Gustafson, I. Broemstrup, D. del-castillo-negrete, W. Dorland and M. Barnes Department of Physics, CSCAMM, University
More informationOn the Nature of ETG Turbulence and Cross-Scale Coupling
J. Plasma Fusion Res. SERIES, Vol. Vol. 6 6 (2004) (2004) 11 16 000 000 On the Nature of ETG Turbulence and Cross-Scale Coupling JENKO Frank Max-Planck-Institut für Plasmaphysik, EURATOM-Association, D-85748
More informationNonlinear quasisteady state benchmark of global gyrokinetic codes
PHYSICS OF PLASMAS 7, 3 Nonlinear quasisteady state benchmark of global gyrokinetic codes X. Lapillonne, B. F. McMillan, T. Görler, S. Brunner, T. Dannert, F. Jenko, F. Merz, and L. Villard Ecole Polytechnique
More informationCritical gradient formula for toroidal electron temperature gradient modes
PHYSICS OF PLASMAS VOLUME 8, NUMBER 9 SEPTEMBER 2001 Critical gradient formula for toroidal electron temperature gradient modes F. Jenko, W. Dorland, a) and G. W. Hammett b) Max-Planck-Institut für Plasmaphysik,
More informationTurbulent Transport due to Kinetic Ballooning Modes in High-Beta Toroidal Plasmas
1 TH/P-3 Turbulent Transport due to Kinetic allooning Modes in High-eta Toroidal Plasmas A. Ishizawa 1, S. Maeyama, T.-H. Watanabe 1, H. Sugama 1 and N. Nakajima 1 1 National Institute for Fusion Science,
More informationMonte-Carlo finite elements gyrokinetic simulations of Alfven modes in tokamaks.
Monte-Carlo finite elements gyrokinetic simulations of Alfven modes in tokamaks. A. Bottino 1 and the ORB5 team 1,2,3 1 Max Planck Institute for Plasma Physics, Germany, 2 Swiss Plasma Center, Switzerland,
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 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 informationGreg Hammett Imperial College, London & Princeton Plasma Physics Lab With major contributions from:
Greg Hammett Imperial College, London & Princeton Plasma Physics Lab With major contributions from: Steve Cowley (Imperial College) Bill Dorland (Imperial College) Eliot Quataert (Berkeley) LMS Durham
More informationLocal gyrokinetic turbulence simulations with realistic tokamak geometries towards ITER and JT-60SA
Local gyrokinetic turbulence simulations with realistic tokamak geometries towards ITER and JT-6SA M. Nakata, A. Matsuyama, N. Aiba, S. Maeyama in collaboration with T. -H. Watanabe, H. Sugama, and M.
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 informationKINETIC DESCRIPTION OF MAGNETIZED TECHNOLOGICAL PLASMAS
KINETIC DESCRIPTION OF MAGNETIZED TECHNOLOGICAL PLASMAS Ralf Peter Brinkmann, Dennis Krüger Fakultät für Elektrotechnik und Informationstechnik Lehrstuhl für Theoretische Elektrotechnik Magnetized low
More informationLow-collisionality density-peaking in GYRO simulations of C-Mod plasmas
Low-collisionality density-peaking in GYRO simulations of C-Mod plasmas D. R. Mikkelsen, M. Bitter, K. Hill, PPPL M. Greenwald, J.W. Hughes, J. Rice, MIT J. Candy, R. Waltz, General Atomics APS Division
More informationComputational Issues in the Continuum Gyrokinetic Code GYRO
Computational Issues in the Continuum Gyrokinetic Code GYRO presented by: Eric Bass GSEP SciDAC project at General Atomics CScADS Workshop Snowbird, UT July 19, 2010 About GYRO Purpose: To predict transport
More informationValidation Study of gyrokinetic simulation (GYRO) near the edge in Alcator C-Mod ohmic discharges
Validation Study of gyrokinetic simulation (GYRO) near the edge in Alcator C-Mod ohmic discharges C. Sung, A. E. White, N. T. Howard, D. Mikkelsen, C. Holland, J. Rice, M. Reinke, C. Gao, P. Ennever, M.
More informationTurbulence spreading and transport scaling in global gyrokinetic particle simulations
PHYSICS OF PLASMAS VOLUME 11, NUMBER 3 MARCH 2004 Turbulence spreading and transport scaling in global gyrokinetic particle simulations Z. Lin Department of Physics and Astronomy, University of California,
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 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 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 informationMulti-scale turbulence, electron transport, and Zonal Flows in DIII-D
Multi-scale turbulence, electron transport, and Zonal Flows in DIII-D L. Schmitz1 with C. Holland2, T.L. Rhodes1, G. Wang1, J.C. Hillesheim1, A.E. White3, W. A. Peebles1, J. DeBoo4, G.R. McKee5, J. DeGrassie4,
More informationFine-Scale Zonal Flow Suppression of Electron Temperature Gradient Turbulence
Fine-Scale Zonal Flow Suppression of Electron Temperature Gradient Turbulence S.E. Parker, J.J. Kohut, Y. Chen, Z. Lin, F.L. Hinton and W.W. Lee Center for Integrated Plasma Studies, University of Colorado,
More informationC-Mod Transport Program
C-Mod Transport Program PAC 2006 Presented by Martin Greenwald MIT Plasma Science & Fusion Center 1/26/2006 Introduction Programmatic Focus Transport is a broad topic so where do we focus? Where C-Mod
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 informationTomo-Hiko Watanabe and Shinya Maeyama Department of Physics, Nagoya University, Japan
Tomo-Hiko Watanabe and Shinya Maeyama Department of Physics, Nagoya University, Japan watanabe.tomohiko@nagoya-u.jp www.p.phys.nagoya-u.ac.jp 2017/1/30 DMC2017@IPAM, UCLA 1 Outline Introduction What is
More informationMacroscopic plasma description
Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion
More informationProgress of Full- f Gyrokinetic Simulation Toward Reactor Relevant Numerical Experiments )
Progress of Full- f Gyrokinetic Simulation Toward Reactor Relevant Numerical Experiments ) Yasuhiro IDOMURA, Motoki NAKATA 1) and Sébastien JOLLIET 2) Japan Atomic Energy Agency, Kashiwanoha, 5-1-5 Kashiwa,
More informationUnderstanding and Predicting Profile Structure and Parametric Scaling of Intrinsic Rotation. Abstract
Understanding and Predicting Profile Structure and Parametric Scaling of Intrinsic Rotation W. X. Wang, B. A. Grierson, S. Ethier, J. Chen, and E. Startsev Plasma Physics Laboratory, Princeton University,
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 informationLong Time Simulations of Microturbulence in Fusion Plasmas
Long Time Simulations of Microturbulence in Fusion Plasmas W. W. Lee, S. Ethier, T. G. Jenkins, W. X. Wang, J. L. V. Lewandowski, G. Rewoldt, and W. M. Tang Princeton Plasma Physics Laboratory, Princeton,
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 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 informationSpace Plasma Physics Thomas Wiegelmann, 2012
Space Plasma Physics Thomas Wiegelmann, 2012 1. Basic Plasma Physics concepts 2. Overview about solar system plasmas Plasma Models 3. Single particle motion, Test particle model 4. Statistic description
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 informationConfinement of toroidal non-neutral plasma
10th International Workshop on Non-neutral Plasmas 28 August 2012, Greifswald, Germany 1/20 Confinement of toroidal non-neutral plasma in magnetic dipole RT-1: Magnetospheric plasma experiment Visualized
More informationGA A25566 COUPLED ITG/TEM-ETG GYROKINETIC SIMULATIONS
GA A25566 COUPLED ITG/TEM-ETG GYROKINETIC SIMULATIONS by J. CANDY and R.E. WALTZ SEPTEMBER 2006 DISCLAIMER This report was prepared as an account of work sponsored by an agency of the United States Government.
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-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 informationUCIrvine. Gyrokinetic Studies of Turbulence Spreading IAEA-CN-116/TH1-4
AEA-CN-116/TH1-4 Gyrokinetic Studies of Turbulence Spreading T.S. Hahm, Z. Lin, a P.H. Diamond, b G. Rewoldt, W.X. Wang, S. Ethier, O. Gurcan, b W. Lee, and W.M. Tang Princeton University, Plasma Physics
More informationtokamak and stellarator geometry, regarding both its physical character and its interaction
THE INFLUENCE OF ZONAL EXB FLOWS ON EDGE TURBULENCE IN TOKAMAKS AND STELLARATORS B. SCOTT, F. JENKO, A. KENDL Max-Planck-Institut fíur Plasmaphysik, Garching, Germany We report on æuid, gyroæuid and gyrokinetic
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 informationTurbulence and Transport The Secrets of Magnetic Confinement
Turbulence and Transport The Secrets of Magnetic Confinement Presented by Martin Greenwald MIT Plasma Science & Fusion Center IAP January 2005 FUSION REACTIONS POWER THE STARS AND PRODUCE THE ELEMENTS
More informationin tokamak plasmas Istvan Pusztai 1 Jeff Candy 2 Punit Gohil 2
Isotope mass and charge effects in tokamak plasmas Istvan Pusztai 1 Jeff Candy 2 Punit Gohil 2 (1) Chalmers University of Technology, Applied Physics, SE-412 96, Göteborg, Sweden (2) General Atomics, P.O.
More informationDevelopment and Validation of a Predictive Model for the Pedestal Height (EPED1)
Development and Validation of a Predictive Model for the Pedestal Height (EPED1) P.B. Snyder 1 with R.J. Groebner 1, A.W. Leonard 1, T.H. Osborne 1, M. Beurskens 3, L.D. Horton 4, A.E. Hubbard 5, J.W.
More informationKinetic Turbulence in Magnetised Plasma
5 March 2013 Kinetic Turbulence in Magnetised Plasma Alexander Schekochihin w i t h S. Cowley (Culham), G. Howes, J. TenBarge (Iowa), W. Dorland, A. Kanekar (Maryland), T. Tatsuno (Tokyo), G. Hammett (Princeton),
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 informationFull-wave Simulations of Lower Hybrid Wave Propagation in the EAST Tokamak
Full-wave Simulations of Lower Hybrid Wave Propagation in the EAST Tokamak P. T. BONOLI, J. P. LEE, S. SHIRAIWA, J. C. WRIGHT, MIT-PSFC, B. DING, C. YANG, CAS-IPP, Hefei 57 th Annual Meeting of the APS
More informationGyrokinetic theory for particle transport in fusion plasmas
Gyrokinetic theory for particle transport in fusion plasmas Matteo Valerio Falessi 1,2, Fulvio Zonca 3 1 INFN - Sezione di Roma Tre, Via della Vasca Navale, 84 (00146) Roma (Roma), Italy 2 Dipartimento
More information