Abstract
The steady heat transfer in a heat-generating fin with simultaneous surface convection and radiation is studied analytically using optimal homotopy asymptotic method (OHAM). The steady response of the fin depends on the convection-conduction parameter, radiation-conduction parameter, heat generation parameter, and dimensionless sink temperature. The heat transfer problem is modeled using two-point boundary value conditions. The results of the dimensionless temperature profile for different values of convection-conduction, radiation-conduction, heat generation, and sink temperature parameters are presented graphically and in tabular form. Comparison of the solution using OHAM with homotopy analysis method (HAM) and Runge-Kutta-Fehlberg fourth-fifth-order numerical method for various values of controlling parameters is presented. The comparison shows that the OHAM results are in excellent agreement with NM.
1. Introduction
Fins (extended surfaces) are widely used to enhance the heat transfer rate between a hot surface and its surrounding fluid. Fin applications have included the cooling of computer processors, air conditioning units, refrigerators, air-cooled engines, and oil carrying pipelines. In the past three decades, fins have gained vast recognition for cooling electronic tools as heat sinks. The subject of extended surface heat transfer is now a fully developed technology but with continuing contributions from numerous researchers. Background information on heat transfer in extended surfaces may be found in the books [1, 2], where the authors have presented wide-ranging coverage of the various facts of this technology.
Numerous mathematical models related to heat transfer in fins of various shapes with different boundary conditions are well documented in the research literature. For instance, the mathematical analysis of convective fins was first provided by Gardener [3] based on the assumption of constant conductivity and a uniform coefficient of convective heat transfer along the fin surface. Khani et al. [4] presented some exact solutions for 1D fin problem with uniform thermal conductivity and heat transfer coefficient. Khani et al. [5] also provided a series solution for 1D fin problem with constant heat transfer coefficient and temperature dependent thermal conductivity.
A variety of approximate analytical methods have been used to study the transient response of fins. Aziz and Na [6] presented a coordinate perturbation expansion for the response of an infinitely long fin due to a step change in the base temperature. Chang et al. [7] used the methods of optimal linearization and variational embedding, and Campo [8] utilized variational techniques to analyze radiative-convective fins under unsteady operating conditions. Solutions for transient heat transfer were constructed for fins by Onur [9]. Aziz and Torabi [10] have presented the numerical analysis of transient heat transfer in fin with temperature dependent heat transfer coefficient.
Exact steady-state solutions of 2D models of fin having constant thermal conductivity and heat transfer coefficient and with no internal heat generation were analyzed in [11–17]. The addition of internal heat generation function based with spatial dependence is discussed [18–20].
In this paper, we used a new approximate method, namely, optimal homotopy asymptotic method [21–27] for steady-state heat transfer with internal heat generation fin, and investigated numerically the effects of the different governing parameters on dimensionless temperature profile in a nonlinear fin-type problem. For comparison purposes the governing highly nonlinear problem is also solved using Runge-Kutta-Fehlberg fourth-fifth-order method and homotopy analysis method (HAM) developed by Liao [28].
The paper is planned as follows: in Section 2 we formulate our nonlinear problem, basic principles of OHAM are discussed in Section 3, solution of the problem via OHAM is presented in Section 4, and Section 5 is reserved for results and discussion. Conclusions are drawn in Section 6.
2. Mathematical Formulation
Consider a straight fin of constant cross-sectional area (rectangular, cylindrical, elliptic, etc.), perimeter of the cross-section , and length as shown in Figure 1. The fin has a thermal conductivity and a thermal diffusivity . The surface of the fin behaves as a gray diffuse surface with an emissivity . The fin is deemed to be initially in thermal equilibrium with the surroundings at temperature . Its tip is insulated. A volumetric internal heat generation rate occurs in the fin. The fin loses heat by simultaneous convection and radiation to its surroundings at temperature . The same sink temperature is used for both convection and radiation to avoid the introduction of an additional parameter in the problem.
For one-dimensional steady conduction in the fin, the energy equation may be written as The initial and boundary conditions are where is measured from the tip of the fin with the introduction of the following definitions: Equations (2) and (3) can be written in dimensionless form as follows: The instantaneous base heat flow is given by: which may be expressed in dimensionless form as follows. The instantaneous convective heat loss from the fin is given by or in dimensionless form as Similarly, the instantaneous radiative heat loss from the fin can be obtained as or in dimensionless form as The instantaneous total surface heat loss in dimensionless form is the sum of convective and radiative losses given by (11) and (13); that is, The instantaneous rate of energy storage in the fin can be calculated from the energy balance as follows: or in dimensionless form as where
3. Basic Principles of OHAM
We review the basic principles of OHAM as expounded by Marinca et al. [21–24] as well as other researchers including [25, 26]. (i) Let us consider the following differential equation: where is problem domain, , where and are linear and nonlinear operators, is an unknown function, and is a known function. (ii) Construct an optimal homotopy equation as where is an embedding parameter, and is auxiliary function on which the convergence of the solution greatly depends. The auxiliary function also adjusts the convergence domain and controls the convergence region.(iii) Expand in Taylor’s series about . One has an approximate solution: Many researchers have observed that the convergence of the series equation (19) depends upon , . If it is convergent, then we obtain (iv) Substituting (20) in (17), we have the following residual: If , then will be the exact solution. For nonlinear problems, generally this will not be the case. For determining , Galerkin’s Method, Ritz Method, or the method of least squares can be used.(v) Finally, substitute these constants in (21), and one can get the approximate solution.
4. OHAM Solution for Heat-Generating Fin
According to the OHAM, (1) can be written as where prime denotes differentiation with respect to .
We consider and as follows: Using (23) in (22) and after some simplifying and rearranging the terms based on the powers of , we obtain the zeroth-, first-, and second-order problems as follows.
The zeroth-order problem is with boundary conditions Its solution is The first-order problem is with boundary conditions having solution The second-order problem is with boundary conditions It is given by The second-order approximate solution by OHAM for is We use the method of least squares to obtain the unknown convergent constant in .
In particular case, and , for , , , and .
By considering the values of in (33) and after simplifying, the second-order approximate analytical OHAM solution can be obtained
5. Results and Discussion
Equation (4) shows that fin temperature is based on four parameters: , , , and which govern this highly nonlinear second-order differential equation. The effect of each parameter on fin temperature is tabulated and graphically presented for different values of the controlling parameters.
In order to validate the accuracy of our approximate solution via OHAM, we have presented a comparative study of OHAM solution with homotopy analysis method (HAM) and numerical solution (Runge-Kutta-Fehlberg fourth-fifth-order method). Table 1 has been prepared to exhibit the comparison of dimensionless temperature obtained by OHAM, homotopy analysis method (HAM), and the numerical method (NM) for several values of heat-generating parameter , when other parameters are fixed. It is observed that, with increasing values of internal heat-generating parameter , the temperature profile gradually increases. Clearly the OHAM solutions are very close to the numerical solution as compared to HAM. This can be seen from the percentage error in the dimensionless temperature obtained by OHAM, HAM, and NM. The increase in dimensionless temperature is also evident in Table 2, in which we have used different values of sink temperature parameter , and other parameters values are predetermined. From Tables 1 and 2, it is observed that our OHAM solutions are more accurate than HAM; this confirms that OHAM is more consistent with approximate analytical method than with HAM. The major factor in HAM is its computational time for finding the (h curve), while in OHAM the ensuring convergence of the solution depends on parameters , which are optimally determined, resultantly HAM in more time consuming than OHAM.
In Table 3, we show the comparison of dimensionless temperature obtained by OHAM and the numerical method (NM) for several values of convection parameter , while other parameters are kept unchanged. It is observed that, with the increase of , the temperature profile shows decrease, and the same phenomena of decrease in dimensionless temperature can be observed in Table 4 for different values of radiation parameter , when the other parameters values are fixed.
In Figures 2, 3, 4, and 5 we depict the dimensionless temperature profile and its variation for different values of parameters. It is important to note that the dimensionless temperature increases with each controlling parameter.
6. Conclusion
We have successfully applied the optimal homotopy asymptotic method for the approximate solution of steady state of heat-generating fin with simultaneous surfaces convection and radiation. The effects of radiation parameter , convection parameter , internal heat-generating parameter , and the sink temperature parameter on temperature profile in the fin are investigated analytically. It is observed that dimensionless fin temperature profile is dependent on the four parameters , , , and . Comparison for the dimensionless temperature has been made between the solutions obtained using OHAM with HAM and Runge-Kutta-Fehlberg fourth-fifth-order method. It is found that the OHAM solution is very close to the numerical solution than HAM, which reveals the reliability and efficiency of OHAM. Approximate analytical solution to highly nonlinear problem was achieved without any assumption of linearization, and we can extend this approach to a variety of nonlinear heat transfer problems.