Abstract
A computational investigation of the dynamics of gas expansion due to intense nanosecond laser evaporation into vacuum has been carried out. The problem is solved in a one-dimensional approximation, which simplifies calculations and at the same time allows one to analyze the main features of the expansion dynamics. For analysis we use three different approaches. Two of them are based on kinetic analysis via the direct simulation Monte Carlo (DSMC) method and numerical solution of the model Bhatnagar-Gross-Krook (BGK) equation. The third one focuses on derivation of an analytical continuum solution. Emphasis is placed on the analysis of the velocity distribution function and the average energy of particles passing through the time-of-flight detector on the normal to the evaporation surface, which is important for interpreting experimental measurements. The formulated problem is quite difficult as the considered flow is time-dependent, contains discontinuities in boundary conditions and involves large variations of local Knudsen numbers as well as steep gradients of the velocity distribution function. Data were obtained on the particle energy in the time-of-flight distribution for the range of regimes from the free molecular flow to continuum one. The maximum attainable average energy of particles in the time-of-flight distribution is determined. The non-monotonicity of the energy increase was found, which is explained based on analysis of the velocity distribution of particles.