Abstract

In this study, the optimum design dimensions and performance analyses of convective-radiative cooling fin subjected to magnetic field are presented using finite element method. The numerical solutions are verified by the exact analytical solution for the linearized models using Laplace transform. The optimum dimensions for the optimum performance of the convection-radiative fin with variable thermal conductivity are investigated and presented graphically. Also, the effects of convective, radiative, and magnetic parameters as well as Biot number on the thermal performance of the cooling fin are analyzed using the numerical solutions. From the results, it is established that the optimum length of the fin and the thermogeometric parameter increases as the nonlinear thermal conductivity term increases. Further analyses also reveal that as the Biot number, convective, radiative, and magnetic parameters, increases, the rate of heat transfer from the fin increases and consequently improves the efficiency of the fin. Additionally, effects of the thermal stability values for the various multiboiling heat transfer modes are established. It is established that, in order to ensure stability and avoid numerical diffusion of the solution by the Galerkin finite element method, the thermogeometric parameter must not exceed some certain values for the different multiboiling heat transfer modes. It is hope that the present study will enhance the understanding of thermal response of solid fin under various factors and fin design considerations.

1. Introduction

The continuous productions of high performance thermal equipment due to the growing demands for the thermal systems require the development of enhanced heat transfer devices for the effective performance and thermal management of the equipment. Moreover, the excessive heat leads to thermal-induced failure in the thermal systems which necessitate the use of fins for heat transfer enhancement. The applications of extended surfaces in thermal systems, electronic and microelectronics components, high-power semiconductor devices, high-power lasers, light emitting diodes (LEDs), computer cooling, sensitive devices, etc. have attracted various research interests in past decades. The thermal analysis of the extended surface involves the development of thermal models for various operating conditions. Different analytical (exact and approximate) and numerical methods have been employed by various researchers to analyze the thermal performance of the extended surfaces under different conditions. Although, different exact analytical methods have been employed [28], they are based on the assumptions of constant thermal properties. Indubitably, the idealization of a uniform or constant heat transfer coefficient is not realistic. This is because in practice, heat transfer coefficients have significantly greater values at the fin tip than more than the fin base. Also, the thermal property is temperature-dependent. Such variation of the heat transfer coefficient as a function of temperature is often governed by a power law. Moreover, the thermal conductivity of the fin is also temperature-dependent. Under these circumstances, the differential equations governing the thermal response of the fin become strictly nonlinear. The development of exact analytical solutions for such nonlinear models is very difficult. Consequently, some of the past studies have developed approximate analytical solution in terms of series solutions for the thermal analyses of fins using different approximate analytical methods [928]. Nevertheless, the series solutions involve large number of terms. In practice, such huge expressions are not convenient for use by designers and engineers. Therefore, over the years, various numerical methods have been applied for thermal analysis of various extended surfaces [2944].

It should be noted that most of the above studies are based on steady state analysis of fin. The transient response of fins is important in a wide range of engineering devices, automobiles, and industrial sectors. In fact, an accurate transient analysis provides insight into the design of fins that would fail in steady state operations but are sufficient for desired operating periods. Consequently, transient analysis of the extended surfaces has been investigated [4566]. However, they are established on the assumption of fins with insulated tips. Although the effects of tip conditions have been analyzed in the past studies [6772], they are based on steady stated heat transfer.

The determinations of the optimum dimension of fins have been one of the central focus in the study of effective performance and applications of the extended surfaces. Consequently, some past studies have focused on determining the optimum dimensions and performance of fins operating under different conditions. On such optimization studies, Hrymak [73] presented a combined analysis and optimization of extended heat transfer surfaces while Sonn and Bar-Cohen [74] analyzed optimum cylindrical of pin fin. Razelos and Imre [75] carried out a study on minimum mass convective fins with variable heat transfer coefficients. The authors [76] also investigated the optimum dimensions of circular fins with variable thermal parameters. Pitchumani and Shenoy [77] submitted a unified approach to determining optimum shapes for cooling fins of various geometries. Bhargava and Duffin [78] applied the nonlinear method of Wilkins for cooling fin optimization. Wilkins [79] determined the minimum mass of thin fins which transfer heat only by radiation to surroundings at absolute zero. Razelos [80] investigated the optimum dimensions of convective pin fins with internal heat generation. Maday [81] explored the minimum weight of one-dimensional straight fin. Laor and Kalman [82] studied the performance and optimum dimensions of different cooling fins with a temperature-dependent heat transfer coefficient. Effects of magnetic field and radiative heat transfer on the extended surfaces under steady state conditions have considered in literature [1, 8386]. Ma et al. [87] adopted spectral collocation method for transient thermal analysis of coupled conductive, convective, and radiative heat transfer in the moving plate with temperature-dependent properties and heat generation. In another work, with the aid of spectral element method, Ma et al. [88] studied the combined conductive, convective, and radiative heat transfer in moving irregular porous fins. Chen et al. [89] utilized least square spectral collocation method to examine the nonlinear heat transfer in moving porous plate with convective and radiative boundary conditions.

The past studies have been based on the optimum analysis of fin under steady state. To the best of the authors’ knowledge, the transient analysis and optimum design analysis of heat transfers in convective-radiative cooling fin with convective tip under the influence of magnetic field using finite element method have not been studied in open literature. Therefore, in this present study, finite element method is used to determine the optimum dimension and study the transient thermal behaviour of convective-radiative fin with convective tip and under the influence of magnetic field is analyzed. The numerical solutions are used to investigate of thermogeometric parameter and nonlinear thermal conductivity on the optimum performance of the fin. Also, effects of convective, radiative, magnetic, and convective tip parameters on the transient thermal performance of the cooling fin are investigated. Additionally, effects the thermal stability values for the various multiboiling heat transfer modes are established.

2. Problem Formulation

Consider a straight fin of length L and thickness which is exposed on both faces to a convective-radiative environment at temperature and subjected to a uniform magnetic field as shown in Figure 1. In order to develop the mathematical model governing the thermal behaviour, the following assumptions are made:(i)The fin material is homogeneous and isotropic and with constant physical properties.(ii)The thermal properties of the surrounding medium and the magnetic field vary with temperature according to the power law.(iii)The temperature of the surrounding fluid and the temperature of the base of the fin are uniform.(iv)The temperature variation inside the fin is one-dimensional. This is because the fin thickness is so small compared to its height that temperature gradients normal to the surface may be neglected.(v)Heat loss through the fin edges is negligible compared to that which passes through the sides.(vi)There is no contact resistance where the base of the fin joins the prime surface. Also, the specific heat capacity of the fin is constant.(vii)There are no heat sources or internal heat generation within the fin.

From the energy balance analysis based on the assumptions stated above, the thermal energy balance could be expressed as where Jc is the conduction current intensity and is given asAfter expansion and simplification of the above equation, using the fact that as , then (1) reduces According to Fourier’s law of heat conduction, the rate of heat conduction in the fin is given by while the radiation heat transfer rate is given asFrom (4) and (5), it could be stated that the total rate of heat transfer is given byIf one substitutes (6) into (3), one arrives atA further simplification of (7) gives the governing differential equation for the fin as The initial and boundary conditions areHowever, if the tip of the fin is assumed insulated or a negligible rate of heat transfer from it, we have It should be noted thatAfter substitution of (11) into (8), If a small temperature difference exists within the fin material during the heat flow, the term T4 can be expressed as a linear function of temperature asAlso, with the aid of Rosseland’s approximation, one can express the nonlinear derivative term in (12) asIf one substitutes (13) and (14) into (12), we arrive atwhereFor most industrial applications the heat transfer coefficient may be given as the power law [3, 20], where the exponents and are constants. The constant may vary between −6.6 and 5. However, in most practical applications it lies between −3 and 3 [20]. So, the power temperature-dependent thermal properties of the surrounding fluid and the magnetic field are defined as Extending the same power temperature-dependent relationship to the fin emissivity and the magnetic field, we haveThe exponent on the heat transfer coefficient represents laminar film boiling or condensation when p = −1/4, laminar natural convection when p = 1/4, turbulent natural convection when p = 1/3, nucleate boiling when p = 2, and radiation when p = 3. p = 0 implies a constant heat transfer coefficient.

Substitution of (17)-(19) gives

3. Finite Element Method for the Transient Analysis

It is very difficult to develop exact analytical solution to the nonlinear equation in (15) or (20). Therefore, Galerkin finite element method is used in this work to solve the nonlinear equation. Using the shape/interpolating function on the governing equation and Integrating over the domain V of the control volume according to Galerkin finite element method, we haveFor the one-dimensional problem of which the dependent variable varies only along x-axis and the boundary integrals turn to be a point value on boundaries, one can replace dV by Acrdx in (21). Here, is the uniform cross-sectional area of the fin and P is the perimeter of the fin from which convection takes place. Equation (23) is a weak formulation of the nonlinear governing differential equation.

Following the usual finite element procedures, the analysis of (23) givesEquation (24) can be written form aswhereandis the shape function matrix andis the vector of unknown temperaturesUsing a 3-node element as shown in Figure 2, one arrives atAfter the substitution of (32)-(34), we arrived at

3.1. Time Discretization Using the Finite Element Method

The above equation is a general representation of a one-dimensional problem with one linear element. All the terms are included irrespective of the boundary condition. Equation (25) is semidiscrete as it is discretized only in space. The differential operator still contains the time-dependent term and it has to be discretized. We now require a method of discretizing the transient terms of the equation. The following subsections give the details of how the transient terms will be discretized. In (25), the temperature is now discretized in the time domain as in Figure 2.

Using Figure 3,Following the similar procedure as done previously, we can derive the linear shape functions asTherefore, the time derivative of the temperature is thus written asSubstituting (36) and (37) into (25) and applying the weighted residual principle (Galerkin method), we obtain for a time interval of ,After expansionSubstituting (37) into (40)Again, after expansion of (41), one arrives atAfter the evaluation of (40), we obtained the characteristic equation over the time interval Δt asThe above equation involves the temperature values at the nth and (n + 1)th level. The boundary conditions and the temperature-dependent parameters are incorporated in the computer program used to solve the system of differential equations. Although a dimensionless form of the governing equation can be derived for the computer program, the handling of physical quantities is simplified. It should be noted that the thermal properties are evaluated directly in each time step from the nodal temperatures. This eliminates any iteration within each time step for the evaluations of the temperature-dependent parameters.

The element equation/matrix has been derived as shown in the previous equations. It should be noted that the whole domain was divided into a set of 50 line elements. Assembling all the elements equation/matrices, a global matrix or a system of equations was obtained. After applying the boundary conditions, the resulting systems of equations are solved numerically. The convergence criterion of the numerical solution along with error estimation has been set to where is the general dependent variable and is the number of iteration.

It should be noted that a steady state is attained when or .

4. Optimization of the Longitudinal Fin

The fin weight and material costs are the primary design considerations in most applications of fin. Therefore, it is highly desirable to obtain the optimum design information of fins.

The optimization of the fin could be achieved either by minimizing the volume (weight) for any required heat dissipation or by maximizing the heat dissipation for any given fin volume [7376]. The later approach is adopted in this work.

The constant fin volume is defined as V=w. We can therefore write the heat dissipation per unit volume as The dimensionless form of (45) is given asEquation (46) could be written as where The maximum heat dissipation value occurs at the condition when the optimum fin characteristics have been achieved. The fin dimensions in this situation represent the optimum fin configuration per unit volume. With the volume constant, the optimization procedure is also carried out to fix the profile area by first expressing as a function of the thermogeometric parameter, M (or fin length, b) and then searching for the optimum value of M [7376] where .

5. Development of an Exact Analytical Solution for the Verification of the Numerical Solution

For constant thermal properties of the surrounding fluid and the magnetic field, we have a linear equation of the form It should be noted that (49) can be solved analytically. Using Laplace transform, it can easily be shown that the exact analytical solution of the equation based on the initial and the boundary conditions in (9) is given asFor the insulated tip, we havewhere are the positive roots of the characteristics equationIt should be noted that a steady state is attained when .

And the steady state solution for the fin with convective tip can be written asWhile, for the fin with insulated tip, For the sake of convenience in subsequent analysis, it should be noted that “b” has been replaced with “L” in the above analytical solution.

Table 1 shows the thermogeometric parameters used for the simulation.

6. Results and Discussion

For the computational domain, numerical solutions are computed and the necessary convergence of the results is achieved with the desired degree of accuracy. Using the numerical solutions, parametric studies are carried out. Also, in order to define the validity of the results of thermal analysis of fin with assumed insulated tip and that of convective tip, effects of the fin tip conditions on the transient thermal response are investigated. The results with the discussion are illustrated through the Figures 416 and Table 2 to substantiate the applicability of the present analysis.

In order to verify the accuracy of the present numerical method, the numerical results are compared with results obtained by exact analytical method for the linearized equation (Table 2). It is inferred from the figure that there are excellent agreements between the FEM results and the analytical results, which testifies to the validity of the FEM code. This validation boosts the confidence in the numerical outcome of the present study. Moreover, it is observed that, in the same domain by increasing the polynomial degree of approximation or the number of nodes in an element, one can achieve the desired accuracy with less DOF.

Table 2 shows the comparison of the results obtained by exact analytical and finite element methods for a conductive-convective fin with constant thermal and physical properties of fin having negligible radiation and magnetic field effects. Very good agreements are found between the exact analytical and finite element solutions. The average percentage error of the numerical solution is 0.133%.

Figure 4 shows the nondimensional heat transfer (for a unit fin volume) varying with M from 1 and 2 for specified values of nonlinear thermal conductivity terms, β, under a given profile area, ; the heat transfer first rises and then falls as the fin length increases. From the figure, it shows that the optimum fin length (at which reaches a maximum value) increases as the nonlinear thermal conductivity term, β, increases. It also shows that the optimum value of M can be obtained based upon the value of nonlinear term. Therefore, from the analysis the optimum dimensions of the convection fin with variable thermal conductivity are established and the relative values of optimum M and β are shown in Figure 5.

Figure 6 depicts the effects of multiboiling parameter on the optimum heat transfer rate. For the case of carrying out the optimization of the fin by minimizing the volume of the for the required heat dissipation, Figure 7 shows the effects of multiboiling parameter on the optimum fin dimensions. It is illustrated from the figures that the optimum heat transfer rate and fin length increase with increasing multiboiling heat transfer mode parameter while the optimum fin width decreases with the multiboiling heat transfer mode parameter.

Figure 8 shows the comparison of result of analytical solution and finite element method. Figures 9 and 10 depict temperature-time history at different four points (0.025 m, 0.050m, 0.075 m, and 0.100 m) of the convective-radiative fin with convective and insulated tips, respectively. Figure 11 shows the temperature profile of different heat transfer modes while Figures 12 and 13 show the temperature profiles of the fin at difference times. The temperature histories at the four points decrease at a faster rate initially, slow down thereafter, and finally tend to reach a constant value showing to be near to steady state. Also, a marginal or slightly higher temperature differences are notice between the convective and insulated tip. However, this temperature differences become appreciable as the length of the fin increases and heat is transferred within a short period of time. It could be inferred from the figures and the preceding discussion that, for a short fin that undergo heat transfer for a prolonged period of time, adiabatic/insulted condition at the tip can be assumed without any significant loss in accuracy.

It has been established that the criterion and errors due to one-dimensional heat transfer analysis is that fin base thickness Biot number should be much smaller than unity (precisely, Bi < 0.1). To this end, a one-dimensional analysis has been carried out and simulated within 0 < Bi <0.1. In this case the error made in the determination of the rate of heat transfer from the fin to the fluid surrounding it is less than 1% [68, 69]. However, when the Biot number is greater than 0.1 (Bi > 0.1), two-dimensional analysis of the fin is recommended as one-dimensional analysis predicts unreliable results for such limit.

Figures 14 and 15 show the effects of Biot number (conduction-convection parameter) on the temperature distribution in the fin with convective and insulated tips, respectively. From the figures, it is shown that, as the Biot number increases, the rate of heat transfer through the fin increases as the temperature in the fin drops faster (becomes steeper reflecting high base heat flow rates) as depicted in the figures.

Effects of heat transfer coefficient on the temperature distribution in the fin are shown in Figure 16. It is shown that the temperature profiles for the various heat transfer coefficient coincide initially but part away as we move towards the tip of the fin. This is due to the fact that coefficient of heat transfer coefficient is a factor/multiplier of the temperature difference between the fin surface and surrounding medium (T-T). It should be noted that the temperature difference between the fin surface and the surrounding decreases as we move away from the fin base to the fin tip despite the increase in the heat transfer coefficient.

It should be noted that, for the fin with heat transfer coefficient which varies according to power law, the hypothetical boundary condition (that is, insulation) at the tip of the fin is taken into account. If the tip is not assumed to be insulated, then the problem becomes overdetermined [69]. This boundary condition is realized for sufficiently long fins. Also, it should be stated that the assumption that the heat transfer coefficient is constant yields incorrect.

Figure 17 presents the impact of emissivity on the temperature distribution. The temperature of the fin decreases with increase of emissivity value. This is because of increase of emissive heat by radiation from the fin surface especially when the distance from the base increases. Therefore, heat transfer rate increases as the emissivity increases. The radiative eat transfer can be neglected if the base temperature of the fin is low and the emissivity of the fin surface is near zero. The important things in fins surface must be emissive because high emissivity gives a great amount of heat radiation transfer from the fin [39]. Figure 18 shows effects of magnetic parameter and Hartman number on the temperature distribution in the porous fin. The figure depicts that the induced magnetic field in the fin can improve heat transfer through the fin. It is shown that increase in magnetic field on the fin increases the rate of heat transfer from the fin and consequently improves the efficiency of the fin. Figure 19 shows the effect of thermal conductivity of the fin materials on the thermal response of the fin. It could be inferred from the figure that more heat is transferred from fin made of copper material than the fins made of stainless steel and aluminum materials.

Effects of the thermal and geometric parameters on the temperature profile of the fin are shown in Figure 20 while Figure 21 shows the influence of thermogeometric parameter (M=(hP/kA)0.5) on the thermal stability of the fin. It was established that the value of M produces physically unsound behaviour for larger values of the thermogeometric parameter. It is shown that for growing values of the thermogeometric parameter the temperature tends to negative values at the tip of the fin which shows thermal instability, contradicting the assumption of (9). Following the assumptions made regarding the numerical solution of the problem, it was realized that these solutions are not only physically unsound but also point towards thermal instability. Therefore, in order for the solution to be physically sound the fin thermogeometric parameter must not exceed a specific value. By extension in order to ensure stability and avoid numerical diffusion of the solution by the Galerkin finite element method, the thermogeometric parameter must not exceed some certain values for the different multiboiling heat transfer modes.

7. Conclusion

The determination of optimum dimensions and performance of convective-radiative cooling fins subjected to magnetic field has been carried out in this study using finite element method. The numerical solutions were verified by the exact solution for the linearized models using Laplace transform. The optimum dimensions of the convection-radiative fin with variable thermal conductivity were investigated and presented graphically. Also, the effects of other operating parameters on the thermal performance of the fin were investigated using the numerical solutions. The study established that increase in Biot number, convective, radiative, and magnetic parameters, increases the rate of heat transfer from the fin and consequently improves the efficiency of the fin. It is hope that the present study will enhance the understanding of thermal response of solid fin under various factors and fin design considerations.

Nomenclature

:Cross-sectional area of the fins,
:Surface area of the fins,
:Length of the fin
:Magnetic field intensity (T)
:Specific heat (J kg−1 K−1)
:Heat transfer coefficient (Jm−2 K−1)
:Heat transfer coefficient at the base of the fin (Wm−2k−1)
:Total current intensity (A)
:Conduction current intensity (A)
k:Thermal conductivity of the fin material (Wm−1k−1)
:Thermal conductivity of the fin material at the base (Wm−1k−1)
K:Permeability (m−1)
L:Length of the fin (m)
M:Dimensionless thermogeometric parameter
P:Perimeter of the fin (m)
:Heat transfer rate, W
Qf:Dimensionless optimum fin parameter
:Thickness of the fin (m)
t:Time
T:Fin temperature (K)
:Ambient temperature, K
:Temperature at the base of the fin, K
u:Fluid velocity (m/s)
V:Volume of the fin (m3)
w:Width of the fin (m)
x:Axial length measured from fin base (m).
Greek Symbols
ε:Emissivity
:Electric conductivity (A/m)
σst:Stefan–Boltzmann constant (W K4)
ρ:Density of the fluid (kgm−3)
β:Thermal conductivity parameter
δ:Thickness of the fin, m
ρ:Density of the fin material (kg/m3).

Data Availability

The data used in this study are in open literature.

Disclosure

The author declares that the present research did not receive specific funding but was performed as part of the employment of the author in the University of Lagos, Nigeria.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The author expresses sincere appreciation to University of Lagos, Nigeria, for providing material supports and good environment for this work.