Heat Generation and Removal in Solid State Lasers

Based on the type of laser gain medium, lasers are mostly divided into four categories; gaseous, liquid (such as dye lasers), semiconductor, and solid state lasers. In recent past years, solid state lasers have been attracted considerable attentions in industry and scientific researches to achieve the high power laser devices with good beam quality. In solid state lasers the gain medium might be a crystal or a glass which is doped with rare earth or transition metal ions. These lasers can be made in the form of bulk [1, 2], fiber [3-7], disk [8, 9] and Microchip lasers [10,11].


Introduction
Based on the type of laser gain medium, lasers are mostly divided into four categories; gaseous, liquid (such as dye lasers), semiconductor, and solid state lasers. In recent past years, solid state lasers have been attracted considerable attentions in industry and scientific researches to achieve the high power laser devices with good beam quality. In solid state lasers the gain medium might be a crystal or a glass which is doped with rare earth or transition metal ions. These lasers can be made in the form of bulk [1,2], fiber [3][4][5][6][7], disk [8,9] and Microchip lasers [10,11].
Optical pumping is associated with the heat generation in solid state laser materials [12]. Moving of heat toward the surrounding medium which is mostly designed for the cooling management causes thermal gradient inside the medium [13]. This is the main reason of appearance of unwanted thermal effects on laser operation. Thermal lensing [14], thermal stress fracture limit [15], thermal birefringence and thus thermal bifocusing [16][17][18][19] are some examples of thermal effects.
Optimizing the laser operation in presence of thermal effects needs to have temperature distribution inside the gain medium. Solving the heat differential equation beside considering boundary conditions gives the temperature field.
Boundary conditions are directly related to the cooling methods which lead to convective or conductive heat transfer from gain to the surrounding medium.
In the case of bulky solid state laser systems (such as rod shape gain medium), water cooling is the most common method which is almost used in high power regime. Design of optimized cooling cavity to achieve the most effective heat transfer is the first step to scale up high power laser devices. Then determination of temperature distribution inside the laser gain medium is essential for evaluation of induced thermo-optic effects on laser operation. This chapter is organized according to the requirements of reader with thermal considerations in solid state laser which are mainly dependent on several kinds of laser materials, pumping procedure and cooling system. We hope the subjects included in this chapter will be interesting for two guilds of scientific and engineering researchers. The first category relates to the laser scientists, who need enough information about the recent cooling methods, their benefits and disadvantages, thermal management and effects of utilizing a specific cooling method on laser operation. And the second one is the mechanical or opto-mechanical engineers who are responsible for designing and manufacturing of the cooling systems. In this chapter our efforts directed such a way to satisfy both the mentioned categories of researchers.
At First, the principle of heat generation process inside the laser gain medium due to the optical pumping are introduced. Then, heat differential equation in laser gain medium and relating boundary conditions are introduced in detail. Formulation of heat problem for a specific form of gain medium such as bulk, disk, fiber and Microchip lasers and details of solution are presented through individual subsections.

Principle of heat generation in solid state laser gain medium 2.1. Laser pumping
A laser device is composed of three essential components which are the "active medium", "pumping source" and the "optical resonator". In the case of solid state lasers, the active medium which is made of a definite glass or crystal, is placed inside the optical resonator and receives energy from another external optical source through the pump beam light. Then it can itself emit an amplified laser beam light delivering a completely modified energy and wavelength [20]. The act of energy transfer from the external source to the active medium is called the laser pumping.
In recent years, diode lasers [21][22] have attracted considerable attention between laser scientists as available, high power and beam quality pumping sources. In this chapter we just concentrate on this kind of pumping sources rather than traditional flash lamp pump sources [15].
The pumping process commonly performs in two methods, which are continuous wave (CW) and pulse pumping laser systems. Furthermore, diode pumped solid state lasers (DPSS lasers) can be divided into side pumped and end pumped configurations. Figure 1, Shows schematically typical solid state lasers, including gain medium, optical resonator and diode pumping in the case of end and side configurations. In the end-pumped geometry, the pump light mostly transfers from diode Laser (DL) to the laser material through either optical system or fiber optics which yields a desirable pump beam shape and size. Then it focuses to the gain medium longitudinally, collinear to the propagation of laser light. In the side-pumped geometry, the diode arrays locate along the laser material in a definite arrangement around it, such that the pump light is perpendicular to the propagation of laser light.

Heat Generation and Removal in Solid State Lasers 343
The pumping geometry and the resultant pump characteristics (such as beam shape and size along the gain medium) play an important role in heat generation and therefore thermal gradient inside the gain medium. The issue will review in details in the following subsections for each kind of laser gain mediums. Simple drawing for two common pumping methods of solid state laser gain medium; a) side pumping b) end pumping. The blue color used for marking the pump light carrying energy from laser diode to the gain medium, and the red beam concerns to the laser resonator mode.

Heat generation
In solid state lasers, a fraction of the pump energy converts to heat which acts as the heat source inside the laser material [23,24]. Spatial and time dependence of the heat source causes important effects on temperature distribution and warming rate of the gain medium, respectively. The spatial form is assumed to be the same shape as pumping light [23,24] and time dependence relates to the pumping procedure, which may perform in CW or pulsed regime. Furthermore, depending on the gain medium configuration and cooling geometry, deposited heat may mostly flow through a preferable direction inside the gain medium and therefore causes thermal gradient. For instance, in traditional rod shape laser mediums with water cooling configuration and also fiber lasers, the main proportion of heat removal occurs through the radial direction which leads to the considerable radial thermal gradient inside the medium. Figure 2 shows a schematic setup of pumping procedure for several kinds of solid state lasers. The dominant directions of heat removal which are associated with the gain medium and cooling system geometry are illustrated. Heat differential equation should be solved for evaluation of temperature and thermal gradient induced by optical pumping in solid state lasers. The general form of heat differential equation in cylindrical coordinate system is given by [25]  The overall thermal load in the laser gain medium due to the optical pumping can be obtained from In which, o P  is the pump power and is the fractional thermal load [12]. In the case of diode pumping, the fractional thermal load, originates from two basic phenomenon, which show the main role in heat generation; quantum defect heating [15] and energy transfer upconversion (ETU) [26]. In most cases, the first is responsible for the heat generation and therefore has the main contribution. However it must be noted that, influence of the second phenomena cannot be ignored in some cases such as Er doped laser materials. The fractional thermal load in the gain medium is due to the quantum defect and related to the pumping and laser wavelength which are shown by p  and L  respectively.

. Side pumping
Pumping configuration performing by one module is illustrated in Figure3.a [27]. The pump beam emitted by diode laser is focused on the rod by the interfacing optics which consist of two lenses. End view of side-pumping geometry is depicted in Figure 3.b. Transversal directions of pump and signal beams are easily observable. In the case of CW pumping of a laser rod, the steady state heat equation can be written as [27] .
Where, h is heat flux and associated with the temperature in the rod by And (,) Qrz is determined according to spatial variation of pump intensity and is given by

End pumping
One of the common pumping configurations which are used in diode-pumped solid state lasers is end-pumping or longitudinal pump scheme. The pumping beam is coaxial with the resonator beam in end-pumped lasers; it leads to highly efficient lasers with good beam quality.In this geometry, the pumping beam of diode laser(s) is delivered to the end of the active medium by optical focusing lenses or optical fibers. In lower power operations (less than a few watts), end-pumping yields more acceptable results [15]. Today's development of diode lasers and new techniques such as using micro-lenses in beam shaping of diode laser bars make end-pumped lasers very promising specially in commercial lasers [28]. Although end-pumping is a common configuration in solid state lasers which include many types of active medium geometries such as slab and microchip, but it is more commonly used in rod shape lasers. Many end-pumping investigations and results have reported about rod lasers and most of commercial end pump systems involve rod lasers [29][30][31][32][33]. Thus, the discussions in this section are focused on end-pumped rod lasers. Schematic diagram of an end-pumped system is shown in figure 4. . End pump systems major elements [15] In order to establish high matching efficiency between resonator beam and pumping modes, in end pumped systems the pumping beam is focused in the active medium with the small beam waist. This issue leads to generation of an intense local heating and then, creation of refracting index gradient inside the laser crystal [15]. As a consequence, the laser rod acts as a thermal lens inside the resonator, which can destroy the beam quality and decrease the output efficiency. Additionally, in contrast to side pumped lasers, heat distribution within the laser material is inhomogeneous in high-power end-pumped lasers, leading to increased in stress and strain [33]. The restricting factors in end-pumped lasers in high power regime are thermal lensing and thermal fracture limit for the laser crystal. The aforementioned restrictions make thermal problems very important in end-pumped lasers especially in highpower systems.
From the thermal point of view, the flat top pumping profile is superior to Gaussian profile in high power end-pumped systems, due to the creation of lower thermal gradient inside the laser crystal leading to lower thermal distortions [28]. However, Gaussian pumping profiles is more investigated because of practical considerations in laser resonator design. One of the earliest thermal analysis in end pumped systems is presented in [15], which relates to the solution of steady state heat differential equation for Nd:YAG crystal with the assumption of Gaussian pumping profile. ( Figure 5). The temperature profile and thus associated thermal effects in bulky solid state lasers had been the subject of consideration in past years by various authors. In the case of CW pumping, a collection of excellent literatures discussing numerical and analytical thermal analysis can be found in [13,23,[34][35][36][37][38][39]. Similar works concerning the transient heat analysis are available in [40].
A famous work which presents analytical expressions for temperature and describes the behavior of temperature inside the rod belongs to Innocenzi [13]. In this work, the laser rod is surrounded by a copper heat sink and exposed to the pump beam with the Gaussian intensity profile as Neglecting axial heat flux, the steady state heat differential equation can be solved analytically. The derived temperature difference is obtained as As would be expected, the temperature decays exponentially along the laser rod and have the highest temperature on the pumping surface (z=0).
A common cooling method for end-pumped systems concerns to utilize of water jacket or copper tube surrounding the laser rod and keep the cylindrical surface at the definite temperature. Heat generated inside the gain medium flows to this surface through the heat conduction process in radial direction. Although this method is considered as a simple efficient technic providing considerable heat removal from gain medium, but the uncooled pumped surface of the rod which is in direct contact with air, performs very week heat transfer. This issue may cause undesirable effects, especially in high power regimes [34]. The thermal load on this surface not only increases thermal lensing effects but also restricts the maximum pumping power because of the fracture limit of the crystal, and damage threshold of optical coatings.
One of the effective strategies to reduce thermal effects is based on the cooling of pumping surface of laser rod. In this respect, three technics to achieve more efficient cooling process have been presented in [33], which are schematically shown in figure 6. In the first method (b), the cooling water is directly in contact with the pumping surface. In the second method (c), a cooling plate cooled by water is mounted in close contact with the rod pumping surface. The cold plate should have large Yang's module and high heat conductivity. The other method utilizes an un-doped cap on the pumping surface (d). The pumping power does not absorb in the undoped cap so there is not heat load in this region, but using this cap increases the effective cooling surface as well as the ratio of cooling surface to heat generation volume.
The influence of thermal effects on laser gain medium in the mentioned methods has been analyzed using FE method [33]. The maximum temperature decreased almost 30%, and 25% using an undoped cap and a sapphire cooling plate in pumping surface respectively. The maximum stress occurred in the configuration with water cooled pumped surface and reduced to half of uncooled system value using an undoped cap or cooling plate. According to recently developed ceramic laser materials using composite rods with undoped cap could be very promising as the best choice for high power end pumped systems. The undoped This not only reduces the thermal lensing effects and thermal stresses but also lowers the maximum temperature of the laser rod and so removes some constraints imposed on the coatings. Prominent role of composite rod in reduction of thermal destructive effects on laser operation have been frequently examined and reflected in literatures [41][42][43]. Figure 7 illustrates the pump model of the dual-end-pumped geometry of composite Nd:YVO4 Figure 7. Pump modeling of the dual-end-pumped geometry [44].
laser [44]. The Nd:YVO4 as the laser gain medium is connected to two YVO4 caps at two ends and the pump energy lunches to it from both ends. As can be seen, every point inside the rod absorbs pump power and experiences heat generation. No absorption inside the caps takes place and therefore, they can show important role in axial heat transfer from end surfaces. Numerical calculations of temperature distribution in composite laser rod can be found in [45] In the case of pulsed pumping laser rod, interesting numerical analysis has been done by Wang published in [46]. In this work the laser rod is surrounded by a cylindrical heat sink which leads to conductive heat transfer from rod surface to the ambient medium. Schematic figure of the rod and cooling system geometry is depicted in Figure 8-a. The laser rod is assumed to be coupled by a fiber-optic to a laser diode; therefore the pump intensity profile with a good approximation has the top-hat shape. Thus the heat source density can be described by Where P(t) is a periodic function of time describing pump power in a repetitively pumping regime given by ,0 In which Po is the pump peak power, rep f is the repetition frequency, and τ is the pump duration time. Figure 8.b tries to simply specify the P(t). Detailed information of solving the transient heat differential equation can be found in [46]. Figure 9 shows the temperature distribution inside a Nd:YAG laser rod exposing to the pump beam with pulse duration of 300μs, and repetition frequency of 50 Hz. According to the results, the temperature at first increases by passing the time and then becomes nearly constant with small fluctuations which lead the temperature to the steady state condition.

Thin disk lasers
Thin disk lasers are one of the recent frontiers in solid state lasers. The most important features which make the thin disk laser distinguishable between solid state lasers are power scalability, good beam quality and minimal thermal lensing [47,48]. These features are related to the thermal characteristics of the thin disk laser. In disk lasers, active medium is cooled from the disk face. The surface to volume ratio of the disk is large due to the disk geometry, therefore cooling is very efficient and as a result thermal distortion of the active medium is very low. Considering axial heat flow in a thin disk laser there is no thermal lensing in a first-order approximation. In fact, however weak thermal lensing occurs because of two residual effects: first, pumped diameter is typically smaller than the diameter of the crystal and second, thermo-mechanical contribution to the thermal lensing from bending of the disk due to thermal expansion [49].
Thermal lensing is important issue in laser design and operation. This factor can be calculated theoretically using thermo-mechanical modeling softwares. In thin disk lasers, the disk is mounted with a cold plate on a heat sink (figure 10). At the same time the other side of disk is radiated by pumping laser, accordingly there is a temperature difference between two faces of the disk. This will cause a temperature distribution through the disk bulk. Generally the refractive index of materials is depended on temperature; accordingly the refractive index of disk will be a function of position. The other effect is the expansion of the disk due to the temperature distribution formed in it. Also mounting the disk on heat sink causes deformation and stress in the disk. The stress itself will affect the refractive index of the disk crystal. To complete analysis of the effects of the disk on the laser and pump beam, one must calculate cumulative effects of expansion and deformation of the disk, also thermo-optical and stress dependent variations of refractive index [50]. Total effect of the disk on laser beam phase can be calculated by [51]: In which 0 n is refractive index of disk at reference temperature 0 T , nT   is the thermooptical coefficient, s n  is the changes of refractive index due to stresses, z  is the strain in direction of thickness of the disk, 0 z is the displacement of back side (High Reflection coated) of the disk and h is the thickness of the disk. As relation (1) shows, optical behavior of the disk is strongly depended to the temperature distribution of disk. The temperature distribution is result of optical pumping and cooling of the disk.

Pump and cooling configurations
There are two conventional methods for pumping disk lasers; first is (quasi) end pumping and the second is edge pumping. Schematic diagram of end pumping is shown in figure 10. Also figure 11 shows a schematic setup of the edge pumped thin disk laser.
In both mentioned methods, the disk is cooled from the face. The disk can be cooled by jet impingement ( figure 12) or cryogenic technique. The disk is mounted with a cold plate on the heat sink. In jet impingement a jet of cooling liquid is sprayed to the cold plate. Different liquids can be considered as coolant which most common is water.
Heat Generation and Removal in Solid State Lasers 353 Figure 11. Schematic setup of edge pumped thin disk laser [53] Figure 12. Schematic diagram of jet impingement cooling system for thin disk laser [54] Laser cooling has been an important problem from the invention of the first practical laser in 1960. After invention of laser, cryogenic cooling of solid state lasers has been interested and first time proposed by Bowness [55] in 1963 and then by McMahon [56] in 1969. In mentioned references the conduction cooling is used and laser element was placed in contact with a material with very high thermal conductivity. That material was, in turn, in contact with a cryogen such as liquid nitrogen near 77 K, liquid Ne near 27 K, or He near 4 K.
At cryogenic temperatures the absorption and emission cross sections increases and the Yb:YAG absorption band near 941 nm narrows at 77 K, however it is still broad enough for pumping with practical diode lasers. At 77 K, Yb:YAG crystal behaves as a four-level active medium however in room temperature Yb:YAG is a quasi three-level material. Significant problems associated with quasi-three-level materials like Yb:YAG, such as the need to provide a significant pump density to reach transparency, a high pump threshold power, and the associated loss of efficiency, disappear at 77 K [57]. When the Yb:YAG is cooled from room to cryogenic temperatures, the lasing threshold decreases and slope efficiency increases. Figure 13 shows drop of lasing threshold from 155 W to near 10 W and increasing slope efficiency from 54% to a 63% for a typical thin disk laser [58]. The inset spectral peak in Figure 13 is the 80 K laser output and is centered near 1029.1 nm. At room temperatures this peak is near 1030.2 nm.

Temperature distribution
Specifications of disk laser beam are tightly related to the active medium geometry. The precise geometry of the active medium geometry is also strongly depended to the thermomechanical and opto-mechanical properties of the disk and the temperature distribution in the disk. Temperature distribution of the disk can be obtained by solving the heat conduction equation using proper boundary conditions.
The flow of heat generated by the pumping diode radiation through a laser gain media in general form is described by a non-homogeneous partial differential equation: In which, Kc is heat conductivity that is assumed to be isotropic and is the heat source density in the laser crystal. In CW regime of output laser, the steady-state temperature distribution obeys the heat diffusion equation Heat Generation and Removal in Solid State Lasers 355 As it is seen the heat conductivity is depended on the temperature and doping concentration (dop) of the crystal. Temperature dependency of YAG crystal in room temperature is not significant and it can be considered as a constant [37] however this approximation would not be valid anymore at cryogenic temperature. Heat conductivity of Yb:YAG crystal which is conventional active medium of thin disk laser can be given by [59]: Characterizing the behavior of the thermally induced lensing effect of the thin disk gain medium is not a trivial task. In order to fully analyze the dynamics of the heat flow and thus the induced nT stresses and strains on the gain medium one must solve the 3-D heat equation with appropriate boundary conditions. This can be accomplished in several ways. The most common is to employ a finite element analysis (FEA) method. Another method is to solve the 3-D heat equation using the Hankel transform. For more details the reader can refer to [60].
Initial estimation of disk thermal behavior can be carried out by calculation of maximum and average temperature of the disk. In thin disk lasers, the thickness of the disk is very low. When the pump spot size is very larger than the disk thickness, one dimensional heat conduction is a good approximation. If pump power of pump P radiates to the disk in a pump spot with radius of p r , the heat load per area can be given by [ In which abs  is absorption efficiency and  is heat generation coefficient in the disk. The heat generation coefficient in the disk is due to the quantum defect and related to the pumping and laser wavelength which are shown by p  and l  respectively.
A parabolic temperature profile will be formed along the axis inside the disk due to the loaded heat which is given by: in which , / th disk c Rh k  is the heat resistance of the disk material and 0 T is the temperature of the disk's cooled face. Also z is the distance along the disk axis in the thickness of the disk and h is the thickness of the disk. In particular, one can calculate maximum temperature from relation (19) which is given by In this way using relations (19) to (21) one can evaluate, temperature distribution and maximum and average temperature in the disk in one dimensional heat conduction approximation.

Introduction to fiber geometry and cooling methods
In recent years, design and manufacturing of high efficient fiber lasers which deliver excellent beam quality, has made them as the main adversary of other types of high power solid state lasers, such as bulk and thin-disk lasers. Achieving to multi-kilowatt output powers [7,62] with diffraction limited laser beam could be considered as the unique record in laser technology. This progress can be attributed directly to the capability of more efficient cooling procedure in fiber lasers, which originates from inherent large surface to volume ratio. In fact, thermal load spreads over meters or tens of meters of fibers, which causes convenient and efficient cooling management and therefore avoids thermo-optic problems.
Fiber lasers are consists of fiber core which is mostly surrounded by two coaxial fiber cladding (double clad fiber lasers) and is pumped by diode bars or diode laser at one or both ends. The laser light can only propagate through the fiber core and doesn't have any role in heat generation inside the fiber.
There are two general methods for lunching pump light into the fiber laser which are called as "core pumping" and "cladding pumping". The conventional core pumping was initially used to achieve single mode output laser, in which the pump light was coupled into the small core. On the other hand, small core causes a serious restriction on pump power level [63]. Furthermore, the core size leads to highly localized pump intensity which usually induces thermal damage at the fiber ends. Therefore, cladding pumping has been developed as the proper solution which ensures lunching high pump power into the double clad fiber lasers. In this method the pump light couples into the inner cladding and propagates through it and gradually absorbs in doped core. In both cases, the pump light only absorbs within the core, where heat generation takes place. Figure 14 shows a simple diagram of cladding pumping of fiber laser geometry [63].
In most cases, cooling procedure in fiber lasers does not need any special cooling system and are called passive cooling, which easily can be performed by the air through the convectional process [62][63][64]65]. However, in modern fiber lasers an active cooling system is Heat Generation and Removal in Solid State Lasers 357 Figure 14. Cladding pumped fiber amplifier [63] considered to scale up high power lasers, which ensures a forced heat dissipation process. Liquid cooled fiber lasers [66] is an example of new cooling methods in which, the whole or some part of the fiber is placed inside a liquid with a definite temperature. Therefore, heat removal occurs through the convection from the fiber periphery to the cooling liquid. This technique is usually applied to the long fiber lasers. Another technique which is often utilized for cooling of short fiber lasers, concerns to the thermoelectric cooling system (TEC). In this case, the fiber medium surrounds by a copper heat sink and therefore heat removal performs through the heat conduction process. The three common methods which imply to the passive and active cooling techniques are introduced in more details in follow.
Example of conductive boundary condition in which a short length fiber is surrounded by a temperature controlled copper heat sink can be found in [67]. Ignoring of axial heat flux as an approximation, heat differential equation can be solved numerically by means of Finite Element (FE) method. Figure 15, shows a drawing of a TEC-cooled fiber assembly. Cooling of long fiber laser based on the conductive heat transfer is reported in [67]. The fiber is placed between aluminum plates with constant temperature caused by water cooling.
Practical models of high power fiber lasers with the unforced convective heat transfer from fiber to air are reported in [62][63][64]65]. Figure16, illustrates the experimental set up for an Er: An Overview of Heat Transfer Phenomena 358 ZBLAN double-clad fiber laser. Fiber is placed inside the resonator and is pumped from one end by a diode laser after passing the pump beam from the designed optics. Figure 16. experimental set up for high power Er: ZBLAN fiber laser [64]. The fiber is pumped by a diode laser at one end. Figure 17 shows another example relating to high power Yb doped Fiber laser (YDFL), which is pumped from both ends [62]. Convectional cooling process from fiber to air is freely established. Pump power is delivered from two diode stacks propagate from the both ends toward the fiber center and cause two individual heat sources inside the fiber. Efficient fiber cooling leads to scale up high power lasers without thermal damage and avoiding destructive thermal effects on laser operation. A new technique for thermal management of fiber was examined in [66], which called direct liquid cooling. In this method the fiber was in direct contact with the fluorocarbon liquid. Furthermore, the both ends of fiber facet are in physical contact with the CaF2 windows. This leads to conductive heat transfer from fiber facet to the window and considerable axial heat removal, which allows increasing pump power without thermal damage. This technique was already used in composite bulky solid state laser mediums [24,[41][42][43] and had found highly operative [68,69]. Figure 18, shows a drawing of system assembly. The fiber is pumped by two fiber coupled diode lasers from both ends.

Continuous Wave (CW) pumping conditions
Using CW pump sources lead to generation of time independent heat source density in fiber core. Therefore, Eq.1 turns to steady state heat differential equation as In which, the azimuthal part of the temperature is omitted due to the cylindrical symmetry of pumping spatial distribution. As we mentioned before, spatial form of heat source density obeys from the spatial form of pump intensity profile lunched to the fiber. "Top hat" and "Gaussian" are two common shapes for the pump beam profile, which are usually considered as the spatial form of heat source density in thermal analysis. Different considerations lead to different differential equations and therefore, needs different solutions. Analytical and Numerical solutions of Eq.1, to specify the temperature behavior inside fiber medium, are reported in various literatures with different approximations and methods.
As we mentioned before, different cooling arrangements lead to different boundary conditions which are conductive or convective heat transfer from fiber periphery to the surrounding medium.
In the case of fiber coupled fiber laser, pump intensity distribution with a good approximation has a top hat profile across the beam. Entering the pump beam inside the fiber core and propagating along the fiber length causes exponential decay in axial direction. Therefore, the heat source density Q(r,z), inside the fiber can be expressed by [70] Where, o P is the pump power, a is radius of the fiber core, b is the radius of outer cladding, is the effective fiber length, L is the geometrical fiber length and  is the effective pump absorption coefficient. Substitution of (,) Qrz from Eq. 23 into Eq. 22, the heat differential equation for two regions can be written as  -a corresponds to the fiber core, where exposes to the pump beam and therefore experiences the heat generation. Equation 24-b relates to fiber cladding that is just responsible for heat transmission to the surrounding medium. Solution of these differential equations give the steady state three dimensional temperature at any point in fiber core T1(r,z) and cladding T2(r,z). Figure 19, illustrates schematically double clad circular fiber geometry under pumping process [70].  Figure 19. End (a) and side (b) view of circular fiber laser. Absorption of pump power inside the fiber core causes heat generation which reduces exponentially along the fiber axis [70].
Assuming the outer surface of fiber is in direct contact with a liquid or gas such as air, the convective boundary condition could be defined according to Newton's law as 2 2 (,) (( , ) ) Where, h is heat transfer coefficient [71],Tc is the coolant temperature andT2 (b,z) represents the temperature along the fiber on its cylindrical surface and the first derivative is taken with respect to the surface normal. This equation expresses that the radial thermal flux which arrives at the fiber periphery via the conduction method, removes by the coolant through the heat convection process. Further information about the analytical solution of Equation and some technical points on numerical approach can be found in [70]. The resultant temperature expressions are Where C, A1 and A2 are constant coefficients which are derived from the boundary conditions as follows The temperature is reduced exponentially along the fiber axis due to the exponential decay of heat deposition. Radial dependence in the core includes zero order Bessel function and inside the cladding, linear combination of the first and second kind of zero order Bessel functions are established. This is illustrated in Figure 20 which shows the calculated temperature distribution in the r-z plane of ZBLAN (ZrF4-BaF2-LaF3-AlF3-NaF) double clad fiber. More information about the fiber geometry and pumping characteristics is given in Table 1.  In the case of short length fiber laser, analytical approach of temperature expressions for an Er 3+/ Yb 3+ co-doped fiber laser can be found in [26]. Schematic drawing of fiber pumping is illustrated in Figure 21. In short fiber lasers, the pumping energy doesn't absorb completely by one passing of light along the fiber length and thus, the reflected part from the coupler acts as a pump beam delivers inside the fiber at the end face. Heat source function is composed of the proportion of the both the left and right traveling pump beams. The heat source density function with the assumption of Gaussian pump beam shape is given by [ wavelength respectively and p  is the Gaussian radius of the pump light. Figure 22 illustrates an end view drawing of the modeled fiber. As it can be seen, the active fiber is surrounded by a ceramic ferrule which is placed inside a copper tube. Heat conduction from fiber cladding to ceramic ferrule occurs through the fiber cylindrical surface.
The steady state heat differential equations are Where, 0 n  is nth zero point of the zero order Bessel function of the first kind, 0 J and 1 J are the zero and first order Bessel function of the first kind, 0 I , 1 I and 0 K , 1 K are the zero and first order modified Bessel function of the first and second kind, respectively [26 ]. m C , m D and nm A are the constants which derived from boundary conditions, via solution of coupled equations and are given in [26]. The above equations are plotted in figure 23, illustrating the temperature distribution in fiber core and cladding in a phosphate fiber. As would be expected, the maximum temperature achieved at the center of the fiber facet and is 479.85K. Amount of some parameters used in calculations are summarized in Table 2.

Pulsed pumping conditions
An example of pulsed pumped fiber laser which causes transient temperature field is presented in [72]. This work relates to short length fiber laser and introduces analytical temperature expressions for solution of transient heat differential equations which are 22 22 ( Where, Pin is the pulse energy and g(t) is temporal shape of pump pulse. The other parameters were previously introduced. The equations (38)- (40) are applicable here as the boundary conditions. Furthermore, the laser system is assumed to be in thermal equilibrium with the ambient air before pumping process which expresses by Solution of equations (36) and (37) In which, n denotes number of pulses, 0 t is the pulse duration and 0 T is pulse period. Figure 24, shows the temperature distribution achieved from the analytical temperature expressions for a typical Er 3+ /Yb 3+ co-doped fiber laser. Characteristics of the fiber laser which was introduced in Table 2 were considered in calculations. The fiber is in encounter of pulses with the pulse pump of 1 in PW  . Entering of each pulse into the fiber core causes definite heat load and thus leads to increase of temperature. Figure 25-a, shows the temperature evolution at fiber end facet during the space time of pulse generation. The curves correspond to the time increment of 1ms.
If the pulse period is longer than the space time of heat removal, the temperature will return to the initial value before loading the next pulse. But if the pulse period is shorter than the time which is necessary for completely remove of generated heat, the heat will gradually deposit inside the fiber medium and leads to increase of the fiber temperature. Figure25-b illustrates reduction of fiber facet temperature after the first pulse. At the time of 15ms temperature profile has nearly flat shape but the temperature in all positions is 2 C  over the initial temperature.

Microchip lasers
In microchip lasers the coated surface of active medium is used as laser's mirror. Popular types of microchip lasers have less than 1 Watt output power. However, recently some reports are demonstrated microchip lasers with several hundred Watts [74][75][76][77][78]. Due to the importance of thermal problems in high power operation, we focus our investigation on this regime. Generally, composite active media are used in high power microchip lasers. The laser material includes a central doped region as the active medium surrounded by an undoped rim as pumping waveguide.
The microchip laser has a thin active medium with geometry very similar to that of thin disk lasers. Considering the uniform absorbed power distribution, the thermal considerations assumed in thin disk lasers are still valid in microchip lasers. Most common microchip lasers in high power operations are composite microchips with an un-doped rim surrounding active region acting as pump wave guide. The main drawback in this geometry is the non-uniformity of absorbed power distribution which leads to non-uniform temperature distribution. Regarding boundary conditions and non-uniform heat deposition, the temperature distribution in the active medium should be investigated.
An example of thermal analysis in high power microchip laser is presented in [79]. Utilizing optimum geometry and dimensions for active medium, one can consider uniform absorbed power distribution in microchip lasers which leads to relatively uniform temperature distribution and can overcome the common disadvantage of microchip lasers. A common geometry of the Yb:YAG/YAG composite microchip lasers is illustrated in figure 27 [79].  [79] In this arrangement, a thin disk Yb:YAG core enclosed by an irregular symmetric eight sided un-doped YAG rim is bonded on a water-cooled Cu-W cold-plate.
Fraction of absorbed pump power in the gain medium produces heat in the core and raises the core temperature. The temperature distribution in the laser material is calculated using the, which is Thermal conductivity is strongly dependent on temperature and atomic doping concentration of Yb ions CYb as follows [80] Heat The absorbed power distribution in active medium is calculated using Mont-carlo ray tracing of pumping photons through the optical elements of the system. Heat generation Q is estimated based on the energy difference between the pump and laser photons, given by [81 ] (1 ) In which Pp is local absorbed power density. The heat transfer between the gain medium and air is week, which is described by In which z axis coincides on microchip optical axis and t denotes the microchip thickness. Moreover, the temperature of the back face of the cold-plate and contact cooling water were considered as the same. The numerical values used in calculations are given in Table 3.

Quantity magnitude
Pumping wavelength 942nm Laser wavelength 1030nm active medium thickness 0.3mm cold plate thickness 1mm thermal conductivity [9-82] cold plate 1.9 W/cm.K cold plate heat transfer coefficient 12 W/cm2.K Table 3. The microchip laser parameters used in numerical calculations.
Heat differential equation (49), was calculated in numerical Finite-Difference method for the described gain medium [79]. This code is used to determine the temperature distribution in the active medium. Figure 28 shows the maximum temperature variation versus a/w ratio for various core diameters with 9 at.% Yb3+ doping concentration. Moreover, the temperature distribution profiles of these cores are also shown in Figure 29.
According to figure 29, the laser gain mediums with the smaller core size experience higher temperatures at the same pump power, which is attributed to the higher thermal load density.