Abstract

One-dimensional steady-state heat transfer in fins of different profiles is studied. The problem considered satisfies the Dirichlet boundary conditions at one end and the Neumann boundary conditions at the other. The thermal conductivity and heat coefficients are assumed to be temperature dependent, which makes the resulting differential equation highly nonlinear. Classical Lie point symmetry methods are employed, and some reductions are performed. Some invariant solutions are constructed. The effects of thermogeometric fin parameter, the exponent on temperature, and the fin efficiency are studied.

1. Introduction

Heat transfer through extended surfaces has been studied quite extensively [1], perhaps because of its frequent applications in engineering. Through the process of mathematical modeling, heat transfer problems are reduced to nonlinear differential equations.

Accurate and efficient exact, analytical, and approximate schemes for solving differential equations have been devised through considerable effort, particularly those arising in heat conduction through one-dimensional fin problems (see, e.g., [28]). The obtained solutions include series solutions [3, 4, 7], homotopy methods [2], and differential transformation methods (approximate analytical methods) [9]. Few exact solutions exist for one-dimensional problems. In fact, the existing solutions are constructed only for constant thermal conductivity and heat transfer coefficient. Recently, Moitsheki et al. [10] constructed the exact solutions of the one-dimensional fin problem given nonlinear thermal conductivity and heat transfer coefficient. This work has been extended in [11] whereby the introduction of the Kirchhoff transformation linearized the one-dimensional fin problem when heat transfer is a differential consequence of thermal conductivity.

Symmetry methods have been used to analyze the one-dimensional fin problems with heat transfer coefficient depending on the spatial variable [1216]. However, these analyses excluded real-world applications. In recent years, many authors have been interested in the steady-state problems [24, 10] describing heat flow in one-dimensional longitudinal rectangular fins. The symmetry analysis, in particular, group classification of the unsteady fin problem, has attracted some interests (see, e.g., [1216]). Recently, Moitsheki and Harley [17] considered fins of various profiles with both heat transfer coefficient and thermal conductivity being given as temperature dependent. An analysis of a steady nonlinear one-dimensional fin of a rectangular profile was given by Moitsheki and Mhlongo [18].

An accurate transient analysis provided insight into the design of fins that failed in steady-state operations but worked well for some operating periods [19]. The transient problem is considered for a fin of arbitrary profile in [20]. However, both thermal conductivity and heat transfer are considered to be constants. Transient response of longitudinal rectangular fins to step change in base temperature and in base heat flow conditions was studied by Mhlongo et al. [21].

In this paper, we determine exact solutions of nonlinear fin problem for steady heat transfer in longitudinal fin of various profiles where the thermal conductivity is related to temperature by a power law. In Section 2, we provide the mathematical formulation of the problem. We determine the exact solution using MAPLE in Section 3. A brief description of symmetry analysis is provided in Section 4. In Section 5, we employ the symmetry techniques to determine, wherever possible, the invariant solutions. Some discussions and concluding remarks are given in Sections 6 and 7, respectively.

2. Mathematical Models

We consider a longitudinal one-dimensional fin with a cross-sectional area as shown in Figure 1. The perimeter of the fin is denoted by and the length of fin by . The fin is attached to a fixed base surface of temperature and extends into a fluid of temperature . The fin profile is given by the function and the fin thickness at the base is .

The energy balance for a longitudinal fin is given by where and are the nonuniform thermal conductivity and heat transfer coefficient depending on the temperature, is the temperature distribution, is the fin profile, is time, and is the spatial variable. The fin length is measured from the tip to the base. The prescribed boundary conditions are given by (see, e.g., [1])

Introducing the dimensionless variables and the dimensionless numbers,

reduce (1) to

The dimensionless boundary conditions are given by

Here, is the thermogeometric fin parameter and is the fin thickness, is the dimensionless temperature, is the dimensionless spatial variable, is the dimensionless fin profile, is the dimensionless thermal conductivity, is the thermal conductivity of the fin at the ambient temperature, is the dimensionless heat transfer coefficient, and is the heat transfer coefficient at the fin base. For most industrial applications, the heat transfer coefficient may be given as the power law [22]: where the exponents and are constants. The constant may vary between and . However, in most practical applications, it lies between and [22]. If the heat transfer coefficient is given by (6), then the hypothetical boundary condition (i.e., insulation) at the tip of the fin is taken into account [22]. If the tip is not assumed to be insulated, then the problem becomes overdetermined (see also [23]). This boundary condition is realized for sufficiently long fins [22]. Besides, the heat transfer through the outermost edge of the fin is negligible compared to that which passes through the side [23]. The exponent represents laminar film boiling or condensation when , laminar natural convection when , turbulent natural convection when , nucleate boiling when , and radiation when , and implies a constant heat transfer coefficient. Exact solutions may be constructed for the steady-state one-dimensional differential equation describing temperature distribution in a straight fin when the thermal conductivity is a constant and , , , and [22].

The thermal conductivity of the fin may be assumed to vary nonlinearly with the temperature; that is, The one-dimensional heat balance equation is then given by

Recently, (8) has been analyzed using the differential transform methods (DTM) [9]. A proposition in the work of Ndlovu and Moitsheki in [9] concluded that , in equations such as (8), needs to be given by an exponential or power law with exponent being strictly for DTM to work successfully. Here, we employ basic integration and Lie point symmetry techniques.

2.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 was at the base temperature and it is given by (see, e.g., [1]) In dimensionless variables, we have

2.2. Heat Flux

Heat flux at the base of the fin is given by Fourier’s law: The total heat flux of the fin is given by [1] In dimensionless variables, we have where the dimensionless parameter is the Biot number.

3. Exact Solutions

In this section, we analyze the governing equation (8), given . In this case, (8) is linearizable by a transformation . Under such a transformation, (8) becomes and the boundary conditions transform to Equation (15) is analyzed for various situations in the next sections and all the solutions in illustrative examples satisfy both the Dirichlet and the Neumann boundary conditions.

3.1. Case: with

As an illustration, we use two examples, and the rest of the exact solutions are listed in Tables 1 and 2.

Example 1. Given , (15) becomes with solution where and , with . Consider the following: The efficiency is given by The temperature distribution along the surface for this profile is depicted in Figures 2 and 3. The fin efficiency as function of the thermogeometric parameter is shown in Figure 4.

Example 2. In case of , then (15) is transformed into
which is solved by Bessel functions, and the solution in terms of original variables is where and
The efficiency is given by
The temperature distribution along the surface for this profile is depicted in Figures 5 and 6. The fin efficiency as function of the thermogeometric fin parameter is shown in Figure 7.

3.2. Case: with

We will use two examples as in the previous case; the rest of the solutions are listed in Tables 3 and 4.

Example 3. Starting with , (15) in its changed form will be
with solution
The efficiency is given by

Example 4. We consider as the second example. This transforms (15) into
with solution where and . Consider the following:
The efficiency is given by

3.3. Case:

The governing equation (8) becomes

After simplification, (32) becomes The solutions to (33) for various are given in Table 9.

Example 5. is considered as an example, and the solution is and satisfying boundary conditions The corresponding fin efficiency is given by Other exact solutions are listed in Tables 5, 6, and 7.

3.4. Case:

Example 6. We consider and and let , which transforms the governing equation (8) into By direct integration using Polyanin and Zaitsev [24], we get
which yields
Using the boundary conditions (16), we get the solution In terms of the original variables, we have
The corresponding fin efficiency is given by Other exact solutions are listed in Tables 8 and 9.

4. Lie Point Symmetry Analysis

The theory and applications of symmetry analysis may be found in excellent text such as those of [2529]. In brief, the symmetry of a differential equation is an invertible transformation of dependent and independent variables which leave the form of the equation in question unchanged [30]. To determine the symmetry for the governing equation (8), one may seek the transformations of the following form: The infinitesimal transformations (43) act on the space with the corresponding infinitesimal generator which leaves the governing equation invariant. We will then apply the boundary condition to the obtained invariant solutions. The action of is extended in the governing equation through the second prolongation given by where

with being the total derivative operator defined by The prime implies differentiation with respect to . The invariance surface condition is given by The coefficients of do not involve derivatives; we can separate (48) with respect to the derivatives of and solve the resulting overdetermined system of linear homogeneous partial differential equations known as determining equations. Further calculations were facilitated by the freely available package Dimsym [31], a subprogram of REDUCE [32].

5. Symmetry Reductions and Invariant Solutions

We employ the direct group classification to calculate the Lie point symmetries admitted by (8). A few cases arise.

5.1. Case:

For and , the symmetries are listed in Table 9.

5.2. Case: , , and

Example 7. As an illustrative example, we consider and . The governing equation (8) becomes
We observe that (49) admits a non-Abelian two-dimensional Lie algebra spanned by the base vectors listed in Table 9. This noncommuting pair of symmetries leads to the canonical variables The corresponding canonical forms of and are Writing transforms (49) to Here, the prime denotes the total derivative with respect to . Three cases arise.

Subcase 1. For , we obtain the constant which is not related to the original problem. Thus, we ignore it.

Subcase 2. If the term in the square bracket vanishes, then we obtain in terms of the original variables the exact “particular” solution which is not physically realistic. Therefore, we ignore this solution.

Subcase 3. Solving the entire equation (52), we obtain the two solutions that satisfy the boundary conditions. Consider the following:
The fin efficiency is given by

Example 8. Given and , , then the governing equation (8) becomes after the substitution . The two-dimensional Lie algebra admitted by (56) is listed in Table 9. These resulting canonical variables are The corresponding canonical forms of and are By writing , (56) is transformed into
As in the previous example, three cases arise.

Subcase 4. For . This is not related to the original problem. Thus, it is not considered.

Subcase 5. If the term in the square bracket vanishes, then we obtain in terms of the original variables the exact “particular” solution which satisfies the boundary condition only at one end.

Subcase 6. Lastly, we solve the entire equation (59) and obtain the solution that satisfies both the Dirichlet and the Neumann boundary conditions. Consider the following:
The fin efficiency is given by
The temperature distribution along the surface for this profile is depicted in Figures 8 and 9. The fin efficiency as function of the thermogeometric fin parameter is shown in Figure 10.

6. Some Discussions

We now analyze fin problem using solutions given in (18) and (21). We observe in Figure 5 that, for the case of laminar film boiling or condensation, the temperature is inversely proportional to the thermogeometric fin parameter. An increase in values of yielded the decrease in values of temperature. Temperature distribution along the surface was studied for varying values of , while was kept constant. The results depicted in Figure 6 show that the temperature is directly proportional to the parameter . The fin efficiency as function of the thermogeometric fin parameter is shown in Figure 7. Similar trends can be observed from the figures showing temperature distribution and efficiency for other profiles.

The Lie commutators or Lie brackets are given in Table 10 and further reductions are provided in Table 11. Solutions for are furnished in Table 12. The solutions in terms of the original variables are listed in Table 13. Most of these exact solutions do not satisfy one of the boundary conditions. Symmetries and further analysis of were ignored in this paper. The solution for for is given in [33]; therefore, we focused on the case where .

7. Concluding Remarks

Exact solutions for fin problem with power law temperature-dependent thermal conductivity and heat transfer coefficient were constructed. Lie symmetry techniques were used in cases where direct integration was not feasible. Results showing longitudinal fin of various profiles were presented. The obtained solutions satisfy the physical boundary conditions. The exact solutions constructed here could be used as benchmarks or validation tests for numerical schemes.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

M. D. Mhlongo would like to thank the Defence, Peace, Safety and Security Landward Sciences in the Council for Scientific and Industrial Research and the University of the Witwatersrand for financial support. R. J. Moitsheki wishes to thank the National Research Foundation of South Africa for the generous financial support.