Abstract
We present the MPI-based implementation of the method of 2D and 3D calculations of the evolution of temperature fields and the phase transitions dynamics in materials irradiated by high-energy heavy ion and by pulsed ion beams. We utilize a modified thermal spike model based on a system of coupled heat conductivity equation describing thermal processes in the electron and ion subsystems of the target sample. Such equations are numerically solved in the cylindrical coordinate system in axially symmetric (2D) and axially nonsymmetrical (3D) cases. The dynamics of phase transitions is realized on the basis of Stephan’s problem in the framework of the enthalpy approach. The mathematical formulation of the problem is given; a numerical scheme is described; a parallel algorithm is presented. The numerical results confirm the efficiency of our approach and the corresponding MPI/C++ computer code. It is shown that the results of numerical simulations are in agreement with experimental estimations of the track sizes which appear in the target samples exposed to heavy ion.
Email this article 
