Small angle reconstruction of plasma luminosity for spherical tokamak

Abstract

Публикация посвящена применению методов вычислительной геометрии, интервального анализа и линейного программирования к задачам физики управляемого термоядерного синтеза. Рассмотрены геометрические аспекты проблемы, получены проекции светимостей различных объемов сферического токамака на плоскость матричного детектора, изучены изображения предполагаемых макроскопических структур и микроскопических включений. Для набора модельных распределений светимости объема токамака поставлена задача восстановления сигнала. Решение получено с использованием задач линейного программирования. The problems of reconstruction of plasma luminosity are important for physics and technology of power plants-tokamaks. The Globus-M research tokamak obtained a large amount of data using a matrix detector in pinhole camera geometry. From the mathematical point of view, finding the luminosity for different regions of the plasma volume according to the matrix detector is an inverse problem related to the field of integral geometry. An essential feature of the particular task is the use of a single fixed camera with a small viewing angle. In this regard, application of methods of harmonic analysis of data is not enough. The paper investigates the geometric aspects of the problem. In the general view, a threedimensional object is projected onto a two-dimensional plane through a diaphragm. Under the assumption of azimuthal symmetry, there is a central projection of the luminosity of the body of rotation onto a flat matrix detector. The initial information for the calculation is the plasma boundary obtained from magnetic sensors. There is no reliable information about the internal structure of the plasma, so its division into regions of the equal luminosity is not unambiguous. The paper presents an algorithm for finding the projections of the luminosity of plasma volumes on the plane of the matrix detector. A set of model direct problems for the construction of algorithms for their recognition according to the detector data was investigated. Images of supposed macroscopic structures and microscopic inclusions were obtained. The methodological basis of the work is the use of interval analysis methods for solving geometric and algebraic problems. This approach allows obtaining qualitative and quantitative results that takes into account the uncertainty of the input data with the minimum amount of computational costs. Algebraic solvability is investigated in the interval formulation using response functionality. Solutions for a set of test problems are obtained, which demonstrate the availability of successful reconstruction for real data. An important result of the study is an information about the presence of uncertainties in geometric data and related calculations by obtaining results about the luminosity of the plasma by solving linear programming problems.