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.

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 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 . 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 h . 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.  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  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).  applied the fractional Adams-Bashforth-Moulton Method for obtaining the numerical solutions of some linear and nonlinear fractional ordinary differential equations.  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.  studied solitary wave solutions for vakhnenko-parkes equation via exp-function and Exp (−φ(ξ ))-expansion method. Also they  studied traveling wave solutions of nonlinear partial differential equation via new extended (G ′ /G)-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.

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 α > 0 (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 Ross where R + is the set of positive real numbers and Ŵ(·) is the gamma function.

Definition
The fractional derivative of f(t) in the Caputo sense is defined by where 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.

Definition
The sumudu transform of a function f(t), defined for all real numbers t ≥ 0, is the function F(u), defined by Watugala (1993) Definition Sumudu transform of function derivatives is defined as Belgacem and Karaballi (2006) Definition Sumudu transform of the Caputo fractional derivative is defined as Belgacem and Karaballi (2006) Theorem Let G(u) be the sumudu transform of f(t) such that

function f(t) is given by
For the proof see Belgacem and Karaballi (2006).

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 condition where d α dζ α denotes without loss of generality the Caputo fraction derivative operator, R is the linear differential operator, N represents the nonlinear differential operator, and g(ζ ) is the source term.
On applying the sumudu transform and the Caputo fractional derivative in Eq. (7), it obtains Sumudu inverse transform of Eq. (9), gives where θ(0) = K is prescribed initial conditions. Using Adomian decomposition method which obtains the approximate solution of Eq. (10) in a series form as Adomian and the nonlinear term can be expressed as a sum of Adomian Polynomials where the Adomian polynomials A n (θ), depends upon the solution components θ 0 , θ 1 , θ 2 , . . . , θ n , can be defined as follows On substituting Eqs. (11) and (12) in Eq. (10), it obtains The 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 as

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 Q . Fin is attached to a surface with constant temperature T b and losses heat to the surrounding medium with temperature T ∞ 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.
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 T ∞ ,that can be defined as Where Q ∞ is the internal heat generation at temperature T ∞ . On introducing the dimensionless variables The 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 (α > 0) 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 obtain Summing these terms, the final temperature field θ (ζ ), is calculated up to five terms as follows Equation (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 θ | ζ =1 = 1.

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 as The mass flow rate of the fluid passing through the porous material can be written as From Darcy's model, the passage velocity ) is On substituting Eqs. (31) and (32) 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 (α > 0) in order to find the fin temperature in porous fins as, with boundary conditions Now, again applying sumudu transform on both sides of Eq. (41), it obtains Applying Inverse sumudu transfer on both sides, we get by applying Adomian Decomposition Method, it obtains the following equation where A n (θ) 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 by On comparing the coefficients of like powers of in Eq. (45), we get decomposition components as

By solving above equations, we get
The approximate solution of θ (ζ ) up to four terms becomes where value of K can be determined at the boundary condition θ | ζ =1 = 1 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 (X, �·�) be a Banach space and H a self-map of X. Let θ n+1 = f (H , θ n ) be particular recursive procedure. Suppose that F(H) the fixed point set of H has at lease one element and that θ n converges to a point p ∈ F (H). Let θ n ⊆ X and define e n = θ n+1 − f (H, θ n ) . If lim n→∞ θ n = p, then the iteration θ n+1 = f (H, θ n ) is said to be H-stable. Without any loss of generality, we must assume that, our sequence θ n has an upper boundary; otherwise we cannot expect the possibility of convergence. If all these conditions are satisfied for θ n+1 = H θ n which is known as Picard's iteration, then the iteration will be H-Stable. Now we state the following theorem (Atangana 2016 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 (n, k) ∈ N × N Using the linearity property of the inverse sumudu transform, we obtain Using the triangle inequality for the norm, we get The above can be further be transformed using the property of norm and integral as follows The evaluation of Eq. (53) can be done as follows �T (θ n (ζ )) − T (θ k (ζ ))� ≤ S −1 u α S M 2 − M 2 GI G (θ n (ζ ) − θ k (ζ )) and where ξ β 2 . Now putting together Eqs. (54) and (55) into Eq. (53), we obtain the following with {β 1 + β 2 κ} < 1, ∀β 1 β 2 . 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 for shows 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.

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. Figure 2 shows the temperature distribution for this state and I G = G = 0.2 , I G = G = 0.4 and I G = 0.4, 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 I G and G because of increasing in heat generation. The comparison of obtained results with numerical results reveals that ADSTM has good efficiency and accuracy.
(54) Values 0.5, 1 and 5 0.1-0.9 0.1-0.9 0.1-0.9 Figure 3 shows the temperature distribution for this state and I G = G = 0.2, I G = G = 0.4 and I G = 0.4, 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. Figure 4 shows the comparison between the numerical solution and the ADSTM solution for temperature distribution when M = 1, G = 0.4, I G = 0.6 and for the different value of ξ and α = 2. 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, I G = 0.6 and for the different value of ξ and for α = 1.75, 1.5, 1.25. The variation of the temperature distribution along the fin for the different values of internal heat generation, I G when M = 1, G = 0.4, ζ = 0.4 and for α = 2 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, ξ = 0.4 and for the different value of I G and for α = 1.75, 1.5, 1.25. Figure 8 shows the comparison between the numerical solution and the ADSTM solution for temperature distribution when M = 1, I G = 0.4, ξ = 0.4 and for the different value of G and for α = 2. Figure 9 shows the variation of temperature distribution of porous fin with temperature dependent internal heat generation when M = 1, I G = 0.4, ξ = 0.4 and for the different value of G and for α = 1.75, 1.5, 1.25. Figure 10 shows the comparison between the numerical solution and ADSTM solution for temperature distribution when G = 0.4, I G = 0.6, ξ = 0.4 and for the different value of M and for α = 2. 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, I G = 0.6, ξ = 0.4 and for the different value of M and for α = 1.75, 1.5, 1.25.