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 iterationmethod (VIM), Adomian decompositionmethod (ADM), homotopy analysis method (HAM), and numerical solution (NS) in order to verify the accuracy of the proposedmethod.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.


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 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 [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 cross-section (SRC) profile subject to combined heat and mass transfer.Cos ¸kun and Atay [4,5] used variational iteration method to analyze convective straight and radial fins with temperature-dependent 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 temperature-dependent thermal conductivity.Khani et al. [10] 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 [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 convective-radiative fin with temperature-dependent 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 seminumerical-analytical 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][16][17][18][19][20][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/control-volume 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 two-dimensional 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 convective-radiative cooling of a lumped system with temperature-dependent specific heat.
More recently, Ganji et al. [29] applied DTM to the problem of convective-radiative 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 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.

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 energy balance equation for a differential element of the fin [2][3][4][5] is given by where () and  are thermal conductivity and the Stefan-Boltzmann constant, respectively.The thermal conductivity of the fin material is assumed to be a linear function of temperature [4][5][6][7] according to where  0 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][9][10][11] reduces to the following equation: with the following boundary conditions:

Fin Efficiency
The heat transfer rate from the surfaces of a fin is found by applying the Stefan-Boltzmann 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

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   = 0 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 T-function).
The differential spectrum of () is confined within the interval  ∈ [0, ], 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, (0) is known as the zero discrete, (1) 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  + 1 is total number of polynomial terms used in the DTM.Mathematical operations performed by differential transform method are listed in Table 1.

Solution Using the Differential Transformation Method
DTM. Now we apply differential transformation method to . Taking the differential transform of ( 4) with respect to , and considering  = 1 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 Θ( + 2) as follows: where  is constant, and we will calculate it with considering another boundary condition in (5b) in point  = 1.
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  = 1.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.
In reference to (4), the BPES is applied to the system through setting the expression where  4 are the 4-order Boubaker polynomials,   are  4 minimal positive roots,  0 is a prefixed integer, and   | =1⋅⋅⋅ 0 are unknown pondering real.The BPES solution is obtained through the following steps.
(i) Integrate, for a given value of  0 , 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   | =1⋅⋅⋅ 0 .Consequently, it comes for (4) and ( 25)-( 26) that with the matrix standard form  The system ( 28) is hence reduced to approximately (64 0 ) 3 arithmetical operations and solved using the Householder algorithm [53,54].(ii) Incremente  0 .
The final result is hence (obtained for  0 = 233)

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  0 , 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.,  = 30) 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  = −0.6 and  = 0.6 when  = 1, more than 10 terms are needed to obtain the value of the .Also, it is seen that if  = −0.6 and  = 100 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  = 30 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).

Validation of the Results.
To check the accuracy of the present solution, by considering   =   =  = 0 our problem is converted to a simpler case which was studied in four pioneering works [4,[30][31][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 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  = 0 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.

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  = 0 represents the case when the thermal conductivity is a constant and its value is  0 .The curves with  > 0 correspond to fin materials whose thermal conductivity increases as temperature increases.The converse is true of curves with  < 0. 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 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 [55] against numerical solution, and results of both methods were in good agreement with each other.

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 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.

Figure 3 :
Figure 3: Comparison for dimensionless temperature variation for various  and  between DTM (solid line), VIM (box symbol), and ADM (solid-circle symbol).

Figure 5 :
Figure 5: Variation of dimensionless fin tip temperature for various  with DTM (solid line) and NS (box symbol).

3 Figure 6 :
Figure 6: Dimensionless temperature in fin for various  with DTM (solid line) and NS (box symbol).

1 Figure 7 :
Figure 7: Dimensionless temperature in fin for various  with DTM (solid line) and NS (box symbol).

Figure 9 :
Figure 9: Dimensionless temperature in fin for various   with DTM (solid line) and NS (box symbol).

Figure 10 :Figure 11 :
Figure 10: Dimensionless temperature in fin for various  with DTM (solid line) and NS (box symbol).

Table 1 :
The fundamental operations of differential transform method.

Table 2 :
The results of DTM, VIM, ADM, and NS for dimensionless temperature.
Figure 8: Dimensionless temperature in fin for various   with DTM (solid line) and NS (box symbol).