Abstract
Exact solutions for models describing heat transfer in a two-dimensional rectangular fin are constructed. Thermal conductivity, internal energy generation function, and heat transfer coefficient are assumed to be dependent on temperature. We apply the Kirchoff transformation on the governing equation. Exact solutions satisfying the realistic boundary conditions are constructed for the resulting linear equation. Symmetry analysis is carried out to classify the internal heat generation function, and some reductions are performed. Furthermore, the effects of physical parameters such as extension factor (the purely geometric fin parameter) and Biot number on temperature are analyzed. Heat flux and fin efficiency are studied.
1. Introduction
Fins are extended surfaces used to increase the heat transfer rate between a hot body and its surroundings. There are a variety of uses such as in air conditioning systems, compressors, and cooling of electronic components. The theory on heat transfer in extended surface may be found in texts such as [1, 2]. Few exact solutions exist even for the one dimensional fin problem with constant thermal conductivity and heat transfer coefficient [3]. Series solutions for one-dimensional fin problem with constant heat transfer coefficient and temperature-dependent thermal conductivity are given in [4]. Furthermore, analytical and exact solutions for one dimensional fins models with temperature-dependent thermal conductivity and heat transfer coefficient were obtained for example, in [5β7]. A compendium of heat transfer in all types of one dimensional fins is given in [8]. Exact steady-state solutions exist for two-dimensional models with constant thermal conductivity and heat transfer coefficient, with no internal heat generation [9β15], and with internal heat generation function depending on a spatial variable [16, 17]. Solutions for transient heat transfer in fins are constructed in [18].
In this paper, we apply the Kirchoff transformation to the two-dimensional steady-state heat balance equation in Section 2. The resulting model contains the arbitrary source term. In Section 3, we construct exact solutions which satisfy the physical boundary conditions. Symmetry analysis is carried out for the single governing equation, wherein group classification is performed in Section 4. Heat transfer results including fin efficiency and heat flux are analyzed in Section 5. Lastly, we provide concluding remarks in Section 6.
2. Mathematical Model
We consider a two-dimensional rectangular fin of length as shown in Figure 1. The fin is mounted to a base surface of temperature and extended into its surrounding of temperature . The heat flow is assumed to be symmetric along the line . We assume that the heat transfer coefficient along the fin is nonuniform and temperature dependent. Also, the internal heat generation is nonzero and temperature dependent.
The two-dimensional heat balance equation is given by (see, e.g., [16, 17]) The imposed boundary conditions are
Here, is the dimensionless temperature, is the fin base temperature, is the heat transfer coefficient, is the longitudinal coordinate, is the transverse coordinate, is the internal heat generation function, and is the thermal conductivity. Several authors have considered the two-dimensional problem with and thermal conductivity being a constant (see, e.g., [19, 20]) and the case with a temperature-dependent thermal conductivity [21].
Introducing the dimensionless variables
we obtain
The corresponding dimensionless boundary conditions are where is the fin extension factor (purely geometric parameter), and and are the Biot numbers. is the reciprocal to aspect ratio (see, e.g., [16]). and are the heat transfer at the base and thermal conductivity of the fin at the ambient temperature, respectively.
3. Exact Solutions
In this section, we construct exact solutions for the boundary value problem (BVP) (2.4)β(2.5). The problem is simplified by the introduction of the Kirchoff transformation. The Kirchoff's transformation (see, e.g., [21]) with being an arbitrary constant, reduces the BVP (2.4)β(2.5) to subject to the conditions
We consider two cases of the thermal conductivity. Note that in (3.3) and (3.4), we required the product of the heat transfer coefficient and temperature to match the integral of thermal conductivity. In fact, in a one-dimensional case, must be the differential consequence of (see, e.g., [7]). is quadratic in when is linear, and is given by the power law when is nonlinear.
3.1. Case 1. Linear Thermal Conductivity
Thermal conductivity is assumed to be a linear function of temperature for many engineering applications [3]. We assume thermal conductivity to be linear function of temperature (see also [21, 22])
is the parameter that describes temperature dependency [3, 4]. In dimensionless variables, where . This case of requires , so that BVP (3.2)β(3.6) hold. Note that renders heat transfer coefficient to be singular. However, one can remove singularity by choosing, without loss of generality, . Assuming the internal heat generation to be linearly dependent on temperature, then the governing equation becomes the modified Helmholtz type equation We seek exact solutions of (3.8) subject to (3.3)β(3.6). Using method of separation of variables, we obtain the nontrivial exact solutions for two cases of the separation constant . Note that leads to trivial solutions.
3.1.1.
where is an arbitrary constant and satisfies . Also,
3.1.2.
where is satisfied by The temperature profile for solution (3.11) is given in Figures 2(a) and 2(b). In Figures 3(a) and 3(b), we plot the temperature profile at the boundaries and , respectively. For simplicity, we allowed and to be equal. We list the first five eigenvalues of (3.12) for various values of the Biot number in Table 2.
(a) Parameters used are from bottom to top sheet, respectively. is fixed at 2 |
(b) Parameters used are from top to bottom sheet, respectively. Both and are fixed at 0.2 |
(a) Temperature profile at . Parameters used are |
(b) Temperature profile at . Parameters used are |
In terms of the original temperature variable, we obtain the exact solution
where .
3.2. Case 2: Nonlinear Thermal Conductivity
Assuming the power law temperature-dependent thermal conductivity
which is given in dimensionless variables as , requires
Again, one may, without loss of generality, assume, . We obtain in terms of original variables the exact solution
where is given in Sections 3.1.1 and 3.1.2.
4. Symmetry Analysis
In this section, we analyze (3.2) with an arbitrary source term. We employ symmetry techniques. A symmetry group of a system of differential equations is a group of transformations mapping any arbitrary solution to another solution of the system. Such groups depend on continuous parameters. Given a continuous one-parameter symmetry group, in most practical cases, one may reduce the number of independent variables by one. The most familiar symmetry is the rotational symmetry that enables one to reduce the variables to a single radial variable . Recent accounts on this theory may be found in many excellent texts such as those of [23, 24]. We use the methods in [24] (which exclude the use of equivalence transformations) to determine possible cases of for which extra symmetries are obtained.
In the initial symmetry analysis where is arbitrary, we obtain the translations of the spatial variables and and the rotational symmetry. Cases for which extra symmetries are obtained are listed in Table 1. Note that reduces (3.2) to the steady two-dimensional thermal diffusion equation (the Laplace equation which has been already analyzed).
4.1. Symmetry Reductions: Illustrative Example
Given the internal heat generation term as the power law (as listed in Table 1), leads to the functional form
where satisfies a nonlinear ordinary differential equation (ODE) The ODE (4.2) is harder to solve exactly. Furthermore, the boundary conditions are not invariant under . In fact, the BVP (2.4)β(2.5) is not invariant under all the admitted symmetries listed. One may seek numerical solutions when internal heat generation term is nonlinear. We omit numerical analysis in this paper. However, the obtained exact solutions in Section 3 may be used as benchmarks for the numerical schemes.
5. Heat Transfer Results
The number of eigenvalues required to calculate the temperature distribution, heat flux, and fin efficiency accurately depends on the Biot number . We observe in Table 2 below that Biot number is directly proportional to the eigenvalues. Similar results are obtained for heat transfer in orthotropic convective pin fin [14]. The expression for the temperature distribution is given explicitly in (3.9) and (3.11). However, in further analysis, we focus on solution (3.11). The temperature distribution depends on a number of variables including , , eigenvalues, and the arbitrary function of describing the temperature at the base of the fin. We are free to choose any function . Note that temperature distribution is proportional to both and .
We observe in Figure 2(a) that the Biot number is directly proportional to the temperature distribution in the fin. Also, in Figure 2(b), we notice that temperature decreases with the increase of the extension factor. Clearly, if the length of the fin is increased, temperature is at its lowest value, or increased width of the fin results in increased temperature distribution. Figures 3(a) and 3(b) depict the temperature distribution at the boundaries of the fin. We note that, in Figure 3(a), there is a significant reduction in temperature along the boundary and toward the tip of the fin (one may recall that we assumed that the fin is measured from the tip to the base). Figure 3(b) shows temperature variation at the tip of the fin.
5.1. Fin Efficiency and Heat Flux
5.1.1. Heat Flux
The heat transfer from the fin base may be constructed by evaluating heat conduction rate at the base (see, e.g., [12])
The dimensionless heat transfer rate from the base of the fin is defined by [12]
5.1.2. Fin Efficiency
Fin efficiency (overall fin performance) is defined as the ratio of the actual heat transfer from the fin rate of heat that would be ideally transferred if the entire fin was at the temperature of the fin base [25]. The local fin efficiency is defined by [13] or simply
5.1.3. Flux and Fin Efficiency Given (3.11)
Given the solution (3.11), we obtain heat flux in terms of , and fin efficiency
Fin efficiency (5.6) and heat flux (5.5) are depicted in Figures 4(a) and 4(b), respectively. Wherein both heat flux and fin efficiency are plotted against the extension factor. Here, is fixed at 0.2. In Figure 4(a), we observe that fin performance decreases with increased extension factor. Moreover, the increased Biot number yields decreased fin efficiency. In our entire analysis, we have assumed a nonuniform internal heat generation. Internal heat generation function is assumed to be a uniform in [26]. In Figure 4(b), we observe that Biot number is inversely proportional to the heat flux. However, heat flux is directly proportional to the extension factor (implying that longer fins result in higher heat flux (see also [14])). The observations in Figures 4(a) and 4(b) are consistent with the results in the literature (see, e.g., [26, 27] or Chapter 15 in [1]). In [26], finite elements methods were used to determine, among others, the effects of uniform internal heat generation. Here, internal energy generation is given as linear function of temperature, and base temperature is quadratic in . Sikka and Iqbal [27] provided the fin efficiency for two-dimensional pin fin. The observations in Figure 4(b) are consistent with those of [27].
(a) Fin efficiency
(b) Heat flux
6. Concluding Remarks
We have successfully applied the Kirchoffβs transformation to reduce the nonlinearity of (2.4) (the source term remained arbitrary). Exact solutions for two-dimensional rectangular fin with temperature-dependent thermal conductivity and heat transfer coefficient, furthermore with internal energy generation function, which depend linearly on temperature are constructed. In the analysis, we allowed the temperature of the fin base to be quadratic in . The forms of the internal energy generation term for which extra symmetries are admitted were obtained. Reduction of the single PDE (given nonlinear source term) to the ODE is achieved. As far as we know, symmetry methods have not yet been employed to two-dimensional fin problems. However, the entire BVP is not invariant under the admitted symmetries. Constructed exact solutions have provided insight into the heat transfer processes in a rectangular straight fin and may be used as benchmarks for the numerical schemes, particularly when thermal conductivity, heat transfer coefficient, and internal thermal energy generation function are all temperature dependent.
Nomenclature
: | Temperature distribution |
: | Fin base temperature |
: | Surrounding temperature |
: | Dimensionless temperature |
: | Transformed temperature variable |
: | Thermal conductivity |
: | Thermal conductivity of the fin at the ambient temperature |
: | Dimensionless thermal conductivity |
: | Heat transfer coefficient |
: | Heat transfer coefficient at the fin base |
: | Dimensionless heat transfer coefficient |
: | Transverse coordinate |
: | Longitudinal coordinate |
: | Spatial variable |
: | Spatial variable |
: | Length of the fin |
: | Half width |
: | Internal energy generation function |
: | Dimensionless internal energy generation function |
: | Biot number |
: | Biot number |
: | Extension factor (purely geometric parameter) |
: | Heat flux at the fin base |
: | Heat flux |
: | Overall heat flux for an ideal fin |
: | Fin efficiency |
: | Constants given in the paper. |