Open access peer-reviewed chapter

Application of Perturbation Theory in Heat Flow Analysis

Written By

Neelam Gupta and Neel Kanth

Submitted: 07 August 2020 Reviewed: 21 December 2020 Published: 12 January 2021

DOI: 10.5772/intechopen.95573

From the Edited Volume

A Collection of Papers on Chaos Theory and Its Applications

Edited by Paul Bracken and Dimo I. Uzunov

Chapter metrics overview

585 Chapter Downloads

View Full Metrics

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.

Advertisement

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
Advertisement

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 A is a general differential operator, linear or nonlinear, fr is a known analytic function, B is a boundary operator and Γ is the boundary of the domain Ω. The operator A can be generally divided into two operators, L and N, where Lis linear and N is 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 p01 is called the homotopy parameter and u0 is 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 p as 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
Advertisement

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 p and 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θθ0n is the taylor series expansion which can be used in Eq. (9).

Advertisement

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 v in 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β2 5L4π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 ƥ=1 in (23) is given by

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

On substituting values of vis in 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 p and 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.

Advertisement

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.

References

  1. 1. N. Gupta and N. Kanth (2019). Study of heat flow in a rod using homotopy analysis method and homotopy perturbation method. AIP Conference Proceedings, 2061(1) 020013 1–8
  2. 2. D. Grover, V. Kumar, and D. Sharma, A Comparative Study of Numerical Techniques and Homotopy Perturbation Method for Solving Parabolic Equation and Nonlinear Equations, Int. J. Comput. Methods Eng. Sci. Mech., vol. 13, no. 6, pp. 403–407, 2012
  3. 3. J.-H. He, Homotopy perturbation technique, Comput. Methods Appl. Mech. Eng., vol. 178, no. 3, pp. 257–262, Aug. 1999
  4. 4. J.-H. He, Comparison of homotopy perturbation method and homotopy analysis method, Appl. Math. Comput., vol. 156, no. 2, pp. 527–539, Sep. 2004
  5. 5. J.-H. He, A coupling method of a homotopy technique and a perturbation technique for non-linear problems, Int. J. Non-Linear Mech., vol. 35, no. 1, pp. 37–43, Jan. 2000
  6. 6. J.-H. He, Some asymptotic methods for strongly nonlinear equations, Int. J. Mod. Phys. B, vol. 20, no. 10, pp. 1141–1199, Apr. 2006
  7. 7. S. Liang and D. Jeffrey, Comparison of Homotopy Analysis Method and Homotopy Perturbation Method through an Evolution Equation, Commun. Nonlinear Sci. Numer. Simul., vol. 14, pp. 4057–4064, Dec. 2009
  8. 8. A. Demir, S. Erman, B. Ozgur, E. Korkmaz, Analysis of the new homotopy perturbation method for linear and nonlinear problems. Boundary Value Problems, 2013(1), 1–11, 2013
  9. 9. J-H. He, Homotopy perturbation method: a new nonlinear analytical technique. Applied Mathematics and Computation, 135(1) 73–79, 2003
  10. 10. J. Biazar, M. Eslami, A new homotopy perturbation method for solving systems of partial differential equations. Computers and Mathematics with Applications, 62(1), 225–234, 2011
  11. 11. J-H. He, Homotopy perturbation method for solving boundary value problems. Physics Letters A. 350(1), 87–88, 2006
  12. 12. M. Mirzazadeh, Z. Ayati, New homotopy perturbation method for system of Burgers equations. Alexandria Engineering Journal. 55(3), 1619–1624, 2016
  13. 13. M. Elbadri, A New Homotopy Perturbation Method for Solving Laplace Equation. Advances in Theoratical and Applied Mathematics. 8(3), 237–242, 2013
  14. 14. M. Elbadri, T.M. Elzaki, New Modification of Homotopy Perturbation Method and the Fourth - Order Parabolic Equations with Variable Coefficients. Pure Appl. Math. J. 4(6) 242–247, 2015
  15. 15. J.H. He, A coupling method of homotopy technique and a perturbation technique for non linear problems, Int. J. Non linear Mech., 35, pp. 37–43, 2000
  16. 16. J.H. He, Application of homotopy perturbation method to nonlinear wave equations, Chaos Solitons Fractals, 26, pp. 295–700,2005
  17. 17. A. Yıldırım, Solution of BVPs for fourth-order integro-differential equations by using homotopy perturbation method, Comput. Math. Appl., 56 (12), 3175–3180, 2008
  18. 18. S. Abbasbandy, Homotopy perturbation method for quadratic Riccati differential equation and comparison with Adomian’s decomposition method, Appl. Math. Comput., 172, pp. 485–490, 2006
  19. 19. S.J. Liao, An approximate solution technique not depending on small parameters: A special example, Int. J. Non linear Mech., 36, 371–180, 1995
  20. 20. N. Gupta and N. Kanth, Study of heat conduction inside rolling calender nip for different roll temperatures, Journal of physics: Conference series, 1276(1), 012044 1–9, 2019
  21. 21. N. Gupta and N. Kanth, Analytical Approximate Solution of Heat Conduction Equation Using New Homotopy Perturbation Method, Matrix Science Mathematic (MSMK), 3(2), 01–07, 2019
  22. 22. N. Gupta and N. Kanth, Numerical Solution of Diffusion Equation Using Method of Lines, Indian Journal of Industrial and Applied Mathematics, 10(2), 194–203, 2019
  23. 23. N. Gupta and N. Kanth, Analysis of heat conduction inside the calender nip used in textile industry, AIP Conference Proceedings, 2214(1), 020008, 2020

Written By

Neelam Gupta and Neel Kanth

Submitted: 07 August 2020 Reviewed: 21 December 2020 Published: 12 January 2021