- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Mathematical Problems in Engineering
Volume 2013 (2013), Article ID 273052, 14 pages
Analytical Solutions for Steady Heat Transfer in Longitudinal Fins with Temperature-Dependent Properties
Center for Differential Equations, Continuum Mechanics and Applications, School of Computational and Applied Mathematics, University of the Witwatersrand, Private Bag 3, Johannesburg 2050, South Africa
Received 30 January 2013; Accepted 24 April 2013
Academic Editor: Ireneusz Zbicinski
Copyright © 2013 Partner L. Ndlovu and Raseelo J. Moitsheki. 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.
Explicit analytical expressions for the temperature profile, fin efficiency, and heat flux in a longitudinal fin are derived. Here, thermal conductivity and heat transfer coefficient depend on the temperature. The differential transform method (DTM) is employed to construct the analytical (series) solutions. Thermal conductivity is considered to be given by the power law in one case and by the linear function of temperature in the other, whereas heat transfer coefficient is only given by the power law. The analytical solutions constructed by the DTM agree very well with the exact solutions even when both the thermal conductivity and the heat transfer coefficient are given by the power law. The analytical solutions are obtained for the problems which cannot be solved exactly. The effects of some physical parameters such as the thermogeometric fin parameter and thermal conductivity gradient on temperature distribution are illustrated and explained.
Fins are surfaces that extend from a hot object (body) to increase the rate of heat transfer to the surrounding fluid. In particular, fins are used extensively in various industrial applications such as the cooling of computer processors, air conditioning, and oil carrying pipe lines. A well-documented review of heat transfer in extended surfaces is presented by Kraus et al. . The problems on heat transfer particularly in fins continue to be of scientific interest. These problems are modeled by highly nonlinear differential equations which are difficult to solve exactly. However, Moitsheki et al. [2–4] have attempted to construct exact solutions for the steady state problems arising in heat flow through fins. A number of techniques, for example, Lie symmetry analysis , He’s variational iteration method , Adomain decomposition methods , homotopy perturbation methods , homotopy analysis methods , methods of successive approximations , and other approximation methods  have been used to determine solutions of the nonlinear differential equations describing heat transfer in fins.
Recently, the solutions of the nonlinear ordinary differential equations (ODEs) arising in extended surface heat transfer have been constructed using the DTM [11–19]. The DTM is an analytical method based on the Taylor series expansion and was first introduced by Zhou  in 1986. The DTM approximates the exact solution by a polynomial, and previous studies have shown that it is an efficient means of solving nonlinear problems or systems with varying parameters . Furthermore, DTM is a computational inexpensive tool for obtaining analytical solution, and it generalizes the Taylor method to problems involving procedures such as fractional derivative (see e.g., [22–24]). Also, this method converges rapidly (see e.g., ).
Models arising in heat transfer through fins may contain temperature-dependent properties such as thermal conductivity and heat transfer coefficient. The dependency of thermal conductivity and heat transfer coefficient on temperature renders such problems highly nonlinear and difficult to solve, particularly exactly. Thermal conductivity may be modeled for many engineering applications by the power law and by linear dependency on temperature. On the other hand, heat transfer coefficient can be expressed as a power law for which values of the exponent represent different phenomena (see e.g., ).
In this paper, the DTM is employed to determine the analytical solutions to the nonlinear boundary value problem describing heat transfer in longitudinal fins of rectangular, exponential, and convex parabolic profiles. Both thermal conductivity and heat transfer coefficient are temperature dependent. We adopt the terminology exact solutions to refer to solutions given in terms of fundamental expressions such as logarithmic, trigonometric, and exponential. However, analytical solutions will be series solutions and in particular, those constructed using the DTM. The mathematical modelling of the problem under consideration is described in Section 2. A brief discussion on the fundamentals of the DTM is be provided in Section 3. The comparison of the exact and analytical solutions constructed by DTM is given in Section 4. In Section 5, we provide analytical solutions for the heat transfer in longitudinal fins of various profiles. Here, heat transfer coefficient is given by the power law, and we consider two cases of the thermal conductivity, namely, the power law and the linear function of temperature. Furthermore, we describe the fin efficiency and the heat flux in Section 6. Some exciting results are discussed in Section 7. Lastly, the concluding remarks are provided in Section 8.
2. Mathematical Models
We consider a longitudinal one dimensional fin of cross-sectional area . The perimeter of the fin is denoted by and its length by . The fin is attached to a fixed prime surface of temperature and extends to an ambient fluid of temperature . The fin thickness at the prime surface is given by and its profile is given by . Based on the one dimensional heat conduction, the energy balance equation is then given by (see e.g., ) where and are nonuniform temperature-dependent thermal conductivity and heat transfer coefficients, respectively, is the temperature distribution, is the fin profile, and is the space variable. The length of the fin is measured from the tip to the prime surface as shown in Figure 1. Assuming that the fin tip is adiabatic (insulated) and the base temperature is kept constant, then the boundary conditions are given by Introducing the following dimensionless variables (see e.g., ): with being defined as the thermal conductivity at the ambient temperature and as the heat transfer at the prime surface (fin base), reduces (1) to Here is the dimensionless temperature, is the dimensionless space variable, is the dimensionless fin profile, is the dimensionless thermal conductivity, is the dimensionless heat transfer coefficient, and is the thermogeometric fin parameter. The dimensionless boundary conditions then become We assume that the heat transfer coefficient is given by the power law where the exponent is a real constant. In fact the values of may vary as between and . However, in most practical applications, it lies between and . In dimensionless variables, we have . We consider the two distinct cases of the thermal conductivity as follows: (a)the power law with being a constant and (b)the linear function The dimensionless thermal conductivity given by the power law and the linear function of temperature is and , respectively. Here the thermal conductivity gradient is . Furthermore, we consider a various fin profiles including the longitudinal rectangular , the longitudinal convex parabolic , and the exponential profile with being the constant (see also ).
3. Fundamentals of the Differential Transform Method
In this section, the basic idea underlying the DTM is briefly introduced. Let be an analytic function in a domain . The Taylor series expansion function of with the center located at is given by 
The particular case of (10) when is referred to as the Maclaurin series expansion of and is expressed as The differential transform of is defined as follows: where is the original analytic function and is the transformed function. The differential spectrum of is confined within the interval , where is a constant. From (11) and (12), the differential inverse transform of is defined as follows: and if is expressed by a finite series, then
Some of the useful mathematical operations performed by the differential transform method are listed in Table 1.
The delta function is given by
4. Comparison of Exact and Analytical Solutions
In this section, we consider a model describing temperature distribution in a longitudinal rectangular fin with both thermal conductivity and heat transfer coefficient being functions of temperature given by the power law (see e.g., ). The exact solution of (4) when both the power laws are given by the same exponent is given by  We use this exact solution as a benchmark or validation of the DTM. The effectiveness of the DTM is determined by comparing the exact and the analytical solutions. We compare the results for the cases and with fixed values of .
Applying the DTM to (4) with the power law thermal conductivity, (rectangular profile) and given , one obtains the following recurrence relation: Exerting the transformation to the boundary condition (6) at a point , The other boundary conditions are considered as follows: where is a constant. Equation (17) is an iterative formula of constructing the power series solution as follows:
To obtain the value of , we substitute the boundary condition (5) into (34) at the point . Thus, we have We omit presenting the tedious process of finding . However, one may obtain the expression for upon substituting the obtained value of into (34).
Following a similar approach given in Section 4.1 and given , one obtains the analytical solution The constant may be obtained using the boundary condition at the fin base. The comparison of the DTM and the exact solutions are reflected in Tables 2 and 3 for different values of . Furthermore, the comparison of the exact and the analytical solutions is depicted in Figures 2(a) and 2(b).
5. Analytical Solutions
It is well known that exact solutions for ODEs such as (4) exist only when thermal conductivity and the term containing heat transfer coefficient are connected by differentiation (or simply if the ODE such as (4) is linearizable) . In this section we determine the analytical solutions for the nonlinearizable (4), firstly when thermal conductivity is given by the power law and secondly as a linear function of temperature. In both cases and throughout this paper, the heat transfer coefficient is assumed to be a power law function of temperature. These assumptions of the thermal properties are physical realistic. We have noticed that DTM runs into difficulty when the exponent of the power law of the thermal conductivity is given by fractional values and also when the function is given in terms of fractional power law. One may follow Moradi and Ahmadikia  by introducing a new variable to deal with fractional powers of , and on the other hand, it is possible to remove the fractional exponent of the heat transfer coefficient by fundamental laws of exponent and binomial expansion.
Proposition 1. A nonlinear ODE such as (4) may admit the DTM solution if is a constant or exponential function. However, if is a power law, then such an equation admits a DTM solution if the product, holds. Here is a real constant.
(a) If , then transforms (4) into
only if .
This example implies that the DTM may only be applicable to problems describing heat transfer in fins with convex parabolic profile. In the next subsections, we present analytical solutions for (4) with various functions and .
5.1. The Exponential Profile and Power Law Thermal Conductivity
In this section, we present solutions for equation describing heat transfer in a fin with exponential profile and power law thermal conductivity and heat transfer coefficient. That is, given (4) with , and both heat transfer coefficient and thermal conductivity being power law functions of temperature, we construct analytical solutions. In our analysis, we consider and indicating the fin subject to nucleate boiling and radiation into free space at zero absolute temperature, respectively. Firstly, given and and applying the DTM, one obtains the following recurrence revelation:
The pervious process is continuous and one may consider as many terms as desired (but bearing in mind that DTM converges quite fast). Substituting (18) to (19) and (42) into (13), we obtain the following closed form of the solution:
Substituting this value of into (43), one finds the expression for . On the other hand, given , one obtains the solution Here, the constant maybe obtained by evaluating the boundary condition . The solutions (43) and (45) are depicted in Figures 3(b) and 3(a), respectively.
5.2. The Rectangular Profile and Power Law Thermal Conductivity
In this section, we provide a detailed construction of analytical solutions for the heat transfer in a longitudinal rectangular fin with a power law thermal conductivity; that is, we consider (4) with and . The analytical solutions are given in the following expressions.
(a) Case (b) Case The constant is obtained by solving the appropriate at the fin base boundary condition. The analytical solutions in (46) and (47) are depicted in Figures 4(a) and 4(b), respectively.
5.3. The Convex Parabolic Profile and Power Law Thermal Conductivity
In this section, we present solutions for the equation describing the heat transfer in a fin with convex parabolic profile and power law thermal conductivity. Equation (40) is considered. Here we consider the values . The final analytical solution is given by
5.4. The Rectangular Profile and Linear Thermal Conductivity
In this section, we present solutions for the equation representing the heat transfer in a fin with rectangular profile and the thermal conductivity depending linearly on temperature. That is, we consider (4) with and . The analytical solutions for this problem for different values of are given by
5.5. The Convex Parabolic Profile and Linear Thermal Conductivity
In this section, we present solutions for the equation describing the heat transfer in a fin with convex parabolic profile and the thermal conductivity depending linearly on temperature. That is, we consider (40) with . The analytical solution for this problem for different values of is given by
5.6. The Exponential Profile and Linear Thermal Conductivity
In this section, we present solutions for equation heat transfer in a fin with convex parabolic profile and the thermal conductivity depending linearly on temperature. That is, we consider (4) with and . Here is a constant. The analytical solutions for this problem for different values of are given by
6. Fin Efficiency and Heat Flux
6.1. Fin Efficiency
The heat transfer rate from a fin is given by Newton’s second law of cooling: Fin efficiency is defined as the ratio of the fin heat transfer rate to the rate that would be if the entire fin were at the base temperature and is given by (see e.g., ) In dimensionless variables, we have We consider the solutions (50), (51), (52), and (53) and depict the fin efficiency (63) in Figure 9.
6.2. Heat Flux
The fin base heat flux is given by the Fourier’s law The total heat flux of the fin is given by  Introducing the dimensionless variable as described in Section 2 implies where the dimensionless parameter is the Biot number. We consider a number of cases for the thermal conductivity and the heat transfer coefficient.
6.2.1. Linear Thermal Conductivity and Power Law Heat Transfer Coefficient
6.2.2. Power Law Thermal Conductivity and Heat Transfer Coefficient
7. Some Discussions
The DTM has resulted in some interesting observations and study. We have observed in Figures 2(a) and 2(b) an excellent agreement between the analytical solutions generated by DTM and the exact solution obtained in . In particular, we considered a fin problem in which both thermal conductivity and heat transfer coefficient are given by the same power law. Furthermore, we notice from Table 2 that an absolute error of approximately is produced by DTM of order . In Table 3, an absolute error of approximately is produced for the same order. This confirms that the DTM converges faster and can provide accurate results with a minimum computation. As such, a tremendous confidence in the DTM in terms of the accuracy and effectiveness was built, and thus we used this method to solve other problems for which exact solutions are harder to construct.
In Figures 3(a), 3(b), 4(a), 4(b), 5(a), and 5(b), we observe that the fin temperature increases with the decreasing values of the thermogeometric fin parameter. Here, the values of the exponents are fixed. Also, we observe that fin temperature is higher when , that is, when heat transfer coefficient is higher than the thermal conductivity. We observe in Figures 6, 7, and 8 that the fin temperature increases with the increasing values of . Furthermore, it appears that the fin with exponential profile performs the least in transferring the heat from the base, since the temperature in such a fin is much higher than that of the rectangular and the convex parabolic profiles. In other words, heat dissipation to the fluid surrounding the extended surface is much faster in longitudinal fins of rectangular and convex parabolic profiles. In Figure 9, fin efficiency decreases with increasing thermogeometric fin parameter. Also, fin efficiency increases with increasing values of . It is easy to show that the thermogeometric fin parameter is directly proportional to the aspect ratio (extension factor) with square root of the Biot number being the proportionality constant. As such, shorter fins are more efficient than longer ones. Else, the increased Biot number results in less efficient fin whenever the space is confined, that is, where the length of the fin cannot be increased. Figure 10 depicts the heat flux at the fin base. The amount of heat energy dissipated from the fin base is of immense interest in engineering . We observe in Figure 10 that the base heat flux increases with the thermogeometric fin parameter for considered values of the exponent (see also ). Figures 11(a), 11(b), and 11(c) display the heat flux across the fin length. We note that the heat flux across the fin length increases with increasing values of the thermogeometric fin parameter.
8. Concluding Remarks
In this study, we have successfully applied the DTM to highly nonlinear problems arising in heat transfer through longitudinal fins of various profiles. Both thermal conductivity and heat transfer coefficient are given as functions of temperature. The DTM agreed well with exact solutions when the thermal conductivity and heat transfer coefficient are given by the same power law. A rapid convergence to the exact solution was observed. Following the confidence in DTM built by the results mentioned, we then solved various exciting problems. The exotic results have been shown in tables and figures listed in this paper.
The results obtained in this paper are significant improvements on the known results. In particular, both the heat transfer coefficient and the thermal conductivity are allowed to be given by the power law functions of temperature, and also we considered a number of fin profiles. We note that exact solutions are difficult if not impossible to construct when the exponents of these properties are distinct.
Perhaps the notable advantage of the DTM is the generalization of the Taylor method to problems involving unusual derivative procedures such as fractional, fuzzy, or q-derivative . Some generalizations have been made by Odibat et al. , and they referred to their new method as the Generalized Differential Transform Method (GDTM). This showed great improvement compared to the Fractional Differential Transform Method (FDTM) introduced by Arikoglu and Ozkol .
We have shown with the help of an example that DTM may only be applied to fin problems involving heat transfer through fins with convex parabolic profile. Note that given an ODE such as (4) with a power law heat transfer coefficient of a fractional exponent, then one can easily remove the fraction by basic exponential rules and employment of the Binomial expansion. However, using the DTM, one runs into difficulty if the power law thermal conductivity in the same equation is given by the fractional exponent. We do not know whether these observations call for the “modified” DTM to solve problems arising in heat flow through fins with other profiles, such as longitudinal triangular and concave parabolic, and also with fractional power law thermal conductivity.
The main results obtained in this paper give insight into heat transfer in boiling liquids where the heat transfer coefficient is temperature dependent and may be given by a power law. The thermal conductivity of some materials such as gallium nitride (GaN) and Aluminium Nitride (AlN) may be modeled by power law temperature dependency [29, 30]. Thus, the solutions constructed here give a better comparison of heat transfer in terms of material used since in many engineering applications thermal conductivity is given as a linear function of temperature. Furthermore, a good study in terms of performance and efficiency of fin with different profiles is undertaken. These finding could help in the design of fins. It is claimed in  that DTM results are more accurate than those constructed by Variational Iteration Methods (VIM) and Homotopy Perturbation Methods (HPM). However, it would be risky to use the DTM approximate solutions as benchmarks for the numerical schemes. Nevertheless, we have also shown that DTM converges rapidly in just fifteen terms to the exact solution.
The authors wish to thank the National Research Foundation of South Africa and the University of the Witwatersrand for the generous financial support. The authors are grateful to Professor D. P. Mason and the anonymous reviewers for their invaluable comments which improved this paper significantly.
- A. D. Kraus, A. Aziz, and J. Welty, Extended Surface Heat Transfer, Wiley, New York, NY, USA, 2001.
- R. J. Moitsheki, T. Hayat, and M. Y. Malik, “Some exact solutions of the fin problem with a power law temperature-dependent thermal conductivity,” Nonlinear Analysis. Real World Applications, vol. 11, no. 5, pp. 3287–3294, 2010.
- R. J. Moitsheki, “Steady one-dimensional heat flow in a longitudinal triangular and parabolic fin,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 10, pp. 3971–3980, 2011.
- R. J. Moitsheki, “Steady heat transfer through a radial fin with rectangular and hyperbolic profiles,” Nonlinear Analysis. Real World Applications, vol. 12, no. 2, pp. 867–874, 2011.
- Mo. Miansari, D. D. Ganji, and Me. Miansari, “Application of He's variational iteration method to nonlinear heat transfer equations,” Physics Letters A, vol. 372, no. 6, pp. 779–785, 2008.
- C. H. Chiu and C. K. Chen, “A decomposition method for solving the convective longitudinal fins with variable thermal conductivity,” International Journal of Heat and Mass Transfer, vol. 45, no. 10, pp. 2067–2075, 2002.
- A. Rajabi, “Homotopy perturbation method for fin efficiency of convective straight fins with temperature-dependent thermal conductivity,” Physics Letters A, vol. 364, no. 1, pp. 33–37, 2007.
- 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.
- A. Campo and R. J. Spaulding, “Coupling of the methods of successive approximations and undetermined coefficients for the prediction of the thermal behaviour of uniform circumferential fins,” Heat and Mass Transfer, vol. 34, no. 6, pp. 461–468, 1999.
- A. Campo and F. Rodrfguez, “Approximate analytic temperature solution for uniform annular fins by adapting the power series method,” International Communications in Heat and Mass Transfer, vol. 25, no. 6, pp. 809–818, 1998.
- R. Chiba, “Application of differential transform method to thermoelastic problem for annular disks of variable thickness with temperature-dependent parameters,” International Journal of Thermophysics, vol. 33, pp. 363–380, 2012.
- S. Sadri, M. R. Raveshi, and S. Amiri, “Efficiency analysis of straight fin with variable heat transfer coefficient and thermal conductivity,” Journal of Mechanical Science and Technology, vol. 26, no. 4, pp. 1283–1290, 2012.
- M. Torabi, H. Yaghoori, and A. Aziz, “Analytical solution for convective-radiative continously moving fin with temperature dependent thermal conductivity,” International Journal of Thermophysics, vol. 33, pp. 924–941, 2012.
- M. Torabi and H. Yaghoobi, “Two dominant analytical methods for thermal analysis of convective step fin with variable thermal conductivity,” Thermal Science. In press.
- A. Moradi, “Analytical solutions for fin with temperature dependant heat transfer coefficient,” International Journal of Engineering and Applied Sciences, vol. 3, no. 2, pp. 1–12, 2011.
- A. Moradi and H. Ahmadikia, “Analytical solution for different profiles of fin with temperature-dependent thermal conductivity,” Mathematical Problems in Engineering, vol. 2010, Article ID 568263, 15 pages, 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.
- S. Ghafoori, M. Motevalli, M. G. Nejad, F. Shakeri, D. D. Ganji, and M. Jalaal, “Efficiency of differential transformation method for nonlinear oscillation: Comparison with HPM and VIM,” Current Applied Physics, vol. 11, no. 4, pp. 965–971, 2011.
- 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.
- J. K. Zhou, Differential Transform Method and Its Applications for Electric circuIts, Huazhong University Press, Wuhan, China, 1986.
- S. Mukherjee, “Reply to comment on ‘solutions of the duffing-van der pol oscillator equation by the differential transform method’,” Physica Scripta, vol. 84, Article ID 037003, 2 pages, 2011.
- C. Bervillier, “Status of the differential transformation method,” Applied Mathematics and Computation, vol. 218, no. 20, pp. 10158–10170, 2012.
- A. Arikoglu and I. Ozkol, “Solution of fractional differential equations by using differential transform method,” Chaos, Solitons and Fractals, vol. 34, no. 5, pp. 1473–1481, 2007.
- Z. Odibat, S. Momani, and V. S. Erturk, “Generalized differential transform method: application to differential equations of fractional order,” Applied Mathematics and Computation, vol. 197, no. 2, pp. 467–477, 2008.
- C. W. Bert, “Application of differential transform method to heat conduction in tapered fins,” Journal of Heat Transfer, vol. 124, no. 1, pp. 208–209, 2002.
- 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.
- H. C. Ünal, “An anlytical study of boiling heat transfer from a fin,” International Journal of Heat and Mass Transfer, vol. 31, no. 7, pp. 1483–1496, 1988.
- M. H. Chang, “A decomposition solution for fins with temperature dependent surface heat flux,” International Journal of Heat and Mass Transfer, vol. 48, no. 9, pp. 1819–1824, 2005.
- A. Jezowski, B. A. Danilchenko, M. Boćkowski et al., “Thermal conductivity of GaN crystals in 4.2–300 K range,” Solid State Communications, vol. 128, no. 2-3, pp. 69–73, 2003.
- S. Vitanov, V. Palankovski, S. Maroldt, and R. Quay, “High-temperature modeling of AlGaN/GaN HEMTs,” Solid-State Electronics, vol. 54, no. 10, pp. 1105–1112, 2010.