Assessment of Homotopy Perturbation and Variational Iteration Methods in Heat Transfer Equations

Assessment of Homotopy Perturbation and Variational Iteration Methods in Heat Transfer Equations

DOI: 10.4018/978-1-5225-2713-8.ch004


We explained VIM and HPM methods in Chapter 1 and Chapter 2. In this chapter, we compare these methods for nonlinear heat transfer problems.
Chapter Preview

The Application Of Vim, Hpm, And Perturbation Method In Heat Transfer

In order to assess the accuracy of VIM for solving nonlinear equations and to compare it with HPM and perturbation method, we will consider the following examples.

Cooling of a Lumped System by Combined Convection and Radiation

Consider the problem of combined convective–radiative cooling of lumped system (Aziz & Na, 1984). Let the system have volume , surface area , density , specific heat emissivity and the initial temperature . At , the system is exposed to an environment with convective heat transfer with the coefficient of and the temperature . The system also loses heat through radiation and the effective sink temperature is . The cooling equation and the initial conditions are as follows:


To solve the equation, we do the following changes of parameters:


After the parameter change, the heat transfer equation will result in the following form:


For simplicity, we assume the case . So we have:


Eqs. (6) and (7) are solved though PM, HPM and VIM as follows:

Variational Iteration Method

To solve Eqs. (6) and (7), we can construct the following correction functional:


Its stationary conditions can be obtained as follows:


The Lagrangian multiplier can therefore be identified as:


As a result, we obtain the following iteration formula:


Now we start with an arbitrary initial approximation that satisfies the initial condition:


Using the above variational formula (11), we have:


Substituting Eq. (12) into Eq. (13) and after some simplifications, we have:


In the same way, we obtain as follows:

(15) and so on. In the same way the rest of the components of the iteration formula can be obtained.

Homotopy Perturbation Method

After separating the linear and nonlinear parts of the equation, we apply HPM to Eqs. (6) and (7), as follows:

16) where is the linear part of the equation and is the initial approximation (Rajabi, Ganji & Taherian, 2006).We consider as:


Substituting Eq. (17) into Eq. (16) and rearranging based on powers of terms, we have:

(19) and
(21) and finally,


Solving Eqs. (18)- (23) results in When , we have:


Perturbation Method

For very small , let us assume a regular perturbation expansion and calculate the first three terms (Nayfeh, 1973), thus we assume:


Substituting Eq. (25) into Eq. (6) and after expansion and rearranging based on coefficients of terms, we have:


The solutions of Eqs. (26)- (28) are:


The three-term expansion in Eq. (25) now becomes:


Complete Chapter List

Search this Book: