NUMERICAL SOLUTIONS OF TIME FRACTIONAL KORTEWEG—DE VRIES EQUATION AND ITS STABILITY ANALYSIS

In this study, the fractional derivative and finite difference operators are analyzed. The time fractional KdV equation with initial condition is considered. Discretized equation is obtained with the help of finite difference operators and used Caputo formula. The inherent truncation errors in the method are defined and analyzed. Stability analysis is explored to demonstrate the accuracy of the method. While doing this analysis, considering conservation law, with the help of using the definition discovered by Lax-Wendroff, von Neumann stability analysis is applied. The numerical solutions of time fractional KdV equation are obtained by using finite difference method. The comparison between obtained numerical solutions and exact solution from existing literature is made. This comparison is highlighted with the graphs as well. Results are presented in tables using the Mathematica software package wherever it is needed.


Introduction
Nowadays, one of the developing conceptions is the fractional di¤erential equations. This notion began to develop since 17 th century with the help of several mathematicians' studies on di¤erential and integration, like Leibniz, Euler, Lagrange, Abel, Liouville etc. [1,2,3]. (0: 5) th order derivative was de…ned by Leibniz in the year 1695. Riemann-Liouville, Hadamard, Grunwald-Letnikov, Riesz and Caputo have given the integral inequalities to the literature. In 2006, by Kilbas, Srivastava, Trujillo and in 1993 by Samko, Kilbas, Marichev de…ned the fractional theory and di¤erent derivatives with developments [4,5].
The exact solutions of the fractional di¤erential equations may not be easily obtained, so we need numerical methods for fractional di¤erential equations. One of them is …nite di¤erent method and it is one of the most popular methods of 354 ASIF YOKUŞ numerical solution of partial di¤erential equations. There are some studies about this method's stability analysis. B.F. Feng, in his study, examined Von Neumann's Stability analysis by linearizing Korteweg-de Vries (in short, KdV) equation.
In this study, classical partial di¤erential equations have been extended to the fractional partial di¤erential equations. There are many applications of this equation in the literature. The fractional partial di¤erential equations have been used in applications such as ‡uid, ‡ow, …nance, hydrology and others [6,21]. In this paper, we investigate …nite di¤erence numerical methods to solve the time fractional KdV equation of the form [22] @ u(x; t) @t where m is an integer that m 1 < m. The function (:) is called as Gamma function.

Analysis of Finite Difference Method
Let us de…ne some notations to describe the …nite forward di¤erence method.
x is the spatial step, t is the time step, x i = a+i x, i = 0; 1; 2; : : : ; N points are the coordinates of mesh and N = b a x , t j = j t, j = 0; 1; 2; : : : ; M and M = T t . The function u(x; t) is the value of the solution at these grid points which are u(x i ; t j ) = u i;j , where we denote by u i;j the numerical estimate of the exact value of u(x; t) at the point (x i ; t j ). Now, we de…ne the di¤erence operators as H t u i;j = u i;j+1 u i;j ; (2.1) Thus, partial derivatives are approximated through the …nite di¤erence operators as According to the shifted Caputo de…nition [23], There are many studies in the literature on fractional derivatives of Taylor Series. The generalized Taylor series which is in these studies has been awarded by Odibat [24]. @ u(x; t) @t In the …nite di¤erence method, substituting Eqs. (2.5), (2.6) and (2.7) into Eq.
(1.1) can be written as indexed

Consistency Analysis and Truncation Error
In this section, we investigate the consistency the Eq. (1.1) by the …nite di¤erence method. At …rst, Taylor series expansions can be given in the form as follows, Now, let us de…ne an operator L, The indexed form of operator L can be written as If we substitute the indexed form Eqs.
Eq. (1.1) can be written as indexed The truncation error is O( t 2 + x 2 ).

Linear Stability Analysis
In this section, we mainly study the stability for the …nite di¤erence method. To describe this method, we consider the …rst-order conservation equation @u @t + @u @x = 0; where p = i x, q = j t and I = p 1. If we substitute the Eq. (4.6) into the Eq. (4.5) yields: (4.7) According to the Von Neumann's Stability analysis; if j j 1, …nite di¤erence method for the KdV equation is stable.
For the Eq. (4.8) the stability depends on the constant . However, due to the nature of the method of …nite di¤erence, stability will be examined with respect to parameter h. For this reason, if we choose = 2 , U = 2 and = 1 in the Eq. (4.8) and have 1 x 1, then the …nite di¤erence method for the KdV equation is stable. By using the Eq. (4.7), neutral stability curve can be drawn [26] for the example Eq. (4.8).
The neutral stability curve is locally a parabola with minimum (0, 0). As shown in the graphs, if we choose the x close to zero, …nite di¤erence methods for the KdV equation is stable. In other words, the …nite di¤erence algorithm is stable if

Numerical Example
We consider the fractional KdV equation of the form Eq. (1.1) with the initial condition as follow: In the following numerical experiments we choose = 0:8. The fractional KdV Eq. (1.1) together with the above initial condition is constructed [22] such that the exact solution is The numerical solutions are obtained from the …nite di¤erence schemes discussed above considering Eq. (2.9). The numerical solutions in the interval 0 x 1: and the numerical solutions in the interval 1 x < 0:  Table 2. Numerical and exact solutions of Eq. (1.1) and absolute errors when x = 0:02 and 1 x < 0 We know that truncation error will be small if x and t choose su¢ ciently small. There are appointed values close to zero indicate that the truncation error

Conclusions
In this study, we considered the numerical solution of fractional dispersion equation by using Finite Di¤erence Method. The method can be applied to many other nonlinear equations. What is more, this method is also computerizable, which allows us to perform complicate and tedious algebraic calculation on a computer. Fractional …nite di¤erence methods are useful to solve the fractional di¤erential equations. In some way, these numerical methods have similar form with the classical equations. Some of them can be seen as the generalizations of the …nite di¤erence methods for the typical di¤erential equations. The numerical method for solving the fractional reaction-dispersion equation has been described and demonstrated. Finally, we point out that, for given equation with initial values, the corresponding analytical and numerical solutions are obtained according to the recurrence Eq. (2.9) using Mathematica software package.