Solution of Temperature Distribution in a Radiating Fin Using Homotopy Perturbation Method
Radiating extended surfaces are widely used to enhance heat transfer between primary surface and the environment. The present paper applies the homotopy perturbation to obtain analytic approximation of distribution of temperature in heat fin radiating, which is compared with the results obtained by Adomian decomposition method (ADM). Comparison of the results obtained by the method reveals that homotopy perturbation method (HPM) is more effective and easy to use.
Most scientific problems and phenomena such as heat transfer occur nonlinearly. Except a limited number of these problems, it is difficult to find the exact analytical solutions for them. Therefore, approximate analytical solutions are searched and were introduced [1–5], among which homotopy perturbation method (HPM) [6–12] and Adomain decomposition method (ADM) [13, 14] are the most effective and convenient ones for both weakly and strongly nonlinear equations.
The analysis of space radiators, frequently provided in published literature, for example, [15–22], is based upon the assumption that the thermal conductivity of the fin material is constant. However, since the temperature difference of the fin base and its tip is high in the actual situation, the variation of the conductivity of the fin material should be taken into consideration and includes the effects of the variation of the thermal conductivity of the fin material. The present analysis considers the radiator configuration shown in Figure 1. In the design, parallel pipes are joined by webs, which act as radiator fins. Heat flows by conduction from the pipes down the fin and radiates from both surfaces.
Here, the fin problem is solved to obtain the distribution of temperature of the fin by homotopy perturbation method and compared with the result obtained by the Adomian decomposition method, which is used for solving various nonlinear fin problems [23–25].
2. The Fin Problem
A typical heat pipe space radiator is shown in Figure 1. Both surfaces of the fin are radiating to the vacuum of outer space at a very low temperature, which is assumed equal to zero absolute. The fin is diffuse-grey with emissivity , and has temperature-dependent thermal conductivity , which depends on temperature linearly. The base temperature of the fin and tube surfaces temperature is constant; the radiative exchange between the fin and heat pipe is neglected. Since the fin is assumed to be thin, the temperature distribution within the fin is assumed to be one-dimensional. The energy balance equation for a differential element of the fin is given as : where and are the thermal conductivity and the Stefan-Boltzmann constant, respectively. The thermal conductivity of the fin material is assumed to be a linear function of temperature according to where is the thermal conductivity at the base temperature of the fin and is the slope of the thermal conductivity temperature curve.
Employing the following dimensionless parameters: the formulation of the fin problem reduces to with boundary conditions
3. Basic Idea of Homotopy Perturbation Method
In this study, we apply the homotopy perturbation method to the discussed problems. To illustrate the basic ideas of the method, we consider the following nonlinear differential equation, where is defined as follows: where stands for the linear and for the nonlinear part. Homotopy perturbation structure is shown as the following equation: Obviously, using (3.3) we have where is an embedding parameter and is the first approximation that satisfies the boundary condition. We Consider and as
4. The Fin Temperature Distribution
By increasing the number of the terms in the solution, higher accuracy will be obtained. Since the remaining terms are too long to be mentioned in here, the results are shown in tables.
The coefficient representing the temperature at the fin tip can be evaluated from the boundary condition given in (2.4c) using the numerical method.
Tables 1 and 2 show the dimensionless tip temperature, that is, coefficient , for different thermal conductivity parameters, . The tables state that the convergence of the solution for the higher thermo-geometric fin parameter, , is faster than the solution with lower fin parameter. It is clear from the tables that the solution is convergent. In order to investigate the accuracy of the homotopy solution, the problem is compared with decomposition solution , also the corresponding results are presented in Figure 2. It should be mentioned that the homotopy results in the tables are arranged for first 5, 7, 10, 14 terms of the solution (). It is seen that the results by homotopy perturbation method, and adomian decomposition method are in good agreement. The results of the comparison show that the difference is 3.1% in the case of the strongest nonlinearity, that is, and .
In this work, homotopy perturbation method has been successfully applied to a typical heat pipe space radiator. The solution shows that the results of the present method are in excellent agreement with those of ADM and the obtained solutions are shown in the figure and tables. Some of the advantage of HPM are that reduces the volume of calculations with the fewest number of iterations, it can converge to correct results. The proposed method is very simple and straightforward. In our work, we use the Maple Package to calculate the functions obtained from the homotopy perturbation method.
A. Aziz and T. Y. Na, Perturbation Methods in Heat Transfer, Series in Computational Methods in Mechanics and Thermal Sciences, Hemisphere, Washington, DC, USA, 1984.View at: MathSciNet
J.-H. He, “Homotopy perturbation method for bifurcation of nonlinear problems,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 6, no. 2, pp. 207–208, 2005.View at: Google Scholar
G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, vol. 60 of Fundamental Theories of Physics, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1994.View at: MathSciNet
J. G. Bartas and W. H. Sellers, “Radiation fin effectiveness,” Journal of Heat Transfer, vol. 82C, pp. 73–75, 1960.View at: Google Scholar
J. E. Wilkins, Jr., “Minimizing the mass of thin radiating fins,” Journal of Aerospace Science, vol. 27, pp. 145–146, 1960.View at: Google Scholar