Research Article  Open Access
Ozlem Ersoy, Idris Dag, "The Exponential Cubic BSpline Algorithm for Kortewegde Vries Equation", Advances in Numerical Analysis, vol. 2015, Article ID 367056, 8 pages, 2015. https://doi.org/10.1155/2015/367056
The Exponential Cubic BSpline Algorithm for Kortewegde Vries Equation
Abstract
The exponential cubic Bspline algorithm is presented to find the numerical solutions of the Kortewegde Vries (KdV) equation. The problem is reduced to a system of algebraic equations, which is solved by using a variant of Thomas algorithm. Numerical experiments are carried out to demonstrate the efficiency of the suggested algorithm.
1. Introduction
The splines consist of piecewise functions defined on the distributed knots on problem domain and have certain continuity inside problem subdomain and at the knots. Until now, some types of splines have been developed and especially polynomial splines. The exponential splines are defined as more general splines by McCartin [1–3]. The basis of the exponential splines known as the exponential Bsplines is also given in the studies of McCartin. Existence of the free parameter in the exponential Bsplines yields the different shapes of the splines functions. He has also showed a reliable algorithm by using the exponential spline functions to solve the hyperbolic conservation laws, McCartin and Jameson [4]. However, McCartin stated that application of the exponential spline/exponential Bspline functions has been neglected in the numerical analysis. So, use of the exponential spline in the numerical methods for finding solutions of the differential equations is not common and few papers exist in the literature. McCartin has shown that the exponential splines admit a basis known as the exponential Bsplines. These Bsplines have been started using to form approximates functions recently which are adapted to set up the numerical methods to find solutions of the differential equations recently. An application of the simple exponential splines is considered for setting up the collocation method to solve the numerical solution of singular perturbation problem [5]. Cardinal exponential Bsplines are applied in solving singularly perturbed boundary problems [6]. A variant of Bspline exponential collocation method was also built up for computing numerical solutions of the singularly perturbed boundary value problem [7]. Very recently, the exponential Bspline collocation method has been applied to obtain the solutions of the onedimensional linear convectiondiffusion equation [8].
Types of spline functions are utilized to form approximate solutions for Kortewegde Vries equation (KdVE). The standard Galerkin formulation using the smooth splines on uniform mesh is set up for 1periodic solutions of KdVE by Baker and his coauthors [9]. The Galerkin finite element method together with the cubic Bsplines is used to solve the KdVE in the paper [10]. The quadratic Bspline Galerkin method is built up to find solutions of the KdVE [11]. A collocation solution of the Kortewegde Vries equation using septic Bsplines is proposed by Soliman [12]. A variant of the Galerkin finite element method is designed for solving the KdVE by Aksan and Özdeş [13]. The collocation method using quintic Bsplines is developed to solve the KdVE [14]. A numerical method is developed for the KdVE by using splitting finite difference technique and quintic Bspline functions [15]. The spline finite element method using quadratic polynomial spline for the numerical solution of the KdVE is given by G. Micula and M. Micula [16]. A cubic Bspline TaylorGalerkin method is developed to find numerical solution of the KdVE by Canıvar et al. in [17]. A study based on cubic Bspline finite element method for the solution of the KdVE is suggested by Kapoor et al. [18]. A BubnovGalerkin finite element method with quintic Bspline functions taken as element shape and weight functions is presented for the solution of the KdVE [19]. The paper deals with the numerical solution of the KdVE using quartic Bsplines Galerkin method as both shape and weight functions over the finite intervals [20]. A blended spline quasiinterpolation scheme is employed to solve the onedimensional nonlinear KdVE [21]. A multilevel quartic spline quasiinterpolation scheme is fulfilled to exhibit a large number of physical phenomena for KdVE [22].
The aim of the present paper is to develop an approximate solution of KdVE by collocation method. In Section 2, the exponential Bspline collocation algorithm is defined for the KdVE. In Section 3, the three numerical experiments are constructed to demonstrate the efficiency of the proposed method and the results are documented in tables and graphs are depicted.
We will solve the KdVE: where , are positive parameters and the subscripts and denote differentiation. The boundary conditions will be chosen as
KdVE is prototypical example of exactly solvable mathematical model of waves on shallow water surface. It arises for evolution, interaction of waves, and generation in physics. Due to the term , (1) is called the evolution equation, the nonlinear term causes the steepness of the wave, and the dispersive term defines the spreading of the wave. It is known that the effect of the steepness and spreading results in soliton solutions for the KdVE.
2. Exponential BSpline Collocation Method
The region is partitioned into equal subintervals by points , . On these points together with additional points , outside the domain, the exponential Bsplines, , can be defined as where
forms a basis for the exponential spline space on the interval . On the four consecutive subintervals, an exponential Bspline is defined and it is secondorder continuously differentiable functions.
In Table 1, the values of , , and at the points , which can be obtained from (3) are listed where denotes differentiation with respect to space variable .

The global approximation to the solution will be searched in terms of the unknown parameters and exponential Bspline function defined on the problem domain: Substitution of the points in (5), in its first and its second derivatives, respectively, yields the numerical solution in terms of parameters: where , , and .
Over the subregion , the local approximation is given by where , , and act as subregion parameters and , , and are known as the subregion shape parameters.
To be able to apply the collocation method formed with the exponential Bsplines, KdV equation is spacesplitted as This system includes the secondorder derivatives so that smooth approximation can be done with the exponential Bsplines. To integrate system (8) in time, discretize by the usual finite difference scheme and , , and by CrankNicolson method and we get where represent the solution at the th time level. Here, and is the time step; superscripts denote th time level, .
One linearizes terms and in (9) as to obtain the timeintegrated KdVE: We approximate and in terms of the element parameters and exponential Bsplines separately as
Putting the approximate solution (12) and its derivatives into (11) and evaluating the resulting equations at the points , yield the following system of equations: where
The system consists of linear equation in unknown parameters . A unique solution of the system can be obtained by imposing the boundary conditions , , , to have the following the equations:
Elimination of the parameters , , , , using (15) from system (13), gives a solvable system of linear equation including unknown parameters. Placing solution parameters in (12) when computed from the system via a variant of the Thomas algorithm gives the approximate solution over the subregion . We need the initial parameter vectors , to start the iteration process for system of (13). To do that, the following requirements help to determine the initial parameters:
3. Numerical Tests
Since the conservation laws remain constant at all time, first three numerical conservations are calculated using the rectangular rule for integrals: The error norm is calculated to show the error between analytical and numerical solutions.
(a) The soliton solution of the KdVE is where and . This solution represents propagation of single soliton, having velocity and amplitude .
The analytical solution (19) is used as the initial condition when . The Dirichlet boundary conditions and are adapted to the system to control numerical solutions at the boundaries. Parameters , , , , space step , and time step on the interval from time to are chosen. At time , numerical magnitude of the single soliton is calculated as , so that the numerical amplitude is obtained to be almost the same as the analytical amplitude. Figure 1 illustrates the amplitudes at some times. The distribution of the absolute values of errors can be observed in Figure 2.
error norms and invariants are presented at the selected times in Table 2, as seen from the table that error norm is found small enough and conservation invariants are excellent throughout the simulation. The method gives good results when the free parameter is used. Invariants and remain constant during the run and remains the constant up to the third decimal digits seen in Table 2.
(b) Wave generation is performed by using the Maxwellian initial condition: and boundary conditions:
and and are taken. We have verified the case in which is some critical parameter and, according to the parameters , initial condition breaks up into a number of solitons and, for values , soliton turns into exhibiting the rapidly oscillating wave packets. When together with parameters , , , and , the solution takes the form of the leading soliton and an oscillating tail. This case is shown in Figure 3.
For , we observe a solitary wave plus an oscillating tail (Figure 3). The actual velocity of the wave has been measured and also computed from the measured amplitude using the formula . We find that and , so the solitary waves are indeed solitons. In Table 3 invariant for Maxwellian and .

When , we find three solitons. We have measured the velocity of the largest solitary wave as and calculated the expected velocity from the observed amplitude 1.54468 as . In Figure 4, Maxwellian initial condition is depicted for , , , , and . The invariants are given in Tables 4 and 5 for , and , , respectively.


For , we observed nine solitons moving to the right in Figure 5. The measured velocity of leading soliton is and the corresponding velocities calculated from their measured amplitudes are . The agreement is good. The initial perturbation breaks up into a number of solitons in the course of time depending on the value of chosen. So, if we decrease the value of , then the number of solitons, amplitude, and the velocity increase.
(c) As a final test example, initial condition together with boundary conditions cause the production of a train of solitons depending on the value of for the KdVE. Computation is done on region [] up to time with parameters , , , and . Visual representation of the solution in Figures 6(a)–6(f) is drawn that 10 solitons have been broken up from the given initial condition.
(a)
(b)
(c)
(d)
(e)
(f)
The first three conservation laws are recorded at some times in Table 6. These are favorably constant. The observed velocity of the leading soliton having the amplitude 1.96342 is which was in close agreement with that calculated from its observed amplitude of .

4. Conclusion
The numerical solution of the KdVE is obtained by the collocation method using the exponential basis functions. Performance of the present method is shown by calculating —the error norm, and conservation laws. The present method gives accurate results and simulations such as the propagation of soliton and generation of waves which are substantiated fairly. Using the exponential cubic Bsplines, alternative numerical methods can be set up for finding numerical solutions of the differential equations with high accuracy when an appropriate free parameter is chosen.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 B. J. McCartin, “Theory computation and application of exponential splines,” Tech. Rep. DOE/ER/03077171, 1981. View at: Google Scholar
 B. J. McCartin, “Computation of exponential splines,” SIAM Journal on Scientific and Statistical Computing, vol. 11, no. 2, pp. 242–262, 1990. View at: Publisher Site  Google Scholar  MathSciNet
 B. J. McCartin, “Theory of exponential splines,” Journal of Approximation Theory, vol. 66, no. 1, pp. 1–23, 1991. View at: Publisher Site  Google Scholar  MathSciNet
 B. J. McCartin and A. Jameson, “Numerical solution of nonlinear hyperbolic conservation laws using exponential splines,” Computational Mechanics, vol. 6, no. 2, pp. 77–91, 1990. View at: Publisher Site  Google Scholar
 M. Sakai and R. A. Usmani, “A class of simple exponential $B$splines and their application to numerical solution to singular perturbation problems,” Numerische Mathematik, vol. 55, no. 5, pp. 493–500, 1989. View at: Publisher Site  Google Scholar  MathSciNet
 D. Radunovic, “Multiresolution exponential Bsplines and singularly perturbed boundary problem,” Numerical Algorithms, vol. 47, no. 2, pp. 191–210, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 S. C. Rao and M. Kumar, “Exponential Bspline collocation method for selfadjoint singularly perturbed boundary value problems,” Applied Numerical Mathematics, vol. 58, no. 10, pp. 1572–1581, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 R. Mohammadi, “Exponential Bspline solution of convectiondiffusion equations,” Applied Mathematics, vol. 4, no. 6, pp. 933–944, 2013. View at: Publisher Site  Google Scholar
 G. A. Baker, V. A. Dougalis, and O. A. Karakashian, “Convergence of Galerkin approximations for the Kortewegde Vries equation,” Mathematics of Computation, vol. 40, no. 162, pp. 419–433, 1983. View at: Publisher Site  Google Scholar  MathSciNet
 G. A. Gardner and L. R. T. Gardner, “A finite element solution for the Korteweg de vries equation using cubic Bspline shape functions,” in Proceedings of the International Conference on Modelling and Simulation, vol. 1, AMSE Conferance Press, TassinlaDemiLune, France, 1988. View at: Google Scholar
 L. R. Gardner, G. A. Gardner, and A. H. Ali, “Simulations of solitons using quadratic spline finite elements,” Computer Methods in Applied Mechanics and Engineering, vol. 92, no. 2, pp. 231–243, 1991. View at: Publisher Site  Google Scholar  MathSciNet
 A. A. Soliman, “Collocation solution of the Kortewegde Vries equation using septic splines,” International Journal of Computer Mathematics, vol. 81, no. 3, pp. 325–331, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 E. N. Aksan and A. Özdeş, “Numerical solution of Kortewegde Vries equation by Galerkin Bspline finite element method,” Applied Mathematics and Computation, vol. 175, no. 2, pp. 1256–1265, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 S. I. Zaki, “A quintic Bspline finite element scheme for the KdVB equation,” Computational and Applied Mathematics, vol. 190, pp. 532–547, 2006. View at: Google Scholar
 P. C. Jain, R. Shankar, and D. Bhardwaj, “Numerical solution of the Kortewegde Vries (KDV) equation,” Chaos, Solitons and Fractals, vol. 8, no. 6, pp. 943–951, 1997. View at: Publisher Site  Google Scholar  MathSciNet
 G. Micula and M. Micula, “On the numerical approach of Kortewegde VriesBurger equations by spline finite element and collocation methods,” Seminar on Fixed Point Theory ClujNapoca, vol. 3, pp. 261–270, 2002. View at: Google Scholar  MathSciNet
 A. Canıvar, M. Sari, and I. Dag, “A TaylorGalerkin finite element method for the KdV equation using cubic Bsplines,” Physica B: Condensed Matter, vol. 405, no. 16, pp. 3376–3383, 2010. View at: Publisher Site  Google Scholar
 S. Kapoor, S. Rawat, and S. Dhawan, “Numerical investigation of separated solitary waves solution for KdV equation through finite element technique,” International Journal of Computer Applications, vol. 40, no. 14, pp. 27–33, 2012. View at: Publisher Site  Google Scholar
 N. K. Amein and M. A. Ramadan, “A small time solutions for the KdV equation using BubnovGalerkin finite element method,” Journal of the Egyptian Mathematical Society, vol. 19, no. 3, pp. 118–125, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 B. Saka and I. Dag, “Quartic Bspline GALerkin approach to the numerical solution of the KdVB equation,” Applied Mathematics and Computation, vol. 215, no. 2, pp. 746–758, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 R. Yu, R. Wang, and C. Zhu, “A numerical method for solving KdV equation with blended bspline quasiinterpolation,” Journal of Information and Computational Science, vol. 10, no. 16, pp. 5093–5101, 2013. View at: Publisher Site  Google Scholar
 R.G. Yu, R.H. Wang, and C.G. Zhu, “A numerical method for solving KdV equation with multilevel Bspline quasiinterpolation,” Applicable Analysis, vol. 92, no. 8, pp. 1682–1690, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 S. I. Zaki, “A quintic Bspline finite elements scheme for the KdVB equation,” Computer Methods in Applied Mechanics and Engineering, vol. 188, no. 1, pp. 121–134, 2000. View at: Publisher Site  Google Scholar
 B. Saka, “Cosine expansionbased differential quadrature method for numerical solution of the KdV equation,” Chaos, Solitons and Fractals, vol. 40, no. 5, pp. 2181–2190, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 Dağ and Y. Dereli, “Numerical solutions of KdV equation using radial basis functions,” Applied Mathematical Modelling, vol. 32, no. 4, pp. 535–546, 2008. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2015 Ozlem Ersoy and Idris Dag. 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.