HPC Simulations of Magnetic Fusion Plasmas
|
|
- Janice Conley
- 5 years ago
- Views:
Transcription
1 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 Centre de Recherches en Physique des Plasmas Association Euratom Suisse, EPFL, Lausanne, Switzerland A. Bottino, T. Dannert, T. Görler, R. Hatzky, F. Jenko, A. Mishchenko, B. Scott Max-Planck Gesellschaft, IPP Garching and Greifswald, Rechenzentrum, Germany S. Jolliet, Y. Idomura, Japan Atomic Energy Agency, Taitou, Tokyo, Japan P. Angelino, G.L. Falchetto, X. Garbet, V. Grandgirard, N. Mellet, Y. Sarazin CEA/DSM/IRFM, Cadarache, France T.H. Watanabe, National Institute for Fusion Science, Toki, Japan T. Stitt, Swiss Center for Supercomputing, CSCS, Manno, Switzerland
2 Outline 1. Introduction Complexity of the physical systems we are studying What is «special» about plasma physics? 2. Direct Numerical Simulations of Turbulence Challenges and progresses 3. Some other physics problems 3D configurations, fast particles, 4. High Performance Computing Massive parallelization 5. Conclusion DEISA PRACE
3 Outline 1. Introduction Complexity of the physical systems we are studying What is «special» in plasma physics? 2. Direct Numerical Simulations of Turbulence Challenges and progresses 3. Some other physics problems 3D configurations, fast particles, 4. High Performance Computing Massive parallelization 5. Conclusion DEISA PRACE
4 ITER: the way to fusion The plasma will be there (see next p) You are here Source:ITER EU+CH, Japan, USA, China, India, South Korea, Russian Federation > 5.3G construction cost Virtually inexhaustible, environmentally benign source of energy: from Deuterium and Lithium; gives Helium. (Tritium is recycled) DEISA PRACE
5 Magnetic confinement: tokamak B r field line helical torsion (rotational transform) of B v field lines is essential to confine particles magneticsurface B and curvature drifts Larmor radius ρ L particle trajectory Trapped particle Passing particle DEISA PRACE
6 Particle motion in a tokamak Passing (top) and trapped (bottom) particle orbits: the guiding centres almost follow the equilibrium magnetic field. 3D view (left) and 2D cross-section (right) In the plasma core, particles are very weakly collisional mean free path along the magnetic field >> system size importance of the geometry of the magnetic field DEISA PRACE Y. Idomura et al., C.R. Physique 7 (2006)
7 Complexity: many nonlinearities coils transformer induction plasma RF current drive current r j, I p operational limits equilibrium configuration heat heat bootstrap macroscopic stability r B( x r ) n T heat heat Neutral Beams fusion α particles transport turbulence microscopic stability Geometry of magnetic configuration is an essential feature of fusion plasma physics DEISA PRACE
8 Timescales in the ITER plasma machine lifetime 1 shot energy confinement turbulence ion cyclotron electron cyclotron t [s] Physics spans several orders of magnitude Direct Numerical Simulation (DNS) of everything is unthinkable Need to separate timescales using approximations How to integrate all phenomena in a consistent manner? DEISA PRACE
9 Kinetic effects: wave-particle interaction f ω / k Surfers with velocity too different from the phase velocity of the wave will not ride the wave Surfers with velocity just below the phase velocity of the wave will be accelerated -> momentum and energy transfer Net energy transfer from the wave to the particles if f / v < 0 Collisionless Landau damping General: distribution function in 6D phase space To be solved with consistent electromagnetic fields v r r f ( x, v; t) r r r r E( x, t), B( x, t) DEISA PRACE
10 Outline 1. Introduction Complexity of the physical systems we are studying What is «special» in plasma physics? 2. Direct Numerical Simulations of Turbulence Challenges and progresses 3. Some other physics problems 3D configurations, fast particles, 4. High Performance Computing Massive parallelization 5. Conclusion DEISA PRACE
11 Turbulence and Transport Hot plasma core, colder edge: => grad T grad T > (grad T) crit => instabilities Small scales, i.e. spatial scales ~ρ s (Larmor radius) large scales, i.e. spatial scales ~ system size a low frequency, i.e. ω << Ω i (cyclotron frequency) Growth and nonlinear saturation of unstable modes => turbulent state => transport (heat, particles, momentum) anomalous, essentially by collisionless processes Turbulent transport >> collisional transport In spite of this, typical grad T ~ 10 6 K / cm in tokamaks Need a theory and numerical tools! DEISA PRACE
12 Gyrokinetic model Assume ω << turbulence ω ion cyclotron Average out the fast motion of the particle around the guiding center Fast parallel motion, slow perpendicular motion (drifts) Strong anisotropy of turbulent perturbations (// vs perp to B) ξ ξ phase space dimension reduction 6D ---> 5D r r r fs( q, p; t) fs( R, u, μ; t) r R = q r r ρ s k// Zonal Flows ion-driven k electron-driven k [m 1 ] DEISA PRACE
13 Gyrokinetic-Maxwell: summary Assuming frequencies << cyclotron frequencies, analytical phase space reduction from 6D to 5D A time evolution partial differential equation (PDE) describing the advection-diffusion in 5D phase space of the distribution function f A set of 5 nonlinear, coupled ordinary differential equations (ODEs) for the characteristics A set of integral partial differential equations for the perturbed fields (potentials φ, A ) in 3D real space DEISA PRACE
14 Gyrokinetic equations r f s ( R, v//, μ) r f s dr fs dv// f s + r + r = C( fs, fs' ) t dt R dt R r dr r dv r // =... fct( φ, A), =...fct( φ, A) dt dt ( φ, A r ) distribution function of species s in 5D phase space solution of Maxwell s equations, with ρ, j obtained as moments of f s advection-diffusion equations of motion (characteristics) (orbits) DEISA PRACE
15 Solving GK: numerical approaches Lagrangian: Particle In Cell (PIC), Monte Carlo Sample phase space with markers and follow their orbits Noise accumulation is difficult to control in long simulations Semi-Lagrangian Fixed grid in phase space, trace orbits back in time Multi-dimensional interpolation is difficult Eulerian Fixed grid in phase space, finite differences for operators CFL condition; overshoot versus dissipation Common to all three: Poisson (+Ampere) field solver Various methods: finite differences, finite elements, FFT, ORB5 (Lagrange-PIC) code: in-house (CRPP) + Max-Planck-IPP GENE (Euler) code: Max-Planck-IPP + CRPP GYSELA (Semi-Lagrange) code: CEA (+CRPP) DEISA PRACE
16 Global linear simulation: unstable mode Linearised δf Contours of perturbed potential CYCLONE base case 1/ρ*=180. GYGLES code [M. Fivaz et al, CPC 111 (1998) 27] DEISA PRACE
17 Nonlinear : Zonal Flows regulate turbulence Color contours of perturbed electrostatic potential hot core turbulent eddies => radial transport cold edge large scale Zonal Flows => breakup of turbulent eddies => self-regulation of turbulence Direct Numerical Simulation of Ion Temperature Gradient driven turbulence ORB5 code, resolution: 128x448x320, N = 83 million, Ω i Δt=40 ~6 h on 1024 processors IBM BG/L@EPFL No sources: final state has (undamped) Zonal Flows and decaying turbulence [S. Jolliet, et al., CPC 2007] DEISA PRACE
18 Nonlinear Direct Numerical Simulation Global gyrokinetic ORB5 code. ITG-driven turbulence. Contours of φ - <φ> (~ perturbed density) Zonal flows breakup eddies and regulate turbulence DEISA PRACE
19 Noise control is essential for Lagrange-PIC (with noise control) (without noise control) McMillan, et al., PoP 2008 When S/N is lost, there is a decrease in the computed transport level Loss of the large bursts DEISA PRACE
20 Noise control is essential to preserve the perturbation structures without noise control with noise control McMillan, et al., PoP 2008 Contours of φ-<φ> (~ perturbed electron density) DEISA PRACE
21 Turbulent processes have a chaotic nature Convergence studies with increasing N Signature of chaos: finite time after which two neighbouring initial conditions diverge from each other How to define a good (converged???) estimate of the transport? B.F. McMillan, et al, PoP 15 (2008) DEISA PRACE
22 Scaling with physical system size TCV DIII-D JET ITER GENE-local GENE ORB5 GENE-local s-α inconsistent geometry [J. Candy, et al., Phys. Plasma, 2004] [X. Lapillonne, PhD Thesis, 2010] Discrepancy between the GYRO and GTC codes: now resolved Importance of finite size effects in existing tokamaks for extrapolation to ITER DEISA PRACE
23 Semi-Lagrange approach Main motivation: get the better of two worlds avoid the statistical noise problem inherent to Lagrange PIC Monte Carlo methods avoid the CFL condition inherent to Euler methods t t k + Δt ( x i, v j ) f=const along orbit GYSELA code: 5D gyrokinetic t k Orbit, backward integrated from grid point Interpolation of foot of orbit Operator splitting Leapfrog scheme, with averaging of integer and half-integer timesteps [V. Grandgirard, et al., J. Comput. Phys. 217 (2006) 395] DEISA PRACE
24 Lagrange (ORB5) vs ½-Lagrange (GYSELA) Heat diffusivity Ion temperature gradient Heat diffusivity Alternative approaches bring useful cross-checks of numerical and physical results Fit from other codes (Dimits) Ion temperature gradient N s N χ Nϕ = , a / ρ = 184 [V. Grandgirard, et al., PPCF 49 (2007) B173] ORB5: N p = 256M, GYSELA : N v L // N μ = 32 8, half - torus DEISA PRACE
25 Outline 1. Introduction Complexity of the physical systems we are studying What is «special» in plasma physics? 2. Direct Numerical Simulations of Turbulence Challenges and progresses 3. Some other physics problems 3D configurations, fast particles, 4. High Performance Computing Massive parallelization 5. Conclusion DEISA PRACE
26 The plasma edge A. Kendl, B.D. Scott, PoP 13, (2006) The plasma crosssection is shaped: elongation, triangularity As we approach the edge, the geometry becomes more complex: separatrix X-point Occurrence of transport barrier : pedestal -like n and T profiles with steep gradients, leading to overall improved confinement. An indispensable ingredient for the success of ITER DEISA PRACE
27 Alternative to the tokamak: stellarator Inside the largest stellarator : LHD (NIFS, Japan) Unlike the tokamak, the stellarator is an intrisically steady-state magnetic confinement concept (does not need a net plasma current, does not need a transformer), but DEISA PRACE
28 Alfvén Eigenmode in W7-X Toroidicity-induced Alfvén Eigenmode Could be destabilised by fast ions Global wave code LEMAN: DNS of Maxwell s equations (finite elements, Fourier) Solves 4-potential (φ,a) in order to avoid spectral pollution f=65.3khz <B mag >=4.2T DEISA PRACE
29 Outline 1. Introduction Complexity of the physical systems we are studying What is «special» in plasma physics? 2. Direct Numerical Simulations of Turbulence Challenges and progresses 3. Some other physics problems 3D configurations, fast particles, 4. High Performance Computing Massive parallelization and scaling with system size 5. Conclusion DEISA PRACE
30 Parallelization scheme ORB5 Aim: solve very large problems Domain decomposition (toroidal direction) => huge number of processors MPI move particles Domain Cloning: field quantities (density, potential) are replicated MPI global sum fields N c = 1 ~ 100 N d = 100 ~ 1000 N pes = N d *N c DEISA PRACE
31 Massive parallelism - Scalability 2x10 9 particles grid 128x512x256=16.8x10 6 BG/P (Babel, Idris) (Deisa Extreme Computing Initiative) 4096 to cores XT5 (Monte Rosa, CSCS) (thanks to Tim Stitt) 1024 to 8192 cores 256k particles/core grid 256x1024x512=134x10 6 DEISA PRACE
32 Massive parallelism - Scalability 2x10 9 particles grid 128x512x256=16.8x10 6 BG/P (Babel, Idris) (Deisa Extreme Computing Initiative) 4096 to cores XT5 (Monte Rosa, CSCS) (thanks to Tim Stitt) 1024 to 8192 cores 256k particles/core grid 256x1024x512=134x10 6 DEISA PRACE
33 ITER global gyrokinetic turbulence simulation: computational requirements Turbulence spatial scales down to gyroradius ρ s = c s /Ω i c s =(T e /m i ) 1/2 Normalized size = system size / gyration radius : 1 / ρ* = a / ρ s ITER: 1 / ρ* = ~ D field solver: min 4 pts/wavelength 7x10 9 grid cells (3D) Velocity space resolution: ~ 300 / cell 2 000x10 9 phase space cells or markers (5D) Time until statistical steady-state: ~ 2000 a / c s ~10 6 time steps This is clearly unfeasible! Cost of algorithm ~ (1 / ρ*) 4 ~ (system size) 4 DEISA PRACE
34 Scalability Large Physical System Size ~ Half ITER size => 5.7x10 9 particles =>grid 512x1024x1024 (537x10 6 grid points) BG/P (JUGENE, Jülich) to cores (EFDA-HPC-FF) A parallel data transpose operation scales very badly (left) and starts to dominate the time (right) when using several racks DEISA PRACE
35 Scalability Large Physical System Size speedup ORB5 BG/P JUGENE strong scaling ideal particle push total pptransp ncores x 10 4 ~ Half ITER size => 5.7x10 9 particles =>grid 512x1024x1024 (537x10 6 grid points) BG/P (JUGENE, Jülich) to cores (EFDA-HPC-FF) A parallel data transpose operation scales very badly (left) and starts to dominate the time (right) when using several racks DEISA PRACE
36 Global large system size linear simulation ITER-size cyclone case 1/ρ* = 1120 m~170 wavelengths, solved with N θ =64 grid points: how is it possible? DEISA PRACE
37 Turbulence in magnetized plasmas hot core is very anisotropic [S. Jolliet et al., Comput. Phys. Commun. 177, 409 (2007)] DEISA PRACE
38 At the heart of the HP2C proposal: Field-aligned coordinates reduction of grid by a factor ~ 1/ρ* (1000 for ITER) requires major refactoring of existing codes (ORB5, GENE) and/or new code development (semi-lagrangian) DEISA PRACE
39 ITER global gyrokinetic turbulence simulation: computational requirements Turbulence spatial scales down to gyroradius ρ s = c s /Ω i c s =(T e /m i ) 1/2 Normalized size = system size / gyration radius : 1 / ρ* = a / ρ s ITER: 1 / ρ* = ~ D field solver: min 4 pts/wavelength 7x10 9 grid cells (3D) 60x10 6 Velocity space resolution: ~ 300 / cell 2 000x10 9 phase space cells or markers (5D) Time until statistical steady-state: ~ 2000 a / c s 18x10 9 ~10 6 time steps 2x10 4 This is clearly unfeasible! Feasible (*) Cost of algorithm ~ (1 / ρ*) 4 ~ (system size) 4 (1 / ρ*) 2 (*) This is one of the goals of the HP2C project. Requires major work! DEISA PRACE
40 Conclusion and outlook Fusion plasma physics is an extraordinary rich and diverse field This diversity is reflected in the various theoretical and numerical approaches Direct Numerical Simulation of Turbulence remains very challenging Future developments will include more physical effects and require further numerical advances Numerical codes for the simulation of fusion plasmas are mature for the current top-end HPC platforms Further developments are required to adapt them to the future generation of HPC mainly a human resource issue! Special thanks to Prof. Thomas Schulthess for the HP2C project Thank you for your attention DEISA PRACE
The gyrokinetic turbulence code GENE - Numerics and applications
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
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 informationGyrokinetics 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationBenchmark of gyrokinetic, kinetic MHD and gyrofluid codes for the linear calculation of fast particle driven TAE dynamics
ITR/P-34 Benchmark of gyrokinetic, kinetic MHD and gyrofluid codes for the linear calculation of fast particle driven TAE dynamics A.Könies,S.Briguglio 2,N.Gorelenkov 3,T.Fehér 4,M.Isaev 5,P.Lauber 4,A.Mishchenko,
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 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 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 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 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 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 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 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 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 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 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 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 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 informationPLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT. Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford
1 PLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford 2 Overview Why do we need plasmas? For fusion, among other things
More informationExploration of Configurational Space for Quasi-isodynamic Stellarators with Poloidally Closed Contours of the Magnetic Field Strength
Exploration of Configurational Space for Quasi-isodynamic Stellarators with Poloidally Closed Contours of the Magnetic Field Strength V.R. Bovshuk 1, W.A. Cooper 2, M.I. Mikhailov 1, J. Nührenberg 3, V.D.
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 informationExperimental Study of the Stability of Alfvén Eigenmodes on JET
IAEA FEC, Paper EX/P-7 Experimental Study of the Stability of Alfvén Eigenmodes on JET D.Testa,, A.Fasoli,, G.Fu 4, A.Jaun 3, D.Borba, P.de Vries 6, and JET-EFDA contributors [] Plasma Science and Fusion
More 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 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 informationA particle-in-cell method with adaptive phase-space remapping for kinetic plasmas
A particle-in-cell method with adaptive phase-space remapping for kinetic plasmas Bei Wang 1 Greg Miller 2 Phil Colella 3 1 Princeton Institute of Computational Science and Engineering Princeton University
More informationTransport and turbulence reduction with negative triangularity : Correlation ECE measurements in TCV
Transport and turbulence reduction with negative triangularity : Correlation ECE measurements in TCV A. Pochelon, M.Rancic, V.S.Udintsev 1, T.P.Goodman, E.Fable 2, B.Labit, O.Sauter, R.Behn, A.Bottino
More informationGyrokinetic simulations of magnetic fusion plasmas
Gyrokinetic simulations of magnetic fusion plasmas Tutorial 1 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr
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 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 informationComparison of Ion Internal Transport Barrier Formation between Hydrogen and Helium Dominated Plasmas )
Comparison of Ion Internal Transport Barrier Formation between Hydrogen and Helium Dominated Plasmas ) Kenichi NAGAOKA 1,2), Hiromi TAKAHASHI 1,2), Kenji TANAKA 1), Masaki OSAKABE 1,2), Sadayoshi MURAKAMI
More informationBeams and magnetized plasmas
Beams and magnetized plasmas 1 Jean-Pierre BOEUF LAboratoire PLAsma et Conversion d Energie LAPLACE/ CNRS, Université Paul SABATIER, TOULOUSE Beams and magnetized plasmas 2 Outline Ion acceleration and
More informationPresented on the 15th International Stellarator Workshop, Madrid, 3-7 Oct Submitted to Fusion Science and Technology Revised 20 March 2006
Presented on the 15th International Stellarator Workshop, Madrid, 3-7 Oct. 2005 Submitted to Fusion Science and Technology Revised 20 March 2006 VENUS+ffif - A BOOTSTRAP CURRENT CALCULATION MODULE FOR
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 informationTurbulence and transport reduction with innovative plasma shapes in TCV - correlation ECE measurements and gyrokinetic simulations
Turbulence and transport reduction with innovative plasma shapes in TCV - correlation ECE measurements and gyrokinetic simulations A. Pochelon, and the TCV team 1 Ecole Polytechnique de Lausanne (EPFL)
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 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 informationMichel Mehrenberger 1 & Eric Sonnendrücker 2 ECCOMAS 2016
for Vlasov type for Vlasov type 1 2 ECCOMAS 2016 In collaboration with : Bedros Afeyan, Aurore Back, Fernando Casas, Nicolas Crouseilles, Adila Dodhy, Erwan Faou, Yaman Güclü, Adnane Hamiaz, Guillaume
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 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 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 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 informationTH/P4-9. T. Takizuka 1), K. Shimizu 1), N. Hayashi 1), M. Hosokawa 2), M. Yagi 3)
1 Two-dimensional Full Particle Simulation of the Flow Patterns in the Scrape-off-layer Plasma for Upper- and Lower- Null Point Divertor Configurations in Tokamaks T. Takizuka 1), K. Shimizu 1), N. Hayashi
More informationFast ion physics in the C-2U advanced, beam-driven FRC
Fast ion physics in the C-2U advanced, beam-driven FRC Richard Magee for the TAE Team 216 US-Japan Workshop on the Compact Torus August 23, 216! High β FRC embedded in magnetic mirror is a unique fast
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 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 informationA global collisionless PIC code in magnetic coordinates
A global collisionless PIC code in magnetic coordinates S. Jolliet a,, A. Bottino b, P. Angelino a, R. Hatzky c, T. M. Tran a, B. F. Mcmillan a, O. Sauter a, K. Appert a, Y. Idomura d, L. Villard a a Ecole
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 informationToroidal confinement devices
Toroidal confinement devices Dr Ben Dudson Department of Physics, University of York, Heslington, York YO10 5DD, UK 24 th January 2014 Dr Ben Dudson Magnetic Confinement Fusion (1 of 20) Last time... Power
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 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 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 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 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 informationIssues of Perpendicular Conductivity and Electric Fields in Fusion Devices
Issues of Perpendicular Conductivity and Electric Fields in Fusion Devices Michael Tendler, Alfven Laboratory, Royal Institute of Technology, Stockholm, Sweden Plasma Turbulence Turbulence can be regarded
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 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 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 informationA fast solver for the gyrokinetic field equation with adiabatic electrons
A fast solver for the gyrokinetic field equation with adiabatic electrons M.Borchardt, R.Kleiber Max-Planck-Institut für Plasmaphysik, EURATOM Association 17491 Greifswald, Germany W.Hackbusch Max-Planck-Institut
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 informationEffect of Ion Orbit Loss on Rotation and the Radial Electric Field in the DIII-D Tokamak
Effect of Ion Orbit Loss on Rotation and the Radial Electric Field in the DIII-D Tokamak by T.M. Wilks 1 with W.M. Stacey 1 and T.E. Evans 2 1 Georgia Institute of Technology 2 General Atomics Presented
More informationWaves in plasma. Denis Gialis
Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.
More 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 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 informationPer Helander. Contributions from: R. Kleiber, A. Mishchenko, J. Nührenberg, P. Xanthopoulos. Wendelsteinstraße 1, Greifswald
Rotation and zonal flows in stellarators Per Helander Wendelsteinstraße 1, 17491 Greifswald Contributions from: R. Kleiber, A. Mishchenko, J. Nührenberg, P. Xanthopoulos What is a stellarator? In a tokamak
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 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 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 informationA simple formula for the trapped fraction in tokamaks including the effect of triangularity
A simple formula for the trapped fraction in tokamaks including the effect of triangularity O. Sauter Centre de Recherches en Physique des plasmas, Ecole Polytechnique Fédérale de Lausanne, CRPP-EPFL,
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 informationIntroduction to Nuclear Fusion. Prof. Dr. Yong-Su Na
Introduction to Nuclear Fusion Prof. Dr. Yong-Su Na What is a stellarator? M. Otthe, Stellarator: Experiments, IPP Summer School (2008) 2 Closed Magnetic System ion +++ ExB drift Electric field, E - -
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 informationElectron Bernstein Wave Heating in the TCV Tokamak
Electron Bernstein Wave Heating in the TCV Tokamak A. Mueck 1, Y. Camenen 1, S. Coda 1, L. Curchod 1, T.P. Goodman 1, H.P. Laqua 2, A. Pochelon 1, TCV Team 1 1 Ecole Polytechnique Fédérale de Lausanne
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 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 informationPhysics and Modelling of a Negative Ion Source Prototype for the ITER Neutral Beam Injection
1 ITR/P1-37 Physics and Modelling of a Negative Ion Source Prototype for the ITER Neutral Beam Injection J.P. Boeuf a, G. Fubiani a, G. Hagelaar a, N. Kohen a, L. Pitchford a, P. Sarrailh a, and A. Simonin
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 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 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 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 informationOn the Physics of the L/H Transition
ASDEX Upgrade On the Physics of the L/H Transition B. Scott Max Planck Institut für Plasmaphysik Euratom Association D-85748 Garching, Germany EFDA-TTG Workshop, Sep 2010, updated Apr 2012 Outline Physical
More informationGyrokinetic neoclassical study of the bootstrap current in the tokamak edge pedestal with fully non-linear Coulomb collisions
Physics of Plasmas Gyrokinetic neoclassical study of the bootstrap current in the tokamak edge pedestal with fully non-linear Coulomb collisions Robert Hager 1, a) 1, b) and C.S. Chang Princeton Plasma
More information