UDC 519.633.6, 519.688 Computational Scheme for Solving Heat Conduction Problem in Multilayer Cylindrical Domain A. S. Ayriyan∗, J. Buˇ H. Grigorian∗§, J. Pribiˇ Technical University of Kosice Letna str. 9, Kosice, Slovak Republic, 04001 ‡ Laboratory of High Energies Joint Institute for Nuclear Research Joliot-Curie str. 6, Dubna, Moscow region, Russia, 141980 § Department of Theoretical Physics Yerevan State University Alek Manoogian str. 1, Yerevan, Armenia, 0025 The computational scheme for solving heat conduction problem with periodic source function in multilayer cylindrical domain is suggested. <...> The domain has a non-trivial geometry and the thermal coefficients are non-linear functions of temperature and have discontinuity of the first kind at the borders of the layers. <...> The computational scheme is based on an algorithm for solving difference problem using the explicit-implicit method. <...> The OpenCL realization of the suggested algorithm for calculations performed on a GPU is also compared to calculations performed using a CPU. <...> It is shown that the scheme can be successfully applied to simulations of thermal processes in pulsed cryogenic cell, which is intended for pulse feeding the working gases into the working space of the ion source within the millisecond range. <...> The results are given for a simulation of one of the particular cell structures, which is assumed to correspond to the practical realization. <...> The computational scheme can be used for the optimization problem of the cell model parameters. <...> Key words and phrases: heat evolution, periodical heating source, multilayer cylindrical structure, finite-difference scheme, OpenCL realization. <...> The main goal of this work is to suggest a computational scheme for solving heat conduction problem with a periodic source in multilayer cylindrical domain. <...> This problem describes the thermal processes inside a pulsed cryogenic cell in frame of previously suggested model [1, 2]. <...> The function of this cell is pulse feeding (in the millisecond range) the working gases into the working space of the ion source [3]. <...> The work was supported by RFBR grants 14-01-00628 and 14-01-31227. z0 zmax r0 * r1 * r2 * 4 rmax * sa Jr.†, E. E. Donets‡, s† ∗ Laboratory of Information Technologies Joint Institute for Nuclear Research <...>