ISRN Applied MathematicsVolume 2012 (2012), Article ID 462731, 12 pagesdoi:10.5402/2012/462731

Research Article

The Spectral Method for Solving Sine-Gordon Equation Using a New Orthogonal Polynomial

Zoleikha Soori and Azim Aminataei

Department of Mathematics, K. N. Toosi University of Technology, P.O. Box 16315-1618, Tehran, Iran

Received 17 September 2011; Accepted 12 October 2011

Academic Editors: P. J. García Nieto and H.-T. Hu

Copyright © 2012 Zoleikha Soori and Azim Aminataei. 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 have presented a numerical scheme for solving one-dimensional nonlinear sine-Gordon equation. We apply the spectral method with a basis of a new orthogonal polynomial which is orthogonal over the interval with weighting function one. The results show the accuracy and efficiency of the proposed method.

1. Introduction

The sine-Gordon equation is a nonlinear hyperbolic partial differential equation (PDE). It was originally considered in the nineteenth century in the course of surface of constant negative curvature. This equation attracted a lot of attention in the 1970s due to the presence of soliton solutions [1, 2]. The sine-Gordon equation appears in a number of physical applications [3–5], including applications in the chain of coupled pendulums and modelling the propagation of transverse electromagnetic (TEM) wave on a superconductor transmission system. Consider the one-dimensional nonlinear sine-Gordon equation with the initial conditions and the boundary conditions In the last decade, several numerical schemes have been developed for solving (1.1), for instance, high-order solution of one-dimensional sine-Gordon equation using compact finite difference and DIRKN methods [3]. The authors of [6] proposed a numerical scheme for solving (1.1) using collocation and radial basis functions. Also, the boundary integral equation approach is used in [7]. Bratsos has proposed a numerical method for solving one-dimensional sine-Gordon equation and a third-order numerical scheme for the two-dimensional sine-Gordon equation in [8, 9], respectively. Also, in [10], a numerical scheme using radial basis function for the solution of two-dimensional sine-Gordon equation is used. In addition, several authors proposed spectral methods and Fourier pseudospectral method for solving nonlinear wave equation using a discrete Fourier series and Chebyshev orthogonal polynomials [11–13].

In this paper, we have proposed a numerical scheme for solving (1.1) using spectral method with a basis of a new orthogonal polynomial. This polynomial is introduced by Chelyshkov in [14].

The outline of this paper is as follows. In Section 2, we introduce a new orthogonal polynomial. In Section 3, we apply spectral method for solving (1.1). In Section 4, we show the numerical results. Finally, a brief conclusion is given in Section 5.

2. The New Orthogonal Polynomial

In this section, we briefly introduce the new orthogonal polynomial, that is introduced in [14].

Let be a fixed whole number, and is a sequence of polynomial, that is, that satisfy the orthogonality relationship and standardization The coefficient of the polynomial is defined by requirements (2.2), (2.3), and the Gram-Schmidt orthogonalization procedure (without normalization).The explicit definition of the polynomial is as follows: This yields Rodrigue’s type representation as follows:

and the orthogonality relationship (2.2) is confirmed by applying (2.5). It also follows from (2.5) thatThis polynomial can be connected to a fixed set of the Jacobi polynomials [15], that is,

3. The Proposed Method

In this section, the introduced orthogonal polynomial is applied as the basis for spectral method to solve sine-Gordon equation.

The approximate solution to the exact solution can be written in the form

where is the time-dependent quantities which must be determined. In this paper, we apply orthogonal polynomial for . From (3.1), we get

Using (3.2) and substituting in (1.1), we obtain

where and denote the second-order derivatives of and with respect to and , respectively, wherein the collocation points are as ,.

Equation (13) is rewritten as follows:

For displaying this equation in matrix form, for , and , we define

The matrices and are positive definite matrices of order . The system given in (3.4) consists of equations in unknowns. To obtain a unique solution for this system, we need two additional equations. For this, we add two boundary conditions, which are given by (1.3), in the following forms:

From (3.6), we get

Therefore, by using (3.6), and (3.7), a new system is obtained as follows:

where , , , and are vectors of order with components. The matrices and are of order , which are defined as follows: Also, the vectors and are defined as follows: The nonlinear second-order system of ODEs (3.8) can be solved numerically using the fourth-order A-stable DIRKN method [3, 16]. From the initial conditions (1.2), we determine the initial vectors and . Using (3.1), we get

Equation (3.12) is system of equations in unknowns. These equations can be written in the matrix form and , where

3.1. The DIRKN Method

For solving (3.8), we use DIRKN method [3, 16]:

where is the internal stages, and represent approximations to and , respectively. Also we definewhere . The RKN method can be denoted in notation by the table of coefficients:

In the DIRKN method, as and are equal. Also, in the above table of coefficients, is the number of stages. We use the four-stage A-stable DIRKN method [16], with algebraic fourth order, which is defined by the following table of coefficients for solving all experiments.

Using Table 1, we can calculate , , separately from (3.14) with Newton-Raphson iteration method. We iterate until the following criterion is satisfied: where is the value of at the kth iteration of the Newton-Raphson method.

Table 1: Four-stage A-stable DIRKN method.

0.6582
0.24

0.9300
0.186484967952
0.24

0.0700
−3.148664086195
2.945628022183
0.24

0.3418
1.491936101409
−1.816304852835
0.103007001689
0.24

0.108572237906
0.012764608972
0.169586947770
0.209076205352

0.317648443258
0.182351556742
0.182351556742
0.317648443258

4. Numerical Experiments

In this section, we present the results of numerical experiments using the method introduced in Section 3. The and error norms are defined as follows:

We choose , for solving all experiments.

The computations associated with the experiments discussed above were performed in Maple 13 on a PC with a CPU of 2.4 GHZ.

Experiment 1. In this experiment [3, 6, 10], we consider (1.1) with and the initial conditions with the boundary conditions (1.3), which can be obtained from the following exact solution: In Table 2, the approximate solution for several final times () with different locations is shown. Figure 1 shows two type plots at of Experiment 1.

Table 2: The approximate solution for of Experiment 1.

0.0
0.3986746100
0.7895822394
1.1658271779
1.5220255084
1.8545904360

0.2
0.3908822100
0.7744385442
1.1441349052
1.4948183070
1.8229915644

0.4
0.3689530692
0.7311732829
1.0827641117
1.4175301638
1.7328129513

0.6
0.3366233389
0.6685450323
0.9914459040
1.3016971244
1.5965446689

0.8
0.2985244887
0.5937600941
0.8826302804
1.1624500568
1.4310481369

1.0
0.2588597329
0.5155692281
0.7680835230
1.0145556439
1.2534069046

(a)

(b)

Figure 1: (a) The exact and approximate solutions for are shown. (b) The error norm for is shown.

Experiment 2. In this experiment [3, 8], we consider (1.1) with and the initial conditions with the boundary conditions (1.3), which can be obtained from the following exact solution: In Tables 3 and 4, we show the and error norms, for several final times (). As we see from Tables 3 and 4, for the given values of , the error norm is less than the error norm. Figure 2 shows four type plots at of Experiment 2.

Table 3: The error norm for different values of of Experiment 2.

T

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

Table 4: The error norm for different values of of Experiment 2.

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

(a)

(b)

(c)

(d)

Figure 2: (a) The exact and approximate solutions for are shown. (b) The error norm for is shown. (c) The error norm for different values of is shown. (d) The error norm for different values of is shown.

Experiment 3. In experiment [3, 9], we consider (1.1) with and the initial conditions with the boundary conditions (1.3), which can be obtained from the following exact solution: In Table 5, we represent the and error norms for different values of at . With increasing , we obtain better results, and these error norms decrease five orders in magnitude. Figure 3 shows four-type plots at of Experiment 3.

Table 5: Error norms for different values of of Experiment 3.

error norm
error norm

4

6

8

10

12

14

16

(a)

(b)

(c)

(d)

Figure 3: (a) The exact and approximate solutions for are shown. (b) The error norm for is shown. (c) The error norm for different values of is shown. (d) The error norm for different values of is shown.

5. Conclusion

In this paper, we have presented the spectral method for solving one-dimensional nonlinear sine-Gordon equation using a new orthogonal polynomial. Numerical experiments show that the spectral method is an efficient one, and the results for value of is more accurate than the other values of . The proposed method has higher order of accuracy even for small values of (i.e., ). We obtain better results when we choose small value of , in comparison to those obtained in the literature.

References

A. D. Polyanin and V. F. Zaitsev, Handbook of Nonlinear Partial Differential Equations, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2004.

R. Rajaraman, Solitons and Instantons, North-Holland, Amsterdam, The Netherlands, 1982.

A. Mohebbi and M. Dehghan, “High-order solution of one-dimensional sine-Gordon equation using compact finite difference and DIRKN methods,” Mathematical and Computer Modelling, vol. 51, no. 5-6, pp. 537–549, 2010.

A. Scott, Nonlinear Science: Emergence and Dynamics of Coherent Structure, vol. 8 of Oxford Texts in Applied and Engineering Mathematics, Oxford University Press, Oxford, UK, 2nd edition, 2003.

T. Dauxois and M. Peyrard, Physics of Solitons, Cambridge University Press, 2006.

M. Dehghan and A. Shokri, “A numerical method for one-dimensional nonlinear sine-Gordon equation using collocation and radial basis functions,” Numerical Methods for Partial Differential Equations, vol. 24, no. 2, pp. 687–698, 2008.

M. Dehghan and D. Mirzaei, “The boundary integral equation approach for numerical solution of the one-dimensional sine-Gordon equation,” Numerical Methods for Partial Differential Equations, vol. 24, no. 6, pp. 1405–1415, 2008.

A. G. Bratsos, “A third order numerical scheme for the two-dimensional sine-Gordon equation,” Mathematics and Computers in Simulation, vol. 76, no. 4, pp. 271–282, 2007.

A. G. Bratsos, “A numerical method for the one-dimensional sine-Gordon equation,” Numerical Methods for Partial Differential Equations, vol. 24, no. 3, pp. 833–844, 2008.

M. Dehghan and A. Shokri, “A numerical method for solution of the two-dimensional sine-Gordon equation using the radial basis functions,” Mathematics and Computers in Simulation, vol. 79, no. 3, pp. 700–715, 2008.

M. Lakestani and M. Dehghan, “Collocation and finite difference-collocation methods for the solution of nonlinear Klein-Gordon equation,” Computer Physics Communications, vol. 181, no. 8, pp. 1392–1401, 2010.

M. Dehghan and F. Fakhar-Izadi, “The spectral collocation method with three different bases for solving a nonlinear partial differential equation arising in modeling of nonlinear waves,” Mathematical and Computer Modelling, vol. 53, no. 9-10, pp. 1865–1877, 2011.

A. H. A. Ali, “Chebyshev collocation spectral method for solving the RLW equation,” International Journal of Nonlinear Science, vol. 7, no. 2, pp. 131–142, 2009.

V. S. Chelyshkov, “Alternative orthogonal polynomials and quadratures,” Electronic Transactions on Numerical Analysis, vol. 25, p. 17–26 (electronic), 2006.

A. Imani, A. Aminataei, and A. Imani, “Collocation method via Jacobi polynomials for solving nonlinear ordinary differential equations,” International Journal of Mathematics and Mathematical Sciences, vol. 2011, Article ID 673085, 11 pages, 2011.

J. M. Franco and I. Gómez, “Accuracy and linear stability of RKN methods for solving second-order stiff problems,” Applied Numerical Mathematics, vol. 59, no. 5, pp. 959–975, 2009.