Open access peer-reviewed chapter - ONLINE FIRST

Application of Perturbation Theory in Heat Flow Analysis

By Neelam Gupta and Neel Kanth

Submitted: August 7th 2020Reviewed: December 21st 2020Published: January 12th 2021

DOI: 10.5772/intechopen.95573

Downloaded: 21

Abstract

Many physical and engineering problems can be modeled using partial differential equations such as heat transfer through conduction process in steady and unsteady state. Perturbation methods are analytical approximation method to understand physical phenomena which depends on perturbation quantity. Homotopy perturbation method (HPM) was proposed by Ji Huan He. HPM is considered as effective method in solving partial differential equations. The solution obtained by HPM converges to exact solution, which are in the form of an infinite function series. Biazar and Eslami proposed new homotopy perturbation method (NHPM) in which construction of an appropriate homotopy equation and selection of appropriate initial approximation guess are two important steps. In present work, heat flow analysis has been done on a rod of length L and diffusivity α using HPM and NHPM. The solution obtained using different perturbation methods are compared with the solution obtained from most common analytical method separation of variables.

Keywords

  • heat conduction equation
  • homotopy perturbation method
  • new homotopy perturbation method
  • specific heat
  • diffusivity

1. Introduction

Partial differential equations play a dominant role in applied mathematics. The classical heat conduction equation is second order linear partial differential equation. The solutions of which are obtained by using various analytical and numerical methods [1, 2, 3]. This equation describes the heat distribution in each domain over some time. Jean-Joseph Fourier was the first to formulate and describe the heat conduction process [1, 4]. Perturbation methods depending upon small/large parameters have been encountered from past few years. Perturbation methods are analytical approximation method to understand physical phenomena which depends on perturbation quantity. But these methods do not provide an easy way to find out the rapid convergence of approximate series. Therefore, this method is simple, suitable and appropriate method to provide the rapid convergence of series [5, 6, 7]. The perturbation method along with the homotopy method has been employed to develop a hybrid method known as homotopy perturbation method (HPM) [1, 2, 3, 4]. Ji-Huan was the first to introduce HPM. Homotopy perturbation method provides analytical approximation to linear/nonlinear problems without linearization or discretization. It helps in formulating simpler equations by breaking down the complex problems, which can be solved easily. Since HPM does not depend on small parameters, therefore drawbacks of the existing perturbation methods can be abolished [8, 9, 10, 11]. The solution obtained by HPM converges to exact solution, which are in the form of an infinite function series. Various problems are modeled by linear and non-linear partial differential equations problems in the fields of physics, engineering etc. To solve such kind of partial differential equations (PDE), many methods are used to find the numerical or exact solutions. Homotopy perturbation method (HPM) is one of the methods used in recent years to solve various linear and non-linear PDE [12, 13, 14, 15]. Initial and boundary value problems can be solved using HPM extensively. Many researchers and scientists show great interest in homotopy perturbation method. Huan was the first who described homotopy perturbation method. He showed that this method is a one of the powerful tools used to investigate various problems which are arising nowadays. HPM is used for solving linear and non-linear ordinary and partial differential equations [16].

In HPM, complex linear or non-linear problem can be continuously distorted into simpler ones. Perturbation theory and homotopy theory in topology is combined to develop homotopy perturbation method [1]. HPM is applicable to linear and non-linear boundary value problems. The solution obtained by HPM gives the solution approximately near to the universally accepted method of separation of variable [17, 18, 19].

Recently, Biazar and Eslami proposed the new homotopy perturbation method (NHPM). Construction of an appropriate homotopy equation and selection of appropriate initial approximation guess are two important steps of NHPM [19, 20]. The study reveals that with less computational work, we can construct proper homotopy by decomposition of source function in a correct way. New homotopy perturbation method is the most powerful tool which can be used to obtain analytical solution of various kinds of linear and nonlinear PDE’s. This method is widely used by researchers to obtain solution of various functional Equations [20, 21, 22].

To develop this new technique, HPM is combined with the decomposition of source function. The decomposition of a source function is the basis of homotopy used in this method because convergence of a solution is affected by the decomposition of source functions [23]. Different kind of homotopy can be formed using various decomposition of a source functions. This study is aimed at constructing suitable homotopy by decomposition of a source function which requires less computational efforts and made calculations in simpler form unlike other perturbation methods. The obtained results directly imply the fact that NHPM is very influential as compared to HPM or any other perturbation technique. To establish exact solution of linear and non-linear problem with boundary and initial condition, new homotopy method is most appropriate method to apply [23].

The two most important steps in application of new homotopy perturbation method to construct a suitable homotopy equation and choose a suitable initial guess, we aim in this work to effectively employ the (NHPM) to establish exact solution for two-dimensional Laplace equation with Dirichlet and Neumann boundary condition, the difference between (NHPM) and standard (HPM) is starts from the form of initial approximation of the solution.

In this chapter, the semi analytic solution of one-dimensional heat conduction equation is obtained by means of homotopy perturbation method and new homotopy perturbation method. These methods are effectively applied to obtain the exact solution for the problem in hand which reveals the effectiveness and simplicity of the method. Numerical results have also been analyzed graphically to show the rapid convergence of infinite series expansion. The obtained analytic solution for one dimensional heat conduction equation with boundary and initial conditions using NHPM is same as the universally accepted exact solution. This tells us about the capability and reliability of this method. The solution obtained using NHPM is considered in the form of an infinite series. The convergence of solution to the exact solution is very rapid.

2. Heat conduction equation

The one-dimensional heat equation

Uθ=β2Uz2E1

with boundary conditions

U0θ=0,U1θ=0,E2

and initial condition

Uz0=hz,0z1.E3

3. Basic idea of Homotopy perturbation method

First, we outline the general procedure of the homotopy perturbation method developed and advanced by He. We consider the differential Eq. [2]

Aufr=0,rΩE4
Bu∂u∂x=0,rΓE5

where Ais a general differential operator, linear or nonlinear, fris a known analytic function, Bis a boundary operator and Γis the boundary of the domain Ω. The operator Acan be generally divided into two operators, Land N, where Lis linear and Nis a nonlinear operator. Eq. (4) can be written as

Lu+Nufr=0E6

Using the homotopy technique, we can construct a homotopy [1,2]

vrp:Ω×01Rwhich satisfies the relation
Hvp=1pLvLu0+pAvfr=0,rΩE7

Here p01is called the homotopy parameter and u0is an initial approximation for the solution of Eq. (4), which satisfies the boundary conditions. Clearly, from Eq. (7), we have

Hv0=LvLu0E8
Hv1=AvfrE9

We assume that the solution of Eq. (7) can be expressed as a series in pas follows:

v=v0+pv0+p2v2+p3v3+E10

On setting p=1, we obtain the approximate solution of Eq. (10) as

u=limp1v=v0+v0+v2+v3+E11

4. Basic idea of new Homotopy perturbation method

First, following homotopy is constructed for solving heat conduction equation using NHPM

1pTθU0+pTθβ2Tz2=0E12

Taking L1=θ0θ.i.e. inverse operator on Eq. (12), then

Tzθ=θ0θU0zθpθ0θU0β2Tz2+Tzθ0E13
WhereTzθ0=Uzθ0

Let the solution of Eq. (13) is given by

T=T0+pT1+p2T2+p3T3+E14

where T0,T1,T2,T3,are to be determined.

Suppose solution given by Eq. (14) is the solution of Eq. (13). On comparing the coefficients of powers of pand equating to zero and using Eq. (14) in Eq. (13), following are obtained:

p0:T0zθ=θ0θU0zθ+Tzθ0
p1:T1zθ=θ0θU0zθβ2T0z2
p2:T2zθ=θ0θβ2T1z2
p3:T3zθ=θ0θβ2T2z2
andsoonE15

Consider the initial approximation of Eq. (1) as

U0zθ=n=0cnzPnθ,Tz0=Uz0,Pkθ=θk,E16

where, P1θ,P2θ,P3θ,and c0z,c1z,c2z,are specified functions and unknown coefficients respectively, depending on the problem.

Using Eq. (16) in (15), following are obtained:

T0zθ=c0zθ+c1zθ22+c2zθ33+c3zθ44++Uz0
T1zθ=c0zβπ2sinπzθ+12c1z+12βc0zθ2+13c2z+13c1zθ3+
andsoonE17

Now solving the above equations in such a manner that, T1zθ=0.

Therefore Eq. (17) reduces to

T1zθ=T2zθ==0.

So Uzθ=T0zθ=n=0cnzPnθis obtained solution which is found to be exactly same as the exact solution obtained through method of separation of variable.

If U0zθis analytic at θ=θ0,

U0zθ=n=0cnzθθ0nis the taylor series expansion which can be used in Eq. (9).

5. Applications of Homotopy perturbation method and new Homotopy perturbation method

For understanding the application of HPM and NHPM, we will solve the one-dimensional heat equation given by

Uθ=β2Uz2E18

with boundary conditions

U0θ=0,U1θ=0,E19

and initial condition

Uz0=sin2πzL,0zL.E20

The homotopy for the diffusion equation given by (18) is obtained as follows [2].

vθu0θ+ƥu0θβ2vz2=0E21

Let u0=sin2πzLcosπ2θbe the initial approximation, which satisfies boundary conditions given by (19).

Let solution of (18) has the following form

v=v0+ƥv1+ƥ2v2+ƥ3v3+ƥ4v4+E22

On substituting the value of vin Eq. (21) and comparing the coefficients of like powers of ƥ we obtain

ƥ0:v0θ=u0θ
ƥ1:v1θ=β2v0z2,v10θ=0=v1Lθ
ƥ2:v2θ=β2v1z2,v20θ=0=v2Lθ
ƥ3:v3θ=β2v2z2,v30θ=0=v3Lθ
ƥn:vnθ=β2vn1z2,vn0θ=0=vnLθE23

On solving the system of Eq. (23) using Mathematica 5.2

v0=u0=sin2πzLcosπ2θ
v1θ=4βπ2L2v0v1=βsin2πzLsinπ2θL2+sinπzLv2θ=4βπ2L2v1v2=βL2π2θ+αcosπ2θsin2πzLL4+L4sin2πzL+β2sin2πzLL4v3θ=4βπ2L2v2v3=βπ2θ2L4+L2π2βθ2β2+2β2sinπ2θsin2πzL2L6+sin2πzL
v4θ=4βπ2L2v3v4=16L8(β(π2θ(6L6+3L4π2βθL2π4θ2β2+3β3π2θ)6β3cos[π2θ])sin[2πzL])+L8sin[2πzL]β4sin[2πzL]L8
v5θ=4βπ2L2v4v5=124L10(β(π2θ(24L812L6π2βθ+4L4π4θ2β2L2π6θ3β3+4(6+π4θ2)β4)+24β4sin[π2θ])sin[2πzL])+sin[2πzL]
v6θ=βπ2L2v5v6=1120L12(β(π2θ(120L1060L8π2βθ+20L6π4θ2β25L4π6θ3β3+L2π8θ4β45π2θ(12+π4θ2)β5)+120β5cos[π2θ])sin[2πzL])+L12sin[2πzL]+β6sin[2πzL]L12
andsoonE24

The approximate solution of (1) by setting ƥ=1in (23) is given by

u=limp1v=v0+v1+v2+v3+v3+E25

On substituting values of visin Eq. (25), solution is obtained in terms of a summation of infinite series which gives results near to the exact solution.

Now we will solve the Eq. (18) using NHPM. First of all, following homotopy is constructed for solving heat conduction equation using NHPM

1pTθU0+pTθβ2Tz2=0E26

Taking L1=θ0θ.i.e. inverse operator on Eq. (26), then

Tzθ=0θU0zθp0θU0β2Tz2+Tz0.E27

Let the solution of the (27) is

T=T0+pT1+p2T2+p3T3+,E28

where, T0,T1,T2,are to be determined.

Suppose Eq. (25) is the solution of Eq. (24). Comparing the coefficients of powers of pand equating to zero and using Eq. (25) in Eq. (24), following are obtained:

p0:T0zθ=0θU0zθ+Tz0p1:T1zθ=0θU0zθβ2T0z2p2:T2zθ=0θβ2T1z2p3:T3zθ=0θβ2T2z2
andsoon.E29

Consider initial approximation of Eq. (18) as

U0zθ=n=0cnzPnθ,Tz0=Uz0,Pkθ=θk,E30

where, P1θ,P2θ,P3θ,and c0z,c1z,c2z,are specified functions and unknown coefficients respectively, depending on the problem.

Using Eq. (30) in (29), following are obtained:

T0zθ=c0zθ+c1zθ22+c2zθ33+c3zθ44++sin2πzL
T1zθ=c0z4βπ2Lsin2πzLθ+12c1z+12βc0zθ2+13c2z+13c1zθ3+
andsoonE31

Now solving the above equations in such a manner that, T1zθ=0.

Therefore Eq. (31) reduces to

c0z=22βπ2Lsin2πzL
c1z=24β2π4L2sin2πzL
c2z=26β3π6L3sin2πzL
U(z,θ)=T0(z,θ)=sin2πzL+c0(z)θ+c1(z)θ22+c2(z)θ33+c3(z)θ44+=sin2πzL[122βπ2Lθ+24β2π4L2θ2226β3π6L3θ33+]=sin2πzLe4βπ2LθE32

which is same as the universally accepted exact solution for the problem which is shown in Figure 1.

Figure 1.

Solution using NHPM.

The solution of one-dimensional heat conduction equation is solved using HPM and NHPM and then compared with the universally accepted exact solution obtained from method of separation of variable. Figure 2 represents the comparison of solution of heat equation using HPM, NHPM and method of separation of variable. It is found that the solution obtained using HPM gives result near to the exact solution whereas solution using NHPM gives same results as the exact solution.

Figure 2.

Comparison of HPM, NHPM and the exact solution.

6. Conclusion

The analytical approximate solutions of one-dimensional heat conduction equation are obtained by applying new homotopy perturbation method and new homotopy perturbation method. It is found that new homotopy perturbation method (NHPM) converges very rapidly as compared to homotopy perturbation method (HPM) and other traditional methods. The exact solutions are obtained up to more accuracy using NHPM. An infinite convergent series solution for particular initial conditions are obtained using these methods which shows the effectiveness and efficiency of NHPM and HPM. The convergence rate of NHPM is much faster than traditional methods which directly indicates that this method is better than other methods. The solution of heat equation obtained by homotopy perturbation method and new homotopy perturbation method are exactly same and very close to the solution obtained by universally accepted and tested analytical method of separation of variables. If the initial guess in homotopy perturbation method is effective and properly chosen which satisfy boundary and initial condition, homotopy perturbation method provides solution with rapid convergence. It is illustrated that NHPM is very prominent, when accuracy has a vital role to play. The numerical results also reflect the remarkable applicability of NHPM to linear and non-linear initial and boundary value problems. NHPM provides the rapid convergence of the series solution for linear as well as non-linear problems with less computational work.

Download for free

chapter PDF

© 2021 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Neelam Gupta and Neel Kanth (January 12th 2021). Application of Perturbation Theory in Heat Flow Analysis [Online First], IntechOpen, DOI: 10.5772/intechopen.95573. Available from:

chapter statistics

21total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us