Моделирование лазерной абляции материалов можно проводить в рамках одномерного нестационарного уравнения теплопроводности в системе координат, связанной с движущимся фронтом испарения. Действие лазера на образец учитывается через функции источника в уравнении теплопроводности, задавая координатную и временную зависимости источника лазера. В предыдущих работах было проведено численное моделирование лазерной абляции материалов, возникающей под действием ультракоротких лазерных импульсов в полуограниченных образцах. Численные расчеты были выполнены с применением метода конечных разностей. Полученные результаты согласовались с результатами работ других авторов.
В настоящей работе проведено аналогичное численное исследование при непрерывном воздействии лазера постоянной интенсивности на образец. При этом нестационарное решение уравнения теплопроводности после некоторого времени переходит к его стационарному решению даже при учете температурной зависимости теплофизических параметров материала образца, в то время, когда аналитическое решение можно получить только при предположении независимости искомых параметров от температуры. Проведен сравнительный анализ полученных результатов с результатами аналитических решений. Динамика перехода к стационарному решению полученных аналитических решений выявляется только численным моделированием.
__Numerical Simulation of Stationary Waves of Laser Ablation of Materials in Semi-Bounded Samples __
Modeling of laser ablation of materials can be carried out within the framework of a one-dimensional non-stationary heat conduction equation in a coordinate system associated with a moving evaporation front. The action of the laser on the sample is taken into account through the source functions in the heat conduction equation, specifying the coordinate and time dependences of the laser source. In previous works, numerical simulation of laser ablation of materials arising under the action of ultrashort laser pulses in semi-bounded samples was carried out. Numerical calculations were performed using the finite difference method. The obtained results were consistent with the results of other authors.
In this paper, a similar numerical study was carried out with a continuous action of a constant-intensity laser on a sample. In this case, the non-stationary solution of the heat conductivity equation after some time passes to its stationary solution even when taking into account the temperature dependence of the thermophysical parameters of the sample material, while an analytical solution can only be obtained under the assumption that the sought parameters are independent of temperature. A comparative analysis of the obtained results with the results of analytical solutions was carried out. The dynamics of the transition to a stationary solution of the obtained analytical solutions are revealed only by numerical modeling.