Neutral Beam Ion Research Articles

Effect of the cation-to-anion mass ratio of ionic liquid on plume neutralization characteristics for novel electrospray thruster

Abstract The ionic liquid electrospray thruster (ILET), a promising technology for space electric propulsion, has been gaining attention due to its theoretical capability to operate without an external neutralizer, which is a key component for traditional ion thrusters. Some experiments have shown that cation-to-anion mass ratio significantly affects the self-neutralization of thruster plumes, but further depth research is still needed. Therefore, numerical simulations of ion emission process both in uni-thruster and bi-thruster modes of ILETs are conducted using particle-in-cell method, to investigate the effect of the cation-to-anion mass ratio in ionic liquid on plume neutralization characteristics. In bi-thruster mode, the anions and cations are emitted simultaneously, while in uni-thruster mode, only anions or cations are emitted under the action of electric field force. Plume neutralization characteristics of four ionic liquid propellants, with different cation-to-anion mass ratio, in bi-thruster mode are obtained and compared to the uni-thruster mode. The results show that the plume profile morphology is significantly dependent on operating mode of the thruster and the cation-to-anion mass ratio of ionic liquid. Unlike the circular profile in uni-thruster mode, the plume in bi-thruster mode exhibit distortion with vertical stratification. The contours of the high-potential region in bi-thruster mode also distort, showing stratification and a V-shape, respectively, influenced by the mass ratio. Ion beam neutralization in bi-thruster mode is primarily achieved through the horizontal displacement of anion and cation beams, as well as their interactions. The cation-to-anion mass ratio influences the oscillation effects of the mass center in anion and cation clouds during neutralization, which in turn affects the neutralization effect. Ionic liquids with a significant mass ratio difference between cations and anions necessitate a longer time and greater distance for neutralization under identical conditions. These findings are instrumental for understanding the plume neutralization process and advancing the design of ILETs.

Comprehensive research facility for negative ion source neutral beam injection at CRAFT: design and first operations

Abstract Neutral beam injection (NBI) has been proven as a reliable heating and current drive method for fusion plasma. For the high-energy NBI system (particle energy > 150 keV) of large-scale fusion devices, the negative ion source neutral beam injection (NNBI) system is inevitable, which can obtain an acceptable neutralization efficiency (> 55%). But the NNBI system is very complex and challengeable. To explore and master the key NNBI technology for future fusion reactor in China, an NNBI test facility is under development in the framework of the Comprehensive Research Facility for Fusion Technology (CRAFT). The initial goal of CRAFT NNBI facility is to achieve a 2 MW hydrogen neutral beam at the energy of 200–400 keV for lasting 100 s. In the first operations of the CRAFT NNBI facility, a negative ion source with dual RF drivers was developed and tested. By using the 50 keV accelerator, the long-pulse and high-current extractions of negative hydrogen ions have been achieved and the typical values were 55.4 keV, 7.3 A (~ 123 A/m2), 105 s and 55.0 keV, 14.7 A (~ 248 A/m2), 30 s, respectively. By using the 200 keV accelerator, the megawatt-class negative hydrogen beam has also been achieved (135.9 keV, 8.9 A, 8 s). The whole process of the gas neutralization of negative ion beam, electric removal of residual ions, and beam transport have been demonstrated experimentally.

Recent Patents for Optical Component Polishing Technology

Background: Fluid jet polishing, Airbag polishing, Magnetorheological polishing, and Ion beam polishing are all emerging polishing technologies commonly used for processing optical components, which can be used to manufacture various high-precision and high-quality optical components. Their processing mechanisms are to remove materials by high-speed jet impacting the surface of the workpiece; using a spherical flexible airbag as a polishing tool to obtain a super-smooth surface by adding polishing fluid containing fine abrasives; controlling the magnetic field to form a flexible polishing belt of Magnetorheological fluid, which undergoes relative motion to achieve material removal; using a neutral ion beam to bombard the surface of the workpiece to remove workpiece surface atoms and obtain a super-smooth surface through energy and momentum transfer between ions and workpiece atoms. The characteristics and removal mechanisms of various polishing methods are briefly described, and the advantages and key issues that need to be addressed for each polishing technology are pointed out, and suggestions and prospects for future research, application, and development are also given. Objective: In order to meet the growing demand for high-precision optical components and mass-manufactured optical components, this paper aims to study the various polishing technologies that are continuously improved and suggestions are made for future research and development directions. Method: This paper reviews the current patents related to various polishing technologies, such as Fluid jet polishing, Airbag polishing, Magnetorheological polishing and Ion beam polishing and other representative polishing technologies. Results: Various types of polishing technologies have been studied, and the main problem at the current stage is the lack of accuracy. The problem of the polishing device, such as fluid jet polishing, only considers the change of fluid velocity and pressure to construct the removal function model for analysis and summary, and suggestions for the latest research progress of various polishing technologies are given. Conclusion: The improvement and optimization of polishing technology are beneficial to improving the polishing accuracy and quality of optical components and have great help for the future development of related patents.

Fast ion relaxation in ITER mediated by Alfvén instabilities

We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

Characterization of an RF-excited broad beam ion source operated with a mixture of CHF3 and O2

The composition and ion energy distributions of the main ion species of an ion beam were recorded and analyzed. The RF-type broad beam ion source was operated with a mixture of CHF3 and O2. A plasma bridge neutralizer operating with Ar was employed for ion beam neutralization. The data were collected with an energy-selective mass spectrometer (ESMS). The mass spectrum showed numerous ion species, beginning with ionized molecules, dissociation products of the process gases and products from reactions with background gas and the plasma discharge vessel, and the extraction system. For a quantification of the ion beam composition, the mass dependent transmission functions for two ESMS were determined. The ion energy distributions show that, in comparison to operation with inert gases, there are additional slower ions present. These ions can be related to dissociation processes outside of the ion beam source. As a result of their typically lower etching yield, these slower ions affect the etching behavior.

Improving the neutralization of a pulsed ion beam by electron cooling and accumulation: A kinetic study

Direct electron injection, such as through electron-emitting filaments, is usually difficult to neutralize ion beams to a very high degree. In this paper, the possibility of a pulsed ion beam achieving very high neutralization through the continuous accumulation of cold electrons is investigated using a two-dimensional particle-in-cell code. Three schemes of electron injection, namely, single-point injection, periodic point-source injection, and periodic line-source injection, are numerically studied and compared. The simulations show that even if an excess of electrons are injected, the single-point electron source is difficult to neutralize the ion beam pulse to exceed 90%, consistent with existing experiments. It is found that the spontaneous cooling mechanism of neutralizing electrons is able to improve the neutralization of the ion beam to a certain extent, but it requires a lot of time. By using a smaller injection current, the latter two injection schemes not only effectively suppress solitary waves, but more importantly, they continuously provide cold electrons that can accumulate inside the ion beam, thereby significantly improving the neutralization of the ion beam in a short period of time. The results show that periodic line-emission sources can neutralize the ion beam to over 99%, but periodic point-emission sources exhibit relatively poor neutralization performance due to their higher virtual-cathode potential. The research results can provide a reference for the design of neutralizing sources in applications that pursue very high neutralization of ion beam pulses, such as heavy ion fusion accelerators.

Accelerated neutral atom beam (ANAB) and gas clustered ion beam (GCIB) treatment of implantable device polymers leads to decreased bacterial attachment in vitro and decreased inflammation in vivo

Accelerated neutral atom beam (ANAB) and gas clustered ion beam (GCIB) treatment of implantable device polymers leads to decreased bacterial attachment in vitro and decreased inflammation in vivo

Influence of positive ions on the beamlet optics for negative-ion neutral beam injectors

Neutral beam injectors are based on the neutralization of ion beams accelerated at the desired energy. In the case of the ITER heating and diagnostic neutral beams, the target heating power translates into stringent requirements on the acceptable beamlet divergence and aiming to allow the beam to reach the fusion plasma. The beamlets composing the accelerated beam are experimentally found to feature a transverse velocity distribution exhibiting two Gaussian components: the well-focused one is referred to as the core component while the rest of the beam, the halo, describes beam particles with much worse optics. The codes that simulate beam extraction and acceleration usually assume that the negative ions move towards the plasma meniscus with a laminar flow (no transverse velocity) or that the transverse velocity distribution can be modelled as a Maxwellian and that the current density is uniformly illuminating the meniscus; under such approximations, the presence of highly divergent components cannot be explained. In this work, we develop a simple test-particle tracing code with Monte Carlo collisions, named ICARO (for Ions Coming Around), to study the transport of negative ions in the extraction region and derive the spatial and velocity distribution of the negative ions at the meniscus (i.e. the plasma boundary where a beamlet is extracted). In particular, the origin of the beamlet halo and its dependence on the source parameters are discussed, highlighting as a key parameter the energy distribution of positive ions in the source plasma.

Calibrating beam fluxes of a low-energy neutral atom beam facility.

Scientific detection and imaging instruments for low-energetic neutral atoms (ENA) onboard spacecraft require thorough pre-flight laboratory calibration against a well-characterized neutral atom beam source. To achieve this requirement, a dedicated test facility is available at the University of Bern, which is equipped with a powerful plasma ion source and an ion beam neutralization stage. Using surface neutralization, low-energy neutral atom beams of any desired gas species can be produced in the energy range from 3keV down as low as 10eV. As the efficiency of the neutralization stage is species and energy dependent, the neutralizer itself needs to be calibrated against an independent reference. We report on the calibration and characterization of this neutral atom beam source using our recently developed Absolute Beam Monitor (ABM) as a primary calibration standard. The ABM measures the absolute ENA flux independent of neutral species in the energy range from 10eV to 3keV. We obtain calibration factors of a few 100cm-2s-1 pA-1, depending on species at beam energies above about 100eV, and a power-law decrease for energies below 100eV. Furthermore, the energy loss of neutralized ions in the surface neutralizer is estimated from time-of-flight measurements using the ABM. The relative energy loss increases with ENA energy from low levels near zero up to 20%-35% at 3keV, depending on atomic species. Having calibrated our neutral beam source allows for accurate calibration of ENA space instruments.

High-order coupling of shear and sonic continua in JET plasmas

A recent model coupling the shear-Alfvén and acoustic continua, which depends strongly on the equilibrium shaping and on elongation in particular, is employed to explain the properties of Alfvénic activity observed on JET plasmas below but close to the typical frequency of toroidicity-induced Alfvén eigenmodes (TAEs). The frequency gaps predicted by the model result from high-order harmonics of the geodesic field-line curvature caused by plasma shaping (as opposed to lower-order toroidicity) and give rise to high-order geodesic acoustic eigenmodes (HOGAEs), their frequency value being close to one-half of the TAEs one. The theoretical predictions of HOGAE frequency and radial location are found to be in fair agreement with measurements in JET experiments, including magnetic, reflectometry and soft x-ray data. The stability of the observed HOGAEs is evaluated with the linear hybrid magnetohydrodynamic/drift-kinetic code CASTOR-K, taking into account the energetic-ion populations produced by neutral beam injection and ion cyclotron resonance heating systems. Wave-particle resonances, along with drive/damping mechanisms, are also discussed in order to understand the conditions leading to HOGAEs destabilisation in JET plasmas.

Final design and manufacturing of the neutralizer for the negative ion-based neutral beam injector test facility of CRAFT

The Comprehensive Research Facility for Fusion Technology (CRAFT) is a highly integrated research and development base. It consists of 20 different demonstrating or testing facilities that address most of the critical technologies and systems for the fusion reactor in China. A negative ion-based neutral beam injector (NNBI) test facility is established in CRAFT. A gas neutralizer is applied to the CRAFT NNBI test facility for the neutralization of the negative hydrogen ion beam. The target of the neutralization efficiency is more than 50%. The length of the neutralizer is 3 m. The neutralizer is divided into two adjacent vertical channels and the total opening is 1.7 m × 0.32 m. The neutralizer is mainly composed of channel plate, support structure, beam scraper, cooling water pipes, and gas inlet pipes. Based on the design requirements of the neutralizer, the support structure and cooling water pipes were optimized, to reduce the total weight. The structural strength of all components was simulated and compared to the preliminary design. The manufacturing of the neutralizer is ongoing. Some manufacturing challenges have been overcome, such as deep hole drilling. Several 1.7 m long cooling channels were successfully drilled in a copper panel and the ratio of depth to diameter was 94.

Excitation of surface waves in 3D ion beam neutralization

Neutralization of beams with 2D and 3D geometries by the electrons emitted from an external source is studied using particle-in-cell simulations. Our work reveals that the high-energy electrons excite Trivelpiece–Gould (TG) surface waves in the beams with 3D axisymmetric geometries. These high-energy electrons are generated because of a large amplitude electrostatic solitary wave (ESW) that forms near the electron source and has an electric potential amplitude more than three times the electron thermal energy. We also find that surface wave excitation only happens when the beam radius is large enough at the ion source to attract enough electrons that could form the large amplitude ESW. A comparison of the 3D TG surface wave dispersion relation with an expression for 2D surface waves reveals that they become excited in 3D axisymmetric but not in 2D planar beam because of a higher phase speed requirement in the latter case.

Tomography of fast ion distribution function under neutral beam injection and ion cyclotron resonance heating on EAST

In magnetic confinement fusion devices, velocity-space tomography of fast-ion velocity distribution function is crucial for investigating fast-ion distribution and transport. In the neutral beam injection (NBI) and ion cyclotron resonance heating (ICRF) synergistic heating experiments in Experimental Advanced Superconducting Tokamak (EAST), high-energy particles with energy exceeding the particle energy in NBI are observed. Simulations of synergistic effect on fast-ion velocity distribution function given by TRANSP also show the existence of particles with energy higher than the particle energy in NBI. To investigate the behaviors of fast ion distribution and calculate the velocity distribution functions under different heating conditions, the first-order Tikhonov regularization tomographic inversion method with higher inversion accuracy is introduced by comparing various regularization techniques. The limitations of the dual-view fast-ion D<sub>α</sub> (FIDA) diagnostic measurements in velocity space are addressed by incorporating prior information such as null measurement and the known peaks and effectively mitigate the occurrence of artifacts. This method is first employed in the case of NBI heating. The NBI peak is successfully reconstructed at the expected location in velocity space, which shows significant improvement in the inversion results. In order to further validate the synergistic effect of NBI-ICRF heating and study the mechanism of fast ion distribution under synergistic heating, the combination of FIDA and neutron emission spectrometer (NES) is applied to the first-order Tikhonov regularization tomographic inversion method for enhancing the coverage of velocity space, through which the issue of artifacts in the inversion results is significantly improved, and thus the precision of the obtained fast-ion velocity distribution functions is enhanced. Based on the benefit described above, the method of combining NES diagnosis and FIDA diagnosis is used to obtain fast-ion velocity distribution functions in the NBI and ICRF synergistic heating discharge. The synergistic heating effect is manifested in the fast-ion velocity distribution. The availability of this inversion method in reconstructing fast-ion velocity distribution functions during high-performance operation of NBI-ICRF synergistic heating in the EAST experiment is confirmed. In the next-step EAST research, high performance discharge will demand more efficiency NBI and ICRF synergistic heating, the present work builds the stage for investigating the underlying mechanism of synergistic heating and the intricate behaviors associated with fast ion distribution and transport.

Kinetic modeling of solitary wave dynamics in a neutralizing ion beam

In this work, we characterize the formation and evolution of electrostatic solitary waves (ESWs) in the space-charge neutralization of ion beams using particle-in-cell simulations. These waves become excited when the electrons emitted from an external filament source initiate a two-stream instability in the beam. We show that such electrostatic waves become excited in both two-dimensional (2D) and three-dimensional (3D) beams with different shapes and sizes. Through a 1D Bernstein–Greene–Kruskal (BGK) analysis of the 2D beam, we find that the non-Maxwellian nature of the beam electrons gives rise to large-sized ESWs that are not predicted by BGK theory since it assumes a Maxwellian electron velocity distribution in the beam. Finally, we show that a 1D BGK theory is inadequate to describe ESWs in 3D beams because of complex electron trajectories.

Overview of initial negative triangularity plasma studies on the ASDEX Upgrade tokamak

An overview of results of negative triangularity (NT) studies from the ASDEX Upgrade tokamak (AUG) is given. Moderate values of the triangularity in the range of δ average ≈ −0.2 have been obtained. Results from both unfavorable and favorable magnetic configuration in terms of high-confinement mode (H-mode) access are presented. In unfavorable configuration, the ion grad B drift points away from the active x-point, which usually results in a factor of 2 higher power thresholds to enter H-mode. In unfavorable configuration, low-confinement mode (L-mode) operation is confirmed at high auxiliary heating input powers MW at low collisionality. Under these conditions, plasmas with positive triangularity (PT) are usually in H-mode. Highest plasma energies have been obtained using a combination of neutral beam injection, electron- and ion cyclotron resonance heating, where values of the plasma beta of β N ≈ 1.7 have been reached in L-mode. While discharges with mixed auxiliary heating power input show moderate confinement on L-mode levels, plasmas heated with electron cyclotron resonance heating show H-mode like confinement. Operation in favorable configuration, however, leads to H-mode with type-I edge localized modes (ELMs). The power required for the L–H transition is comparable to the typical L–H power threshold in PT cases ( MW), and good confinement of typical PT H-modes is obtained. The ELM frequency is demonstrated to reach 700 Hz. These results show a clear effect of the NT configuration on the ELM behavior in H-modes, and that on AUG, unfavorable configuration is an essential ingredient if H-mode operation is to be avoided in NT plasmas. So far, in all AUG NT plasmas, power degradation is observed. Initial gyrokinetic simulations of the edge plasma in unfavorable configuration reveal the ion-temperature gradient mode as the fastest growing mode with a subdominant trapped-electron-mode.

Fusion by beam ions in a low collisionality, high mirror ratio magnetic mirror

The classical problem of neutral beam ions slowing down in a magnetic mirror geometry is revisited to provide predictive capability for the new Wisconsin HTS Axisymmetric Mirror under development at the University of Wisconsin. A Fokker–Planck model named fast beam ion solver (FBIS) is developed to include the spatial non-uniformity of a physical mirror geometry. The mathematical framework allows for efficient orbit averaging of the pitch-angle scattering operator, and permits a determination of the axial profile of the ambipolar potential confining the electrons. The numerical results from FBIS are consistent with earlier work, but further show how mirror-ratio and a near square magnetic well optimizes the fusion gain. The numerical results are also applied to inform the conceptual design of WHAM++, a low capital-cost breakeven-class magnetic mirror device.

Simulations of neutral beam injection and ion cyclotron resonance heating synergy in high power EAST scenarios.

The EAST plasmas heated with deuterium neutral beam injection and ion cyclotron resonance heating (ICRH) have been simulated by the TRANSP code. The analysis has been conducted using the full wave solver TORIC5, the radio frequency (RF)-kick operator, and NUBEAM to model the RF heating effects on fast ion velocity distribution. In this work, we present several simulated results compared with experiments for high power EAST scenarios, indicating that the interactions between ICRH and fast ions can significantly accelerate fast ions, which are confirmed by the increased neutron yield and broadened neutron emission spectrum measurements.

Excitation of long electrostatic solitary waves in ion beam neutralization process

Unusually long electrostatic solitary waves (ESWs) are discovered in a particle-in-cell simulation study of the process of ion beam neutralization by electron emission from a filament. These ESWs are long because the density perturbation responses to the potential wells created by the ESWs are very small. The density perturbation is small because the trapped (positive) and untrapped (negative) electron density perturbations nearly compensate each other because of a non-Maxwellian electron velocity distribution in the beam.

Doublet splitting of fusion alpha particle driven ion cyclotron emission

Ion cyclotron emission (ICE) originating from confined populations of fast ions in toroidal fusion plasmas is an important non-invasive, passive diagnostic for current and next-generation devices. The ability to model ICE signals accurately is an essential step towards inferring the characteristics of confined energetic alpha particles or fast neutral beam ions and background plasma. In this paper, the linear growth rates of the magnetoacoustic cyclotron instability, which is the leading explanation for ICE, are calculated in high-resolution 2D space for parameters corresponding to JET Pulse No. 26148. These calculations further highlight the origin of doublet splitting of the cyclotron harmonics and the perceived need to invoke nonlinear interactions with 1D3V particle-in-cell studies to account for lower harmonics with substantial amplitudes. The inclusion of signals of ICE from a well-resolved range of propagation angles can also account for these effects.

Integrative simulation of a 2 cm electron cyclotron resonance ion source with full particle-in-cell method

A two-dimensional full particle-in-cell model is developed to simulate the operation process of a 2 cm electron cyclotron resonance ion thruster. The simulation regions include the discharge chamber, grid system, and plume, and the processes of neutral gas injection, plasma discharge, and ion beam extraction are simulated continuously. Neutrals are treated as particles and managed through an adaptive particle management algorithm. The collisions between charged and neutral particles are approached by the method of combining the Monte Carlo and directly simulation Monte Carlo collisions. The integrative simulation provides a global perspective to learn the gridded ion thruster, such as the Child-Langmuir sheath, acceleration grid erosion, plasma bridge, and so on. In this paper, the response of the plasma system to the accelerating voltage is analyzed, which is related to the bulk plasma, antenna floating voltage, and Child-Langmuir sheath. Besides, an ion density wave is observed in plume, which will increase the current of acceleration grid and cause ion backflow.