A Study of the Characteristics of Plasma Generated by Infrared Pulse Laser-Induced Fused Silica

: When high energy infrared laser pulses are incident on fused silica, the surface of the fused silica is damaged and a laser-induced plasma is produced. Based on the theory of ﬂuid mechanics and gas dynamics, a two-dimensional axisymmetric gas dynamic model was established to simulate the plasma generation process of fused silica induced by a millisecond pulse laser. The results show that the temperature of the central region irradiated by the laser is the highest, and the plasma is ﬁrst produced in this region. When the laser energy density is 1.0 × 10 4 J/cm 2 and the pulse width is 0.2 ms, the maximum expansion velocity of the laser-induced plasma is 17.7 m/s. Under the same experimental conditions, the results of the simulation and experiment are in good agreement. With an increase in pulse width, the plasma expansion rate gradually decreases.


Introduction
With the popularization and application of a strong laser in various fields [1][2][3][4], optical elements as indispensable components of the laser system have received extensive attention. Because fused silica material has good light transmittance [5], it is widely used as an optical element in high power laser systems [6].
Fused silica material is one of the weakest components in a high power laser system. Its damage resistance directly affects the output performance of the laser. Even if there are subtle flaws, the performance of the laser will decrease. Therefore, it is of great significance to study the damage process of the interaction between the laser and fused silica [7,8]. R.W. Hopper et al. [9] established a theoretical model of laser damage induced by defects and studied the thermal stress damage caused by impurity defects in the target. F. Bolmeau et al. [10,11] used numerical simulation to study the laser energy deposition process caused by defects in fused silica.
The process of the interaction between the laser and fused silica is not only limited to material damage but also causes a series of processes such as phase transition, ablation, splashing, and plasma expansion. As a followup phenomenon of fused silica damage, plasma expansion has a high research value. However, the research on millisecond pulse laser-induced plasma generation of fused silica is quite scarce. Therefore, it is necessary to carry out related research on the plasma generation of fused silica induced by an infrared millisecond pulse laser.
As a high efficiency and commonly used numerical calculation method, the finite element analysis method is widely used in the treatment of various complex physical problems. M. Courtois et al. [12] used the finite element analysis method to numerically simulate and analyze the coupled heat transfer, fluid heat transfer, and thermal-mechanical coupling processes in laser welding. In this paper, the finite element method is used to simulate the plasma generation process of fused silica induced by the infrared millisecond

Materials and Methods
In the process of laser heating the target, the state of the material changed from solid to liquid and then changed to a gaseous state. The laser interacted with the ejected gaseous material, the atoms were ionized, the plasma was formed, and ablation craters appeared in the laser-irradiated target area. The condensed target material spontaneously splashed into the action of atmosphere under the high pressure generated by the high temperature, which was the physical process of plasma expansion.
In the simulation process, considering plasma generation and wave propagation at geometric junctions, we used free triangular meshes to divide the computational domain. We refined the grids in the handover area and the upper area. The maximum cell size of the grid in the dense area was 0.02 mm, the smallest cell was 1.5 × 10 −4 mm, the maximum cell growth rate was 1.08, and the curvature factor was 0.25. The time step used was the range of 0, 0.001, and 2, and the unit was millisecond. Figure 1 is a schematic diagram of the numerical simulation of the laser-induced plasma. The initial plasma was generated on the surface of the target and propagated against the direction of the laser under the action of the laser. simulate and analyze the coupled heat transfer, fluid heat transfer, and thermal-mechanical coupling processes in laser welding. In this paper, the finite element method is used to simulate the plasma generation process of fused silica induced by the infrared millisecond pulse laser; the results of the simulation and experiment are in good agreement, and the results have good convergence.

Materials and Methods
In the process of laser heating the target, the state of the material changed from solid to liquid and then changed to a gaseous state. The laser interacted with the ejected gaseous material, the atoms were ionized, the plasma was formed, and ablation craters appeared in the laser-irradiated target area. The condensed target material spontaneously splashed into the action of atmosphere under the high pressure generated by the high temperature, which was the physical process of plasma expansion.
In the simulation process, considering plasma generation and wave propagation at geometric junctions, we used free triangular meshes to divide the computational domain. We refined the grids in the handover area and the upper area. The maximum cell size of the grid in the dense area was 0.02 mm, the smallest cell was 1.5 × 10 −4 mm, the maximum cell growth rate was 1.08, and the curvature factor was 0.25. The time step used was the range of 0, 0.001, and 2, and the unit was millisecond. Figure 1 is a schematic diagram of the numerical simulation of the laser-induced plasma. The initial plasma was generated on the surface of the target and propagated against the direction of the laser under the action of the laser. The whole physical process of plasma generation and expansion satisfies conservation of mass, conservation of momentum, and conservation of energy. In the numerical simulation, the parameters of plasma changed with temperature T (where the units of T are K ) are shown in Table 1.  The whole physical process of plasma generation and expansion satisfies conservation of mass, conservation of momentum, and conservation of energy. In the numerical simulation, the parameters of plasma changed with temperature T (where the units of T are K) are shown in Table 1.
The following needed to be clear during the simulation: (1) The propagation velocity of the plasma is much lower than the speed of sound, which can be regarded as subsonic laminar flow. (2) Plasma generation and expansion during a small span of time without heat exchange with the outside and inside is a local thermal equilibrium. Throughout the physical process, the total mass of matter remained unchanged, which satisfied the law of conservation of mass: In the formula: u is the fluid velocity, ρ is the density. Momentum conservation equation: In the formula: p is the pressure, µ is the material viscosity. Energy conservation equation: In the formula: C P is the specific heat capacity, K is the thermal conductivity, Q Laser is the laser heat source, and Q Loss is the heat source that maintains the plasma expansion loss.
In the energy conservation equation, Q Laser is one of the source terms of the laser heat source; it is simplified to the Gaussian surface heat source distribution that varies with radius and temperature: where I 0 is the power density, r 0 is the laser spot radius, and A(T) is the dynamic absorption rate. τ(t) is the time distribution function of the laser pulse width; the expression is as follows: τ P is the laser pulse width. The vaporization and expansion process included the vapor evaporation process of the material leaving the target surface and the process of spraying outwards under continued heating. The mass loss of the surface layer of the fused silica material during the outward spraying led to the further formation of molten pits and continued to induce plasma. The quantity that affects the quality loss is the mass mobility [12]: where β is the diffusion coefficient, the value is assumed to be 1 at the beginning of the evaporation, and the value is 0 when it is stabilized, k B is Boltzmann's constant, P sat (T) is the saturated vapor pressure of the target steam [12,13], and the expression is: where P amb is the atmospheric pressure, L v is the latent heat of evaporation, T v is the evaporation temperature, and M SiO 2 is the molar mass of the fused silica. In Equation (3), the expression of Q Loss is: ρ l is the density of the material after melting, L V is the latent heat of vaporization, v r is the expansion speed of the plasma, and k is the energy loss rate in the system.
Computational domain boundary conditions:

Principle of the Experiment
The experimental measurement principle diagram of the plasma expansion of fused silica induced by the millisecond pulse laser is shown in Figure 2. In order to analyze the expansion law of the plasma generated by laser-induced f silica, the high-time-resolution optical shadow method was used to monitor the evol of the laser-induced plasma flow field. The time delay between the high-speed ca and the laser system was adjusted by using the DG645 digital delay generator. The h speed camera received the image. By analyzing the optical shadow images under diff time series conditions, the expansion and propagation process of the laser-ind plasma was obtained. In order to analyze the expansion law of the plasma generated by laser-induced fused silica, the high-time-resolution optical shadow method was used to monitor the evolution of the laser-induced plasma flow field. The time delay between the high-speed camera and the laser system was adjusted by using the DG645 digital delay generator. The high-speed camera received the image. By analyzing the optical shadow images under different time series conditions, the expansion and propagation process of the laser-induced plasma was obtained.

Analysis of the Plasma Temperature Field Induced by a Pulsed Laser in Fused Silica
Under millisecond pulse laser irradiation, the laser energy densities used to simulate the plasma generated by laser irradiation of fused silica were 8.0 × 10 3 J/cm 2 , the pulse width was 0.2 ms, and the irradiation radius was 0.25 mm.
The temperature field results of the numerical simulation of the plasma generated by the millisecond pulse laser induced fused silica were as follows: When the laser with an energy density of 8.0 × 10 3 J/cm 2 irradiated the surface of the fused silica material, the main form of plasma generated in the laser irradiation area was thermally induced plasma. Due to the laser heating, the surface of the fused silica material melted, vaporized, and splashed, forming a numerical simulation of the low-temperature thermally induced plasma temperature field as shown in the figure. Figure 3 shows the temperature field distribution of the laser-induced plasma at different times of the plasma. At the initial stage of laser action 0-0.06 ms, the temperature of the irradiated area rose rapidly and plasma began to be generated; the temperature gradually stabilized within the time range of 0.10-0.15 ms, the plasma slowly expanded outward, and thermally induced plasma was generated. In the time range of 0.25-0.30 ms, since the internal energy of the thermally induced plasma was converted into kinetic energy, the plasma expanded outwards and the temperature decreased accordingly.  Figure 4 is the curve of radial temperature change with time when the laser energy density was 8.0 × 10 3 J/cm 2 in the period of 0.005-0.030 ms. From Figure 4, it can be seen that the temperature change in the center area of the laser irradiation on the surface of the target was much greater than that of other areas. This is because the pulse laser has a   Figure 4 is the curve of radial temperature change with time when the laser energy density was 8.0 × 10 3 J/cm 2 in the period of 0.005-0.030 ms. From Figure 4, it can be seen that the temperature change in the center area of the laser irradiation on the surface of the target was much greater than that of other areas. This is because the pulse laser has a Gaussian distribution, and the laser energy is concentrated in the center, which makes the highest temperature rise in the irradiated center area the fastest. Figure 5 is the curve of radial temperature change with time when the laser energy density was 8.0 × 10 3 J/cm 2 in the period of 0.035-0.055 ms. Compared with Figure 4, it can be seen that the temperature in the central area of laser irradiation was stable at 0.030 ms, and plasma began to be generated. At this time, since the action time of the pulsed laser did not stop, the generation area of the plasma began to expand.       Figure 6 shows the flow field distribution of the plasma at different times. In the initial stage of pulse laser action, the velocity of the flow field in the action area increased gently (0-0.14 ms). In the time range of 0.14-0.20 ms, the plasma flow field velocity increased rapidly and reached maximum velocity at 0.2 ms, with a value of 17.7 m/s. The velocity of the plasma flow field decreased gradually between 0.24 ms and 0.28 ms. This was because the laser action stopped, and the laser irradiation did not continue to provide energy; the plasma could only rely on its own internal energy to be converted into kinetic energy to continue to maintain its outward expansion.
The plasma expansion velocity as a function of time t was found using where the expanding plasma front is measured to reach distance () Lt at time t , and t is a small interval of time typically; the unit of t we used here is μs [14]. Figure 7A shows the axial velocity distribution of plasma generated by the millisecond laser-induced fused silica with a pulse width of 0.2 ms; the dotted line is the experimental results, the solid line is the simulation results. Within 0-0.2 ms, the laser induced plasma continuously absorbed laser energy and expanded outward rapidly. Figure 7B shows that the axial plasma expansion rate was much greater than the radial expansion rate. The reason for this phenomenon is that the plasma expansion speed is determined by the laser energy. When the laser was incident perpendicular to the target surface, the laser energy distributed in the axial direction was much greater than the laser energy in the radial direction. The axial position made it easier for the plasma to move faster than  Figure 6 shows the flow field distribution of the plasma at different times. In the initial stage of pulse laser action, the velocity of the flow field in the action area increased gently (0-0.14 ms). In the time range of 0.14-0.20 ms, the plasma flow field velocity increased rapidly and reached maximum velocity at 0.2 ms, with a value of 17.7 m/s. The velocity of the plasma flow field decreased gradually between 0.24 ms and 0.28 ms. This was because the laser action stopped, and the laser irradiation did not continue to provide energy; the plasma could only rely on its own internal energy to be converted into kinetic energy to continue to maintain its outward expansion.
The plasma expansion velocity as a function of time t was found using where the expanding plasma front is measured to reach distance L(t) at time t, and t is a small interval of time typically; the unit of t we used here is µs [14]. Figure 7A shows the axial velocity distribution of plasma generated by the millisecond laser-induced fused silica with a pulse width of 0.2 ms; the dotted line is the experimental results, the solid line is the simulation results. Within 0-0.2 ms, the laser induced plasma continuously absorbed laser energy and expanded outward rapidly. Figure 7B shows that the axial plasma expansion rate was much greater than the radial expansion rate. The reason for this phenomenon is that the plasma expansion speed is determined by the laser energy. When the laser was incident perpendicular to the target surface, the laser energy distributed in the axial direction was much greater than the laser energy in the radial direction. The axial position made it easier for the plasma to move faster than in the radial position. Therefore, there was a large gap between the expansion speed of the axial plasma and the expansion speed of the radial plasma.  Figure 8 shows the evolution process of the plasma expansion based on the shadow method using a high-speed camera with a pulse width of 0.2 ms under the same laser energy density (which was 50 μs/frame). Through calculation, we obtained the maximum expansion velocity of the combustion wave in this process, about 18.2 m/s. The observed experimental phenomena were in good agreement with our calculation results.     Figure 8 shows the evolution process of the plasma expansion based on the shadow method using a high-speed camera with a pulse width of 0.2 ms under the same laser energy density (which was 50 µs/frame). Through calculation, we obtained the maximum expansion velocity of the combustion wave in this process, about 18.2 m/s. The observed experimental phenomena were in good agreement with our calculation results.  Figure 8 shows the evolution process of the plasma expansion based on the shadow method using a high-speed camera with a pulse width of 0.2 ms under the same laser energy density (which was 50 μs/frame). Through calculation, we obtained the maximum expansion velocity of the combustion wave in this process, about 18.2 m/s. The observed experimental phenomena were in good agreement with our calculation results.     Figure 9 shows the relationship between the plasma velocity field and the pulse width at three different energy densities of 8.0 × 10 3 J/cm 2 , 9.0 × 10 3 J/cm 2 , and 1.0 × 10 4 J/cm 2 ( Figure 9A-D), respectively, corresponding to the pulse widths of 0.2 ms, 0.5 ms, 0.8 ms, and 1.0 ms.  Figure 10 shows the relationship between the plasma expansion velocity field and the laser pulse width under three different laser energy densities. From Figure 10, it can be seen that when the laser energy density was 8.0 × 10 3 J/cm 2 , the maximum speeds corresponding to 0.   Figure 10 shows the relationship between the plasma expansion velocity field and the laser pulse width under three different laser energy densities. From Figure 10, it can be seen that when the laser energy density was 8.0 × 10 3 J/cm 2 , the maximum speeds corresponding to 0.2 ms, 0.5 ms, 0.8 ms, and 1.0 ms were 11.545 m/s, 10.169 m/s, 8.0608 m/s, and 6.3985 m/s. When the energy density was 9.0 × 10 3 J/cm 2 , the maximum speeds corresponding to the simulation results were 14.76 m/s, 11.848 m/s, 9.8884 m/s, and 8.159 m/s, respectively. When the energy density was 1.0 × 10 4 J/cm 2 , the corresponding maximum speeds were 17.7 m/s, 15.769 m/s, 13.453 m/s, and 10.046 m/s, respectively.
Comparing the plasma flow field velocity in Figures 8 and 9, it can be seen that the time for the plasma flow field velocity to reach the maximum value increased with the increase in the pulse width, and the maximum velocity decreased with the increase in the pulse width. After the laser ended, the plasma flow field velocity rapidly decreased to zero. Velocity (m/s) Pulse width (ms) 8.0×10 3 J/cm 2 (Simulation) 9.0×10 3 J/cm 2 (Simulation) 1.0×10 4 J/cm 2 (Simulation) 8.0×10 3 J/cm 2 (Experiment) 9.0×10 3 J/cm 2 (Experiment) 1.0×10 4 J/cm 2 (Experiment) Figure 10. The relationship between the velocity field and pulse width under different laser energy densities.
Comparing the plasma flow field velocity in Figures 8 and 9, it can be seen that the time for the plasma flow field velocity to reach the maximum value increased with the increase in the pulse width, and the maximum velocity decreased with the increase in the pulse width. After the laser ended, the plasma flow field velocity rapidly decreased to zero.

Conclusions
This paper considered the laser energy transmission processes including thermal radiation, heat conduction, and convection, and simulated the whole physical process of the millisecond pulse laser induced fused silica to generate plasma, and the simulation results have good convergence. The simulation of laser-induced fused silica plasma can easily provide some details, such as the initial plasma temperature, which is not easy to measure in the experiment. It can deepen people's understanding of the initial details of plasma generation and make the plasma generation process less abstract. The plasma simulation model still has some shortcomings, and people need to continue to improve it in future work, so that it can represent a more realistic physical process. The results were as follows: (1) The temperature field mainly discussed the relationship of temperature distribution in the radial position over time. The simulation results showed that the temperature in the center region irradiated by the laser was the highest, and plasma was first generated in this region. (2) In the velocity field, the plasma flow field under different laser energy densities and pulse widths were analyzed and studied, and the relationship between the plasma expansion law and the laser energy density and pulse width were obtained. The simulation results showed that there was a large gap between the axial plasma expansion speed and the radial plasma expansion speed; the plasma expansion speed increased with the increase in the laser energy density and decreased with the increase in the pulse width. The simulation results were in good agreement with the experimental results, and it also provided a theoretical and experimental basis for studying the laser-induced plasma process of fused silica materials.

Conclusions
This paper considered the laser energy transmission processes including thermal radiation, heat conduction, and convection, and simulated the whole physical process of the millisecond pulse laser induced fused silica to generate plasma, and the simulation results have good convergence. The simulation of laser-induced fused silica plasma can easily provide some details, such as the initial plasma temperature, which is not easy to measure in the experiment. It can deepen people's understanding of the initial details of plasma generation and make the plasma generation process less abstract. The plasma simulation model still has some shortcomings, and people need to continue to improve it in future work, so that it can represent a more realistic physical process. The results were as follows: (1) The temperature field mainly discussed the relationship of temperature distribution in the radial position over time. The simulation results showed that the temperature in the center region irradiated by the laser was the highest, and plasma was first generated in this region. (2) In the velocity field, the plasma flow field under different laser energy densities and pulse widths were analyzed and studied, and the relationship between the plasma expansion law and the laser energy density and pulse width were obtained.
The simulation results showed that there was a large gap between the axial plasma expansion speed and the radial plasma expansion speed; the plasma expansion speed increased with the increase in the laser energy density and decreased with the increase in the pulse width. The simulation results were in good agreement with the experimental results, and it also provided a theoretical and experimental basis for studying the laser-induced plasma process of fused silica materials.