New techniques for calculating heat and particle source rates due to neutral beam injection in axisymmetric tokamaks

A set of numerical techniques for calculating heat and particle source rates due to neutral beam injection in axisymmetric tokamaks is described. While these techniques consume a substantial amount of computer time, they take into account a number of significant, and normally neglected, effects. Examples of these effects are reionization of escaping charge exchanged beam particles, finite fast ion orbit excursions, beam deposition through collisions of beam neutrals with circulating beam ions, and the transport of thermal neutrals in the plasma due to charge changing collisions with beam ions.