Analysis of Radiative Radial Fin with Temperature-Dependent Thermal Conductivity Using Nonlinear Differential Transformation Methods
Radiative radial fin with temperature-dependent thermal conductivity is analyzed. The calculations are carried out by using differential transformation method (DTM), which is a seminumerical-analytical 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 radiation-conduction fin parameter, radiation sink temperature, heat generation, and thermal conductivity parameters.
Extended surfaces are extensively used in various industrial applications. An extensive review on this topic is presented by Kraus et al. . 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 well-established 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  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  analytically analyzed the performance and optimization of longitudinal and pin fins of step reduction in local cross-section (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 temperature-dependent thermal conductivity. Sharqawy and Zubair  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  solved nonlinear straight fin differential equations by homotopy analysis method (HAM) to evaluate the temperature distribution within the fin. Arslanturk  obtained correlation equations for optimum design of annular fins with temperature-dependent thermal conductivity. Kulkarni and Joglekar  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 temperature-dependent thermal conductivity. Khani et al.  used HAM to evaluate the analytical approximate solutions and efficiency of the nonlinear fin problem with temperature-dependent thermal conductivity and heat transfer coefficient. Kundu  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  introduced a new concept called the double optimal linearization method (DOLM) to derive simple and accurate expressions for predicting the thermal performance of a convective-radiative fin with temperature-dependent thermal conductivity. Khani and Aziz  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 temperature-dependent.
The differential transformation method (DTM) is a seminumerical-analytical method. DTM, which is based on the Taylor series expansion, was first proposed by Zhou  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  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  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  proposed an alternative numerical method to investigate hyperbolic heat conduction problems using the hybrid differential transfer/control-volume method. Joneidi et al.  provided analytical solution to fin efficiency of convective straight fins with temperature-dependent thermal conductivity by DTM. Jang et al.  investigated and characterized a two-dimensional thermal conductive boundary value problem with discontinuous boundary and initial conditions. Rashidi et al.  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  applied DTM to solve the problems of convective and convective-radiative cooling of a lumped system with temperature-dependent specific heat.
More recently, Ganji et al.  applied DTM to the problem of convective-radiative straight fins with temperature-dependent 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) , homotopy analysis method (HAM) , 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 radiation-conduction 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 temperature-dependent thermal conductivity, , which depends on temperature linearly, and fin is diffuse-grey 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 one-dimensional, because the fin is assumed to be thin. Hence, only fin tip length is considered as the computational domain.
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:
3. Fin Efficiency
The heat transfer rate from the surfaces of a fin is found by applying the Stefan-Boltzmann law, namely , 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 :
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  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 T-function). 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
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 Newton-Raphson iterative technique for determination of unknown tip temperature . We employed the Maple’s built-in Roots command which numerically approximates the roots of an algebraic function using the specified method, such as Newton-Raphson, bisection, secant, and fixed-point iteration methods, and returns the specified outputs. This command uses the Newton-Raphson 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 : And the ADM solution of (4), with , is obtained in the following form :
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 applied-physics 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 , heat transfer , homodynamic [44, 45], material characterization , 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, , radiation-conduction 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 high-order 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 Runge-Kutta-Fehlberg method which is a well-tested numerical solution. Moreover, since the problem is a highly nonlinear radiative equation, an exact analytical solution cannot be found for (4).
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  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 radiation-conduction fin parameter, results of the present analysis are tabulated against the ADM, VIM, and again with numerical solution (NS), of the fourth-fifth order Runge-Kutta-Fehlberg 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  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 radiation-conduction 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 radiation-conduction 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 radiation-conduction 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 convective-radiative fins  against numerical solution, and results of both methods were in good agreement with each other.
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 temperature-dependent 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 radiation-conduction 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.
|:||Fin tip length, m|
|:||Temperature-dependent thermal conductivity,|
|:||Thermal conductivity at the base temperature,|
|:||Volumetric heat generation,|
|:||Dimensionless volumetric heat generation|
|:||Heat transfer rate from the surfaces of a fin|
|:||Fin’s base temperature,|
|:||Radiation sink temperature,|
|:||Transformed analytical function|
|:||Original analytical function|
|:||Semithickness of the fin, m.|
|:||Thermal conductivity parameter|
|:||Slope of the thermal conductivity-temperature curve,|
|:||Dimensionless radiation sink temperature|
|:||Radiation-conduction fin parameter.|
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
A. D. Kraus, A. Aziz, and J. R. Welty, Extended Surface Heat Transfer, John Wiley, New York, NY, USA, 2002.
S. B. Coşkun and M. T. Atay, “Analysis of convective straight and radial fins with temperature-dependent 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
F. Khani, M. A. Raji, and H. H. Nejad, “Analytical solutions and efficiency of the nonlinear fin problem with temperature-dependent 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
J. K. Zhou, Differential Transform and Its Applications for Electrical CircuIts, Huarjung University Press, Wuhan, China, 1986.
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
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
R. Siegel and J. R. Howell, Thermal Radiation Heat Transfer, McGraw-Hill, 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, “Meixner-type 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, 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
M. Benhaliliba, C. E. Benouis, K. Boubaker, M. Amlouk, and A. Amlouk, “A new guide to thermally optimized doped oxides monolayer spray-grown solar cells: the Amlouk-boubaker optothermal expansivity ψab in the book,” in Solar Cells-New Aspects and Solutions, L. A. Kosyachenko, Ed., pp. 27–41, 2011.View at: 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 vacuum-grown ZnS1-xSex compounds using Boubaker polynomials expansion scheme (BPES),” Journal of Alloys and Compounds, vol. 479, no. 1-2, 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 4q-order m-boubaker polynomials,” International Journal of Heat and Technology, vol. 26, no. 1, pp. 49–53, 2008.View at: Google Scholar
H. Rahmanov, “A solution to the non linear Korteweg-De-Vries equation in the particular case dispersion-adsorption 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. Mohyud-Din, and D. H. Zhang, “Analytical solutions to the pulsed Klein-Gordon 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