Design of the “half-size” ITER neutral beam source for the test facility ELISE

Long-pulse diagnostic calorimeter for the negative ion source testbed BATMAN upgrade

The RF-driven negative ion source testbed BATMAN upgrade is being developed at IPP Garching in the framework of the ion source development for ITER and DEMO neutral beam injection systems. The testbed has recently been enhanced to allow for steady state operation with a focus on beam optics studies. The previous titanium sublimation pumps and inertial calorimeter limited the beam pulse length to about 6 s every 3 min. The upgrade comprises a long-pulse compatible, actively cooled diagnostic calorimeter. This has been designed and is currently being manufactured to substitute the inertially cooled calorimeter that has limited diagnostic capabilities. The new diagnostic calorimeter consists of a copper plate with dimensions of 910 × 660 × 25 mm³ placed about 2 m from the ion source extraction grids, and through a novel solution, it will provide a 2D profile of beam power density with a 20 mm spatial resolution. Water flowing through cooling channels embedded in the copper plate will actively cool the calorimeter, which is loaded with about 160 kW beam power at ITER-relevant current density, but 45 kV acceleration. A fraction of the beam will pass through many small apertures (ø2 mm) positioned in the calorimeter plate and will be collected by thin (0.2 mm) copper foils attached to the calorimeter back side. Evaluation of power density will be performed by measuring the temperature of the heat flux foils with a high-resolution infrared camera observing the calorimeter from the back side and calibrated by thermocouples.

Simulation of the gas density distribution in the accelerator of the ELISE test facility

In multiaperture electrostatic accelerators of negative ion sources, the plasma discharge is sustained by injecting gas in the plasma source, in a dynamic equilibrium with the gas flowing out through the accelerator. In this work, we present a three-dimensional numerical simulation of the gas flow inside the accelerator system of the large negative ion source ELISE at Max-Planck-Institut für Plasmaphysik Garching. ELISE has 640 apertures per electrode and lateral gaps between the electrode support structures that also contribute to the total gas conductance. Assuming molecular regime, we estimated the gas conductance, the gas density profile along the path of the ion beams from upstream of the plasma grid to downstream of the ground grid, and the transverse nonuniformities in the accelerator. The simulation included the most relevant geometrical features, while the results are compared to analytical results.

Formation of large negative deuterium ion beams at ELISE

Negative ion sources for neutral beam injection (NBI) in fusion experiments are based on the surface production of H^- or D^- on cesiated low work function surfaces. In the recent years, it was demonstrated at the large RF driven ion source of the ELISE (Extraction from a Large Ion Source Experiment) test facility that the requirements for the ITER NBI systems can be fulfilled by hydrogen. This is a big step toward the first operational period of ITER, planned for up to 2035. However, for the following operational period, neutral beam systems working in deuterium are needed. Operation of negative hydrogen ion sources in deuterium is significantly more demanding than in hydrogen: the amount of coextracted electrons is much higher and their increase during pulses is much more pronounced, limiting the achievable performance. This paper presents the results of investigations aimed to improve the insight into the physics related to this isotope effect. Due to the higher atomic mass of deuterium, cesium is removed much more effectively from reservoirs at the walls, resulting in a depletion of these reservoirs and a strongly increased cesium density in the plasma. Additionally, a correlation between the fluxes of charged particles toward the inner ion source surfaces and the coextracted electrons is identified.

Structure design and analysis of RF ion source for negative ion source test facility

A negative ion source acts as a critical part in a neutral beam injector (NBI). A high current ion source is required for the high-power NBI. In this paper, a prototype radio frequency (RF) ion source and its test facility are developed in the Institute of Plasma Physics, Chinese Academy of Sciences, to demonstrate the key technology of the high power negative ion source. The structure design of the RF negative ion source is presented, involving the designs of the ion source plasma generator and accelerator. The detailed structure design and analysis of the key parts of the ion source are also presented, such as the Faraday shield (FS) and accelerator grids. The fluid-thermal-structural coupling characteristics of the FS and grid are explored with different mechanisms of fluid pressure, RF power, and the structure type on the thermal stress. Then, the processing and manufacturing scheme of the FS and grids are also given. Finally, the results were presented with a manufactured three cooling channel FS. The experimental results prove that the developed structure design of the RF ion source is effective and reliable, and the correctness of finite element analysis is also verified by experimental data comparison.

Review of particle-in-cell modeling for the extraction region of large negative hydrogen ion sources for fusion

Particle-in-cell (PIC) codes are used since the early 1960s for calculating self-consistently the motion of charged particles in plasmas, taking into account external electric and magnetic fields as well as the fields created by the particles itself. Due to the used very small time steps (in the order of the inverse plasma frequency) and mesh size, the computational requirements can be very high and they drastically increase with increasing plasma density and size of the calculation domain. Thus, usually small computational domains and/or reduced dimensionality are used. In the last years, the available central processing unit (CPU) power strongly increased. Together with a massive parallelization of the codes, it is now possible to describe in 3D the extraction of charged particles from a plasma, using calculation domains with an edge length of several centimeters, consisting of one extraction aperture, the plasma in direct vicinity of the aperture, and a part of the extraction system. Large negative hydrogen or deuterium ion sources are essential parts of the neutral beam injection (NBI) system in future fusion devices like the international fusion experiment ITER and the demonstration reactor (DEMO). For ITER NBI RF driven sources with a source area of 0.9 × 1.9 m² and 1280 extraction apertures will be used. The extraction of negative ions is accompanied by the co-extraction of electrons which are deflected onto an electron dump. Typically, the maximum negative extracted ion current is limited by the amount and the temporal instability of the co-extracted electrons, especially for operation in deuterium. Different PIC codes are available for the extraction region of large driven negative ion sources for fusion. Additionally, some effort is ongoing in developing codes that describe in a simplified manner (coarser mesh or reduced dimensionality) the plasma of the whole ion source. The presentation first gives a brief overview of the current status of the ion source development for ITER NBI and of the PIC method. Different PIC codes for the extraction region are introduced as well as the coupling to codes describing the whole source (PIC codes or fluid codes). Presented and discussed are different physical and numerical aspects of applying PIC codes to negative hydrogen ion sources for fusion as well as selected code results. The main focus of future calculations will be the meniscus formation and identifying measures for reducing the co-extracted electrons, in particular for deuterium operation. The recent results of the 3D PIC code ONIX (calculation domain: one extraction aperture and its vicinity) for the ITER prototype source (1/8 size of the ITER NBI source) are presented.

MS-CASPT2 study of the ground and low lying states of CsH. J Mol Model 2017;23:339. [PMID: 29124408 DOI: 10.1007/s00894-017-3503-9]

Correlated ab initio methods [CASPT2 and R-CCSD(T)] in conjunction with the ANO-RCC basis sets in large contraction were used to calculate potential energy curves (PECs) of the ground and excited electronic states of CsH⁺ (doublets and quartets) with the inclusion of the scalar relativistic effects and spin-orbit interaction. The ground X²Σ⁺ state is a rather fragile van der Waals molecular ion. The binding energy of this X²Σ⁺ state provided by both computational methods is estimated to be 0.02-0.04 eV, and is compared with the reported experimental binding energy (0.51-0.77 eV). This large binding energy can be attributed to the A²Σ⁺ state, and can thus explain the apparent disagreement between theory and experiment. The spectroscopic constants of all bound states were calculated from the PECs and compared with previous published data for X²Σ⁺ and A²Σ⁺ states. Graphical abstract Low-lying Ω states of cesium hydride cation.

Caesium hydride: MS-CASPT2 potential energy curves and A1Σ+→X1Σ+ absorption/emission spectroscopy

Correlated ab initio methods (CASPT2 and CCSD(T)) in conjunction with the ANO-RCC basis sets were used to calculate potential energy curves (PECs) of the ground, valence, and Rydberg electronic states of CsH with the inclusion of the scalar relativistic effects. The spectroscopic constants of bound states were calculated from the PECs and compared with previous theoretical and/or available experimental data. Absorption and emission spectra arising from the transition between X¹Σ⁺ and A¹Σ⁺ states were modelled using vibrational and rotational energy levels and corresponding nuclear wave functions obtained via the direct numerical integration of one-dimensional rovibrational Schrödinger equation in the CASPT2/ANO-RCC electronic potentials. The anharmonic shape of the A¹Σ⁺ potential and the shape of the pertinent vibrational wave functions have an interesting impact on the final shape of the spectrum and result in the complicated fine structure of individual emission bands.

Final design of the beam source for the MITICA injector

The megavolt ITER injector and concept advancement experiment is the prototype and the test bed of the ITER heating and current drive neutral beam injectors, currently in the final design phase, in view of the installation in Padova Research on Injector Megavolt Accelerated facility in Padova, Italy. The beam source is the key component of the system, as its goal is the generation of the 1 MeV accelerated beam of deuterium or hydrogen negative ions. This paper presents the highlights of the latest developments for the finalization of the MITICA beam source design, together with a description of the most recent analyses and R&D activities carried out in support of the design.

Towards 20 A negative hydrogen ion beams for up to 1 h: Achievements of the ELISE test facility (invited)

The large-scale RF-driven ion source of the test facility extraction from a large ion source experiment is aimed to deliver an accelerated ion current of 20 A D(-) (23 A H(-)) with an extracted electron-to-ion ratio below one for up to 1 h. Since the first plasma pulses for 20 s in volume operation in early 2013, followed by caesiation of the ion source, substantial progress has been achieved in extending the pulse length and the RF power. The record pulses in hydrogen are stable 400 s pulses with an extracted ion current of 18.3 A at 180 kW total RF power and 9.3 A at 80 kW stable for 1 h. For deuterium pulse, length and RF power are limited by the amount of co-extracted electrons.

3D numerical simulations of negative hydrogen ion extraction using realistic plasma parameters, geometry of the extraction aperture and full 3D magnetic field map

Decreasing the co-extracted electron current while simultaneously keeping negative ion (NI) current sufficiently high is a crucial issue on the development plasma source system for ITER Neutral Beam Injector. To support finding the best extraction conditions the 3D Particle-in-Cell Monte Carlo Collision electrostatic code ONIX (Orsay Negative Ion eXtraction) has been developed. Close collaboration with experiments and other numerical models allows performing realistic simulations with relevant input parameters: plasma properties, geometry of the extraction aperture, full 3D magnetic field map, etc. For the first time ONIX has been benchmarked with commercial positive ions tracing code KOBRA3D. A very good agreement in terms of the meniscus position and depth has been found. Simulation of NI extraction with different e/NI ratio in bulk plasma shows high relevance of the direct negative ion extraction from the surface produced NI in order to obtain extracted NI current as in the experimental results from BATMAN testbed.

First results of the ITER-relevant negative ion beam test facility ELISE (invited)

An important step in the European R&D roadmap towards the neutral beam heating systems of ITER is the new test facility ELISE (Extraction from a Large Ion Source Experiment) for large-scale extraction from a half-size ITER RF source. The test facility was constructed in the last years at Max-Planck-Institut für Plasmaphysik Garching and is now operational. ELISE is gaining early experience of the performance and operation of large RF-driven negative hydrogen ion sources with plasma illumination of a source area of 1 × 0.9 m(2) and an extraction area of 0.1 m(2) using 640 apertures. First results in volume operation, i.e., without caesium seeding, are presented.

Optical emission spectroscopy at the large RF driven negative ion test facility ELISE: instrumental setup and first results

One of the main topics to be investigated at the recently launched large (A(source) = 1.0 × 0.9 m(2)) ITER relevant RF driven negative ion test facility ELISE (Extraction from a Large Ion Source Experiment) is the connection between the homogeneity of the plasma parameters close to the extraction system and the homogeneity of the extracted negative hydrogen ion beam. While several diagnostics techniques are available for measuring the beam homogeneity, the plasma parameters are determined by optical emission spectroscopy (OES) solely. First OES measurements close to the extraction system show that without magnetic filter field the vertical profile of the plasma emission is more or less symmetric, with maxima of the emission representing the projection of the plasma generation volumes, and a distinct minimum in between. The profile changes with the strength of the magnetic filter field but under all circumstances the plasma emission in ELISE is much more homogeneous compared to the smaller IPP prototype sources. Planned after this successful demonstration of the ELISE OES system is to combine OES with tomography in order to determine locally resolved values for the plasma parameters.

The development of the radio frequency driven negative ion source for neutral beam injectors (invited)

Large and powerful negative hydrogen ion sources are required for the neutral beam injection (NBI) systems of future fusion devices. Simplicity and maintenance-free operation favors RF sources, which are developed intensively at the Max-Planck-Institut für Plasmaphysik (IPP) since many years. The negative hydrogen ions are generated by caesium-enhanced surface conversion of atoms and positive ions on the plasma grid surface. With a small scale prototype the required high ion current density and the low fraction of co-extracted electrons at low pressure as well as stable pulses up to 1 h could be demonstrated. The modular design allows extension to large source dimensions. This has led to the decision to choose RF sources for the NBI of the international fusion reactor, ITER. As an intermediate step towards the full size ITER source at IPP, the development will be continued with a half-size source on the new ELISE testbed. This will enable to gain experience for the first time with negative hydrogen ion beams from RF sources of these dimensions.

Physics design of a 100 keV acceleration grid system for the diagnostic neutral beam for international tokamak experimental reactor

This paper describes the physics design of a 100 keV, 60 A H(-) accelerator for the diagnostic neutral beam (DNB) for international tokamak experimental reactor (ITER). The accelerator is a three grid system comprising of 1280 apertures, grouped in 16 groups with 80 apertures per beam group. Several computer codes have been used to optimize the design which follows the same philosophy as the ITER Design Description Document (DDD) 5.3 and the 1 MeV heating and current drive beam line [R. Hemsworth, H. Decamps, J. Graceffa, B. Schunke, M. Tanaka, M. Dremel, A. Tanga, H. P. L. De Esch, F. Geli, J. Milnes, T. Inoue, D. Marcuzzi, P. Sonato, and P. Zaccaria, Nucl. Fusion 49, 045006 (2009)]. The aperture shapes, intergrid distances, and the extractor voltage have been optimized to minimize the beamlet divergence. To suppress the acceleration of coextracted electrons, permanent magnets have been incorporated in the extraction grid, downstream of the cooling water channels. The electron power loads on the extractor and the grounded grids have been calculated assuming 1 coextracted electron per ion. The beamlet divergence is calculated to be 4 mrad. At present the design for the filter field of the RF based ion sources for ITER is not fixed, therefore a few configurations of the same have been considered. Their effect on the transmission of the electrons and beams through the accelerator has been studied. The OPERA-3D code has been used to estimate the aperture offset steering constant of the grounded grid and the extraction grid, the space charge interaction between the beamlets and the kerb design required to compensate for this interaction. All beamlets in the DNB must be focused to a single point in the duct, 20.665 m from the grounded grid, and the required geometrical aimings and aperture offsets have been calculated.