Numerical model of high-power transient heating of tungsten with considering of various erosion effects

Surface melting of tungsten under exposure to a pulsed electron beam was simulated numerically, the evaporation process taken into account. The calculation is based on the experimental time dependence of the total beam power. The model of the tungsten heating process is based on solving the two-phase Stefan problem. The position of the phase boundary depends on discontinuous time- and space-nonlinear coefficients and boundary conditions. The aim of the study is to provide a detailed resolution of the heat flow deep into the material with a fine spatial grid step. As compared with the size of the tungsten plate, the heating depth is very small. The problem statement under consideration is multiscale. Further expansion of the model involves taking into account microcracks. Micro-cracks occur during the cooling process after exposure and affect the temperature of the tungsten surface during the subsequent heating process. The article presents a modeling of cracks of different geometries typical for this process. The results of the calculations correlate with the experimental data obtained on the experimental test facility BETA at BINP SB RAS.