Dynamic transport simulation code including plasma rotation and radial electric field

A new one-dimensional transport code named TASK/TX, which is able to describe dynamic behavior of tokamak plasmas, has been developed. It solves simultaneously a set of flux-surface averaged equations composed of Maxwell's equations, continuity equations, equations of motion, heat transport equations, fast-particle slowing-down equations and two-group neutral diffusion equations. The set of equations describes plasma rotations in both toroidal and poloidal directions through momentum transfer and evaluates the radial electric field self-consistently. The finite element method with a piecewise linear interpolation function is employed with a fine radial mesh near the plasma surface. The Streamline Upwind Petrov-Galerkin method is also used for robust calculation. We have confirmed that the neoclassical properties are well described by the poloidal neoclassical viscous force. The modification of density profile during neutral beam injection is presented. In the presence of ion orbit loss, the generation of the inward radial electric field and torque due to radial current is self-consistently calculated. (c) 2007 Elsevier Inc. All rights reserved.