Chinese Journal of EngineeringVolume 2013 (2013), Article ID 470696, 12 pageshttp://dx.doi.org/10.1155/2013/470696

Research Article

Analysis of Radiative Radial Fin with Temperature-Dependent Thermal Conductivity Using Nonlinear Differential Transformation Methods

Mohsen Torabi,¹ Hessameddin Yaghoobi,¹ Andrea Colantoni,² Paolo Biondi,² and Karem Boubaker³

¹Young Researchers and Elite Club, Central Tehran Branch, Islamic Azad University, Tehran 14174, Iran²Department of Agriculture, Forest, Nature and Energy (DAFNE), University of Tuscia, Via S. Camillo de Lellis snc, 01100 Viterbo, Italy³Unité de Physique des Dispositifs à Semi-Conducteurs, Faculté des Sciences de Tunis, Université de Tunis El Manar, 2092 Tunis, Tunisia

Received 8 August 2013; Accepted 1 September 2013

Academic Editors: B.-Y. Cao and J.-w. Zhou

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.

Abstract.

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.

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 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. 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 [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 temperature-dependent 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 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 [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/control-volume method. Joneidi et al. [25] provided analytical solution to fin efficiency of convective straight fins with temperature-dependent 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 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) [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.

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.

Figure 1: A heat pipe/fin radiating element.

The energy balance equation for a differential element of the fin [2–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–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 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

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

Table 1: The fundamental operations of differential transform method.

Original function
Transformed function

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 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 [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 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 [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, , 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).

(a)

(b)

Figure 2: Convergence test of the dimensionless temperature with DTM, for (dash-dot line), (dash line), (solid line), and NS (box symbol).

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.

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

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.

Table 2: The results of DTM, VIM, ADM, and NS for dimensionless temperature.

DTM
VIM [4]
ADM [31]
NS
DTM
VIM [4]
ADM [31]
NS

0
0.8133693583
0.8145587262
0.8122259487
0.8133695812
0.8679116912
0.8677934850
0.8678227415
0.8679120939

0.05
0.8137822645
0.8149299130
0.8126321774
0.8137824427
0.8682139355
0.8680864947
0.8681245987
0.8682140666

0.10
0.8150223528
0.8160457268
0.8138530912
0.8150225764
0.8691214181
0.8689665327
0.8690309943
0.8691217405

0.15
0.8170937494
0.8179129642
0.8158953383
0.8170939665
0.8706363930
0.8704366340
0.8705444008
0.8706367938

0.20
0.8200033838
0.8205430743
0.8187698837
0.8200036083
0.8727626384
0.8725018855
0.8726689432
0.8727630106

0.25
0.8237610702
0.8239523423
0.8224918404
0.8237612930
0.8755054875
0.8751694686
0.8754104067
0.8755058710

0.30
0.8283796228
0.8281621511
0.8270802450
0.8283798480
0.8788718740
0.8784487205
0.8787762466
0.8788722707

0.35
0.8338750084
0.8331993221
0.8325577851
0.8338752350
0.8828703906
0.8823512130
0.8827756075
0.8828707739

0.40
0.8402665394
0.8390965420
0.8389504951
0.8402667643
0.8875113635
0.8868908496
0.8874193466
0.8875117414

0.45
0.8475771116
0.8458928803
0.8462874349
0.8475773321
0.8928069419
0.8920839837
0.8927200706
0.8928073243

0.50
0.8558334915
0.8536344030
0.8546003708
0.8558337131
0.8987712062
0.8979495558
0.8986921830
0.8987715631

0.55
0.8650666613
0.8623748942
0.8639234800
0.8650668754
0.9054202946
0.9045092538
0.9053519505
0.9054206440

0.60
0.8753122271
0.8721766901
0.8742931088
0.8753124399
0.9127725525
0.9117876937
0.9127175865
0.9127729094

0.65
0.8866109039
0.8831116398
0.8857476040
0.8866111038
0.9208487049
0.9198126293
0.9208093621
0.9208490358

0.70
0.8990090866
0.8952622030
0.8983272589
0.8990092815
0.9296720572
0.9286151827
0.9296497443
0.9296723561

0.75
0.9125595243
0.9087226999
0.9120743995
0.9125597131
0.9392687258
0.9382301051
0.9392635667
0.9392690380

0.80
0.9273221155
0.9236007249
0.9270336526
0.9273222778
0.9496679059
0.9486960683
0.9496782441
0.9496681990

0.85
0.9433648484
0.9400187469
0.9432524344
0.9433650219
0.9609021785
0.9600559865
0.9609240279
0.9609023769

0.90
0.9607649133
0.9581159096
0.9607817046
0.9607650173
0.9730078657
0.9723573730
0.9730343160
0.9730082160

0.95
0.9796100255
0.9780500563
0.9796770290
0.9796103729
0.9860254404
0.9856527331
0.9860460201
0.9860261076

1
1.0000000000
1.0000000000
1.0000000000
1.0000000000
1.0000000000
1.0000000000
1.0000000000
1.0000000000

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.

Figure 4: Comparison for dimensionless temperature variation for and using DTM (solid line), BPES (triangle), NS (box symbol), and HAM (solid-box symbol).

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 5: Variation of dimensionless fin tip temperature for various with DTM (solid line) and NS (box symbol).

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.

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

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 7: Dimensionless temperature in fin for various with DTM (solid line) and NS (box symbol).

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 8: Dimensionless temperature in fin for various with DTM (solid line) and NS (box symbol).

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

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

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.

Figure 11: Variation of the fin efficiency with and with DTM (solid line) and NS (box symbol).

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

Nomenclature

: | Fin tip length, m
: | Domain
: | Constant
: | 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
: | Temperature,
: | Fin’s base temperature,
: | Radiation sink temperature,
: | Transformed analytical function
: | Original analytical function
: | Semithickness of the fin, m.

Greek Symbols

: | Thermal conductivity parameter
: | Emissivity
: | Fin efficiency
: | Slope of the thermal conductivity-temperature curve,
: | Stefan-Boltzmann constant,
: | Dimensionless temperature
: | 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.

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. 5-6, pp. 976–987, 2007.

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. 7-8, pp. 1545–1558, 2007.

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.

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. 17-18, pp. 2345–2352, 2008.

M. H. Sharqawy and S. M. Zubair, “Efficiency and optimization of straight fins with combined heat and mass transfer-an analytical solution,” Applied Thermal Engineering, vol. 28, no. 17-18, pp. 2279–2288, 2008.

G. Domairry and M. Fazeli, “Homotopy analysis method to determine the fin efficiency of convective straight fins with temperature-dependent thermal conductivity,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 2, pp. 489–499, 2009.

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.

D. B. Kulkarni and M. M. Joglekar, “Residue minimization technique to analyze the efficiency of convective straight fins having temperature-dependent thermal conductivity,” Applied Mathematics and Computation, vol. 215, no. 6, pp. 2184–2191, 2009.

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.

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. 11-12, pp. 2646–2659, 2009.

M. N. Bouaziz and A. Aziz, “Simple and accurate solution for convective-radiative fin with temperature dependent thermal conductivity using double optimal linearization,” Energy Conversion and Management, vol. 51, no. 12, pp. 2776–2782, 2010.

F. Khani and A. Aziz, “Thermal analysis of a longitudinal trapezoidal fin with temperature-dependent thermal conductivity and heat transfer coefficient,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 3, pp. 590–601, 2010.

J. K. Zhou, Differential Transform and Its Applications for Electrical CircuIts, Huarjung University Press, Wuhan, China, 1986.

I. H. Abdel-Halim Hassan, “Application to differential transformation method for solving systems of differential equations,” Applied Mathematical Modelling, vol. 32, no. 12, pp. 2552–2559, 2008.

A. S. V. Ravi Kanth and K. Aruna, “Solution of singular two-point boundary value problems using differential transformation method,” Physics Letters A, vol. 372, no. 26, pp. 4671–4673, 2008.

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.

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.

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.

H. Yaghoobi and M. Torabi, “Analytical solution for settling of non-spherical particles in incompressible Newtonian media,” Powder Technology, vol. 221, pp. 453–463, 2012.

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.

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.

H. P. Chu and C. Y. Lo, “Application of the hybrid differential transform-finite difference method to nonlinear transient heat conduction problems,” Numerical Heat Transfer A, vol. 53, no. 3, pp. 295–307, 2008.

C. Y. Lo and B. Y. Chen, “Application of hybrid differential transform/control-volume method to hyperbolic heat conduction problems,” Numerical Heat Transfer B, vol. 55, no. 3, pp. 219–231, 2009.

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.

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.

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 DTM-Padé,” International Journal of Thermal Sciences, vol. 49, no. 12, pp. 2405–2412, 2010.

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.

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.

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.

C. Arslanturk, “Optimum design of space radiators with temperature-dependent thermal conductivity,” Applied Thermal Engineering, vol. 26, no. 11-12, pp. 1149–1157, 2006.

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.

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.

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.

A. Belhadj, O. F. Onyango, and N. Rozibaeva, “Boubaker polynomials expansion scheme-related heat transfer investigation inside keyhole model,” Journal of Thermophysics and Heat Transfer, vol. 23, no. 3, pp. 639–640, 2009.

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.

K. B. Ben Mahmoud and M. Amlouk, “The 3D Amlouk-Boubaker expansivity-energy gap-Vickers hardness abacus: a new tool for optimizing semiconductor thin film materials,” Materials Letters, vol. 63, no. 12, pp. 991–994, 2009.

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.

M. Dada, O. B. Awojoyogbe, and K. Boubaker, “Heat transfer spray model: an improved theoretical thermal time-response to uniform layers deposit using Bessel and Boubaker polynomials,” Current Applied Physics, vol. 9, no. 3, pp. 622–624, 2009.

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.

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.

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.

T. Ghrib, K. Boubaker, and M. Bouhafs, “Investigation of thermal diffusivitymicrohardness correlation extended to surface-nitrured steel using Boubaker polynomials expansion,” Modern Physics Letters B, vol. 22, no. 29, pp. 2893–2907, 2008.

C. Khélia, K. Boubaker, T. Ben Nasrallah, M. Amlouk, and S. Belgacem, “Morphological and thermal properties of β-SnS₂ sprayed thin films using Boubaker polynomials expansion,” Journal of Alloys and Compounds, vol. 477, no. 1-2, pp. 461–467, 2009.

A. S. Kumar, “An analytical solution to applied mathematics-related Love's equation using the Boubaker polynomials expansion scheme,” Journal of the Franklin Institute, vol. 347, no. 9, pp. 1755–1761, 2010.

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.

A. Milgram, “The stability of the Boubaker polynomials expansion scheme (BPES)-based solution to Lotka-Volterra problem,” Journal of Theoretical Biology, vol. 271, no. 1, pp. 157–158, 2011.

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.

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.

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.

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.

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.

A. H. Hassan, “Differential transformation technique for solving higher-order initial value problems,” Applied Mathematics and Computation, vol. 154, no. 2, pp. 299–311, 2004.

M. Torabi, H. Yaghoobi, and A. Aziz, “Analytical solution for convective-radiative continuously moving fin with temperature-dependent thermal conductivity,” International Journal of Thermophysics, vol. 33, pp. 924–941, 2012.