Numerical study of alpha particle loss with toroidal field ripple based on CFETR steady-state scenario

Effects of plasma equilibrium parameters on the alpha particle loss with the toroidal field ripple based on the CFETR steady-state scenario have been numerically investigated by the orbit-following code GYCAVA. It is found that alpha particle losses decrease and loss regions become narrower with the plasma current increasing or with the magnetic field decreasing. It is because the ripple stochastic transport and the ripple well loss of alpha particle are reduced with the safety factor decreasing. Decrease of the plasma density and temperature can reduce alpha particle losses due to enhancement of the slowing-down effect. The direction of the toroidal magnetic field can significantly affect heat loads induced by lost alpha particle. The vertical asymmetry of heat loads induced by the clockwise and counter-clockwise toroidal magnetic fields are due to the fact that the ripple distribution is asymmetric about the mid-plane, which can be explained by the typical orbits of alpha particle. The maximal heat load of alpha particle for the clockwise toroidal magnetic field is much smaller than that for the counter-clockwise one.

Gyrokinetic simulations of the kinetic electron effects on the electrostatic instabilities on the ITER baseline scenario

The linear and nonlinear simulations are carried out using the gyrokinetic code NLT for the electrostatic instabilities in the core region of a deuterium plasma based on the International Thermonuclear Experimental Reactor(ITER)baseline scenario.The kinetic electron effects on the linear frequency and nonlinear transport are studied by adopting the adiabatic electron model and the fully drift-kinetic electron model in the NLT code,respectively.The linear simulations focus on the dependence of linear frequency on the plasma parameters,such as the ion and electron temperature gradientsκ_(Ti,e)≡R=L_(Ti,e),the density gradientκ_(n)≡R/L_(n)and the ion-electron temperature ratioτ=T_(e)=T_(i).Here,is the major radius,and T_(e)and T_(i)denote the electron and ion temperatures,respectively.L_(A)=-(δ_(r)lnA)^(-1)is the gradient scale length,with denoting the density,the ion and electron temperatures,respectively.In the kinetic electron model,the ion temperature gradient(ITG)instability and the trapped electron mode(TEM)dominate in the small and large k_(θ)region,respectively,wherek_(θ)is the poloidal wavenumber.The TEMdominant region becomes wider by increasing(decreasing)κ_(T_(e))(κ_(T_(i)))or by decreasingκ_(n).For the nominal parameters of the ITER baseline scenario,the maximum growth rate of dominant ITG instability in the kinetic electron model is about three times larger than that in the adiabatic electron model.The normalized linear frequency depends on the value ofτ,rather than the value of T_(e)or T_(i),in both the adiabatic and kinetic electron models.The nonlinear simulation results show that the ion heat diffusivity in the kinetic electron model is quite a lot larger than that in the adiabatic electron model,the radial structure is finer and the time oscillation is more rapid.In addition,the magnitude of the fluctuated potential at the saturated stage peaks in the ITGdominated region,and contributions from the TEM(dominating in the higher k_(θ)region)to the nonlinear transport can be neglected.In the adiabatic electron model,the zonal radial electric field is found to be mainly driven by the turbulent energy flux,and the contribution of turbulent poloidal Reynolds stress is quite small due to the toroidal shielding effect.However,in the kinetic electron model,the turbulent energy flux is not strong enough to drive the zonal radial electric field in the nonlinear saturated stage.The kinetic electron effects on the mechanism of the turbulence-driven zonal radial electric field should be further investigated.

SOLPS-ITER drift modeling of neon impurity seeded plasmas in EAST with favorable and unfavorable toroidal magnetic field direction

To better understand divertor detachment and asymmetry in the Experimental Advanced Superconducting Tokamak(EAST),drift modeling via the comprehensive edge plasma code SOLPS-ITER of neon impurity seeded plasmas in favorable/unfavorable toroidal magnetic field(BT)has been performed.Firstly,electrostatic potential/field(f/E)distribution has been analyzed,to make sure that f and E are correctly described and to better understand drift-driven processes.After that,drift effects on divertor detachment and asymmetry have been focused on.In accordance with the corresponding experimental observations,simulation results demonstrate that in favorable BTthe onset of detachment is highly asymmetric between the inner and outer divertors;and reversing BT can significantly decrease the magnitude of in-out asymmetry in the onset of detachment,physics reasons for which have been explored.It is found that,apart from the well-known E×B drift particle flow from one divertor to the other through the private flux region,scrape-off layer(SOL)heat flow,which is much more asymmetrically distributed between the high field side and low field side for favorable BTthan that for unfavorable B_T,is also a critical parameter affecting divertor detachment and asymmetry.During detachment,upstream pressure(P_u)reduction occurs and tends to be more dramatical in the colder side than that in the hotter side.The convective SOL heat flow,emerging due to in-out asymmetry in P_u reduction,is found to be critical for understanding divertor detachment and asymmetry observed in EAST.To better understand the calculated drastic power radiation in the core and upstream SOL,drift effects on divertor leakage/retention of neon in EAST with both BTdirections have been addressed for the first time,by analyzing profile of poloidal neon velocity and that of neon ionization source from atoms.This work can be a reference for future numeric simulations performed more closely related to experimental regimes.

Simulations of NBI fast ion loss in the presence of toroidal field ripple on EAST

NBI fast ion losses in the presence of the toroidal field ripple on EAST have been investigated by using the orbit code GYCAVA and the NBI code TGCO.The ripple effect was included in the upgraded version of the GYCAVA code.It is found that loss regions of NBI fast ions are mainly on the low field side near the edge in the presence of ripple.For co-current NBIs,the synergy effect of ripple and Coulomb collision on fast ion losses is dominant,and fast trapped ions located on the low field side are easily lost.The ripple well loss and the ripple stochastic loss of fast ions have been identified from the heat loads of co-current NBI fast ions.The ripple stochastic loss and the collisioninduced loss are much larger than the ripple well loss.Heat loads of lost fast ions are mainly localized on the right side of the radio frequency wave antennas from the inside view toward the first wall.For counter-current NBIs,the first orbit loss due to the magnetic drift is the dominant loss channel.In addition,fast ion loss fraction with ripple and collision for each NBI linearly increases with the effective charge number,which is related to the pitch angle scattering effect.

Simulations of first-orbit losses of neutral beam injection(NBI)fast ions on EAST

Simulations of first-orbit losses of neutral beam injection(NBI)fast ions in the EAST tokamak have been studied in detail by using the orbit-following code GYCAVA and the NBI code TGCO.Beam ion losses with the wall boundary are smaller than those with the last closed flux surface boundary.In contrast to heat loads on the wall without radio frequency wave(RFW)antennas,heat loads on the wall with RFW antennas are distributed more locally near the RFW antennas.The direction of the toroidal magnetic field dramatically affects the final positions of lost fast ions,which is related to the magnetic drift.The numerical results on heat loads of beam ions corresponding to different toroidal magnetic fields are qualitatively consistent with the experimental results.Beam ion losses increase with the beam energy for the co-current NBIs and the counter-perpendicular NBI.We have studied the behavior of fast ions produced by a small section neutral beam(beamlet)by using the numerical tool NBIT.The distributions of the loss fraction of beamlet fast ions peaked near the edge of the beam section for the counter-current NBIs,and they are related to the injection angle.This indicates that the first-orbit losses can be reduced by changing the shape of beam cross section.

Effects of fast ions produced by ICRF heating on the pressure at EAST

Ion cyclotron resonance heating(ICRH),which can produce fast ions,is an important auxiliary heating method at EAST.To analyze the effect of ICRH-induced fast ions on the plasma pressure at EAST,simulations are performed using TRANSP and TORIC codes.It is found that the ICRF-induced fast ion pressure cannot be negligible when the ICRF power is sufficiently high.The magnitude of the total ion pressure can be raised up to 60%of the total pressure as the input power rises above 3 MW.The pressure profile is also significantly modified when the resonant layer is changed.It is shown that by changing the wave frequency and antenna position,the total ion pressure profile can be broadened,which might provide an option for profile control at EAST.

Implementation of field-aligned coordinates in a semi-Lagrangian gyrokinetic code for tokamak turbulence simulation

Field-aligned coordinates have been implemented in the gyrokinetic semi-Lagrangian code NLT, Ye et al （2016 J. Comput. Phys. 316 180）, to improve the computational efficiency for the numerical simulations of tokamak turbulence and transport. 4D B-spline interpolation in field- aligned coordinates is applied to solve the gyrokinetic Vlasov equation. A fast iterative algorithm is proposed for efficiently solving the quasi-neutrality equation. A pseudo transform method is used for the numerical integration of the gyro-average operator for perturbations with a high toroidal mode number. The new method is shown to result in an improved code performance for reaching a given accuracy. Some numerical tests are presented to illustrate the new methods.

Simulations of energetic alpha particle loss in the presence of toroidal field ripple in the CFETR tokamak

Energetic alpha particle losses with the toroidal field ripple and the Coulomb collision in the CFETR tokamak have been simulated by using the orbit-following code GYCAVA for the steady-state and hybrid scenarios.The effects of the outer boundary and the ripple amplitude on alpha particle losses have been investigated.The loss fractions and heat loads of alpha particles in the hybrid scenario are much smaller than those in the steady-state scenario for a significant ripple amplitude.Some alpha particles in the plasma core are lost due to the ripple stochastic transport for a large ripple amplitude parameter.The heat loads with the last closed flux surface boundary are different from those with the wall boundary for the CFETR tokamak,which can be explained by typical alpha particle orbits.Discrete heat load spots have been observed in alpha particle loss simulations,which is due to the ripple well loss.The transition of the lost alpha particle behavior from the ripple stochastic diffusion to the ripple well trapping has been identified in our CFETR simulations.The Coulomb collision effect is responsible for this transition.

China's Export-led Growth Strategy:An International Comparison

A pressing challenge for China is determining where to accommodate millions of migrant workers displaced by the closing of many export-oriented factories. The current global financial crisis＂ has exposed the fragility of the export-led growth strategy China has adopted over the past 30years. Is there a better alternative for providing non-agricultural jobs than the sweatshops of cheap export production？ In the present paper, international experience is reviewed to shed light on China＇s situation. Using pooled regression models, we analyze data from the World Bank for 209 economies. We investigate the experience of other economies to answer the following questions： What is the common process of expanding the nonagricultural economy？ How is that process affected by the level of the real exchange rate？ Is export production a common way of absorbing surplus rural labor？ Finally, what are the ways that domestic demand and service employment can be expanded？