Mathematical Problems in Engineering

Volume 2014, Article ID 108560, 10 pages

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

## Application of Hybrid Cubic B-Spline Collocation Approach for Solving a Generalized Nonlinear Klien-Gordon Equation

^{1}School of Mathematical Sciences, Universiti Sains Malaysia, 11800 Pulau Pinang, Malaysia^{2}Institute of Engineering Mathematics, Universiti Malaysia Perlis, 02600 Pauh, Perlis, Malaysia^{3}Department of Mathematics, University of Sargodha, Sargodha 40100, Pakistan

Received 24 June 2014; Accepted 11 December 2014; Published 30 December 2014

Academic Editor: Anuar Ishak

Copyright © 2014 Shazalina Mat Zin 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

The generalized nonlinear Klien-Gordon equation is important in quantum mechanics and related fields. In this paper, a semi-implicit approach based on hybrid cubic B-spline is presented for the approximate solution of the nonlinear Klien-Gordon equation. The usual finite difference approach is used to discretize the time derivative while hybrid cubic B-spline is applied as an interpolating function in the space dimension. The results of applications to several test problems indicate good agreement with known solutions.

#### 1. Introduction

Consider a generalized nonlinear Klien-Gordon (KG) equation in the form of [1] subject to the initial conditions and the Dirichlet boundary conditions where denotes the wave displacement at position and time , is a nonlinear function in , and are constants, and, , , , and are known functions.

In the last decade or so, spline functions have been utilized to solve differential equations. For example, Caglar et al. [2] have introduced a cubic B-spline interpolation method to solve two-point boundary value problems. The results obtained were compared to finite difference, finite element, and finite volume method. Caglar et al. [2] concluded that the B-spline interpolation is a better method to interpolate any smooth functions than others. Hamid et al. [3] have developed an alternative cubic trigonometric B-spline interpolation method for the same problem. They have found that the trigonometric B-spline gives better approximation compared to technique used by Caglar et al. Goh et al. [4] have presented a comparison between cubic B-spline and extended cubic B-spline collocation method for solving heat equation. It was concluded that the extended cubic B-spline gives better results. By using the same method, Abbas et al. [5] have solved coupled reaction diffusion system. They found that the B-spline function approximates the system very well and the results are in good agreement with known solutions.

Great deals of research on solving KG equation have been carried out and the results can be found in [6–13]. Dehghan and Shokri [14] have approximated the numerical solution of the nonlinear KG equation using thin plate splines (TPS) radial basis functions. The implementation of the method has been claimed to be simple as the finite difference method and the numerical results obtained were more accurate than others in literature. Khuri and Sayfy [1] have solved the generalized nonlinear KG equation using a finite element collocation approach based on third degree B-spline polynomials. Six examples of nonlinear KG equation including Sine-Gordon equation have been analyzed. The proposed method gives compatible results and better approximation compared to Dehghan and Shokri’s method in [14]. Rashidinia et al. [15] have presented a cubic B-spline collocation method for solving linear KG equation. The results show the proposed scheme is effective and accurate.

In this paper, a new approach by combining hybrid cubic B-spline function and central finite difference is proposed to solve the KG equation. The finite difference approach is used for the time derivative and the hybrid cubic collocation method is applied to interpolate the solutions at space dimension. The scheme obtained is analyzed by Von Neumann stability analysis. To show the feasibility and accuracy of the method, three problems are considered. Numerical solutions, absolute errors, maximum error, and order of convergence are calculated.

#### 2. Temporal Discretization

Consider a uniform mesh with grid points to discretize the grid region with and , where and . The values of and denote mesh space size and time step size, respectively. The Klien-Gordon equation is approximated at time level as follows [16]: where , and the subscripts and are successive time levels. Central difference approach has been used to discretize the time derivative. In order to produce Crank-Nicolson scheme, is chosen to be 0.5. Hence, the scheme becomes At , there is term, , that is outside of domain. Therefore, initial condition (2b) is approximated by the following central difference approach. Thus, The whole scheme is solved numerically by substituting hybrid cubic B-spline function discussed in next section for at each time level .

#### 3. Hybrid Cubic B-Spline Collocation Method

Through this section, hybrid cubic B-spline (HCuBS) is used to solve nonlinear KG equation. The approximate solution, , to the analytical solution, , is considered as where are time-dependent unknowns to be determined and is hybrid cubic B-spline basis function of order 4 as where is cubic B-spline basis function given as and is cubic trigonometric B-spline basis function [18, 19] given as with , , and . The value of plays an important role in the hybrid cubic basis function. If , the basis function is equal to cubic trigonometric B-spline basis function and if , the basis function is equal to cubic B-spline basis function. Hence, this work just considers the value of .

Due to local support properties of B-spline basis function, there are only three nonzero basis functions; namely, , , and are included over subinterval . Thus, the approximate solution and its derivatives with respect to at are where with where These approximations are substituted into (4) to produce the matrix system of order with unknown. In order to generate a unique solution, two additional equations are needed in the system. Hence, boundary conditions (2c) and (2d) are approximated as follows: Thus, the resulting system can be written as where with , , and .

This tridiagonal matrix system can be solved using Thomas Algorithm repeatedly for .

#### 4. Initial State,

The initial vector, , is obtained from initial condition (2a) and boundary values of the derivatives of the initial condition as follows [2, 16]:(i) for ,(ii) for ,(iii) for .

This operation yields matrix system: The solution of the tridiagonal system is obtained by using the Thomas Algorithm [20].

#### 5. Von Neumann Stability Analysis

The growth of error in single Fourier mode is considered as where and is the mode number. It is known that this method is applicable to linear scheme. Hence, (1) is linearized by assuming all nonlinear terms equal zero [15]. The following equation is obtained after substituting (10) into the linear scheme: where By substituting (18) into (19), the following characteristic equation is generated: where , and . Based on Routh-Hurwitz criterion, the transformation, , is applied to the characteristic equation [15, 21]. Then, the equation becomes The necessary and sufficient conditions for are , and . Thus, the following terms have been proved: Hence, this scheme is concluded to be unconditionally stable.

#### 6. Numerical Results and Discussions

In this section, two problems involving KG equation with initial conditions and Dirichlet boundary conditions are tested. In order to measure the accuracy of the method, absolute errors and maximum error are calculated using [22] where and are analytical solution and approximate solution, respectively. The numerical order of convergence, , is obtained by using [1] where and are the at number of partitions and , respectively.

*Problem 1. *Consider the following nonlinear Klien-Gordon equation as [1, 14]
subject to the initial conditions
and boundary conditions
The analytical solution is given by . Figure 1 shows the space-time plot for this analytical solution.