|
1.IntroductionHigh-power all-solid-state 2-μm laser systems have played an increasingly important role in the technological landscape during the past few years. Applications have ranged from an abundance of medical uses, such as laser surgery and therapy, to more technical applications, such as laser ranging and remote sensing.1–4 It is important to develop a laser system with high efficiency, long lifetime, and high stability; however, the laser output window may compromise the system’s performance because of the thermal and mechanical effects caused by the window absorption in the process of high-power laser-beam outputting, which not only limits the service life, but also affects the power scaling of the system.5 The thermal-mechanical effects and the damage mechanism of the laser output window have to be studied for further improving the power levels of the system before the onset of window damage.6 These issues have received a great deal of attention, and various models have been proposed including the heating conduction model, the multiphoton ionization model, and cumulative ionization breakdown. Studies have shown7,8 that the multiphoton ionization model and cumulative ionization breakdown dominate only in the case of high-power pulsed output, such situations have difficulties in quantitative analysis, and the heating conduction model is frequently used in the interpretation of continuous wave (cw) laser damage of output windows.9 In this article, we establish a three-dimensional (3-D) thermal model of the laser output window. Expressions for the transient thermal and mechanical distributions are deduced in the case of continuous laser output using the integral-transform method. Then, taking the infrared quartz window material as an example, the temperature and mechanical field distributions of a high-power all-solid-state 2-μm laser system window are simulated. Finally, the laser-induced damage mechanism is analyzed. 2.Model Building2.1.2-μm Laser SystemThe 2-μm Tm:YAG laser system contains a two-rod scheme, in which either rod employs a laser diode (LD) side-pumping laser module. The module contains five diode arrays, and each array consists of 12 diodes operating in cw-mode with a central wavelength of 785 nm. The rod is 4 mm in diameter and 69 mm in length with a -ion doping concentration of 3.5%. The Tm:YAG rod is cooled to 8°C by deionized water. M1 is the flat rear mirror coated to obtain a reflectivity more than 99.5% at 2 μm, and M2 is the output coupler with a transmission of 5% around 2 μm; in this instance, the wavelength of the output laser is 2.07 μm, and the maximum cw output power is 115 W. The schematic of the laser system is shown in Fig. 1 (Ref. 10). The selection of the laser window is extremely important during the design of lasers. Materials such as , , and are commonly used for optical windows,11 and quartz glass, specifically, exhibits outstanding thermodynamic properties, which makes it an attractive selection as the window material of a high-power 2-μm laser system. 2.2.Thermal Model DescriptionHeat deposition arises due to the absorption of the laser-beam energy when a high-energy laser beam passes through the output window, resulting in temperature rise as well as thermal stress. The thermal model of the output window is defined and illustrated in Fig. 2. A laser beam with a radius is incident along the positive -direction for the window M2, which has a radius and a thickness . The 3-D heat transfer equation of the window can be expressed as where , , and indicate the density, heat capacity, and thermal conductivity of the window material, respectively. High-power laser systems generally operate in a multimode state, and thus the laser beam can be assumed to be uniformly distributed. The heat source , therefore, can be described as follows: where is the absorption coefficient of the material. is the peak power density, which can be described as where is the transmittance of the window, and represents the total beam power. The initial and boundary conditions are where is the ambient temperature.3.Theoretical AnalysisThere are many analytical methods for solving Eq. (1) based on the initial and boundary conditions in Eq. (4).12 Numerical methods, including the finite-element method13,14 and the finite-difference method,15 do not involve derivation of formulas or transformations, and so the solution procedure is convenient, especially in simulating multidimensional systems.16 However, numerical solutions cannot provide an intuitive understanding of the relationships among the various physical data, whereas analytical solutions can clearly depict how the temperature-stress field distribution evolves with the parameters of the laser and window material.17,18 Therefore, we employ the approach of using the integral-transform method to obtain such analytical solutions in this article. Based on the model established above, we obtain the following formulas through utilizing the integral-transform method to solve Eqs. (1)–(4) and performing the positive and inverse integral transforms for variables and :19 Here, is the temperature distribution of the window, and Eq. (6) is the inverse transform of Eq. (5). The eigenfunctions and the norm are obtained as where the eigenvalues are the positive roots of , refers to the ’th Bessel function, is the radial coordinate, and is the window radius.The eigenfunctions and the norm are obtained as where the eigenvalues are the roots of .The 3-D temperature distribution of the output window is given by Eq. (11) according to the inferential method from Eqs. (5)–(10): The nonuniform temperature distribution experienced by the output window causes mechanical stresses, which can be obtained from the stress-strain equation and Hooke equation20 where , , and are the thermal stresses in the radius, hoop, and axial directions, respectively; is the thermal coefficient of expansion; is the Young’s modulus; and is the Poisson’s ratio. Introducing the temperature distribution of Eq. (11) into Eqs. (12)–(14), we get the stress distributions of the output window4.Numerical Computation and DiscussionsIn this section, some numerical calculations concerning the thermal–mechanical effects are performed, and the damage mechanism of the laser output window is investigated in the case of taking infrared quartz glass as the output window. Infrared quartz glass is used as the 2-μm laser optical window material owing to its fine performance, whose application wavelength is in the range of 260 to 3500 μm and the infrared transmittance is 85%. Some parameters of the glass in the calculation are listed in Table 1. We assume that the diameter and the thickness of the window are, respectively, 25 and 6 mm, the laser beam radius is 2 mm, and the ambient temperature is 20°C. Table 1Physical parameters of infrared quartz glass.
4.1.Analysis of Temperature and Stress FieldsFigure 3 illustrates the 3-D temperature distribution of the window in the case of a laser system producing continuous output for 8 s when the laser output power density is . Figure 4 displays the temperature distributions of the window surface () for different laser-beam radii and different durations when the laser operates at a maximum output power of 115 W. As can be seen in Fig. 3, in the radius direction, the nearer the sample is to the center of the window, the higher the temperature is, whereas in the axial direction, the temperature does not change significantly. At the edge of the laser beam, the temperature gradient increases sharply, reflecting the thermal expansion of the actively heated area against the passively heated area of the window. Figure 4 shows that the temperature rises higher with longer run times and larger laser-beam radii. It can also be seen that the maximum temperature rise is 289°C when the laser operates at maximum output power. Figures 5(a), 5(b), and 5(c) depict the stress distributions of the window in the radial, hoop, and axial directions, respectively. Figures 6(a), 6(b), and 6(c) illustrate the situation for different run times and different laser-beam radii, all of which depict how the radial, hoop, and axial stresses evolve at the window surface () when the laser operates at the maximum output power. From Fig. 5, it can be seen that all stresses are compressive on the centerline of the window, where the maximum heat deposition occurs. Radial stress is always compressive; the hoop and axial stresses, however, turn tensile at the outer rim of the window. Figure 6 illustrates that the stress gradients become greater with longer run times or larger laser-beam radii. 4.2.Damage Mechanism AnalysisThe infrared quartz window may be damaged because of either the temperature reaching the melting point (1780°C) or the tensile (compressive) stress exceeding the tensile (compressive) strength of 48 MPa (1100 MPa). From the analysis mentioned above, it can be seen that the temperature of window would exceed the melting point while the thermal stresses would not reach the limits, providing that the laser power density reaches and the continuous laser output lasts for 8 s. The thermal stress on the quartz material caused by the temperature gradient are relatively small due to the good thermal stability of the quartz; thus the window is primarily damaged by melting in the area where the temperature exceeds the melting point of the material, so that melting, ablation, or even perforation arises as a result. 5.Laser OutputThe wavelength of the 2-μm laser system can switch depending on the transmission of the output window. When the coupler M2 has a transmission of 10%, the wavelength of the output laser is 2.02 μm, and the maximum cw output power is 77.1 W. When the coupler has a transmission of 5%, the wavelength of the laser output is shifted to 2.07 μm, and the maximum cw output power is 115 W. The laser spectra at the maximum output power were recorded by an optical spectrometer (NIRQuest256-2.5), as shown in Fig. 7. Damage to the output window can be avoided by controlling the output power of the laser system. The dependence of the laser output power on the LD pump power is obtained with a power meter (Ophir F300A-SH), as illustrated in Fig. 8. With a coupler transmission of 5%, the maximum output power is 115 W under an LD pump power of 1188 W, while the maximum power of 77.1 W is achieved with a coupler transmission of 10% under an LD pump power of 1188 W. 6.ConclusionsIn summary, we have defined a thermal model of a laser window and deduced expressions for the transient thermal and mechanical distributions. Taking the infrared quartz glass as an example, the temperature and mechanical field distributions of a high-power all-solid-state 2-μm laser system window were simulated, and the laser-induced damage mechanism was analyzed. The results showed that the window damage mechanism is mainly due to the melting in the areas where the temperature exceeds the melting point of the material, when a 2-μm high-power laser continuously outputs the optical power. We believe that the research methods and results presented in this article will be helpful for the selection of laser windows. AcknowledgmentsThis work was supported by the National Natural Science Foundation under Grant No. 61074158 (China). ReferencesQ. L. Maet al.,
“Light scattering and 2 μm laser performance of Tm:YAG ceramic,”
Opt. Commun., 284
(6), 1645
–1674
(2011). http://dx.doi.org/10.1016/j.optcom.2010.11.051 OPCOB8 0030-4018 Google Scholar
B. Q. Yaoet al.,
“Continuous-wave operation of a room-temperature Tm:YAG-pumped Ho:YAG laser,”
Chin. Opt. Lett., 6
(7), 520
–522
(2008). Google Scholar
J. Yuet al.,
“ Q-switched 2 μm solid-state laser,”
Opt. Lett., 31
(4), 462
–464
(2006). http://dx.doi.org/10.1364/OL.31.000462 OPLEDP 0146-9592 Google Scholar
L. Wanget al.,
“Resonantly pumped monolithic nonplanar Ho:YAG ring laser with high-power single-frequency laser output at 2122 nm,”
Opt. Express, 21
(8), 9541
–9546
(2013). http://dx.doi.org/10.1364/OE.21.009541 OPEGAR 0091-3286 Google Scholar
H. J. Liuet al.,
“Experimental comparison of damage performance induced by nanosecond laser between K9 and fused silica optics,”
Acta Phys. Sin., 61
(7), 076103
(2012). http://dx.doi.org/10.7498/aps.61.076103 WLHPAR 1000-3290 Google Scholar
M. E. Innocenziet al.,
“Thermal modeling of continuous-wave end-pumped solid-state lasers,”
Appl. Phys. Lett., 56
(19), 1831
–1833
(1990). http://dx.doi.org/10.1063/1.103083 APPLAB 0003-6951 Google Scholar
Y. X. Panet al.,
“Effect of inclusion matrix model on temperature and thermal stress fields of K9-glass damaged by long-pulse laser,”
Opt. Eng., 52
(4), 044302
(2013). http://dx.doi.org/10.1117/1.OE.52.4.044302 OPEGAR 0091-3286 Google Scholar
A. S. EpifanovA. A. ManenkovA. M. Prokhorov,
“Theory of avalanche ionization induced in transparent dielectrics by an electromagnetic field,”
Sov. Phys. JETP, 43
(2), 377
–381
(1976). SPHJAR 0038-5646 Google Scholar
C. L. Wanget al.,
“Wavelength switchable high-power diode-side-pumped rod Tm: YAG laser around 2 μm,”
Opt. Express, 21
(6), 7156
–7161
(2013). http://dx.doi.org/10.1364/OE.21.007156 OPEGAR 0091-3286 Google Scholar
F. HuangY. X. NiuY.F. Wang,
“Calculation of thermal and mechanical effect induced by laser in optical window materials,”
Acta Opt. Sin., 26
(4), 576
–580
(2006). GUXUDC 0253-2239 Google Scholar
C. A. Klein,
“High-energy laser windows: case of fused silica,”
Opt. Eng., 49
(9), 091006
(2010). http://dx.doi.org/10.1117/1.3484946 OPEGAR 0091-3286 Google Scholar
Y. F. Penget al.,
“Laser-induced temperature distributions and thermal deformations in sapphire, silicon, and calcium fluoride substrates at 1.315 μm,”
Opt. Eng., 40
(12), 2822
–2829
(2001). http://dx.doi.org/10.1117/1.1418712 OPEGAR 0091-3286 Google Scholar
Y. X. Niuet al.,
“Thermal shock effect on diamond-like carbon thin films induced by pulsed-laser,”
Acta Phys. Sin., 54
(10), 4816
–4821
(2005). WLHPAR 1000-3290 Google Scholar
X. F. DuanY. X. NiuC. Zhang,
“Calculation of laser irradiation effect and analysis of laser-induced damage threshold in semiconductor,”
Acta Opt. Sin., 24
(8), 1057
–1061
(2004). GUXUDC 0253-2239 Google Scholar
P. H. LuR. W. Wang,
“Three-dimensional temperature distribution analysis and thermal lens effect calculation for high power laser windows,”
Acta Opt. Sin., 21
(8), 965
–969
(2001). GUXUDC 0253-2239 Google Scholar
Y. S. FengX. X. Li,
“ANSYS numerical simulation of temperature field in silica glass irradiated by pulsed laser,”
Phys. Exp., 32
(2), 35
–37
(2012). Google Scholar
J. Liet al.,
“Numerical analysis of CW laser damage in GaAs,”
Acta Photon. Sin., 41
(5), 571
–574
(2012). http://dx.doi.org/10.3788/gzxb 1004-4213 Google Scholar
J. Wuet al.,
“Thermal effect in pulsed laser diode dual-end pumped Tm:YAG laser,”
Proc. SPIE, 8904 890413
(2013). http://dx.doi.org/10.1117/12.2034333 PSISDG 0277-786X Google Scholar
M. N. Ozisik, Heat Conduction, High Education Press, Beijing
(1983). Google Scholar
P. Xiuer, Thermal Stress and Thermal Fatigue, 12
–30 National Defense Industry Press, Beijing
(1984). Google Scholar
BiographyWenwen Liu received her BS degree in optical engineering from Beihang University, Beijing, China, in 2010. She is currently pursuing her MS degree in optical engineering at Beihang University. Her research is focused on laser technology. |