Abstract and Applied Analysis

Volume 2014 (2014), Article ID 705204, 8 pages

http://dx.doi.org/10.1155/2014/705204

Research Article

Petrov-Galerkin Method for the Coupled Schrödinger-KdV Equation

M. S. Ismail, Farida M. Mosally, and Khadeejah M. Alamoudi

Department of Mathematics, College of Science, King Abdulaziz University, P. O. Box 80203, Jeddah 21589, Saudi Arabia

Received 27 March 2014; Accepted 24 May 2014; Published 15 June 2014

Academic Editor: Mohammad T. Darvishi

Copyright © 2014 M. S. Ismail et al. 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.

Petrov-Galerkin method is used to derive a numerical scheme for the coupled Schrödinger-KdV (SKdV) equations, where we have used the cubic B-splines as a test functions and a linear B-splines as a trial functions. Product approximation technique is used to deal with the nonlinear terms. An implicit midpoint rule and the Runge-Kutta method of fourth-order (RK4) are used to discretize in time. A block nonlinear pentadiagonal system is obtained. We solve this system by the fixed point method. The resulting scheme has a fourth-order accuracy in space direction and second-order in time direction in case of the implicit midpoint rule and it is unconditionally stable by von Neumann method. Using the RK4 method the scheme will be linear and fourth-order in time and space directions, and it is also conditionally stable. The exact soliton solution and the conserved quantities are used to assess the accuracy and to show the robustness and the efficiency of the proposed schemes.

1. Introduction

The coupled nonlinear Schrödinger-KdV equation has attracted extensive interest in physics and mathematics. The nonlinear Schrödinger-KdV equation describing the nonlinear dynamics.

Describing the nonlinear dynamics of one-dimensional Langmuir and ion-acoustic waves in a system of coordinates moving at the ion-acoustic speed has the following form [1, 2]:

and the initial condition

where and are constants, together with the homogenous Dirichlet boundary conditions

Analytical and numerical methods are very important tools to understand the feature and the behavior of the nonlinear wave equations. Many numerical methods have been used to solve numerically the single nonlinear Schrödinger and the single KdV equation using finite element and finite difference methods [3–6]. Ismail et al. [7–12] solved numerically the coupled nonlinear Schrodinger equation and the coupled KdV equation using the finite difference and finite element methods. Recently, Bhatt and Khaliq [13] solve the coupled nonlinear Schrodinger using high order exponential differencing method.

In this work we are going to study numerical the SKdV equation. Analytical solutions for this system using different methods are given in [14–16]. For numerical methods for the SKdV equation the literature is not rich. Appert and Vaclavik [17] solved the SKdV equations using a Crank-Nicolson scheme. A very good numerical works by Bai and Zhang [1, 2], they have solved this system using a split-step quadratic B-spline finite element method.

The exact solution of coupled Schrödinger-KdV equations (1) is

where is a free positive parameter. The SKdV equations has the following conserved quantities [1, 2, 18]:

To avoid complex computation, we assume

where , and are real functions. This will reduce Schrödinger-KdV equations to the coupled system

2. Numerical Method

Application of the numerical method requires truncation of the infinite interval to a finite interval . We assume that the solution of the coupled Schrödinger-KdV equation is negligible outside this interval. Also we assume all space derivatives at the boundaries approaches to zero in the region .

A standard weak formulation [3–5] of this problem is obtained by multiplying (7) by a twice differentiable test function and integrating by parts to obtain

where denotes the usual inner product

The space interval is discretized by uniform grid points

where the grid spacing is given by . We introduce finite elements in space in (7) and approximate the exact solution of the SKdV equation by

and the product approximation technique is used for treating the nonlinear terms in the following manner:

The trial functions , are chosen to be the piecewise linear functions

The unknown functions are determined from the numerical solution of the first-order ordinary differential system

where .

We choose the test functions ,  , to be the cubic B-spline with compact support

Direct calculations of the inner products in (14), using (11)–(13) and (15), will produce the first-order ordinary differential system

where

The first order ordinary differential system (16) can be solved by many ordinary differential equation solver. In this work, we will choose two methods, the implicit midpoint rule and the explicit Runge-Kutta method of fourth-order (RK4). These methods will be discussed in the following subsections.

2.1. The Implicit Midpoint Rule

By making use of the following substitution for the implicit midpoint rule

where is the time step size, into the system (16), this will lead us to the nonlinear block pentadiagonal system

The numerical solution of the nonlinear system (19) can be obtained by many iterative methods like Newton’s method, fixed point method, and in this work we will adopt the fixed point method and this will be given next.

2.2. Fixed Point Method

The fixed point method for solving the nonlinear system (19) can be given as follows:

where

We apply the iterative scheme (20) until the condition

is satisfied. Concerning the accuracy of the scheme, the scheme is of second-order in time and fourth-order in space.

2.3. Stability of the Scheme

To study the stability of the resulting scheme (19), we use von Neumann stability method. This method can only be applied for linear schemes, so we consider the linear version of the proposed scheme (19) by freezing all terms which make the scheme nonlinear. The linear version of the proposed method can be displayed as follows:

where ,  , and is assumed to be constant on the whole range.

We consider the first two equations and the third equation can be done in the similar way. We assume that the solution of (23) to be of the form

By substituting this into (23), this will lead us after some manipulations to the following system:

where

The system (25) can be written in a matrix-vector form as

where and is the matrix

where

The von Neumann stability condition for the system (27) is the maximum modulus of the eigenvalues of the matrix are to be less than or equal to one. The eigenvalues of the matrix are

with modulus equal to one. This means the scheme which we have derived is unconditionally stable according to von Neumann stability analysis.

The same methodology of stability can be applied for the third equation.

2.4. Runge-Kutta Method of Fourth Order (RK4)

Another approach for solving the semidiscrete system (16) is to use the Runge-Kutta method of fourth-order. The semidiscrete system (16) can be written in a matrix vector form as

where is a pentadiagonal matrix of the form

which can be in a compact form as

where

The RK4 method for the block system (33) can be given as follows:

where

So in order to calculate , and , we need to solve a linear pentadiagonal system, and this can be done by Crout’s method as follows.

At the beginning of the calculations, we factor the matrix into , where and are lower and upper triangular matrices one time only, and then we left only with a lower and an upper triangular systems which can be solved easily by forward and backward substitution process, respectively. This can be easily accomplished with minimum cost. The RK4 method is conditionally stable. In order to get stable results we choose the time step as . Regarding the accuracy in this case the scheme is of fourth-order in time and space directions.

3. Numerical Results

In this section, we present some numerical results for the proposed schemes. The accuracy of our methods is tested by calculating the error norms

The conserved quantities of the coupled Schrödinger-KdV equation are calculated using trapezoidal rule. In all calculations, we choose the initial condition

where ; and are free positive parameters. The following tests will be discussed

(i) Implicit Midpoint Rule with . To compute the numerical solution using the system obtained by using the implicit midpoint rule (19), the following parameters are used:

In Tables 1 and 2, we display the errors and the conserved quantities. In Figures 1 and 2, we display the numerical solution of and for . The maximum of and is constant and equal, respectively, to for . It is very easy to see the accuracy and how the conserved quantities are almost conserved.

Table 1: Errors using implicit midpoint rule.

Time

1

2

5

8

Table 2: Conservative quantities.

Time

0
2.749231

2
2.749229

5
2.749226

8
2.749223

Figure 1: The evolution of the numerical solution of with .

Figure 2: The evolution of the numerical solution of with .

(ii) Implicit Midpoint Rule with . We compute the numerical solution using the implicit midpoint rule. We choose in this case , we display the numerical solution of and for . In Figures 3 and 4, we display the real ( and imaginary ( parts of the numerical solution, the same maximum values obtained as the case of . This is compatible with the exact solution, and we have noticed some differences between our results in Figures 3 and 4 with the one in [2]. We think this occurrence is due to their selected interval where boundary values are not zero. In Figures 5 and 6, we display the simulation of the numerical solution of the modulus of and , respectively.

Figure 3: The evolution of the numerical solution of real part of and .

Figure 4: The evolution of the numerical solution of imaginary l part of and .

Figure 5: The evolution of the numerical solution of with .

Figure 6: The evolution of the numerical solution of with .

(iii) Runge-Kutta Method of Fourth-Order. Now we will test the method obtained by using the RK4. We choose the following set of parameters:

The errors and conserved quantities are displayed in Tables 3 and 4; we can easily see the high accuracy of the proposed method and how the conserved quantities are exactly conserved.

Table 3: Errors using Runge-Kutta method.

Time

0.2

0.6

0.8

1.0

Table 4: Conservative quantities.

Time

0
2.749231

2
2.749231

5
2.749231

8
2.749231

(iv) Comparison of Some Existing Methods. A comparative study has been conducted with some existing methods and the results are displayed in Table 5. From the numerical results displayed for different methods, we can see the superiority of our proposed method comparing to the previous methods which are mentioned in [1, 2].

Table 5: The comparison of maximum errors using different schemes at .1 with .

Method

Present (with midpoint rule)

Present (with RK4)

SSQBS FEM

SSCBS FEM

Semidiscrete FEM

Crank Nicolson

4. Conclusion

In this work, we derived highly accurate numerical schemes for the SKdV equation, using finite element method (Petrov-Galerkin) with product approximation technique, where we have used the linear hat functions as trial function and the cubic B-spline as test functions. The differential systems obtained are solved by using implicit midpoint rule. Also we have solved the first-order ordinary differential system using the RK4, and in this case the numerical scheme is conditionally stable and of fourth-order in space and time directions. We held a comparison with some existing methods, and we found that our scheme produced highly accurate results and conserved the conserved quantities almost exactly.

Conflict of Interests

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

References

D. Bai and L. Zhang, “Numerical studies on a novel split-step quadratic B-spline finite element method for the coupled Schrödinger-KdV equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 3, pp. 1263–1273, 2011.

D. Bai and L. Zhang, “The finite element method for the coupled Schrödinger-KdV equations,” Physics Letters A, vol. 373, no. 26, pp. 2237–2244, 2009.

I. Christie, D. F. Griffiths, A. R. Mitchell, and J. M. Sanz-Serna, “Product approximation for nonlinear problems in the finite element method,” IMA Journal of Numerical Analysis, vol. 1, no. 3, pp. 253–266, 1981.

J. M. Sanz-Serna and I. Christie, “Petrov-Galerkin methods for nonlinear dispersive waves,” Journal of Computational Physics, vol. 39, no. 1, pp. 94–102, 1981.

S. W. Schoombie, “Spline Petrov-Galerkin methods for the numerical solution of the Korteweg-de Vries equation,” IMA Journal of Numerical Analysis, vol. 2, no. 1, pp. 95–109, 1982.

H. Wang, “Numerical studies on the split-step finite difference method for nonlinear Schrödinger equations,” Applied Mathematics and Computation, vol. 170, no. 1, pp. 17–35, 2005.

M. S. Ismail, “Numerical solution of a coupled Korteweg-de Vries equations by collocation method,” Numerical Methods for Partial Differential Equations, vol. 25, no. 2, pp. 275–291, 2009.

M. S. Ismail, “A fourth-order explicit schemes for the coupled nonlinear Schrödinger equation,” Applied Mathematics and Computation, vol. 196, no. 1, pp. 273–284, 2008.

M. S. Ismail, “Numerical solution of coupled nonlinear Schrödinger equation by Galerkin method,” Mathematics and Computers in Simulation, vol. 78, no. 4, pp. 532–547, 2008.

M. S. Ismail and T. R. Taha, “A linearly implicit conservative scheme for the coupled nonlinear Schrödinger equation,” Mathematics and Computers in Simulation, vol. 74, no. 4-5, pp. 302–311, 2007.

M. S. Ismail and S. Z. Alamri, “Highly accurate finite difference method for coupled nonlinear Schrödinger equation,” International Journal of Computer Mathematics, vol. 81, no. 3, pp. 333–351, 2004.

M. S. Ismail and T. R. Taha, “Numerical simulation of coupled nonlinear Schrödinger equation,” Mathematics and Computers in Simulation, vol. 56, no. 6, pp. 547–562, 2001.

H. P. Bhatt and A. Q. M. Khaliq, “Higher order exponential time differencing scheme for system of coupled nonlinear Schrödinger equations,” Applied Mathematics and Computation, vol. 228, pp. 271–291, 2014.

R. Abardeh, M. Ramezanpour, A. Doosthoseini, and E. Rezaie, “New method for solving coupled Schrödinger KdV equation,” Applied Mathematical Sciences, vol. 7, no. 86, pp. 4273–4280, 2013.

D. Kaya and S. M. El-Sayed, “On the solution of the coupled Schrödinger-KdV equation by the decomposition method,” Physics Letters A, vol. 313, no. 1-2, pp. 82–88, 2003.

M. Labidi, G. Ebadi, E. Zerrad, and A. Biswas, “Analytical and numerical solutions of the Schrödinger-KdV equation,” Pramana—Journal of Physics, vol. 78, no. 1, pp. 59–90, 2012.

K. Appert and J. Vaclavik, “Dynamics of coupled solitons,” Physics of Fluids, vol. 20, no. 11, pp. 1845–1849, 1977.

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