Research Article  Open Access
Analysis of Radiative Radial Fin with TemperatureDependent Thermal Conductivity Using Nonlinear Differential Transformation Methods
Abstract
Radiative radial fin with temperaturedependent thermal conductivity is analyzed. The calculations are carried out by using differential transformation method (DTM), which is a seminumericalanalytical solution technique that can be applied to various types of differential equations, as well as the Boubaker polynomials expansion scheme (BPES). By using DTM, the nonlinear constrained governing equations are reduced to recurrence relations and related boundary conditions are transformed into a set of algebraic equations. The principle of differential transformation is briefly introduced and then applied to the aforementioned equations. Solutions are subsequently obtained by a process of inverse transformation. The current results are then compared with previously obtained results using variational iteration method (VIM), Adomian decomposition method (ADM), homotopy analysis method (HAM), and numerical solution (NS) in order to verify the accuracy of the proposed method. The findings reveal that both BPES and DTM can achieve suitable results in predicting the solution of such problems. After these verifications, we analyze fin efficiency and the effects of some physically applicable parameters in this problem such as radiationconduction fin parameter, radiation sink temperature, heat generation, and thermal conductivity parameters.
1. Introduction
Extended surfaces are extensively used in various industrial applications. An extensive review on this topic is presented by Kraus et al. [1]. Fins are very frequently encountered in many engineering applications to enhance heat transfer. Numerous contributions have been made in the heat transfer analysis of the fins. Constant thermophysical properties and uniform heat transfer coefficient are often assumed in the determination of the temperature distribution along an extended surface. The mathematical complexity of the conservation energy equation is reduced by this assumption and therefore wellestablished closed form analytical solutions can be obtained for a number of cases. If a large temperature difference exists within a fin, the thermal conductivity may not be constant. Furthermore, in general, the heat transfer coefficient may vary along a fin. The heat transfer coefficient may be a function of the spatial coordinate only along a fin or may depend on the local temperature difference between the fin surface and the surrounding fluid. Kundu [2] analytically carried out the thermal analysis and optimization of longitudinal and pin fins of uniform thickness subject to fully wet, partially wet, and fully dry surface conditions. Moreover, Kundu [3] analytically analyzed the performance and optimization of longitudinal and pin fins of step reduction in local crosssection (SRC) profile subject to combined heat and mass transfer. Coşkun and Atay [4, 5] used variational iteration method to analyze convective straight and radial fins with temperaturedependent thermal conductivity. Sharqawy and Zubair [6] carried out an analysis to study the efficiency of straight fins with different configurations when subjected to simultaneous heat and mass transfer mechanisms. Domairry and Fazeli [7] solved nonlinear straight fin differential equations by homotopy analysis method (HAM) to evaluate the temperature distribution within the fin. Arslanturk [8] obtained correlation equations for optimum design of annular fins with temperaturedependent thermal conductivity. Kulkarni and Joglekar [9] proposed and implemented a numerical technique based on residue minimization to solve the nonlinear differential equation, which governs the temperature distribution in straight convective fins having temperaturedependent thermal conductivity. Khani et al. [10] used HAM to evaluate the analytical approximate solutions and efficiency of the nonlinear fin problem with temperaturedependent thermal conductivity and heat transfer coefficient. Kundu [11] described an analytical method for temperature and heat transfer characteristics of an annular step fin (ASF) with the simultaneous heat and mass transfer mechanisms. Bouaziz and Aziz [12] introduced a new concept called the double optimal linearization method (DOLM) to derive simple and accurate expressions for predicting the thermal performance of a convectiveradiative fin with temperaturedependent thermal conductivity. Khani and Aziz [13] used HAM to develop analytical solutions for the thermal performance of a straight fin of trapezoidal profile when both thermal conductivity and heat transfer coefficient are temperaturedependent.
The differential transformation method (DTM) is a seminumericalanalytical method. DTM, which is based on the Taylor series expansion, was first proposed by Zhou [14] in 1986 for the solution of linear and nonlinear initial value problems that appear in electrical circuits. This method obtains a solution in the form of a polynomial. Indeed, this fact can be seen in Section 4, where the concept of differential transforms is briefly described. It may be pointed out that, later, this method has been successfully used in a series of literature [15–21] dealing with many engineering problems.
There are recent studies on the application of DTM to the heat transfer problems in the literature. Chu and Chen [22] applied a hybrid method of differential transformation and finite difference method to solve a transient heat conduction problem which had complex nonlinear terms. Chu and Lo [23] applied the differential transformation technique for transforming and discretizing the governing equations as well as the boundary conditions and provided two numerical examples. Lo and Chen [24] proposed an alternative numerical method to investigate hyperbolic heat conduction problems using the hybrid differential transfer/controlvolume method. Joneidi et al. [25] provided analytical solution to fin efficiency of convective straight fins with temperaturedependent thermal conductivity by DTM. Jang et al. [26] investigated and characterized a twodimensional thermal conductive boundary value problem with discontinuous boundary and initial conditions. Rashidi et al. [27] applied DTM to find the analytic solution for the problem of mixed convection about an inclined flat plate embedded in a porous medium. Yaghoobi and Torabi [28] applied DTM to solve the problems of convective and convectiveradiative cooling of a lumped system with temperaturedependent specific heat.
More recently, Ganji et al. [29] applied DTM to the problem of convectiveradiative straight fins with temperaturedependent thermal conductivity. But they considered zero dimensionless convective and radiative sink temperatures and without heat generation.
This paper is an analytical study of the thermal performance of a radiative fin with variation of thermal conductivity with temperature. The problem considers nonzero dimensionless radiation sink temperature, , and dimensionless heat generation, , which is the novelty of the present work. The resulting nonlinear differential equation is solved by DTM to evaluate the temperature distribution within the fin. Accordingly, the appropriate convergence study and comparison with previously published related articles, the results obtained using variational iteration method (VIM) [4, 30], Adomian decomposition method (ADM) [31], homotopy analysis method (HAM) [32], and numerical solution (NS), were employed in order to verify the accuracy of the proposed method. Using the temperature distribution, we express the efficiency of the fin in terms of radiationconduction fin parameter, , and thermal conductivity parameter, . Because a broad range of governing parameters are investigated, the results should be useful in a number of engineering applications.
2. Description of the Problem
A typical heat pipe/fin space radiator is shown in Figure 1. Both surfaces of the fin are radiating to the outer space at a very low temperature, which is assumed equal to zero absolute. The fin has temperaturedependent thermal conductivity, , which depends on temperature linearly, and fin is diffusegrey with emissivity . The tube surfaces’ temperature and the base temperature of the fin are constant, and the convective exchange between the fin and the heat pipe is neglected. The temperature distribution within the fin is assumed to be onedimensional, because the fin is assumed to be thin. Hence, only fin tip length is considered as the computational domain.
The energy balance equation for a differential element of the fin [2–5] is given by where and are thermal conductivity and the StefanBoltzmann constant, respectively.
The thermal conductivity of the fin material is assumed to be a linear function of temperature [4–7] according to where is the thermal conductivity at the temperature of the fin and is the measure of variation of the thermal conductivity with temperature.
We introduce the following dimensionless parameters:
The formulation of the fin problem [8–11] reduces to the following equation: with the following boundary conditions:
3. Fin Efficiency
The heat transfer rate from the surfaces of a fin is found by applying the StefanBoltzmann law, namely [31], Fin efficiency is defined as the ratio of energy radiated away by the fin to the energy that would be radiated if the entire fin was at the base temperature [33]:
Employing the dimensionless parameters in (3), fin efficiency is expressed as
4. Fundamentals of Differential Transformation Method
Let be analytic in a domain , and let represent any point in . The function is then represented by a power series whose center is located at . A Taylor series expansion function of about takes the form The particular case of (9) when is the Maclaurin series of and is expressed as As explained in [34] the differential transformation of the function is defined as follows: where is the original function and is the transformed function (commonly referred to as the Tfunction). The differential spectrum of is confined within the interval , where is a constant. The differential inverse transform of is defined as follows:
It is clear that the concept of differential transformation is based upon the Taylor series expansion. Values of the function are referred to as discretes; that is, is known as the zero discrete, is the first discrete, and is the th discrete. The more discretes are available, the more precise it is possible to restore the unknown function. The function consists of function , and its value is given by the sum of the function with as its coefficient. In real applications, with the right choice of constant , the larger the values of argument are, the more rapid the discretes of spectrum are reduced. The function is expressed by a finite series and (12) can be written as where is total number of polynomial terms used in the DTM. Mathematical operations performed by differential transform method are listed in Table 1.

5. Solution Protocols
5.1. Solution Using the Differential Transformation Method DTM
Now we apply differential transformation method to (4). Taking the differential transform of (4) with respect to , and considering according to Table 1, we gives
From boundary condition in (5a), and employing the second formulation from Table 1, we obtain Supposing that and using (14)–(16) one can obtain as follows:where is constant, and we will calculate it with considering another boundary condition in (5b) in point .
The above process is continuous. Substituting (17a), (17b), (17c), and (17d) into the main equation based on DTM, it can be obtained that the closed form of the solutions is To obtain the value of , we substitute the boundary condition from (5b) into (18) in point . So, we have Solving (19) gives the value of . The resultant equation can be solved by using NewtonRaphson iterative technique for determination of unknown tip temperature . We employed the Maple’s builtin Roots command which numerically approximates the roots of an algebraic function using the specified method, such as NewtonRaphson, bisection, secant, and fixedpoint iteration methods, and returns the specified outputs. This command uses the NewtonRaphson method by default. Also, the estimated initial value was picked as a starting number near the solution using the graph to locate sufficiently close number to the root.
As an example, let us assume that , , , , and . Therefore, the value of , applying which will be used in this paper, will be obtained as follows: Substituting this obtained parameter in (18), the temperature profile of fin for this special case will be as follows: The VIM solution of (4), with , is obtained in the following form [4]: And the ADM solution of (4), with , is obtained in the following form [31]:
5.2. Solution Using the Boubaker Polynomials Expansion Scheme (BPES)
The Boubaker polynomials expansion scheme (BPES) [34–52] is a resolution protocol which has been successfully applied to several appliedphysics and mathematics problems. The BPES protocol ensures the validity of the related boundary conditions regardless of main equation features. The Boubaker polynomials expansion scheme (BPES) is based on the Boubaker polynomials first derivatives properties: with .
Several solutions have been proposed through the BPES in many fields such as numerical analysis [34–37], theoretical physics [38–41], mathematical algorithms [42], heat transfer [43], homodynamic [44, 45], material characterization [46], fuzzy systems modeling [47–50], and biology [51, 52].
In reference to (4), the BPES is applied to the system through setting the expression where are the order Boubaker polynomials, are minimal positive roots, is a prefixed integer, and are unknown pondering real. The BPES solution is obtained through the following steps. (i) Integrate, for a given value of , the whole expressions given in the left side of (4) along the given domain. By introducing expression (26) in (25), boundary conditions become redundant since they already verified by the proposed expansion (26). Consecutively, and by majoring the quadrature terms, the problem is transformed in a linear system with unknown real variables . Consequently, it comes for (4) and (25)(26) that with the matrix standard form The system (28) is hence reduced to approximately arithmetical operations and solved using the Householder algorithm [53, 54]. (ii) Incremente . (iii) Test the convergence of the coefficients until stability of the solution.
The final result is hence (obtained for )
6. Results and Discussion
In this section, we divide our study into three subsections. Firstly, the convergence rate of the DTM and BPES is checked. Secondly, comparison with previously published related articles is employed in order to verify the accuracy of the proposed method. Finally, some figures are introduced regarding the effects of physically applicable parameters such as thermal conductivity parameter, , dimensionless temperature whose is constant, , dimensionless radiation sink temperature, , radiationconduction parameter, , and dimensionless heat generation parameter, , on the temperature distribution within the fin. Moreover, for all numerical results which are reported here, the following values of variables are used unless otherwise indicated by graphs or the table:
6.1. Convergence Study
Opposite to the convergence of the Boubaker Polynomials Expansion Scheme (BPES), which is intrinsically tested (via the value of , Section 5.2), the proposed differential transformation method DTM adopts an iterative procedure to obtain the highorder Taylor series. Since the Taylor series is an infinite series, the differential transformation should theoretically comprise an infinite series. However, the present results indicate that in practice a small value (i.e., ) is sufficient to provide an accurate solution.
Figure 2 depicts the convergence of the dimensionless temperature for six different cases. From Figure 2 it is clearly visible that for two cases of and when , more than 10 terms are needed to obtain the value of the . Also, it is seen that if and are chosen to plot the temperature within the fin, at least 30 terms are required to obtain an accurate value of the . As more terms are taken, the converges to its exact value. Therefore, the numerical results from the DTM approach which are presented in the next sections were obtained by taking sufficient terms to the temperature solution. It should be noted that the final results from DTM are in good agreement with the RungeKuttaFehlberg method which is a welltested numerical solution. Moreover, since the problem is a highly nonlinear radiative equation, an exact analytical solution cannot be found for (4).
(a)
(b)
6.2. Validation of the Results
To check the accuracy of the present solution, by considering our problem is converted to a simpler case which was studied in four pioneering works [4, 30–32] in this field. The results obtained from DTM analysis are compared with VIM [4, 30] and ADM [31] in Figure 3. As it can be seen, this method leads to acceptable results compared with those methods.
For the case of different values for thermal conductivity and radiationconduction fin parameter, results of the present analysis are tabulated against the ADM, VIM, and again with numerical solution (NS), of the fourthfifth order RungeKuttaFehlberg method using the Maple package, in Table 2. A very interesting agreement between the results is observed, which confirms the validity of the DTM.

In Figure 4, the efficiency of the DTM and BPES solution technique for can be illustrated with comparison respect to HAM [32] and NS. Also by means of these comparisons, it can be shown that DTM and BPES are a better alternatives in the solution of such problems.
6.3. Results of Present Study
Figure 5 shows the behavior of fin tip temperature, , relative to the thermal conductivity parameters, , and the radiationconduction fin parameter, . Figure 5 clearly demonstrates that increasing in the values of thermal conductivity parameter produces increase in values of fin tip temperature.
Figure 6 illustrates the dimensionless temperature distributions along the fin surface with varying from −0.3 to 0.3. The curve marked represents the case when the thermal conductivity is a constant and its value is . The curves with correspond to fin materials whose thermal conductivity increases as temperature increases. The converse is true of curves with . As the parameter increases, the average thermal conductivity of the material increases, and as expected, the result is a gradual increase in the local temperature.
In Figure 7 we illustrate the effect of the radiationconduction parameter, , on the temperature distribution in the fin. As the radiative transport becomes stronger, the radiative cooling becomes more effective, which in turn causes the lowering of temperatures in the fin.
Figure 8 shows that as the increases, the temperature within the fin decreases. This is because of the fact that, as becomes lower, the becomes more sensitive to the temperature. Moreover, Figure 9 displays the effect of the radiation sink temperature, , on the temperature distribution in the fin. As the radiation sink temperature increases, the radiative heat loss decreases, resulting in higher temperatures in the material. Figure 10 allows us to assess the effect of the dimensionless heat generation in (4), which varies from 0.2 to 0.6. As expected, an increase in the value of causes an increase in the value of within the fin.
Figure 11 shows the behavior of the fin efficiency relative to and . Figure 11 clearly demonstrates that increasing in the values of radiationconduction parameter produces a decrease in the value of the fin efficiency. In addition, the fin efficiency increases as the thermal conductivity parameter, , increases. To explain the effect of parameter , we note that while the temperature increases as increases (see Figure 6), the temperature of the fin becomes closer to , and therefore, from (6), increases resulting in an increase in the value of . For all results reported in Figures 5–11 the DTM results were checked against NS which were in excellent agreement with each other. It would be useful to mention that DTM was tested for moving convectiveradiative fins [55] against numerical solution, and results of both methods were in good agreement with each other.
7. Conclusion
In this study, the DTM and BPES have been utilized to derive approximate explicit analytical solution for nonlinear radiative radial fin heat transfer problem with temperaturedependent thermal conductivity and heat generation. The figures and table clearly show that these methods provide excellent approximations to the solution of this nonlinear equation with high accuracy. As the radiationconduction parameter increases the effect is to lower the fin temperature. As expected, it was observed that increasing in the value of or increases the temperature distribution through the fin. The collection of temperature graphs should be useful in the study and design of a variety of engineering systems where fins are attached to realize heat transfer enhancement.
Nomenclature
:  Fin tip length, m 
:  Domain 
:  Constant 
:  Temperaturedependent thermal conductivity, 
:  Thermal conductivity at the base temperature, 
:  Volumetric heat generation, 
:  Dimensionless volumetric heat generation 
:  Heat transfer rate from the surfaces of a fin 
:  Temperature, 
:  Fin’s base temperature, 
:  Radiation sink temperature, 
:  Transformed analytical function 
:  Original analytical function 
:  Semithickness of the fin, m. 
:  Thermal conductivity parameter 
:  Emissivity 
:  Fin efficiency 
:  Slope of the thermal conductivitytemperature curve, 
:  StefanBoltzmann constant, 
:  Dimensionless temperature 
:  Dimensionless radiation sink temperature 
:  Radiationconduction fin parameter. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 A. D. Kraus, A. Aziz, and J. R. Welty, Extended Surface Heat Transfer, John Wiley, New York, NY, USA, 2002.
 B. Kundu, “Performance and optimum design analysis of longitudinal and pin fins with simultaneous heat and mass transfer: unified and comparative investigations,” Applied Thermal Engineering, vol. 27, no. 56, pp. 976–987, 2007. View at: Publisher Site  Google Scholar
 B. Kundu, “Performance and optimization analysis of SRC profile fins subject to simultaneous heat and mass transfer,” International Journal of Heat and Mass Transfer, vol. 50, no. 78, pp. 1545–1558, 2007. View at: Publisher Site  Google Scholar
 S. B. Coşkun and M. T. Atay, “Analysis of convective straight and radial fins with temperaturedependent thermal conductivity using variational iteration method with comparison with respect to finite element analysis,” Mathematical Problems in Engineering, vol. 2007, Article ID 42072, 15 pages, 2007. View at: Publisher Site  Google Scholar
 S. B. Coşkun and M. T. Atay, “Fin efficiency analysis of convective straight fins with temperature dependent thermal conductivity using variational iteration method,” Applied Thermal Engineering, vol. 28, no. 1718, pp. 2345–2352, 2008. View at: Publisher Site  Google Scholar
 M. H. Sharqawy and S. M. Zubair, “Efficiency and optimization of straight fins with combined heat and mass transferan analytical solution,” Applied Thermal Engineering, vol. 28, no. 1718, pp. 2279–2288, 2008. View at: Publisher Site  Google Scholar
 G. Domairry and M. Fazeli, “Homotopy analysis method to determine the fin efficiency of convective straight fins with temperaturedependent thermal conductivity,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 2, pp. 489–499, 2009. View at: Publisher Site  Google Scholar
 C. Arslanturk, “Correlation equations for optimum design of annular fins with temperature dependent thermal conductivity,” Heat and Mass Transfer, vol. 45, no. 4, pp. 519–525, 2009. View at: Publisher Site  Google Scholar
 D. B. Kulkarni and M. M. Joglekar, “Residue minimization technique to analyze the efficiency of convective straight fins having temperaturedependent thermal conductivity,” Applied Mathematics and Computation, vol. 215, no. 6, pp. 2184–2191, 2009. View at: Publisher Site  Google Scholar
 F. Khani, M. A. Raji, and H. H. Nejad, “Analytical solutions and efficiency of the nonlinear fin problem with temperaturedependent thermal conductivity and heat transfer coefficient,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 8, pp. 3327–3338, 2009. View at: Publisher Site  Google Scholar
 B. Kundu, “Analysis of thermal performance and optimization of concentric circular fins under dehumidifying conditions,” International Journal of Heat and Mass Transfer, vol. 52, no. 1112, pp. 2646–2659, 2009. View at: Publisher Site  Google Scholar
 M. N. Bouaziz and A. Aziz, “Simple and accurate solution for convectiveradiative fin with temperature dependent thermal conductivity using double optimal linearization,” Energy Conversion and Management, vol. 51, no. 12, pp. 2776–2782, 2010. View at: Publisher Site  Google Scholar
 F. Khani and A. Aziz, “Thermal analysis of a longitudinal trapezoidal fin with temperaturedependent thermal conductivity and heat transfer coefficient,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 3, pp. 590–601, 2010. View at: Publisher Site  Google Scholar
 J. K. Zhou, Differential Transform and Its Applications for Electrical CircuIts, Huarjung University Press, Wuhan, China, 1986.
 I. H. AbdelHalim Hassan, “Application to differential transformation method for solving systems of differential equations,” Applied Mathematical Modelling, vol. 32, no. 12, pp. 2552–2559, 2008. View at: Publisher Site  Google Scholar
 A. S. V. Ravi Kanth and K. Aruna, “Solution of singular twopoint boundary value problems using differential transformation method,” Physics Letters A, vol. 372, no. 26, pp. 4671–4673, 2008. View at: Publisher Site  Google Scholar
 H. S. Yalcin, A. Arikoglu, and I. Ozkol, “Free vibration analysis of circular plates by differential transformation method,” Applied Mathematics and Computation, vol. 212, no. 2, pp. 377–386, 2009. View at: Publisher Site  Google Scholar
 B. Jang, “Solving linear and nonlinear initial value problems by the projected differential transform method,” Computer Physics Communications, vol. 181, no. 5, pp. 848–854, 2010. View at: Publisher Site  Google Scholar
 H. Yaghoobi, P. Khoshnevisrad, and A. Fereidoon, “Application of the differential transformation method to a modified van der pol oscillator,” Nonlinear Science Letter A, vol. 2, pp. 171–180, 2011. View at: Google Scholar
 H. Yaghoobi and M. Torabi, “Analytical solution for settling of nonspherical particles in incompressible Newtonian media,” Powder Technology, vol. 221, pp. 453–463, 2012. View at: Publisher Site  Google Scholar
 M. Torabi and H. Yaghoobi, “Accurate solution for acceleration motion of a vertically falling spherical particle in incompressible newtonian media,” Canadian Journal of Chemical Engineering, 2012. View at: Publisher Site  Google Scholar
 H. P. Chu and C. L. Chen, “Hybrid differential transform and finite difference method to solve the nonlinear heat conduction problem,” Communications in Nonlinear Science and Numerical Simulation, vol. 13, no. 8, pp. 1605–1614, 2008. View at: Publisher Site  Google Scholar
 H. P. Chu and C. Y. Lo, “Application of the hybrid differential transformfinite difference method to nonlinear transient heat conduction problems,” Numerical Heat Transfer A, vol. 53, no. 3, pp. 295–307, 2008. View at: Publisher Site  Google Scholar
 C. Y. Lo and B. Y. Chen, “Application of hybrid differential transform/controlvolume method to hyperbolic heat conduction problems,” Numerical Heat Transfer B, vol. 55, no. 3, pp. 219–231, 2009. View at: Publisher Site  Google Scholar
 A. A. Joneidi, D. D. Ganji, and M. Babaelahi, “Differential transformation method to determine fin efficiency of convective straight fins with temperature dependent thermal conductivity,” International Communications in Heat and Mass Transfer, vol. 36, no. 7, pp. 757–762, 2009. View at: Publisher Site  Google Scholar
 M. Jang, Y. Yeh, C. Chen, and W. Yeh, “Differential transformation approach to thermal conductive problems with discontinuous boundary condition,” Applied Mathematics and Computation, vol. 216, no. 8, pp. 2339–2350, 2010. View at: Publisher Site  Google Scholar
 M. M. Rashidi, N. Laraqi, and S. M. Sadri, “A novel analytical solution of mixed convection about an inclined flat plate embedded in a porous medium using the DTMPadé,” International Journal of Thermal Sciences, vol. 49, no. 12, pp. 2405–2412, 2010. View at: Publisher Site  Google Scholar
 H. Yaghoobi and M. Torabi, “The application of differential transformation method to nonlinear equations arising in heat transfer,” International Communications in Heat and Mass Transfer, vol. 38, no. 6, pp. 815–820, 2011. View at: Publisher Site  Google Scholar
 D. D. Ganji, M. Rahimi, M. Rahgoshay, and M. Jafari, “Analytical and numerical investigation of fin efficiency and temperature distribution of conductive, convective, and radiative straight fins,” Heat Transfer, vol. 40, no. 3, pp. 233–245, 2011. View at: Publisher Site  Google Scholar
 H. Tari, D. D. Ganji, and H. Babazadeh, “The application of He's variational iteration method to nonlinear equations arising in heat transfer,” Physics Letters A, vol. 363, no. 3, pp. 213–217, 2007. View at: Publisher Site  Google Scholar
 C. Arslanturk, “Optimum design of space radiators with temperaturedependent thermal conductivity,” Applied Thermal Engineering, vol. 26, no. 1112, pp. 1149–1157, 2006. View at: Publisher Site  Google Scholar
 S. Abbasbandy, “The application of homotopy analysis method to nonlinear equations arising in heat transfer,” Physics Letters A, vol. 360, no. 1, pp. 109–113, 2006. View at: Publisher Site  Google Scholar
 R. Siegel and J. R. Howell, Thermal Radiation Heat Transfer, McGrawHill, New York, NY, USA, 1972.
 M. Agida and A. S. Kumar, “A Boubaker polynomials expansion scheme solution to random Love's equation in the case of a rational Kernel,” Electronic Journal of Theoretical Physics, vol. 7, no. 24, pp. 319–326, 2010. View at: Google Scholar
 P. Barry and A. Hennessy, “Meixnertype results for Riordan arrays and associated integer sequences,” Journal of Integer Sequences, vol. 13, no. 9, pp. 1–34, 2010. View at: Google Scholar
 A. Belhadj, O. F. Onyango, and N. Rozibaeva, “Boubaker polynomials expansion schemerelated heat transfer investigation inside keyhole model,” Journal of Thermophysics and Heat Transfer, vol. 23, no. 3, pp. 639–640, 2009. View at: Publisher Site  Google Scholar
 A. Belhadj, J. Bessrour, M. Bouhafs, and L. Barrallier, “Experimental and theoretical cooling velocity profile inside laser welded metals using keyhole approximation and Boubaker polynomials expansion,” Journal of Thermal Analysis and Calorimetry, vol. 97, no. 3, pp. 911–915, 2009. View at: Publisher Site  Google Scholar
 K. B. Ben Mahmoud and M. Amlouk, “The 3D AmloukBoubaker expansivityenergy gapVickers hardness abacus: a new tool for optimizing semiconductor thin film materials,” Materials Letters, vol. 63, no. 12, pp. 991–994, 2009. View at: Publisher Site  Google Scholar
 M. Benhaliliba, C. E. Benouis, K. Boubaker, M. Amlouk, and A. Amlouk, “A new guide to thermally optimized doped oxides monolayer spraygrown solar cells: the Amloukboubaker optothermal expansivity ψab in the book,” in Solar CellsNew Aspects and Solutions, L. A. Kosyachenko, Ed., pp. 27–41, 2011. View at: Google Scholar
 M. Dada, O. B. Awojoyogbe, and K. Boubaker, “Heat transfer spray model: an improved theoretical thermal timeresponse to uniform layers deposit using Bessel and Boubaker polynomials,” Current Applied Physics, vol. 9, no. 3, pp. 622–624, 2009. View at: Publisher Site  Google Scholar
 S. Fridjine and M. Amlouk, “A new parameter: an abacus for optimizing PVT hybrid solar device functional materials using the boubaker polynomials expansion scheme,” Modern Physics Letters B, vol. 23, no. 17, pp. 2179–2191, 2009. View at: Publisher Site  Google Scholar
 S. Fridjine, K. B. Ben Mahmoud, M. Amlouk, and M. Bouhafs, “A study of sulfur/selenium substitution effects on physical and mechanical properties of vacuumgrown ZnS1xSex compounds using Boubaker polynomials expansion scheme (BPES),” Journal of Alloys and Compounds, vol. 479, no. 12, pp. 457–461, 2009. View at: Publisher Site  Google Scholar
 J. Ghanouchi, H. Labiadh, and K. Boubaker, “An attempt to solve the heat transfer equation in a model of pyrolysis spray using 4qorder mboubaker polynomials,” International Journal of Heat and Technology, vol. 26, no. 1, pp. 49–53, 2008. View at: Google Scholar
 T. Ghrib, K. Boubaker, and M. Bouhafs, “Investigation of thermal diffusivitymicrohardness correlation extended to surfacenitrured steel using Boubaker polynomials expansion,” Modern Physics Letters B, vol. 22, no. 29, pp. 2893–2907, 2008. View at: Publisher Site  Google Scholar
 C. Khélia, K. Boubaker, T. Ben Nasrallah, M. Amlouk, and S. Belgacem, “Morphological and thermal properties of βSnS_{2} sprayed thin films using Boubaker polynomials expansion,” Journal of Alloys and Compounds, vol. 477, no. 12, pp. 461–467, 2009. View at: Publisher Site  Google Scholar
 A. S. Kumar, “An analytical solution to applied mathematicsrelated Love's equation using the Boubaker polynomials expansion scheme,” Journal of the Franklin Institute, vol. 347, no. 9, pp. 1755–1761, 2010. View at: Publisher Site  Google Scholar
 S. Lazzez, K. B. Ben Mahmoud, S. Abroug, F. Saadallah, and M. Amlouk, “A Boubaker polynomials expansion scheme (BPES)related protocol for measuring sprayed thin films thermal characteristics,” Current Applied Physics, vol. 9, no. 5, pp. 1129–1133, 2009. View at: Publisher Site  Google Scholar
 A. Milgram, “The stability of the Boubaker polynomials expansion scheme (BPES)based solution to LotkaVolterra problem,” Journal of Theoretical Biology, vol. 271, no. 1, pp. 157–158, 2011. View at: Publisher Site  Google Scholar
 H. Rahmanov, “A solution to the non linear KortewegDeVries equation in the particular case dispersionadsorption problem in porous media using the spectral Boubaker polynomials expansion scheme (BPES),” Studies in Nonlinear Sciences, vol. 2, no. 1, pp. 46–49, 2011. View at: Google Scholar
 S. Slama, M. Bouhafs, and K. B. Ben Mahmoud, “A boubaker polynomials solution to heat equation for monitoring A3 point evolution during resistance spot welding,” International Journal of Heat and Technology, vol. 26, no. 2, pp. 141–145, 2008. View at: Google Scholar
 A. Yildirim, S. T. MohyudDin, and D. H. Zhang, “Analytical solutions to the pulsed KleinGordon equation using Modified Variational Iteration Method (MVIM) and Boubaker Polynomials Expansion Scheme (BPES),” Computers and Mathematics with Applications, vol. 59, no. 8, pp. 2473–2477, 2010. View at: Publisher Site  Google Scholar
 E. J. Kontoghiorghes, “Parallel strategies for computing the orthogonal factorizations used in the estimation of econometric models,” Algorithmica, vol. 25, no. 1, pp. 58–74, 1999. View at: Google Scholar
 M. Benhamadou, “New method to solve large linear systems: the algorithm of ‘recursive reduction’,” Advances in Engineering Software, vol. 31, no. 7, pp. 481–491, 2000. View at: Publisher Site  Google Scholar
 A. H. Hassan, “Differential transformation technique for solving higherorder initial value problems,” Applied Mathematics and Computation, vol. 154, no. 2, pp. 299–311, 2004. View at: Publisher Site  Google Scholar
 M. Torabi, H. Yaghoobi, and A. Aziz, “Analytical solution for convectiveradiative continuously moving fin with temperaturedependent thermal conductivity,” International Journal of Thermophysics, vol. 33, pp. 924–941, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Mohsen Torabi 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.