NUMERICAL MODELLING OF TURBULENCE AND TRANSPORT IN THE EDGE PLASMA OF TOKAMAKS: A MANDATORY BUT CHALLENGING TOOL FOR ITER
|
|
- Louise Holland
- 5 years ago
- Views:
Transcription
1 NUMERICAL MODELLING OF TURBULENCE AND TRANSPORT IN THE EDGE PLASMA OF TOKAMAKS: A MANDATORY BUT CHALLENGING TOOL FOR ITER Patrick Tamain on behalf of the SOLEDGE2D and TOKAM3X teams Séminaire Maison de la Simulation April 11 th 2016 CEA 10 AVRIL 2012 PAGE 1
2 Producing energy with nuclear fusion Energetically favourable to fuse light nuclei (mass < Fe) The most attractive reaction: 17.6 MeV Deuterium abundant Tritium produced from Lithium Li + n He + T Deuterium: 115ppm=33mg/L
3 What is necessary to make fusion? Condition 1: heat up hot plasma (>10keV) Condition 2: confine due to energy balance n e τ E = 12T E α L T P ext P fus P losses ~E plasma /τ E Plasma ~30keV
4 Making a star on Earth? Not exactly Every single star = (very) successful fusion reactor gravitational confinement => not reproducible on Earth sun: kg/s vs 250 kg/year for power plant On Earth, 2 possible ways: large n e (~10 31 m -3 ), small τ E (~10-11 s): inertial confinement (Mégajoule) small n e (~10 20 m-3), large τ E (~1s): magnetic confinement
5 Magnetic confinement in tokamaks Cyclotron motion in magnetic field: helicoidal trajectory, radius ρ L = mv = Larmor radius qb for D + in 30keV and B = 8T : ρ L ~2mm ρ L centre guide Tokamak configuration = toroidal field (coils) + plasma current (transfo) helicoidal field around ~toroidal surfaces = magnetic (flux) surfaces Poloidal direction θ Field lines Toroidal direction φ Magnetic flux surfaces Radial direction
6 Confinement is determined by turbulence Tokamak plasma = out of thermodynamic equilibrium transverse transport => confinement degradation Measured transport >> collisional transport transverse transport dominated by turbulence ~1m n~10 29 m 3 T~0.03eV n~10 20 m 3 T~10keV n~10 18 m 3 T~10 100eV
7 ITER in sight ITER = next generation machine under construction in Cadarache Objectives: achieve for the first time Q>1 (target Q=10) assess technologies for reactors (e.g. Li breeding blankets) Dramatic change of scale with respect to current machines extrapolation based on scaling laws very uncertain need theoretical understanding to back-up
8 The edge plasma of tokamaks open field lines = Scrape-Off Layer (SOL) transition = Last Closed Flux Surface (LCFS) Two types of configurations: LCFS / separatrix JET Tore Supra limiter divertor
9 The wall, key actor of edge plasma physics Solid material = perfect sink for the plasma sink for particles, energy boundary conditions on fluxes v i = c s = T e + T i m i Stationnary shock wave Recycling of neutral particles complex atomic and molecular reactions when re-entering the plasma important particle source for the plasma (>95%) energy and momentum sink (up to 100%) sputtering of impurities limiter
10 The wall, key actor of edge plasma physics Solid material = perfect sink for the plasma sink for particles, energy boundary conditions on fluxes v i = c s = T e + T i m i Stationnary shock wave Recycling of neutral particles complex atomic and molecular reactions when re-entering the plasma important particle source for the plasma (>95%) energy and momentum sink (up to 100%) sputtering of impurities limiter
11 The specific properties of the edge Multiple consequences to the presence of the wall: 2D equilibrium (r, θ) centre exponential decay of n e et T e/i => (relatively) cold and thin plasma with strong gradients different dynamics of turbulence Poloidal equilibrium of parallel velocity in limiter geometry δn n bord ~10 100% δn n centre < 1%
12 The edge plasma determines the life expectancy of the reactor Power flux per unit of surface determined by decay length λ (a.k.a. SOL width) λ q wall = P edge 2πR 0 λf geo
13 The edge plasma determines the life expectancy of the reactor Power flux per unit of surface determined by decay length λ (a.k.a. SOL width) q wall = 6m P edge 2πR 0 λf geo 1mm (?) 100MW ~200MW. m 2 ~10 Surface of the sun q sun ~60MW. m 2 Technological limits: q max ~10 20MW. m 2 Two axes of study: understanding of λ: equilibrium between and transport => turbulence? dissipate P edge before it reaches the wall by harnessing atomic and molecular processes
14 The edge plasma sets performances H-mode = improved confinement mode (High-confinement) bifurcation beyond some heating power threshold Large additional gradients in the edge plasma => n e and T e larger in the whole core regime expected for reactor pedestal
15 The edge plasma sets performances H-mode = improved confinement mode (High-confinement) bifurcation beyond some heating power threshold Large additional gradients in the edge plasma => n e and T e larger in the whole core regime expected for reactor Very limited understanding related to drop of turbulence + sheared rotation velocity in the edge = transport barrier? recover mode H = challenge for turbulence models
16 How to model the edge plasma? Main difficulties: multi-scale physics τ turb ~μs τ E ~s l turb ~ρ l ~mm r ~m multi-physics system: balance between plasma transport atomic and molecular processes complex geometry (plasma shape + anisotropy) Transport models Closure assumptions real time hours weeks years Scaling laws Operational codes Transport codes Full fluid codes Kinetic codes Luck : edge plasma cold => collisional (ω col T 3 2) fluid approach possible
17 How to model the edge plasma? Main difficulties: multi-scale physics τ turb ~μs τ E ~s l turb ~ρ l ~mm r ~m multi-physics system: balance between plasma transport atomic and molecular processes complex geometry (plasma shape + anisotropy) Transport models Closure assumptions real time hours weeks years Scaling laws Operational codes Transport codes Full fluid codes Kinetic codes SolEdge2D TOKAM3X Luck : edge plasma cold => collisional (ω col T 3 2) fluid approach possible
18 SOLEDGE2D-EIRENE: a coarse description of transverse transport SOLEDGE2D = edge plasma mean-field transport solver Solves only the fluctuation averaged fluid equations n t +. nu = S n, u, n = n + n u = u + u n t +. nu + nu t,r = S n, u, + S n, u, t,r + n +. nu + nu t = S
19 SOLEDGE2D-EIRENE: a coarse description of transverse transport SOLEDGE2D = edge plasma mean-field transport solver Solves only the fluctuation averaged fluid equations n t +. nu = S n, u, n = n + n u = u + u n t +. nu + nu t,r = S n, u, + S n, u, t,r + n +. nu + nu t = S gradient-diffusion closure to model turbulent transport nu t,r = D n
20 SOLEDGE2D-EIRENE: a coarse description of transverse transport SOLEDGE2D = edge plasma mean-field transport solver Solves only the fluctuation averaged fluid equations n t +. nu = S n, u, 3D, mm m μs s n = n + n u = u + u n t +. nu + nu t,r = S n, u, + S n, u, t,r + n +. nu + nu t = S 2D, cm m steady-state gradient-diffusion closure to model turbulent transport nu t,r = D n
21 but very detailed plasma-wall/neutral interaction physics 4 coupled tools for a precise treatment of the plasma-wall interaction extremely rich description of plasma solid interaction and neutral particles physics Direct coupling via STYX interface
22 Multi-fluid equations Multi-fluid mean-field model for arbitrary # of species follows each degree of ionization for each element described (ex: for carbon, 6 ionized species) 3 fluid equations per ion species: particle balance, momentum balance, energy balance + 3 equations for electrons n α t + n αu α b + n α v α = S nα t m αn α u α + m α n α u α u α b + v α = p α q α n α φ + ν α u α b + m α n α ν α u α + R α + S uα t 3 2 n αt α m αn α u α 2 + with 5 2 n αt α m αn α u α 2 u α b n αt α m αn α u α 2 v α = κ α T α b + n α χ α T α q α u b + v α = D α n α n α n T v q n u φ + R u + Q + S Parallel advection Perpendicular advection Collisional closure Gradient-diffusion closure 1 2 ν α u α 2 b m αn α ν α u α
23 SOLEDGE2D-EIRENE numerics Time discretization: 1 st order (steady-state!) semi-implicit operator splitting to separate fast (implicit) and slow terms (explicit) Spatial discretization: conservative finite differences structured mesh of quadrangles Parallelization: OpenMP (sub-domains) for SOLEDGE2D plasma solver MPI for EIRENE Monte-Carlo kinetic neutrals code typical run = 2 nodes of 12 cores for 10 to 30 days (mésocentre AMU) room for improvements although EIRENE = very old code!
24 Arbitrary 2D geometry with fluxsurface aligned grid Geometry: arbitrary 2D plasma equilibrium complex geometries made possible by domain decomposition COMPASS geometry DEMO geometry (zoom on divertor)
25 A penalization method for plasma-wall boundary conditions Problem: geometry of the wall (BC) often complex and different from the plasma magnetic geometry Solution: adapt to plasma equations the penalization method developed for neutral fluids [L. Isoardi, JCP 2010] [A. Paredes, JCP 2014] Allows treatment of very complex wall geometries with little impact on the mesh of the plasma solver
26 SOLEDGE applied to main present machines and involved in the European road-map
27 Example of application of the SOLEDGE2D-EIRENE package Confrontation between ASDEX-U discharge and SOLEDGE2D simulation
28 Example of application of the SOLEDGE2D-EIRENE package Confrontation between ASDEX-U discharge and SOLEDGE2D simulation Reproduces very well the divertor heat and particle loads captures accurately most of the main processes at play
29 The main limit of the mean-field approach No theoretical drive for transverse transport coefficients should be considered as free parameters depend on the machine, the plasma and even the position in the plasma => huge number of degrees of freedom
30 Poloidal asymmetries as a textbook example A marking example: the ballooning of transverse transport poloidal asymmetries of plasma flows observed in all machines observation of turbulence only on the outboard side of the poloidal section interpreted as enhanced turbulent transport on the outboard side b
31 Mean-field codes unable to capture poloidal asymmetries self-consistently Mean-field codes fail at reproducing the basic feature of poloidal flow asymmetries must enhance artificially transport coefficients by x100 on the outboard side, and trends not reproduced Parallel Mach number Radial direction [Erents et al., PPCF (2000)] [Erents et al., PPCF (2004)]
32 The promises of 3D global turbulence codes Turbulent transport is fundamentally not a diffusive process 2D turbulence codes data cannot be fitted with a Fick law Necessity to solve for turbulent transport selfconsistently Direct Numerical Simulations [Garcia et al., JNM (2007)] But tackling issues such as poloidal asymmetries or divertor heat loads requires to go 3D => TOKAM3X code
33 What does it take to treat electrostatic turbulence? Adjunctions to the SOLEDGE2D model to treat turbulent transport self-consistently: 1) 2 perpendicular directions (lift the toroidal symmetry) must be 3D if geometrical features are to be adressed B 2) drift velocities for perpendicular advection i/e i/e B Φ i/e 2T u i/e = ue + u B + up i u E = B 2 u B = ± B 3) a current balance equation to solve for the electrostatic potential B B B 2 t W +. W Γ N b +. Wu i i E + Wu B =. N u B e u B +. J b + D W 2 W with W =. 1 B 2 Φ + ln N
34 TOKAM3X vs SOLEDGE2D in a nutshell: geometry and physics SOLEDGE2D 2D, realistic magnetic geometry, realistic wall geometry Transport model: mean-field diffusion, no drift Multi-fluid (impurities) (new) A&M physics by coupling with EIRENE TOKAM3X 3D (full torus), realistic magnetic geometry (axisymmetric) Transport model: fluid-drift turbulence Two-fluid (single ion species) A&M physics by coupling with EIRENE (new)
35 TOKAM3X vs SOLEDGE2D in a nutshell: numerics SOLEDGE2D Explicit / implicit operator splitting Domain decomposition TOKAM3X Explicit / implicit operator splitting Multi-level domain decomposition
36 TOKAM3X vs SOLEDGE2D in a nutshell: numerics SOLEDGE2D Explicit / implicit operator splitting Domain decomposition fdsfdsfdsfds OpenMP parallelization dsqdsqdsqdqs Shock-capturing (Roe-Marquina) schemes for advection terms Coupled to EIRENE kinetic neutrals Monte-Carlo code (MPI) via unified interface STYX PASTIX sparse LU solver for implicit operators Penalization method for boundary conditions typical run: 24 CPUs, 2 weeks TOKAM3X Explicit / implicit operator splitting Multi-level domain decomposition Hybrid MPI+openMP parallelization Shock-capturing (Roe-Marquina) schemes for advection terms Coupled to EIRENE kinetic neutrals Monte-Carlo code (MPI) via unified interface STYX MUMPS or PASTIX sparse LU solver for implicit operators Standard boundary conditions typical run: 256 CPUs, 1 month
37 TOKAM3X applied to poloidal asymmetries (1) Parameters: ρ = 1/256, D = , η = , q parabolic 3-6 φ N A priori observations from simulations output: change of turbulence characteristics across separatrix filamentary (field-aligned) structures poloidal asymmetry (ballooning)
38 TOKAM3X applied to poloidal asymmetries (2) Poloidal inhomogeneities of turbulent transport recovered consistently (no free parameter) Turbulence properties outboard inboard and core SOL Density N Potential fluctuations Φ
39 Density (au) TOKAM3X applied to poloidal asymmetries (2) Poloidal inhomogeneities of turbulent transport recovered consistently (no free parameter) Turbulence properties outboard inboard and core SOL Turbulence amplitude outboard >> inboard N N t,φ Time (10 5 /ω ci ) 0.05
40 TOKAM3X applied to poloidal asymmetries (2) Poloidal inhomogeneities of turbulent transport recovered consistently (no free parameter) Turbulence properties outboard inboard and core SOL Turbulence amplitude outboard >> inboard Transverse particle flux strongly ballooned (70% in ±60 sector) Γ r E B t,φ %
41 TOKAM3X applied to poloidal asymmetries (2) Poloidal inhomogeneities of turbulent transport recovered consistently (no free parameter) Turbulence properties outboard inboard and core SOL Turbulence amplitude outboard >> inboard Transverse particle flux strongly ballooned (70% in ±60 sector) Parallel flow asymmetries quantitatively recovered 1.5 M > 1 M t,φ Where the stagnation point should be
42 TOKAM3X catching up with SOLEDGE2D When they have reached full maturity: SOLEDGE2D TOKAM3X until recently, TOKAM3X lagging behind SOLEDG2D on several aspect but progressively catching up Recent TOKAM3X advances: inclusion of neutrals physics via coupling with EIRENE Electron density (au) Microturbulence filaments Ionization source (log 10 scale, au) Wall recycling of filaments Limiter recycling
43 TOKAM3X catching up with SOLEDGE2D When they have reached full maturity: SOLEDGE2D TOKAM3X until recently, TOKAM3X lagging behind SOLEDG2D on several aspect but progressively catching up Recent TOKAM3X advances: inclusion of neutrals physics via coupling with EIRENE turbulent simulations in realistic X-point geometry φ (au)
44 TOKAM3X catching up with SOLEDGE2D on neutrals physics When they have reached full maturity: SOLEDGE2D TOKAM3X until recently, TOKAM3X lagging behind SOLEDG2D on several aspect but progressively catching up Recent TOKAM3X advances: inclusion of neutrals physics via coupling with EIRENE turbulent simulations in realistic X-point geometry non-isothermal model implemented and currently being tested but very costly (x5!) T e (au)
45 TOKAM3X s numerical knot Current balance equation contains extremely fast anisotropic operator (due to parallel dynamics) implicit scheme chosen to avoid crippling too much the time step t. 1 B 2 φ +. 1 η φ = RHS η ~ sparse linear system (~10-5 ) with very large size (1 to 250 millions DF) and very ill-conditionned (κ>10 6 ) numerical error critical for physics Current method = direct solver (MUMPS / PASTIX) main limitation to parallel scalability (>1000 CPU) and to up-scaling of cases need to explore other solutions but no luck for now (EoCoE project)
46 Bridging the gap between mean-field and turbulence approaches Are 3D global turbulence codes going to replace 2D mean-field codes? Solve the main limitation of 2D codes Still lag behind in term of description of the plasma-wall and plasma-neutral interaction but should catch up within 2-3 years BUT: computing costs 3D turbulence codes >> 2D mean-field codes For the near/mid term, meanfield codes are the only ones able to tackle ITER/DEMO size machines for routine usage and/or design purposes 2D mean-field range 3D turbulence range
47 Can the solution come from other branches of physics? CFD community facing same type of issue for aeronautic and automotive industries 2D mean-field codes RANS models Plethora of solutions developed and widely used for 50 years to avoid the cost of direct numerical simulation (DNS)
48 The k-epsilon model philosophy Example: the k-epsilon turbulence model 2 additional equations: for turbulence intensity k = 1 2 u2 and damping rate ε = u + T u 2 T t nk + nku = μ t u + u T t nε + nεu = c 1ε μ t u + u μ t k σ k μ t ε σ ε nε c 2ε n ε2 k n u u = μ t u + u T nki μ t = c μ n k2 ε New model parameters introduced and tuned from experiments Remain valid for most applications => universal parameters (vs free parameters)
49 First application to the ballooning issue Attempt to integrate similar turbulence model in SOLEDGE2D Instability and damping mechanisms adapted to the specific case of magnetized plasmas and tuned from outputs of TOKAM3X Very promising first results Spontaneous ballooning of transverse transport recovered Amplitude and trends of parallel flow asymmetries recovered
50 Summary (1) The edge plasma of tokamaks is at the centre of key issues for the success of ITER and fusion in general: Life-expectancy of the walls of the reactor Plasma performances (i.e., energetic and economical yield of the machine) Numerical simulations with predictive capabilities mandatory for design and operation purposes in sight of ITER and future IRFM and its partners have strong position in the field with the current development of 2 complementary numerical tools: SOLEDGE2D-EIRENE, a 2D mean-field code TOKAM3X(-EIRENE), a 3D global turbulence code 3D turbulence codes solve the major physical limitation of 2D mean-field codes e.g. no need for arbitrary twicking of transport coefficients to match experiments
51 Summary (2) Tremendous progress of 3D global turbulence codes means that 2D mean-field model might soon just be a subset of 3D tools However, 3D turbulent simulations are not THE solution Computational cost hinders their use for design or routine usage Unlikely to change in the next few years Currently exploring a third way inspired from the history of CFD: RANS turbulence models, e.g. k-epsilon Very promising first application to the case of poloidal asymmetries Need to be pushed further For all these approaches (2D / 3D / advanced turbulence closure), very strong need for algorithmic and numerical optimization
52 Thank you for your attention! The SOLEDGE2D and TOKAM3X teams: C. Baudoin, H. Bufferand, C. Colin, G. Ciraolo, R. Futtersack, D. Galassi, Ph. Ghendrih, G. Latu, N. Nace, F. Schwander, E. Serre, P. Tamain PAGE 52 CEA 10 AVRIL 2012 Commissariat à l énergie atomique et aux énergies alternatives Centre de Cadarache Saint Paul Lez Durance Cedex T. +33 (0) F. +33 (0) DSM IRFM SIPP Etablissement public à caractère industriel et commercial RCS Paris B
53 Back-up slides PAGE 53 CEA 10 AVRIL 2012 Commissariat à l énergie atomique et aux énergies alternatives Centre de Cadarache Saint Paul Lez Durance Cedex T. +33 (0) F. +33 (0) DSM IRFM SIPP Etablissement public à caractère industriel et commercial RCS Paris B
54 Fundamentally anisotropic physics // direction free directions constrained by B Characteristic length Characteristic time Dir. L = 2πqR~100m τ = L v th ~10 5 s Dir. L = r~1m τ E ~100ms 1s Main consequence: flux surfaces = iso-n e /T e / plasma equilibrium determined by 1D radial profiles
55 Inputs / Outputs of a mesh generation tool Minimal input: a ψ R, Z map wall line coordinates and material Output: HDF5 format For each subdomain: position of nodes of the grid in (R,Z) plane ψ for each mesh flux surf. B R, B Z, B φ at each node the number of neighbouring subdomains in each direction EIRENE mesh in EIRENE specific format
56 A new mesh generation tool for SOLEDGE A grid tool has been developed in Matlab (user-friendly GUI) to generate grids for SOLEDGE2D-EIRENE able to mesh arbitrary number of X-points and generate domain decomposition automatically (inputs: psi map & wall geometry) can be fully automated or partially hand driven Limits: generated mesh can be quite irregular (problem for TOKAM3X)
57 Careful verification process using MMS Both codes verified using Method of Manufactured Solutions (MMS) Ex: TOKAM3X
58 Computing ressources: growing needs not simple to satisfy Typical production run in 2015, SOLEDGE 24 CPUs, 1 to 2 weeks => 4 to 10 kh.cpu TOKAM3X 200 CPUs, 1 to 2 weeks => 30 to 50 kh.cpu likely to triple for SOLEDGE and quintuple for TOKAM3X in 2016 Growing allocation and consumption but 2016 will be difficult Alloc SOLEDGE Conso SOLEDGE Alloc TOKAM3X Conso TOKAM3X Key contribution of mesocentre AMU Difficulty of juggling with too many different machines (porting,reporting )
59 A strongly collaborative work Both codes developed in strong collaboration with Aix-Marseille University A fast-growing team: IRFM: H. Bufferand (SOLEDGE RO), C. Baudoin*, G. Ciraolo, J. Denis*, R. Futtersack*, Ph. Ghendrih, R. Mao*, N. Nace*, P. Tamain (TOKAM3X RO), M. Valentinnuzzi* M2P2: D. Galassi*, G. Giorgiani*, F. Schwander, E. Serre, J.-A. Soler* PIIM: Y. Marandet (STYX-EIRENE RO) + 7 PhDs & 1 post-doc already finished
60 3D edge turbulence codes: an unprecedented effort BOUT++ Drift-fluid / MHD SOL+core Flexible geometry Gradient or fluxdriven ATTEMPT Drift-fluid EM, isothermal Limiter SOL+core Gradient-driven GEMR Drift-gyrofluid EM, 8 fields Limiter SOL+core TOKAM3X Drift-fluid ES, SOL+core Flexible geometry Flux-driven GBS Drift-fluid EM, T i =0 Limiter SOL Flux-driven PAGE 60
61 Kinetic codes on the global scale: the BIT1 code Kinetic codes required to investigate assumptions made in MHD & fluid models for kinetic effects: Sheath boundary conditions Parallel transport and fluid closures Impact of non-maxwellian distribution functions on atomic physics Most advanced code: BIT1 [Tskhakaya, Contrib. Plasma Phys. 2012] 1D plasma, 2D neutrals and quasi-2d impurity models (3V) full size SOL with finest resolution (up to cells in poloidal direction) Massively parallel runs up to 8096 processors Rich atomic physics and PSI database PAGE 61
62 Example of application: heat flux limiters must be revisited for consistent parallel transport Heat-flux limiters required in fluid codes to avoid non-realistic divergence of Spitzer-Harm heat fluxes at low collisionality: Characterized by free parameter α q 1 q SH 1 nv T T 1 q SH T s 0.1?? Application of BIT1 to heat flux limiters: Limiters are strongly nonuniform along the field lines and sensitive to plasma conditions in the SOL Need for a more precise model 10-1 α e α i Perspectives of development: Enrich atomic data base 2D3V BIT2 code under test 10-2 [Tskhakaya, JNM 2009] PAGE 62
63 Edge gyrokinetic codes already available but cost will remain prohibitive in the near future 2 main motivations: GK = only approach able to combine selfconsistently the fore-mentioned physics Suspected interdependence of core and edge physics [Holland C., PoP 2009] [XGC1] On-going effort towards extending GK codes to the edge and SOL: Main existing projects: TEMPEST, XGC But edge GK comes at a huge cost: Non-adiabatic electrons necessary for sheath and interchange Extremely greedy in numerical resources (N proc ~ for ~months) Prohibitive for parameter scans [Chang, PoP 16 (2009)] PAGE 63
64 Summary: where are we standing w/r to the heat exhaust issue? The multiple facets of heat exhaust issues commend the use of a broad range of numerical tools Scaling laws Fluid diffusive transport MHD routine week(s) / 1-2 years Fluid turbulence kinetic gyrokinetic 0D 2D 3D 3D 2D 4D 5D 1 st ppleness PAGE 64
65 Summary: where are we standing w/r to the heat exhaust issue? The multiple facets of heat exhaust issues commend the use of a broad range of numerical tools 3D fluid turbulence codes: cross-field transport, L-H transition Limiter configuration: readily available (GEMR) or within reach Divertor configuration: within reach for 2 codes (BOUT++, TOKAM3X) Scaling laws Fluid diffusive transport MHD routine week(s) / 1-2 years Fluid turbulence week(s) / 1-5 years kinetic gyrokinetic 0D 2D 3D 3D 2D 4D 5D 1 st ppleness PAGE 65
66 Summary: where are we standing w/r to the heat exhaust issue? The multiple facets of heat exhaust issues commend the use of a broad range of numerical tools 3D fluid turbulence codes: cross-field transport, L-H transition Limiter configuration: readily available (GEMR) or within reach Divertor configuration: within reach for 2 codes (BOUT++, TOKAM3X) 1D / 2D (3V) PIC codes: parallel transport and kinetic effects Full SOL simulations available => feedback for fluid modelling 2D under test but costly Scaling laws Fluid diffusive transport MHD routine week(s) / 1-2 years Fluid turbulence week(s) / 1-5 years kinetic week(s) / 1-5 years gyrokinetic 0D 2D 3D 3D 2D 4D 5D 1 st ppleness PAGE 66
67 Summary: where are we standing w/r to the heat exhaust issue? The multiple facets of heat exhaust issues commend the use of a broad range of numerical tools 3D fluid turbulence codes: cross-field transport, L-H transition Limiter configuration: readily available (GEMR) or within reach Divertor configuration: within reach for 2 codes (BOUT++, TOKAM3X) 1D / 2D (3V) PIC codes: parallel transport and kinetic effects Full SOL simulations available => feedback for fluid modelling 2D under test but costly 3D (2V) gyrokinetic codes: fully consistent transport On-going but extremely challenging for edge issues Scaling laws Fluid diffusive transport MHD routine week(s) / 1-2 years Fluid turbulence week(s) / 1-5 years kinetic week(s) / 1-5 years gyrokinetic 0D 2D 3D 3D 2D 4D 5D month(s) / >5 years 1 st ppleness PAGE 67
68 A multi-facets issue The unavoidable EIRENE 3D fluid turbulence codes 3D (2V) gyrokinetic codes Pedestal Separatrix 2D + EIRENE 1D/2D kinetic 3D fluid 5D gyrokinetic Perpendicular transport 1D (3V) kinetic codes Parallel transport Atomic physics ~ ~ PAGE 68
69 A multi-facets issue The unavoidable EIRENE 3D fluid turbulence codes 3D (2V) gyrokinetic codes Pedestal Separatrix Empirical codes + EIRENE 2D + EIRENE 1D/2D kinetic 3D fluid 5D gyrokinetic Perpendicular transport 1D (3V) kinetic codes Parallel transport Atomic physics ~ ~ PAGE 69
70 A multi-facets issue The unavoidable EIRENE 3D fluid turbulence codes 3D (2V) gyrokinetic codes Pedestal Separatrix Empirical codes + EIRENE 2D + EIRENE 1D/2D kinetic 3D fluid 5D gyrokinetic Perpendicular transport 1D (3V) kinetic codes Parallel transport Atomic physics ~ ~ PAGE 70
71 The limits of the diffusive description Local Fick s law does not fit with 2D turbulence results TOKAM2D SOL interchange turbulence code Perpendicular transport better represented as convective process but what is the thermodynamic force? => non local operator?? need for new model or self-consistent simulations PAGE 71
72 Pre-study: How can we improve transport description in empirical codes? What influence has a sampled diffusion coefficient on the large-scale system? Solving
73 La turbulence de bord : un phénomène 2D? Expérience : structures turbulentes alignées sur B ~ 2D Hypothèse : transport plus rapide que => découplage / & k = 0 équations moyennées le long de B instabilité de courbure B Conservation matière + charge t D N + φ, N = σne Λ φ + S t ν φ + φ, φ = σ 1 e Λ φ g y lnn advection puits (paroi) courbure ( B) T 1 ~ Rayleigh Bénard (g = courbure de B) T 2 > T 1
74 y θ Simulations numériques avec le code TOKAM2D Simulation d une boîte 2D dans les lignes de champ ouvertes au-delà d un seuil en gradient de pression, transport turbulent 250 B x r
75 Les modèles 2D permettent de retrouver l essentiel des propriétés expérimentales Comparaison qualitative avec l expérience : forte intermittence (45% flux par évènements > 1rms) Flux (a.u.) PDF time (ms) flux (u.a.)
76 Les modèles 2D permettent de retrouver l essentiel des propriétés expérimentales Comparaison qualitative avec l expérience : forte intermittence (45% flux par évènements > 1rms) profils exponentiels => détermination largeur de décroissance λ [Sarazin, JNM (2003)]
77 Propagation balistique de structures cohérentes B a r s Phénomène d avalanches : structures perturbent l équilibre => conditions de maintien et de propagation fondamentalement multiéchelles (r-a) r s n (a.u.) t 1 t 2 = t / W i <n> ~ 40.2 r s Dr~20r s Transport = balistique diffusif M ~ (r - a) / r s n / n ~ 4.6 r s
Impact of diverted geometry on turbulence and transport barrier formation in 3D global simulations of tokamak edge plasma
1 Impact of diverted geometry on turbulence and transport barrier formation in 3D global simulations of tokamak edge plasma D. Galassi, P. Tamain, H. Bufferand, C. Baudoin, G. Ciraolo, N. Fedorczak, Ph.
More informationDrive of parallel flows by turbulence and large-scale E B transverse transport in divertor geometry
Drive of parallel flows by turbulence and large-scale E B transverse transport in divertor geometry D Galassi, Patrick Tamain, Hugo Bufferand, G Ciraolo, Philippe Ghendrih, C Baudoin, C Colin, Nicolas
More informationParallel Kelvin-Helmholtz instability in edge plasma
Journal of Physics: Conference Series OPEN ACCESS Parallel Kelvin-Helmholtz instability in edge plasma To cite this article: H Guillard et al 2014 J. Phys.: Conf. Ser. 561 012009 View the article online
More informationMean-field and turbulent transport in divertor geometry Davide Galassi
Mean-field and turbulent transport in divertor geometry Davide Galassi In collaboration with: Ph. Ghendrih, P. Tamain, C. Baudoin, H. Bufferand, G. Ciraolo, C. Colin and E. Serre Our goal: quantify turbulence
More informationImpact on Divertor Operation of the Pattern of Edge and SOL Flows Induced by Particle Sources and Sinks
1 TH/P4-26 Impact on Divertor Operation of the Pattern of Edge and SOL Flows Induced by Particle Sources and Sinks Ph. Ghendrih 1, T. Auphan 2, B. Bensiali 3, M. Bilanceri 4, K. Bodi 3,5, J. Bucalossi
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 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 informationMaths pavingthe pathto Fusion Energy & ITER. Marseille, Maths master class, May 28 th 2018, Philippe GHENDRIH
Maths pavingthe pathto Fusion Energy & ITER Marseille, Maths master class, May 28 th 2018, Philippe GHENDRIH Maths & Physics A story of Fusion Energy What is Fusion Energy Research made of? MATHS & PHYSICS
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 informationMAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT
MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT ABSTRACT A. G. Tarditi and J. V. Shebalin Advanced Space Propulsion Laboratory NASA Johnson Space Center Houston, TX
More informationEdge Momentum Transport by Neutrals
1 TH/P3-18 Edge Momentum Transport by Neutrals J.T. Omotani 1, S.L. Newton 1,2, I. Pusztai 1 and T. Fülöp 1 1 Department of Physics, Chalmers University of Technology, 41296 Gothenburg, Sweden 2 CCFE,
More informationCurrent density modelling in JET and JT-60U identity plasma experiments. Paula Sirén
Current density modelling in JET and JT-60U identity plasma experiments Paula Sirén 1/12 1/16 Euratom-TEKES Euratom-Tekes Annual Seminar 2013 28 24 May 2013 Paula Sirén Current density modelling in JET
More informationThe role of statistical noise in edge plasma transport codes based on kinetic Monte Carlo solvers for neutrals: an analogy with turbulent fluctuations
EUROFUSION WP15ER-CP(16) 15192 Y Marandet et al. The role of statistical noise in edge plasma transport codes based on kinetic Monte Carlo solvers for neutrals: an analogy with turbulent fluctuations Preprint
More informationOverview of edge modeling efforts for advanced divertor configurations in NSTX-U with magnetic perturbation fields
Overview of edge modeling efforts for advanced divertor configurations in NSTX-U with magnetic perturbation fields H. Frerichs, O. Schmitz, I. Waters, G. P. Canal, T. E. Evans, Y. Feng and V. Soukhanovskii
More informationImplementation of drift velocities and currents in SOLEDGE2D-EIRENE
Implementation of drift velocities and currents in SOLEDGE2D-EIRENE Hugo Bufferand, C Baudoin, J Bucalossi, G Ciraolo, J Denis, Nicolas Fedorczak, D Galassi, Philippe Ghendrih, R Leybros, Y Marandet, et
More informationModeling of ELM Dynamics for ITER
Modeling of ELM Dynamics for ITER A.Y. PANKIN 1, G. BATEMAN 1, D.P. BRENNAN 2, A.H. KRITZ 1, S. KRUGER 3, P.B. SNYDER 4 and the NIMROD team 1 Lehigh University, 16 Memorial Drive East, Bethlehem, PA 18015
More informationTriggering Mechanisms for Transport Barriers
Triggering Mechanisms for Transport Barriers O. Dumbrajs, J. Heikkinen 1, S. Karttunen 1, T. Kiviniemi, T. Kurki-Suonio, M. Mantsinen, K. Rantamäki 1, S. Saarelma, R. Salomaa, S. Sipilä, T. Tala 1 Euratom-TEKES
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 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 informationDrift-Driven and Transport-Driven Plasma Flow Components in the Alcator C-Mod Boundary Layer
Drift-Driven and Transport-Driven Plasma Flow Components in the Alcator C-Mod Boundary Layer N. Smick, B. LaBombard MIT Plasma Science and Fusion Center PSI-19 San Diego, CA May 25, 2010 Boundary flows
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 informationUnderstanding Edge Harmonic Oscillation Physics Using NIMROD
Understanding Edge Harmonic Oscillation Physics Using NIMROD J. King With contributions from S. Kruger & A. Pankin (Tech-X); K. Burrell, A. Garofalo, R. Groebner & P. Snyder (General Atomics) Work supported
More 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 informationRotation and Neoclassical Ripple Transport in ITER
Rotation and Neoclassical Ripple Transport in ITER Elizabeth J. Paul 1 Matt Landreman 1 Francesca Poli 2 Don Spong 3 Håkan Smith 4 William Dorland 1 1 University of Maryland 2 Princeton Plasma Physics
More informationSimple examples of MHD equilibria
Department of Physics Seminar. grade: Nuclear engineering Simple examples of MHD equilibria Author: Ingrid Vavtar Mentor: prof. ddr. Tomaž Gyergyek Ljubljana, 017 Summary: In this seminar paper I will
More 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 informationThe Levitated Dipole Experiment: Towards Fusion Without Tritium
The Levitated Dipole Experiment: Towards Fusion Without Tritium Jay Kesner MIT M.S. Davis, J.E. Ellsworth, D.T. Garnier, M.E. Mauel, P.C. Michael, P.P. Woskov MCP I3.110 Presented at the EPS Meeting, Dublin,
More informationOperational Phase Space of the Edge Plasma in Alcator C-Mod
Operational Phase Space of the Edge Plasma in B. LaBombard, T. Biewer, M. Greenwald, J.W. Hughes B. Lipschultz, N. Smick, J.L. Terry, Team Contributed talk RO.00008 Presented at the 47th Annual Meeting
More informationThe Path to Fusion Energy creating a star on earth. S. Prager Princeton Plasma Physics Laboratory
The Path to Fusion Energy creating a star on earth S. Prager Princeton Plasma Physics Laboratory The need for fusion energy is strong and enduring Carbon production (Gton) And the need is time urgent Goal
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 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 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 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 informationMultiscale modelling of sheath physics in edge transport codes
EUROFUSION WPPFC-CP(6) 6599 N Mellet et al. Multiscale modelling of sheath physics in edge transport codes Preprint of Paper to be submitted for publication in Proceedings of 26th IAEA Fusion Energy Conference
More informationImpact of neutral atoms on plasma turbulence in the tokamak edge region
Impact of neutral atoms on plasma turbulence in the tokamak edge region C. Wersal P. Ricci, F.D. Halpern, R. Jorge, J. Morales, P. Paruta, F. Riva Theory of Fusion Plasmas Joint Varenna-Lausanne International
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 informationINTRODUCTION TO MAGNETIC NUCLEAR FUSION
INTRODUCTION TO MAGNETIC NUCLEAR FUSION S.E. Sharapov Euratom/CCFE Fusion Association, Culham Science Centre, Abingdon, Oxfordshire OX14 3DB, UK With acknowledgments to B.Alper for use of his transparencies
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 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 informationProgress in characterization of the H-mode pedestal
Journal of Physics: Conference Series Progress in characterization of the H-mode pedestal To cite this article: A W Leonard 2008 J. Phys.: Conf. Ser. 123 012001 View the article online for updates and
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 informationNon-linear MHD Simulations of Edge Localized Modes in ASDEX Upgrade. Matthias Hölzl, Isabel Krebs, Karl Lackner, Sibylle Günter
Non-linear MHD Simulations of Edge Localized Modes in ASDEX Upgrade Matthias Hölzl, Isabel Krebs, Karl Lackner, Sibylle Günter Matthias Hölzl Nonlinear ELM Simulations DPG Spring Meeting, Jena, 02/2013
More informationDynamics of Zonal Shear Collapse in Hydrodynamic Electron Limit. Transport Physics of the Density Limit
Dynamics of Zonal Shear Collapse in Hydrodynamic Electron Limit Transport Physics of the Density Limit R. Hajjar, P. H. Diamond, M. Malkov This research was supported by the U.S. Department of Energy,
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 informationArbiTER studies of filamentary structures in the SOL of spherical tokamaks
ArbiTER studies of filamentary structures in the SOL of spherical tokamaks D. A. Baver, J. R. Myra, Research Corporation F. Scotti, Lawrence Livermore National Laboratory S. J. Zweben, Princeton Plasma
More informationPlasma shielding during ITER disruptions
Plasma shielding during ITER disruptions Sergey Pestchanyi and Richard Pitts 1 Integrated tokamak code TOKES is a workshop with various tools objects Radiation bremsstrahlung recombination s line s cyclotron
More informationTokamak Fusion Basics and the MHD Equations
MHD Simulations for Fusion Applications Lecture 1 Tokamak Fusion Basics and the MHD Equations Stephen C. Jardin Princeton Plasma Physics Laboratory CEMRACS 1 Marseille, France July 19, 21 1 Fusion Powers
More informationTurbulence bursts probing of transport barriers analyzed in terms of competing stochastic processes
Author manuscript, published in "Plasma Phys. Control. Fusion 55, 9 (013) 09501" DOI : 10.1088/071-3335/55/9/09501 Turbulence bursts probing of transport barriers analyzed in terms of competing stochastic
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 informationProgressing Performance Tokamak Core Physics. Marco Wischmeier Max-Planck-Institut für Plasmaphysik Garching marco.wischmeier at ipp.mpg.
Progressing Performance Tokamak Core Physics Marco Wischmeier Max-Planck-Institut für Plasmaphysik 85748 Garching marco.wischmeier at ipp.mpg.de Joint ICTP-IAEA College on Advanced Plasma Physics, Triest,
More informationIntegrated Simulation of ELM Energy Loss Determined by Pedestal MHD and SOL Transport
1 Integrated Simulation of ELM Energy Loss Determined by Pedestal MHD and SOL Transport N. Hayashi, T. Takizuka, T. Ozeki, N. Aiba, N. Oyama Japan Atomic Energy Agency, Naka, Ibaraki-ken, 311-0193 Japan
More informationCharacteristics of the H-mode H and Extrapolation to ITER
Characteristics of the H-mode H Pedestal and Extrapolation to ITER The H-mode Pedestal Study Group of the International Tokamak Physics Activity presented by T.Osborne 19th IAEA Fusion Energy Conference
More informationDer Stellarator Ein alternatives Einschlusskonzept für ein Fusionskraftwerk
Max-Planck-Institut für Plasmaphysik Der Stellarator Ein alternatives Einschlusskonzept für ein Fusionskraftwerk Robert Wolf robert.wolf@ipp.mpg.de www.ipp.mpg.de Contents Magnetic confinement The stellarator
More informationOn the physics of shear flows in 3D geometry
On the physics of shear flows in 3D geometry C. Hidalgo and M.A. Pedrosa Laboratorio Nacional de Fusión, EURATOM-CIEMAT, Madrid, Spain Recent experiments have shown the importance of multi-scale (long-range)
More informationSimulation of alpha particle current drive and heating in spherical tokamaks
Simulation of alpha particle current drive and heating in spherical tokamaks R. Farengo 1, M. Zarco 1, H. E. Ferrari 1, 1 Centro Atómico Bariloche and Instituto Balseiro, Argentina. Consejo Nacional de
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 informationA kinetic neutral atom model for tokamak scrape-off layer tubulence simulations. Christoph Wersal, Paolo Ricci, Federico Halpern, Fabio Riva
A kinetic neutral atom model for tokamak scrape-off layer tubulence simulations Christoph Wersal, Paolo Ricci, Federico Halpern, Fabio Riva CRPP - EPFL SPS Annual Meeting 2014 02.07.2014 CRPP The tokamak
More informationMagnetically Confined Fusion: Transport in the core and in the Scrape- off Layer Bogdan Hnat
Magnetically Confined Fusion: Transport in the core and in the Scrape- off Layer ogdan Hnat Joe Dewhurst, David Higgins, Steve Gallagher, James Robinson and Paula Copil Fusion Reaction H + 3 H 4 He + n
More informationA neoclassical model for toroidal rotation and the radial electric field in the edge pedestal. W. M. Stacey
A neoclassical model for toroidal rotation and the radial electric field in the edge pedestal W. M. Stacey Fusion Research Center Georgia Institute of Technology Atlanta, GA 30332, USA October, 2003 ABSTRACT
More informationModelling of JT-60U Detached Divertor Plasma using SONIC code
J. Plasma Fusion Res. SERIES, Vol. 9 (2010) Modelling of JT-60U Detached Divertor Plasma using SONIC code Kazuo HOSHINO, Katsuhiro SHIMIZU, Tomonori TAKIZUKA, Nobuyuki ASAKURA and Tomohide NAKANO Japan
More informationThermodynamical and microscopic properties of turbulent transport in the edge plasma
Home Search Collections Journals About Contact us My IOPscience Thermodynamical and microscopic properties of turbulent transport in the edge plasma This content has been downloaded from IOPscience. Please
More informationNon-linear MHD Modelling of Rotating Plasma Response to Resonant Magnetic Perturbations.
Non-linear MHD Modelling of Rotating Plasma Response to Resonant Magnetic Perturbations. M. Becoulet 1, F. Orain 1, G.T.A. Huijsmans 2, P. Maget 1, N. Mellet 1, G. Dif-Pradalier 1, G. Latu 1, C. Passeron
More informationTH/P6-14 Integrated particle simulation of neoclassical and turbulence physics in the tokamak pedestal/edge region using XGC a)
1 TH/P6-14 Integrated particle simulation of neoclassical and turbulence physics in the tokamak pedestal/edge region using XGC a) 1 Chang, C.S., 1 Ku, S., 2 Adams M., 3 D Azevedo, G., 4 Chen, Y., 5 Cummings,
More informationAnalysis and modelling of MHD instabilities in DIII-D plasmas for the ITER mission
Analysis and modelling of MHD instabilities in DIII-D plasmas for the ITER mission by F. Turco 1 with J.M. Hanson 1, A.D. Turnbull 2, G.A. Navratil 1, C. Paz-Soldan 2, F. Carpanese 3, C.C. Petty 2, T.C.
More 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 informationConnections between Particle Transport and Turbulence Structures in the Edge and SOL of Alcator C-Mod
Connections between Particle Transport and Turbulence Structures in the Edge and SOL of Alcator C-Mod I. Cziegler J.L. Terry, B. LaBombard, J.W. Hughes MIT - Plasma Science and Fusion Center th 19 Plasma
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 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 informationDOPPLER RESONANCE EFFECT ON ROTATIONAL DRIVE BY ION CYCLOTRON MINORITY HEATING
DOPPLER RESONANCE EFFECT ON ROTATIONAL DRIVE BY ION CYCLOTRON MINORITY HEATING V.S. Chan, S.C. Chiu, Y.A. Omelchenko General Atomics, San Diego, CA, U.S.A. 43rd Annual APS Division of Plasma Physics Meeting
More informationCorrelations of ELM frequency with pedestal plasma characteristics
cpp header will be provided by the publisher Correlations of ELM frequency with pedestal plasma characteristics G. Kamberov 1 and L. Popova 2 1 Stevens Institute of Technology, Hoboken NJ, USA 2 Institute
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 informationA THEORETICAL AND EXPERIMENTAL INVESTIGATION INTO ENERGY TRANSPORT IN HIGH TEMPERATURE TOKAMAK PLASMAS
A THEORETICAL AND EXPERIMENTAL INVESTIGATION INTO ENERGY TRANSPORT IN HIGH TEMPERATURE TOKAMAK PLASMAS Presented by D.P. SCHISSEL Presented to APS Centennial Meeting March 20 26, 1999 Atlanta, Georgia
More informationAlcator C-Mod. Particle Transport in the Scrape-off Layer and Relationship to Discharge Density Limit in Alcator C-Mod
Alcator C-Mod Particle Transport in the Scrape-off Layer and Relationship to Discharge Density Limit in Alcator C-Mod B. LaBombard, R.L. Boivin, M. Greenwald, J. Hughes, B. Lipschultz, D. Mossessian, C.S.
More informationElectron temperature barriers in the RFX-mod experiment
Electron temperature barriers in the RFX-mod experiment A. Scaggion Consorzio RFX, Padova, Italy Tuesday 5 th October 2010 ADVANCED PHYSICS LESSONS 27/09/2010 07/10/2010 IPP GARCHING JOINT EUROPEAN RESEARCH
More informationC-Mod Core Transport Program. Presented by Martin Greenwald C-Mod PAC Feb. 6-8, 2008 MIT Plasma Science & Fusion Center
C-Mod Core Transport Program Presented by Martin Greenwald C-Mod PAC Feb. 6-8, 2008 MIT Plasma Science & Fusion Center Practical Motivations for Transport Research Overall plasma behavior must be robustly
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 informationHybrid Simulations: Numerical Details and Current Applications
Hybrid Simulations: Numerical Details and Current Applications Dietmar Krauss-Varban and numerous collaborators Space Sciences Laboratory, UC Berkeley, USA Boulder, 07/25/2008 Content 1. Heliospheric/Space
More informationFusion Nuclear Science Facility (FNSF) Divertor Plans and Research Options
Fusion Nuclear Science Facility (FNSF) Divertor Plans and Research Options A.M. Garofalo, T. Petrie, J. Smith, V. Chan, R. Stambaugh (General Atomics) J. Canik, A. Sontag, M. Cole (Oak Ridge National Laboratory)
More informationTokamak Divertor System Concept and the Design for ITER. Chris Stoafer April 14, 2011
Tokamak Divertor System Concept and the Design for ITER Chris Stoafer April 14, 2011 Presentation Overview Divertor concept and purpose Divertor physics General design considerations Overview of ITER divertor
More informationD.J. Schlossberg, D.J. Battaglia, M.W. Bongard, R.J. Fonck, A.J. Redd. University of Wisconsin - Madison 1500 Engineering Drive Madison, WI 53706
D.J. Schlossberg, D.J. Battaglia, M.W. Bongard, R.J. Fonck, A.J. Redd University of Wisconsin - Madison 1500 Engineering Drive Madison, WI 53706 Concept Overview Implementation on PEGASUS Results Current
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 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 informationNonlocal heat flux application for Scrape-off Layer plasma
Nonlocal heat flux application for Scrape-off Layer plasma Hugo Bufferand, G. Ciraolo, P Di Cintio, N Fedorczak, Ph Ghendrih, S. Lepri, R. Livi, Yannick Marandet, E. Serre, P. Tamain To cite this version:
More informationOVERVIEW OF THE ALCATOR C-MOD PROGRAM. IAEA-FEC November, 2004 Alcator Team Presented by Martin Greenwald MIT Plasma Science & Fusion Center
OVERVIEW OF THE ALCATOR C-MOD PROGRAM IAEA-FEC November, 2004 Alcator Team Presented by Martin Greenwald MIT Plasma Science & Fusion Center OUTLINE C-Mod is compact, high field, high density, high power
More informationShear Flow Generation in Stellarators - Configurational Variations
Shear Flow Generation in Stellarators - Configurational Variations D. A. Spong 1), A. S. Ware 2), S. P. Hirshman 1), J. H. Harris 1), L. A. Berry 1) 1) Oak Ridge National Laboratory, Oak Ridge, Tennessee
More informationDriving Mechanism of SOL Plasma Flow and Effects on the Divertor Performance in JT-60U
EX/D-3 Driving Mechanism of SOL Plasma Flow and Effects on the Divertor Performance in JT-6U N. Asakura ), H. Takenaga ), S. Sakurai ), G.D. Porter ), T.D. Rognlien ), M.E. Rensink ), O. Naito ), K. Shimizu
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 informationNon-Solenoidal Plasma Startup in
Non-Solenoidal Plasma Startup in the A.C. Sontag for the Pegasus Research Team A.C. Sontag, 5th APS-DPP, Nov. 2, 28 1 Point-Source DC Helicity Injection Provides Viable Non-Solenoidal Startup Technique
More informationELMs and Constraints on the H-Mode Pedestal:
ELMs and Constraints on the H-Mode Pedestal: A Model Based on Peeling-Ballooning Modes P.B. Snyder, 1 H.R. Wilson, 2 J.R. Ferron, 1 L.L. Lao, 1 A.W. Leonard, 1 D. Mossessian, 3 M. Murakami, 4 T.H. Osborne,
More informationImpact of Toroidal Flow on ITB H-Mode Plasma Performance in Fusion Tokamak
Impact of oroidal Flow on I H-Mode Plasma Performance in Fusion okamak oonyarit Chatthong 1,*, hawatchai Onjun 1, Roppon Picha and Nopporn Poolyarat 3 1 School of Manufacturing Systems and Mechanical Engineering,
More informationGyrokinetic Theory and Dynamics of the Tokamak Edge
ASDEX Upgrade Gyrokinetic Theory and Dynamics of the Tokamak Edge B. Scott Max Planck Institut für Plasmaphysik D-85748 Garching, Germany PET-15, Sep 2015 these slides: basic processes in the dynamics
More informationJoint ITER-IAEA-ICTP Advanced Workshop on Fusion and Plasma Physics October Introduction to Fusion Leading to ITER
2267-1 Joint ITER-IAEA-ICTP Advanced Workshop on Fusion and Plasma Physics 3-14 October 2011 Introduction to Fusion Leading to ITER SNIPES Joseph Allan Directorate for Plasma Operation Plasma Operations
More informationOverview of Tokamak Rotation and Momentum Transport Phenomenology and Motivations
Overview of Tokamak Rotation and Momentum Transport Phenomenology and Motivations Lecture by: P.H. Diamond Notes by: C.J. Lee March 19, 2014 Abstract Toroidal rotation is a key part of the design of ITER
More informationfor the French fusion programme
The ITER era : the 10 year roadmap for the French fusion programme E. Tsitrone 1 on behalf of IRFM and Tore Supra team 1 : CEA, IRFM, F-13108 Saint-Paul-lez-Durance, France Association EURATOM-CEA TORE
More informationMHD. Jeff Freidberg MIT
MHD Jeff Freidberg MIT 1 What is MHD MHD stands for magnetohydrodynamics MHD is a simple, self-consistent fluid description of a fusion plasma Its main application involves the macroscopic equilibrium
More informationTheory Work in Support of C-Mod
Theory Work in Support of C-Mod 2/23/04 C-Mod PAC Presentation Peter J. Catto for the PSFC theory group MC & LH studies ITB investigations Neutrals & rotation BOUT improvements TORIC ICRF Mode Conversion
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 informationTurbulent Transport Analysis of JET H-mode and Hybrid Plasmas using QuaLiKiz, TGLF and GLF23
EFDA JET CP(1)/ B. Baiocchi, J. Garcia, M. Beurkens, C. Bourdelle, F. Crisanti, C. Giroud, J. Hobirk, F. Imbeaux, I. Nunes, EU-ITM ITER Scenario Modelling group and JET EFDA contributors Turbulent Transport
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 informationEdge Impurity Dynamics During an ELM Cycle in DIII D
Edge Impurity Dynamics During an ELM Cycle in by M.R. Wade 1 in collaboration with K.H. Burrell, A.W. Leonard, T.H. Osborne, P.B. Snyder, J.T. Hogan, 1 and D. Coster 3 1 Oak Ridge National Laboratory General
More information