Research Article  Open Access
Iliycho Petkov Iliev, Snezhana Georgieva GochevaIlieva, Krassimir Angelov Temelkov, Nikolay Kirilov Vuchkov, Nikola Vassilev Sabotinov, "Modeling of the Radial Heat Flow and Cooling Processes in a Deep Ultraviolet NeCuBr Laser", Mathematical Problems in Engineering, vol. 2009, Article ID 582732, 16 pages, 2009. https://doi.org/10.1155/2009/582732
Modeling of the Radial Heat Flow and Cooling Processes in a Deep Ultraviolet NeCuBr Laser
Abstract
An improved theoretical model of the gas temperature profile in the crosssection of an ultraviolet copper ion excited copper bromide laser is developed. The model is based on the solution of the onedimensional heat conduction equation subject to special nonlinear boundary conditions, describing the heat interaction between the laser tube and its surroundings. It takes into account the nonuniform distribution of the volume power density along with the radius of the laser tube. The problem is reduced to the boundary value problem of the first kind. An explicit solution of this model is obtained. The model is applied for the evaluation of the gas temperature profiles of the laser in the conditions of free and forced aircooling. Comparison with other simple models assumed constant volume power density is made. In particular, a simple expression for calculating the average gas temperature is found.
1. Introduction
Deep ultraviolet (DUV) gas laser sources have been objects of great interest in the recent years, because of a wide variety of applications, such as highprecision processing of different materials, highresolution laser lithography in microelectronics, highdensity optical recording of information, laserinduced modification in various materials newly developed, as well as laserinduced fluorescence in plasma and widegap semiconductors [1–7]. These applications require a DUV laser source in order to achieve the necessarily high resolution or to induce some phenomena in different materials, such as material and its surface modification, image recording, and fluorescence. The photon energy is sufficient first, to produce photochemical ablation rather than melting and subsequent vaporization, or second, to initiate the induced processes. The DUV laser ablation produces clean precision cutting and drilling with a minimum of thermal or mechanical damages on the target. DUV beams, in addition to heatless etching, offer the advantage of focusing with less limitation due to diffraction. The good spatial structure of the laser beam is generally important for the goodquality image projected, while narrow linewidth would reduce the chromatic aberrations. Chromatic aberrations are very difficult to correct because of the limited number of optical materials that transmit in the DUV range. The excimer lasers, which are generally used as light sources in the DUV range, can neither reach the spatial structure nor the linewidth required. Their wide application is also impeded by their cost and considerable overheads. Lately, there has been widespread interest in the ion metal vapor lasers, operating below 300 nm in the DUV spectral region. High beam quality and narrowlinewidth are inherent for metal vapor lasers. For the first time laser oscillation on four DUV Cu^{+} lines: 248.6 nm, 252.9 nm, 260.0 nm, and 270.3 nm, in a nanosecond, pulsed longitudinal NeCuBr discharge was obtained in [1]. In [2] an active zone diameter scanning of the DUV Cu^{+} NeCuBr laser was carried out and the optimal discharge conditions were found for each active zone diameter. A record average output power of 1.3 W was obtained at multiline operation. A record specific average output power was 57 mW/cm^{3} at an active volume of 23 cm^{3}. The highest peak pulse power and average laser power on the 248.6 nm laser line for the DUV Cu^{+} lasers were 3.25 W and 0.85 W, respectively.
For the metal and metal halide vapor lasers, the thermal mode, and in particular the radial temperature distribution, is of great importance for the stability of the laser operation and for the achievement of high output characteristics as well [8, 9]. In addition, crosssections and rate constants for heavy particle collisions, which thoroughly determine the inverse population and laser output parameters, such as asymmetric charge transfer, Penning ionization, depend strongly on the gas temperature. That is why the radial temperature distribution in the active zone of the DUV Cu^{+} NeCuBr laser is object of a detailed investigation. The gas temperature problem is also topical in radiofrequency discharges, widely encountered in gas lasers, electron and ion technologies, nanotechnology, plasmachemical technology, and more [8].
It is well known that the experimental techniques for gas temperature measurement, using spectral lines broadening and thermocouple are definitely imprecise. That is why mathematical and computer models are applied to carry out the temperature distribution. These models are used not only to evaluate the temperature but also for complex modeling of plasma kinetics processes [9–14].
The determination of the gas temperature profile in metal vapor lasers has been firstly considered in [9], where selfconsistent mathematical models have been developed based on an exact solution of the onedimensional steadystate heat conduction equation subject to the simple boundary conditions of the first and second kind. Till now, the same method has been used in all other papers in this area, when the buffer gas temperature is calculated, see for instance, [9–16]. In these studies, the expected average temperature and profile of the gas are obtained using the following simplifying assumptions: the average power input is considered uniform in the discharge—no radial or longitudinal variations in power density are considered (); the wall temperature is assumed to be constant.
However, in practice, and in computer modeling and simulations, these assumptions are not fully applicable. That is so because the real values of the volume power density are not uniformly distributed and vary strongly along the tube radius taking its maximum in the central axis of the discharge. In addition, the temperature of the outer surface of the laser tube under insulation is unknown and will change with variation of the laser geometry, input electric power, and other laser parameters including temperature of the surroundings.
In [17, 18], the temperature profile in the case of a copper bromide vapor laser with wavelength 510.6 nm and 578.2 nm was determined by using a new approach. It is based on solving the heat conduction equation with at nonlinear boundary value conditions for a given temperature of the surroundings.
In this study, for the first time, the analytical investigation of the temperature profile in the crosssection of the laser tube is performed with the assumption of a specified qualitative distribution of , dependent on the tube radius, namely, . Using the approach similar to this in [17, 18], in the case of UV Cu^{+} CuBr lasers, a new improved analytical model consisting of the onedimensional heat conduction equation, subject to nonlinear boundary conditions of the third and fourth kind is derived. At a given air temperature of the surroundings, due to the heat convection and heat radiation, the proposed model allows to take into account the heat exchange processes between the outer surface of the laser tube and its surroundings. The gas temperature profile in the tube and the wall temperature are expressed by an explicit solution of the obtained problem and are directly dependent on the basic input laser parameters. The model is applied for evaluating the natural and forced aircooling of the laser tube.
2. Experimental Setup
The construction of the gasdischarge laser tube described by the model is presented in Figure 1. The basic tube with an 18 mm inside diameter and 24.5 mm outside diameter is made of fused quartz. The active length is 86 cm. The CuBr powder is placed in five quartz sidearm reservoirs. A ceramic tube insert with an inside diameter of 5.2 mm is sleeved in the basic tube. In the ceramic insert, five holes are made over each reservoir for the CuBr vapor diffusion into the active zone. The optic cavity of the laser studied is formed by two dielectriccoated mirrors. Mirror separation is 1.8 m.
The laser is excited by a pulsed electrical scheme with Interacting Circuits (IC). The IC excitation of CuBr lasers, operating on selfterminating copper atom transitions, was described in detail in [19].
3. Description of the Mathematical Model
The aforementioned multiline copperbromide laser operates in the UVregion [1–5]. The total input electric power is 1300 W. Taking into account the losses in power supply, the laser tube is fed with power W. The output multiline laser power is 500 mW. The geometric design of the crosssection of the laser tube in the active zone is shown in Figure 2. The laser source is manufactured from quartz (3.2), in which a ceramic tube of (3.1) is inserted along with the active laser volume, and the quartz tube is covered from the outside active volume with extra heatinsulating wadding (3.3) made of felt—glass, mineral material, or zircon oxide.
The model is developed with the following assumptions: (i) the temperature profile is determined in a quasistationary regime; (ii) the gas temperature does not change substantially in the interimpulse period; (iii) the total input electric power W in the active volume is transformed into heat, the power transferred to the walls as a result of the discharge radiation and the deactivation of the excited and charged particles is not taken into account; (iv) the thermal radiation of the heated gas in the active volume is ignored.
The temperature distribution in the crosssection of the laser tube is governed by the following quasi steadystate twodimensional heat conduction equation: where is the thermal conductivity of the gas (here neon) and is the volume density of the discharge. Due to the radial symmetry, depends only on the variable along with the radius of the tube. Consequently, in cylindrical configuration, (3.1) is reduced to the form
Usually, as it was mentioned earlier, equation (3.1) (resp., (3.2)) is solved in publications under the boundary conditions: where is the measured temperature of the outer wall of the tube (under insulation) [9–13].
Commonly, is in the form .
In the case , (3.2) possesses an exact solution [9], written in the form
For solving (3.2), we need to obtain the correct boundary conditions corresponding to (3.3) for (see Figure 2):
For that purpose, we will apply the distribution of the radial heat flow through the composite laser tube. We consider the following mixed boundary conditions of the third and fourth kind in cylindrical configuration [20–22]: where denotes the temperatures at the boundaries of the tubes, , respectively (see Figure 2).
Boundary conditions (3.6) express the equation of the continuity of the heat flow at the borders of the three mediums of the composed tube. Here is the power per unit length, ; is the active length [2, 5, 6]; , are the thermal conductivities of the tube, quartz tube and the thermal insulation, respectively; , are the diameters of the composite tubes (see Figure 2 and Table 1).

The boundary condition (3.7) shows the heat exchange between the outer surface of the laser tube and the surroundings. The first term on the righthand side of (3.7) evolves from NewtonRiemann’s law for heat exchange by convection. The second term represents the StefanBoltzmann law for heat exchange by radiation. The value of is equal to the electric power of W, in accordance with assumption (iii), as it was stated earlier, is the heat transfer coefficient, is the outside area of the insulation, is the integral emissivity of the material, c is the black body radiation coefficient, and is the temperature of the air.
The two unknown values and in boundary condition (3.7) have to be determined.
The values of the constants used in this study are given in Table 1.
In this way we obtain the temperature model described by (3.2) and boundary conditions (3.6)(3.7), equivalent to (3.2), (3.5). Our aim is to find an analytical formula for the solution of this model at .
4. Determination of the Gas Temperature at Radial Distribution of the Volume Power Density
In this section, we will obtain an explicit solution for the gas temperature satisfying the proposed theoretical model (3.2), (3.6), (3.7) and will discuss its application.
4.1. Determining the Variable Radial Distribution
Due to the lack of experimental data for , we will derive it as a qualitative theoretical dependency. From and , we have , where is the electric field intensity and is the electrical conductivity of the medium. In [24] the distribution of the field intensity in the crosssection of the tube is represented by the expression , where is a Bessel function of the first kind, . In this way, we have , where is a constant, which is found below. The Bessel function is well known and usually represented in [23, a table] in what follows. In this form, it is not suitable for direct engineeringphysics calculations. For this reason, we will approximate the term by a polynomial of the thirddegree , where . Based on tabular data from [23] and using the least squares method, we find
For , we obtain: where The constant can be found by using the equality of areas, bounded between the graphics of each of the functions and , and the abscissa (see Figure 3): After integrating in (4.4) and substituting the values of the constants, we find
In Figure 3, the distribution of the volume power densities and , according to (4.2) and (4.5), are illustrated in relative units, assuming here, in order to simplify that . From Figure 3 and (4.5), it can be observed that in the center of the discharge the local electrical volume power density for (4.2) is over two times larger than . This suggests a difference in the distribution of in the two cases being examined: and .
4.2. Determining the Gas Temperature
The solution to (3.2) at mixed boundary conditions (3.6)(3.7) and radial distribution of type (4.2) has the following form: where the constants were introduced in (4.3).
Detailed determination of (4.6) is given in the appendix.
4.3. Application of the Mathematical Model
The obtained explicit solution (4.6) can be used when the value of the temperature of the internal tube is known. There are two cases as follows.
(1)The temperature of the outside surface of the laser tube (i.e., quartz tube) under the insulation is known (see Figure 2). For existing laser devices, it can be measured, for instance, by a thermocouple. Then and can be calculated by means of the corresponding boundary conditions from (3.6).(2)The temperature is unknown. This problem can arise in the development of new laser sources or implementation of different computer simulations. In this case, the temperature of the surroundings must be specified, usually K. To use boundary condition (3.7) we need to find and . In the following section, we discuss the procedure for determining the heat transfer coefficient and obtaining a nonlinear algebraic equation for the temperature . Then, applying (3.6), we calculate , , and .5. Evaluation of Cooling and Discussion
We will apply the derived temperature model (3.2), (3.6)(3.7) for determining the gas temperature in the cases of free and forced convection.
5.1. Cooling of the Laser Tube by Free Convection
In [18], a simplified temperature model in the case of free convection at has been used. Here we will compare the results obtained by our new model for the general case with those in [18].
In the case of free convection, the heat transfer coefficient in (3.7) was calculated in [18] as
The substitution of (5.1) in boundary condition (3.7) with consequent representation in relation to the power per unit length results in [18]
In the previous expressions (5.1)(5.2), the numerical values of , , , and are given in Table 1. The data is correct for air temperature K [22].
However (5.2) is a nonlinear equation with respect to the outside temperature of the laser device insulation—. Also (5.2) can be easily solved by any computer algebra system, for instance by Mathematica [25]. Once the temperature is calculated, the values of , , and can be evaluated from (3.6), and the gas temperature in the internal tube is determined by (4.6).
In Figure 4, on the same coordinate system, the distributions of the gas temperature in the crosssection along with the radius of the laser tube for and are presented. In Table 2, special characteristic temperatures , , , , and the maximum temperature in the center of the laser tube are given (see also Figure 1).

Table 2 shows that temperatures , , , and are equal. Their values are determined by the total electrical power emitted within the active volume and are independent from its radial distribution. In both cases, this power is the same—1000 W. Table 2 and Figure 4 show that when is 90°C higher than the corresponding value when . The results for the gas temperature in the case of have the same behavior as the calculated values in [9, 11].
5.2. Cooling of the Laser Tube by Forced Convection
For all types of convection, the Nusselt criterion, , holds [20–22], from which for and , we find
In the case of forced convection, the Reynolds criterion is represented by [20, 21] where is the velocity of the moving fluid, is the length of the laser tube and is the kinematical viscosity of the air. However (5.4) is valid for [22].
For horizontal tubes with forced air cooling the following equality holds [22]:
From (5.3), (5.4), and (5.5), the heat transfer coefficient is Substituting (5.6) in boundary condition (3.7) and representing it with respect to the power per unit length , we obtain a nonlinear algebraic equation for in the form In this way by solving (5.7), we determine . Then using (3.6) and (4.6), we find the required gas temperature profile in the crosssection of the laser tube.
The obtained values of some characteristic temperatures are given in Table 3, including the maximum value along the center of the tube. The results of the calculated values of in the two cases and are shown in Figure 5 for air flow m/s.

As it is expected the cooling process causes a decrease of the buffer gas temperature in relation with the case of free convection (compare Tables 2 and 3, and Figures 4 and 5).
It can be noted that although the maximum local electric power at the center of the tube is twice higher for (see Figure 3), the difference between the corresponding maximum temperatures is only 95°C. This result is almost the same as in the case of free convection. The deviation is on average around 6%. We can conclude that in principle, solution (3.4) can successfully be used to analyze temperature conditions of existing laser sources, when the temperature is known.
As an absolute quantity, the difference of °C at the center of the discharge should not be overlooked. The model discussed in this paper can better explain and predict the occurrence of a number of negative phenomena connected with the overheating of the laser medium. The increase of °C in the temperature at the center of the discharge can lead to a contraction of the gas discharge, thermal ionizing instability, and thermochemical gas degradation, additional thermal population of lower laser levels. In the end, this leads to decreased laser power and deterioration in mode composition. In some cases, the overheating of the discharge at the center of the tube may lead to a cessation of laser generation along its axis and the appearance of dark spots at the center of laser beam. For this reason, regardless the complexity of the new model, its use is fully advisable.
6. Average Temperature in the Active Volume
During the analysis of the temperature condition of existing or new laser devices, the average gas temperature in the active volume is a characteristic of great importance. It is defined as where is the radial distribution of temperature in the active volume.
The average value of temperature for equal configurations depends only on the electric power, supplied to the active volume, and is independent of the radial distribution of . In our case, the electrical power is W. The average temperature should not change for the temperature distribution of the type in (3.4) and (4.6).
All subsequent calculations have been made using the Mathematica software system [25].
For the radial distribution of given by formula (3.4), the result is The numerical value of (6.2) is K.
For the radial distribution of from (4.6) the quantity (6.1) does not have an exact algebraic solution and the integral in (6.1) is solved numerically. The numerical value is K. The approximate calculation leads to an insignificant deviation of values for and , estimated by a relative error of 0.6%.
The presence of a hypergeometric function of type in the solution of (6.2) makes it practically difficult to use. For this reason, we represent the function by its Gauss series [25, 26], limited only to the first two terms of the expansion The result is
In this way by using (6.2), we obtain
The average value of the temperature from (6.5) is K, with a relative error of 1%. This shows that (6.5) can be used with sufficient accuracy to determine the average temperature in the active volume: .
7. Conclusion
A theoretical mathematical model for evaluating the buffer gas temperature in UV Cu^{+} NeCuBr laser is developed. It takes into account the nonuniformly distributed electrical power along the crosssection of the laser tube. Based on common theoretical dependencies, a suggestion is made for the qualitative distribution of volume power density . The model includes a heat conduction equation subject to nonlinear boundary conditions. An explicit solution with these conditions is obtained. The model is applied in the cases of free and forced convection. A simple expression is established for calculating the average gas temperature in the active volume.
An evaluation of the previous existing solution has been presented, describing the distribution of the gas temperature under the assumption . It has been established that, for such an assumption, the error when determining at the center of the tube is about 6%.
The obtained results at have been compared with similar calculated results by using simple mathematical gas temperature models [9, 11].
A comparison has been made between the obtained temperature profiles of at and . It has been established that at the center of the tube the temperature when is about 90–95°C higher, both for free and forced convection.
It has to be noted that the simple model when cannot be used to evaluate radial buffer gas temperature when the temperature of the wall is unknown. The advantage of the model presented here is that when specific geometric dimensions have been chosen, the temperature of the surroundings and other parameters, the temperature of the outer wall of the laser tube is calculated, after which the values for the gas temperature within the tube are calculated. Therefore, the new model can be applied not only to precise the existing models but to design new laser devices as well.
The proposed theoretical model and methodology could be useful for solving similar engineering problems where the calculation of the radial heat flow with an internal heat source is required. It gives the opportunity to carry out further computer simulations in order to optimize laser generation, while changing the geometrical design, tube materials, heat insulation, input electric power, and laser operating conditions.
Appendix
Let us consider the problem (3.2), (3.5), assuming that the temperature is determined. For , and , (3.2) becomes
By the change of variable we have . After some simplification in (A.1) and (3.3), we obtain the problem
By introducing a new variable However (A.3) can be rearranged in the form
The integration of (A.6) and substitution of by the expression (4.2) yields , where is a constant of integration. We find
Coming back to the variable according to (A.5), we obtain
The second boundary condition in (A.4) requires that , so
The integration of (A.9) yields where is a constant. Using the boundary condition (A.4), we find
The substitution of into (A.10) gives
Finally, by using (A.2), we obtain the required solution of (3.2) in the form (4.6):
Acknowledgments
This work is partially supported by the National Scientific Fund of Bulgarian Ministry of Education and Science, Grant no. VUMI205/2006.
References
 N. K. Vuchkov, K. A. Temelkov, and N. V. Sabotinov, “UV lasing on ${\text{Cu}}^{\text{+}}$ in NeCuBr pulse longitudinal discharge,” IEEE Journal of Quantum Electronics, vol. 35, pp. 1799–1804, 1999. View at: Publisher Site  Google Scholar
 N. K. Vuchkov, K. A. Temelkov, P. V. Zahariev, and N. V. Sabotinov, “Influence of the active zone diameter on the UV ion NeCuBr laser performance,” IEEE Journal of Quantum Electronics, vol. 37, pp. 1538–1546, 2001. View at: Publisher Site  Google Scholar
 K. Beev, K. Temelkov, N. Vuchkov et al., “Optical properties of polymer films for near UV recording,” Journal of Optoelectronics and Advanced Materials, vol. 7, pp. 1315–1318, 2005. View at: Google Scholar
 M. Ilieva, V. Tsakova, N. K. Vuchkov, K. A. Temelkov, and N. V. Sabotinov, “UV copper ion laser treatment of poly3,4 ethylenedioxythiophene,” Journal of Optoelectronics and Advanced Materials, vol. 9, pp. 303–306, 2007. View at: Google Scholar
 N. K. Vuchkov, K. A. Temelkov, and N. V. Sabotinov, “Effect of hydrogen on the average output power of the UV ${\text{Cu}}^{\text{+}}$ NeCuBr laser,” IEEE Journal of Quantum Electronics, vol. 41, pp. 62–65, 2005. View at: Publisher Site  Google Scholar
 N. Vuchkov, “High discharge tube resource of the UV ${\text{Cu}}^{\text{+}}$ NeCuBr laser and some applications,” in New Developments in Lasers & ElectroOptics Research, W. T. Arkin, Ed., pp. 41–47, Nova Science, Lancaster, UK, 2006. View at: Google Scholar
 Yu. G. Gradoboev, Yu. V. Gulyaev, M. A. Kazaryan et al., “UV radiation source based on a copper vapour laser with acoustooptically controlled spectral and temporal parameters,” Quantum Electronics, vol. 34, pp. 1133–1137, 2004. View at: Publisher Site  Google Scholar
 Y. P. Raizer, M. N. Shneide, and N. A. Yatsenko, RadioFrequency Capacitive Discharges, CRC Press, New York, NY, USA, 1995.
 M. J. Kushner and B. E. Warner, “Largebore coppervapor lasers: kinetics and scaling issues,” Applied Physics, vol. 54, pp. 2970–2982, 1983. View at: Publisher Site  Google Scholar
 A. M. Boichenko, G. S. Evtushenko, O. V. Zhdaneev, and S. I. Yakovlenko, “Theoretical analysis of the mechanisms of influence of Hydrogen additions on the emission parameters of a copper vapour laser,” Quantum Electronics, vol. 33, pp. 1047–1058, 2003. View at: Publisher Site  Google Scholar
 B. A. Ghani and M. Hammadi, “Modeling the plasma kinetics mechanisms of CuBr laser with neon–hydrogen additives,” Optics & Laser Technology, vol. 38, no. 2, pp. 67–76, 2006. View at: Google Scholar
 R. J. Carman, D. J. W. Brown, and J. A. Piper, “A selfconsistent model for the discharge kinetics in a highrepetitionrate coppervapor laser,” IEEE Journal of Quantum Electronics, vol. 30, no. 8, pp. 1876–1895, 1994. View at: Publisher Site  Google Scholar
 R. J. Carman, R. P. Mildren, M. J. Withford, D. J. W. Brown, and J. A. Piper, “Modeling the plasma kinetics in a kinetically enhanced copper vapor laser utilizing HCl + ${\text{H}}_{2}$ admixtures,” IEEE Journal of Quantum Electronics, vol. 36, no. 4, pp. 438–449, 2000. View at: Publisher Site  Google Scholar
 B. L. Pan, G. Chen, B. N. Mao, and Z. X. Yao, “Kinetic process of UV ${\text{Cu}}^{\text{+}}$ laser in NeCuBr longitudinal pulsed discharge,” Optics Express, vol. 14, no. 19, pp. 8644–8653, 2006. View at: Publisher Site  Google Scholar
 D. N. Astadjov, N. K. Vuchkov, and N. V. Sabotinov, “Parametric study of the CuBr laser with hidrogen additives,” IEEE Journal of Quantum Electronics, vol. 24, pp. 1926–1935, 1988. View at: Publisher Site  Google Scholar
 K. A. Temelkov, N. K. Vuchkov, B. L. Pan, N. V. Sabotinov, B. Ivanov, and L. Lyutov, “Strontium bromide vapor laser excited by nanosecond pulsed longitudinal discharge,” in 14th International School on Quantum Electronics: Laser Physics and Applications, vol. 6604 of Proceedings of SPIE, pp. 1–5, March 2007. View at: Google Scholar
 I. P. Iliev, S. G. GochevaIlieva, and N. V. Sabotinov, “Analytic study of the temperature profile in a copper bromide laser,” Quantum Electronics, vol. 38, pp. 338–342, 2008. View at: Publisher Site  Google Scholar
 S. G. GochevaIlieva, I. P. Iliev, K. A. Temelkov, N. K. Vuchkov, and N. V. Sabotinov, “Analytical model of the temperature in UV ${\text{Cu}}^{\text{+}}$ CuBr laser,” in Applications of Mathematics in Engineering and Economics, M. D. Todorov, Ed., vol. CP1067, pp. 114–121, American Institute of Physics, Melville, NJ, USA, 2008. View at: Google Scholar
 N. K. Vuchkov, D. N. Astadjov, and N. V. Sabotinov, “Influence of the excitation circuits on the CuBr laser performance,” IEEE Journal of Quantum Electronics, vol. 30, no. 3, pp. 750–758, 1994. View at: Publisher Site  Google Scholar
 M. N. Özişik, Heat Transfer. A Basic Approach, McGrawHill, Boston, Mass, USA, 1985.
 V. A. Milchev, D. Z. Uzunov, V. Y. Jordanov, and D. K. Palov, Heating Technology, Tehnika, Sofia, Bulgaria, 1989.
 M. Oprev, Tz. Batov, and D. Uzunov, Heating Technology, Tehnika, Sofia, Bulgaria, 1978.
 M. Abramowitz and I. Stegun, Handbook of Mathematical Functions, with Formulas, Graphs, and Mathematical Tables, Dover, New York, NY, USA, 9th edition, 1964. View at: MathSciNet
 P. Blau, “Impedance matching and electric field penetration in metal vapour lasers,” in Pulsed Metal Vapour Lasers, E. Little and N. V. Sabotinov, Eds., vol. 5 of Nato Science Partnership SubSeries: 1, pp. 215–220, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1996. View at: Google Scholar
 http://reference.wolfram.com/mathematica/guide/Mathematica.html.
 L. J. Slater, Generalized Hypergeometric Functions, Cambridge University Press, Cambridge, UK, 1966. View at: MathSciNet
Copyright
Copyright © 2009 Iliycho Petkov Iliev et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.