Modelling 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 passes to its stationary solution after some time 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 modelling.

(based on the PhD thesis)