Search for a Critical Electron Temperature Gradient in DIII-D L-mode Discharges

Similar documents
The geometric construction of Ewald sphere and Bragg condition:

LC transfer of energy between the driving source and the circuit will be a maximum.

Supplementary Figure 1. Circular parallel lamellae grain size as a function of annealing time at 250 C. Error bars represent the 2σ uncertainty in

Progress of HMGC Nonlinear Simulation

Stellar Structure and Evolution

= 4 3 π( m) 3 (5480 kg m 3 ) = kg.

Pulse Neutron Neutron (PNN) tool logging for porosity Some theoretical aspects

Basic Bridge Circuits

Chapter 3 Optical Systems with Annular Pupils

Profile Formation and Sustainment of Autonomous Tokamak Plasma with Current Hole Configuration

Tests of Profile Stiffness Using Modulated Electron Cyclotron Heating

Numerical solution of diffusion mass transfer model in adsorption systems. Prof. Nina Paula Gonçalves Salau, D.Sc.

OSCILLATIONS AND GRAVITATION

15 B1 1. Figure 1. At what speed would the car have to travel for resonant oscillations to occur? Comment on your answer.

On the Sun s Electric-Field

11) A thin, uniform rod of mass M is supported by two vertical strings, as shown below.

Determining solar characteristics using planetary data

EXAM NMR (8N090) November , am

Electron Density Distribution in HSX

Experiment I Voltage Variation and Control

Magnetic Fluctuation-Induced Particle Transport. and Zonal Flow Generation in MST

Physics 221 Lecture 41 Nonlinear Absorption and Refraction

Section 11. Timescales Radiation transport in stars

B. Spherical Wave Propagation

Multipole Radiation. February 29, The electromagnetic field of an isolated, oscillating source

MATH 415, WEEK 3: Parameter-Dependence and Bifurcations

1) Emits radiation at the maximum intensity possible for every wavelength. 2) Completely absorbs all incident radiation (hence the term black ).

LINEAR AND NONLINEAR ANALYSES OF A WIND-TUNNEL BALANCE

Introduction to Nuclear Forces

A scaling-up methodology for co-rotating twin-screw extruders

THEORY AND THEORY-BASED MODELS FOR THE PEDESTAL, EDGE STABILITY AND ELMs IN TOKAMAKS

you of a spring. The potential energy for a spring is given by the parabola U( x)

Plasma heating in reversed field pinches at the fundamental ion cyclotron frequency

Lecture 24 Stability of Molecular Clouds

Gravitational Radiation from Oscillating Gravitational Dipole

Extra notes for circular motion: Circular motion : v keeps changing, maybe both speed and

Phys101 Lectures 30, 31. Wave Motion

Ultrafast effects in 3D Metamaterials

ASTR415: Problem Set #6

Black Body Radiation and Radiometric Parameters:

Electron Transport Stiffness and Heat Pulse Propagation on DIII-D

ATMO 551a Fall 08. Diffusion

Graphs of Sine and Cosine Functions

SAMPLING DELAY AND BACKLASH IN BALANCING SYSTEMS

ME 3600 Control Systems Frequency Domain Analysis

Coupled Electromagnetic and Heat Transfer Simulations for RF Applicator Design for Efficient Heating of Materials

1) Consider an object of a parabolic shape with rotational symmetry z

ANTENNAS. Vector and Scalar Potentials. Maxwell's Equations. D = εe. For a linear, homogeneous, isotropic medium µ and ε are contant.

Numerical Integration

Collision Frequency of Adsorbed Particles

INFLUENCE OF GROUND INHOMOGENEITY ON WIND INDUCED GROUND VIBRATIONS. Abstract

Impedance-based Biosensors. Pittsburgh, PA 15213, U.S.A.

GENERAL RELATIVITY: THE GEODESICS OF THE SCHWARZSCHILD METRIC

CHAPTER 3. Section 1. Modeling Population Growth

EM Boundary Value Problems

Evolution of ETG-mode scale turbulence spectra and anomalous electron transport in dynamic experiments at FT-2 tokamak

V G. In this class, we will look at a possible hypothesis for way the time dependence is t

Field emission of Electrons from Negatively Charged Cylindrical Particles with Nonlinear Screening in a Dusty Plasma

Modeling of Non-stationary Gas Flow in Annular Nozzle

FEMA 451B Topic 3 Handouts Structural Dynamics 1

Astronomy 111, Fall October 2011

Chapter 13 Gravitation

Magnetohydrodynamics (MHD) I

Conventional Paper-I (a) Explain the concept of gradient. Determine the gradient of the given field: ( )

Solutions. V in = ρ 0. r 2 + a r 2 + b, where a and b are constants. The potential at the center of the atom has to be finite, so a = 0. r 2 + b.

1 kev Plasma. Theoretical model of a collisionally ionized plasma kt=1 kev with solar abundances The lines are 'narrow' Notice dynamic range of 10 5

One-Dimensional, Steady-State. State Conduction with Thermal Energy Generation

Sensor and Simulation Notes. Note 525. Oct Lens Design for a Prolate-Spheroidal Impulse radiating Antenna (IRA)

Core Fluctuations and Current Profile Dynamics in the MST Reversed-Field Pinch

Analysis of high speed machining center spindle dynamic unit structure performance Yuan guowei

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

General Railgun Function

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s

Modeling of Noble Gas Injection into Tokamak Plasmas

Objective Notes Summary

Molecular dynamics simulation of ultrafast laser ablation of fused silica

Mathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden

GA A22703 POLOIDAL GRADIENTS, DIVERTOR DETACHMENT AND MARFES

FI 2201 Electromagnetism

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 4: Toroidal Equilibrium and Radial Pressure Balance

Tutorial Exercises: Central Forces

Do not turn over until you are told to do so by the Invigilator.

Problem 1. Part b. Part a. Wayne Witzke ProblemSet #1 PHY 361. Calculate x, the expected value of x, defined by

Physics 181. Assignment 4

Central Force Motion

LECTURER: PM DR MAZLAN ABDUL WAHID PM Dr Mazlan Abdul Wahid

, and the curve BC is symmetrical. Find also the horizontal force in x-direction on one side of the body. h C

Rydberg-Rydberg Interactions

2017 US/EU Transport Task Force Workshop April 26 th 2017 Williamsburg, VA

Gauss s Law: Circuits

ECE 6340 Intermediate EM Waves. Fall Prof. David R. Jackson Dept. of ECE. Notes 4

7.2. Coulomb s Law. The Electric Force

Mobility of atoms and diffusion. Einstein relation.

10. Force is inversely proportional to distance between the centers squared. R 4 = F 16 E 11.

Conservative Averaging Method and its Application for One Heat Conduction Problem

Physics 161 Fall 2011 Extra Credit 2 Investigating Black Holes - Solutions The Following is Worth 50 Points!!!

PES 3950/PHYS 6950: Homework Assignment 6

24.2 Naive Backtracking

20th Century Atomic Theory - Hydrogen Atom

Magnetic Field and Inductance Calculations in Theta-Pinch and Z-Pinch Geometries

Transcription:

Seach fo a Citical Electon Tempeatue Gadient in DIII-D L-mode Dischages.C. Deoo, S. Ciant 1, T.C. Luce, A. Manini, C.C. Petty, F. Ryte, M.E. Austin 3, D.R. ake, K.W. Gentle 3, C.M. Geenfield,.E. Kinsey and G.M. Staeble DIII-D National Fusion Facility 1 Instituto di Fisica del Plasma - CNR IPP-Gaching 3 Univesity of Texas at Austin Lehigh Univesity th IAEA Fusion Enegy Confeence Novembe 1-, Vilamoua, Potugal DIII D NATIONAL FUSION FACILITY SAN DIEGO Deoo IAEA

Intoduction Pevious expeiments on DIII-D L-mode and H-mode dischages showed no evidence of nonlinea behavio in the T e esponse to modulated ECH. Othe tokamaks, notably ASDEX-U and FTU, have found that a tanspot model based on a citical tempeatue gadient scale length agees well with both steady state powe balance analysis and the esponse fom modulated ECH. Two expeiments on DIII-D have been pefomed with the pupose of seaching fo evidence of χ e dependent on a citical value of T e. ECH swing expeiment (based on FTU expeiments) systematic heat flux scan (based on AUG expeiments) DIII D NATIONAL FUSION FACILITY SAN DIEGO Deoo IAEA

ECH Swing Expeiment Fist caied out in FTU (Ciant, IAEA 3). Citical tempeatue gadient model χ e = χ o + f(t e )( T e HP (q/n) k cit )H k χ e T e T e = χ o + f(t e )( T e T e k cit )H k If the dive is stong enough to cause the local T e to go fom below a citical value to above it then the nonlinea change in χ e HP will poduce a nonlinea esponse in T e. Utilize 1 o out-of-phase, squae wave modulation of pais of gyotons absobed at slightly diffeent adii in ode to localize the heat flux and T e modulation at constant total input powe. DIII D NATIONAL FUSION FACILITY SAN DIEGO Deoo IAEA

L MODE TARGET PLASMA Inside wall limited L mode. MA,. T 1.5-.9 x 1 19 m 3 Sawtooth fee ECH esonance at ρ ECH =.-.3 and.-.5 fmod = 5 Hz and PECH1 ~ PECH ~ 1 MW Thee heat flux conditions:.,., MW NI R (m) = 1.-1.73 a (m) =. κ = 1.5-1. ECH Resonance Laye ECH Launch 3..5 1153 DIII D NATIONAL FUSION FACILITY SAN DIEGO T e (kev). 1.5 1..5 ECH1 ECH.... 1 ρ Deoo IAEA

3 Compaison of Pofiles with P NI = MW and MW ECH deposition at ECH =.-.3 Cental T i is educed without NI Lowe density causes an incease in Zeff T e (kev) 115 MW NI 115 No NI 3 T i (kev) 3 T i 1 1 1 1 W (1 ad/sec) n e (1 19 m -3 ) 5 Z eff 3 1.... 1.... 1.... 1 Deoo IAEA

-ST e /T e in the ECH Absoption Region -ST e vaies fom.5 to 5 m -1 at =.9 whee T e is nealy constant T esponse is clealy not vey stiff as ST vaied by facto e e T aveaged ove 9 modulation peiods fom 11-1 msec on shot e 115 1 115 3.5 ECH on ECH1 on Rho.31.91.5.39 3..5. 1.5 1..5 1137 msec 1117 msec ECH1 ECH 11 111 11 113 11 115 Time (msec).15..5.3.35 Deoo IAEA

The Relative Change in ST e Is Lagest Fo The Ohmic Case The Ohmic case is most likely case to obseve nonlinea behavio Changes in ST e ae small outside the ECH deposition since in this egion the ECH heat flux is essentially constant ST max -ST min ST min 1 Ohmic, ECH ~. MW, ECH ~. MW, ECH ~. 115 115 1153.... 1.... 1.... 1 Deoo IAEA

No Evidence of a Non-lineaity in c e HP is Obseved Fo Case with no NI and ECH =.-.3 o Any Othe Case The linea combination of the Fouie amplitude and phase of the Te oscillations due to ECH1 and ECH applied individually is the same as the esult with ECH1 and ECH modulated simultaneously and out of phase with espect to each othe 115 Simultaneous modulation of ECH1 and ECH 15 Calculated linea combination of ECH1 and ECH applied individually 1 5.... 1 3 1 ECH1 ECH.... 1 Deoo IAEA

Heat Flux Scan Expeiments Expeimental appoach (follow expeiments on AUG, Ryte NF 3 3) scan T e locally by vaying the electon heat flux while holding the flux at the edge (total powe) 3 constant ECH at two locations, ρ 1 and ρ vay P ECH1 and P ECH at constant 1 P ECH1 + P ECH... modulate one gyoton at ρ to ρ obtain χ HP. 1 T e (kev) Taget plasma (simila to low density ECH swing plasmas) L-mode limited on inside wall of vessel I p =. MA, T =.T, n e = 1.7x1 19 m -3, q a = 7 ealy NI to delay onset of sawteeth, then ECH only DIII D NATIONAL FUSION FACILITY SAN DIEGO Deoo IAEA

T e is Systematically Vaied etween ECH1 and ECH Locations EC powe is moved fom ECH1 to ECH shot by shot T e pofile outside ECH is constant since PECH1+PECH is held fixed T e (kev) 3 ECH1 ECH 1 / % / 1 1 / 39 7 / 73 / 1 1 ECH1 ECH.... 1 Deoo IAEA

The Electon Heat Flux is Systematically Vaied at the Half Radius The total powe in the electons is compised of P ~ 1.5 MW ECH and P ~.5 MW. OH The dominant loss is by conduction, convection is negligible. At =.5 adiation and collisional tansfe to ions is small..5 Suface Integated Electon Heat Flux (MW). 1.5 1..5..... 1 Deoo IAEA

Modulate ECH Powe on One Gyoton at Outside Location to Obtain Heat Pulse Diffusivity Modulate at Hz with 5% duty cycle fequency chosen to avoid matching sawtooth fequency Fouie analyze ECE channels Spatial deivative of pulse amplitude and phase yields an estimate of the heat pulse diffusivity χ HP Use 3 d hamonic pulse amplitude and phase to estimate χ HP DIII D NATIONAL FUSION FACILITY SAN DIEGO Deoo IAEA

No Evidence of a Citical Gadient Scale Length Was Obseved ut the Existence of One Cannot e Ruled Out A citical gadient scale length, k cit, could exist below T e = 3. m -1 at =.5. This coesponds to L Te. m o ove % of the plasma adius. If k exists it would be at values typically found nea the plasma coe, cit indicating most of the tempeatue pofile would be above it. c e 3/ (m /sec/kev 3/ ) 1 1 1 1 Heat Pulse Powe alance Data at =.5 1 3 5 7 - T e (m -1 ) c e nomed to T e 3/ to take out gyo- ohm dependence Deoo IAEA

Results Of This Expeiment Ae Consistent With ASDEX Upgade Expeiments The diffusivity levels and tends with 1/L Te ae simila ASDEX Upgade found k cit =.3 m -1 fit thei data (Ryte, NF 3 3) c e 3/ (m /sec/kev 3/ ) 1 1 1 1 AUG DIII-D H Heat Pulse P Powealance H H H H H P PP P H P Data at =.5 P 1 3 5 7 - T e (m -1 ) Deoo IAEA

TEMs Dominate Dift Wave Spectum in Region of Inteest Analysis pefomed with the linea gyokinetic stability code GKS TEMs ae unstable ove most of plasma mino adius Lagest gowth ates fo most unstable modes occu at..5 1 1 1 11.... 1 7 5 3 1.... 1 Deoo IAEA

Spectum of Wavenumbes Whee Gowth Rate Is Maximum Spectal width at =. is dk ~ 3 cm -1 Fequency at the lagest gowth ate is f ~ 1 khz 1 1 11 =. 1 3 5 7 9 1 k (cm -1 ) Deoo IAEA

ETG Modes Ae Stable In The Region Of Inteest ETG modes unstable only nea plasma edge, >. Citical gadient fo ETG modes calculated fom linea gyokinetic stability code GKS - T e (m -1 ) 5 15 1 5 11 Citical level fo ETG modes.... 1 Deoo IAEA

Thee Possible Dependencies of c HP on T e Offset linea: c HP becomes negative fo T e < m -1, implying heat flux deceases as the gadient becomes steepe (pehaps due to a heat pinch) Nonlinea: ( T e ) a, fit shows a ~ 5/3. Exponent can not be too lage in ode to emain consistent with ealie expeiments whee no obvious nonlinea behavio was obseved. Citical gadient: If a citical gadient exists it must be at T e < 3. m -1 at =.5 (abitaily chose 3. m -1 fo illustation). AUG has epoted a citical value of k cit ~.3 m -1. 1 Offset Linea Nonlinea Citical Gadient 1 1 1 3 5 7 - T e (m -1 ) 1 3 5 7 - T e (m -1 ) 1 3 5 7 - T e (m -1 ) Deoo IAEA

Summay No evidence of a citical gadient scale length dependence of χ HP was obseved but the existence of one cannot be uled out by the expeimental esults. The T e pofile was not vey stiff since the gadient was vaied by a facto of at locations in the pofile that wee above k cit, if k cit exists. If a citical gadient scale length exists, all expeimental esults to date ae consistent with the obsevations being made at values above it. The esults indicate that if k cit exists then k cit < 3. m -1 at ρ=.5 and k cit <.5 m -1 at ρ=.9 coesponding to a citical gadient scale length lage than 3% and 5% of the plasma mino adius, espectively. Models othe than one based on k cit ae also consistent with the expeimental obsevations such as an offset linea o powe law dependence of χ HP on 1/L Te. DIII D NATIONAL FUSION FACILITY SAN DIEGO Deoo IAEA