Numerical study of beam propagation and plasma properties in the neutralizer and the E-RID of the ITER Neutral Beam Injector

Non-ohmic heating will be used in the experimental nuclear fusion reactor ITER to reach thermonuclear temperatures. Two heating mechanism will be implemented, i.e. microwaves resonant with ion and electron cyclotron frequencies and energetic neutral beam injection, which contributes also to the current drive. Each one of the two neutral beam injector planned for ITER will deliver 16MW of 1MeV D-0 beam. In the injector, negative ions D-coming from a 40A negative ion source are electrostatically accelerated to 1 MeV, and stripped of their extra electron by collision with a target gas in a structure known as the neutralizer. Residual charged particles are deflected after the neutralizer in an electrostatic ion dump (E-RID). The ionization of the deuterium buffer gas filling the neutralizer induced by the D-beam creates a rarefied plasma which is expected to efficiently screens the Coulomb repulsion of the beam. Moreover, this plasma can eventually escape from the neutralizer and move back in the accelerator, towards the accelerating grids and the negative ion source. The transport of the beam through the neutralizer and the RID and the related plasma properties were studied using a 3D electrostatic particle-in-cell code called OBI-3 (Orsay Beam Injector 3 dimensional). Particle-particle and particle-wall collisions are treated using the Monte Carlo collision approach. Simulations show that the secondary plasma effectively screens the beam space charge preventing beam transverse expansion. Plasma ions created in the neutralizer form an upstream current with a magnitude of similar to 0.5% of the negative ion current. Gas breakdown leading to arc formation in the RID was not observed. Finally, results for the propagation of non-ideal beams coming from simulations of the extraction and consecutive acceleration taken from Revel et al 2013 Nucl. Fusion 53 073027 are presented.