Research Article  Open Access
M. G. Sobamowo, G. A. Oguntala, A. A. Yinusa, "Nonlinear Transient Thermal Modeling and Analysis of a ConvectiveRadiative Fin with Functionally Graded Material in a Magnetic Environment", Modelling and Simulation in Engineering, vol. 2019, Article ID 7878564, 16 pages, 2019. https://doi.org/10.1155/2019/7878564
Nonlinear Transient Thermal Modeling and Analysis of a ConvectiveRadiative Fin with Functionally Graded Material in a Magnetic Environment
Abstract
Nonlinear transient thermal analysis of a convectiveradiative 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 inhomogeneity index improve the thermal performance of the fin. Also, the transient responses reveal that the FGM fin with quadraticlaw and exponentiallaw 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.
1. Introduction
The increasing demands for highperformance 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), highpower semiconductor devices, highpower lasers, lightemitting 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 AlNimr [1] and Kiwan [2–4] on the thermal performance analysis of porous fin in the natural convection environment. ShalchiTabrizi 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 crosssectional 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 crosscut 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 doublelayered 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 platepin 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 fanshaped 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 hightemperature 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 indepth 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 nonFourier 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 convectiveradiative 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 length b and thickness t exposed on both faces to a convectiveradiative environment at temperature as shown in Figure 1, assuming the fin medium is homogeneous, isotropic, and saturated with a singlephase fluid. Moreover, the physical properties of solid as well as fluid are considered as constant, and the temperature variation inside the fin is onedimensional, 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.
Following the assumptions, the governing differential equation is developed as
But the Lorentz force can be written as
After substitution of equation (2) into equation (1), taking the magnetic term as a linear function of temperature, we have
Initial condition is stated aswhile the boundary conditions are
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:
Quadratic variation of thermal conductivity with the fin length:
Exponential variation of thermal conductivity with the fin length:
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:
Case 2. Quadratic variation of thermal conductivity with the fin length:
Case 3. Exponential variation of thermal conductivity with the fin length:
In order to nondimensionalize the developed governing equations, the following dimensionless parameters are used in equations (9)–(11):
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:
Case 2. Quadratic variation of thermal conductivity with the fin length:
Case 3. Exponential variation of thermal conductivity with the fin length:and the dimensionless initial and boundary conditions are
Case 1. Linear variation of thermal conductivity with the fin length:
Case 2. Quadratic variation of thermal conductivity with the fin length:
Case 3. Exponential variation of thermal conductivity with the fin length:
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
Using the recursive relation in equation (21), one arrives at
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
With the aid of the recursive relation in equation (23), we have
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
From the recursive relation in equation (25), one arrives at
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 temperatureinvariant physical and thermal properties of the fin and expression T^{4} as a linear function of temperature can be done without any loss of generality and accuracy:
On substituting equation (27) into equation (9), we arrived at the following equation:
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:
Case 2. Quadratic variation of thermal conductivity with the fin length:
Case 3. Exponential variation of thermal conductivity with the fin length:
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:
Case 2. Quadratic variation of thermal conductivity with the fin length:
Case 3. Exponential variation of thermal conductivity with the fin length:and the dimensionless initial and boundary conditions are
4.1. Method of Solution Using Laplace Transform Method
In order to provide closedform 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 function and its inversion formulas are defined as follows:where 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:
Case 2. Quadratic variation of thermal conductivity with the fin length:
Case 3. Exponential variation of thermal conductivity with the fin length:and the dimensionless initial and boundary conditions in Laplace domain are
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:
Case 2. Quadratic variation of thermal conductivity with the fin length:where
Case 3. Exponential variation of thermal conductivity with the fin length:
It should be noted that the first and second kinds of Bessel, modified Bessel, and Legendre functions are given as follows:
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
Lee et al. [50] suggested values of 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]
Case 1. Linear variation of thermal conductivity with the fin length:
Case 2. Quadratic variation of thermal conductivity with the fin length:where
Case 3. Exponential variation of thermal conductivity with the fin length:
For homogenous material (HM), β = 0. The solution for the homogeneous material using Laplace transform is given as
5. Results and Discussion
The results of the exact analytical solutions of the fin with the functionally graded material are shown in Figures 2–8. The influences of magnetic field, thermal radiation, inhomogeneity index, and the temperature distribution and history are investigated.
(a)
(b)
(c)
(d)
(a)
(b)
Figures 2–4 show the effects of varying time on the temperature distribution on the FGM fin with variable thermal conductivity according to linear, quadratic, and exponentiallaw 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 quadraticlaw and exponentiallaw 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 exponentiallaw function requires small time to reach steady state than the linear and quadraticlaw functions.
Figures 5(a)–5(d) show the variation of temperature with length at difference dimensionless time for linear, quadratic, and exponentiallaw 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 exponentiallaw 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.
Figures 6 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, Figures 8(a) and 8(b) show the effects of the inhomogeneity index on the dimensionless temperature distribution and consequently on the rate of heat transfer. From the figures, it is shown that as the inhomogeneity 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 exponentiallaw functions. Moreover, increasing the inhomogeneity 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 inhomogeneity index, β, of the FGM fin. It is therefore established from these results that the temperature profiles of the fin are highly sensitive in linearlaw function with FGM than that of powerlaw 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 powerlaw class FGM heat sink is significantly higher than the linear class FGM heat sink.
Table 1 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.

6. Conclusion
Transient thermal analysis of a convectiveradiative 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 inhomogeneity 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 inhomogeneity index improve the thermal performance of the fin. Also, it was established that the transient responses of the functionally graded material (FGM) fin with quadraticlaw and exponentiallaw 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.
Nomenclature
a_{r}:  Aspect ratio of the porous fin base area to the surface area 
A:  Crosssectional area of the fins (m^{2}) 
A_{b}:  Porous fin base area 
A_{s}:  Porous fin surface area 
Bi:  Biot number 
h:  Heat transfer coefficient (W·m^{−2}·k^{−1}) 
h_{b}:  Heat transfer coefficient at the base of the fin (W·m^{−2}·k^{−1}) 
c_{p}:  Specific heat of the fluid passing through the porous fin (J/kg·K) 
Da:  Darcy number 
:  Gravity constant (m/s^{2}) 
h:  Heat transfer coefficient over the fin surface (W/m^{2}·K) 
H:  Dimensionless heat transfer coefficient at the base of the fin (W·m^{−2}·k^{−1}) 
k:  Thermal conductivity of the fin material (W·m^{−1}·k^{−1}) 
k_{b}:  Thermal conductivity of the fin material at the base of the fin (W·m^{−1}·k^{−1}) 
k_{eff}:  Effective thermal conductivity ratio 
K:  Permeability of the porous fin (m^{2}) 
L:  Length of the fin (m) 
M:  Dimensionless thermogeometric parameter 
m:  Mass flow rate of fluid passing through the porous fin (kg/s) 
P:  Perimeter of the fin (m) 
Q:  Dimensionless heat transfer rate per unit area 
q_{b}:  Heat transfer rate per unit area at the base (W/m^{2}) 
Q_{b}:  Dimensionless heat transfer rate at the base in the porous fin 
Q_{s}:  Dimensionless heat transfer rate at the base in the solid fin 
Ra:  Rayleigh number 
S_{h}:  Porosity parameter 
t:  Thickness of the fin 
T_{b}:  Base temperature (K) 
T:  Fin temperature (K) 
T_{a}:  Ambient temperature (K) 
T_{b}:  Temperature at the base of the fin (K) 
:  Average velocity of fluid passing through the porous fin (m/s) 
x:  Axial length measured from the fin tip (m) 
X:  Dimensionless length of the fin 
:  Width of the fin. 
β:  Inhomogeity index 
δ:  Thickness of the fin (m) 
δ_{b}:  Fin thickness at its base 
γ:  Dimensionless internal heat generation parameter 
θ:  Dimensionless temperature 
θ_{b}:  Dimensionless temperature at the base of the fin 
η:  Efficiency of the fin 
β′:  Coefficient of thermal expansion (K^{−1}) 
ɛ:  Porosity or void ratio 
υ:  Kinematic viscosity (m^{2}/s) 
ρ:  Density of the fluid (kg/m^{3}). 
s:  Solid properties 
f:  Fluid properties 
eff:  Effective porous properties. 
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors express their sincere appreciation to University of Lagos, Nigeria, for providing material supports and good environment for this work. This research was performed as part of the employment of the authors under the University of Lagos, Nigeria.
References
 S. Kiwan and M. A. AlNimr, “Using porous fins for heat transfer enhancement,” Journal of Heat Transfer, vol. 123, no. 4, pp. 790–795, 2000. View at: Publisher Site  Google Scholar
 S. Kiwan, “Effect of radiative losses on the heat transfer from porous fins,” International Journal of Thermal Sciences, vol. 46, no. 10, pp. 1046–1055, 2007. View at: Publisher Site  Google Scholar
 S. Kiwan, “Thermal analysis of natural convection porous fins,” Transport in Porous Media, vol. 67, no. 1, pp. 17–29, 2006. View at: Publisher Site  Google Scholar
 S. Kiwan and O. Zeitoun, “Natural convection in a horizontal cylindrical annulus using porous fins,” International Journal of Numerical Methods for Heat & Fluid Flow, vol. 18, no. 5, pp. 618–634, 2008. View at: Publisher Site  Google Scholar
 A. ShalchiTabrizi and H. R. Seyf, “Analysis of entropy generation and convective heat transfer of Al_{2}O_{3} nanofluid flow in a tangential micro heat sink,” International Journal of Heat and Mass Transfer, vol. 55, no. 1516, pp. 4366–4375, 2012. View at: Publisher Site  Google Scholar
 S. M. H. Hashemi, S. A. Fazeli, H. Zirakzadeh, and M. Ashjaee, “Study of heat transfer enhancement in a nanofluidcooled miniature heat sink,” International Communications in Heat and Mass Transfer, vol. 39, no. 6, pp. 877–884, 2012. View at: Publisher Site  Google Scholar
 T.C. Hung, W.M. Yan, X.D. Wang, and C.Y. Chang, “Heat transfer enhancement in microchannel heat sinks using nanofluids,” International Journal of Heat and Mass Transfer, vol. 55, no. 910, pp. 2559–2570, 2012. View at: Publisher Site  Google Scholar
 H. R. Seyf and M. Feizbakhshi, “Computational analysis of nanofluid effects on convective heat transfer enhancement of micropinfin heat sinks,” International Journal of Thermal Sciences, vol. 58, pp. 168–179, 2012. View at: Publisher Site  Google Scholar
 S. A. Fazeli, S. M. H. Hashemi, H. Zirakzadeh, and M. Ashjaee, “Experimental and numerical investigation of heat transfer in a miniature heat sink utilizing silica nanofluid,” Superlattices and Microstructures, vol. 51, no. 2, pp. 247–264, 2012. View at: Publisher Site  Google Scholar
 S.M. Kim and I. Mudawar, “Analytical heat diffusion models for different microchannel heat sink crosssectional geometries,” International Journal of Heat and Mass Transfer, vol. 53, no. 1920, pp. 4002–4016, 2010. View at: Publisher Site  Google Scholar
 P. Naphon, S. Klangchart, and S. Wongwises, “Numerical investigation on the heat transfer and flow in the minifin heat sink for CPU,” International Communications in Heat and Mass Transfer, vol. 36, no. 8, pp. 834–840, 2009. View at: Publisher Site  Google Scholar
 P. Naphon and O. Khonseur, “Study on the convective heat transfer and pressure drop in the microchannel heat sink,” International Communications in Heat and Mass Transfer, vol. 36, no. 1, pp. 39–44, 2009. View at: Publisher Site  Google Scholar
 T. Y. Kim and S. J. Kim, “Fluid flow and heat transfer characteristics of crosscut heat sinks,” International Journal of Heat and Mass Transfer, vol. 52, no. 2324, pp. 5358–5370, 2009. View at: Publisher Site  Google Scholar
 T.C. Hung, W.M. Yan, and W.P. Li, “Analysis of heat transfer characteristics of doublelayered microchannel heat sink,” International Journal of Heat and Mass Transfer, vol. 55, no. 1112, pp. 3090–3099, 2012. View at: Publisher Site  Google Scholar
 Z. M. Wan, G. Q. Guo, K. L. Su, Z. K. Tu, and W. Liu, “Experimental analysis of flow and heat transfer in a miniature porous heat sink for high heat flux application,” International Journal of Heat and Mass Transfer, vol. 55, no. 1516, pp. 4437–4441, 2012. View at: Publisher Site  Google Scholar
 D. Lelea, “Effects of inlet geometry on heat transfer and fluid flow of tangential microheat sink,” International Journal of Heat and Mass Transfer, vol. 53, no. 1718, pp. 3562–3569, 2010. View at: Publisher Site  Google Scholar
 X. Yu, J. Feng, Q. Feng, and Q. Wang, “Development of a platepin fin heat sink and its performance comparisons with a plate fin heat sink,” Applied Thermal Engineering, vol. 25, no. 23, pp. 173–182, 2005. View at: Publisher Site  Google Scholar
 L. Chai, G. Xia, M. Zhou, and J. Li, “Numerical simulation of fluid flow and heat transfer in a microchannel heat sink with offset fanshaped reentrant cavities in sidewall,” International Communications in Heat and Mass Transfer, vol. 38, no. 5, pp. 577–584, 2011. View at: Publisher Site  Google Scholar
 Wikipedia, “Functionally graded material,” https://en.wikipedia.org/wiki/functionally_graded_material. View at: Google Scholar
 R. Hassanzadeh and M. Bilgili, “Improvement of thermal efficiency in computer heat sink using functionally graded materials,” Communications on Advanced Computational Science with Applications, vol. 2014, Article ID cacsa00018, 13 pages, 2014. View at: Publisher Site  Google Scholar
 C. Bilgili, “Thermal management technologies for electronic packaging: current capabilities and future challenges for modelling tools,” in Proceedings of 2008 10th Electronics Packaging Technology Conference, pp. 527–532, Singapore, December 2008. View at: Google Scholar
 D. P. Kulkarni and D. K. Das, “Analytical and numerical studies on microscale heat sinks for electronic applications,” Applied Thermal Engineering, vol. 25, no. 1415, pp. 2432–2449, 2005. View at: Publisher Site  Google Scholar
 T. Dang, J.T. Teng, and J.C. Chu, “A study on the simulation and experiment of a microchannel counterflow heat exchanger,” Applied Thermal Engineering, vol. 30, no. 1415, pp. 2163–2172, 2010. View at: Publisher Site  Google Scholar
 L. Li, M. Nagar, and X. Jie, “Effect of thermal interface materials on manufacturing and reliability of Flip Chip PBGA and SiP packages,” in Proceedings of 2008 58th Electronic Components and Technology Conference, pp. 973–978, Lake Buena Vista, FL, USA, May 2008. View at: Google Scholar
 X. Luo, W. Xiong, T. Cheng, and S. Liu, “Design and optimization of horizontallylocated plate fin heat sink for high power LED street lamps,” in Proceedings of 2009 59th Electronic Components and Technology Conference, pp. 854–859, San Diego, CA, USA, May 2009. View at: Google Scholar
 R. Kandasamy, X.Q. Wang, and A. S. Mujumdar, “Transient cooling of electronics using phase change material (PCM)based heat sinks,” Applied Thermal Engineering, vol. 28, no. 89, pp. 1047–1057, 2008. View at: Publisher Site  Google Scholar
 K.T. Chiang, “Modeling and optimization of designing parameters for a parallelplain fin heat sink with confined impinging jet using the response surface methodology,” Applied Thermal Engineering, vol. 27, no. 1415, pp. 2473–2482, 2007. View at: Publisher Site  Google Scholar
 D. Gottlieb and S. A. Orszag, Numerical Analysis of Spectral Methods: Theory and Applications (CBMSNSF Regional Conference Series in Applied Mathematics), Society for Industrial and Applied Mathematics, Philadelphia, PA, USA.
 Y. H. Claudio Canuto, A. Quarteroni, and T. A. Zang, Spectral Methods in Fluid Dynamics, Springer Berlin Heidelberg, Berlin, Germany, 1988.
 R. Peyret, Spectral Methods for Incompressible Viscous Flow, SpringerVerlag, New York, NY, USA, 2002.
 F. B. Belgacem and M. Grundmann, “Approximation of the wave and electromagnetic diffusion equations by spectral method,” SIAM Journal on Scientific Computing, vol. 20, pp. 13–32, 1998. View at: Google Scholar
 X. Shan, D. Montgomery, and H. Chen, “Nonlinear magnetohydrodynamics by Galerkinmethod computation,” Physical Review A, vol. 44, no. 10, pp. 6800–6818, 1991. View at: Publisher Site  Google Scholar
 J. P. Wang, “Fundamental problems in spectral methods and finite spectral method,” Acta Aerodynamica Sinica, vol. 19, pp. 161–171, 2001. View at: Google Scholar
 E. M. E. Elbarbary and M. ElKady, “Chebyshev finite difference approximation for the boundary value problems,” Applied Mathematics and Computation, vol. 139, no. 23, pp. 513–523, 2003. View at: Publisher Site  Google Scholar
 Z. J. Huang and Z. J. Zhu, “Chebyshev spectral collocation method for solution of Burgers’ equation and laminar natural convection in twodimensional cavities,” University of Science and Technology of China, Hefei, China, 2009, B.Sc. thesis. View at: Google Scholar
 N. T. Eldabe and M. E. M. Ouaf, “Chebyshev finite difference method for heat and mass transfer in a hydromagnetic flow of a micropolar fluid past a stretching surface with Ohmic heating and viscous dissipation,” Applied Mathematics and Computation, vol. 177, no. 2, pp. 561–571, 2006. View at: Publisher Site  Google Scholar
 A. H. Khater, R. S. Temsah, and M. M. Hassan, “A Chebyshev spectral collocation method for solving Burgers’type equations,” Journal of Computational and Applied Mathematics, vol. 222, no. 2, pp. 333–350, 2008. View at: Publisher Site  Google Scholar
 E. H. Doha, A. H. Bhrawy, and S. S. EzzEldien, “Efficient Chebyshev spectral methods for solving multiterm fractional orders differential equations,” Applied Mathematical Modelling, vol. 35, no. 12, pp. 5662–5672, 2011. View at: Publisher Site  Google Scholar
 B. Kundu, R. Das, P. A. Wankhade, and K.S. Lee, “Heat transfer improvement of a wet fin under transient response with a unique design arrangement aspect,” International Journal of Heat and Mass Transfer, vol. 127, pp. 1239–1251, 2018. View at: Publisher Site  Google Scholar
 P. A. Wankhade, B. Kundu, and R. Das, “Establishment of nonFourier heat conduction model for an accurate transient thermal response in wet fins,” International Journal of Heat and Mass Transfer, vol. 126, pp. 911–923, 2018. View at: Publisher Site  Google Scholar
 R. Das and B. Kundu, “Direct and inverse approaches for analysis and optimization of fins under sensible and latent heat load,” International Journal of Heat and Mass Transfer, vol. 124, pp. 331–343, 2018. View at: Publisher Site  Google Scholar
 R. Das and B. Kundu, “Estimation of internal heat generation in a fin involving all modes of heat transfer using golden section search method,” Heat Transfer Engineering, vol. 39, no. 1, pp. 58–71, 2017. View at: Publisher Site  Google Scholar
 R. Das and B. Kundu, “Prediction of heat generation in a porous fin from surface temperature,” Journal of Thermophysics and Heat Transfer, vol. 31, no. 4, pp. 781–790, 2017. View at: Publisher Site  Google Scholar
 S. S. EzzEldien, “New quadrature approach based on operational matrix for solving a class of fractional variational problems,” Journal of Computational Physics, vol. 317, pp. 362–381, 2016. View at: Publisher Site  Google Scholar
 A. H. Bhrawy and S. S. EzzEldien, “A new Legendre operational technique for delay fractional optimal control problems,” Calcolo, vol. 53, no. 4, pp. 521–543, 2015. View at: Publisher Site  Google Scholar
 S. S. EzzEldien and E. H. Doha, “Fast and precise spectral method for solving pantograph type Volterra integrodifferential equations,” Numerical Algorithms, pp. 1–21, 2018, In press. View at: Publisher Site  Google Scholar
 A. H. Bhrawy, S. S. EzzEldien, E. H. Doha, M. A. Abdelkawy, and D. Baleanu, “Solving fractional optimal control problems within a ChebyshevLegendre operational technique,” International Journal of Control, vol. 90, no. 6, pp. 1230–1244, 2017. View at: Publisher Site  Google Scholar
 M. G. Sobamowo, “Nonlinear thermal and flowinduced vibration analysis of fluidconveying carbon Nanotube resting on Winkler and Pasternak foundations,” Thermal Science and Engineering Progress, vol. 4, pp. 133–149, 2017. View at: Publisher Site  Google Scholar
 M. G. Sobamowo, O. Adeleye, and A. A. Yinusa, “Analysis of convectiveradiative porous fin with temperaturedependent internal heat Generation and magnetic field using Homotopy Perturbation method,” Journal of Computational and Applied Mechanics, vol. 12, no. 2, pp. 127–145, 2017. View at: Google Scholar
 S. T. Lee, M. C. H. Chien, and W. E. Culham, “Vertical singlewell pulse testing of a threelayer stratified reservoir,” in Proceedings of SPE Annual Technical Conference and Exhibition, Society of Petroleum Engineers, Houston, TX, USA, 1984. View at: Google Scholar
Copyright
Copyright © 2019 M. G. Sobamowo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.