Research Article

# Variational Iteration Method and Homotopy-Perturbation Method for Solving Different Types of Wave Equations

A. Barari, Abdoul R. Ghotbi, F. Farrokhzad and D.D. Ganji

ABSTRACT

Due to wide range of interest in use of wave equations to gain insight in to vibrating problems, Homotopy Perturbation Method (HPM) and Variational Iteration Method (VIM) are employed to approximate the solution of the three types of wave equations including one-dimensional wave equation, kinematic wave equation and non-linear homogeneous wave equation. The final results obtained by HPM and VIM are compared with those results obtained from the exact solution. The comparison shows a precise agreement between the results and introduces these new methods as the applicable methods which they need less computations and are much easier and more convenient than other approximate methods, so they can be widely used in engineering.

 Services Related Articles in ASCI Search in Google Scholar View Citation Report Citation

 How to cite this article: A. Barari, Abdoul R. Ghotbi, F. Farrokhzad and D.D. Ganji, 2008. Variational Iteration Method and Homotopy-Perturbation Method for Solving Different Types of Wave Equations. Journal of Applied Sciences, 8: 120-126. DOI: 10.3923/jas.2008.120.126 URL: https://scialert.net/abstract/?doi=jas.2008.120.126

INTRODUCTION

In this study, we consider the linear and non-linear wave equation (Inc et al., 2004):

 (1)

On the finite x-interval [0, π] with Dirichlet boundary condition:

 (2)

and F (x,t) is a given function of x and t.

This equation describes the propagation of a wave (or disturbance) and it arises in a wide variety of physical problems. some of these problems include a vibrating string, vibrating membrane, longitudinal vibrations of an elastic rod or beam, shallow water waves, acoustic problems for the velocity potential for a fluid flow through which sound can be transmitted, transmission of electric signals along a cable, shock waves, chemical exchange processes in chromatography, sediment transport in rivers and waves in plasmas and both electric and magnetic fields in the absence of charge and dielectric (Debnath, 1997).

We will apply the Variational Iteration Method (VIM) (He, 1999a, 2000; Momani and Abuasad, 2006; Ganji et al., 2007; Sweilam and Khader, 2007; Bildik and Konuralp, 2006) and Homotopy Perturbation Method (HPM) (He, 1999b, 2006; Zhang and He, 2006; Ganji and Sadighi, 2007) for the three types of partial differential wave equations.

Variational Iteration Method (VIM) and Homotopy-Perturbation Method (HPM) are the most effective and convenient ones for both linear and nonlinear equations. The VIM is to construct correction functional using general Lagrange multipliers identified optimally via the variational theory and the initial approximation can be freely chosen with unknown constants. The HPM deforms a difficult problem in to a simple problem which can be easily solved using ordinary methods.

Extensive studies have been done regarding this matter. Such as investigation of Helmholtz equation and fifth-order KdV equation or Nonlinear Coupled Systems of reaction-diffusion Equations using Homotopy perturbation method (Ganji and Sadighi, 2006).

BASIC IDEA OF HOMOTOPY-PERTURBATION METHOD

Linear and Nonlinear phenomena are of fundamental importance in various fields of science and engineering. Most models of real-life problems are still very difficult to solve. Therefore, approximate analytical solutions such as Homotopy Perturbation Method (HPM) were introduced. This method is the most effective and convenient ones for both linear and nonlinear equations.

Perturbation method is based on assuming a small parameter. The majority of nonlinear problems, especially those having strong nonlinearity, have no small parameters at all and the approximate solutions obtained by the perturbation methods, in most cases, are valid only for small values of the small parameter.

Generally, the perturbation solutions are uniformly valid as long as a scientific system parameter is small. However, we cannot rely fully on the approximations, because there is no criterion on which the small parameter should exists. Thus, it is essential to check the validity of the approximations numerically and/or experimentally. To overcome these difficulties, HPM have been proposed recently.

To explain this method, let us consider the following function:

 (3)

With the boundary conditions of:

 (4)

Where, A,B, f (r) and Γ are a general differential operator, a boundary operator, a known analytical function and the boundary of the domain Ω, respectively.

Generally speaking the operator can be divided in to a linear part L and a nonlinear part N (u). Equation 3 can therefore, be written as:

 (5)

By the homotopy technique, we construct a homotopy v (r, p): Ωx[0,1]→R Which satisfies:

 (6)

Or

 (7)

Where, p ε [0,1] is an embedding parameter, while u0 is an initial approximation of Eq. 3, which satisfies the boundary conditions. Obviously, from Eq. 6 and 7 we will have:

 (8)

 (9)

The changing process of p from zero to unity is just that of v (r, p) from u0 to u (r). In topology, this is called deformation, while L (v)-L (u0) and A (v) -f (r) are called homotopy.

According to the HPM, we can first use the embedding parameter p as a small parameter and assume that the solutions of Eq. 6 and 7 can be written as a power series in p:

 (10)

Setting p = 1 yields in the approximate solution of Eq. 3 to:

 (11)

The combination of the perturbation method and the homotopy method is called the HPM, which eliminates the drawbacks of the traditional perturbation methods while keeping all its advantage.

The series (11) is convergent for most cases. However, the convergent rate depends on the nonlinear operator A (v). Moreover, He (1999b) made the following suggestions:

 • The second derivative of N (v) with respect to v must be small because the parameter may be relatively large, i.e., p→1. • The norm of must be smaller than one so that the series converges.

BASIC IDEA OF VARIATIONAL ITERATION METHOD

To clarify the basic ideas of VIM, we consider the following differential equation:

 (12)

Where:
 L = A linear operator N = A nonlinear operator g (t) = A homogeneous term

According to VIM, we can write down a correction functional as follows:

 (13)

Where, λ is a general lagrangian multiplier which can be identified optimally via the variational theory. The subscript n indicates the nth approximation and un is considered as a restricted variation, i.e., δ " én = 0.

Example 1: We consider a homogeneous linear wave equation:

 (14)

With the initial and boundary conditions posed are:

 (15)

Exact solution of this equation is:

 (16)

APPLICATION OF HOMOTOPY-PERTURBATION METHOD

To solve Eq. 14 by means of HPM, we consider the following process after separating the linear and nonlinear parts of the equation; we apply Homotopy-Perturbation to Eq. 6.

A homotopy-perturbation method can be constructed as follows:

 (17)

Substituting Eq. 10 into 17 and rearranging the resultant equation based on powers of p-terms, one has:

 (18)

 (19)

 (20)

With the following conditions:

 (21)

With the effective initial approximation for v0 from the conditions (21) and solutions of Eq. 18-20 may be written as follows:

 (22)

 (23)

 (24)

In the same manner, the rest of components were obtained using the maple package.

According to the HPM, we can conclude that:

 (25)

Therefore, substituting the values of v0 (x, t), v1 (x, t) and v2 (x, t) from Eq. 22-24 into 25 yields:

 (26)

APPLICATION OF VARIATIONAL ITERATION METHOD

 (27)

Its stationary conditions can be obtained as follows:

 (28)

We obtain the lagrangian multiplier:

 (29)

As a result, we obtain the following iteration formula:

 (30)

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

 (31)

Using the above variational formula (30), we have

 (32)

Substituting Eq. 31 into 32 and after simplifications, we have:

 (33)

In the same way, we obtain u2 (x, t) as follows:

 (34)

And so on. In the same way the rest of the components of the iteration formula can be obtained.

As it can be seen, using these two approximate methods lead to the exact solution (Fig. 1) using more iteration.

Example 2: We consider the first-order non-linear wave equation. The equation of the form:

 (35)

With the exact solution of:

 (36)

APPLICATION OF HOMOTOPY-PERTURBATION METHOD

A homotopy-perturbation method can be constructed as follows:

 (37)

 Fig. 1: 3D obtained results of u(x,t) by the HPM, VIM and the exact solution

Substituting Eq. 10 into 37 and rearranging the resultant equation based on powers of p-terms, one has:

 (38)

 (39)

 (40)

 (41)

With the following conditions:

 (42)

With the effective initial approximation for from the conditions (42) and solutions of Eq. 38-41 may be written as follows:

 (43)

 (44)

 (45)

 (46)

In the same manner, the rest of components were obtained using the maple package.

According to the HPM, we can conclude that:

 (47)

Therefore, substituting the values of v0 (x, t), v1 (x, t), v2 (x, t) and v3 (x, t) from Eq. 43-46 into 47 yields:

 (48)

APPLICATION OF VARIATIONAL ITERATION METHOD

 (49)

Its stationary conditions can be obtained as follows:

 (50)

The lagrangian multiplier can therefore be identified as:

 (51)

As a result, we obtain the following iteration formula:

 (52)

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

 (53)

Using the above variational formula (52), we have

 (54)

Substituting Eq. 53 into 54 and after simplification, we have:

 (55)

 Fig. 2: 3D obtained results of u(x,t) by the HPM, VIM and the exact solution

In the same way, we obtain u2 (x, t), u3 (x, t) as follows:

 (56)

 (57)

And so on. In the same way the rest of the components of the iteration formula can be obtained.

As it can be seen, using these two approximate methods in solving this equation also lead to the exact solution (Fig. 2) using more iteration.

Example 3: We consider a second-order non-linear wave equation:

 (58)

With the initial and boundary conditions posed are:

 (59)

Exact solution of this equation is:

 (60)

APPLICATION OF HOMOTOPY-PERTURBATION METHOD

A homotopy can be constructed as follows:

 (61)

Substituting Eq. 10 into 61 and rearranging the resultant equation based on powers of p-terms, one has:

 (62)

 (63)

 (64)

With the following conditions:

 (65)

With the effective initial approximation for v0 from the conditions (65) and solutions of Eq. 62-64 may be written as follows:

 (66)

 (67)

 (68)

In the same manner, the rest of components were obtained using the maple package.

According to the HPM, we can conclude that:

 (69)

Therefore, substituting the values of v0 (x, t), v1 (x, t) and v2 (x, t) from Eq. 66-69 yields:

 (70)

Application of variational iteration method

 (71)

Its stationary conditions can be obtained as follows:

 (72)

The lagrangian multiplier can therefore be identified as:

 (73)

As a result, we obtain the following iteration formula:

 (74)

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

 (75)

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

 (76)

Substituting Eq. 75-76 and after simplification, we have:

 (77)

In the same way, we obtain u2 (x, t), u3 (x, t) as follows:

 (78)

 Fig. 3: 3D obtained results of u(x,t) by the HPM, VIM and the exact solution

 (79)

And so on. In the same way the rest of the components of the iteration formula can be obtained.

As it can be seen, using these two approximate methods in solving this equation also lead to the exact solution (Fig. 3) using more iteration.

CONCLUSION

The homotopy perturbation method and variational iteration method have been successfully used to study three types of partial differential wave equations. These equations describe the propagation of a wave (disturbance) and it arises in a wide variety of physical problems. The results obtained here were compared with the exact solutions. It can be easily seen that applying these two methods to the mentioned equations were led to the exact solution in all examples, so the results revealed that the homotopy perturbation method and variational iteration method are powerful mathematical tools for solutions of differential equations in terms of accuracy and efficiency.

REFERENCES

1:  Bildik, N. and A. Konuralp, 2006. The use of variational iteration method, differential transform method and adomian decomposition method for solving different types of nonlinear partial differential equations. Int. J. Nonlinear Sci. Numer. Simul., 7: 65-70.

2:  Debnath, L., 1997. Nonlinear partial differential equations for scientists and engineering. Birkhauser Boston, 18: 141-230.

3:  Ganji, D.D. and A. Sadighi, 2006. Application of He’s homotopy-perturbation method to nonlinear coupled systems of reaction-diffusion equations. Int. J. Nonlinear Sci. Numer. Simul., 7: 411-418.

4:  Ganji, D.D. and A. Sadighi, 2007. Application of homotopy-perturbation and variational iteration methods to nonlinear heat transfer and porous media equations. J. Comput. Applied Math., 207: 24-34.

5:  Ganji, D.D., M. Jannatabadi and E. Mohseni, 2007. Application of He's variational iteration method to nonlinear Jaulent-Miodek equations and comparing it with ADM. J. Comput. Applied Math., 207: 35-45.

6:  He, J.H., 1999. Variational iteration method: A kind of non-linear analytical technique. Some examples. Int. J. Nonlinear Mech., 344: 699-699.

7:  He, J.H., 1999. Homotopy perturbation technique. Comput. Methods Applied Mech. Eng., 178: 257-262.

8:  He, J.H., 2000. Variational iteration method for autonomous ordinary differential systems. Applied Math. Comput., 114: 115-123.

9:  He, J.H., 2006. New Interpretation of homotopy perturbation method. Int. J. Mod. Phys. B, 20: 2561-2568.

10:  Inc, M., Y. Cherrualt and K. Abbaoui, 2004. A computational approach to the wave equations: An application of the decomposition method. Kybernates, 33: 80-97.

11:  Sweilam, N.H. and M.M. Khader, 2007. Variational iteration method for one dimensional nonlinear thermoelasticity. Chaos Solitons Fract., 32: 145-149.

12:  Zhang, L.N. and J.H. He, 2006. Homotopy perturbation method for the solution of the electrostatic potential differential equation. Math. Problems Eng., 2006: 1-6.

13:  Momani, S. and S. Abuasad, 2006. Application of He's variational iteration method to Helmholtz equation. Chaos Soliton Fract., 27: 1119-1123.