post-ablation evolution of the tungsten vuv/xuv spectra in...

1
Introduction We suggest that the introduction of tungsten into fusion plasmas significantly changes their behaviour due to considerable cooling of the pedestal. Study of this hypothesis depends on understanding both the transport of tungsten in the fusion plasma and its radiative properties. Experiment overview The presence of sawteeth can dominate the transport and introduce non-local T e effects [1], so in this experiment tungsten is ablated into JET L-mode, low temperature (T e < 2 keV), non-sawtoothing plasmas. Four JET pulses: #92284, #92285, #92286, #92293, ablation at 3.1 s in octant 6 #92286 pulse data: B = 3T, I = 1.55 MA at 3.1 s, 1.7 MA from 3.5 s. NBI 1 MW, Z eff ~1.3 (no visible change due to ablation) Measurement setup Four spectrometers used: Vertical l-o-s 15–140 nm KT7/1 14–45 nm KT7/2 4–7 nm, KT7/3 Horizontal l-o-s 10–104 nm, KT2 Results – spectroscopy Spectra preparation - average of five frames before ablation subtracted to show only ablated tungsten features Evolution: VUV – maximum amplitude just after ablation, decay: horizontal – time constant 160 ms vertical – t c 200 ms XUV features: increase – first 100 ms decay – t c 470 ms Results – bolometry reconstruction Just after ablation radiation localised along the separatrix, originates from many species, esp. deuterium after ~100 ms (time of increase of the XUV feature) maximum radiation in the plasma core (apart from divertor) and then start of the decay, t c ~500 ms Modelling- overview Transport is modelled with neoclassical edge convection and turbulent core diffusion (resulting core n W insensitive to model selection) Measurement of line-integrated radiation is used to constrain the transport model Synthetic modelling of the total radiated power and spectral radiance can be used to assess the W cooling function L W Radiation modelling W cooling function between 0.8 – 2 keV from baseline (Pütterich [2]) atomic data is smoothly increasing with temperature With updates to the total line power (PLT) coefficients [3] and an (over-)estimate of the dielectronic rate (DR) coefficient enhancements [4], the cooling function can feature structure within this temperature region By using inverted bolometry measurements to determine P rad /n e across the plasma radius, L W n W can be assessed as a function of temperature n W is modelled using JETTO [5] to understand the contribution of n W to the P rad /n e measurement Spectroscopic features – discussion XUV feature comes from W 27+ - 35+ [2] VUV features are emitted by W 16+ - 30+ ions and are sensitive to the configurations included in the collisional radiative population model and to the ionization balance Modelling of W ion emission between Z=23-26 has been calculated previously [2, 6] Work is ongoing to assess the contribution to the emission feature from W 18+ - 22+ Summary We have introduced tungsten by ablation in a L-mode low temperature (<2 keV) plasma without sawteeth and performed measurements and data treatment to obtain VUV and XUV spectra of ablated tungsten and reconstruction of the spatial distribution of the radiated power Modelling results of P rad /n e suggest that the current baseline set of ionisation, recombination, and line power coefficients do not recover the experimental results between T e =0.6 - 1.3 keV An enhanced baseline data set better recovers the flatter radiation profile at lower temperature but work is ongoing to deliver new dielectronic rate coefficients Spectral fitting is an ongoing work for W 18+ - 22+ with data already available for W 23+ - 26+ A finalised set of atomic data producing L w should provide consistency over both the total radiation and VUV emission profile. References [1] E. R. Solano et al, 43rd EPS Conference on Plasma Physics, P2.005 (2016); [2] T. Pütterich et al, Plasma Physics and Controlled Fusion 50, 085016 (2008); [3] S. Henderson et al., Plasma Physics and Controlled Fusion 59, 055010 (2017); [4] N R Badnell et al, Phys. Rev. A, 85, 052716 (2012); ; [5] G. Cenacchi and A. Taroni, Rapporto ENEA RT/TIB 88(5), (1988) [6] T. Pütterich et al, AIP Conf. Proc. 1545, 132 (2013). E. Pawelec 1 , S. S. Henderson 2 , M. O’Mullane 2 , E. R. Solano 3 , A. Huber 4 , A. Baciero 3 , J. Boom 5 , I. Coffey 6 , J. Flanagan 5 , S. Jachmich 4 , K. Lawson 5 , M. Maslov 2 , L. Meneses 7 , S. Menmuir 5 , T. Pütterich 8 , S. Schmuck 5 , G. Sergienko 3 , M. Stamp 5 , H. Summers 2 and JET Contributors* EUROfusion Consortium, JET, Culham Science Centre, Abingdon, OX14 3DB, UK 1 Institute of Physics, Opole University, Oleska 48, 45-052 Opole, Poland; 5 UKAEA-EURATOM Association, Culham Science Centre, Abingdon, OX14 3DB, UK; 2 Department of Physics, University of Strathclyde, Glasgow G4 0NG, UK; 6 Queens University Belfast, BT71NN, UK; 3 Laboratorio Nacional de Fusión, CIEMAT, Avda Complutense 40, 280040 Madrid, Spain; 7 Instituto de Plasmas e Fusao Nuclear, Universidade de Lisboa, Portugal; 4 Forschungszentrum Jülich GmbH, 52425 Jülich, Germany; 8 Max-Planck-Institut für Plasmaphysik, D-85748 Garching, Germany * See Litaudon et al, Overview of the JET results in support to ITER, accepted for publication in Nuclear Fusion Post-ablation evolution of the tungsten VUV/XUV spectra in JET 10 15 20 25 30 35 40 45 0 2 4 6 8 0 1 2 3 4 5 10 15 20 25 30 35 40 45 0 1 2 3 0 1 2 3 4 10 15 20 25 30 35 40 45 0 1 2 0 1 2 10 15 20 25 30 35 40 45 0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 Wavelength (nm) 3.1 s 3.15 s KT2 intensity (a.u.) Wavelength (nm) KT7/1 intensity (a. u.) 3.2 s Wavelength (nm) 3.3 s Wavelength (nm) 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 -1 0 1 2 3 4 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 -1 0 1 2 3 4 5 Intensity (a.u.) Wavelength (nm) 3.15 s 3.05 s 3.3 s 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 -1 0 1 2 3 3.5 s 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 -1 0 1 2 3 3.75 s VUV XUV Temporal evolution of VUV and XUV spectra. Left – VUV spectra from vertical (black) and horizontal (magenta) line of sight. Wavelength dependent calibration sensitivities have been applied to each spectrometer. Right – XUV spectra from vertical line of sight. The sensitivity response is expected to be flat as a function of wavelength for this instrument. n e (ρ,t)/10 19 [m -3 ] T e (ρ,t) [keV] P rad (ρ,t) [kW/m -3 ] P rad /n e n e T e P rad 0.2 0 0 0.38 0.2 0.2 0.55 0.4 0.4 0.73 0.6 0.6 0.91 0.8 0.8 1.09 1.0 1.0 1.27 1.2 1.2 1.44 1.4 1.4 1.62 1.6 1.6 1.8 1.8 1.8 3.1 3.2 3.3 3.4 Time (s) 2.972 3.072 3.122 3.172 3.222 3.272 3.322 3.372 3.422 3.472 Plasma parameters for JPL #92286 Ablation time Ablation time Note that P rad /n e is proportional to the cooling factor, n W L W (T e ), and shows the penetration of tungsten to the plasma Reconstruction of the radiated power from the bolometer, averaged over 50 ms time frames. The ablation is at 3.1s Bolometry VUV/XUV spectroscopy W insertion point KT2 line of sight (VUV horizontal) KT7 lines of sight (VUV and XUV vertical) Horizontal bolometer, all spectrometers and ablation insert in octant 6, vertical bolometer in octant 3 Up: Modeled n W profiles for different transport coefficients Down: comparison of the synthetic and experiment radiation along the bolometer LoS The spectrometer (blue) and bolometer (red) line-of-sight used for synthetic diagnostic measurements Different tungsten transport coefficients used in JETTO simulations

Upload: others

Post on 15-Aug-2020

0 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Post-ablation evolution of the tungsten VUV/XUV spectra in JETfusionwiki.ciemat.es/fusionwiki/images/e/eb/Pawelec-poster-EPS201… · Summary We have introduced tungsten by ablation

Introduction

We suggest that the introduction of tungsten into fusion plasmas significantly changes their behaviour due to considerable cooling of the pedestal. Study of this hypothesis depends on understanding both the transport of tungsten in the fusion plasma and its radiative properties.

Experiment overview

The presence of sawteeth can dominate the transport and introduce non-local T

e

effects [1], so in this experiment tungsten is ablated into JET L-mode, low temperature (T

e < 2 keV), non-sawtoothing plasmas.

• Four JET pulses: #92284, #92285, #92286, #92293, ablation at 3.1 s in octant 6

• #92286 pulse data: – B = 3T, I = 1.55 MA at 3.1 s,

1.7 MA from 3.5 s.

– NBI 1 MW, Zeff ~1.3 (no visible change due to ablation)

Measurement setup

Four spectrometers used:

● Vertical l-o-s – 15–140 nm KT7/1– 14–45 nm KT7/2– 4–7 nm, KT7/3

● Horizontal l-o-s– 10–104 nm, KT2

Results – spectroscopy

● Spectra preparation - average of five framesbefore ablation subtractedto show only ablatedtungsten features

Evolution:

● VUV – maximum amplitude just after ablation, decay:– horizontal – time

constant 160 ms

– vertical – tc 200 ms

● XUV features:– increase – first 100 ms

– decay – tc 470 ms

Results – bolometry reconstruction

● Just after ablation radiation localised along the separatrix,originates from many species,esp. deuterium

● after ~100 ms (time of increaseof the XUV feature) maximumradiation in the plasma core(apart from divertor) and thenstart of the decay, t

c ~500 ms

Modelling- overview● Transport is modelled with

neoclassical edge convection and turbulent core diffusion (resulting core n

W insensitive

to model selection)● Measurement of line-integrated

radiation is used to constrain the transport model

● Synthetic modelling of the total radiated power and spectral radiance can be used to assess the W cooling function L

W

Radiation modelling● W cooling function between 0.8 – 2 keV

from baseline (Pütterich [2]) atomic data is smoothly increasing with temperature

● With updates to the total line power (PLT) coefficients [3] and an (over-)estimate of the dielectronic rate (DR) coefficient enhancements [4], the cooling function can feature structure within this temperature region

● By using inverted bolometry measurements to determine P

rad/n

e across the plasma radius,

LWn

W can be assessed as a function of temperature

● nW is modelled using JETTO [5] to understand

the contribution of nW to the P

rad/n

e measurement

Spectroscopic features – discussion● XUV feature comes from W27+ - 35+ [2]● VUV features are emitted by W16+ - 30+ ions

and are sensitive to the configurations included in the collisional radiative population model and to the ionization balance

● Modelling of W ion emission between Z=23-26 has been calculated previously [2, 6]

● Work is ongoing to assess the contribution to the emission feature from W18+ - 22+

Summary● We have introduced tungsten by ablation in a L-mode low temperature (<2 keV)

plasma without sawteeth and performed measurements and data treatment to obtain VUV and XUV spectra of ablated tungsten and reconstruction of the spatial distribution of the radiated power

● Modelling results of Prad

/ne suggest that the current baseline set of ionisation,

recombination, and line power coefficients do not recover the experimental results between T

e=0.6 - 1.3 keV

● An enhanced baseline data set better recovers the flatter radiation profile at lower temperature but work is ongoing to deliver new dielectronic rate coefficients

● Spectral fitting is an ongoing work for W18+ - 22+ with data already available for W23+ - 26+

● A finalised set of atomic data producing Lw should provide consistency over both the

total radiation and VUV emission profile.

References[1] E. R. Solano et al, 43rd EPS Conference on Plasma Physics, P2.005 (2016); [2] T. Pütterich et al, Plasma Physics and Controlled Fusion 50, 085016 (2008); [3] S. Henderson et al., Plasma Physics and Controlled Fusion 59, 055010 (2017); [4] N R Badnell et al, Phys. Rev. A, 85, 052716 (2012); ; [5] G. Cenacchi and A. Taroni, Rapporto ENEA RT/TIB 88(5), (1988) [6] T. Pütterich et al, AIP Conf. Proc. 1545, 132 (2013).

E. Pawelec1, S. S. Henderson2, M. O’Mullane2, E. R. Solano3, A. Huber4, A. Baciero3, J. Boom5, I. Coffey6, J. Flanagan5, S. Jachmich4, K. Lawson5, M. Maslov2, L. Meneses7, S. Menmuir5, T. Pütterich8, S. Schmuck5, G. Sergienko3, M. Stamp5, H. Summers2

and JET Contributors*

EUROfusion Consortium, JET, Culham Science Centre, Abingdon, OX14 3DB, UK1Institute of Physics, Opole University, Oleska 48, 45-052 Opole, Poland; 5UKAEA-EURATOM Association, Culham Science Centre, Abingdon, OX14 3DB, UK; 2Department of Physics, University of Strathclyde, Glasgow G4 0NG, UK; 6Queens University Belfast, BT71NN, UK;3Laboratorio Nacional de Fusión, CIEMAT, Avda Complutense 40, 280040 Madrid, Spain; 7Instituto de Plasmas e Fusao Nuclear, Universidade de Lisboa, Portugal;4Forschungszentrum Jülich GmbH, 52425 Jülich, Germany; 8Max-Planck-Institut für Plasmaphysik, D-85748 Garching, Germany* See Litaudon et al, Overview of the JET results in support to ITER, accepted for publication in Nuclear Fusion

Post-ablation evolution of the tungsten VUV/XUV spectra in JET

10 15 20 25 30 35 40 4502468

012345

10 15 20 25 30 35 40 450

1

2

3

01234

10 15 20 25 30 35 40 450

1

2

0

1

2

10 15 20 25 30 35 40 45

0.0

0.5

1.0

1.5

0.00.51.01.5

Wavelength (nm)

3.1 s

3.15 sKT2

inte

nsity

(a.

u.)

Wavelength (nm)

KT7/

1 in

tens

ity (a.

u.)

3.2 s

Wavelength (nm)

3.3 s

Wavelength (nm)

4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5-101234

4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5-1012345

Inte

nsity

(a.

u.)

Wavelength (nm)

3.15 s 3.05 s

3.3 s

4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5-1

0

1

2

3

3.5 s

4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5

-1

0

1

2

3 3.75 s

VUV XUV

Temporal evolution of VUV and XUV spectra. Left – VUV spectra from vertical (black) and horizontal (magenta) line of sight. Wavelength dependent calibration sensitivities have been applied to each spectrometer. Right – XUV spectra from vertical line of sight. The sensitivity response is expected to be flat as a function of wavelength for this instrument.

ne(ρ,t)/1019 [m-3]

Te(ρ,t) [keV]

Prad(ρ,t) [kW/m-3]

Prad/ne

ne

Te

Prad

0.2 0 0

0.38 0.2 0.2

0.55 0.4 0.4

0.73 0.6 0.6

0.91 0.8 0.8

1.09 1.0 1.0

1.27 1.2 1.2

1.44 1.4 1.4

1.62 1.6 1.6

1.8 1.8 1.8

3.1 3.2 3.3 3.4

Time (s)

2.972 3.072 3.122 3.172 3.222

3.272 3.322 3.372 3.422 3.472

Plasma parameters for JPL #92286

Ablation time

Ablation time

Note that Prad/ne is proportional to the cooling factor, nW LW(Te), and shows the penetration of tungsten to the plasma

Reconstruction of the radiated power from the bolometer, averaged over 50 ms time frames. The ablation is at 3.1s

Bolometry VUV/XUV spectroscopy

W insertionpoint

KT2 line of sight (VUV horizontal)

KT7 lines of sight (VUV and XUV vertical)

Horizontal bolometer, all spectrometers andablation insert in octant 6, vertical bolometerin octant 3

Up: Modeled nW profiles for different

transport coefficientsDown: comparison of the synthetic and experiment radiation along the bolometer LoS

The spectrometer (blue) and bolometer (red) line-of-sight used for synthetic diagnostic measurements

Different tungsten transport coefficients used in JETTO simulations