Numerical simulation of the impurity peak positions in TFR 400 discharges

A numerical code which solves evolution equations for the densities of each ionization level of oxygen and molybdenum is used for simulating the stationary state of T.F.R. discharge. Using neo classical (Pfirsch-Schluter) fluxes and anomalous 'spreading' fluxes, gives the positions of oxygen peaks inwardly displaced from those calculated using the coronal equilibrium model, in reasonable agreement with the experimental measurements. Reasonable agreement with the experiment for the emitted photon is also found. In addition the positions of the maxima of the molybdenum peaks are in agreement with both the coronal equilibrium model and experimental observations