Laplace Decomposition Method to Study Solitary Wave Solutions of Coupled Nonlinear Partial Differential Equation
Arun Kumar1and Ram Dayal Pankaj2
Academic Editor: L. S. Heath, P. Amodio
Received02 May 2012
Accepted13 Jun 2012
Published04 Sept 2012
Abstract
Analytical and numerical solutions are obtained for coupled nonlinear partial differential equation by the well-known Laplace decomposition method. We combined Laplace transform and Adomain decomposition method and present a new approach for solving coupled SchrΓΆdinger-Korteweg-de Vries (Sch-KdV) equation. The method does not need linearization, weak nonlinearity assumptions, or perturbation theory. We compared the numerical solutions with corresponding analytical solutions.
1. Introduction
Systems of partial differential equations have attracted much attention in a variety of applied sciences because of their wide applicability. These systems were formally derived to describe wave propagation to model the shallow water waves [1β5] and to examine some chemical reaction-diffusion model of Brusselator [4β6]. Some of the commonly used methods to solve these equations are the method of characteristics, the Riemann invariants, and Adomian decomposition method [6].
In this work, we used Laplace decomposition method introduced by Khuri [7, 8] which is further used by Yusufoglu to solve Duffing equation [9] and Elgasery for Falkner-Skan equations [10]. This technique, modified by Hussain and Khan [11], illustrates how the Laplace transform may be used to approximate the solutions of the nonlinear partial differential equations by extending the decomposition method [12, 13].
2. Laplace Decomposition Method (LDM)
In this section, we outline the main steps of the method. We consider the nonlinear partial differential equations in an operator form:
With initial data
where is a first-order partial differential operator, , and , are linear and nonlinear operators, respectively, and are the source terms. By applying the Laplace transform to both sides of (1) and using initial conditions (2), we have
By the differentiation property of Laplace transform, we get
where ββ ββis Laplace domain function. Solutions and in LDM are defined as
The nonlinear terms , represented by infinite series,
are the Adomian polynomials [14], generated for all forms of nonlinearity and they are determined by the following relations:
Substituting (5) and (6) into (4), gives
Applying the linearity of the Laplace transform, we deduct the following recursive relations:
In general, for , the recursive relations are given by
Applying the inverse Laplace transform, we can evaluate and (). In some cases the exact solution in the closed form can be obtained.
3. Application
At the classical level, a set of coupled nonlinear wave equations describes the interaction between high-frequency and low-frequency waves [15], and the calculation of exact and numerical solutions of the equations, in particular, travelling wave solutions, plays an important role in wave-wave interaction and soliton theory [1, 16].
We consider the SchrΓΆdinger-KdV (Sch-KdV) equation as a model for the interaction of long and short nonlinear waves:
With initial conditions
where , are arbitrary constant.
With Laplace decomposition method on (11) and using the differentiation property of Laplace transform, initial conditions and the inverse Laplace transforms are
where ,, and β ββ β are adomian polynomials that represent nonlinear terms.
So the recursive relation is deduced as
By this recursive relation we can find other components of the solution as
The other components of the decomposition series can be determined in a similar way; we can obtain the expression of in a Taylor series, which gives the closed form solutions as:
4. Numerical Description of the Solution
The Laplace decomposition method is used for finding the exact and approximate travelling-waves solutions of the Sch-KdV equation. Both the exact and approximate solutions obtained for n = 2 using LDM are plotted in Figure 1. It is evident that when compute more terms for the decomposition series the numerical results are getting much closer to the corresponding analytical solutions.
(a)
(b)
(c)
(d)
5. Conclusion
The Laplace decomposition method is a powerful method which has provided an efficient potential for the solution of physical applications modeled by nonlinear differential equations. The algorithm can be used without any need to complex calculations except for simple and elementary operations.
References
L. Debnath, Nonlinear Partial Differential Equations for Scientists and Engineers, Birkhäauser, Boston, Mass, USA, 1997.
J. D. Logan, An Introduction to Nonlinear Partial Differential Equations, Wiley-Inter Science, New York, NY, USA, 1994.
G. B. Whitham, Linear and Nonlinear Waves, Wiley, New York, NY, USA, 1974.
C. Lubich and A. Ostermann, βMulti-grid dynamic iteration for parabolic equations,β BIT Numerical Mathematics, vol. 27, no. 2, pp. 216β234, 1987.
S. Vandewalle and R. Piessens, βNumerical experiments with nonlinear multigrid waveform relaxation on a parallel processor,β Applied Numerical Mathematics, vol. 8, no. 2, pp. 149β161, 1991.
S. A. Khuri, βA Laplace decomposition algorithm applied to a class of nonlinear differential equations,β Journal of Applied Mathematics, vol. 1, no. 4, pp. 141β155, 2001.
E. Yusufoglu, βNumerical solution of Duffing equation by the Laplace decomposition algorithm,β Applied Mathematics and Computation, vol. 177, no. 2, pp. 572β580, 2006.
A. M. Wazwaz, Partial Differential Equations Methods and Applications, Balkema, Rotterdam, The Netherlands, 2002.
A. M. Wazwaz, βThe variational iteration method for solving linear and nonlinear systems of PDEs,β Computers and Mathematics with Applications, vol. 54, no. 7-8, pp. 895β902, 2007.
A. Kumar, βA variational description of non-linear Schrödinger Equation,β International Journal of Computational and Applied Mathematics, vol. 4, no. 2, pp. 159β164, 2009.