Thermoelastic study of a functionally graded annular fin with variable thermal parameters using semiexact solution

Abstract

In this paper, the thermoelastic behavior of a functionally graded material (FGM) annular fin is investigated. The material properties of the annular fin are assumed to vary radially. The heat transfer coefficient and internal heat generation are considered to be functions of temperature. A closed form solution of nonlinear heat transfer equation for the FGM fin is obtained using the homotopy perturbation method (HPM) which leads to nonuniform temperature distributions within the fin. The temperature field is then coupled with the classical theory of elasticity and the associated thermal stresses are derived analytically. For the correctness of the present closed form solution for the stress field, the results are compared with the ANSYS-based finite element method (FEM) solution. The present HPM-based closed form solution of the stress field exhibits a good agreement with the FEM results. The effect of various thermal parameters such as the thermogeometric parameter, conduction-radiation parameter, internal heat generation parameter, coefficient of variation of thermal conductivity, and the coefficient of thermal expansion on the thermal stresses are discussed. The results are presented in both nondimensional and dimensional form. The dimensional stress analysis discloses the suitability of FGM as the fin material in practical applications.     

Determination of temperature distribution for annular fins with temperature‐dependent thermal conductivity

In this paper, homotopy analysis method (HAM) has been used to evaluate the temperature distribution of annular fin with temperature‐dependent thermal conductivity and to determine the temperature distribution within the fin. This method is useful and practical for solving the nonlinear heat transfer equation, which is associated with variable thermal conductivity condition. HAM provides an approximate analytical solution in the form of an infinite power series. The annular fin heat transfer rate with temperature‐dependent thermal conductivity has been obtained as a function of thermo‐geometric fin parameter and the thermal conductivity parameter describing the variation of the thermal conductivity. © 2011 Wiley Periodicals, Inc. Heat Trans Asian Res; Published online in Wiley Online Library ( wileyonlinelibrary.com/journal/htj ). DOI 10.1002/htj.20353     

Entropy generation in a hollow cylinder with temperature dependent thermal conductivity and internal heat generation with convective–radiative surface cooling

This article investigates entropy generation in an asymmetrically cooled hollow cylinder with temperature dependent thermal conductivity and internal heat generation. The inside surface of the cylinder is cooled by convection on its inside surface while the outside surface experiences simultaneous convective–radiative cooling. The thermal conductivity of the cylinder as well as the internal heat generation within the cylinder are linear functions of temperature, introducing two nonlinearities in the one-dimensional steady state heat conduction equation. A third nonlinearity arises due to radiative heat loss from the outside surface of the cylinder. The nonlinear system is solved analytically using the differential transformation method (DTM) to obtain the temperature distribution which is then used to compute local and total entropy generation rates in the cylinder. The accuracy of DTM is verified by comparing its predictions with the analytical solution for the case of constant thermal conductivity and constant internal heat generation. The local and total entropy generations depend on six dimensionless parameters: heat generation parameter Q, thermal conductivity parameter β, conduction–convection parameters Nc₁ and Nc₂, conduction–radiation parameter Nr, convection sink temperature δ and radiation sink temperature η.     

A decomposition method for fin efficiency of convective straight fins with temperature-dependent thermal conductivity

The Adomian decomposition method (ADM) has been used to evaluate the efficiency of straight fins with temperature-dependent thermal conductivity and to determine the temperature distribution within the fin. The method is useful and practical for solving the nonlinear heat diffusion equation, which is associated with variable thermal conductivity condition. The ADM provides an analytical solution in the form of an infinite power series. The fin efficiency of the straight fins with temperature-dependent thermal conductivity has been obtained as a function of thermo-geometric fin parameter and the thermal conductivity parameter describing the variation of the thermal conductivity. It has been observed that the thermal conductivity parameter has a strong influence over the fin efficiency. The data from the present solutions has been correlated for a wide range of thermo-geometric fin parameter and the thermal conductivity parameter. The resulting correlation equations can assist thermal design engineers for designing of straight fins with temperature-dependent thermal conductivity.     

Thermal performance of straight porous fin with variable thermal conductivity under magnetic field and radiation effects

The present numerical study reports the thermal performance of the straight porous fin with temperature-dependent thermal conductivity, radiation, and magnetic field effects. The heat transfer model comprising the Darcy's law for simulating flow with solid-fluid interactions in porous medium, Rosseland approximation for heat transfer through radiation, Maxwell equations for magnetic field effect and linearly varying temperature dependent thermal conductivity, results into highly nonlinear ordinary differential equation. The governing equation is solved using a finite difference scheme with suitable boundary conditions. The obtained solutions are physically interpreted by considering the impact of different nondimensional parameters on thermal performance, efficiency, and effectiveness of the system through plotted graphs. A detailed result with regard to the Nusselt number at the fin base is calculated. The results obtained are observed to be in excellent agreement with previous studies. From the study, it is observed that there is a significant effect on the thermal performance of the fin in the presence of porous constraints; also, results reveal that the nonlinear thermal conductivity parameter strengthens the thermal performance, efficiency, and effectiveness of the fin. Furthermore, the results of the study reveal that the rate of heat transfer of the fin increases with the increase in the magnetic parameter and radiation parameter.     

Free Convective Flow of Electrically Conducting and Viscous Immiscible Fluid Flow in a Vertical Channel in the Presence of First‐Order Chemical Reaction

In this paper, the effects of chemical reaction on free convective flow of electrically conducting and viscous incompressible immiscible fluids are analyzed. The coupled nonlinear equations governing the heat and mass transfer are solved analytically and numerically with appropriate boundary conditions for each fluid and the solutions have been matched at the interface. The analytical solutions are solved by using regular perturbation method valid for small values of perturbation parameter and numerically by using finite difference method. The numerical results for various values of thermal Grashof number, mass Grashof number, Hartman number, viscosity ratio, width ratio, conductivity ratio, and chemical reaction parameter have been presented graphically in the presence and in the absence of electric field load parameter. In addition, the closed form expression for volumetric flow rate, Nusselt number, species concentration, and total heat rate added to the flow is also analyzed. The solutions obtained by finite difference method and perturbation method agree very well to the order of 10^?4 for small values of perturbation parameter.     

Inverse and efficiency of heat transfer convex fin with multiple nonlinearities

In this article, we first propose the novel semi‐analytical technique—modified Adomian decomposition method (MADM)—for a closed‐form solution of the nonlinear heat transfer equation of convex profile with singularity where all thermal parameters are functions of temperature. The longitudinal convex fin is subjected to different boiling regimes, which are defined by particular values of n (power index) of heat transfer coefficient. The energy balance equation of the convex fin with several temperature‐dependent properties are solved separately using the MADM and the spectral quasi‐linearization method. Using the values obtained from the direct heat transfer method, the unknown parameters of the profile, such as thermal conductivity, surface emissivity, heat generation number, conduction‐convection parameter, and radiation‐conduction parameter are inversely predicted by an inverse heat transfer analysis using the simplex search method. The effect of the measurement error and the number of measurement points has been presented. It is found that present measurement points and reconstruction of the exact temperature distribution of the convex fin are fairly in good agreement.     

Hybrid differential transformation and finite difference method to annular fin with temperature-dependent thermal conductivity

A hybrid numerical technique which combines the differential transformation and finite difference method is utilized to investigate the annular fin with temperature-dependent thermal conductivity. The exposed surfaces of the fin dissipate heat to the surroundings by convection and radiation. The influences of the convective heat transfer coefficient, absorptivity, emissivity and thermal conductivity parameter on the temperature distribution are examined. The results show that the convective heat transfer plays a dominant role for heat dissipation under the convection–radiation condition. The optimum radii ratio of fin which maximizes the heat transfer rate and fin efficiency is also discussed.     

Thermal performance of continuously moving radiative–convective fin of complex cross-section with multiple nonlinearities

Spectral collocation method (SCM) is adopted to predict the temperature distribution in the fin with temperature dependent thermal conductivity, heat transfer coefficient and surface emissivity. These temperature dependent properties or parameters cause multiple nonlinearities of energy equation. In order to reduce these multiple nonlinearities, a local linearization approach is adopted. The spatial distribution of dimensionless temperature is discretized by Lagrange interpolation polynomials. Accordingly, the differential form of energy equation is transformed to the matrix form of algebraic equation. The accuracy of the SCM model is verified by comparing SCM results with available data in references. Meanwhile, compared with analytical solutions, it can be found that the convergence rate of SCM approximately follows exponential law. In addition, effects of various physical parameters, such as Peclet number, thermal conductivity parameter, emissivity parameter, parameter of heat transfer coefficient, convective–conductive parameter and radiative–conductive parameter on the dimensionless temperature, the dimensionless fin-tip temperature and the volume adjusted fin efficiency are comprehensively analyzed.     

Homotopy Perturbation Method for Thermal Stresses in an Annular Fin with Variable Thermal Conductivity

An approximate analytical solution method for thermal stresses in an annular fin with variable thermal conductivity is presented. Homotopy perturbation method (HPM) is employed to estimate the non-dimensional temperature field by solving nonlinear heat conduction equation. The closed-form solutions for the thermal stresses are formulated using the classical thermoelasticity theory coupled with HPM solution for temperature field. The plane state of stress conditions are considered in this study. The effects of thermal parameters such as variable thermal conductivity parameter (β), thermogeometric parameter (K), and the non-dimensional coefficient of thermal expansion (χ) on the temperature field and stress field are studied. The results for temperature field and stress field obtained from HPM-based solution are found to be in very close agreement with the results available in literature. Furthermore, the HPM solution is found to be very efficient and handles nonlinear heat transfer equation with greater convenience.     

Analytical Solution and Numerical Simulation for One‐Dimensional Steady Nonlinear Heat Conduction in a Longitudinal Radial Fin with Various Profiles

In this article we consider a model describing the temperature profile in a longitudinal fin with rectangular, concave, triangular, and convex parabolic profiles. Both thermal conductivity and the heat transfer coefficient are assumed to be temperature‐dependent, and given by a linear function and by power laws, respectively. In addition, the effects of the thermal conductivity gradient have been investigated. Optimal homotopy analysis method (OHAM) is employed to analyze the problem. The effects of the physical applicable parameters such as thermo‐geometric fin, thermal conductivity, and heat transfer mode are analyzed. The OHAM solutions are obtained and validity of obtained solutions is verified by the Runge–Kutta fourth‐order method and numerical simulation. A very good agreement is found between analytical and numerical results. Also for investigation of lateral effects on the accuracy of results, numerical simulation (by Ansis software) is compared with the homotopy analysis method (HAM) and numerical solution (by Runge–Kutta) of the energy balance equation. © 2013 Wiley Periodicals, Inc. Heat Trans Asian Res; Published online in Wiley Online Library ( wileyonlinelibrary.com/journal/htj ). DOI 10.1002/htj.21104     

Thermal distribution analysis of rectangular fin considering multiboiling heat transfer modes

In this investigation, a numerical method is used to compute the thermal distribution analysis of a rectangular fin with surface emissivity and internal heat generation. Here, the thermal conductivity, heat generation, emissivity at the surface, and coefficient of heat transfer depend on temperature linearly. The role of four distinct multiboiling heat transfer modes such as laminar film boiling (condensation), laminar convection, turbulent convection, and nucleate boiling are discussed in detail and the corresponding outcomes are displayed graphically. Isolated (insulated) and convective tip boundary conditions for the fin tip are employed in this study. The solution is obtained using shooting technique involving Runge Kutta Fehlberg method. It is emphasized that the thermal distribution shows a diminishing trend for the convective tip condition compared to the insulated tip. In addition to this, it is illustrated that laminar film boiling and laminar convection are two effective modes of heat transfer in comparison with turbulent convection and nucleate boiling for a finned surface in boiling liquids. The study on fin efficiency shows that fin efficiency increases with the increase in internal heat generation number.     

Stresses in radiative annular fin under thermal loading and its inverse modeling using sine cosine algorithm (SCA)

This work presents a novel mathematical model for the analysis of thermal stresses in a radiative annular fin with temperature-dependent thermal conductivity and radiative parameter. An approximate analytical solution for thermal stresses is derived using a homotopy perturbation method (HPM)-based closed-form solution of steady-state nonlinear heat transfer equation, coupled with classical elasticity theory. The effect of thermal parameters on the temperature field and the thermal stress fields are discussed. The various thermal parameters, such as a parameter describing the temperature-dependent thermal conductivity, coefficient of thermal expansion, coefficient of radiative parameter, and the variable radiative parameter, are inversely estimated for a given stress field. For inverse modeling, a population-based sine cosine algorithm (SCA) was employed to estimate the thermal parameters. The inverse modeling is verified by using the estimated thermal parameters in the closed-form solution of stress field. The reconstructed stress fields obtained from the inversely estimated parameters are then compared with the reference stress field. Results show a very good agreement between the reference stress field and the inversely estimated stress fields.     

Efficiency of the Horizontal Single Pin Fin Subjected to Free Convection and Radiation Heat Transfer

In this communication, the results of numerical calculations of the heat transfer coefficient and temperature distribution along the horizontal single pin fin in still ambient air, as well as the fin efficiency, are presented and compared with classical analytical results in the case of the constant heat transfer coefficient fin theory. The measured temperature distributions along the two low carbon steel pin fins having a length-to-diameter ratio of 35—one covered with the polished nickel and the other painted mat black—agree very well with the numerical results and are higher than the classical results. The analytically calculated fin efficiency does not differ significantly from the results of the numerical calculations if they are compared for the same dimensionless fin parameter in which the heat transfer coefficient is determined for the fin base temperature. More extended numerical calculations showed that beyond the fin parameter of five, the analytical results of the fin efficiency are higher than the numerical results by no more than about 1%. The largest difference between the classical and numerical fin base efficiencies is about 3.5%, and it was observed at a fin parameter of about 1, where the length of the pin fin has the optimal value based on the classical theory.     

Application of homotopy perturbation method and inverse prediction of thermal parameters for an annular fin subjected to thermal load

In this work, application of the homotopy perturbation method (HPM) and an inverse solution for estimating unknown thermal parameters such as the variable thermal conductivity parameter (β), the thermogeometric parameter (K), and the nondimensional coefficient of thermal expansion (χ) in an annular fin subjected to thermal stresses is presented. Initially, to obtain the nondimensional temperature distribution from the heat equation, the forward method is employed using an approximate analytical solution based on HPM. Thereafter, a closed form solution for the temperature-dependent thermal stresses is obtained using the classical theory of thermoelasticity coupled with HPM solution containing the temperature distribution. Next, for satisfying a particular stress criterion which makes relevance in selecting appropriate configurations for selecting the finned system, unknown thermal parameters are obtained using an inverse approach based on the Nelder–Mead simplex search minimization technique. The objective function is taken as the sum of square of the residuals between the measured stress field and an initially guessed value which is updated iteratively. It is found that more than one type of temperature distribution may yield a given stress distribution, thereby giving rise to different fin efficiencies. The agreement between the actual and the predicted results was found to be satisfactory.     

Convective heat transfer and second law analysis of non-Newtonian fluid flows with variable thermophysical properties in circular channels

This study presents an analytical solution, for fully developed non-Newtonian fluid flows in circular channels under isoflux thermal boundary conditions based on perturbation techniques. Since the physical properties are generally a function of temperature and may not be assumed constant under certain circumstances, the change in viscosity and thermal conductivity with temperature was taken into account. Viscous dissipation term was also included in the performed analysis. In this study, first closed form expressions for velocity, temperature distributions, and Nusselt numbers corresponding to constant thermophysical properties were given in terms of governing parameters. Then, numerical calculation was performed to obtain the values of Nusselt number and global entropy generation for variable thermophysical properties. The results revealed that neglecting the property variation significantly affects heat transfer characteristics and entropy generation, in which the deviation from the constant physical property assumption may reach up to about 32.6%.     

Series solution of convective radiative conduction equation of the nonlinear fin with temperature dependent thermal conductivity

In this study, differential transform method (DTM) is used to evaluate the analytical solution of the nonlinear fin problem with temperature dependent thermal conductivity. Results are presented for the dimensionless temperature distribution and fin efficiency for a range of values of the problem parameters. Using DTM, the differential equation and related boundary conditions are transformed into a recurrence set of equations and finally, the coefficients of power series are obtained based on the solution of this set of equations. Results of this method are compared with the variational iteration method (VIM) and numerical solutions. Here, it is shown that the results of DTM are more accurate than the VIM with lesser calculation cost, which is the main innovation of the current study.     

Thermal Analysis of Orthotropic Pin Fins With Contact Resistance: A Closed-Form Analytical Solution

Light weight composite fins are considered to deal with thermal management problems for many microelectronic components. These composite fins are inherently anisotropic, therefore cannot be handled by a traditional one-dimensional approach; however, these materials can be designed to provide high thermal conductivity values in the desired direction to handle application-specific demands. In this article, we present analytical solutions for temperature distribution and heat transfer rate for orthotropic two-dimensional pin fins subject to convective-tip boundary condition and the contact resistance at the fin base. The generalized results are presented in terms of fin aspect ratio (fin length-to-radius ratio) and three dimensionless fin parameters that relate the internal conductive resistance to three convective resistances discussed in terms of dimensionless variables such as contact, tip, and axial Biot numbers, in addition to the axial-to-radial conductivity ratio. Several special cases including the insulated tip boundary condition are presented. It is demonstrated that the temperature distribution and heat transfer rate from the two-dimensional isotropic annular fin introduced earlier in the literature, can easily be recovered from the benchmark solutions presented in this article. Furthermore, dimensionless heat transfer rates are presented for the pin fins with contact resistance that can help to solve design and optimization problems of many natural-to-forced convection composite fins that are typically encountered in many microelectronic applications.     

Application of Simplex Search Method for Predicting Unknown Parameters in an Annular Fin Subjected to Thermal Stresses

This article presents an inverse approach aimed at estimating the unknown parameters such as the coefficient of thermal expansion and the Biot number for satisfying a prescribed thermal stresses field in radial fin geometry. The variation of temperature with the radius of fin is obtained by solving the heat conduction-convection equation using regular perturbation method and applying proper boundary conditions. A closed form solution for the temperature dependent stress field has been derived by employing the classical elasticity theory coupled with semi-analytical solution for the temperature field. Using the data obtained from a forward method based on the analytical solution, two unknown parameters such as the coefficient of thermal expansion and the Biot number of the fin are simultaneously estimated by an inverse technique using the Nelder-Mead simplex search method. Effects of measurement errors and number of measurement points have been analyzed in detail. It is found that even with 10 measurement points a fairly good reconstruction of the stress field can be achieved.     

Minimization of thermal resistance in an air cooled porous matrix made up of solid spheres with heat generation

In this paper an analytical model was developed to minimize the thermal resistance of an air cooled porous matrix made up of solid spheres with internal heat generation. This was done under the assumption of local thermal equilibrium. The analytical solution of the optimum sphere diameter was found to be independent of the heat generation rate of the solid spheres, but was dependent on the applied pressure drop and fluid properties. The analytical model compared very well to a numerical model found in a computational fluid dynamics code when air and liquid water properties were used for the fluid phase and wood and silica/sand properties were used for the solid phase.