Advances in Mathematical Physics
Volume 2016 (2016), Article ID 7909632, 8 pages
http://dx.doi.org/10.1155/2016/7909632

Research Article

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

Saeed Ahmed and Muhammad Kalim

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

where

In (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

where

It 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 write

where

The higher order deformation equation (12) governs and so is valid for any nonlinear operator . So, we can let

which then becomes

Substituting (15) into (12), we have

Now, 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, let

Now, 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 as

Consider that the inverse operator exists and is defined as

Applying 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 formula

where 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) as

Substituting 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.

Table 1: The values of at , = 0.1–0.9.

Exact

= 0.1
0.01040032
0.01042119
0.01011233
0.01042119

= 0.2
0.04321024
0.04356014
0.04090019
0.04356076

= 0.3
0.10087776
0.10274115
0.09304346
0.10275264

= 0.4
0.18592772
0.19214213
0.16722674
0.19223594

= 0.5
0.30100029
0.31705000
0.26413946
0.31754163

= 0.6
0.44888979
0.4841644
0.3844759
0.48612181

= 0.7
0.63258414
0.70194985
0.5289352
0.70843802

= 0.8
0.85530539
0.98103624
0.69822148
1.00000000

= 0.9
1.12055235
1.33466772
0.89304379
1.38562172

Table 2 shows the absolute error of these solutions at different values of . Again, we see that the approximate solution of ADM has less percentage error than others. The percentage error is increasing as is getting close to 1.25 due to the domain of .

Table 2: The absolute error of at , = 0.1–0.9.

= 0.1
0.00002087171
0.00000000451
0.00030886334

= 0.2
0.00035052253
0.00000062180
0.00266056974

= 0.3
0.00187487820
0.00001149230
0.00029081445

= 0.4
0.00187480000
0.00001150000
0.00970910000

= 0.5
0.00630820000
0.00009380000
0.02500920000

= 0.6
0.01654130000
0.00049160000
0.05340210000

= 0.7
0.03723200000
0.00195740000
0.10164590000

= 0.8
0.07585390000
0.00648810000
0.17950280000

= 0.9
0.14469460000
0.01896380000
0.30177850000

In Figure 1, the graph shows the approximate solutions by the iterations methods together with the exact solution (purple line). We see that the approximate solution by ADM (red line) is closer to the exact solution followed by HAM (blue line) and VIM (green line). In Figure 2, the error plot is given and we again see that the absolute error by ADM (red line) is the smallest, while the absolute error of VIM (green line) is the highest.

Figure 1: The exact solution is in purple line, HAM is in blue, ADM is in red, and VIM is in green.

Figure 2: The error graph of HAM is in blue, ADM is in red, and VIM is in green.

5. Conclusion

A comparative study has been made to seek the semianalytical solution of the ZK equation with initial and boundary conditions. The solution of the ZK equation is obtained in the form of an infinite series that converges rapidly in its domain. This form of solution gives the acoustic pressure of the beam in terms of propagation coordinate and time more explicitly. Since such form of the solution has not been obtained before, we believe that this will open new ways to understand the propagation of the confined beam in the nonlinear medium. Different iteration methods like HAM, ADM, and VIM have been applied and compared. The result revealed that ADM generated surprisingly effective results.

Conflict of Interests

The authors declare that they there is no conflict of interests regarding the publication of this paper.

References

O. V. Rudenko, “The 40th anniversary of the Khokhlov-Zabolotskaya equation,” Acoustical Physics, vol. 56, no. 4, pp. 457–466, 2010.

O. V. Bessonova, V. A. Khokhlova, M. R. Bailey, M. S. Canney, and L. A. Crum, “Focusing of high power ultrasound beams and limiting values of shock wave parameters,” Acoustical Physics, vol. 55, no. 4-5, pp. 463–473, 2009.

M. Kumar, R. Kumar, and A. Kumar, “On similarity solutions of Zabolotskaya-Khokhlov equation,” Computers & Mathematics with Applications, vol. 68, no. 4, pp. 454–463, 2014.

M. Tajiri, “Similarity reductions of the Zabolotskaya-Khokhlov equation with a dissipative term,” Journal of Nonlinear Mathematical Physics, vol. 2, no. 3-4, pp. 392–397, 1995.

I. Kostin and G. Panasenko, “Khokhlov-Zabolotskaya-Kuznetsov-type equation: nonlinear acoustics in heterogeneous media,” SIAM Journal on Mathematical Analysis, vol. 40, no. 2, pp. 699–715, 2008.

A. Rozanova-Pierrat, “Qualitative analysis of the Khokhlov-Zabolotskaya-KUZnetsov (KZK) equation,” Mathematical Models & Methods in Applied Sciences, vol. 18, no. 5, pp. 781–812, 2008.

Z. Xiaofeng and R. J. McGough, “The Khokhlov—Zabolotskaya—Kuznetsov (KZK) equation with power law attenuation,” in Proceedings of the IEEE International Ultrasonics Symposium (IUS '14), pp. 2225–2228, Chicago, Ill, USA, September 2014.

T. Y. Cheung, Application of the Khokhlov-Zabolotskaya-Kuznetsov equation to modeling high-intensity focused ultrasound beams [M.S. thesis], Boston University, Boston, Mass, USA, 2008.

W. I. Fushchich, W. M. Shtelen, and N. I. Serov, Symmetry Analysis and Exact Solutions of Equations of Nonlinear Mathematical Physics, vol. 246 of Mathematics and its Applications, Kluwer Academic, Dordrecht, The Netherlands, 1993.

N. Euler and W.-H. Steeb, Continuous Symmetries, Lie Algebras and Differential Equations, Bibliographisches Institut, Mannheim, Germany, 1992.

J. Biazar and F. Azimi, “He's homotopy perturbation method for solving Helmholtz equation,” International Journal of Contemporary Mathematical Sciences, vol. 3, no. 13–16, pp. 739–744, 2008.

R. Courant and D. Hilbert, Methods of Mathematical Physics, Vol 2: Partial Differential Equations, Interscience Publishers, New York, NY, USA, 1962.

L. Debnath, Nonlinear Partial Differential Equations for Scientists and Engineers, Birkhäuser, Boston, Mass, USA, 1997.

L. C. Evans, Partial Differential Equations, Graduate Studies in Mathematics, vol. 19, American Mathematical Society, 1999.

F. John, Partial Differential Equations, Springer, New York, NY, USA, 1982.

M. Khan, M. A. Gondal, and S. Kumar, “A novel homotopy perturbation algorithm for linear and nonlinear system of partial differential equations,” World Applied Sciences Journal, vol. 12, no. 12, pp. 2352–2357, 2011.

J.-H. He, “Variational principles for some nonlinear partial differential equations with variable coefficients,” Chaos, Solitons & Fractals, vol. 19, no. 4, pp. 847–851, 2004.

G. Adomian and R. Rach, “Noise terms in decomposition solution series,” Computers & Mathematics with Applications, vol. 24, no. 11, pp. 61–64, 1992.

G. Adomian, “A review of the decomposition method in applied mathematics,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 501–544, 1988.

G. Adomian and R. Rach, “Modified decomposition solution of linear and nonlinear boundary-value problems,” Nonlinear Analysis: Theory, Methods & Applications, vol. 23, no. 5, pp. 615–619, 1994.

J. Biazar and H. Ebrahimi, “An approximation to the solution of hyperbolic equations by Adomian decomposition method and comparison with characteristic method,” Applied Mathematics and Computation, vol. 163, no. 2, pp. 633–638, 2005.

H. Gu and Z.-B. Li, “A modified Adomian method for system of nonlinear differential equations,” Applied Mathematics and Computation, vol. 187, no. 2, pp. 748–755, 2007.

R. Jebari, I. Ghanmi, and A. Boukricha, “Adomian decomposition method for solving nonlinear heat equation with exponential nonlinearity,” International Journal of Mathematical Analysis, vol. 7, no. 15, pp. 725–734, 2013.