Next Article in Journal
Interactions between Artificial Channel Protein, Water Molecules, and Ions Based on Theoretical Approaches
Next Article in Special Issue
Identification of Material Properties of Elastic Plate Using Guided Waves Based on the Matrix Pencil Method and Laser Doppler Vibrometry
Previous Article in Journal
Robust Nonlinear Non-Referenced Inertial Frame Multi-Stage PID Controller for Symmetrical Structured UAV
Previous Article in Special Issue
Fault Diagnosis of Crack on Gearbox Using Vibration-Based Approaches
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploration of Temperature Distribution through a Longitudinal Rectangular Fin with Linear and Exponential Temperature-Dependent Thermal Conductivity Using DTM-Pade Approximant

by
Ravikumar Shashikala Varun Kumar
1,
Rangaswamy Naveen Kumar
1,
Ganeshappa Sowmya
2,
Ballajja Chandrappa Prasannakumara
1 and
Ioannis E. Sarris
3,*
1
Department of Studies in Mathematics, Davangere University, Davangere 577002, India
2
Department of Mathematics, M S Ramaiah Institute of Technology, Bangalore 560054, India
3
Department of Mechanical Engineering, University of West Attica, 12244 Athens, Greece
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(4), 690; https://doi.org/10.3390/sym14040690
Submission received: 2 March 2022 / Revised: 19 March 2022 / Accepted: 23 March 2022 / Published: 26 March 2022
(This article belongs to the Special Issue Solid Mechanics and Mechanical Mechanics)

Abstract

:
The present study elaborates on the thermal distribution and efficiency of a longitudinal rectangular fin with exponentially varying temperature-dependent thermal conductivity and heat transfer coefficient concerning internal heat generation. Also, the thermal distribution of a fin is comparatively studied for both exponentially varying temperature-dependent thermal conductivity and linearly varying temperature-dependent thermal conductivity. Further, the thermal distribution of a longitudinal fin is examined by using ANSYS software with different fin materials. Many physical mechanisms can be explained by ordinary differential equations (ODEs) with symmetrical behavior, the significance of which varies based on the perspective. The governing equation of the considered problem is reduced to a non-linear ODE with the assistance of dimensionless terms. The resultant equation is solved analytically using the DTM-Pade approximant and is also solved numerically using Runge-Kutta Fehlberg’s fourth-fifth (RKF-45) order method. The features of dimensionless parameters influencing the fin efficiency and temperature profile are discussed through graphical representation for exponentially and linearly varying temperature-dependent thermal conductivity. This study ensures that the temperature field enhances for the higher magnitude of thermal conductivity parameter, whereas it diminishes for diverse values of the thermo-geometric parameter. Also, greater values of heat generation and heat transfer parameters enhance the temperature profile. Highlight: Thermal distribution through a rectangular profiled straight fin is examined. Linear and non-linear thermal properties are considered. The combined impact of conduction, convection, and internal heat generation is taken for modeling the energy equation of the fin. Thermal simulation is performed for Aluminum Alloy 6061 (AA 6061) and Cast Iron using ANSYS.

1. Introduction

The heat transfer process is significant when the temperature of different bodies or parts of the same body varies. Convection of fluids, conduction in solids, and radiation are the three modes of heat transfer. A higher heat transfer rate is generally essential in industrial manufacturing processes. Traditional heat transfer methods are ineffective in generating the essential heat transfer rate for industrial applications. The use of cooling fluid, nanofluid, and hybrid nanofluids in the industrial sector is one of the techniques for achieving a high heat transfer rate. Numerous investigations explained the role of these kinds of liquids in the heat transfer process [1,2,3,4,5,6,7]. Usage of extended surfaces in various mechanical and engineering applications seems to be another approach to enhance heat transfer rate. Extended surfaces or fins are frequently used in various engineering systems to expand the heat transference area, dissipate heat generated inside the system, and extend the system’s functioning. Fins are used in oil-carrying pipelines, computer processor cooling, air conditioning, refrigeration, and air-cooled craft engines. Fins also play a prominent role in solidification augmentation in the energy storage system [8,9]. Fins are divided into longitudinal, radial, and pin fins. To improve heat dissipation from a heated main surface, a longitudinal fin of different profiles with a constant cross-sectional area is extensively utilized in practice. This fin is most effective when used in a natural convection setting with a low convection heat transfer coefficient. Motivated by these properties, several researchers examined the heat transfer properties of longitudinal fins with several influencing factors. Kundu [10] elaborated the comparative and unified research on the longitudinal fin and discussed the thermal analysis with optimal design. For a fin having a constant area with temperature variance, internal heat production, and thermal conductivity, Aziz and Bouaziz [11] derived the analytically approximate with very accurate solutions for temperature thermal distribution, fin efficiency, and optimal fin parameter. Najafabadi et al. [12] investigated analytical solutions for the non-linear differential equation (NLDE) model of a longitudinal fin with thermal conductivity and heat production. In the presence of internal heat generation, Sowmya et al. [13,14] conferred the thermal properties of longitudinal fin and moving rod owing to natural convection and radiation.
On the other hand, due to differences in thickness, the longitudinal fin has variable profiles depending on the shape. There are two techniques for solving fin optimization issues. The initial step in solving an optimization issue is to choose a basic profile (triangular or rectangular) and to calculate the fin dimensions. In reality, the rectangular profile is commonly utilized in the heat exchanger system due to its ease of production. As a result, additional research is being done to find the best size for these sorts. Recently, numerous researchers deliberated the heat transference analysis via different fins with rectangular profiles [15,16,17,18,19]. The study of fins for improving heat transmission has piqued the interest of many scholars during the last decade. When a fin’s cross-section is tiny, the temperature variation in the cross-section is correspondingly minimal. Then there is just a considerable temperature variation in the direction of the fin’s length. Put another way, the temperature distribution in the fin is a one-dimensional temperature field running parallel to the fin’s length. Several researchers explored thermal distribution analysis in past decades by considering different fins. Tarobi et al. [20] inspected the aspects of thermal distribution in a radiative-convective fin having various profiles by employing the methodology of DTM and elaborated on the fin’s efficiency. In recent years, many investigators have manifested significant reports to discuss the features of temperature distribution through different kinds of fins [21,22,23,24,25].
The traditional fin analysis assumes constant conductivity and homogeneous convective heat transference coefficient over the fin surface. But in some instances, the fin creates internal heat due to the presence of an electric current, such as in an electric filament, or as a result of a chemical or atomic reaction that can be observed in an atomic reactor. The internal heat generation may be invariable, but it is a function of temperature in a more realistic situation. Two non-linear factors are included in the governing equation for a fin in which both thermal conductivity and internal heat production are temperature dependent. Even with one non-linear factor, the effective analytical solution for this framed differential equation (DE) cannot be obtained easily. As a result, the non-linear fin equation has been solved numerically or semi-analytically. Initially, Minkler and Rouleau [26] inspected a convective fin with homogeneous internal heat production and established mathematical explanations for the thermal distribution. Later, numerous researchers explored the numerical solutions for heat transfer analysis of different fins with internal heat generation [27,28,29,30]. Several researchers have recently derived more feature representations of real-world problems using newly designed simulation studies. Symmetry is recognized as such a method, and researchers obtained the symmetrical features of such a model utilizing mathematical norms. As a result, many investigators presented new operators and models. Moreover, numerous physical processes in fluid dynamics, bio-modeling, and physics are described using symmetrical ordinary differential equations (ODEs). The DTM is an analytic approach used for solving the considered ODEs. In this technique, the solution is expressed in Taylor’s series form. As a unique feature, this procedure can be used effectively to the problem without requiring any linearization, perturbation, or discretization. In addition, this technique delivers better accurate or precise responses. In electric circuit analysis, Zhou [31] presented the DTM as a technique for addressing both linear and non-linear initial value problems. Kundu and Lee [32] employed DTM to assess the excessive heat transference aspects of annular fins when radiation effects were considered. Moradi et al. [33] used the fundamental algorithm of DTM to investigate the thermal performance of triangular fin with radiation and convection phenomenon. Several studies have recently used DTM to investigate the heat transferal aspects of various liquids through various geometries [34,35]. The DTM approach generates an analytical solution in a power series for differential equations. Moreover, for high values of a space variable, power series are typically insufficient. Pade approximants [36,37] are widely recognized for manipulating polynomial approximation into rational polynomial functions. Several researchers have recently used this DTM-Pade approximant method to solve different fluid flow models [38,39,40].
According to the abovementioned literature, DTM is an effective methodology for solving non-linear fin problems. The constant term in the boundary condition that appears after the transformation is specified and determined by employing boundary conditions in most research articles. In the present investigation, the Pade approximant technique is used to achieve this value, and the DTM-Pade approximant technique offers desirable solutions to heat transfer problems. The novelty of this investigation is to examine the temperature distribution in a longitudinal rectangular fin with internal heat generation for linearly and nonlinearly varying temperature-dependent thermal conductivity using the technique of DTM-Pade approximant. More explicitly, exponentially varying temperature-dependent thermal conductivity and linearly temperature-dependent thermal conductivity are taken into account in this inspection for analyzing the nature of temperature distribution within a straight fin. Also, the thermal distribution through a fin of material Aluminum Alloy 6061 (AA 6061) and Cast Iron is examined using ANSYS software.

2. Mathematical Formulation

A rectangular straight fin of length L , width W , thickness t with temperature-dependent internal heat generation q * ( T ) , thermal conductivity k * ( T ) whose faces are subjected to convection heat transfer of coefficient h * ( T ) at temperature T is considered in this analysis as shown in Figure 1. Since the cross-section’s temperature distribution is negligible and changes significantly in the longitudinal direction, the nature of the current problem is steady-state heat transfer. The thickness of a longitudinal fin t is comparatively lesser than its width W and length L , so it is presumed that the thermal attribute over the fin’s thickness and heat transference from its boundaries is insignificant compared to the heat-dissipation at its adjacent region.
With all these norms, the governing equation for the thermal distribution of a fin is given as (Aziz and Bouaziz [11]):
d d x [ k * ( T ) A c r * d T d x ] h * ( T ) P ( T T ) + A c r * q * ( T ) = 0
The corresponding boundary conditions are specified as:
x = 0   ;     d T d x = 0 , T = T b   ;     x = L .
The temperature dispersion through a fin is investigated in this study by determining two scenarios:
  • k * ( T ) and h * ( T ) are taken to be exponentially varying with temperature change.
  • k * ( T ) varies linearly with temperature and h * ( T ) exponentially varies with temperature change.
Case 1: Fin with exponential temperature-dependent k * ( T )
Here, k * ( T ) and h * ( T ) depend exponentially on the temperature change and are mathematically expressed as (see Moradi and Ahmadikia [41]):
k * ( T ) = k 0 e a 0 T   and   h * ( T ) = h 0 e b 0 T ,
and q * ( T ) is signified as:
q * ( T ) = q a [ 1 + ψ ( T T ) ]
For the sake of ease, the following non-dimensional terms are used.
θ = T T T b T ,   α 2 = P h b L 2 k a A c r * ,   γ = ψ ( T b T ) ,   X = x L , Q = q a A c r * P h b ( T b T ) ,   β = κ ( T b T ) .
And Equation (3) is indicated as:
k * ( T ) = k a e a θ   and   h * ( T ) = h b e b θ .
where,
a = a 0 ( T b T ) ,   k a = k 0 e a 0 T , b = b 0 ( T b T ) ,   h b = h 0 e h 0 T .
Using the above Equation from (3) to (7), Equation (1) together with the boundary conditions (2) can be derived as:
d d X [ e a θ d θ d X ] α 2 θ e b θ + α 2 Q ( 1 + γ   θ ) = 0 ,
Subjected to
d θ d X ( 0 ) = 0 ,   θ ( 1 ) = 1 .
Case 2: Fin with linear temperature-dependent k * ( T )
Here, the thermal conductivity k * ( T ) depends linearly on the temperature change, and the h * ( T ) depends exponentially on the temperature change and are mathematically expressed as:
k * ( T ) = k a [ 1 + κ ( T T ) ]   and   h * ( T ) = h 0 e b 0 T .
And Equation (10) is indicated as:
h * ( T ) = h b e b θ
where,
b = b 0 ( T b T ) ,   h b = h 0 e h 0 T .
Using the above Equation from (10) to (12) with Equation (5), Equation (1) together with the boundary conditions (2) can be derived as:
d 2 θ d X 2 + β θ d 2 θ d X 2 + β ( d θ d X ) 2 α 2 θ e b θ + α 2 Q ( 1 + γ   θ ) = 0
Subjected to
d θ d X ( 0 ) = 0 , θ ( 1 ) = 1 .

3. The Fundamental Concept of DTM-Pade Approximant

Zhou [31] introduced and applied the idea of DTM to solve linear and non-linear initial value problems in electric circuit investigation. This technique was later implemented to achieve analytical solutions to the partial differential equation (PDE), ODE, and other kinds of equations (see Moradi and Ahmadikia [41], Jawad [42]). Therefore, to derive an analytical solution for the prescribed differential equation, the basic definitions and preliminary results of DTM and Pade approximants are required and are found in the literature Wang et al. [18], Boyd [43], Rashidi et al. [38], and Ismail et al. [44].

4. Solution Procedure with DTM-Pade Approximant

Case 1: Fin with Exponential k * ( T )
Equation (8) is transformed using the DTM technique as:
s = 0 Κ F [ s ] ( Κ s + 1 ) ( Κ s + 2 )   Θ [ Κ s + 2 ] + a s = 0 Κ     ( s + 1 )   Θ [ s + 1 ]     l = 0 Κ s   F [ l ] ( Κ s l + 1 ) Θ [ Κ s l + 1 ] α 2 s = 0 Κ G [ s ] Θ [ Κ s ] + α 2 Q ( Δ [ Κ ] + γ   Θ [ Κ ] ) = 0
The transformed function F   [ s ] , F   [ l ] and G   [ s ] of e a θ and e b θ in the Equation (15) are evaluated using the basic properties of DTM (see Moradi and Ahmadikia [41]) and Θ [ Κ ] is the differential transform of θ ( X ) .
After applying DTM to Equation (9), the transformed relation is given as;
Θ [ 0 ] = A ,     Θ [ 1 ] = 0 , r = 0 Θ [ r ] = 1
Substituting Equation (16) and Κ = 0 ,   1 ,   2 ,   3 in Equation (15), The consecutive approximants are derived as follows:
Θ [ 2 ] = 1 2 e a A     α 2 ( A γ Q + A e b A Q ) Θ [ 3 ] = 0
Θ [ 4 ] = 1 24 ( e a A ) 2       α 4 ( 3 A 2 γ 2 Q 2 a 6 A 2 γ e b A Q a + A 2 γ e b A Q b + 3 A 2 ( e b A ) 2 a A 2 ( e b A ) 2 b A γ 2 Q 2 + 6 A γ Q 2 a + 2 A γ e b A Q 6 A e b A Q a + A e b A Q b A ( e b A ) 2 γ Q 2 + 3 Q 2 a + e b A Q )
Θ [ 5 ] = 0
Θ [ 6 ] = 1 720 ( e a A ) 3       α 6     ( 90 A 3 γ 2 e b A Q 2 a 2 3 A 3 γ 2 e b A Q 2 b 2 + 90 A 3 γ ( e b A ) 2 Q a 2 + 7 A 3 γ ( e b A ) 2 Q b 2 + 54 A 2 γ 2 e b A Q 2 a 8 A 2 γ 2 e b A Q 2 b 180 A 2 γ e b A Q 2 a 2 6 A 2 γ e b A Q 2 b 2 54 A 2 γ ( e b A ) 2 Q a + 16 A 2 γ ( e b A ) 2 Q b 36 A 2 ( e b A ) 2 Q a b + 72 A γ e b A Q 2 a 14 A γ e b A Q 2 b + 18 A e b A Q 2 a b + 90 A γ Q 3 a 2 3 A γ 2 e b A Q 2 90 A e b A Q 2 a 2 3 A e b A Q 2 b 2 + 3 A γ ( e b A ) 2 Q 36 A ( e b A ) 2 Q a + 14 A ( e b A ) 2 Q b + 18 A 3 γ 2 e b A Q 2 a b 36 A 3 γ ( e b A ) 2 Q a b + 36 A 2 γ e b A Q 2 a b + 30 Q 3 a 2 + γ 2 Q 3 A ( e b A ) 3 + ( e b A ) 2 Q 30 A 3 ( e b A ) 3 a 2 4 A 3 ( e b A ) 3 b 2 + A γ 3 Q 3 + 18 A 2 ( e b A ) 3 a 8 A 2 ( e b A ) 3 b 18 γ Q 3 a 2 γ e b A Q 2 + 18 e b A Q 2 a 6 e b A Q 2 b + 30 A 3 γ 3 Q 3 a 2 18 A 2 γ 3 Q 3 a + 90 A 2 γ 2 Q 3 a 2 + 18 A 3 ( e b A ) 3 a b + 90 A 2 ( e b A ) 2 Q a 2 + 7 A 2 ( e b A ) 2 Q b 2 36 A γ 2 Q 3 a )
And so on.
Where the constant A , which has to be evaluated by employing Equation (9). We have developed a power series of order 6 utilizing Equations (16)–(20), which is expressed as:
θ ( X )         =       A +   1 2 e a A     α 2 ( A γ Q + A e b A Q ) X 2 1 24 ( e a A ) 2       α 4 ( 3 A 2 γ 2 Q 2 a 6 A 2 γ e b A Q a + A 2 γ e b A Q b + 3 A 2 ( e b A ) 2 a A 2 ( e b A ) 2 b A γ 2 Q 2 + 6 A γ Q 2 a + 2 A γ e b A Q 6 A e b A Q a + A e b A Q b A ( e b A ) 2 γ Q 2 + 3 Q 2 a + e b A Q ) X 4   1 720 ( e a A ) 3       α 6     ( 90 A 3 γ 2 e b A Q 2 a 2 3 A 3 γ 2 e b A Q 2 b 2 + 90 A 3 γ ( e b A ) 2 Q a 2 + 7 A 3 γ ( e b A ) 2 Q b 2 + 54 A 2 γ 2 e b A Q 2 a 8 A 2 γ 2 e b A Q 2 b 180 A 2 γ e b A Q 2 a 2 6 A 2 γ e b A Q 2 b 2 54 A 2 γ ( e b A ) 2 Q a + 16 A 2 γ ( e b A ) 2 Q b 36 A 2 ( e b A ) 2 Q a b + 72 A γ e b A Q 2 a 14 A γ e b A Q 2 b + 18 A e b A Q 2 a b + 90 A γ Q 3 a 2 3 A γ 2 e b A Q 2 90 A e b A Q 2 a 2 3 A e b A Q 2 b 2 + 3 A γ ( e b A ) 2 Q 36 A ( e b A ) 2 Q a + 14 A ( e b A ) 2 Q b + 18 A 3 γ 2 e b A Q 2 a b 36 A 3 γ ( e b A ) 2 Q a b + 36 A 2 γ e b A Q 2 a b + 30 Q 3 a 2 + γ 2 Q 3 A ( e b A ) 3 + ( e b A ) 2 Q 30 A 3 ( e b A ) 3 a 2 4 A 3 ( e b A ) 3 b 2 + A γ 3 Q 3 + 18 A 2 ( e b A ) 3 a 8 A 2 ( e b A ) 3 b 18 γ Q 3 a 2 γ e b A Q 2 + 18 e b A Q 2 a 6 e b A Q 2 b + 30 A 3 γ 3 Q 3 a 2 18 A 2 γ 3 Q 3 a + 90 A 2 γ 2 Q 3 a 2 + 18 A 3 ( e b A ) 3 a b + 90 A 2 ( e b A ) 2 Q a 2 + 7 A 2 ( e b A ) 2 Q b 2 36 A γ 2 Q 3 a ) X 6
The A value is extracted by applying the Pade approximant to Equation (21) with Equation (9). The constant value A = 0.9954320629 , Q = 0.2 ,   a = 1 ,   α = 0.1 ,   b = 1 ,   γ = 0.1 are substituted in the series, we get
θ ( X ) = 0.9954320629 + 0.004570809020 X 2 0.000002873685247 X 4 +                               5.745191892 × 10 9 X 6 +
Case 2: Fin with Linear k * ( T )
Equation (13) is transformed using the DTM technique as:
( Κ + 1 ) ( Κ + 2 ) Θ [ Κ + 2 ] + β     s = 0 Κ ( s + 1 ) Θ [ s + 1 ]     ( Κ s + 1 ) Θ [ Κ s + 1 ] + β   s = 0 Κ   Θ [ Κ s ] ( s + 1 ) ( s + 2 ) Θ [ s + 2 ] α 2   s = 0 Κ   G [ s ] Θ [ Κ s ] + α 2   Q     ( γ   Θ [ Κ ] + Δ [ Κ ] ) = 0
applying DTM to Equation (14) yields
Θ [ 0 ] = B ,     Θ [ 1 ] = 0 ,   r = 0 Θ [ r ] = 1
Substituting Equation (24) and Κ = 0 ,   1 ,   2 ,   3 in Equation (23), The consecutive approximants are derived as follows:
Θ [ 2 ] = 1 2 ( B β + 1 ) α 2 ( B γ Q + B e b B Q )
Θ [ 3 ] = 0
Θ [ 4 ] = 1 24 ( B β + 1 ) 3 [ α 4 ( B γ Q + B e b B Q ) ( B 2 e b B b β + 2 B γ Q β + B e b B b 2 B e b B β Q γ + 3 Q β + e b B ) ]
Θ [ 5 ] = 0
Θ [ 6 ] = 1 720 ( B β + 1 ) 5   { α 6 ( B γ Q B e b B + Q ) ( ( e b B ) 2 + 3 B 3 e b B Q b 2 β 2 + 3 B 2 γ e b B Q b 2 + 56 B 2 γ e b B Q β 2 + 6 B 2 e b B Q b 2 β 12 B 2 e b B Q b β 2 + 8 B γ e b B Q b 32 B γ e b B Q β 6 B e b B Q b β 4 B 4 ( e b B ) 2 b 2 β 2 8 B 3 ( e b B ) 2 b 2 β + 10 B 3 ( e b B ) 2 b β 2 28 B 2 γ 2 Q 2 β 2 + 2 B 2 ( e b B ) 2 b β + 16 B γ 2 Q 2 β 72 B γ Q 2 β 2 + 3 B e b B Q b 2 + 72 B e b B Q β 2 + 3 B 4 γ e b B Q b 2 β 2 + 6 B 3 γ e b B Q b 2 β 10 B 3 γ e b B Q b β 2 2 B 2 γ e b B Q b β 4 B 2 ( e b B ) 2 b 2 28 B 2 ( e b B ) 2 β 2 8 B ( e b B ) 2 b + 16 B ( e b B ) 2 β + 18 γ Q 2 β + 2 γ e b B Q + 6 e b B Q b 18 e b B Q β 45 Q 2 β 2 γ 2 Q 2 ) }
And so on.
Where the constant B , which has to be evaluated by employing Equation (14). We have developed a power series of order 6 utilizing Equations (24)–(29), which is expressed as:
θ ( X )         =       A + 1 2 ( B β + 1 ) α 2 ( B γ Q + B e b B Q )   X 2 + 1 24 ( B β + 1 ) 3 [ α 4 ( B γ Q + B e b B Q ) ( B 2 e b B b β + 2 B γ Q β + B e b B b 2 B e b B β Q γ + 3 Q β + e b B ) ] X 4 + 1 720 ( B β + 1 ) 5   { α 6 ( B γ Q B e b B + Q ) ( ( e b B ) 2 + 3 B 3 e b B Q b 2 β 2 + 3 B 2 γ e b B Q b 2 + 56 B 2 γ e b B Q β 2 + 6 B 2 e b B Q b 2 β 12 B 2 e b B Q b β 2 + 8 B γ e b B Q b 32 B γ e b B Q β 6 B e b B Q b β 4 B 4 ( e b B ) 2 b 2 β 2 8 B 3 ( e b B ) 2 b 2 β + 10 B 3 ( e b B ) 2 b β 2 28 B 2 γ 2 Q 2 β 2 + 2 B 2 ( e b B ) 2 b β + 16 B γ 2 Q 2 β 72 B γ Q 2 β 2 + 3 B e b B Q b 2 + 72 B e b B Q β 2 + 3 B 4 γ e b B Q b 2 β 2 + 6 B 3 γ e b B Q b 2 β 10 B 3 γ e b B Q b β 2 2 B 2 γ e b B Q b β 4 B 2 ( e b B ) 2 b 2 28 B 2 ( e b B ) 2 β 2 8 B ( e b B ) 2 b + 16 B ( e b B ) 2 β + 18 γ Q 2 β + 2 γ e b B Q + 6 e b B Q b 18 e b B Q β 45 Q 2 β 2 γ 2 Q 2 ) } X 6 +
The B value is extracted by applying the Pade approximant to Equation (30) with Equation (14). The constant value B = 0.9888654429 , Q = 0.2 ,   β = 0.1 ,   α = 0.1 ,   b = 1 ,   γ = 0.1 are substituted in the series, we get
θ ( X ) = 0.9888654429 + 0.01109520280 X 2 + 0.00003921518999 X 4 +                           1.73779481310 7 X 6 +

5. Fin Efficiency

Fin efficiency is stated as the ratio of the fin heat transference rate to the rate that it would be if the complete fin were maintained at the base temperature and is mathematically expressed as:
η = k * A c r * d T d x x = L P h * L ( T b T )
The dimensionless form of fin efficiency is expressed as:
η = e a b θ ( 1 ) α 2

6. Result and Discussion

The significant consequence of several dimensionless variables, namely thermal conductivity parameter β ,   thermo-geometric parameter α , internal heat generation variable γ , and heat transfer parameter Q on the thermal dispersal of a longitudinal rectangular fin with exponential temperature-dependent heat transfer coefficient is examined by considering two cases of thermal conductivity. In the first case, it is taken to be exponential temperature-dependent, whereas in the other case, it is linearly dependent on temperature. The significance of these pertinent parameters on the non-dimensional temperature profile θ is explained in this section via graphical representation for both the aforementioned cases by using the DTM-Pade and RKF-45 methods.
Also, the effect of parameters on fin efficiency η is debriefed graphically. Table 1 is developed to compare the results of the current DTM-Pade solution for temperature profile with the existing literature (Languri et al. [45]) to authenticate the correctness of the solution. As per the tabulated result, there is an appropriate consistency between the result of the present work and those in Languri et al. [45], signifying that the solution obtained by the proposed technique is accurate. Table 2 and Table 3 are constructed to endorse the results achieved from the DTM-Pade and RKF-45 schemes for exponential temperature-dependent and linearly temperature-dependent thermal conductivity, respectively, and are close to each other. The present work is validated graphically with the existing work (Sun and Li [46]) as exhibited in Figure 2 by considering the values of the parameters a = 1 , b =   0 , α =   1.6 , Q   =   γ   = 0 and the graph shows the excellent agreement with the existing work.
Figure 3 illustrates the nature of thermal distribution θ for diverse α   ( 1 ,   1.1 ,   1.2 ,   1.3 ) values. This figure manifests that θ declines significantly with the impact of various α values. When the magnitude of α   upsurges, the convective heat transferal through the fin escalates, causing the fin temperature to drop. This results in a high heat transmission rate at the base due to the augmentation in the heat transfer ratio of convection to conduction ( h b / k b ). As an impact, the temperature through the fin decreases as α   augments. Thus, when the fin convective heat transfer improves, more heat transfer occurs by conduction through the fin, and as a result, the rate of heat transfer is augmented. The thermo-geometric parameter α   is substantial in assessing the amount of heat transmission from the fin since it accounted for the consequences of temperature reduction on heat transfer. The impact of dimensionless parameters γ   (   0.1 ,   0.3 ,   0.5 ,   0.7 ) on thermal distribution θ is portrayed in Figure 4. Here, θ enhances gradually for enhanced values of γ in the case of exponential temperature-dependent thermal conductivity. As these variable upsurges, the heat generation develops as a stronger function of the local fin temperature, resulting in an elevated temperature in the fin. Figure 5 exemplifies the behavior of θ with an increase of dimensionless heat transfer parameter Q   ( 0.2 ,   0.4 ,   0.6 ,   0.8 ) values. It is noted from this figure that enrichment in the scale of Q upsurges the thermal distribution θ .
Case 1: Exponential thermal conductivity and Exponential heat transfer coefficient.
Since the fin must disperse a large quantity of heat to the ambient fluid, excess heat production causes an improvement in the non-dimensional temperature. As demonstrated by the equations derived in this research, the consequence of internal heating on the heat transfer features of fins is substantially determined by the generation number Q , with the lower the generation number signifying higher heat transfer. Figure 6 manifests the deviation in the thermal distribution θ for diverse values of a   ( 2 ,   1 ,   0 ,   1 ) . As the a   ( 2 ,   1 ,   0 ,   1 ) increases, θ enhances remarkably for exponential k * ( T ) . Figure 7 displays the thermal dispersal deviations with the impact of b   (   1 ,   0 ,   1 ,   2 ) . The temperature distribution decreases for higher values of b   ( 1 ,   0 ,   1 ,   2 ) . Physically, with the increase of b   ( 2 ,   1 ,   0 ,   1 ) , the heat conduction coefficient of the fin upsurges, and consequently, the thermal resistance in the fin declines. Thus, even for the small temperature gradient, the large heat flow is generated by a fin. The features of non-dimensional parameters and their impact on fin efficiency η are represented via Figure 8, Figure 9 and Figure 10. The significant values of the fin efficiency show the performance of the heat exchanger, and better values of fin efficiency represent the better-quality performance of the fin. The major impact of α   (   0.6 ,   0.7 ,   0.8 ) and Q on η is shown in Figure 8. Here, the efficiency of fin diminutions with an upsurge of α with increased Q values in the case of exponential temperature-dependent thermal conductivity. Figure 9 reveals the consequence of Q   (   0.70 ,   0.75 ,   0.80 ) on η for increased values of α . In this figure, fin efficiency enriches for the enhanced magnitude of Q   (   0.70 ,   0.75 ,   0.80 ) . Figure 10 shows the influence of b   (   1 ,   0 ,   1 ) on η with diverse values of Q . Fin efficiency enhances for the greater magnitude of Q , and it decreases for different upsurge values of b   (   1 ,   0 ,   1 ) . Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15 elucidate the behavior of the temperature field with the significant influence of several non-dimensionless coefficients for linearly temperature-dependent thermal conductivity. Figure 11 is plotted to describe the impact of α   ( 0.8 ,   1.0 ,   1.2 ,   1.4 ) on θ . Here, higher α values decline θ for linearly temperature-dependent thermal conductivity. As the value of α increases, the convective heat transfer through the fin rises monotonically, while the temperature of the fin declines. As a result, the base possesses a high rate of heat transmission. The heat transfer ratio of convection to conduction ( h b / k b ) improves towards the fin base. Consequently, when the fin convective heat transfer improves, more heat is transferred by conduction in the fin, boosting the rate of heat transmission. The aspect of thermal distribution θ for enhancing values of dimensionless internal heat generation parameter γ   ( 0 ,   0.2 ,   0.4 ,   0.6 ) is depicted in Figure 12. As the magnitude of γ heightens, θ of the fin improves remarkably. The effect of dimensionless heat transfer parameter Q on the behavior of θ is illustrated via Figure 13. Here, the thermal distribution θ upsurges for improved values of Q   ( 0.1 ,   0.3 ,   0.5 ,   0.7 ) in the case of linearly temperature-dependent k * ( T ) . Furthermore, the lower generation number causes the higher heat transfer. Figure 14 reveals the deviation in the thermal distribution θ of the fin with the influence of the thermal conductivity parameter β . This figure denotes that rising values of β   ( 0 ,   0.2 ,   0.4 ,   0.6 ) is accountable for the augmentation of θ . Enhancement in β enriches heat conduction from the fin base, and the temperature inside the fin increases. Physically, the tendency of the fin to conduct heat significantly upsurge as this parameter intensifies so does the local fin temperature. The thermal distribution deviation for increasing magnitude of b   (   1 ,   0 ,   1 ,   2 ) is exhibited in Figure 15. As the b   (   1 ,   0 ,   1 ,   2 ) values increase, the temperature inside the fin decreases for linearly temperature-dependent k * ( T ) . Figure 16 portrays the comparison of temperature distribution in a longitudinal fin for exponential temperature-dependent and linearly temperature-dependent thermal conductivity with the impact of Q   ( 0.2 ,   0.3 ,   0.4 ,   0.5 ) . Thermal distribution is more in a fin with exponential temperature-dependent. compared to a fin with linearly temperature-dependent k * ( T ) for a higher magnitude of Q   ( 0.2 ,   0.3 ,   0.4 ,   0.5 ) .
Case 2: Linear thermal conductivity and Exponential heat transfer coefficient.

7. Inspecting the Thermal Behavior of a Longitudinal Fin Using ANSYS

ANSYS is a general-purpose software package that affords a cost-effective approach to inspecting the better performance of mechanical and fluid models or processes in a virtual environment. This software implements governing equations for corresponding problems to study the behavior of problems and solves them efficiently. The obtained results are presented in tabular or graphical forms.
To inspect the thermal behavior, a longitudinal fin is modeled with the following assumptions:
  • Since aluminum is an excellent thermal and electrical conductor, Aluminum Alloy 6061 (AA 6061) and Cast Iron with constant thermal conductivity 300 W / m   K and 55 W / m   K are taken as fin materials.
  • One-dimensional heat conduction is considered along the longitudinal direction.
  • The convective heat transfer coefficient (39.9 W / m 2   K ) is considered over the complete fin surface.
  • The temperature at the fin base is 550 K, and the ambient temperature is 283 K.
The temperature distribution of a longitudinal fin made up of AA 6061, and Cast Iron can be assessed with an ANSYS workbench. For the longitudinal fin made up of AA 6061, the maximum temperature at the prime surface is 550 K. By contrast, the minimum temperature is 545.58 K observed at the end of the fin, and 548.9 K is the average temperature. From this analysis, it is found that the extreme temperature of 550 K is at the source, and the minimum temperature of 527.28 K will be at the end of the fin for the fin of material Cast Iron. It is also noticed that the average temperature through the fin is 544.33 K. The temperature gradually decreases from one end to another end, that is from base to tip end based on the given heat transfer coefficient as shown in Figure 17, Figure 18, Figure 19 and Figure 20.

8. Final Remarks

The temperature dispersal of a longitudinal rectangular fin with exponential temperature-dependent heat transfer coefficient and internal heat generation is inspected by considering exponential temperature-dependent k * ( T ) with the aid of the DTM-Pade approximant. The significant results of the present scrutiny are as follows:
  • Enhancement in the scale of thermo-geometric parameters reduces temperature dispersal in a fin for both cases.
  • Temperature distribution enriches for a larger magnitude of thermal conductivity parameter in the case of linear temperature-dependent thermal conductivity.
  • Larger values of the internal heat generation and heat transfer parameter upsurge the thermal distribution in both cases.
  • The efficiency of a fin varies with prescribed non-dimensional thermal parameters under internal heat generation.
  • The thermal distribution of a longitudinal fin is studied using ANSYS software by considering the material of the fin body as AA 6061 and Cast Iron. The temperature is higher at the base, decreasing monotonically towards the fin tip.
  • The analytical solution and numerical results obtained by the DTM-Pade approximant afford higher accuracy than other techniques.

Author Contributions

Conceptualization, B.C.P. and I.E.S.; methodology, B.C.P. and I.E.S.; software, G.S. and R.S.V.K.; validation, B.C.P. and I.E.S.; formal analysis, G.S. and R.S.V.K.; investigation, G.S.; resources, I.E.S.; data curation, B.C.P.; writing original draft preparation, G.S. and R.S.V.K.; writing review and editing, I.E.S. and B.C.P.; visualization, G.S. and R.N.K.; supervision, B.C.P. and I.E.S.; project administration, I.E.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

t Fin’s thickness ( m )
X Length (dimensionless)
A c r * Fin’s cross-sectional area ( m 2 )
k a Thermal conductivity at ambient temperature ( W   m 1   K 1 )
Q Dimensionless heat transfer
b 0 Exponential indexes of convection heat transfer coefficient
β Variable thermal conductivity(dimensionless)
h * Convective heat transfer coefficient ( W   m 2   K 1 )
α Thermo-geometric parameter
θ Non-dimensional temperature
a 0 Exponential indexes of thermal conductivity
x Fin axial distance ( m )
κ Thermal conductivity variation parameter ( K 1 )
T b Base temperature ( K )
W Width ( m )
L Length ( m )
T Ambient temperature ( K )
q * ( T ) Uniform internal heat generation ( W   m 3   )
h b Heat transfer coefficient at the fin’s base ( W   m 2   K 1 )
η Fin efficiency
k * ( T ) Thermal conductivity ( W   m 1   K 1 )
h 0 Reference value of convection heat transfer coefficient ( W   m 2   K 1 )
γ Dimensionless internal heat generation parameter
k 0 Reference values of thermal Conductivity ( W   m 1   K 1 )
P Perimeter ( m )
T Temperature ( K )

References

  1. Khan, S.U.; Tlili, I. Significance of Activation Energy and Effective Prandtl Number in Accelerated Flow of Jeffrey Nanoparticles With Gyrotactic Microorganisms. J. Energy Resour. Technol. 2020, 142, 112101. [Google Scholar] [CrossRef]
  2. Hosseinzadeh, K.; Moghaddam, M.E.; Asadi, A.; Mogharrebi, A.; Ganji, D. Effect of internal fins along with Hybrid Nano-Particles on solid process in star shape triplex Latent Heat Thermal Energy Storage System by numerical simulation. Renew. Energy 2020, 154, 497–507. [Google Scholar] [CrossRef]
  3. Khan, S.U.; Al-Khaled, K.; Aldabesh, A.; Awais, M.; Tlili, I. Bioconvection flow in accelerated couple stress nanoparticles with activation energy: Bio-fuel applications. Sci. Rep. 2021, 11, 3331. [Google Scholar] [CrossRef] [PubMed]
  4. Hamid, A.; Kumar, R.N.; Gowda, R.J.P.; Khan, S.U.; Khan, M.I.; Prasannakumara, B.C.; Muhammad, T. Impact of Hall current and homogenous–heterogenous reactions on MHD flow of GO-MoS2/water (H2O)-ethylene glycol (C2H6O2) hybrid nanofluid past a vertical stretching surface. Waves Random Complex Media 2021, 1–18. [Google Scholar] [CrossRef]
  5. Khan, M.I.; Qayyum, S.; Chu, Y.-M.; Khan, N.B.; Kadry, S. Transportation of Marangoni convection and irregular heat source in entropy optimized dissipative flow. Int. Commun. Heat Mass Transf. 2021, 120, 105031. [Google Scholar] [CrossRef]
  6. Madhukesh, J.K.; Kumar, R.N.; Gowda, R.J.P.; Prasannakumara, B.C.; Ramesh, G.K.; Khan, M.I.; Khan, S.U.; Chu, Y.-M. Numerical simulation of AA7072-AA7075/water-based hybrid nanofluid flow over a curved stretching sheet with Newtonian heating: A non-Fourier heat flux model approach. J. Mol. Liq. 2021, 335, 116103. [Google Scholar] [CrossRef]
  7. Raza, A.; Khan, S.U.; Al-Khaled, K.; Khan, M.I.; Haq, A.U.; Alotaibi, F.; Allah, A.M.A.; Qayyum, S. A fractional model for the kerosene oil and water-based Casson nanofluid with inclined magnetic force. Chem. Phys. Lett. 2021, 787, 139277. [Google Scholar] [CrossRef]
  8. Hosseinzadeh, K.; Mogharrebi, A.; Asadi, A.; Paikar, M.; Ganji, D. Effect of fin and hybrid nano-particles on solid process in hexagonal triplex Latent Heat Thermal Energy Storage System. J. Mol. Liq. 2020, 300, 112347. [Google Scholar] [CrossRef]
  9. Hosseinzadeh, K.; Montazer, E.; Shafii, M.B.; Ganji, A. Solidification enhancement in triplex thermal energy storage system via triplets fins configuration and hybrid nanoparticles. J. Energy Storage 2021, 34, 102177. [Google Scholar] [CrossRef]
  10. Kundu, B. Performance and optimum design analysis of longitudinal and pin fins with simultaneous heat and mass transfer: Unified and comparative investigations. Appl. Therm. Eng. 2007, 27, 976–987. [Google Scholar] [CrossRef]
  11. Aziz, A.; Bouaziz, M. A least squares method for a longitudinal fin with temperature dependent internal heat generation and thermal conductivity. Energy Convers. Manag. 2011, 52, 2876–2882. [Google Scholar] [CrossRef]
  12. Najafabadi, M.F.; Rostami, H.T.; Hosseinzadeh, K.; Ganji, D.D. Thermal analysis of a moving fin using the radial basis function approximation. Heat Transf. 2021, 50, 7553–7567. [Google Scholar] [CrossRef]
  13. Sowmya, G.; Sarris, I.E.; Vishalakshi, C.S.; Kumar, R.S.V.; Prasannakumara, B.C. Analysis of Transient Thermal Distribution in a Convective–Radiative Moving Rod Using Two-Dimensional Differential Transform Method with Multivariate Pade Approximant. Symmetry 2021, 13, 1793. [Google Scholar] [CrossRef]
  14. Sowmya, G.; Kumar, R.S.V.; Alsulami, M.D.; Prasannakumara, B.C. Thermal stress and temperature distribution of an annular fin with variable temperature-dependent thermal properties and magnetic field using DTM-Pade approximant. Waves Random Complex Media 2022, 1–29. [Google Scholar] [CrossRef]
  15. Zhang, H.J.; Sun, K.F. Conduction from Longitudinal Fin of Rectangular Profile with Exponential Vary Heat Transfer Coefficient. Adv. Mater. Res. 2012, 614–615, 311–314. [Google Scholar]
  16. Buikis, A.; Pagodkina, I. Comparison of Analytical and Numerical Solutions for a Two-Dimensinal Longitudal Fin of Rectangular Profile. Latv. J. Phys. Tech. Sci. 1996. [Google Scholar] [CrossRef]
  17. Kezzar, M.; Tabet, I.; Eid, M.R. A new analytical solution of longitudinal fin with variable heat generation and thermal conductivity using DRA. Eur. Phys. J. Plus 2020, 135, 120. [Google Scholar] [CrossRef]
  18. Wang, F.; Kumar, R.V.; Sowmya, G.; El-Zahar, E.R.; Prasannakumara, B.; Khan, M.I.; Khan, S.U.; Malik, M.; Xia, W.-F. LSM and DTM-Pade approximation for the combined impacts of convective and radiative heat transfer on an inclined porous longitudinal fin. Case Stud. Therm. Eng. 2022, 101846. [Google Scholar] [CrossRef]
  19. Alhejaili, W.; Kumar, R.V.; El-Zahar, E.R.; Sowmya, G.; Prasannakumara, B.; Khan, M.I.; Yogeesha, K.; Qayyum, S. Analytical solution for temperature equation of a fin problem with variable temperature-dependent thermal properties: Application of LSM and DTM-Pade approximant. Chem. Phys. Lett. 2022, 793, 139409. [Google Scholar] [CrossRef]
  20. Torabi, M.; Aziz, A.; Zhang, K. A comparative study of longitudinal fins of rectangular, trapezoidal and concave parabolic profiles with multiple nonlinearities. Energy 2013, 51, 243–256. [Google Scholar] [CrossRef]
  21. Mallick, A.; Ghosal, S.; Sarkar, P.K.; Ranjan, R. Homotopy Perturbation Method for Thermal Stresses in an AnnularFinwith Variable Thermal Conductivity. J. Therm. Stress. 2014, 38, 110–132. [Google Scholar] [CrossRef]
  22. Atouei, S.; Hosseinzadeh, K.; Hatami, M.; Ghasemi, S.E.; Sahebi, S.; Ganji, D. Heat transfer study on convective–radiative semi-spherical fins with temperature-dependent properties and heat generation using efficient computational methods. Appl. Therm. Eng. 2015, 89, 299–305. [Google Scholar] [CrossRef]
  23. Turkyilmazoglu, M. Heat transfer from moving exponential fins exposed to heat generation. Int. J. Heat Mass Transf. 2018, 116, 346–351. [Google Scholar] [CrossRef]
  24. Padmanabhan, S.; Thiagarajan, S.; Kumar, A.D.R.; Prabhakaran, D.; Raju, M. Investigation of temperature distribution of fin profiles using analytical and CFD analysis. Mater. Today: Proc. 2021, 44, 3550–3556. [Google Scholar] [CrossRef]
  25. Hosseinzadeh, S.; Hasibi, A.; Ganji, D. Thermal analysis of moving porous fin wetted by hybrid nanofluid with trapezoidal, concave parabolic and convex cross sections. Case Stud. Therm. Eng. 2022, 30, 101757. [Google Scholar] [CrossRef]
  26. Minkler, W.S.; Rouleau, W.T. The Effects of Internal Heat Generation on Heat Transfer in Thin Fins. Nucl. Sci. Eng. 1960, 7, 400–406. [Google Scholar] [CrossRef]
  27. Venkitesh, V.; Mallick, A. Thermal analysis of a convective–conductive–radiative annular porous fin with variable thermal parameters and internal heat generation. J. Therm. Anal. 2020, 147, 1519–1533. [Google Scholar] [CrossRef]
  28. Majhi, T.; Kundu, B. New Approach for Determining Fin Performances of an Annular Disc Fin with Internal Heat Generation. In Advances in Mechanical Engineering; Springer: Singapore, 2020; pp. 1033–1043. [Google Scholar] [CrossRef]
  29. Das, R.; Kundu, B. Prediction of Heat-Generation and Electromagnetic Parameters from Temperature Response in Porous Fins. J. Thermophys. Heat Transf. 2021, 35, 761–769. [Google Scholar] [CrossRef]
  30. Jayaprakash, M.; Alzahrani, H.A.; Sowmya, G.; Kumar, R.V.; Malik, M.; Alsaiari, A.; Prasannakumara, B. Thermal distribution through a moving longitudinal trapezoidal fin with variable temperature-dependent thermal properties using DTM-Pade approximant. Case Stud. Therm. Eng. 2021, 28, 101697. [Google Scholar] [CrossRef]
  31. Zhou, J.K. Differential Transformation and Its Applications for Electrical Circuits; Huazhong University Press: Wuhan, China, 1986. [Google Scholar]
  32. Kundu, B.; Lee, K.-S. Analytical tools for calculating the maximum heat transfer of annular stepped fins with internal heat generation and radiation effects. Energy 2014, 76, 733–748. [Google Scholar] [CrossRef]
  33. Moradi, A.; Hayat, T.; Alsaedi, A. Convection-radiation thermal analysis of triangular porous fins with temperature-dependent thermal conductivity by DTM. Energy Convers. Manag. 2014, 77, 70–77. [Google Scholar] [CrossRef]
  34. Christopher, A.J.; Magesh, N.; Gowda, R.J.P.; Kumar, R.N.; Kumar, R.S.V. Hybrid nanofluid flow over a stretched cylinder with the impact of homogeneous–heterogeneous reactions and Cattaneo–Christov heat flux: Series solution and numerical simulation. Heat Transf. 2021, 50, 3800–3821. [Google Scholar] [CrossRef]
  35. Kundu, B.; Yook, S.-J. An accurate approach for thermal analysis of porous longitudinal, spine and radial fins with all nonlinearity effects—Analytical and unified assessment. Appl. Math. Comput. 2021, 402, 126124. [Google Scholar] [CrossRef]
  36. Baker, G.A., Jr. Essentials of Padé Approximants; Academic Press: New York, NY, USA, 1975; 317p, ISBN 0-12-074885-X. [Google Scholar]
  37. Baker, G.A.; Graves-Morris, P.R. Padé Approximants: Extensions and Applications; Addison-Wesley: Boston, MA, USA, 1981; Volume 2. [Google Scholar]
  38. Rashidi, M.M.; Keimanesh, M. Using Differential Transform Method and Padé Approximant for Solving MHD Flow in a Laminar Liquid Film from a Horizontal Stretching Surface. Math. Probl. Eng. 2010, 2010, 491319. [Google Scholar] [CrossRef]
  39. Rashidi, M.M.; Laraqi, N.; Sadri, S.M. A novel analytical solution of mixed convection about an inclined flat plate embedded in a porous medium using the DTM-Padé. Int. J. Therm. Sci. 2010, 49, 2405–2412. [Google Scholar] [CrossRef]
  40. Saha, D.; Sengupta, S. Dual DTM-Padé approximations on free convection MHD mass transfer flow of nanofluid through a stretching sheet in presence of Soret and Dufour Phenomena. WSEAS Trans. Fluid Mech. 2020, 15, 23–40. [Google Scholar] [CrossRef]
  41. Moradi, A.; Ahmadikia, H. Investigation of Effect Thermal Conductivity on Straight Fin Performance with DTM. Int. J. Eng. Appl. Sci. 2011, 3, 42–54. [Google Scholar]
  42. Jawad, A.J.M. Differential Transformation Method for Solving Nonlinear Heat Transfer Equations. Al-Rafidain Univ. Coll. Sci. 2013, 32, 180–194. [Google Scholar]
  43. Boyd, J.P. Padé approximant algorithm for solving nonlinear ordinary differential equation boundary value problems on an unbounded domain. Comput. Phys. 1997, 11, 299. [Google Scholar] [CrossRef]
  44. Ismail, T.; Kezzar, M.; Touafe, K.; Bellel, N.; Gherieb, S.; Khelifa, A.; Adouane, M. Adomian Decomposition Method And Padé Approximation To Determine Fin Efficiency Of Convective Solar Air Collector In Straight Fins. Int. J. Math. Model. Comput. 2015, 5, 335–346. [Google Scholar]
  45. Languri, E.M.; Ganji, D.D.; Jamshidi, N. Variational Iteration and Homotopy perturbation methods for fin efficiency of convective straight fins with temperature dependent thermal conductivity. In Proceedings of the 5th WSEAS international conference On FLUID MECHANICS (Fluids 08), Acapulco, Mexico, 25–27 January 2008; Volume 25. [Google Scholar]
  46. Sun, S.-W.; Li, X.-F. Exact solution of the nonlinear fin problem with exponentially temperature-dependent thermal conductivity and heat transfer coefficient. Pramana 2020, 94, 94. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of a longitudinal fin.
Figure 1. Schematic representation of a longitudinal fin.
Symmetry 14 00690 g001
Figure 2. Validation of the present analysis with existing work (Sun and Li [46]) for a = 1 , b = 0 , α = 1.6 , Q = γ = 0 .
Figure 2. Validation of the present analysis with existing work (Sun and Li [46]) for a = 1 , b = 0 , α = 1.6 , Q = γ = 0 .
Symmetry 14 00690 g002
Figure 3. Impact of α on θ .
Figure 3. Impact of α on θ .
Symmetry 14 00690 g003
Figure 4. Impact of γ on θ .
Figure 4. Impact of γ on θ .
Symmetry 14 00690 g004
Figure 5. Impact of Q on θ .
Figure 5. Impact of Q on θ .
Symmetry 14 00690 g005
Figure 6. Impact of a on θ .
Figure 6. Impact of a on θ .
Symmetry 14 00690 g006
Figure 7. Impact of b on θ .
Figure 7. Impact of b on θ .
Symmetry 14 00690 g007
Figure 8. Influence of α on fin efficiency.
Figure 8. Influence of α on fin efficiency.
Symmetry 14 00690 g008
Figure 9. Impact of Q on fin efficiency.
Figure 9. Impact of Q on fin efficiency.
Symmetry 14 00690 g009
Figure 10. Impact of b on fin efficiency.
Figure 10. Impact of b on fin efficiency.
Symmetry 14 00690 g010
Figure 11. Impact of α on θ .
Figure 11. Impact of α on θ .
Symmetry 14 00690 g011
Figure 12. Impact of γ on θ .
Figure 12. Impact of γ on θ .
Symmetry 14 00690 g012
Figure 13. Impact of Q on θ .
Figure 13. Impact of Q on θ .
Symmetry 14 00690 g013
Figure 14. Impact of β on θ .
Figure 14. Impact of β on θ .
Symmetry 14 00690 g014
Figure 15. Impact of b on θ .
Figure 15. Impact of b on θ .
Symmetry 14 00690 g015
Figure 16. Comparing the temperature distribution of a fin with exponential temperature-dependent and linearly temperature-dependent thermal conductivity.
Figure 16. Comparing the temperature distribution of a fin with exponential temperature-dependent and linearly temperature-dependent thermal conductivity.
Symmetry 14 00690 g016
Figure 17. Temperature distribution for Aluminum Alloy 6061.
Figure 17. Temperature distribution for Aluminum Alloy 6061.
Symmetry 14 00690 g017
Figure 18. Temperature varying along the length of the fin for Aluminum Alloy 6061.
Figure 18. Temperature varying along the length of the fin for Aluminum Alloy 6061.
Symmetry 14 00690 g018
Figure 19. Temperature distribution for Cast Iron.
Figure 19. Temperature distribution for Cast Iron.
Symmetry 14 00690 g019
Figure 20. Temperature varying along the length of the fin for Cast Iron.
Figure 20. Temperature varying along the length of the fin for Cast Iron.
Symmetry 14 00690 g020
Table 1. Comparison of θ ( X ) obtained by DTM-Pade with existing literature (Languri et al. [45]) with constant values of β = 0 , Q = 0 ,   α = 0.5 ,   b = 0 ,     γ = 0 .
Table 1. Comparison of θ ( X ) obtained by DTM-Pade with existing literature (Languri et al. [45]) with constant values of β = 0 , Q = 0 ,   α = 0.5 ,   b = 0 ,     γ = 0 .
θ ( X )
X HPM
(Languri et al. [45])
VIM
(Languri et al. [45])
DTM-PadeError
00.8868190.8868190.8868180.000001
0.20.8912570.8912570.8912560.000001
0.40.9046140.9046140.9046140.000000
0.60.9270260.9270260.9270260.000000
0.80.9587150.9587150.9587150.000000
1.01.0000001.0000001.0000000.000000
Table 2. Comparison of θ ( X ) for a = 1 , a = 0 and a = 1 with constant values of Q = 0.8 ,   α = 1 ,   b = 1 ,     γ = 0.5 obtained by RKF-45 and DTM-Pade.
Table 2. Comparison of θ ( X ) for a = 1 , a = 0 and a = 1 with constant values of Q = 0.8 ,   α = 1 ,   b = 1 ,     γ = 0.5 obtained by RKF-45 and DTM-Pade.
θ ( X )
X a = 1 a = 0 a = 1
RKF-45Present Result % Error RKF-45Present Result % Error RKF-45Present Result % Error
00.6512680.6470210.42470.7268590.7246230.22360.8259040.8253020.0602
0.10.6530870.6487190.43680.7289280.7266580.22700.8275600.8269540.0606
0.20.6586600.6539210.47390.7352030.7328240.23790.8325390.8319200.0619
0.30.6683570.6629620.53950.7458850.7433210.25640.8408720.8402340.0638
0.40.6828460.6764430.64030.7613220.7584890.28330.8526140.8519480.0666
0.50.7031900.6952900.79000.7820320.7788280.32040.8678380.8671360.0702
Table 3. Comparison of θ ( X ) for b = 1 , b = 0 and b = 1 with constant values of Q = 0.8 ,   α = 1 ,       γ = 0.5 ,   β = 0.1 obtained by RKF-45 and DTM-Pade.
Table 3. Comparison of θ ( X ) for b = 1 , b = 0 and b = 1 with constant values of Q = 0.8 ,   α = 1 ,       γ = 0.5 ,   β = 0.1 obtained by RKF-45 and DTM-Pade.
θ ( X )
X b = 1 b = 0 b = 1
RKF-45Present Result % Error RKF-45Present Result % Error RKF-45PresentResult % Error
01.4459731.4457630.02101.0737941.0738620.00680.7357370.7336810.2056
0.11.4414401.4412310.02091.0730901.0731580.00680.7377960.7357100.2086
0.21.4278511.4276420.02091.0709761.0710450.00690.7440360.7418570.2179
0.31.4052331.4050260.02071.0674381.0675080.00700.7546360.7523000.2336
0.41.3736311.3734270.02041.0624571.0625280.00710.7699120.7673460.2566
0.51.3331091.3329090.02001.0560021.0560750.00730.7903230.7874440.2879
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kumar, R.S.V.; Kumar, R.N.; Sowmya, G.; Prasannakumara, B.C.; Sarris, I.E. Exploration of Temperature Distribution through a Longitudinal Rectangular Fin with Linear and Exponential Temperature-Dependent Thermal Conductivity Using DTM-Pade Approximant. Symmetry 2022, 14, 690. https://doi.org/10.3390/sym14040690

AMA Style

Kumar RSV, Kumar RN, Sowmya G, Prasannakumara BC, Sarris IE. Exploration of Temperature Distribution through a Longitudinal Rectangular Fin with Linear and Exponential Temperature-Dependent Thermal Conductivity Using DTM-Pade Approximant. Symmetry. 2022; 14(4):690. https://doi.org/10.3390/sym14040690

Chicago/Turabian Style

Kumar, Ravikumar Shashikala Varun, Rangaswamy Naveen Kumar, Ganeshappa Sowmya, Ballajja Chandrappa Prasannakumara, and Ioannis E. Sarris. 2022. "Exploration of Temperature Distribution through a Longitudinal Rectangular Fin with Linear and Exponential Temperature-Dependent Thermal Conductivity Using DTM-Pade Approximant" Symmetry 14, no. 4: 690. https://doi.org/10.3390/sym14040690

APA Style

Kumar, R. S. V., Kumar, R. N., Sowmya, G., Prasannakumara, B. C., & Sarris, I. E. (2022). Exploration of Temperature Distribution through a Longitudinal Rectangular Fin with Linear and Exponential Temperature-Dependent Thermal Conductivity Using DTM-Pade Approximant. Symmetry, 14(4), 690. https://doi.org/10.3390/sym14040690

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop