Математическое моделирование плавления вольфрама при воздействии лазерного импульса

Работа посвящена математическому моделированию процесса плавления в образце под воздействием импульсной тепловой нагрузки на основе решения двухфазной задачи Стефана. Численная модель основана на подходе Самарского, что позволяет не выделять свободную границу во время расчета. Учет аксиально-симметричной геометрии позволил показать, что на испарение расходуется около четверти падающей энергии в центре области расплава. Это в пять раз больше, чем дают оценки, основанные на решении одномерного уравнения теплопроводности. В случае учета испарения вещества получено хорошее соответствие расчетных и экспериментальных температуры остывающей поверхности и скорости сужения области расплава. Результаты математического моделирования подтвердили существование режима охлаждения испарением при нагреве вольфрама электронным пучком существенно выше порога плавления.

The paper is devoted to mathematical modeling of the melting process in a sample under the influence of a pulsed thermal load based on the solution of the two-phase Stefan problem. The free boundary is ignoring during the calculation, since the numerical model is based on the Samarsky approach. The calculation in axially symmetric geometry allowed us to show that about a quarter of the incident energy is consumed in the center of the melt region. This is five times more than estimates based on the solution of the one-dimensional heat equation give. Considering the evaporation of the substance a good correspondence between the calculated and experimental temperatures of the cooling surface and the rate of narrowing of the melt region is obtained. The results of mathematical modeling confirmed the existence of an evaporation cooling mode when tungsten is heated by an electron beam significantly above the melting threshold.