# Nonlinear Transient Thermal Modeling and Analysis of a Convective-Radiative Fin with Functionally Graded Material in a Magnetic Environment

Authors: M. G. Sobamowo; G. A. Oguntala; A. A. Yinusa
Journal: Modelling and Simulation in Engineering
 (2019)

---

## Abstract

Nonlinear transient thermal analysis of a convective-radiative fin with functionally graded materials (FGMs) under the influence of magnetic field is presented in this study. The developed nonlinear thermal models of linear, quadratic, and exponential variation of thermal conductivity are solved approximately and analytically using the differential transformation method (DTM). In order to verify the accuracies of the nonlinear solutions, exact analytical solutions are also developed with the aids of Bessel, Legendre, and modified Bessel functions. Good agreements are established between the exact and the approximate analytical solutions. In the parametric studies, effects of heat enhancement capacity of fin with functionally graded material as compared to fin with homogeneous material are investigated. Also, influence of the Lorentz force and radiative heat transfer on the thermal performance of the fin are analyzed. From the results, it is shown that increase in radiative and magnetic field parameters as well as the in-homogeneity index improve the thermal performance of the fin. Also, the transient responses reveal that the FGM fin with quadratic-law and exponential-law function shows the slowest and fasted thermal responses, respectively. This study will provide a very good platform for the design and optimization of an improved heat transfer enhancement in thermal systems, where the surrounding fluid is influenced by a magnetic field.


---

## Body

## 1. Introduction

The increasing demands for high-performance engineering systems come with inherent increased heat generation in the thermal and electronics systems which consequently lead to increased thermal damages. Effective thermal management and dissipation of excess heat from heat transfer components have been a major concern and challenge in the design of the thermal and electronic systems. In the quest for the thermal management of these systems, the heat sink has proved to be an effective, passive means of reducing the thermally induced failures by enhancing heat dissipation from thermal systems and electronics systems. Heat sinks are heat exchangers used in dissipating heat from functional thermal systems to the environment to ensure whether the device operates within safe temperature limits. The wide range of applications of heat sinks in cooling different electronic and micro electronics components such as the central processing unit (CPU), high-power semiconductor devices, high-power lasers, light-emitting diodes (LEDs), and sensitive devices affirm its effectiveness as a passive mode of cooling of thermal systems. The improvement of the thermal performance of heat sink and consequently energy saving in the thermal system is timely owing to miniaturization in size of thermal systems, acoustic control due to its destructive effects, and energy saving using smaller fan size.To improve the thermal performance of cooling devices using porous fin as a passive method for heat transfer enhancement in thermal equipment, different authors have proposed different innovative analysis and approach, succeeding the pioneer of work of Kiwan and Al-Nimr [1] and Kiwan [2–4] on the thermal performance analysis of porous fin in the natural convection environment. Shalchi-Tabrizi and Seyf [5] employed nanofluid and presented the entropy generation and convective heat transfer in a micro heat sink. Hashemi et al. [6] investigated the heat transfer enhancement in a miniature heat sink using nanofluid, while Hung et al. [7] presented a study on heat transfer enhancement in micro channel heat sinks using nanofluids. Seyf and Feizbakhshi [8] carried out a computational analysis of nanofluid effects on convective heat transfer enhancement of micro pin fin heat sinks, while Fazeli et al. [9] investigated experimentally and numerically the heat transfer in a miniature heat sink utilizing silica nanofluid. However, other studies focus on the properties of conventional airflow and features of channel cross section as an approach to improve heat transfer enhancement in heat sinks. In these studies, Kim and Mudawar [10] developed analytical heat diffusion models for different micro channel heat sink cross-sectional geometries. Naphon et al. [11] presented a numerical investigation of the heat transfer flow in the minifin heat sink for CPU. Naphon and Khonseur [12] studied the convective heat transfer and pressure drop in the micro channel heat sink whilst Kim and Kim [13] investigated fluid flow and the characteristics of heat transfer in cross-cut heat sinks. Furthermore, the effects of various geometries and airflow paths on the heat sink performance have also been investigated in the previous literature. Hung et al. [14] analyzed the heat transfer characteristics of the double-layered micro channel heat sink. Wan et al. [15] presented an experimental analysis of flow and heat transfer in a miniature porous heat sink for high heat flux applications. Lelea [16] studied the effects of inlet geometry on heat transfer and fluid flow of tangential micro heat sink. Yu et al. [17] presented an enhanced heat transfer study using a plate-pin fin heat sink. Chai et al. [18] carried out a numerical simulation of fluid flow and heat transfer in a micro channel heat sink with offset fan-shaped reentrant cavities in sidewall. The previous studies explored the convective parameters around the heat sink while some past research works employed the properties of the conventional airflow, features of channel cross section, variable geometries, and airflow paths to enhanced the heat sinks thermal performance However, the use of materials with changing composition, microstructure, or porosity across the material volume has been proven to be an effective way of improving the thermal performance of heat sink. Such materials with varying electrical, chemical, mechanical, magnetic, and thermal properties over the volume of the bulk material are referred to as functionally graded materials (FGMs). The continuous variations in properties in the FGM along a specific axis could be based on the porosity and pore size gradient structure, chemical gradient structure, and microstructural gradient structure of the FGM.FGM was initially developed and used for high-temperature applications in Japan during a space plane project [19]. However, their effective thermal, physical, chemical, and magnetic properties have paved ways for their wide range of applications such as engine components, heat shielding of satellites, thermal stress control of cylindrical and spherical bodies, and vessels and several applications in aerospace, automobile, and nuclear industries. Although there are numerous studies in literature on cooling of electronics components using heat sink [20–27], an in-depth examination and review of the past works indicate to the best of our knowledge that FGM has not been applied in heat sinks except in [20]. Recent works on the applications of the non-Fourier heat transfer in extended surfaces, transient analysis of fins under dehumidification, and other optimization studies as well as the recent advancements in Legendre polynomials have been well presented [28–47]. Thus, considering the significant capabilities and benefits of FGM coupled with established enhanced heat transfer characteristics of porous fins. However, different studies on the thermal performance of fins with and without functionally graded materials have been carried out. However, to the best of the authors’ knowledge, a study on the transient thermal analysis of fin with functionally graded materials is scarce. Therefore, in this work, nonlinear transient thermal analysis of a convective-radiative fin with functionally graded materials under the influence of magnetic field is presented. The developed nonlinear thermal models of linear, quadratic, and exponential variation of thermal conductivity are solved approximately and analytically using the differential transformation method. This method converges very fast and is very efficient for the handling of both ordinary and partial differential equations. However, when the domain of convergence is at infinity, there is always a need to augment the classical differential transform method with PADE approximants, cosine and sine after treatments, domain transformation, or multistep DTM. In order to verify the accuracies of the nonlinear solutions, exact analytical solutions are also developed with the aids of Bessel, Legendre, and modified Bessel functions. The effects of heat enhancement capacity of fin with the functionally graded material as compared to fin with the homogeneous material are investigated. Also, influence of the Lorentz force and radiative heat transfer on the thermal performance of the fin are analyzed.


## 2. Problem Formulation

Consider a heat sink made of FGM having lengthb and thickness t exposed on both faces to a convective-radiative environment at temperature T∞ as shown in Figure 1, assuming the fin medium is homogeneous, isotropic, and saturated with a single-phase fluid. Moreover, the physical properties of solid as well as fluid are considered as constant, and the temperature variation inside the fin is one-dimensional, i.e., temperature varies along the length only and remain constant with time and there is no thermal contact resistance at the fin base and the fin tip is an adiabatic type.Figure 1
Schematic diagram of heat transfer in a longitudinal fin with FGM.Following the assumptions, the governing differential equation is developed as(1)ρcp∂T∂t=∂∂xkx∂T∂x−hPAT−Ta−σεPAT4−Ta4−Jc×JcσA.But the Lorentz force can be written as(2)Jc×Jcσ=σB02u2.After substitution of equation (2) into equation (1), taking the magnetic term as a linear function of temperature, we have(3)ρcp∂T∂t=∂∂xkx∂T∂x−hPAT−Ta−σεPAT4−Ta4−σB02u2AT−Ta.Initial condition is stated as(4)t=0,T=Ta,x>0,while the boundary conditions are(5)t>0,x=0,∂T∂x=0,t>0,x=L,T=Tb.In most practical situations, the fin thickness at the end is small that the heat transfer from the fin tip can be neglected and solution can be obtained assuming fin is insulated at the tip. Therefore, in this work, we assumed that the tip of the fin is adiabatic/insulated.In the present analysis, cases of linear, quadratic, and exponential variation of thermal conductivity with the fin length will be considered.Linear variation of thermal conductivity with the fin length:(6)kx=k01+γlx.Quadratic variation of thermal conductivity with the fin length:(7)kx=k01+γqx2.Exponential variation of thermal conductivity with the fin length:(8)kx=k0eγix.It should be noted that when rapid cooling is required in the domain of the FGM fin (such as in the condenser, electronic devices, and CPU), FGM materials with quadratic thermal conductivity variation may be applied. However, if slow heat enhancement is required as in the fluidized bed, exponential is advised. For an intermediate and moderate heat transfer or enhancement process, linear thermal conductivity variation will be useful.Case 1. Linear variation of thermal conductivity with the fin length:(9)ρcpk0∂T∂t=∂2T∂x2+γlx∂2T∂x2+γl∂T∂x−hPk0AT−Ta−σεPAT4−Ta4−σB02u2k0AT−Ta.Case 2. Quadratic variation of thermal conductivity with the fin length:(10)ρcpk0∂T∂t=∂2T∂x2+γqx2∂2T∂x2+2γqx∂T∂x−hPk0AT−Ta−σεPAT4−Ta4−σB02u2k0AT−Ta.Case 3. Exponential variation of thermal conductivity with the fin length:(11)ρcpk0∂T∂t=eγix∂2T∂x2+γieγix∂T∂x−hPk0AT−Ta−σεPAT4−Ta4−σB02u2k0AT−Ta.In order to nondimensionalize the developed governing equations, the following dimensionless parameters are used in equations (9)–(11):(12)X=xL,θ=T−TaTb−Ta,τ=k0tρcpL2,β=γlL=γiL,λ=βqL2,Nc2=PLhAk0,Nr=4σεPTa3k0A,Ma2=σB02u2k0A,CT=T∞Tb−T∞.Then, the dimensionless forms of the governing equations (9)–(11) can be written as follows.Case 1. Linear variation of thermal conductivity with the fin length:(13)∂θ∂τ=∂2θ∂X2+βX∂2θ∂X2+β∂θ∂X−Nc2θ−Nrθ+CT4−CT4−Ma2θ=0.Case 2. Quadratic variation of thermal conductivity with the fin length:(14)∂θ∂τ=∂2θ∂X2+λX2∂2θ∂X2+2λX∂θ∂X−Nc2θ−Nrθ+CT4−CT4−Ma2θ=0.Case 3. Exponential variation of thermal conductivity with the fin length:(15)∂θ∂τ=eβX∂2θ∂X2+βeβX∂θ∂X−Nc2θ−Nrθ+CT4−CT4−Ma2θ=0,and the dimensionless initial and boundary conditions are(16)τ=0,θ=0,X>0,τ>0,X=0,∂θ∂X=0,(17)τ>0,X=1,θ=1.Case 1. Linear variation of thermal conductivity with the fin length:(18)∂θ∂τ=∂2θ∂X2+βX∂2θ∂X2+β∂θ∂X−Nc2θ−Nrθ4−4NrCTθ3−6NrCT2θ2−4NrCT3θ−Maθ=0.Case 2. Quadratic variation of thermal conductivity with the fin length:(19)∂θ∂τ=∂2θ∂X2+λX2∂2θ∂X2+2λX∂θ∂X−Nc2θ−Nrθ4−4NrCTθ3−6NrCT2θ2−4NrCT3θ−Maθ=0.Case 3. Exponential variation of thermal conductivity with the fin length:(20)∂θ∂τ=eβX∂2θ∂X2+βeβX∂θ∂X−Nc2θ−Nrθ4−4NrCTθ3−6NrCT2θ2−4NrCT3θ−Maθ=0.


## 3. Method of Solution for the Nonlinear Thermal Models: Differential Transform Method

The developed nonlinear thermal models in equations (18)–(20) are very difficult to solve exactly and analytically. Therefore, a recourse was made to an approximation analytical method, differential transform method. The basic definitions and the operational properties of the method can be found in our previous study [48].Case 1. Linear variation of thermal conductivity with the fin length.The recursive relation for the linear variation of thermal conductivity with the fin length given in equation (20) is (21)h+1θk,h+1=k+1k+2θk+2,h+β∑r=0k∑u=0hδr−1,h−uk−r+1k−r+2θk−r+2,u+βk+1θk+1,h−Nc2+Ma2θk,h−Nr∑r=0k∑s=0k−r∑t=0k−r−s∑u=0h∑v=0h−u∑w=0h−u−vθr,h−u−v−wθs,uθt,vθk−r−s−t,w+4CT∑r=0k∑s=0k−r∑u=0h∑v=0h−uθr,h−u−vθs,uθk−r−s,w+6CT2∑r=0k∑u=0hθr,h−uθk−r,w+4CT3θk,h.Using the recursive relation in equation (21), one arrives at(22)θ2,1=σ+Ma2σ2+Nc2σ2+2NrσCT3,θ3,1=−16βσ2+Ma2+Nc2+4NrCT3,θ4,1=σ4+Ma2σ6+Ma4σ24+Nc2σ6+112Ma2Nc2σ+Nc4σ24+β2σ12+124Ma2β2σ+124Nc2β2σ+12Nrσ2CT2+23NrσCT3+13Ma2NrσCT3+13Nc2NrσCT3+16Nrβ2σCT3+23Nr2σCT6,θ2,2=2NrCT3σ+3CT2Nrσ2+12Ma2σ+12Nc2σ+32,θ3,2=−16βσ3+Nc2+Ma2+6NrσCT2+4NrCT3,θ5,1=−βσ10−115Nc2βσ−160Nc4βσ−β3σ60−1120Nc2β3σ−115βσMa2−130Nc2βσMa2−1120β3σMa2−160βσMa4−15Nrβσ2CT2−415NrβσCT3−215Nc2NrβσCT3−130Nrβ3σCT3−215NrβσMa2CT3−415Nr2βσCT6.Case 2. Quadratic variation of thermal conductivity with the fin length.The recursive relation for the quadratic variation of thermal conductivity with the fin length given in equation (19) is(23)h+1θk,h+1=k+1k+2θk+2,h+λ∑r=0k∑u=0hδr−2,h−uk−r+1k−r+2θk−r+2,u+2λ∑r=0k∑u=0hδr−1,h−uk−r+1θk−r+1,u−Nc2+Ma2θk,h−Nr∑r=0k∑s=0k−r∑t=0k−r−s∑u=0h∑v=0h−u∑w=0h−u−vθr,h−u−v−wθs,uθt,vθk−r−s−t,w+4CT∑r=0k∑s=0k−r∑u=0h∑v=0h−uθr,h−u−vθs,uθk−r−s,w+6CT2∑r=0k∑u=0hθr,h−uθk−r,w+4CT3θk,h.With the aid of the recursive relation in equation (23), we have(24)θ2,1=σ+Ma2σ2+Nc2σ2+2NrσCT3,θ3,1=0,θ2,2=3σ2+Ma2σ2+Nc2σ2+3Nrσ2CT2+2NrσCT3,θ3,2=0,θ4,1=σ4+Ma2σ6+Ma4σ24+Nc2σ6+112Ma2Nc2σ+Nc4σ24+12Nrσ2CT2+23NrσCT3+13Ma2NrσCT3+13Nc2NrσCT3+23Nr2σCT6,θ5,1=0.Case 3. Exponential variation of thermal conductivity with the fin length.The recursive relation for the exponential variation of thermal conductivity with the fin length given in equation (20) is(25)h+1θk,h+1=∑r=0k∑u=0hβrr!δh−uk−r+1k−r+2θk−r+2,u+β∑r=0k∑u=0hβrr!δh−uk−r+1θk−r+1,u−Nc2+Ma2θk,h−Nr∑r=0k∑s=0k−r∑t=0k−r−s∑u=0h∑v=0h−u∑w=0h−u−vθr,h−u−v−wθs,uθt,vθk−r−s−t,w+4CT∑r=0k∑s=0k−r∑u=0h∑v=0h−uθr,h−u−vθs,uθk−r−s,w+6CT2∑r=0k∑u=0hθr,h−uθk−r,w+4CT3θk,h.From the recursive relation in equation (25), one arrives at(26)θ2,1=σ+Ma2σ2+Nc2σ2+2NrσCT3,θ3,1=−2βσ3−13Ma2βσ−13Nc2βσ−43NrβσCT3,θ2,2=3σ2+Ma2σ2+Nc2σ2+3Nrσ2CT2+2NrσCT3,θ3,2=−βσ−13Ma2βσ−13Nc2βσ−2Nrβσ2CT2−43NrβσCT3,θ4,1=σ4+Ma2σ6+Ma4σ24+Nc2σ6+112Ma2Nc2σ+Nc4σ24+β2σ4+18Ma2β2σ+18Nc2β2σ+12Nrσ2CT2+23NrσCT3+13Ma2NrσCT3+13Nc2NrσCT3+12Nrβ2σCT3+23Nr2σCT6,θ5,1=−3βσ10−15Ma2βσ−120Ma4βσ−15Nc2βσ−110Ma2Nc2βσ−120Nc4βσ−β3σ15−130Ma2β3σ−130Nc2β3σ−35Nrβσ2CT2−45NrβσCT3−25Ma2NrβσCT3−25Nc2NrβσCT3−215Nrβ3σCT3−45Nr2βσCT6.


## 4. Development of Exact Analytical Solution for the Linearized Thermal Models

In order to verify the accuracies of the nonlinear solutions as developed in the previous sections, exact analytical solutions are also developed. In order to do this, we considered a case where the temperature difference within the material during the heat flow is small. Under such scenario, the adoption of temperature-invariant physical and thermal properties of the fin and expressionT4 as a linear function of temperature can be done without any loss of generality and accuracy:(27)T4≅4Ta3T−3Ta4.On substituting equation (27) into equation (9), we arrived at the following equation:(28)ρcp∂T∂t=∂∂xkx∂T∂x−hPAT−Ta−4σεPTa3AT−Ta−σB02u2AT−Ta.Upon substitution of equations (6)–(8) into equation (28) and expanding the resulting equations, we have the following.Case 1. Linear variation of thermal conductivity with the fin length:(29)ρcpk0∂T∂t=∂2T∂x2+γlx∂2T∂x2+γl∂T∂x−hPk0AT−Ta−4σεPTa3k0AT−Ta−σB02u2k0AT−Ta.Case 2. Quadratic variation of thermal conductivity with the fin length:(30)ρcpk0∂T∂t=∂2T∂x2+γqx2∂2T∂x2+2γqx∂T∂x−hPk0AT−Ta−4σεPTa3k0AT−Ta−σB02u2k0AT−Ta.Case 3. Exponential variation of thermal conductivity with the fin length:(31)ρcpk0∂T∂t=eγix∂2T∂x2+γieγix∂T∂x−hPk0AT−Ta−4σεPTa3k0AT−Ta−σB02u2k0AT−Ta.Applying the nondimensionalize parameters in equation (12) to equations (29)–(31), we obtain the dimensionless forms of the linearized governing equations as given below.Case 1. Linear variation of thermal conductivity with the fin length:(32)∂θ∂τ=∂2θ∂X2+βX∂2θ∂X2+β∂θ∂X−Nc2θ−Nrθ−Ma2θ.Case 2. Quadratic variation of thermal conductivity with the fin length:(33)∂θ∂τ=∂2θ∂X2+λX2∂2θ∂X2+2λX∂θ∂X−Nc2θ−Nrθ−Ma2θ.Case 3. Exponential variation of thermal conductivity with the fin length:(34)∂θ∂τ=eβX∂2θ∂X2+βeβX∂θ∂X−Nc2θ−Nrθ−Ma2θ,and the dimensionless initial and boundary conditions are(35)τ=0,θ=0,X>0,τ>0,X=0,∂θ∂X=0,τ>0,X=1,θ=1.


### 4.1. Method of Solution Using Laplace Transform Method

In order to provide closed-form solutions to the developed equations, we adopted integral transforms using Laplace transform. We first apply Laplace transform over the time and then solve the resulting equation using Bessel, modified Bessel, Legendre, and sign and gamma functions. The procedures are given as follows:The Laplace transform of a real functionft and its inversion formulas are defined as follows:(36)Fs=∫0∞e−stftdt,ft=12πi∫s−i∞s+i∞e−stFsdt,where s=a+iba,b∈R is a complex number.Applying Laplace transform on equations (15)–(17), we have the equation in Laplace domain as follows.Case 1. Linear variation of thermal conductivity with the fin length:(37)d2θ˜dX2+βXd2θ˜dX2+βdθ˜dX−s+Nc2+Nr+Ma2θ˜.Case 2. Quadratic variation of thermal conductivity with the fin length:(38)d2θ˜dX2+λX2d2θ˜dX2+2λXdθ˜dX−s+Nc2+Nr+Ma2θ˜.Case 3. Exponential variation of thermal conductivity with the fin length:(39)eβXd2θ˜dX2+βeβXdθ˜dX−s+Nc2+Nr+Ma2θ˜,and the dimensionless initial and boundary conditions in Laplace domain are(40)s=0,θ˜=0,X>0,(41)s>0,X=0,∂θ˜∂X=0,s>0,X=1,θ=1s.With the aids of Bessel, modified Bessel, Legendre, and sign and gamma functions, it could be shown that the solutions of equations (37)–(39) are as follows.Case 1. Linear variation of thermal conductivity with the fin length:(42)θ˜X,s=J02/β−s+Nc2+Nr+Ma21+βXY12csgn1/β1/β−s+Nc2+Nr+Ma2−Y02/β−s+Nc2+Nr+Ma21+βXJ12csgn1/β1/β−s+Nc2+Nr+Ma2×J02/β−s+Nc2+Nr+Ma2Y12csgn1/β1/β−s+Nc2+Nr+Ma2−Y02/β−s+Nc2+Nr+Ma2J12csgn1/β1/β−s+Nc2+Nr+Ma2−1.Case 2. Quadratic variation of thermal conductivity with the fin length:(43)θ˜X,s=πQn−λX−Γ−n/2Γp/2Qm0Pn−λXπQn−λ−Γ−n/2Γp/2Qm0Pn−λ,where(44)m=λ+4s+Nc2+Nr+Ma2+λ2λ,n=λ+4s+Nc2+Nr+Ma2−λ2λ,p=λ+4s+Nc2+Nr+Ma2+5λ2λ.Case 3. Exponential variation of thermal conductivity with the fin length:(45)θ˜X,s=e−βX−1/2K02/βs+Nc2+Nr+Ma2I12/βs+Nc2+Nr+Ma2e−βX/2+I02/βs+Nc2+Nr+Ma2K12/βs+Nc2+Nr+Ma2e−βX/2×K02/βs+Nc2+Nr+Ma2I12/βs+Nc2+Nr+Ma2e−β/2+I02/βs+Nc2+Nr+Ma2K12/βs+Nc2+Nr+Ma2e−β/2−1.It should be noted that the first and second kinds of Bessel, modified Bessel, and Legendre functions are given as follows:(46)Jνz=∑r=0∞−1rz/22r+νr!Γν+r+1,Yνz=2π∑r=0∞−1rz/22r+νr!Γν+r+1lnz2−1π∑r=0∞ν−r−1!z/22r−νr!−1π∑r=0∞−1rz/22r+nr!r+n!ψr+ν+1+ψr+1,Iνz=∑r=0∞z/22r+νr!Γν+r+1,Kνz=−1ν+1∑r=0∞z/22r+νr!Γν+r+1lnz2+12∑r=0∞ν−r−1!z/22r−νr!+12−1v∑r=0∞−1rz/22r+nr!r+n!ψr+ν+1+ψr+1,Pνz=12νν!ddxνx−1ν=2ν∑r=0νxrνrν+r−12r,Qνz=12Pνzln1+z1−z=2ν−1∑r=0νxrνrν+r−12rln1+z1−z,ψr+1=−γ+∑r=1n1m,csgnz=1,ifRez>0−1,ifRez<0sgnImz,ifRez=0,sgnz=zz.It should be noted thatψ is a digamma function while γ is the Euler–Mascheroni constant and c sgn is a sign function for real and complex expressions.The determination of the inverse Laplace transforms of equations (41)–(43) and (45) is very complex. However, they can be numerically evaluated using Simon’s approach [49] given as(47)θX,τ=eapττ12θ˜X,ap+∑n=1NReθ˜X,ap+inπτ−1n.Lee et al. [50] suggested values of apτ ranging between 4 and 5. Equations (41)–(43) and (45) converge more quickly because it does not contain oscillating cosine and sine functions.The optimally value is [50](48)apτ=4.7⇒ap=4.7τ.Case 1. Linear variation of thermal conductivity with the fin length:(49)θX,τ=eapττ12J02/β−4.7/τ+Nc2+Nr+Ma21+βXY12csgn1/β1/β−4.7/τ+Nc2+Nr+Ma2−Y02/β−4.7/τ+Nc2+Nr+Ma21+βXJ12csgn1/β1/β−4.7/τ+Nc2+Nr+Ma2×J02/β−4.7/τ+Nc2+Nr+Ma2Y12csgn1/β1/β−4.7/τ+Nc2+Nr+Ma2−Y02/β−4.7/τ+Nc2+Nr+Ma2J12csgn1/β1/β−4.7/τ+Nc2+Nr+Ma2−1+∑n=1NReJ02/β−1/τ4.7+inπ+Nc2+Nr+Ma21+βXY12csgn1/β1/β−1/τ4.7+inπ+Nc2+Nr+Ma2−Y02/β−1/τ4.7+inπ+Nc2+Nr+Ma21+βXJ12csgn1/β1/β−1/τ4.7+inπ+Nc2+Nr+Ma2×J02/β−1/τ4.7+inπ+Nc2+Nr+Ma2Y12csgn1/β1/β−1/τ4.7+inπ+Nc2+Nr+Ma2−Y02/β−1/τ4.7+inπ+Nc2+Nr+Ma2J12csgn1/β1/β−1/τ4.7+inπ+Nc2+Nr+Ma2−1−1n.Case 2. Quadratic variation of thermal conductivity with the fin length:(50)θX,τ=eapττ12θ˜X,ap+∑n=1NReθ˜X,ap+inπτ−1n,where(51)θ˜X,ap=πQn′−λX−Γ−n′/2Γp′/2Qm′0Pn′−λXπQn′−λ−Γ−n′/2Γp′/2Qm′0Pn′−λ,m′=λ+44.7/τ+Nc2+Nr+Ma2+λ2λ,n′=λ+44.7/τ+Nc2+Nr+Ma2−λ2λ,p′=λ+44.7/τ+Nc2+Nr+Ma2+5λ2λ.Case 3. Exponential variation of thermal conductivity with the fin length:(52)θX,τ=eapττ12e−βX−1/2K02/β4.7/τ+Nc2+Nr+Ma2I12/β4.7/τ+Nc2+Nr+Ma2e−βX/2+I02/β4.7/τ+Nc2+Nr+Ma2K12/β4.7/τ+Nc2+Nr+Ma2e−βX/2×K02/β4.7/τ+Nc2+Nr+Ma2I12/β4.7/τ+Nc2+Nr+Ma2e−β/2+I02/β4.7/τ+Nc2+Nr+Ma2K12/β4.7/τ+Nc2+Nr+Ma2e−β/2−1+∑n=1NRee−βX−1/2K02/β1/τ4.7+inπ+Nc2+Nr+Ma2I12/β1/τ4.7+inπ+Nc2+Nr+Ma2e−βX/2+I02/β1/τ4.7+inπ+Nc2+Nr+Ma2K12/β1/τ4.7+inπ+Nc2+Nr+Ma2e−βX/2×K02/β1/τ4.7+inπ+Nc2+Nr+Ma2I12/β1/τ4.7+inπ+Nc2+Nr+Ma2e−β/2+I02/β1/τ4.7+inπ+Nc2+Nr+Ma2K12/β1/τ4.7+inπ+Nc2+Nr+Ma2e−β/2−1−1n.For homogenous material (HM),β = 0. The solution for the homogeneous material using Laplace transform is given as(53)θX,τ=coshNc2+Nr+Ma2XcoshNc2+Nr+Ma2+4π∑n=1∞−1n2n−1e−2n−12π2/4+Nc2+Nr+Ma2τcos2n−1πX2.


## 4.1. Method of Solution Using Laplace Transform Method

In order to provide closed-form solutions to the developed equations, we adopted integral transforms using Laplace transform. We first apply Laplace transform over the time and then solve the resulting equation using Bessel, modified Bessel, Legendre, and sign and gamma functions. The procedures are given as follows:The Laplace transform of a real functionft and its inversion formulas are defined as follows:(36)Fs=∫0∞e−stftdt,ft=12πi∫s−i∞s+i∞e−stFsdt,where s=a+iba,b∈R is a complex number.Applying Laplace transform on equations (15)–(17), we have the equation in Laplace domain as follows.Case 1. Linear variation of thermal conductivity with the fin length:(37)d2θ˜dX2+βXd2θ˜dX2+βdθ˜dX−s+Nc2+Nr+Ma2θ˜.Case 2. Quadratic variation of thermal conductivity with the fin length:(38)d2θ˜dX2+λX2d2θ˜dX2+2λXdθ˜dX−s+Nc2+Nr+Ma2θ˜.Case 3. Exponential variation of thermal conductivity with the fin length:(39)eβXd2θ˜dX2+βeβXdθ˜dX−s+Nc2+Nr+Ma2θ˜,and the dimensionless initial and boundary conditions in Laplace domain are(40)s=0,θ˜=0,X>0,(41)s>0,X=0,∂θ˜∂X=0,s>0,X=1,θ=1s.With the aids of Bessel, modified Bessel, Legendre, and sign and gamma functions, it could be shown that the solutions of equations (37)–(39) are as follows.Case 1. Linear variation of thermal conductivity with the fin length:(42)θ˜X,s=J02/β−s+Nc2+Nr+Ma21+βXY12csgn1/β1/β−s+Nc2+Nr+Ma2−Y02/β−s+Nc2+Nr+Ma21+βXJ12csgn1/β1/β−s+Nc2+Nr+Ma2×J02/β−s+Nc2+Nr+Ma2Y12csgn1/β1/β−s+Nc2+Nr+Ma2−Y02/β−s+Nc2+Nr+Ma2J12csgn1/β1/β−s+Nc2+Nr+Ma2−1.Case 2. Quadratic variation of thermal conductivity with the fin length:(43)θ˜X,s=πQn−λX−Γ−n/2Γp/2Qm0Pn−λXπQn−λ−Γ−n/2Γp/2Qm0Pn−λ,where(44)m=λ+4s+Nc2+Nr+Ma2+λ2λ,n=λ+4s+Nc2+Nr+Ma2−λ2λ,p=λ+4s+Nc2+Nr+Ma2+5λ2λ.Case 3. Exponential variation of thermal conductivity with the fin length:(45)θ˜X,s=e−βX−1/2K02/βs+Nc2+Nr+Ma2I12/βs+Nc2+Nr+Ma2e−βX/2+I02/βs+Nc2+Nr+Ma2K12/βs+Nc2+Nr+Ma2e−βX/2×K02/βs+Nc2+Nr+Ma2I12/βs+Nc2+Nr+Ma2e−β/2+I02/βs+Nc2+Nr+Ma2K12/βs+Nc2+Nr+Ma2e−β/2−1.It should be noted that the first and second kinds of Bessel, modified Bessel, and Legendre functions are given as follows:(46)Jνz=∑r=0∞−1rz/22r+νr!Γν+r+1,Yνz=2π∑r=0∞−1rz/22r+νr!Γν+r+1lnz2−1π∑r=0∞ν−r−1!z/22r−νr!−1π∑r=0∞−1rz/22r+nr!r+n!ψr+ν+1+ψr+1,Iνz=∑r=0∞z/22r+νr!Γν+r+1,Kνz=−1ν+1∑r=0∞z/22r+νr!Γν+r+1lnz2+12∑r=0∞ν−r−1!z/22r−νr!+12−1v∑r=0∞−1rz/22r+nr!r+n!ψr+ν+1+ψr+1,Pνz=12νν!ddxνx−1ν=2ν∑r=0νxrνrν+r−12r,Qνz=12Pνzln1+z1−z=2ν−1∑r=0νxrνrν+r−12rln1+z1−z,ψr+1=−γ+∑r=1n1m,csgnz=1,ifRez>0−1,ifRez<0sgnImz,ifRez=0,sgnz=zz.It should be noted thatψ is a digamma function while γ is the Euler–Mascheroni constant and c sgn is a sign function for real and complex expressions.The determination of the inverse Laplace transforms of equations (41)–(43) and (45) is very complex. However, they can be numerically evaluated using Simon’s approach [49] given as(47)θX,τ=eapττ12θ˜X,ap+∑n=1NReθ˜X,ap+inπτ−1n.Lee et al. [50] suggested values of apτ ranging between 4 and 5. Equations (41)–(43) and (45) converge more quickly because it does not contain oscillating cosine and sine functions.The optimally value is [50](48)apτ=4.7⇒ap=4.7τ.Case 1. Linear variation of thermal conductivity with the fin length:(49)θX,τ=eapττ12J02/β−4.7/τ+Nc2+Nr+Ma21+βXY12csgn1/β1/β−4.7/τ+Nc2+Nr+Ma2−Y02/β−4.7/τ+Nc2+Nr+Ma21+βXJ12csgn1/β1/β−4.7/τ+Nc2+Nr+Ma2×J02/β−4.7/τ+Nc2+Nr+Ma2Y12csgn1/β1/β−4.7/τ+Nc2+Nr+Ma2−Y02/β−4.7/τ+Nc2+Nr+Ma2J12csgn1/β1/β−4.7/τ+Nc2+Nr+Ma2−1+∑n=1NReJ02/β−1/τ4.7+inπ+Nc2+Nr+Ma21+βXY12csgn1/β1/β−1/τ4.7+inπ+Nc2+Nr+Ma2−Y02/β−1/τ4.7+inπ+Nc2+Nr+Ma21+βXJ12csgn1/β1/β−1/τ4.7+inπ+Nc2+Nr+Ma2×J02/β−1/τ4.7+inπ+Nc2+Nr+Ma2Y12csgn1/β1/β−1/τ4.7+inπ+Nc2+Nr+Ma2−Y02/β−1/τ4.7+inπ+Nc2+Nr+Ma2J12csgn1/β1/β−1/τ4.7+inπ+Nc2+Nr+Ma2−1−1n.Case 2. Quadratic variation of thermal conductivity with the fin length:(50)θX,τ=eapττ12θ˜X,ap+∑n=1NReθ˜X,ap+inπτ−1n,where(51)θ˜X,ap=πQn′−λX−Γ−n′/2Γp′/2Qm′0Pn′−λXπQn′−λ−Γ−n′/2Γp′/2Qm′0Pn′−λ,m′=λ+44.7/τ+Nc2+Nr+Ma2+λ2λ,n′=λ+44.7/τ+Nc2+Nr+Ma2−λ2λ,p′=λ+44.7/τ+Nc2+Nr+Ma2+5λ2λ.Case 3. Exponential variation of thermal conductivity with the fin length:(52)θX,τ=eapττ12e−βX−1/2K02/β4.7/τ+Nc2+Nr+Ma2I12/β4.7/τ+Nc2+Nr+Ma2e−βX/2+I02/β4.7/τ+Nc2+Nr+Ma2K12/β4.7/τ+Nc2+Nr+Ma2e−βX/2×K02/β4.7/τ+Nc2+Nr+Ma2I12/β4.7/τ+Nc2+Nr+Ma2e−β/2+I02/β4.7/τ+Nc2+Nr+Ma2K12/β4.7/τ+Nc2+Nr+Ma2e−β/2−1+∑n=1NRee−βX−1/2K02/β1/τ4.7+inπ+Nc2+Nr+Ma2I12/β1/τ4.7+inπ+Nc2+Nr+Ma2e−βX/2+I02/β1/τ4.7+inπ+Nc2+Nr+Ma2K12/β1/τ4.7+inπ+Nc2+Nr+Ma2e−βX/2×K02/β1/τ4.7+inπ+Nc2+Nr+Ma2I12/β1/τ4.7+inπ+Nc2+Nr+Ma2e−β/2+I02/β1/τ4.7+inπ+Nc2+Nr+Ma2K12/β1/τ4.7+inπ+Nc2+Nr+Ma2e−β/2−1−1n.For homogenous material (HM),β = 0. The solution for the homogeneous material using Laplace transform is given as(53)θX,τ=coshNc2+Nr+Ma2XcoshNc2+Nr+Ma2+4π∑n=1∞−1n2n−1e−2n−12π2/4+Nc2+Nr+Ma2τcos2n−1πX2.


## 5. Results and Discussion

The results of the exact analytical solutions of the fin with the functionally graded material are shown in Figures2–8. The influences of magnetic field, thermal radiation, in-homogeneity index, and the temperature distribution and history are investigated.Figure 2
Temperature history in the fin at various locations when the thermal conductivity varies according to linear-law function.Figure 3
Temperature history in the fin at various locations when the thermal conductivity varies according to quadratic-law function.Figure 4
Temperature history in the fin at various locations when the thermal conductivity varies according to exponential-law function.Figure 5
Temperature distribution in the fin for different variable thermal conductivity at (a)τ = 0.1, (b) τ = 0.15, (c) τ = 0.30, and (d) τ = 0.35.
(a)
(b)
(c)
(d)Figure 6
Effects of the magnetic field parameter on the temperature distribution.Figure 7
Effects of the radiation number on the temperature distribution.Figure 8
Temperature distribution in the fin for varying thermogeometric parameter under (a) linear-law function and (b) exponential-law functions.
(a)
(b)Figures2–4 show the effects of varying time on the temperature distribution on the FGM fin with variable thermal conductivity according to linear-, quadratic-, and exponential-law functions, respectively. We observe that the temperature increases with an increase in time. The solution seems to converge to a steady state solution as time evolves. We also notice that temperature at the tip of the fin increases with time.Also, it is depicted in the figure, the time required to cool down to the ambient temperature for all the FGM fins of the variable thermal conductivities considered. The FGM fin with quadratic-law and exponential-law functions shows the slowest and fasted thermal responses, respectively. It is illustrated that, in such cooling down process of the fin, the FGM fin with exponential-law function requires small time to reach steady state than the linear- and quadratic-law functions.Figures5(a)–5(d) show the variation of temperature with length at difference dimensionless time for linear-, quadratic-, and exponential-law functions. It is depicted in the figure that the variation of FGM fin thermal conductivity according to the linear law shows an enhanced performance as compared to the quadratic- and exponential-law functions. These enhancements in the energy saving of FGM heat sinks in both linear, quadratic, and exponential laws depict some advantages in the cooling of thermal systems as this causes the lower airflow rate around the heat sinks which results to the smaller fan size and lower electrical energy consumption [16]. The other benefit of smaller fan size is a reduction on the airflow rate within such sensitive devices and controlling the acoustic level, which results to lower distortions, and increasing the electrical or computer component lifetimes.Figures6 and 7 illustrate the impacts of magnetic field (Lorentz force) and radiation number on the temperature distribution in the FGM fin. From the figures, as the magnetic field and radiation parameters increase, the temperature decreases rapidly and the rate of heat transfer through the fin increases as the temperature in the fin drops faster as depicted in the figures. The occurrence of rapid decrease in fin temperature as the magnetic parameter is increased is because the augmentation of the Hartmann number causes the force created by the magnetic field to increase which consequently increases the magnetic field strength and created a reduction in the fin temperature. The illustration shows the thermal efficiency of the fin is enhanced by the magnetic force and thermal radiation is induced in the heat transfer process.As it is displayed, Figures8(a) and 8(b) show the effects of the in-homogeneity index on the dimensionless temperature distribution and consequently on the rate of heat transfer. From the figures, it is shown that as the in-homogeneity index increases, the rate of heat transfer through the fin increases. For all the values of convective, magnetic field, and radiative parameters, the temperature gradient along the fin with FGM is smaller than the fin with homogeneous material for both linear- and exponential-law functions. Moreover, increasing the in-homogeneity index, β, the temperature gradient decreases as shown in the figures. Furthermore, the rate of heat transfer enhancement for FGM longitudinal fin increases by increasing the in-homogeneity index, β, of the FGM fin. It is therefore established from these results that the temperature profiles of the fin are highly sensitive in linear-law function with FGM than that of power-law function. In addition, the application of fin with FGM is favourable at low thermogeometric convective, magnetic field, and radiative parameters since the difference between fin with FGM and fin with HM fin temperature profiles slightly decreases as the dimensionless thermogeometric parameters increase. It could be inferred from the figures that the application of fin with FGM decreases the thermal resistance along the fin, and consequently, a fin with FGM has higher temperature at the fin tip than the fin with HM. The rate of energy saving in power-law class FGM heat sink is significantly higher than the linear class FGM heat sink.Table1 shows the comparison of results between the DTM and the exact analytical method. From the table, it is established that the temperature predictions in the heat sink using the differential transformation method are in excellent agreement with the results of the numerical methods.Table 1
Comparison of results when Nc = 0.01, Ma = 0.1, and Nr = 0.5.
X
Exact
DTM
Error
0.00
0.863499158
0.863499231
0.000000073
0.05
0.863828540
0.863828568
0.000000028
0.10
0.864817031
0.864817090
0.000000059
0.15
0.866465671
0.866465743
0.000000072
0.20
0.868776195
0.868776261
0.000000066
0.25
0.871751037
0.871751104
0.000000067
0.30
0.875393336
0.875393404
0.000000068
0.35
0.879706946
0.879707010
0.000000064
0.40
0.884696438
0.884696500
0.000000062
0.45
0.890367120
0.890367181
0.000000061
0.50
0.896725040
0.896725096
0.000000056
0.55
0.903777007
0.903777060
0.000000053
0.60
0.911530606
0.911530658
0.000000052
0.65
0.919994212
0.919994259
0.000000047
0.70
0.929177015
0.929177056
0.000000041
0.75
0.939089039
0.939089079
0.00000004
0.80
0.949741166
0.949741203
0.000000037
0.85
0.961145166
0.961145189
0.000000023
0.90
0.973313722
0.973313764
0.000000042
0.95
0.986260463
0.986260549
0.000000086
1.00
1.000000000
1.000000000
0.000000000


## 6. Conclusion

Transient thermal analysis of a convective-radiative fin with functionally graded materials under the influence of magnetic field has been presented in this work. Different variations of thermal conductivity according to linear, quadratic, and exponential variations considered in the thermal models are solved analytically with the aids of Bessel, Legendre, and modified Bessel functions, respectively. The impacts of Lorentz force and radiative as well as in-homogeneity index on the thermal performance of the fin have been investigated and discussed. The results show that increase in radiative and magnetic field parameters as well as in-homogeneity index improve the thermal performance of the fin. Also, it was established that the transient responses of the functionally graded material (FGM) fin with quadratic-law and exponential-law function show the slowest and fastest thermal responses, respectively. It is therefore hoped that the study will aid in the design and optimization approach of the fin for improved heat transfer enhancements of thermal systems under the influence of magnetic field.



---

*Source: 7878564-2019-02-14.xml*