Mathematical Methods and Models in the Natural to the Life Sciences
View this Special IssueResearch Article  Open Access
Comparison of Exact Solutions for Heat Transfer in Extended Surfaces of Different Geometries
Abstract
We consider a steady state problem for heat transfer in fins of various geometries, namely, rectangular, radial, and spherical. The nonlinear steady state problem is linearizable provided that the thermal conductivity is the differential consequence of the term involving the heat transfer coefficient. As such, one is able to construct exact solutions. On the other hand, we employ the Lie point symmetry methods when the problem is not linearizable. Some interesting results are obtained and analyzed. The effects of the parameters such as thermogeometric fin parameter and the exponent on temperature are studied. Furthermore, fin efficiency and heat flux along the fin length of a spherical geometry are also studied.
1. Introduction
Heat transfer rate from a hot body to the surrounding may be increased by surfaces which extended into that surrounding. These extended surfaces are referred to as fins. Extended surfaces are found in many engineering appliances. Thus, mathematical modeling of the heat transfer through these surfaces and the solution are of continued interest. The heat transfer in fins is governed by boundary value problems (BVPs) which are rendered highly nonlinear by the dependency of thermal properties on temperature. In this study, both the heat transfer coefficient and thermal conductivity are given as a power law function of temperature.
The interest in solutions of fin problems continues unabated. Many symmetry analysts [1–5] analysed the fin equation when the heat transfer coefficient is given as a function of a space variable. Such a function was classified by direct methods [2] and the extended analysis was done in [4]. Only general solutions were provided in this case. It was claimed that exact solutions of steady fin problems exist only when thermal conductivity and heat transfer coefficients are constants [6]. However Moitsheki et al. [7] have shown that solutions may exist even when these thermal parameters are temperature dependent. In recent years Moitsheki [8, 9] and Moitsheki and Mhlongo [10] constructed exact solutions for the convective heat transfer in fins of different profiles. Furthermore, Ndlovu and Moitsheki [11] provided the approximate analytical solutions to steady state heat transfer in fins of different profiles which could not be solved exactly. In their studies an excellent comparison between exact and approximate solutions was established. One may also refer to the work by Moradi [12] and many other scholars.
In this study, we consider heat conduction problem in fins of different geometries and in particular the spherical fin which has never been studied before. We compare the exact solutions of heat transfer in rectangular, radial, and spherical fins. We further compare the fin efficiencies and effectiveness of these fins and determine the effects of thermal parameters in a spherical fin. This paper is arranged as follows. In Section 2, we present the description of the models considered. In Section 3 we briefly discuss the Lie point symmetry methods. Following linearization, the exact solution is provided in Section 4. In Section 5, we analyze the problem when linearization fails. Conclusions are provided in Section 6.
2. Mathematical Description of a Fin Problem
We consider a fin of an arbitrary geometry with the length (or radius) and a crosssectional area . The perimeter of the fin is given by . The fin is attached to a fixed prime surface of temperature and extends to an ambient fluid of temperature . The energy balance equation is given by and the relevant boundary conditions are The Greek letter represents different geometries; for example , and represent the longitudinal rectangular, the rectangular radial, and the spherical fin as shown in Figures 1, 2, and 3.
Introducing the nondimensional variables and numbers, and then (1) and the boundary conditions (2) become Two main cases may be analyzed, namely, and . One may construct exact solution when since the problem is linearizable and employ symmetry methods given .
3. A brief Account on Lie Point Symmetry Analysis
In this section we provide a brief theory of Lie point symmetry techniques. In short, a symmetry of a differential equation is an invertible transformation of the dependent and independent variables that does not change the original differential equation. Detailed theory and applications of Lie symmetry groups may be found in the texts such as those of [13–17]. Since in this study we deal with nonlinear second order ODEs therefore we will restrict our discussion to the determination of symmetries for such equations.
We seek transformations of the form called the infinitesimal transformations generated by the vector which leave the given differential equation invariant. Here is a group parameter. The group generated by transformations (6) is called a oneparameter group of transformations. If the given equation is of second order, for instance, then we extend the symmetry generator (7) to where with being the total derivative operator defined as
The invariance surface condition yields the overdetermined system of linear equations called the determining equation which may be solved to obtain the admitted symmetry generators (or equivalently symmetry transformation groups). In our analysis we determine the symmetries admitted by the single governing equation rather than the boundary value problem (BVP). Usually the symmetry algebra for the BVP is less in dimensions than that of the governing equation (see also [13]).
4. Linearization and Exact Solutions
It has been proven in [9] that equation such as (4) is linearizable provided that . Thus assuming and letting then (4)
Several subcases arise, namely, and arbitrary, given and .
Case 1 (, and ). This case has been solved in [7]. In this case the solutions are given by The solution for is given in terms of sine and cosine functions.
Case 2 (, and ). This case has been solved in [8]. In this case the exact solutions are given: The solutions for are given in terms of Bessel functions.
Case 3 (, and ). In this case we obtain the exact solutions
The solutions (14), (16), and (17) are depicted in Figure 4. Figures 5 and 6 depict the plot of solution (17) for varying and , respectively.
Case 4 (, and ). Given an arbitrary , we obtain the general solutions
Note that one may obtain exact solutions which satisfy the boundary conditions only when but this will coincide with solution (16). One may also construct exact solution when . In this case the solution satisfying the boundary condition is given by
The solution (20) is depicted in Figure 7.
5. Symmetry Reductions
In this section we consider the case , with and being arbitrary. In this case (4) is not linearizable and as such we employ the Lie point symmetry analysis.
5.1. Case
In this case (4) admits onedimensional Lie algebra spanned by the base vector Note that an obvious extra symmetry when is a translation in . We employ method of differential invariants to reduce the order of (4) by one. The first prolongation of the generator is given by and the corresponding characteristic equation is given by Solving the above characteristics gives the invariants Now we let , and write . From the definition of and by using chain rule we obtain Substituting into (8) we have We notice that the above equation may not be solved exactly.
5.1.1. Subcase
This choice of is not physical since there is no relationship between the geometry of the fin and the exponents of the thermal conductivity and heat transfer coefficient. However, this case is mathematically interesting since (4) admits two Lie point symmetries which implies [18] that the ODE in question is integrable or reducible to the one with cubic degree in first derivative. The admitted Lie algebra is spanned by the base vectors We omit further analysis since the initial assumption is not physically realistic.
5.1.2. Subcase
In this subcase (4) is integrable and admits an eightdimensional Lie algebra spanned by the base vectors Equation (4) is equivalent to the simple motion equation [18]. We adopt method of canonical coordinate to demonstrate this claim. We introduce the method of finding the solutions using and from the above dimensional Lie algebra vectors. The two symmetries lead to the canonical variables and the corresponding canonical forms of and are Writing transforms (4) into Integrating the latter equation and writing it into its original variables we obtain Imposing the boundary conditions, we obtain
5.1.3. Subcase
In this subcase the admitted eightdimensional Lie symmetry algebra includes
Since (4) admits eight symmetries it is linearizable (see also [18]). Using the symmetry generators and , we obtain the canonical variables The corresponding generators in canonical variables are Writing transforms (4) into
6. Fin Efficiency and Heat Flux
The fin efficiency is defined as the ratio of the actual heat transfer from the fin surface to the surrounding fluid while the whole fin is kept at the same temperature. On the other hand, heat flux is the total amount of heat flowing per unit area per unit time. The fin efficiency and the heat flux in dimensionless variables are given by respectively. Here the dimensionless parameter is the Biot number.
Given solution (17), we obtain where is the exponential integral [19]. Fin efficiency (41) is depicted in Figure 8.
And the heat flux becomes Heat flux (42) is depicted in Figures 9 and 10.
7. Some Discussions and Concluding Remarks
In Figure 4, we observe that heat is transferred much slower in spherical fins than in radial and rectangular fins. This is also confirmed by the values in Table 1. Now we focus only on spherical fins and observe in Figure 5 that temperature decreases with increase in the values of . Recall that the thermogeometric fin parameter is directly proportional to the aspect ratio of the fin. Thus longer fins ( larger) release heat much more efficiently that shorter ones. In Figure 6, temperature increases with an increase in the values of . Finally Figure 7 depicts the heat transfer where . Figure 8 is a plot of the fin efficiency with varying values of . Figures 9 and 10 depict the heat flux. In this paper we focused on the comparison of temperature distribution (or heat transfer) in fins of different geometries. One may observe from Table 1 that at any given point on the radius the temperature values are much higher in a spherical fin than radial and rectangular geometries. It turns out that the spherical fin is not as effective in transferring heat as the radial or rectangular fins.

Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
R. J. Moitsheki is grateful to the National Research Foundation of South Africa for the generous financial support.
References
 M. Pakdemirli and A. Z. Sahin, “Group classification of fin equation with variable thermal properties,” International Journal of Engineering Science, vol. 42, no. 1718, pp. 1875–1889, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. H. Bokhari, A. H. Kara, and F. D. Zaman, “A note on a symmetry analysis and exact solutions of a nonlinear fin equation,” Applied Mathematics Letters, vol. 19, no. 12, pp. 1356–1360, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Pakdemirli and A. Z. Sahin, “Similarity analysis of a nonlinear fin equation,” Applied Mathematics Letters, vol. 19, no. 4, pp. 378–384, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 O. O. Vaneeva, A. G. Johnpillai, R. O. Popovych, and C. Sophocleous, “Group analysis of nonlinear fin equations,” Applied Mathematics Letters, vol. 21, no. 3, pp. 248–253, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. O. Popovych, C. Sophocleous, and O. O. Vaneeva, “Exact solutions of a remarkable fin equation,” Applied Mathematics Letters, vol. 21, no. 3, pp. 209–214, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. Khani, M. A. Raji, and H. HamediNejad, “Analytical solutions and efficiency of the nonlinear fin problem with temperaturedependent thermal conductivity and heat transfer coefficient,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 8, pp. 3327–3338, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. J. Moitsheki, T. Hayat, and M. Y. Malik, “Some exact solutions of the fin problem with a power law temperaturedependent thermal conductivity,” Nonlinear Analysis, Real World Applications, vol. 11, no. 5, pp. 3287–3294, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. J. Moitsheki, “Steady onedimensional heat flow in a longitudinal triangular and parabolic fin,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 10, pp. 3971–3980, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. J. Moitsheki and M. D. Mhlongo, “Classical lie point symmetry analysis of a steady nonlinear onedimensional fin problem,” Journal of Applied Mathematics, vol. 2012, Article ID 671548, 13 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. L. Ndlovu and R. J. Moitsheki, “Analytical solutions for steady heat transfer in longitudinal fins with temperature dependent properties,” Mathematical Problems in Engineering, vol. 2013, Article ID 273052, 14 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 A. Moradi, “Analytical solutions for fin with temperature dependant heat transfer coefficient,” International Journal of Engineering & Applied Sciences, vol. 3, no. 2, pp. 1–12, 2011. View at: Google Scholar
 G. W. Bluman and S. Kumei, Symmetries and Differential Equations, vol. 81 of Applied Mathematical Sciences, Springer, New York, NY, USA, 1989. View at: MathSciNet
 G. W. Bluman and S. Anco, Symmetry and Integration Methods for Differential Equations, vol. 154 of Applied Mathematical Sciences, Springer, New York, NY, USA, 2002. View at: MathSciNet
 G. W. Bluman, A. Cheviakov, and S. Anco, Applications of Symmetry Methods to Partial Differential Equations, vol. 168 of Applied Mathematical Sciences, Springer, New York, NY, USA, 2010. View at: Publisher Site  MathSciNet
 P. J. Olver, Applications of Lie Groups to Differential Equations, vol. 107 of Graduate Texts in Mathematics, Springer, New York, NY, USA, 1986. View at: Publisher Site  MathSciNet
 N. H. Ibragimov, Elementary Lie Group Analysis and Ordinary Differential Equations, vol. 4 of Mathematical Methods in Practice, John Wiley & Sons, Chichester, UK, 1999. View at: MathSciNet
 F. M. Mahomed, “Symmetry group classification of ordinary differential equations: survey of some results,” Mathematical Methods in the Applied Sciences, vol. 30, no. 16, pp. 1995–2012, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Abramowitz and I. A. Stegun, Handbook of Mathematicsl Functions, Dover, New York, NY, USA, 1972.
Copyright
Copyright © 2014 K. J. Moleofane and R. 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.