Advances in Mathematical Physics

Volume 2016 (2016), Article ID 7909632, 8 pages

http://dx.doi.org/10.1155/2016/7909632

## A Comparative Approach to the Solution of the Zabolotskaya-Khokhlov Equation by Iteration Methods

National College of Business Administration and Economics, Gulberg-III, 54660 Lahore, Pakistan

Received 9 October 2015; Revised 26 December 2015; Accepted 27 December 2015

Academic Editor: Soheil Salahshour

Copyright © 2016 Saeed Ahmed and Muhammad Kalim. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We employed different iteration methods like Homotopy Analysis Method (HAM), Adomian Decomposition Method (ADM), and Variational Iteration Method (VIM) to find the approximate solution to the Zabolotskaya-Khokhlov (ZK) equation. Iteration methods are used to solve linear and nonlinear PDEs whose classical methods are either very complex or too limited to apply. A comparison study has been made to see which of these methods converges to the approximate solution rapidly. The result revealed that, amongst these methods, ADM is more effective and simpler tool in its nature which does not require any transformation or linearization.

#### 1. Introduction

In science, physics, and engineering, many problems are modeled by the means of nonlinear evolutions equations especially in plasma physics, fluid dynamics, biology, and nonlinear acoustics. In this paper, we consider one such nonlinear partial differential equation which is (1 + 1) dimensional Zabolotskaya-Khokhlov (ZK) equation (ZK for short) This equation is one of the basic equations in nonlinear acoustic and nonlinear wave theory. It is named after two Russian mathematicians R. V. Khokhlov and E. A. Zabolotskaya who derived it for the first time almost four decades ago and presented an approximate solution which described some interesting features of nonlinear wave beams [1]. In recent years, an extensive research is being made in nonlinear acoustic due to the development of new medical devices for nonlinear diagnostic ultrasound imaging, for noninvasive destruction of tumors in acoustic surgery and for kidney stone comminution. All of the above applications rely on the focusing acoustic waves in nonlinear medium which can be studied by the ZK equation [2]. The ZK equation describes the propagation of a confined wave beam or sound beam in nonlinear medium without dispersion or absorption and it investigates the beam deformation associated with the properties of the nonlinear medium [3, 4]. This equation has applications in many fields of life such as it is used to simulate the estimation of fish stock abundance, discrimination between fish species, the effect of excess attenuation which occurs due to the nonlinear sound propagation in water, absorption, and diffraction in the focused sound beams [5–8].

The solution of the ZK equation will shed light on some new features in the behavior of nonlinear beam. Many authors [9, 10] have found the infinitesimal symmetries and exact solution of the ZK equation by means of Lie Algebras, but the analytical or approximate solution of the ZK equation in the form of an infinite series has not been obtained before. Keeping this limitation in our mind, in this work, we will find the analytical solution of the ZK equation in the form of an infinite series with the help of iteration methods. This analytical solution would be extremely useful in the areas of theoretical physics, applied mathematics, and engineering as it will express the acoustic pressure of the beam in terms of propagation coordinate and time explicitly.

As we know that many physical problems are modeled by means of partial differential equations (PDEs) like Helmholtz equation, heat equation, wave equation, gas equation, and so forth [11–15]. Of these PDEs, some are linear, but many are quasi-linear or nonlinear [16] whose analytical solution techniques are limited to apply. There is no general method of solving such PDEs; however an extensive research is being made to find exact or approximate solutions of nonlinear partial differential equations with the aid of iteration methods like variational principal [17], decomposition solution [18, 19], and Adomian Decomposition Method [20–23]. These iteration methods are analytical tool that approximate the solution of such PDEs in terms of an infinite series which converges rapidly to the exact solution of the PDE in just few iterations. The accuracy and the speed of convergence of the resulting series are surprisingly high and one needs only a few iterations to approximate the exact solution as shown by the example given in Section 4.

Consider the initial and boundary conditions for where and are two known functions.

Integrating twice and applying the initial conditions (1), we can write whereIn (2), we assume that is bounded function for all in and in , where . The functions and are Lipschitz continuous such that

#### 2. Existence and Uniqueness

In this section, we will analyze uniqueness and existence of the solution to the ZK equation.

Theorem 1. *Equation has a unique solution for .*

*Proof. *We suppose that and are two solutions of (2), and then Since , which proves that has a unique solution.

#### 3. Analysis of the Methods

In this section, we will introduce the iteration methods which we will apply on the ZK equation. We will apply Homotopy Analysis Method (HAM), Adomian Decomposition Method (ADM), and Variational Iteration Method (VIM) to solve the said equation.

##### 3.1. Homotopy Analysis Method (HAM)

To illustrate the basic idea of the Homotopy Analysis Method, we consider the following nonlinear differential equation:where is any differential or integral operator and is the unknown function. Using Homotopy Analysis Method, we can construct the so-called zeroth-order deformation equation where is an initial guess to the exact solution of the given differential equation and is the embedding parameter, whereas is an auxiliary parameter, is an auxiliary linear operator, is an auxiliary function, and is the unknown function. Observe that when , relation (7) turns out to be , and when it becomes . So, as the embedding parameter changes from 0 to 1, we see that continuously changes from the initial condition to the exact solution of the differential equations.

Using Taylor theorem, we can express as a power series of whereIt is obvious that the convergence of the above relation depends mainly on the auxiliary parameter . With the proper choice of the initial guess , the power series (8) converges to the exact solution at Next, we substitute (8) into (7) to get Differentiating (11) times with respect to and dividing both sides by and finally applying , we can writewhere The higher order deformation equation (12) governs and so is valid for any nonlinear operator . So, we can letwhich then becomes Substituting (15) into (12), we haveNow, taking as initial guess , , , and , together with using the value of , we can write the recurrence relation from the last equation as follows: The set of equations in (17) can be solved by simple integration. The resulting components , will then be substituted in (8) with as unity, which will lead towards the solution of in the form of an infinite series whose convergence is sure in most cases in its own domain.

As we see, the initial and boundary conditions play vital role in determining the series solution and different set of conditions will result in different series solutions. Though we do not have general criteria of imposing restrictions at these conditions at the start, one should be careful in choosing initial conditions.

##### 3.2. Adomian Decomposition Method (ADM)

Adomian Decomposition Method was first introduced by Adomian in 1984 [18, 19]. Like HAM, in this method, we split the nonlinear partial differential equation into two parts, namely, linear and nonlinear. Applying the inverse linear operator on both sides of the equation and decomposing the unknown function in the form of a series, we get a solution. In ADM, we also decompose the nonlinear terms in the form of Adomian polynomials and use a special formula to find these unknown terms due to which, no doubt, we need to do extra calculation than we do in HAM.

Consider the following equation: Equation (18) can be split into the following form: where and are linear and nonlinear differential operators, respectively. Consider that the inverse operator exists, and applying it on (19), we get where is the initial guess of the equation.

By the Adomian Decomposition Method [16], we can define the unknown function The nonlinear term of the equation is given by the formula [18]where , are the Adomian polynomials which can be retrieved by the formula Substituting (21) and (22) into (20), we can write Relation (24) can be used to determine the solution of (18).

##### 3.3. Variational Iteration Method (VIM)

Under Variational Iteration Method, we break the given (6) into the following form: where and are nonlinear and linear operators, respectively. According to VIM, a correctional functional can be constructed as follows:where is the Lagrange multiplier which can be obtained optimally by using variation theory. The function is the restricted variation; that is , which helps us to determine the Lagrange multiplier together with integration by parts. The exact solution may be obtained by Applying (26) on (2), we can write where .

Next we find the value of Lagrange multiplier by taking variation on (28)From (29), due to stationary conditions, we get, after a bit of calculations,which means that the Lagrange multiplier is , so from (28) we have We can find the successive values of from relation (31).

##### 3.4. Theorem

The series solution using VIM converges with , .

*Proof. *Consider Removing the subscripts, we can again write Subtracting (33) from (32), we get Now, letNow, observe that is decreasing with respect to , so applying mean value theorem, we can write from (34) the following:Therefore,Since , , so the series converges and hence the proof.

#### 4. Application and Results

In this section, we solve the ZK equation by the above-mentioned methods. We also compare the results obtained to see which method is more sufficient and converges rapidly to the exact solution.

##### 4.1. Example

Consider ZK equation with together with initial and boundary conditions given as follows:Here is the unknown function to be determined which represents the acoustic pressure and is the propagation coordinate. The term denotes the time evolution of the wave, while the nonlinear terms and describe the steepening of the wave.

*Using HAM*. We apply HAM on (38) to get the following iteration scheme due to (17):Using the above relation, we can find the other components as follows:Substituting these values in (8) and making unity, we get the solution of (38) in the form of the following infinite series: This represents the acoustic pressure of the beam in terms of propagation coordinate and time.

*Using ADM*. Let and so that (38) can be written asConsider that the inverse operator exists and is defined asApplying on (42) and using the initial conditions, we get By the Adomian Decomposition Method, we can define the unknown function The nonlinear term of the equation is given by the formulawhere are the Adomian polynomials which can be retrieved by the formula Substituting (45) and (46) into (44), we can write From (48), we get the recurrence relation as Now, using formula (47), we can find the first few Adomian polynomials as Applying the recurrence relation (49) together with the inverse operator (43) on the set of equations given in (50) to get the components of and inserting these components in (45), we get the solution to (38) asSubstituting these values in (8) and making unity, we get the series The exact solution to (38) can be written in a closed form as follows:

*Using VIM*. Using VIM, we can find the recurrence relation given as follows:Since we can do few more iterations and have the successive components of as follows:Consequently, the third-order approximate solution of (38) using VIM is given by Table 1 shows the comparison between the approximate solutions (41), (51), and (57) and the exact solution (53) of (38). From the table, we see that the approximate solution obtained by ADM is closer to the exact solution than others. After ADM, the approximate solution by HAM is better than that of VIM. Since we take the third-order solution from VIM, it is suggested that higher order solution obtained by VIM may rapidly lead to convergence towards the solution. If we look at the table carefully, we will see that these approximate solutions are close to the exact one for the values and after this the approximate solutions seem to be moving away to the exact solution. Since the exact solution of the equation given in (53) has the domain and by fixing , we get that the restriction on is . So, as the values of are getting close to 1.25, the solution is diverging away.