Tungsten control in type-I ELMy H-mode plasmas on EAST

Zong Xu ? Zhen-Wei Wu ? Ling Zhang ? Yun-Xing Cheng,4 ?Wei Gao ? Xiao-Dong Lin ? Xiang Gao, ? Ying-Jie Chen ? Yue-Heng Huang ?Juan Huang ? Qing Zang, ? Yin-Xian Jie, ? Yan-Min Duan ?Ying-Ying Li ? Hai-Qing Liu ? Guo-Sheng Xu ? EAST team

Abstract The first experimental investigation of the tungsten behavior in ELMy H-mode plasmas with co-/counter neutral beam injection (NBI) and unfavorable/favorable Bt was performed on EAST.Tungsten was found to accumulate easily in ELMy H-mode plasma with co-NBI heating and unfavorable Bt.Thus,in this case the tungsten concentration can exceed 10-4,resulting in degradation of the plasma confinement and periodic H—L transitions.To reduce the tungsten concentration in steadystate type-I ELMy H-mode operation,counter-NBI is applied to modify the density and temperature and brake the plasma toroidal rotation.The applied counter-NBI decreases the PHZ + Er inward pinch velocity and reverses the direction of neoclassical inward convection,thus decreasing the tungsten concentration from～7 × 10-5 to～2 × 10-5 in type-I ELMy H-mode plasma with favorable Bt.A comparison of the effects of different Bt directions on the tungsten behavior also shows that favorable Bt is beneficial for reducing the tungsten concentration in the core plasma.These results imply that counter-NBI with favorable Bt can effectively prevent tungsten accumulation and expand the operating window for exploring steady-state type-I ELMy H-mode operation of EAST.

Keywords Tungsten accumulation · Counter-NBI ·Favorable Bt · ELMy H-mode · EAST

1 Introduction

Tungsten is a good candidate for plasma-facing materials in future fusion devices such as ITER and CFETR.Its intrinsic benefits include a high melting point,a low sputtering yield,and a low tritium retention.However,tungsten accumulation in the core plasma can result in serious radiation losses,the degradation of the overall energy confinement [1],magnetohydrodynamic (MHD)instabilities [2],or even disruptions [3—5].For ITER,the tungsten concentration (CW=nW/ne,wherenWandneare the tungsten and electron densities,respectively) must be less than 10-4to avoid unacceptable radiative cooling [6].It is therefore important to avoid tungsten accumulation in the high-confinement regime in devices with tungsten divertor plates,where tungsten atoms can be sputtered from the divertor plates by giant edge-localized modes (ELMs)[7,8].Methods of reducing the stored energy radiated by tungsten in the central plasma are a concern for the economical operation of future commercial reactors [9].

Previous investigations have shown that core impurity accumulation is driven mainly by neoclassical inward convection[10],and the negative radial electric field(Er)is responsible for inward impurity transport [11,12].In addition,the inward pinch of high-Zimpurities (PHZ)attributed to the atomic processes of ionization and recombination also plays a role in core tungsten accumulation in a toroidally rotating plasma [12—14].Thus,to suppress and avoid tungsten accumulation,these inward pinch velocities must be reduced and well controlled.Theoretical predictions and experimental observations have shown that neoclassical convection can be controlled by changing the density and temperature gradients[15—17],and the PHZ andErpinch depend strongly on the toroidal rotation velocity of the background plasma[12].Turbulent transport due to,for example,the ion temperature gradient(ITG) mode and trapped electron mode (TEM),can also affect tungsten transport [18—20],but it is much weaker than neoclassical transport in the central region of a toroidally rotating plasma,as has been found in JET and the ASDEX Upgrade [21—23].Turbulent transport can be enhanced by decreasing the ion-to-electron temperature ratio (Ti/Te) [24,25] so as to increase the impurity diffusivity and suppress impurity accumulation.Neutral beam injection (NBI),which can alter the density and temperature and change the toroidal rotation velocity,will be used to study the tungsten behavior in the core of an EAST type-I ELMy H-mode plasma,because the ELMy H-mode with a low ELM frequency (fELM) of dozens of hertz has inevitably been found to result in strong impurity accumulation in EAST [26—28].Here,the effects of neoclassical transport [10] and the PHZ andErpinch on the tungsten behavior in the type-I ELMy H-mode regime are studied.In addition,the effect of turbulent transport on tungsten transport is briefly discussed.

The remainder of the paper is structured as follows.Section 2 describes the NBI systems and relevant diagnostics in EAST,and Sect.3 reports the experimental investigation of the tungsten behavior in ELMy H-mode plasmas with co-/counter-NBI and different toroidal field(TF)directions.In Sect.4,a summary and discussion are presented.

2 Experimental setup

This work is based on the EAST tokamak device.EAST is a full superconducting device that was designed to operate with major radiiRof 1.7—1.9 m,minor radiiaof 0.4—0.45 m,a maximum TF strengthBtof 3.5 T,a maximum plasma currentIpof 1 MA,an elongationkof 1.2—2,and a triangularity δ of ＜0.8.Since the 2014 campaign,EAST has been upgraded to facilitate progress in achieving long-pulse stable high-performance plasma operation[pulse length ＞1000 s,central electron temperatureTe0-＞10 keV,normalized beta βN＞2,andne～nGW,wherenGW=Ip/πa2(1020m-3)].For these purposes,several auxiliary heating systems have been developed to provide effective heating and current drive for EAST [29,30].As the auxiliary heating power is increased,the materials of plasma-facing components(PFCs)are replaced at different positions of the inner wall to improve the heat exhaust ability [31].The main PFC materials in EAST are molybdenum tiles mounted on the first wall of the vacuum chamber,SiC-coated graphite tiles fixed in the lower divertor,and W/Cu mono-blocks installed in the upper divertor [32].

Two sets of NBI systems have been installed on EAST since the 2014 and 2015 campaigns [33,34].The NBI system in EAST consists of two sets of neutral beam injectors.The maximum beam power of each injector is 4 MW,and the maximum beam energy ranges from 50 to 80 keV.One injector has two ion sources which with an independent beam channel and can be operated independently.The layout of the EAST NBI system is shown in Fig.1.For anticlockwiseIp,the tangential direction of each injector (the co-Ipdirection is denoted as co-NBI,whereas the counter-Ipdirection is denoted as counter-NBI)and the unfavorable direction ofBt(counter-Ipdirection)are shown in the figure.

Fig.1 (Color online) Directions of plasma current Ip and toroidal field Bt,and layout of two sets of NBI systems and EUV spectrometers on EAST

To investigate the tungsten behavior in core plasmas with different NBI andBtdirections,three newly developed extreme ultraviolet (EUV) spectrometers(EUV_Short,EUV_Long,and EUV_Long2 in Fig.1) will be used.EUV_Short and EUV_Long [35,36] operate in wavelength ranges of 20—50 A° and 20—500 A°,respectively,and both have good spectral resolution for measuring the chord-integrated intensity of line emission from tungsten.Other relevant diagnostics include a Thomson scattering system for observing the profiles of the electron temperature and density,a polarimeter—interferometer for monitoring the profiles of the electron density and safety factor,a fast CCD camera for detecting the plasma position[37],and charge exchange recombination spectroscopy for monitoring the ion temperature profile and toroidal rotation velocity of the background plasma.

3 Experimental results

3.1 Tungsten accumulation at high toroidal rotation velocity

Since 2014,EAST has been used to access advanced regimes with steady-state high-performance plasma in ITER-like tungsten divertor operation [38—41].However,the accumulation of tungsten that is transported from the divertor by regular collisions or large ELMy impacts[42,43] always degrades the performance of H-mode plasma.In this section,the behavior of tungsten under a combination of NBI and lower hybrid wave(LHW)heating in the type-I ELMy H-mode regime is investigated to explore suitable heating scenarios without high-Zimpurity accumulation and to develop the steady-state long-pulse and high-performance regime.

Type-I ELMy H-mode discharges with the ITER-like tungsten divertor are selected to investigate the tungsten behavior in the core plasma.In the discharge of shot#63723,type-I ELMs withfELM～80 Hz were obtained when additional co-NBI was performed during the LHW heating phase,as shown in Fig.2.This discharge was performed in the upper single-null (USN) configuration with an unfavorableBtdirection (counter-Ipdirection,where the direction ofBt×?Btpoints away from the divertor X-point) atBt=2.3T,Ip=0.5 MA,a line-averaged electron density,and a safety factor at 95% of the magnetic flux,q95,of 5,where the source power of LHW and co-NBI heating is 1.5 and 4.0 MW,respectively.In this shot,the toroidal rotation velocity of the background plasma can be driven to a quite high level [a central toroidal rotation velocityVt(0) of～140 km/s] because co-NBI heating is applied.

Fig.2 (Color online)Waveforms of a typical tungsten accumulation shot,#63723.a Plasma current Ip and line-averaged electron density,b total LHW and NBI heating power,c intensity of divertor Dα emission,and plasma stored energy WMHD,d central toroidal rotation velocity Vt(0) and energy confinement factor H98,and e external Mirnov coil signal and normalized beta βN

The emission from high ionization stages of Li,C,Fe,Cu,Mo,and W is detected by EUV spectrometers at wavelengths of 20—150 A° [35,36].Figure 3a—c shows the EUV spectra and total radiative power (Prad) at 3.1,3.3,and 4.1 s,respectively.The stronger emission lines,such as Li III (134.99 A°),C VI (33.73 A°),and Fe XXIII(132.91 A°),and the first-and second-order tungsten emission in the unresolved transition array (W-UTA) are labeled.CWin core plasma is evaluated as the W-UTA intensity at 45—70 A°,which is composed of W24+—W45+.TheCWcalculation method is analyzed and introduced in literature [44].The W-UTA intensity is found to increase sharply with increase inPrad.In particular,in Fig.3b and c,whenPradis changed from 1.0 to 1.3 MW,the intensity of the first-order W-UTA becomes saturated,and the secondorder W-UTA emission appears and becomes more intense,whereas the intensity of emission from other impurities remains relatively low.This result confirmed that tungsten accumulates easily in the core plasma when co-NBI is performed and makes a large contribution to this important radiation.

Fig.3 (Color online) EUV spectra from shot #63723 in wavelength range of 20—150 A° at times of (a) 3.1 s,(b) 3.3 s,and (c) 4.1 s

The process of tungsten accumulation is shown in Fig.4a and b.CWincreases continually during each tungsten accumulation phase,and the radiation profile becomes progressively peaked correspondingly.Noticeably,in the first phase of the H-mode,becauseq(0)is slightly less than 1[Fig.5d],core MHD activity(m/n=1/1)begins and is observed on the spectrogram of the core soft X-ray emission signals from 3.0 to 3.6 s,as shown in Fig.4c.This core MHD activity may enhance the impurity diffusivity and decrease the neoclassical convection[45,46],which suppresses further increases inCWin the core plasma.Although this MHD activity vanishes after 3.6 s,CWcontinues to grow.In this case,tungsten accumulation results in a larger radiated power loss;thus,a transition back to the L-mode phase occurs after approximately one energy confinement time τE[τE=WMHD/(Pto-tal-dWMHD/dt)].These sequences of tungsten accumulation and H—L transitions can repeat several times with a period of approximately 3—6τE.In each tungsten accumulation phase,the effective charge numberZeffincreases from～3 to～3.5,andCWincreases from～0.23 × 10-4to～1.6 ×10-4;thus,the power radiation fractionPrad/Ptotalincreases from～15% to a peak value of～50%.The energy confinement factorsH98and βNare clearly attenuated during each tungsten accumulation process,as shown in Fig.2.

Fig.4 (Color online) Analysis of tungsten accumulation from 3.0 to 4.8 s in shot#63723.a Tungsten concentration CW,b radiation profile measured by 64-channel fast bolometer system and W-UTA emission intensity,and c spectrogram of central chord of soft X-ray emission signals

Figure 5 a—d shows the electron density,electron temperature,ion temperature,and safety factor profiles,respectively,during the L-mode phase (at～3.1 s,blue line) and before (at～3.3 s,red line) and during (at～4.0 s,black line) tungsten accumulation in the H-mode phase of shot #63723.The strong radiative power loss caused by tungsten accumulation changes the profiles and results in an aggravated plasma performance and steadystate operation,as follows.(i) The electron temperature profile inside ρ=0.2 becomes hollow owing to radiative power loss in the tungsten accumulation phase,as shown in Fig.5b.This change in the electron temperature profile in the core region can reduce the turbulent outward pinch driven by the TEM [47],which decreases the outward impurity transport.(ii)The increase in electron density andZeffduring the tungsten accumulation phase can reduce the ITG mode growth rate and particle diffusion [48,49],producing a peak in ion temperature in the central plasma region,as shown in Fig.5c.(iii)The reason for the increase in the safety factor in the core region[Fig.5d]may be that plasma cooling and the increase inZeffin the accumulation phase alternate the plasma electric conductivity and the current redistribution.Overall,Ti/Teincreases as co-NBI is added,indicating that the turbulent transport is reduced and tungsten accumulation is accelerated.However,there have been no reports that the turbulent transport exhibitsZdependence [50] (that the impurity accumulation increases with charge numberZ[51]);thus,another mechanism needs to be considered to explain the observation in Fig.3.

Fig.5 (Color online) Radial profiles in shot #63723:a electron density ne,b electron temperature Te,c ion temperature Ti,and d safety factor q.Blue,red,and black lines indicate times in the L-mode phase(3.1 s)and before(3.3 s) and during (4.1 s) the tungsten accumulation phase.ρ is the normalized poloidal flux coordinate,and the dashed black line in d indicates q=1

The contribution of neoclassical transport to tungsten accumulation can be indicated by the ratio of the neoclassical radial drift velocityvneoand the diffusion coefficientDneo(vneo/Dneo).In this analysis,because the core collisionality is greater than 1,vneo/Dneocan be obtained using the simplified expression [10]

where the subscript CL denotes the classical regime,and PS denotes the Pfirsch—Schlu¨ter regime.According to neoclassical theory,in each regime [10],v/Dcan be predicted using the expression

wherenDandTDare the main ion density and temperature,respectively;ZIis the impurity charge number;andZDis the main ion charge number.In the CL regime,H=0.5,and in the PS regime,His determined by the impurity concentration and collisionality of the main ions.The values ofvneo/Dneobefore and during co-NBI in shot#63723 are estimated assumingnD=neandTD=Ti,and are shown in Fig.6a.The value of |vneo/Dneo| is clearly smaller during the co-NBI phase than before co-NBI,which does not explain the tungsten accumulation observed in this shot.

Fig.6 (Color online) Estimates of a vneo/Dneo and b vPHZ+vEr before co-NBI [Vt(0)～20 km/s] and during co-NBI[Vt(0)～140 km/s] in shot#63723.The shaded areas in a are based on the H values,and the dashed lines indicate the averaged vneo/Dneo

The PHZ must be taken into account for plasma with a high toroidal rotation velocity [12].The PHZ velocity is given by the expression

wherevd0is the magnetic drift velocity,vd0=/Z0eRB;Z0is the impurity ion charge state in ionization equilibrium;ω is the angular frequency of the poloidal motion of impurity ions,ω=（Vt-Er/Bθ）/（qR0）;CT=ne?（γk-αk）/?Te;CZ=ne?（γk-αk）/?Z;andC?T=dT/dρ.In addition,mis the impurity mass,Bθis the poloidal magnetic field,γkis the ionization rate,and αkis the recombination rate(γkand αkare obtained from[52]).Eris evaluated using the ion radial force balance equation of the background plasma[53]:

wherePi=ni×Ti(niis the ion density),andVpis the poloidal rotation velocity.BecauseEris a function ofVt,theErpinch due to the effect ofErresulting from Coulomb collisions [12] is also considered in this analysis.TheErpinch is given by the equation.

To simplify the estimation ofvPHZ+vErbefore and during co-NBI,Vt(0)=20 km/s before co-NBI andVt(r)=Vt(0)e-r/athroughout the plasma radial region are assumed.In the ion radial force balance equation,the second term is not included becauseVpis small in the core region,and the third term can be simplified to ?（ne×Ti）/enefor bulk plasma.As the toroidal rotation velocity increases from 20 to 140 km/s,the estimatedvPHZ+vEr(ρ～0.5)is found to increase from approximately— 0.02 m/s to approximately -1.0 m/s,as shown in Fig.6b,which results in increased inward transport of tungsten and causes tungsten accumulation.Because the PHZ velocity depends strongly on the impurity mass,the increase in inward pinch velocityvPHZ(from approximately -0.017 m/s to approximately -0.8 m/s around ρ=0.5)can explain why tungsten accumulates more easily than other impurities(Fig.3).

3.2 Tungsten control using counter-NBI and favorable Bt

To suppress or avoid tungsten accumulation in NBIheated ELMy H-mode plasmas,the values ofTi/Te,vneo/Dneo,andvPHZ+vErshould be decreased.For this purpose,counter-NBI is applied to modify the density and temperature and brake the plasma toroidal rotation.Normally,the counter-NBI is not good for plasma confinement other than the co-NBI because counter-NBI introduces impurities into the main plasma and degrades the confinement[54].To rule out the possibility that the decreased tungsten content is due to the loss of confinement,another mechanism is needed for reducing the penetration of impurities into the main plasma and compensating for plasma confinement.It has been found that a favorableBt(in the co-Ipdirection,where the direction ofBt×?Btpoints toward the X-point) can improve plasma confinement[55](compared with an unfavorableBt,a favorableBtcan improve theH98by～20% in shots with similar plasma parameters).In addition,favorableBtis found to have a positive effect on tungsten control,which will be discussed in the next section.For these reasons,the tungsten behavior was investigated using counter-NBI with a favorableBtdirection,as in shot #71605.As shown in Fig.7,type-I ELMs withfELM～40 Hz were obtained when additional co-NBI was performed during the LHW heating phase.This shot was the last shot before lithium wall conditioning and was performed in the USN configuration withBt=2.4 T,Ip=0.6 MA,～4.2 × 1019m-3,q95=4.8,and maximum source powers of LHW and NBI of 1.6 and 4.8 MW,respectively.

Fig.7 (Color online) Waveforms of EAST shot #71605:a plasma current Ip and line-averaged electron density ,b source power of LHW and NBI (co-and counter-injection directions),c Dα emission signal from upper divertor and plasma stored energy WMHD,d central toroidal rotation velocity Vt(0) and XUV emission signals at r/a=0,e energy confinement factor H98 and normalized beta βN,f tungsten influx at inner target plate,where the black line denotes its tendency.The tungsten flux is indicated by the intensity of the WI emission line at 400.9 nm from the filterscope diagnostic [42]

The central radiation (as indicated by the central XUV emission intensity)clearly decreases rapidly when counter-NBI is performed,although counter-NBI slightly enhances the tungsten source (as indicated by the intensity of W I line emission at 400.9 nm) near the upper divertor region.Further analysis shows that the evolution ofCWis directly correlated with the control ofVt(0)by co-and counter-NBI.As shown in Fig.7d,when co-NBI heating with a 2 MW source power is applied at 3.0 s,Vt(0) begins to increase and reaches a peak value of approximately 150 km/s.As shown in Fig.8a and b,the increase in plasma toroidal rotation velocity results in a rapid increase inCWand a gradual peak in the W-UTA emission intensity measured by EUV_Long2.EUV_Long2 [56] is a newly developed space-resolved EUV spectrometer that can provide the radial profile of impurity line emission covering a wide vertical range of -0.085 m ≤Z≤0.4 m [Fig.12b].These rapid variations inCWand the W-UTA profile in the core region suggest tungsten accumulation.However,CWdecreases and the W-UTA emission intensity in the core plasma suddenly dissipates when counter-NBI heating with a 1.4 MW source power is superimposed at 4.0 s to reduceVt(0) to 45 km/s.When the counter-NBI source power is increased to 2.8 MW,Vt(0) can be further reduced and finally maintained at～20 km/s.When the counter-NBI source power decreases and vanishes,CWincreases again with the increase inVt(0).In addition,additional peaks in the W-UTA emission intensity in the core region indicate that tungsten accumulation again occurs.In this process,when counter-NBI with a maximum power of 2.8 MW is added during the co-NBI phase,Vt(0)changes from 150 to 20 km/s.Consequently,CWdecreases from～7 × 10-5to～2 × 10-5;thus,the power radiation fractionPrad/Ptotalis decreased from～32% to～8%.

Fig.8 (Color online)Time evolution of a tungsten concentration CW,b contour plot of W-UTA emission intensity measured by a newly developed space-resolved EUV spectrometer(EUV_Long2 in Fig.1),c frequency of ELMs fELM,and d normalized amplitude of ELMs AELM (black lines indicate their tendency) in shot #71605

In the ELMy H-mode phase,the superimposed counter-NBI also changes the characteristics of ELMs[57](such as frequencyfELMand amplitudeAELM),as shown in Fig.8(c) and (d).It has been found that high-frequency ELMs can provide continuous transport of impurities in the edge plasma region,which prevents the accumulation of impurities [58—60].This effect of ELMs on the elimination of impurities is called impurity flushing.The correlation analysis ofVt(0),fELM,andCWfrom 4 to 4.55 s shows thatCWis strongly dependent onVt(0)but weakly dependent onfELM.This weak dependence ofCWonfELMresults from the fact that the low-frequency ELMs (which increase from～50 to～65 Hz during that time period) are insufficiently effective in flushing out impurities[58].Note also that no core MHD activity is observed on the spectrogram of the central chord of the soft X-ray emission signals,and the energy confinement factorH98remains nearly constant during the counter-NBI heating phase.This result implies that the counter-NBI added here does not degrade the plasma confinement and can effectively suppress tungsten accumulation in the type-I ELMy H-mode on EAST.

However,the decrease inCWduring the added counter-NBI phase is attributed not only to the reduced toroidal rotation velocity,but also to the variations in density and temperature and their gradients.Figure 9a—c shows the density and temperature profiles at time points before (at 3.5 s,red line),during (at～5.0 s,blue line),and after(at～7.0 s,black line) counter-NBI in shot #71605.Under the assumption described above,the estimatedvPHZ+vEr(ρ～ 0.4) decreases from approximately -0.57 m/s to approximately -0.02 m/s when the counter-NBI changesVt(0) from～80 to～20 km/s,as shown in Fig.10(b).The figure also shows that the applied counter-NBI increases the density and flattens its profile in the region 0.3 ＜ρ ＜0.9 and clearly enhances the ion temperature and its gradient inside ρ～0.7.These changes will increase the impurity diffusivity owing to the enhanced collisionality [61] and may decrease the pinch velocity of impurities in the central region.

Fig.9 (Color online)Radial profiles of shot #71605:a electron density ne,b electron temperature Te,and c ion temperature Ti.Red,blue,and black lines indicate times before (3.5 s),during (5.0 s),and after (7.0 s) counter-NBI,respectively

Fig.10 (Color online)Estimated a vneo/Dneo and b vPHZ+vEr during co-NBI[Vt(0)～80 km/s] and added counter-NBI [Vt(0)～20 km/s]in shot #71605.The shaded areas in a are based on the H values,and the dashed lines indicate the average vneo/Dneo

The values ofvneo/Dneoin the co-NBI and counter-NBI phases in shot #71605 were also estimated and are shown in Fig.10a.When counter-NBI is added during the co-NBI phase,vneo/Dneooutside of ρ～0.2 can be reduced,and the direction ofvneo/Dneocan even be reversed from inward to outward around ρ=0.6.The reversed direction is determined by the ratio of the normalized temperature gradient to the normalized main ion density gradient,.It is inferred that the added counter-NBI in this shot reduces the PHZ +Erinward pinch velocity and weakens the contribution of neoclassical inward transport.Note,however,thatTi/Teincreases slightly when counter-NBI is added,which should reduce the turbulent impurity transport and increase tungsten accumulation.Thus,it can be concluded that the decrease inCWduring the counter-NBI phase is not caused by turbulent transport.

3.3 Effect of favorable Bt on tungsten control

In this section,the effect ofBtdirection on tungsten control is discussed.It was found that the impurity confinement time(τimp)was significantly reduced in favorableBtexperiments in recent EAST campaigns.Typically,τimpcan be reduced by nearly half for favorableBtH-mode plasmas (e.g.,τimp～112.0 ms for unfavorableBtin shot#79062 and τimp～61.2 ms for favorableBtin shot#80303).Shots #79062 and #80303 are similar shots that have oppositeBtdirections [55].The effect was further studied by investigatingCWand the tungsten flux in the divertor.To decrease the effect of toroidal rotation velocity on tungsten transport,the shots with lowVt(0)(＜30 km/s)are chosen for comparison.CWversusVt(0)in steady-state ELMy H-mode discharges with unfavorable and favorableBtdirections is compared in Fig.11.The selected shots have the following parameters:Ip=0.45—0.55 MA,=2.8—4.5 × 1019m-3,and the total heating powers are similar.This figure shows thatCWcan be easily controlled below 8 × 105in the ELMy H-mode plasmas with reduced toroidal rotation velocity [Vt(0)＜30 km/s],whereCWis generally lower in the shots with favorableBtdirection than in those with unfavorableBtdirection.These twoBtdirections produceEr×Bdrift with opposite poloidal directions,as shown in Fig.11c.For the unfavorableBtdirection,the divertorEr×Bdrift increases the backflow(PS flow) of impurity particles from the upper outer divertor to the upstream region.For the favorableBtdirection,the backflow of impurity particles is directed to the lower divertor and can be exhausted by the closed divertor structure [62].The lower divertor structure is closed because it is W-shaped with a dome.This effect partially explains the lowerCWin the shots with the favorableBtdirection.In addition,theEr×Bdrift in the favorableBtdirection increases the tungsten flux on the inboard side[as indicated by the increasedIIn/IOutshown in Fig.11b],which may affect tungsten transport in the core region owing to the poloidal asymmetry of the tungsten density when the toroidal rotation velocity is low.The in—out asymmetry effect is beneficial for removing the tungsten from the core plasma [63—65],for example,by reversing the direction ofvneo/Dneofrom inward to outward.

Fig.11 (Color online) a Tungsten concentration versus toroidal rotation velocity in ELMy H-mode plasmas with parameters of Ip=0.45—0.6 MA and =2.8—4.5 × 1019 m-3,b ratio of tungsten flux density near the inboard and outboard sides(IIn/IOut)in shots withunfavorable(red)and favorable Bt(black)directions,and c schematic diagram of drifts and parallel flows in the corresponding USN configuration

The W-UTA emission intensity profiles in ELMy H-mode plasmas with differentBt×?Btdrift directions are presented in Fig.12a.To eliminate the effect of electron temperature on tungsten transport,two shots with similar magnetic axes (nearZ=0.024 m) and central electron temperatures are compared.Figure 12a shows that the W-UTA emission intensity profile is flatter for the favorableBtdirection than for the unfavorableBtdirection.The flatter W-UTA emission intensity profile suggests that the favorableBtdirection is beneficial for avoiding tungsten accumulation and can decrease the risk of intense radiated power in the core region.From these facts,it can be deduced that lowerCWwith favorableBtwill expand the operating window toward steady-state type-I ELMy H-mode operation.

Fig.12 (Color online)a Comparison of W-UTA emission intensity profiles in type-I ELMy H-mode plasmas with unfavorable and favorable Bt directions,b viewing range of space-resolved EUV spectrometer (blue) and plasma shape in the shots with unfavorable (red) and favorable(black) Bt directions in EAST

4 Summary and discussion

The tungsten behavior in ELMy H-mode plasmas with co-/counter-NBI and unfavorable/favorableBton EAST was investigated.The results showed that tungsten always accumulates in ELMy H-mode plasma with co-NBI owing to the increased inward pinch velocityvPHZ+vErand decreased turbulent transport.CWcan exceed 10-4in the type-I ELMy H-mode withVt(0)＞120 km/s and unfavorableBt,ultimately resulting in the degradation of plasma confinement and periodic H—L transitions.Thus,CWmust be maintained below 10-4for future steady-state type-I ELMy H-mode operation on EAST.

To suppress and avoid tungsten accumulation in ELMy H-mode plasmas,both the neoclassical inward convection and PHZ inward pinch velocities should be decreased,whereas turbulent transport needs to be increased.For this purpose,counter-NBI was applied to modify the density and temperature and brake the plasma toroidal rotation.The result shows thatCWcan be decreased from～7 × 10-5to～2 × 10-5in the type-I ELMy H-mode plasma with favorableBtwhen counter-NBI is superimposed during the co-NBI heating phase.Theanalysis indicated that the added counter-NBI not only decreases the PHZ +Erinward pinch velocity but also reverses thevneo/Dneodirection from inward to outward;both of these changes contribute to the suppression of tungsten accumulation.It was preliminarily found that,because NBI directly heats ions,turbulent transport in the core region driven by co-/counter-NBI is disadvantageous for the suppression of tungsten accumulation.Increasing the turbulent transport during the co-/counter-NBI phases by increasing the power of LHW heating or adopting onaxis electron cyclotron resonance heating is required.Further detailed studies remain as future work.

It was also found thatCWis lower in ELMy H-mode plasmas with favorableBtthan in similar discharges with unfavorableBt.Possible reasons for the effect ofBtdirection onCWwere discussed.On the one hand,in discharges with favorableBtin the USN configuration,the backflow of impurity particles is directed to the lower divertor and can be exhausted by the closed divertor structure.On the other hand,in low-rotation plasmas,theEr×Bdrift in plasmas with favorableBtincreases the tungsten flux on the inboard side,which can weaken the contribution of neoclassical inward convection because of in—out asymmetry.

These results demonstrate the potential of counter-NBI in the favorableBtdirection for suppressing tungsten accumulation in ELMy H-mode plasma.The application of counter-NBI to suppress tungsten accumulation and avoid plasma confinement degradation must ensure thatVtandIphave the same direction(Vt＞0),thatVtis sufficiently low,and thatBtis favorable.Otherwise,impurities can accumulate more easily during counter-NBI than during co-NBI[12,54],which destroys the plasma confinement.The lowerCWin the core plasma caused by counter-NBI in the favorableBtdirection will extend the operating window for exploring steady-state type-I ELMy H-mode operation in EAST in the near future.

Author contributionsAll authors contributed to the study conception and design.Material preparation,data collection,and analysis were performed by Zong Xu,Zhen-Wei Wu,Ling Zhang,Yun-Xing Cheng,Wei Gao,Xiao-Dong Lin,Xiang Gao,Ying-Jie Chen,Yue-Heng Huang,Juan Huang,Qing Zang,Yin-Xian Jie,Yan-Min Duan,Ying-Ying Li,Hai-Qing Liu,and Guo-Sheng Xu.The first draft of the manuscript was written by Zong Xu,and all authors commented on previous versions of the manuscript.All authors read and approved the final manuscript.