Improved equilibrium reconstructions by advanced statistical weighting of the internal magnetic measurements

Abstract

In a Tokamak the configuration of the magnetic fields remains the key element to improve performance and to maximise the scientific exploitation of the device. On the other hand, the quality of the reconstructed fields depends crucially on the measurements available. Traditionally in the least square minimisation phase of the algorithms, used to obtain the magnetic field topology, all the diagnostics are given the same weights, a part from a corrective factor taking into account the error bars. This assumption unduly penalises complex diagnostics, such as polarimetry, which have a limited number of highly significant measurements. A completely new method to choose the weights, to be given to the internal measurements of the magnetic fields for improved equilibrium reconstructions, is presented in this paper. The approach is based on various statistical indicators applied to the residuals, the difference between the actual measurements and their estimates from the reconstructed equilibrium. The potential of the method is exemplified using the measurements of the Faraday rotation derived from JET polarimeter. The results indicate quite clearly that the weights have to be determined carefully, since the inappropriate choice can have significant repercussions on the quality of the magnetic reconstruction both in the edge and in the core. These results confirm the limitations of the assumption that all the diagnostics have to be given the same weight, irrespective of the number of measurements they provide and the region of the plasma they probe.

Investigating the thermal stability of highly radiative discharges on JET with a new tomographic method

“…To obtain reasonable results, it is therefore indispensable to profit from the information about the magnetic topology. Unfortunately, the identification of the magnetic equilibrium is another ill-posed problem and therefore also the magnetic reconstructions have uncertainties [15][16][17]. On the other hand, the calculation of the uncertainties in the tomographic reconstructions, with the iterative formula (5), considers only the propagation of the instrumental errors in the inversion process.…”