Influence of Beta, Shape and Rotation on the H-mode Pedestal Height

Similar documents
EFFECT OF EDGE NEUTRAL SOUCE PROFILE ON H-MODE PEDESTAL HEIGHT AND ELM SIZE

Development and Validation of a Predictive Model for the Pedestal Height (EPED1)

Characteristics of the H-mode H and Extrapolation to ITER

ELMs and Constraints on the H-Mode Pedestal:

Edge Rotational Shear Requirements for the Edge Harmonic Oscillation in DIII D Quiescent H mode Plasmas

GA A26891 A FIRST PRINCIPLES PREDICTIVE MODEL OF THE PEDESTAL HEIGHT AND WIDTH: DEVELOPMENT, TESTING, AND ITER OPTIMIZATION WITH THE EPED MODEL

TH/P8-4 Second Ballooning Stability Effect on H-mode Pedestal Scalings

DEPENDENCE OF THE H-MODE PEDESTAL STRUCTURE ON ASPECT RATIO

Pedestal Stability and Transport on the Alcator C-Mod Tokamak: Experiments in Support of Developing Predictive Capability

Progress in characterization of the H-mode pedestal

The EPED Pedestal Model: Extensions, Application to ELM-Suppressed Regimes, and ITER Predictions

Integrated Simulation of ELM Energy Loss Determined by Pedestal MHD and SOL Transport

MHD Analysis of the Tokamak Edge Pedestal in the Low Collisionality Regime Thoughts on the Physics of ELM-free QH and RMP Discharges

UCLA. Broadband Magnetic and Density Fluctuation Evolution Prior to First ELM in DIII-D Edge Pedestal. Presented by G. Wang a, In collaboration with

Global stabilization effect of Shafranov shift on the edge pedestal plasmas in JET and JT-60U

ELM Suppression in DIII-D Hybrid Plasmas Using n=3 Resonant Magnetic Perturbations

H-mode Pedestal Enhancement and Improved Confinement in DIII-D with Lithium Injection

GA A27433 THE EPED PEDESTAL MODEL: EXTENSIONS, APPLICATION TO ELM-SUPPRESSED REGIMES, AND ITER PREDICTIONS

ITER Predictions Using the GYRO Verified and Experimentally Validated TGLF Transport Model

Local Plasma Parameters and H-Mode Threshold in Alcator C-Mod

Towards the construction of a model to describe the inter-elm evolution of the pedestal on MAST

Predicting the Rotation Profile in ITER

Operational Phase Space of the Edge Plasma in Alcator C-Mod

H-mode performance and pedestal studies with enhanced particle control on Alcator C-Mod

Dependences of Critical Rotational Shear in DIII-D QH-mode Discharges

Scaling of divertor heat flux profile widths in DIII-D

GA A27857 IMPACT OF PLASMA RESPONSE ON RMP ELM SUPPRESSION IN DIII-D

GA A25592 STABILITY AND DYNAMICS OF THE EDGE PEDESTAL IN THE LOW COLLISIONALITY REGIME: PHYSICS MECHANISMS FOR STEADY-STATE ELM-FREE OPERATION

New bootstrap current formula valid for steep edge pedestal, and its implication to pedestal stability

Edge Impurity Dynamics During an ELM Cycle in DIII D

Summary of CDBM Joint Experiments

Enhanced Energy Confinement Discharges with L-mode-like Edge Particle Transport*

Advancing the predictive capability for pedestal structure through experiment and modeling

Progress in Modeling of ARIES ACT Plasma

Modeling of ELM Dynamics for ITER

ITER operation. Ben Dudson. 14 th March Department of Physics, University of York, Heslington, York YO10 5DD, UK

Dependence of Achievable β N on Discharge Shape and Edge Safety Factor in DIII D Steady-State Scenario Discharges

W.M. Solomon 1. Presented at the 54th Annual Meeting of the APS Division of Plasma Physics Providence, RI October 29-November 2, 2012

Scaling of H-mode Pedestal and ELM Characteristics in the JET and DIII-D Tokamaks

GA A26565 LIMITS TO H-MODE PEDESTAL PRESSURE GRADIENT IN DIII-D

PSI meeting, Aachen Germany, May 2012

HIGH PERFORMANCE EXPERIMENTS IN JT-60U REVERSED SHEAR DISCHARGES

GA A22443 STUDY OF H MODE THRESHOLD CONDITIONS IN DIII D

Characteristics of the H-Mode Pedestal and Extrapolation to ITER

Effect of Variation in Equilibrium Shape on ELMing H Mode Performance in DIII D Diverted Plasmas

Role of Magnetic Configuration and Heating Power in ITB Formation in JET.

Correlation between the edge and the internal transport barriers in JT-60U

GA A27437 SCALING OF THE DIVERTOR HEAT FLUX WIDTH IN THE DIII-D TOKAMAK

Particle transport results from collisionality scans and perturbative experiments on DIII-D

Validation of Theoretical Models of Intrinsic Torque in DIII-D and Projection to ITER by Dimensionless Scaling

Relating the L-H Power Threshold Scaling to Edge Turbulence Dynamics

GA A25229 FINAL REPORT ON THE OFES ELM MILESTONE FOR FY05

Integrated Simulation of ELM Energy Loss and Cycle in Improved H-mode Plasmas

Low-collisionality density-peaking in GYRO simulations of C-Mod plasmas

GA A23114 DEPENDENCE OF HEAT AND PARTICLE TRANSPORT ON THE RATIO OF THE ION AND ELECTRON TEMPERATURES

Development of a Systematic, Self-consistent Algorithm for the K-DEMO Steady-state Operation Scenario

GA A23403 GAS PUFF FUELED H MODE DISCHARGES WITH GOOD ENERGY CONFINEMENT ABOVE THE GREENWALD DENSITY LIMIT ON DIII D

THE DIII D PROGRAM THREE-YEAR PLAN

Evolution of the pedestal on MAST and the implications for ELM power loadings

Correlations of ELM frequency with pedestal plasma characteristics

arxiv: v1 [physics.plasm-ph] 24 Nov 2017

The H-mode pedestal structure and its role on confinement in JET with a carbon and metal wall

DIII D Research in Support of ITER

Advanced Tokamak Research in JT-60U and JT-60SA

Direct drive by cyclotron heating can explain spontaneous rotation in tokamaks

GA A26874 ITER PREDICTIONS USING THE GYRO VERIFIED AND EXPERIMENTALLY VALIDATED TGLF TRANSPORT MODEL

QTYUIOP LOCAL ANALYSIS OF CONFINEMENT AND TRANSPORT IN NEUTRAL BEAM HEATED DIII D DISCHARGES WITH NEGATIVE MAGNETIC SHEAR D.P. SCHISSEL.

GA A23736 EFFECTS OF CROSS-SECTION SHAPE ON L MODE AND H MODE ENERGY TRANSPORT

GA A22993 EFFECTS OF PLASMA SHAPE AND PROFILES ON EDGE STABILITY IN DIII D

Minimal Model Study for ELM Control by Supersonic Molecular Beam Injection and Pellet Injection

MHD stability analysis of diagnostic optimized configuration shots in JET

Studies of H Mode Plasmas Produced Directly by Pellet Injection in DIII D

Results From Initial Snowflake Divertor Physics Studies on DIII-D

Understanding Edge Harmonic Oscillation Physics Using NIMROD

ABSTRACT, POSTER LP1 12 THURSDAY 11/7/2001, APS DPP CONFERENCE, LONG BEACH. Recent Results from the Quiescent Double Barrier Regime on DIII-D

Overview of Tokamak Rotation and Momentum Transport Phenomenology and Motivations

Partially Coherent Fluctuations in Novel High Confinement Regimes of a Tokamak

Analysis and modelling of MHD instabilities in DIII-D plasmas for the ITER mission

QTYUIOP ENERGY TRANSPORT IN NEUTRAL BEAM HEATED DIII D DISCHARGES WITH NEGATIVE MAGNETIC SHEAR D.P. SCHISSEL. Presented by. for the DIII D Team*

First Observation of ELM Suppression by Magnetic Perturbations in ASDEX Upgrade and Comparison to DIII-D Matched-Shape Plasmas

OV/2-5: Overview of Alcator C-Mod Results

PROGRESS TOWARDS SUSTAINMENT OF ADVANCED TOKAMAK MODES IN DIIIÐD *

Predictive Study on High Performance Modes of Operation in HL-2A 1

High fusion performance at high T i /T e in JET-ILW baseline plasmas with high NBI heating power and low gas puffing

Critical Physics Issues for DEMO

Theory, Simulation and Modelling for the FY2011 Pedestal Milestone

Impact of High Field & High Confinement on L-mode-Edge Negative Triangularity Tokamak (NTT) Reactor

STATIONARY, HIGH BOOTSTRAP FRACTION PLASMAS IN DIII-D WITHOUT INDUCTIVE CURRENT CONTROL

Self-consistent modeling of ITER with BALDUR integrated predictive modeling code

L-H transition and pedestal studies on MAST

Characterization of neo-classical tearing modes in high-performance I- mode plasmas with ICRF mode conversion flow drive on Alcator C-Mod

THE ADVANCED TOKAMAK DIVERTOR

Dimensionless Identity Experiments in JT-60U and JET

Confinement and edge studies towards low ρ* and ν* at JET

Recent Development of LHD Experiment. O.Motojima for the LHD team National Institute for Fusion Science

Connections between Particle Transport and Turbulence Structures in the Edge and SOL of Alcator C-Mod

Extension of High-Beta Plasma Operation to Low Collisional Regime

GA A26887 ADVANCES TOWARD QH-MODE VIABILITY FOR ELM-FREE OPERATION IN ITER

Scaling of divertor plasma effectiveness for reducing target-plate heat flux

Pedestals and Fluctuations in C-Mod Enhanced D α H-modes

Transcription:

Influence of Beta, Shape and Rotation on the H-mode Pedestal Height by A.W. Leonard with R.J. Groebner, T.H. Osborne, and P.B. Snyder Presented at Forty-Ninth APS Meeting of the Division of Plasma Physics Orlando, Florida November 12 16, 27

Pedestal Height is Greatest Uncertainty in Predicting Tokamak Performance 15 1 5 2 1 6 P ped (kpa) T e,ped (kev) p Ped. Width SOL Pedestal sets boundary condition for core transport For stiff transport, core profiles strongly dependent on pedestal Pedestal uncertainty for ITER corresponds to a wide range of performance; Q=2-15 3.75.85.95 1.5 Separatrix ψ N Uncertainty

Pedestal Width is key to Pedestal Height Scaling Validated edge stability model can predict pedestal top pressure if width is known Pedestal pressure height increases with improved edge stability Pedestal pressure increase with global and stronger shaping No change with toroidal rotation Stability model utilized to extract pedestal width Greater accuracy than directly measured width Pedestal width increases with higher pedestal pressure Pedestal pressure rise double expected from stability with fixed width Pedestal models predict an ion gyro-radius, *, width scaling, yet Width correlated with pressure ( 12 ) ped No T i ( *) dependence

Validated MHD Stability Model a Tool for Extracting Pedestal Width from Pedestal Pressure DIII-D Stability Study: Fixed parameters; I p, B t, Shape, pedestal width, etc. Pedestal Te (kev) The peeling-ballooning model has been validated over a range of conditions Stability model predicts pedestal pressure for a given width Stability model can extract width from a measured pedestal pressure Pedestal top more accurately measured than gradient and width Pedestal Density (1 19 m -3 ) [P.B. Snyder, H.R. Wilson, et al., PPCF 46 (24) A131 ]

Pedestal Pressure Increases with Global Beta I p, P inj (MA, MW) Density (1 19 m -3 ) Pe,ped(kPa) 1.4 Plasma Current 1.2 1..8.6.4.2. 4 β N =β/ 3 I/ab β N 2 1 5 4 3 2 1 1 8 6 4 2 Density Pedestal Electron Pressure Injected Power/1 1 2 3 4 5 Time (msec) Power scan for N =2. to N =3. Other parameters fixed; shape, q 95, density, etc. Pedestal pressure increases with N

Pedestal Profiles Fit Just Prior to an ELM.6.4.2. 2.5 2. 1.5 1..5. 2.5 2. 1.5 1. n e (1 2 m -3 ) T e (kev) T i (kev).5..5.6.7.8.9 1. 1.1 ψ Profiles fit to data combined from the last 2% of many ELM cycles Hyperbolic tangent fit to profiles Pedestal height more certain, 5-1% Width/gradient less certain, ~2% T i pedestal is usually significantly wider than T e or n e Ion pressure width set by density width Pedestal pressure characterizes pedestal height for comparison with stability analysis

Pedestal Pressure Increases with N Pressure (kpa) 4 3 2 1 Total pedestal pressure increases ~5% as N increases from 2.-3. b N =2. b N =2.5 b N =2.8 b N =3..6.7.8.9 1. 1.1 ψ N Pedestal Pressure (kpa) 2. 15. 1. 5. Both electron and ion pedestal pressure increase Total Pressure Electron Pressure Ion Pressure. 1.8 2. 2.2 2.4 2.6 2.8 3. Global Beta ( β N )

Both Pressure Gradient and Width Increase with N Pressure Gradient ( dp/dψ n ) Pedestal Width (Δψ n ) 25 2 15 1 5.7.6.5.4.3.2.1 Pressure Gradient Pedestal Width Electrons Ions Ions Electrons Both gradients and widths increase with N Larger uncertainty, ~2%, for gradient and width measurements Stability constraint applied to pressure can reduce uncertainty in pedestal width scaling. 1.8 2. 2.2 2.4 2.6 2.8 3. Global Beta ( β N )

Stability Space is Mapped for Each Configuration Average Pedestal Current Density 1.4 1.2 1..8.6 Instability Growth Rate Operating Point Unstable Stable 2 4 6 8 1 12 Normalized Peak Pressure Gradient (α) γ E /ω*/2 5 4 3 2 1 MHD equilibria constructed by measured pressure profile and Sauter bootstrap model for edge current Create 2D array of equilibria with varying pedestal current and pressure Fixed pedestal width Fixed total stored energy Stability analysis of each equilibrium using ELITE, n=5-25 Stability threshold: E *2 Variation in stability threshold analyzed for varying global, shape, etc.

Pressure Gradient Limit for Fixed Width Increases with Global Stability Contours for N Variation with Fixed Width Average Pedestal Current Density 1.2 1..8.6.4 β N =2. β N =2.5 β N =3. β N =3.5 Stable β N =2. Operating Point Unstable 2 4 6 8 1 12 Normalized Peak Pressure Gradient (α) Stability dependence on calculated using fixed pedestal width but varying total stored energy Calculated gradient limit increases with due to Shafranov shift

Improved Stability at Fixed Pedestal Width Can Account for Half of Pedestal Pressure Rise with N Pedestal Pressure (kpa) 2. 18. 16. 14. 12. Pedestal Height Measured Pedestal Pressure Fixed Width Stability Limit Measured pedestal pressure proportional to global N Increase in stability limit at fixed pedestal width accounts for only half of observed increase in P ped with N Width must increase with to account for pedestal pressure increase 1. 1.5 2. 2.5 3. 3.5 4. Global Beta (β N )

Pedestal Gradient Expected to be Weakly Reduced with Increasing Width Pped Limit (kpa) Pedestal Pressure Previous numerical studies indicate width ( ) scaling of stability limit: p crit 14 or P ped 34 All other parameters fixed; Shape,, etc. Width variation of p due to radial extent of finite n modes Pedestal Width (% of Flux) [P.B. Snyder, H.R. Wilson, et al., PPCF 46 (24) A131 ]

Stability Analysis Implies Pedestal Width Increase with Pedestal Pressure (kpa) 2. 18. 16. 14. 12. Pedestal Height Measured Pedestal Pressure 1. 1.5 2. 2.5 3. 3.5 4. Global Beta ( β N ) 1/2 Dµb ped Stability Limit Fixed Width Stability Limit Pedestal pressure proportional to global N Width must increase with to account for measured pressure increase Width scaling consistent with previous statistical study 12 ped T.H. Osborne, et. al., Jour. Nucl. Mater.,266-269, 131 (1999)

Similar Dependence of P ped on global N Observed in High Advanced Tokamak Regime Pressure (kpa) Temperature (kev).7.8.9 1. 1.1 2..7.8.9 1. 1.1 High scenario with early heating Increase in N from 2.2 to 3.7, ~68% 3 2 1 1.5 1..5. Total Pressure β N =2.2 β N =3.7 Pedestal T e..7.8.9 1. 1.1 Pedestal pressure rises from 9.8 kpa to 15.5 kpa, ~58% Density (1 2 m -3 ) Temperature (kev).6.4.2..7.8.9 1. 1.1 2. 1.5 1..5 Pedestal n e Pedestal T i

Pedestal Width Scales with in Advanced Regimes Average Pedestal Current Density 1.2 1..8.6 Fixed Pedestal Width Stability Limit β N =2.2 Operating Point β N =2. β N =2.5 β N =3. β N =3.5 β N =4..4 4 6 8 1 12 14 16 Normalized Peak Pressure Gradient (α) Pedestal Pressure (kpa) 18. 16. 14. 12. 1. Pedestal Height 1/2 Dµb ped Stability Limit Measured Pedestal Pressure Fixed Width Stability Limit 8. 1.5 2. 2.5 3. 3.5 4. 4.5 Global Beta ( β N ) Pedestal stability continues to improve over wide range of 12 Pedestal width scaling,, consistent with pressure increase ped

Shape an Important Factor for Pedestal Height Shape dependence for pedestal and core confinement long observed Stronger shaping improves edge stability Elongation, triangularity, squareness, second separatrix, etc. Contribution of stability to increased pedestal for stronger shaping not generally quantified Stronger shaping leads to higher stability limit and wider pedestal

Shaping Can Produce Wide Range of Pedestal Heights 128198 12825 Pressure (kpa) Temperature (kev) Constant conditions ELMing H-mode N =2. q 95 =4.6 Ip=1.1 MA Bt= 1.9T, ±1% to match q95 25 2 15 1 5.7.8.9 1. 1.1 2.5 2. 1.5 1..5. Total Pressure <δ>=.25 <δ>=.48 Pedestal T e.7.8.9 1. 1.1 Density (1 2 m -3 ) Temperature (kev).5.4.3.2.1..7.8.9 1. 1.1 2.5 2. 1.5 1. Pedestal n e Pedestal T i.5..7.8.9 1. 1.1 < >=.25 : P ped =8.1 kpa < >=.48 : P ped =14. kpa Global Energy Confinement: ~5% higher

Pedestal Width Increases with Stronger Shaping Average Pedestal Current Density 1.2 1..8.6.4 Stability Contours for Fixed Width <δ>=.25 <δ>=.48 <δ>=.25 Operating Point Low δ High δ 2 4 6 8 1 Normalized Peak Pressure Gradient (α) Stronger shaping raises pedestal pressure 72% Stability limit increases ~33% assuming fixed width Width increase of 31% for 12 scaling yields ped Expected increase in P ped of 75%, in agreement with experiment

Moderate Changes to ITER Shape Can Significantly Affect Pedestal Squareness scan about the proposed ITER shape kpa kev 25 2 15 1 5 Measured Pedestal Profiles in ELMing H-mode..85.9.95 1. 1.5.85.9.95 1. 1.5 2. 1.5 1. Total Pressure ζ=.32 ζ=.35 ζ=.41 Pedestal Te kev 1 2 m -3.6.4.2 2. 1.5 1. Pedestal ne Pedestal Ti.5.5 121559 121562 121567..85.9.95 1. 1.5..85.9.95 1. 1.5

Pedestal Width Decreases with Squareness Pedestal Pressure (kpa) 16. 14. 12. 1. Pedestal Height versus Squareness Measured Pedestal Pressure Fixed Width Limit Width ~ β Ped 1/2 Limit 8..3.32.34.36.38.4.42 Squareness (z) Pedestal pressure increases ~5% from highest to lowest squareness Improved stability for fixed width can account for half of pedestal pressure increase 12 ped Pedestal width consistent with pressure variation

Toroidal Rotation does not Significantly Affect Pedestal Height km/s 6 4 2 Rotation Scan in Hybrid Configuration Carbon Toroidal Rotation Cntr 6.5 MW 1 Cntr 7.8 MW 2 Cntr 9.2 MW kpa 2 15 1 5 Total Pressure -2.8.85.9.95 1. 1.5 Ratio of co to counter NBI to control rotation While central confinement degraded with counter injection pedestal pressure not significantly affected Previous stability studies in this parameter regime indicate edge stability limit not modified by toroidal rotation [P.B. Snyder, et al., Nucl. Fusion 47 (27) 961 ].8.85.9.95 1. 1.5

Implications for Pedestal Width Scaling Significant gyro-radius dependence expected if width set where ExB velocity shearing exceeds driftwave growth rate.5 2. [ f( S,,q)gZ ( eff )h( * )] Complicated dependence on other dimensionless parameters 12 ped Possible pedestal width scaling mechanisms: ped correlation with *: Improved stability leads to higher T i ped correlation with other parameters: S,, etc.

Gyro-Radius Scaling Not Yet Observed Pressure (kpa) 3 2 1.7.8.9 1. 1.1 2. 1.5 1..5..7.8.9 1. 1.1 Temperature (kev)2.5 Total Pressure Low Density High Density Pedestal T e Density (1 2 m -3 ).8.6.4.2..7.8.9 1. 1.1 2. 1.5 1..5..7.8.9 1. 1.1 Temperature (kev)2.5 Pedestal n e Pedestal T i * reduced, but pedestal not narrower Other dimensionless parameters change; collisionality, shear, etc. The * pedestal width scaling is critical If *, ITER likely not to achieve its mission Current experiments cannot match ITER * Experiments in DIII-D, JT-6U and C-Mod have not found significant * scaling A 12 pedestal scaling ped consistent with single parameter scans in DIII-D and JT-6U T.H. Osborne, et. al., Jour. Nucl. Mater.,266-269, 131 (1999)

Summary Pedestal pressure increases with improved edge stability Increase twice that expected for fixed pedestal width 12 Pedestal width increase, ped, consistent with stability analysis and observed pressure variation Pedestal height and width not dependent on toroidal rotation Consistent with stability analysis Future work to isolate underlying physics of width scaling Universality of 1/2 ped scaling for pedestal width Ion gyro-radius, * Other dimensionless parameters; S,, q, etc.