Estimation of thermal fracture limits in quasi-continuous-wave end-pumped lasers through a time-dependent analytical model

doi:10.1364/OE.16.011115

« Show journal navigation

Estimation of thermal fracture limits in quasi-continuous-wave end-pumped lasers through a time-dependent analytical model

E. H. Bernhardi, A. Forbes, C. Bollig, and M.J. D. Esser »View Author Affiliations

Optics Express, Vol. 16, Issue 15, pp. 11115-11123 (2008)
http://dx.doi.org/10.1364/OE.16.011115

View Full Text Article

Acrobat PDF (428 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Full Text
Article Info
References (17)
Cited By
Figures (5)
Multimedia (2)
Metrics

Abstract

A time-dependent analytical thermal model of the temperature and the corresponding induced thermal stresses on the pump face of quasi-continuous wave (qcw) end-pumped laser rods is derived. We apply the model to qcw diode-end-pumped rods and show the maximum peak pump power that can be utilized without fracturing the rod. To illustrate an application of the model, it is applied to a qcw pumped Tm:YLF rod and found to be in very good agreement with published experimental results.The results indicate new criteria to avoid fracture when operating Tm:YLF rods at low qcw pump duty cycles.

1. Introduction

The power scaling of diode-end-pumped solid-state lasers is a very active area of research. The main problem that limits the power scaling of these lasers is the generation of heat inside the laser gain medium; the generated heat causes steep temperature gradients inside the crystal, which in turn produce stress, leading to fracture. Fracture of the laser material occurs when the thermally induced stress exceeds the ultimate strength of the material [1

W. Koechner Solid-State Laser Engineering 4th ed., (Springer-Verlag Berlin, Heidelberg, Germany, 1996).

]. When the continuous wave (cw) pump power exceeds the power at which crystal fracture occurs, the pump source is often modulated in time, creating a so called quasi-continuous wave (qcw) pump, with the effect of reducing the average pump power to below the fracture limit, while maintaining a high output power during the on-time of the pump pulse.

In order to investigate the thermally induced stresses and the power limitations due to fracture, a thermal model of the laser gain medium is required. Existing analytical thermal models that describe the temperature and stresses in laser crystals are restricted to special cases and approximations, such as cw pump sources and steady-state conditions [1–5

W. Koechner Solid-State Laser Engineering 4th ed., (Springer-Verlag Berlin, Heidelberg, Germany, 1996).

]. In this paper an analytical thermal model that determines the transient behaviour of the temperature and the corresponding induced stresses on the pump face of an isotropic laser rod is derived from first principles. To the best of our knowledge this is the first time that such a time-dependent analytical model has been reported. We validate the model through finite element analysis, and apply the model to qcw pumped Tm:YLF laser rods, and find favourable agreement between the calculated fracture limits and the experimentally determined values reported in the literature.

We have particularly chosen to apply the model to a YLF rod because power scaling with this material is limited by the relatively low fracture limit of 40 MPa, which is ∼5 times lower than that of YAG [6

X. Peng, L. Xu, and A. Asundi, “High-power efficient continous-wave TEM00 intracavity frequencydoubled diode-pumped Nd:YLF laser,” Appl. Opt. 44, 800–807 (2005). [CrossRef] [PubMed]

]. Despite being derived for isotropic rods, we show that the model may be applied to anisotropic rods on condition that the highest linear expansion and the lowest thermal conductivity of the respective a- and c-axis of the crystal are used in the calculation. Finally, we illustrate how the model may be used to estimate the peak pump power that can safely be used to qcw pump a laser rod at a given duty cycle, opening the way to fracture-free power scaling with qcw pump sources.

2. Theory

In this section the transient temperature and stress profiles on the pump face of a longitudinally pumped isotropic laser rod are derived. We assume that the Rayleigh range of the pump beam is much longer than the length of the rod so that the pump beam is treated as perfectly collimated inside the rod. This implies that only the pump face need be considered since the pump light in the gain medium follows an exponential decay, with the steepest temperature gradient (and thus stress) on the pump face [2

C. Pfistner, R. Weber, H. P. Weber, S Merazzi, and R. Gruber, “Thermal Beam Distortions in End-Pumped Nd:YAG, Nd:GSGG and Nd:YLF,” IEEE J. Quantum Electron. 30, 1605–1615 (1994). [CrossRef]

L. Yan and C. H. Lee, “Thermal Effects in End-Pumped Nd:phosphate Glasses,” J. Appl. Phys. 75, 1286–1292 (1994). [CrossRef]

]. A closed form solution for the time dependent temperature profile on the surface of the crystal rod, u(r,t), may be found by solving the non-homogeneous heat diffusion equation [7

G. Barton, Elements of Green’s Functions and Propagation , (Oxford University Press, Oxford,1995).

H. S. Carslaw and J. C. Jaeger, Conduction of Heats in Solids , (Oxford University Press, Oxford, 1959).

\frac{\partial u (r, t)}{\partial t} - D \nabla^{2} u (r, t) = Q (r, t),

(1)

for a generalized source term Q(r,t). Here D=k/ρC_p is the diffusivity, k is the thermal conductivity, ρ is the density and C_p is the heat capacity of the laser material respectively,while all other terms have their usual meaning. We will assume that the rod has a length l and a radius R, and is homogenous and isotropic. Furthermore, we will assume that the boundary of the rod is at a constant temperature, u(R,t)=0, with no initial temperature profile on its pump face: u(r,0)=0. With these boundary conditions, Eq. (1) may be solved directly by use of an appropriate Green’s function, with the solution given in integral form as:

u (r, t) = \int_{0}^{t} \int_{0}^{R} Q (ξ, τ) G (r, ξ, t - τ) dξdτ,

(2)

with the Green’s function G(r,ξ,t) given by [7

G. Barton, Elements of Green’s Functions and Propagation , (Oxford University Press, Oxford,1995).

G (r, ξ, t) = \sum_{m = 1}^{\infty} \frac{2 ξ}{R^{2} J_{1}^{2} (μ_{m})} J_{0} (\frac{μ_{m} r}{R}) J_{0} (\frac{μ_{m} r}{R}) \exp (- \frac{D μ_{m}^{2} t}{R^{2}}) .

(3)

The summation is over the positive roots (μ_m ) of the zeroth order Bessel function, J ₀. For the common experimental configuration where the pump beam is imaged onto the rod face from a multimode fibre, the pump intensity may be approximated as having a top-hat spatial intensity profile. Figure 1 shows an example of a measured top-hat transverse intensity profile as produced by a fibre-coupled diode pump. If in addition the pump is not assumed to be cw but rather the more general qcw, consisting of a pulse train of on-off pulses with on-time τ_on (with τ_on usually longer than the upper-state lifetime of the laser gain medium) and period T,we may write the source term as:

Q (r, t) = {\begin{matrix} \frac{αηE}{π w^{2} ρ C_{p} τ_{on}} \\ 0 \end{matrix} \begin{matrix} ; nT \leq t \leq nT + τ_{on} \\ ; nT + τ_{on} \leq t \leq (n + 1) T \end{matrix},

(4)

where n is the number of pulses, E is the energy in each pulse and w is the radius of the top-hat beam. α is the absorption coefficient in units of inverse length of the crystal, while η is a parameter to account for the fact that not all the absorbed pump light is converted into heat.This model can be applied under lasing or non-lasing conditions by choosing an appropriate value for η. Typical values for the heat load efficiency are η=0.32 (lasing) and η=0.4 (non-lasing) [9

T. Y. Fan, “Heat Generation in Nd:YAG and Yb:YAG,” IEEE J. Quantum Electron. 29, 1457–1459 (1993). [CrossRef]

Fig. 1. An example of a measured top-hat transverse intensity profile produced by a fibre coupled diode laser pump (own experimental results).

Substituting Eqs. (3) and (4) into Eq. (2

C. Pfistner, R. Weber, H. P. Weber, S Merazzi, and R. Gruber, “Thermal Beam Distortions in End-Pumped Nd:YAG, Nd:GSGG and Nd:YLF,” IEEE J. Quantum Electron. 30, 1605–1615 (1994). [CrossRef]

) and solving both the time and spatial integrals separately yields the following analytical expression for the temperature on the pump face of the rod:

u (r, pT + t) = \frac{2 αηER}{kπw τ_{on}} \sum_{m = 1}^{\infty} \frac{J_{0} (\frac{μ_{m} r}{R}) J_{1} (\frac{μ_{m} w}{R}) f (p, t, μ_{m})}{μ_{m}^{3} J_{1}^{2} (μ_{m})},

(5)

with the time dependence given by f(p,t,μ_m ) as:

f (p, t, μ_{m}) = \exp (- μ_{m}^{2} \frac{t}{τ_{D}}) {\frac{[\exp (μ_{m}^{2} \frac{τ_{on}}{τ_{D}}) - 1] [\exp (- μ_{m}^{2} \frac{pT}{τ_{D}}) - 1]}{1 - \exp (μ_{m}^{2} \frac{T}{τ_{D}})} - [1 - \exp (μ_{m}^{2} \frac{τ}{τ_{D}})]},

(6)

with τ=min[t,τ_on ]. Here we have introduced a new variable, τ_D =R ²/D, which we refer to as the diffusion relaxation time of the system, and for convenience the temperature is calculated after p complete pulses plus some time t into the p+1 pulse, so that the total elapsed time from the start of the pumping process is pT+t. While the summation in Eq. (5) is carried to infinity, in practice one finds that 30 terms or above leads to very good convergence of the series. Equation (5) may easily be rewritten in terms of peak pump power (P_p ) or average pump power (P_av ) rather than pump energy (E) by noting that for the source term in this study P_p =E/τ_on and P_av =E/T (note: here we have assumed that the time pulse envelope is a square pulse; for other cases the peak pump power expression may have to be modified in an appropriate manner). In the special case of a cw pump source, T=τ_on so that the peak and average pump powers are identical.

By making use of a plane-strain approximation and assuming that the stress in the axial direction is zero, we may calculate the radial and tangential stresses from the temperature profile from [10

S. P. Timoshenko and J. N. Goodier, Theory of Elasticity 3rd ed., (McGraw-Hill, New York, 1970).

,11

B. A. Boley and J. H. Weiner, Theory of Thermal Stresses , (Courier Dover Publications, New York, 1997).

σ_{r} (r, t) = c [\frac{1}{R^{2}} \int_{0}^{R} u (r, t) rdr - \frac{1}{r^{2}} \int_{0}^{r} u (r, t) rdr];

(7a)

σ_{θ} (r, t) = C [\frac{1}{R^{2}} \int_{0}^{R} u (r, t) rdr + \frac{1}{r^{2}} \int_{0}^{r} u (r, t) rdr - u (r, t)],

(7b)

where C=γY/(1- V ), with γ the linear coefficient of expansion, Y is Young’s modulus and V is Poisson’s ratio. The plain-strain approximation is valid for l/R ≫ 1 (a long rod) [11

B. A. Boley and J. H. Weiner, Theory of Thermal Stresses , (Courier Dover Publications, New York, 1997).

]. For the case where l/R≪1 (a thin disk), the plane-stress approximation is used where Eqs. (7) and (8) also hold with C=γY [11

B. A. Boley and J. H. Weiner, Theory of Thermal Stresses , (Courier Dover Publications, New York, 1997).

]. One can readily show that Eqs. (7) and (8) can be solved analytically to yield:

σ_{r} (r, pT + t) = \frac{2 CαηER}{kπw τ_{on}} \sum_{m = 1}^{\infty} \frac{J_{1} (\frac{μ_{m} w}{R})}{μ_{m}^{3} J_{1}^{2} (μ_{m})} [\frac{J_{1} (μ_{m})}{μ_{m}} - \frac{R J_{1} (\frac{μ_{m} r}{R})}{{rμ}_{m}}] f (p, t, μ_{m});

(8a)

σ_{θ} (r, pT + t) = \frac{2 CαηER}{kπw τ_{on}} \sum_{m = 1}^{\infty} \frac{J_{1} (\frac{μ_{m} w}{R})}{μ_{m}^{3} J_{1}^{2} (μ_{m})} [\frac{J_{1} (μ_{m})}{μ_{m}} + \frac{R J_{1} (\frac{μ_{m} r}{R})}{{rμ}_{m}} - J_{0} (\frac{μ_{m} r}{R})] f (p, t, μ_{m}) .

(8b)

Due to the fact that the stress tensor alone does not provide enough information regarding crystal fracture, we use the maximum shear stress to predict fracture [12

LASCAD Manual, (LASCAD, 2008), https://www.las-cad.com/files/PP_FEA.pdf.

]. This is also known as the stress intensity or the Tresca failure criterion, which in the plain-strain approximation reduces to:

σ_{T} (r, pT + t) = ∣ σ_{θ} - σ_{r} ∣

= ∣ \frac{2 CαηER}{kπw τ_{on}} \sum_{m = 1}^{\infty} \frac{J_{1} (\frac{μ_{m} w}{R})}{μ_{m}^{3} J_{1}^{2} (μ_{m})} J_{2} (\frac{μ_{m} r}{R}) f (p, t, μ_{m}) ∣,

(9)

where we have made use of the well known relation J ₂(x)=2J ₁(x)/x-J ₀(x). Equations (5), (8)and (9) allow the temperature, stresses and fracture limit to be calculated for any qcw pulse train, as a function of both time and position in the crystal rod.

3. Model validation

As a verification of the analytical thermal model, a time-dependent three dimensional coupled thermal-stress finite element analysis was implemented in a commercial software package,ABAQUS [13

ABAQUS Finite Element Software Package (ABAQUS, 2008) http://www.simulia.com/products/abaqus_fea.html.

,14

E. H. Bernhardi, C. Bollig, L. Harris, M. J. D. Esser, and A. Forbes, “Investigating thermal stresses in quasicw pumped Tm:YLF laser crystals,” in Proceedings of Advanced Solid-State Photonics , (Nara, Japan,2008), Poster WB11.

]. The average computation time of each finite element simulation (5560 discretization) was approximately 5.5 hours on a dual-core 2.4 GHz processor (the analytical model takes ∼1-2 s). In analogy with convective heat transfer, we have specified a Neumann boundary condition with a heat transfer coefficient at the interface between the crystal and the copper heat sink. We have used a heat transfer coefficient of 0.9 Wcm^-2K^-1,which corresponds to a layer of indium foil between the crystal and the heat sink [5

S. Chénais, S. Forget, F. Druon, F. Balembois, and P. Georges, “Direct and Absolute Temperature Mapping and Heat Transfer Measurements in Diode-End-Pumped Yb:YAG,” Appl. Phys. B 79, 221–224 (2004). [CrossRef]

]. In the remaining text the finite element analysis will be referred to as the numerical model [14

The analytical and numerical models have been applied to reported results [15

S. So, J. I. Mackenzie, D. P. Shepherd, W. A. Clarkson, J. G. Betterson, and E. K. Gorton, “A power-scaling strategy for longitudinally diode-pumped Tm:YLF lasers,” Appl. Phys. B 84, 389–393 (2006). [CrossRef]

] of a cw pumped, 4% doped, Tm:YLF laser rod with l=12 mm and R=1.5 mm (l/R=8 so that we make use the plain-strain approximation as discussed in section (2)). The pump beam was a near top-hat profile with M ²∼100 and a pump radius of w=470 μm in the middle of the crystal. The onset of fracture was reported at an incident power of 47.2 W (α=1.43 cm^-1;η=0.33). While the actual stress values at fracture are not reported, it is known that the fracture limit of YLF crystals is in the 33-40 MPa range [1

W. Koechner Solid-State Laser Engineering 4th ed., (Springer-Verlag Berlin, Heidelberg, Germany, 1996).

X. Peng, L. Xu, and A. Asundi, “High-power efficient continous-wave TEM00 intracavity frequencydoubled diode-pumped Nd:YLF laser,” Appl. Opt. 44, 800–807 (2005). [CrossRef] [PubMed]

]. The parameters that were used in the numerical and analytical thermal models are shown in Table 1.

The analytical model predicts a maximum Tresca stress of 42 MPa which agrees very well with the reported 33-40 MPa fracture limit range [1

W. Koechner Solid-State Laser Engineering 4th ed., (Springer-Verlag Berlin, Heidelberg, Germany, 1996).

X. Peng, L. Xu, and A. Asundi, “High-power efficient continous-wave TEM00 intracavity frequencydoubled diode-pumped Nd:YLF laser,” Appl. Opt. 44, 800–807 (2005). [CrossRef] [PubMed]

]. The numerical model results in a maximum Tresca stress of 41 MPa. This shows that for a cw pump, the analytical thermal model is consistent with both the experimental fracture data and with the numerical solutions of a three-dimensional finite element analysis. Since the analytical thermal model proved to be accurate for a cw pump beam, the transient behaviour of the temperature and the induced stresses were determined for various qcw pump duty cycles, defined as τ_on /T. Figure 2(a) and 2(b) show the predicted time-dependence of the temperature in the centre of the pump face of the Tm:YLF rod for a 10% and a 50% pump duty cycle respectively (τ_on =10 ms). The upper and lower boundaries of the shaded red region in Fig. 2 indicate the analytical model’s predictions of the temperature when the thermal conductivity of the c- and the a-axis of Tm:YLF were used respectively. It is clear from the graphs that there is very good agreement between the analytical and numerical models when the lowest thermal conductivity is used in the calculations. Figure 3 illustrates that the thermally induced stresses as calculated by the two models show very good agreement when the highest linear expansion coefficient is used in the analytical model. The upper and lower boundaries of the shaded red region in Fig. 3 indicate the analytical model’s predictions of the maximum stress on the pump face when the two respective linear expansion coefficients of Tm:YLF were used along with the lowest thermal conductivity.

By considering the transient stress distribution on the entire pump face with a qcw pump source, it is evident that the maximum stress does not always occur on the edge-surface (r=R) of the rod; this is contrary to the case of a cw pump beam [4

L. Yan and C. H. Lee, “Thermal Effects in End-Pumped Nd:phosphate Glasses,” J. Appl. Phys. 75, 1286–1292 (1994). [CrossRef]

]. The position of the maximum stress changes during a single qcw pump pulse as well as between qcw pulses.Figure 4(a) shows an animation of the analytically predicted σ_T on the pump face of the Tm:YLF rod when subjected to a 90 W peak power pump beam at 50 Hz (T=20 ms; τ_on =10 ms), while Fig. 4(b) shows an animation of the numerically predicted σ_T throughout the bulk of the Tm:YLF rod. A cross-section of the stress is overlaid on the analytical animation, and clearly shows the “peak” where the stress is maximum and its movement in time as the pulses accumulate. Note that the numerical solution does not exhibit the same symmetrical stress distribution on the pump face as the analytical solution since the anisotropic characteristics of Tm:YLF were accounted for in the numerical model and not in the analytical model.

Table 1. Parameter values of the pumped Tm:YLF rod that were implemented in the simulations.

Parameter	Thermal Model	Reference
Pump beam radius (w) [mm]	0.47	[15 S. So, J. I. Mackenzie, D. P. Shepherd, W. A. Clarkson, J. G. Betterson, and E. K. Gorton, “A power-scaling strategy for longitudinally diode-pumped Tm:YLF lasers,” Appl. Phys. B 84, 389–393 (2006). [CrossRef] ]
Rod radius (R) [mm]	1.5	[15 S. So, J. I. Mackenzie, D. P. Shepherd, W. A. Clarkson, J. G. Betterson, and E. K. Gorton, “A power-scaling strategy for longitudinally diode-pumped Tm:YLF lasers,” Appl. Phys. B 84, 389–393 (2006). [CrossRef] ]
Absorption coefficient (α) [cm^-1]	1.43	[15 S. So, J. I. Mackenzie, D. P. Shepherd, W. A. Clarkson, J. G. Betterson, and E. K. Gorton, “A power-scaling strategy for longitudinally diode-pumped Tm:YLF lasers,” Appl. Phys. B 84, 389–393 (2006). [CrossRef] ]
Thermal conductivity (k) [W.m^-1.K^-1]	7.2 (a-axis), 5.8 (c-axis)	[1 W. Koechner Solid-State Laser Engineering 4th ed., (Springer-Verlag Berlin, Heidelberg, Germany, 1996). ,16 M. Pollnau, P. J. Hardman, M. A. Kern, W. A. Clarkson, and D. C. Hanna, “Upconversion-induced heat generation and thermal lensing in Nd:YLF and Nd:YAG,” Phys. Rev. B 58, 16076–16092 (1998). [CrossRef] ]
Linear expansion coefficient (γ) [10^-6 K^-1]	13 (a-axis), 8.0 (c-axis)	[1 W. Koechner Solid-State Laser Engineering 4th ed., (Springer-Verlag Berlin, Heidelberg, Germany, 1996). ,17 B. M. Walsh, N. P. Barnes, M. Petros, J. Yu, and U. N. Singh, “Spectroscopy and modeling of solid state lanthanide lasers: Application to trivalent Tm³⁺ and Ho³⁺ in YLiF₄ and LuLiF₄ ,” J. Appl. Phys. 95, 3255–3271 (2004). [CrossRef] ]
Fractional heat load (η)	0.33	estimated
Poisson’s ratio ( V )	0.33	[1 W. Koechner Solid-State Laser Engineering 4th ed., (Springer-Verlag Berlin, Heidelberg, Germany, 1996). ,16 M. Pollnau, P. J. Hardman, M. A. Kern, W. A. Clarkson, and D. C. Hanna, “Upconversion-induced heat generation and thermal lensing in Nd:YLF and Nd:YAG,” Phys. Rev. B 58, 16076–16092 (1998). [CrossRef] ]
Young’s modulus (Y) [GPa]	75	[1 W. Koechner Solid-State Laser Engineering 4th ed., (Springer-Verlag Berlin, Heidelberg, Germany, 1996). ,2 C. Pfistner, R. Weber, H. P. Weber, S Merazzi, and R. Gruber, “Thermal Beam Distortions in End-Pumped Nd:YAG, Nd:GSGG and Nd:YLF,” IEEE J. Quantum Electron. 30, 1605–1615 (1994). [CrossRef] ]
Density (ρ) [g.cm^-3]	3.9	[1 W. Koechner Solid-State Laser Engineering 4th ed., (Springer-Verlag Berlin, Heidelberg, Germany, 1996). ,17 B. M. Walsh, N. P. Barnes, M. Petros, J. Yu, and U. N. Singh, “Spectroscopy and modeling of solid state lanthanide lasers: Application to trivalent Tm³⁺ and Ho³⁺ in YLiF₄ and LuLiF₄ ,” J. Appl. Phys. 95, 3255–3271 (2004). [CrossRef] ]
Specific heat capacity (C_p ) [J.g^-1.K^-1]	0.79	[1 W. Koechner Solid-State Laser Engineering 4th ed., (Springer-Verlag Berlin, Heidelberg, Germany, 1996). ,17 B. M. Walsh, N. P. Barnes, M. Petros, J. Yu, and U. N. Singh, “Spectroscopy and modeling of solid state lanthanide lasers: Application to trivalent Tm³⁺ and Ho³⁺ in YLiF₄ and LuLiF₄ ,” J. Appl. Phys. 95, 3255–3271 (2004). [CrossRef] ]

Fig. 2. The analytically (red) and numerically (black) predicted temperature in the centre of the Tm:YLF rod as a function of time while the rod is subjected to a qcw pump with a peak power of (a) 200 W at 10 Hz (τ_on =10 ms) and, (b) 90 W at 50 Hz (τ_on =10 ms).

Fig. 3. The maximum stress on the pump face of the Tm:YLF rod as a function of time while the rod is subjected to a qcw pump with a peak power of (a) 200 W at 10 Hz (τ_on = 10 ms), and (b) 90 W at 50 Hz (τ_on = 10 ms). The analytical (red) and numerical (black) solutions are shown.

Fig. 4. (0.75 MB and 0.33 MB respectively) Animations of (a) the analytical stress distribution on the pump face [Media 1] and (b) the numerical stress distribution in volume of the Tm:YLF rod while it is subjected to a 90 W peak power qcw pump beam at 50 Hz (τ_on = 10 ms) [Media 2].

4. Power scaling of Tm:YLF rods by qcw pumping

To achieve ever higher output powers from a diode-end-pumped solid state laser, the pump power itself must be increased in a concomitant manner. When the pump is a cw source, the high average powers required in some applications leads to high thermal loads, and the onset of fracture. A standard solution to this problem is to employ a qcw pump source, with the advantages that: (i) the average thermal load is reduced through a reduced duty cycle, and (ii) the qcw pulsing leads to higher peak pump power, resulting in much higher laser output power during the qcw pump pulses.

During the qcw operation of a solid-state laser, the generally accepted criteria to avoid thermal fracture is to pump the crystal with an average power (P_av ) that is below the cw fracture pump power (P_cw ):

P_{av \leq P_{cw}} .

(10)

The thermal model developed in section (2) provides a more comprehensive criterion for safe qcw operation, through direct application of Eq. (9), and the results are shown in Fig. 5. The green shaded area in Fig. 5(a) shows the analytically predicted average power with which the Tm:YLF crystal can be pumped without causing thermal fracture, while the red region indicate the predicted onset of fracture using Eq. (10). The yellow region then indicates the region of disagreement between the two models. It is clear that for qcw pump duty cycles of 40% and higher, fracture of the Tm:YLF rod will indeed occur at average pump powers which are equal or greater to the cw fracture pump power, so that Eqs. (9) and (10) are in agreement.At pump duty cycles that are lower than 40%, the analytical model predicts that crystal fracture will occur at average powers that are significantly lower than the cw fracture pump power. By way of example, consider a qcw pump duty cycle of 10% (T=100 ms; τ_on =10 ms) where the fracture limit of Tm:YLF is at P_cw =47 W. For this qcw duty cycle, fracture would occur at P_av >26 W, with the region 47 W >P_av >26 W indicating the error of using Eq. (10). This error is indicated for all duty cycles as the yellow shaded area in Fig. 5(a). The implication is that Eq. (10) is a necessary but not sufficient condition for fracture-free qcw operation in Tm:YLF lasers.

An alternative approach to illustrate the results of the analytical thermal model is to consider the peak power during a qcw pump pulse and to note the peak power at which fracture occurs. The green shaded area in Fig. 5(b) indicates the qcw peak power that can be used to pump the Tm:YLF rod without fracturing it as predicted by the analytical model,while the red region indicates the fracture limit using Eq. (10); the yellow region is once again the region of discrepancy.

Fig. 5. (a). The average pump power (as a fraction of the cw fracture power P_cw ) at which fracture of the Tm:YLF rod occurs as a function of qcw pump duty cycle (τ_on =10 ms). The green shaded region indicates the average pump power at which the Tm:YLF rod can be pumped without fracturing according to the analytical model. The yellow shaded region indicates the difference between the analytical model and P_cw . (b) The same notation as in (a) but for the peak pump power (in units of P_cw ) at which fracture of the Tm:YLF rod occurs as a function of qcw pump duty cycle.

Considering the qcw pump duty cycle of 10%, we note that the Tm:YLF rod can be pumped with a peak pump power of 5.5× higher (270 W) than the cw fracture power before fracture will occur, while Eq. (10) would predict a 10× higher value (470 W).

The model predictions can thus be summarised as follows: the rule of thumb given by Eq. (10) is a necessary but not sufficient condition for fracture-free power scaling of Tm:YLF through the use of qcw pump sources. This is not surprising given that this criteria is not derived from an analysis of the onset of fracture, but is based rather on intuition. Equation (9) is derived from a full thermal analysis, and we suggest that this provides the necessary criterion for any laser rod to be pumped without fracture. In the limit that the duty cycle approaches the cw case, Eq. (9) correctly converges to Eq. (10) as expected.

5. Conclusion

A time-dependent analytical thermal model was developed to investigate the transient behaviour of thermally induced stresses in qcw end-pumped laser rods. The versatility of such a model is that all the material and pump laser parameters may be varied analytically, thereby aiding physical insight. For example, the model confirms exactly that the temperature increases linearly with pump power, as expected. Similarly, one can ‘instantaneously’ determine how the rod temperature varies with any of the key parameters, enabling one to easily probe the underlying physics with questions such as: how does the temperature profile vary with pump size w? Determining this numerically would be extremely time consuming. Furthermore, the resulting analytical expression for the temperature of the rod may be used to determine other quantities, such as the optical aberrations likely to be imparted to a propagating wave through the crystal rod. The key advantages of the analytical model are ease of computation, as highlighted above, and time of computation. We reported that the average computation time of each finite element simulation was approximately 5.5 hours on a dual-core 2.4 GHz processor; this is in stark contrast to the analytical model in which the computational time is less than a couple of seconds on a standard PC, i.e., for all practical purposes it is ‘instantaneous’ in comparison.

The analytical model was used to investigate the thermal stress in a Tm:YLF rod at various qcw pump duty cycles, and was found to be in very good agreement with that of a time-dependent coupled thermal-stress finite element analysis, and with published experimental data. We have applied the analytical model to determine the maximum peak power at which a Tm:YLF rod can be pumped before it will fracture. We show that at all qcw duty cycles the crystal will fracture at average power levels that are lower than predicted by using the corresponding cw fracture limit. The discrepancy is largest at low duty cycles,converging to a perfect agreement in the limiting case of a 100% duty cycle (cw). Thus we suggest a new criteria to be applied for safe power scaling of Tm:YLF rods. While the implementation of the analytical model has concentrated on YLF due to its low fracture limit,the model presented here may be applied to any end-pumped laser rod pumped using qcw or cw sources as the assumptions used to develop the analytical model are not material specific.It needs to be verified for other gain materials and other pumping conditions that the same discrepancy (between the predicted average power level at which fracture will occur and the cw fracture limit) exists at low duty cycles.

References and links

1.	W. Koechner Solid-State Laser Engineering 4th ed., (Springer-Verlag Berlin, Heidelberg, Germany, 1996).
2.	C. Pfistner, R. Weber, H. P. Weber, S Merazzi, and R. Gruber, “Thermal Beam Distortions in End-Pumped Nd:YAG, Nd:GSGG and Nd:YLF,” IEEE J. Quantum Electron. 30, 1605–1615 (1994). [CrossRef]
3.	A. K. Cousins, “Temperature and Thermal Stress Scaling in Finite-Length End-Pumped Lasers Rods,” IEEE J. Quantum Electron. 28, 1057–1069 (1992). [CrossRef]
4.	L. Yan and C. H. Lee, “Thermal Effects in End-Pumped Nd:phosphate Glasses,” J. Appl. Phys. 75, 1286–1292 (1994). [CrossRef]
5.	S. Chénais, S. Forget, F. Druon, F. Balembois, and P. Georges, “Direct and Absolute Temperature Mapping and Heat Transfer Measurements in Diode-End-Pumped Yb:YAG,” Appl. Phys. B 79, 221–224 (2004). [CrossRef]
6.	X. Peng, L. Xu, and A. Asundi, “High-power efficient continous-wave TEM00 intracavity frequencydoubled diode-pumped Nd:YLF laser,” Appl. Opt. 44, 800–807 (2005). [CrossRef] [PubMed]
7.	G. Barton, Elements of Green’s Functions and Propagation , (Oxford University Press, Oxford,1995).
8.	H. S. Carslaw and J. C. Jaeger, Conduction of Heats in Solids , (Oxford University Press, Oxford, 1959).
9.	T. Y. Fan, “Heat Generation in Nd:YAG and Yb:YAG,” IEEE J. Quantum Electron. 29, 1457–1459 (1993). [CrossRef]
10.	S. P. Timoshenko and J. N. Goodier, Theory of Elasticity 3rd ed., (McGraw-Hill, New York, 1970).
11.	B. A. Boley and J. H. Weiner, Theory of Thermal Stresses , (Courier Dover Publications, New York, 1997).
12.	LASCAD Manual, (LASCAD, 2008), https://www.las-cad.com/files/PP_FEA.pdf.
13.	ABAQUS Finite Element Software Package (ABAQUS, 2008) http://www.simulia.com/products/abaqus_fea.html.
14.	E. H. Bernhardi, C. Bollig, L. Harris, M. J. D. Esser, and A. Forbes, “Investigating thermal stresses in quasicw pumped Tm:YLF laser crystals,” in Proceedings of Advanced Solid-State Photonics , (Nara, Japan,2008), Poster WB11.
15.	S. So, J. I. Mackenzie, D. P. Shepherd, W. A. Clarkson, J. G. Betterson, and E. K. Gorton, “A power-scaling strategy for longitudinally diode-pumped Tm:YLF lasers,” Appl. Phys. B 84, 389–393 (2006). [CrossRef]
16.	M. Pollnau, P. J. Hardman, M. A. Kern, W. A. Clarkson, and D. C. Hanna, “Upconversion-induced heat generation and thermal lensing in Nd:YLF and Nd:YAG,” Phys. Rev. B 58, 16076–16092 (1998). [CrossRef]
17.	B. M. Walsh, N. P. Barnes, M. Petros, J. Yu, and U. N. Singh, “Spectroscopy and modeling of solid state lanthanide lasers: Application to trivalent Tm³⁺ and Ho³⁺ in YLiF₄ and LuLiF₄ ,” J. Appl. Phys. 95, 3255–3271 (2004). [CrossRef]

OCIS Codes
(140.3480) Lasers and laser optics : Lasers, diode-pumped
(140.6810) Lasers and laser optics : Thermal effects

ToC Category:
Lasers and Laser Optics

History
Original Manuscript: May 15, 2008
Revised Manuscript: July 2, 2008
Manuscript Accepted: July 9, 2008
Published: July 10, 2008

Citation
E. H. Bernhardi, A. Forbes, C. Bollig, and M.J. D. Esser, "Estimation of thermal fracture limits in quasi-continuous-wave end-pumped lasers through a time-dependent analytical model," Opt. Express 16, 11115-11123 (2008)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-15-11115

Sort: Author | Year | Journal | Reset

References

W. Koechner, Solid-State Laser Engineering 4th ed., (Springer-Verlag Berlin, Heidelberg, Germany, 1996).
C. Pfistner, R. Weber, H. P. Weber, S Merazzi and R. Gruber, "Thermal Beam Distortions in End-Pumped Nd:YAG, Nd:GSGG and Nd:YLF," IEEE J. Quantum Electron. 30, 1605-1615 (1994). [CrossRef]
A. K. Cousins, "Temperature and Thermal Stress Scaling in Finite-Length End-Pumped Lasers Rods," IEEE J. Quantum Electron. 28, 1057-1069 (1992). [CrossRef]
L. Yan, and C. H. Lee, "Thermal Effects in End-Pumped Nd:phosphate Glasses," J. Appl. Phys. 75, 1286-1292 (1994). [CrossRef]
S. Chénais, S. Forget, F. Druon, F. Balembois and P. Georges, "Direct and Absolute Temperature Mapping and Heat Transfer Measurements in Diode-End-Pumped Yb:YAG," Appl. Phys. B 79, 221-224 (2004). [CrossRef]
X. Peng, L. Xu and A. Asundi, "High-power efficient continous-wave TEM00 intracavity frequency-doubled diode-pumped Nd:YLF laser," Appl. Opt. 44, 800-807 (2005). [CrossRef] [PubMed]
G. Barton, Elements of Green???s Functions and Propagation, (Oxford University Press, Oxford, 1995).
H. S. Carslaw and J. C. Jaeger, Conduction of Heats in Solids, (Oxford University Press, Oxford, 1959).
T. Y. Fan, "Heat Generation in Nd:YAG and Yb:YAG," IEEE J. Quantum Electron. 29, 1457-1459 (1993). [CrossRef]
S. P. Timoshenko and J. N. Goodier, Theory of Elasticity 3rd ed., (McGraw-Hill, New York, 1970).
B. A. Boley and J. H. Weiner, Theory of Thermal Stresses, (Courier Dover Publications, New York, 1997).
LASCAD Manual, (LASCAD, 2008), https://www.las-cad.com/files/PP_FEA.pdf.
ABAQUS Finite Element Software Package (ABAQUS, 2008) http://www.simulia.com/products/abaqus_fea.html.
E. H. Bernhardi, C. Bollig, L. Harris, M. J. D. Esser and A. Forbes, "Investigating thermal stresses in quasi-cw pumped Tm:YLF laser crystals," in Proceedings of Advanced Solid-State Photonics, (Nara, Japan, 2008), Poster WB11.
S. So, J. I. Mackenzie, D. P. Shepherd, W. A. Clarkson, J. G. Betterson and E. K. Gorton, "A power-scaling strategy for longitudinally diode-pumped Tm:YLF lasers," Appl. Phys. B 84, 389-393 (2006). [CrossRef]
M. Pollnau, P. J. Hardman, M. A. Kern, W. A. Clarkson, and D. C. Hanna, "Upconversion-induced heat generation and thermal lensing in Nd:YLF and Nd:YAG," Phys. Rev. B 58, 16076-16092 (1998). [CrossRef]
B. M. Walsh, N. P. Barnes, M. Petros, J. Yu, and U. N. Singh, "Spectroscopy and modeling of solid state lanthanide lasers: Application to trivalent Tm3+ and Ho3+ in YLiF4 and LuLiF4," J. Appl. Phys. 95, 3255-3271 (2004). [CrossRef]

Asundi, A.
Balembois, F.
Barnes, N. P.
Betterson, J. G.
Chénais, S.
Clarkson, W. A.
Cousins, A. K.
Druon, F.
Fan, T. Y.
Forget, S.
Georges, P.
Gorton, E. K.
Gruber, R.
Hanna, D. C.
Hardman, P. J.
Kern, M. A.
Lee, C. H.
Mackenzie, J. I.
Merazzi, S
Peng, X.
Petros, M.
Pfistner, C.
Pollnau, M.
Shepherd, D. P.
Singh, U. N.
So, S.
Walsh, B. M.
Weber, H. P.
Weber, R.
Xu, L.
Yan, L.
Yu, J.

Appl. Opt.

X. Peng, L. Xu and A. Asundi, "High-power efficient continous-wave TEM00 intracavity frequency-doubled diode-pumped Nd:YLF laser," Appl. Opt. 44, 800-807 (2005). [CrossRef] [PubMed]

Appl. Phys. B

S. So, J. I. Mackenzie, D. P. Shepherd, W. A. Clarkson, J. G. Betterson and E. K. Gorton, "A power-scaling strategy for longitudinally diode-pumped Tm:YLF lasers," Appl. Phys. B 84, 389-393 (2006). [CrossRef]
S. Chénais, S. Forget, F. Druon, F. Balembois and P. Georges, "Direct and Absolute Temperature Mapping and Heat Transfer Measurements in Diode-End-Pumped Yb:YAG," Appl. Phys. B 79, 221-224 (2004). [CrossRef]

IEEE J. Quantum Electron.

C. Pfistner, R. Weber, H. P. Weber, S Merazzi and R. Gruber, "Thermal Beam Distortions in End-Pumped Nd:YAG, Nd:GSGG and Nd:YLF," IEEE J. Quantum Electron. 30, 1605-1615 (1994). [CrossRef]
A. K. Cousins, "Temperature and Thermal Stress Scaling in Finite-Length End-Pumped Lasers Rods," IEEE J. Quantum Electron. 28, 1057-1069 (1992). [CrossRef]
T. Y. Fan, "Heat Generation in Nd:YAG and Yb:YAG," IEEE J. Quantum Electron. 29, 1457-1459 (1993). [CrossRef]

J. Appl. Phys.

L. Yan, and C. H. Lee, "Thermal Effects in End-Pumped Nd:phosphate Glasses," J. Appl. Phys. 75, 1286-1292 (1994). [CrossRef]
B. M. Walsh, N. P. Barnes, M. Petros, J. Yu, and U. N. Singh, "Spectroscopy and modeling of solid state lanthanide lasers: Application to trivalent Tm3+ and Ho3+ in YLiF4 and LuLiF4," J. Appl. Phys. 95, 3255-3271 (2004). [CrossRef]

Phys. Rev. B

M. Pollnau, P. J. Hardman, M. A. Kern, W. A. Clarkson, and D. C. Hanna, "Upconversion-induced heat generation and thermal lensing in Nd:YLF and Nd:YAG," Phys. Rev. B 58, 16076-16092 (1998). [CrossRef]

Other

G. Barton, Elements of Green???s Functions and Propagation, (Oxford University Press, Oxford, 1995).
H. S. Carslaw and J. C. Jaeger, Conduction of Heats in Solids, (Oxford University Press, Oxford, 1959).
W. Koechner, Solid-State Laser Engineering 4th ed., (Springer-Verlag Berlin, Heidelberg, Germany, 1996).
S. P. Timoshenko and J. N. Goodier, Theory of Elasticity 3rd ed., (McGraw-Hill, New York, 1970).
B. A. Boley and J. H. Weiner, Theory of Thermal Stresses, (Courier Dover Publications, New York, 1997).
LASCAD Manual, (LASCAD, 2008), https://www.las-cad.com/files/PP_FEA.pdf.
ABAQUS Finite Element Software Package (ABAQUS, 2008) http://www.simulia.com/products/abaqus_fea.html.
E. H. Bernhardi, C. Bollig, L. Harris, M. J. D. Esser and A. Forbes, "Investigating thermal stresses in quasi-cw pumped Tm:YLF laser crystals," in Proceedings of Advanced Solid-State Photonics, (Nara, Japan, 2008), Poster WB11.

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper