Adomian decomposition sumudu transform method for solving a solid and porous fin with temperature dependent internal heat generation.

In this paper, Adomian decomposition sumudu transform method is introduced and used to solve the temperature distribution in a solid and porous fin with the temperature dependent internal heat generation for a fractional order energy balance equation. In this study, we assume heat generation as a variable of fin temperature for solid and porous fin and the heat transfer through porous media is simulated by using Darcy's model. The results are presented for the temperature distribution for the range of values of parameters appeared in the mathematical formulation and also compared with numerical solutions in order to verify the accuracy of the proposed method. It is found that the proposed method is in good agreement with direct numerical solution.

Background

Fins are commonly used to facilitate the dissipation of heat from a heated wall to the surrounding environment. Examples of fin are the radiator in vehicles and heat exchangers in power plants. In electrical devices like motors and transformers, the generated heat can be efficiently transferred. In the study of heat transfer, fin is a surface, made by metallic material which is used to increase the rate of heat transfer to the environment. The rate of heat transfer depends on the surface area of the fin. Fins are extensively used to improve the rate of heat dissipation from a hot surface, especially in thermal engineering applications (Bergman et al.; Nield and Bejan).

In many thermal engineering applications, convective flow through porous media is mandatory for the investigation. Several numerical and analytical revisions has been conducted so far to afford a profound understanding of the transport system of the heat transfer inside the porous medium. Generally, high thermal conductivity porous substrates are employed to improve the rate of forced convection heat transfer in many engineering applications such as reactor heat exchangers, solar collectors and in cooling process (Alkam and Al-Nimr 1999). However, heat transfer in porous fins has attracted a lot of attention of researchers with a wide range of it’s applications, especially in recent years. Kiwan and Al-Nimr (2001) was the first person who introduced the concept of fins made of porous materials by introducing Darcy’s model (Kiwan 2007; Kiwan and Zeitoun 2008).

Now a days, Heat exchanger industries are looking for more compact and cost–effective heat exchanger manufacturing techniques which leads to use porous fins in enhance heat transfer (Kiwan 2007). The heat-transfer enhancement between two parallel-plate channels was investigated by adding porous fin through the channel (Hamdan and Moh’d 2010) and by adding porous insert to one side of the duct walls (Hamdan et al. 2000). Alkam et al. (2002) investigated the thermal analysis of natural convection porous fins. They studied all the geometric flow parameters that influence the temperature distribution in to a single parameter specified S. They considered three cases: the infinite fin, a finite fin with an insulated tip and a finite fin with uninsulated tip. Similarly Gorla and Bakier (2011) discussed the thermal analysis of natural convection and radiation in the porous fin and showed that the radiation transfers more heat than a similar model without radiation. Hatami et al. (2013) studied the heat transfer through porous fin with different porous material and compared their results with the Differential Transform Method, Collocation Method and Least Square Method. They Hatami and Ganji (2013) also studied the thermal performance of circular convective–radiative porous fins with different section, shapes and materials. Ghasemi et al. (2014) used the Differential Transform Method for solving the nonlinear temperature distribution in solid and porous fin with temperature dependent internal heat generation. Patel and Meher (2015a, b) studied the fractional solution of longitudinal porous fin for the case of temperature distribution, efficiency and effectiveness and also analysed the variation of temperature distribution for a straight rectangular fin with power-law temperature dependent surface heat flux by using Adomian decomposition sumudu transform method.

It is revealed that, the concept of fractional derivative is more suitable for modeling real world problem than the local derivative. Many researchers have devoted their attention in developing new definition of fractional derivative (Atangana 2016). Baskonus and Bulut (2015) applied the fractional Adams–Bashforth–Moulton Method for obtaining the numerical solutions of some linear and nonlinear fractional ordinary differential equations. Baskonus and Bulut (2015) studied it to obtain some new analytical solutions to the (1 + 1)-dimensional nonlinear Dispersive Modified Benjamin Bona Mahony equation by using modified exp-function method. Roshid et al. (2014) studied solitary wave solutions for vakhnenko-parkes equation via exp-function and Exp —expansion method. Also they Roshid et al. (2014) studied traveling wave solutions of nonlinear partial differential equation via new extended —expansion method.

In the present study, we fractionalize the energy balance equation in order to understand the anomalous behavior of this system and to find the temperature distribution in solid and porous fin by using Adomian decomposition sumudu transform method.

Preliminaries

Definitions of Caputo fractional derivative

In this part of the paper it would be useful to introduce some definitions and properties of the fractional calculus theory. There are several definitions of fractional derivatives of order (Miller and Ross; Srivastava et al. 2014). The two most commonly used definitions are Riemann-Liouville and Caputo.

Definition

The Riemann-Liouville fractional integral of fractional order is defined as Miller and Rosswhere is the set of positive real numbers and is the gamma function.

The fractional derivative of f(t) in the Caputo sense is defined bywhere the parameter α is the order of the derivative and is allowed to be real or even complex. Here N is the set of natural numbers. In this paper only real and positive α will be considered. The properties underpinning the use of the Caputo derivative can be found in Atangana (2016), Atangana and Alqahtani (2016) and Atangana and Baleanu (2016).

Basics of sumudu transform method

The sumudu transform is a new integral transform (Kadem and Baleanu 2012; Atangana and Baleanu 2013; Jarad et al. 2012) which is a little known and not widely used whose defined for the functions of exponential order.

The sumudu transform of a function f(t), defined for all real numbers , is the function F(u), defined by Watugala (1993)

Sumudu transform of function derivatives is defined as Belgacem and Karaballi (2006)

Sumudu transform of the Caputo fractional derivative is defined as Belgacem and Karaballi (2006)

Theorem

LetG(u) be the sumudu transform off(t) such thatFor the proof see Belgacem and Karaballi (2006).

is a meromorphic function, with singularities having

There exist a circular regionwith radiusRand positive constantsMandKwith, then the functionf(t) is given by

Formulation of Adomian decomposition sumudu transform method (ADSTM)

Consider a general fractional nonlinear nonhomogeneous differential equation with satisfying the initial condition of the form, as,subject to the initial conditionwhere denotes without loss of generality the Caputo fraction derivative operator, R is the linear differential operator, N represents the nonlinear differential operator, and is the source term.

On applying the sumudu transform and the Caputo fractional derivative in Eq. (7), it obtainsSumudu inverse transform of Eq. (9), giveswhere is prescribed initial conditions.

Using Adomian decomposition method which obtains the approximate solution of Eq. (10) in a series form as Adomianand the nonlinear term can be expressed as a sum of Adomian Polynomialswhere the Adomian polynomials , depends upon the solution components , can be defined as followsOn substituting Eqs. (11) and (12) in Eq. (10), it obtainsThe resulting Eq. (14) is the coupling of the Adomian decomposition method and the sumudu transform. On comparing the coefficients of like powers of , the recursive relation of Eq  (14) can be written asSince is a rapidly converging series, the partial sum is our approximant to the solution (Adomian).

Problem description

Here we considered two cases, namely (1) solid fin and (2) porous fin to study the fin temperature distribution in longitudinal fin with rectangular profile.

Case 1: solid fin with temperature dependent internal heat generation and constant thermal conductivity

Consider a longitudinal fin with a constant rectangular profile, section area A, length L, perimeter P, thermal conductivity k, and heat generation . Fin is attached to a surface with constant temperature and losses heat to the surrounding medium with temperature through a constant convective heat transfer coefficient h. Here we assumed that the temperature variation in the transfer direction is negligible, so heat conduction occurs only in the longitudinal direction (x direction). A schematic diagram of the described fin is shown in Fig. 1.

Fig. 1

Schematic of fin geometry with the internal heat generation source

The governing differential equation and boundary condition for this problem can be written as Ghasemi et al. (2014)Here it is assumed that the temperature heat generation in the solid fin varies with temperature ,that can be defined asWhere is the internal heat generation at temperature . On introducing the dimensionless variablesThe dimensionless form of Eqs. (16)–(18) can be written as

To understand the anomalous behavior of this system, we fractionalize the energy balance Eq. (21) into fractional order () in order to find the fin temperature in solid fins as,With boundary conditions

Now, applying sumudu transform on both sides of Eq. (24), we obtaintaking inverse sumudu transform on both side of Eq. (26), we get

Now on applying Adomian decomposition method it obtains,

Comparing the coefficients of in Eq. (28), we haveSumming these terms, the final temperature field , is calculated up to five terms as followsEquation (29) represents the expression for finding the temperature distribution in solid fin of fractional order energy balance Eq. (24), where the coefficient K denotes the fin tip temperature, and it can be determined at the boundary condition .

Case 2: porous fin with the temperature dependent internal heat generation

Here the energy balance equation for rectangular porous profile fin with temperature dependent internal heat generation can be written asThe mass flow rate of the fluid passing through the porous material can be written asFrom Darcy’s model, the passage velocity (Hatami et al. 2013) isOn substituting Eqs. (31) and (32) into Eq. (30), it givesAs , Eq. (33) becomesAlso, from Fourier’s law of conductionwhere is the effective thermal conductivity of the porous fin, that can be obtained from the following equation (Hatami et al. 2013)Where is the porosity of the porous fin. On substituting Eq. (35) into Eq. (34), it givesOn introducing the dimensionless variables and numbersThe dimensionless form of Eq. (37) can be written aswith its boundary condition where M is a convection parameter that indicates the effect of surface convecting of the fin and is a porous parameter that indicates the effect of the permeability of the porous medium as well as the buoyancy effect so higher the value of indicates higher permeability of the porous medium or higher buoyancy forces.

To understand the anomalous behaviour of this system, we fractionalize the energy balance Eq. (39) into fractional order () in order to find the fin temperature in porous fins as,with boundary conditionsNow, again applying sumudu transform on both sides of Eq. (41), it obtainsApplying Inverse sumudu transfer on both sides, we getby applying Adomian Decomposition Method, it obtains the following equationwhere is the nonlinear term which can be determined by using Eq. (13). The first few components of the adomian’s polynomial for corrosponding nonlinear terms are given byOn comparing the coefficients of like powers of in Eq. (45), we get decomposition components asBy solving above equations, we getThe approximate solution of up to four terms becomeswhere value of K can be determined at the boundary condition using Eq. (47). Since a constant K is assumed as an initial guess, it automatically satisfies the given boundary condition.

Stability analysis via fixed point theorem

Let be a Banach space and H a self-map of X. Let be particular recursive procedure. Suppose that F(H) the fixed point set of H has at lease one element and that converges to a point . Let and define . If , then the iteration is said to be H-stable. Without any loss of generality, we must assume that, our sequence has an upper boundary; otherwise we cannot expect the possibility of convergence. If all these conditions are satisfied for which is known as Picard’s iteration, then the iteration will be H-Stable. Now we state the following theorem (Atangana 2016).

Theorem 1

Letbe a Banach space andHa self-map ofXsatisfyingfor allx, yinXwhereThenHis Picard H-Stable.

Let the following succession correlate to the nonlinear fractional order energy balance equation (Eq. 46),

Theorem

LetTbe a self-map defined asis T-stable inif

Proof

The first step of the proof will consist on showing that T has a fixed point. To achieve this, we evaluate the following for all Using the linearity property of the inverse sumudu transform, we obtainUsing the triangle inequality for the norm, we getThe above can be further be transformed using the property of norm and integral as followsThe evaluation of Eq. (53) can be done as followsandwhere Now putting together Eqs. (54) and (55) into Eq. (53), we obtain the followingwith Hence, the nonlinear T-self mapping has a fixed point. This completes the proof. Further we show that, T satisfies the condition in Theorem 1, Now forshows that conditions of this theorem holds for the nonlinear mapping T. Since all condition in Theorem 1 hold for the defined non-linear mapping T. Hence, T is Picard’s T-stable. This completes the proof of this theorem.

Results and discussion

The purpose of this work is to observe the simultaneous effects of the governing parameters and the different fractional values on solid and porous fins. The range of thermal and physical parameters selected the present work is shown in Table 1.

Table 1

Range of values for physical and thermal parameters

Parameters	M	G	I_G	ξ
Values	0.5, 1 and 5	0.1–0.9	0.1–0.9	0.1–0.9

Solid fin with temperature dependent internal heat generation and constant thermal conductivity

Temperature distribution for this case (temperature dependent heat generation and constant thermal conductivity) is shown in Figs. 2 and 3 where M = 1 that is common in fin design.

Fig. 2

Temperature variation obtained by ADSTM in comparison with the numerical solution for and M = 1

Fig. 3

Temperature distribution in solid fins for M = 1 and a , b , c

Figure 2 shows the temperature distribution for this state and , and , G = 0.6. This choice of parameters represents a fin with moderate temperature dependent heat generation and the thermal conductivity variation of 20 % between the base and the surrounding coolant temperatures that are often used in nuclear rods.

It is shown in Figs. 2 and 3 that temperature of the fin increases by increasing the value of and G because of increasing in heat generation. The comparison of obtained results with numerical results reveals that ADSTM has good efficiency and accuracy.

Figure 3 shows the temperature distribution for this state and and , G = 0.6 and for the different fractional order value of  = 1.75, 1.5 and 1.25. Further, the nature of the graphs depicts that, considered value of represents the point of convergence under the given range of interval between 1 and 2.

Porous fin with temperature dependent internal heat generation

Figure 4 shows the comparison between the numerical solution and the ADSTM solution for temperature distribution when M = 1, G = 0.4, and for the different value of and . When the porous parameter increases, it can be noticed from Fig. 5 that, a declines in fin temperature which causes stronger cooling results decrease in temperature distribution. Figure 5 shows the variation of temperature distribution of porous fin with temperature dependent internal heat generation when M = 1, G = 0.4, and for the different value of and for  = 1.75, 1.5, 1.25.

Fig. 4

Comparison of dimensionless temperature distribution in porous fins with temperature dependent international heat generation for the different values of porous parameter and

Fig. 5

Temperature distribution of porous fin with temperature dependent internal heat generation for the different values of porous parameter and for a , b , c

The variation of the temperature distribution along the fin for the different values of internal heat generation, when M = 1, G = 0.4, and for is illustrated in Fig. 6. It is observed from figure that, if the internal heat generation increases then fin temperature be increases. Figure 7 shows the variation of temperature distribution of porous fin with temperature dependent internal heat generation when M = 1, G = 0.4, and for the different value of and for .

Fig. 6

Comparison of dimensionless temperature distribution in porous fins with temperature dependent international heat generation for the different values of internal heat generation parameter and

Fig. 7

Temperature distribution of porous fin with temperature dependent internal heat generation for the different values of internal heat generation parameter and for a , b , c

Figure 8 shows the comparison between the numerical solution and the ADSTM solution for temperature distribution when M = 1, and for the different value of G and for .

Fig. 8

Comparison of dimensionless temperature distribution in porous fins with temperature dependent international heat generation for the different values of heat generation parameter and

Figure 9 shows the variation of temperature distribution of porous fin with temperature dependent internal heat generation when M = 1, and for the different value of G and for .

Fig. 9

Temperature distribution of porous fin with temperature dependent internal heat generation for the different values of heat generation parameter and for a , b , c

Figure 10 shows the comparison between the numerical solution and ADSTM solution for temperature distribution when G = 0.4, and for the different value of M and for . At last, it can be concluded that the analytical results correspond exactly with the numerical results. This means that the ADSTM has a high aptitude in solving highly nonlinear initial and boundary value problems without involving linearization. Figure 11 shows the variation of temperature distribution of porous fin with temperature dependent internal heat generation when G = 0.4, and for the different value of M and for .

Fig. 10

Comparison of dimensionless temperature distribution in porous fins with temperature dependent international heat generation for the different values of convective parameter and

Fig. 11

Temperature distribution of porous fin with temperature dependent internal heat generation for the different values of convective parameter and for a , b , c

Conclusion

In this study, the heat transfer in rectangular solid and porous fin with the temperature dependent internal heat generation is analyzed by using ADSTM and used the concept of T-stable mapping and the fixed point theorem to prove the stability of ADSTM. Here, it is shown that, ADSTM provide a simple, accurate and appropriate technique for simulating the heat transfer in solid and porous fin in a fractional order energy balance equation. The results shows that the temperature distribution strongly depends on different parameter in solid fin as well as on Darcy’s number in porous fin and also on the fractional parameter.