The numerical methods used in the full particle-orbit following SPIRAL code are described and a number of physics studies performed with the code are presented to illustrate its capabilities. The SPIRAL code is a test-particle code and is a powerful numerical tool to interpret and plan fast-ion experiments in tokamaks. Gyro-orbit effects are important for fast ions in low-field machines such as NSTX and to a lesser extent in DIII-D. A number of physics studies are interlaced between the description of the code to illustrate its capabilities. Results on heat loads generated by a localized error-field on the DIII-D wall are compared with measurements. The enhanced Triton losses caused by the same localized error-field are calculated and compared with measured neutron signals. Magnetohydrodynamic (MHD) activity such as tearing modes and toroidicity-induced Alfven eigenmodes (TAEs) have a profound effect on the fast-ion content of tokamak plasmas and SPIRAL can calculate the effects of MHD activity on the confined and lost fast-ion population as illustrated for a burst of TAE activity in NSTX. The interaction between ion cyclotron range of frequency (ICRF) heating and fast ions depends solely on the gyro-motion of the fast ions and is captured exactly in the SPIRAL code. A calculation of ICRF absorption on beam ions in ITER is presented. The effects of high harmonic fast wave heating on the beam-ion slowing-down distribution in NSTX is also studied.

/pdf/a-description-of-the-full-particle-orbit-following-spiral-4v6odu1244.pdf

A Description of the Full Particle Orbit Following SPIRAL Code for Simulating Fast-ion Experiments in Tokamaks

A fast-ion D-alpha (FIDA) diagnostic has been developed for the fully tungsten coated ASDEX Upgrade (AUG) tokamak using 25 toroidally viewing lines of sight and featuring a temporal resolution of 10 ms. The diagnostic's toroidal geometry determines a well-defined region in velocity space which significantly overlaps with the typical fast-ion distribution in AUG plasmas. Background subtraction without beam modulation is possible because relevant parts of the FIDA spectra are free from impurity line contamination. Thus, the temporal evolution of the confined fast-ion distribution function can be monitored continuously. FIDA profiles during on- and off-axis neutral beam injection (NBI) heating are presented which show changes in the radial fast-ion distribution with the different NBI geometries. Good agreement has been obtained between measured and simulated FIDA radial profiles in MHD-quiescent plasmas using fast-ion distribution functions provided by TRANSP. In addition, a large fast-ion redistribution with a drop of about 50% in the central fast-ion population has been observed in the presence of a q = 2 sawtooth-like crash, demonstrating the capabilities of the diagnostic.

/pdf/fast-ion-d-alpha-measurements-at-asdex-upgrade-4zxcfaz1in.pdf

Fast-ion D-alpha measurements at ASDEX Upgrade

We compute tomographies of 2D fast-ion velocity distribution functions from synthetic collective Thomson scattering (CTS) and fast-ion Dα (FIDA) 1D measurements using a new reconstruction prescription. Contradicting conventional wisdom we demonstrate that one single 1D CTS or FIDA view suffices to compute accurate tomographies of arbitrary 2D functions under idealized conditions. Under simulated experimental conditions, single-view tomographies do not resemble the original fast-ion velocity distribution functions but nevertheless show their coarsest features. For CTS or FIDA systems with many simultaneous views on the same measurement volume, the resemblance improves with the number of available views, even if the resolution in each view is varied inversely proportional to the number of views, so that the total number of measurements in all views is the same. With a realistic four-view system, tomographies of a beam ion velocity distribution function at ASDEX Upgrade reproduce the general shape of the function and the location of the maxima at full and half injection energy of the beam ions. By applying our method to real many-view CTS or FIDA measurements, one could determine tomographies of 2D fast-ion velocity distribution functions experimentally.

/pdf/tomography-of-fast-ion-velocity-space-distributions-from-37se00atlb.pdf

Tomography of fast-ion velocity-space distributions from synthetic CTS and FIDA measurements

The velocity-space observation regions and sensitivities in fast-ion Dα (FIDA) spectroscopy measurements are often described by so-called weight functions. Here we derive expressions for FIDA weight functions accounting for the Doppler shift, Stark splitting, and the charge-exchange reaction and electron transition probabilities. Our approach yields an efficient way to calculate correctly scaled FIDA weight functions and implies simple analytic expressions for their boundaries that separate the triangular observable regions in (v||, v⊥)-space from the unobservable regions. These boundaries are determined by the Doppler shift and Stark splitting and could until now only be found by numeric simulation.

/pdf/on-velocity-space-sensitivity-of-fast-ion-d-alpha-2wsxzf5u18.pdf

On velocity-space sensitivity of fast-ion D-alpha spectroscopy

The neutral-beam induced D(α) emission spectrum contains a wealth of information such as deuterium ion temperature, toroidal rotation, density, beam emission intensity, beam neutral density, and local magnetic field strength magnitude |B| from the Stark-split beam emission spectrum, and fast-ion D(α) emission (FIDA) proportional to the beam-injected fast ion density. A comprehensive spectral fitting routine which accounts for all photoemission processes is employed for the spectral analysis. Interpretation of the measurements to determine physically relevant plasma parameters is assisted by the use of an optimized viewing geometry and forward modeling of the emission spectra using a Monte-Carlo 3D simulation code.

/pdf/active-spectroscopic-measurements-of-the-bulk-deuterium-fr5k0ya67d.pdf

Active spectroscopic measurements of the bulk deuterium properties in the DIII-D tokamak (invited).

A code that models signals produced by charge-exchange reactions between fast ions and injected neutral beams in tokamak plasmas is described. With the fast- ion distribution function as input, the code predicts the efflux to a neutral particle analyzer (NPA) diagnostic and the photon radiance of Balmer-alpha light to a fast- ion Dα (FIDA) diagnostic. Reactions with both the primary injected neutrals and with the cloud of secondary "halo" neutrals that surround the beam are treated. Accurate calculation of the fraction of neutrals that occupy excited atomic states (the collisional- radiative transition equations) is an important element of the code. Comparison with TRANSP output and other tests verify the solutions. Judicious selection of grid size and other parameters facilitate efficient solutions. The output of the code has been validated by FIDA measurements on DIII-D but further tests are warranted.

/pdf/a-code-that-simulates-fast-ion-da-and-neutral-particle-399cvxkl14.pdf

Y. Luo

Papers

A code that simulates fast-ion Dα and neutral particle measurements