Homotopy Analysis Method for Conformable Burgers-Korteweg-de Vries Equation

The main goal of this paper is finding the approximate analytical solution of BurgersKorteweg-de Vries with newly defined conformable derivative by using homotopy analysis method (HAM). Then the approximate analytical solution is compared with the exact solution and comparative tables are given. Introduction The adventure of fractional calculus started with the L’Hospital’s letter to Leibnitz in 1695. In this letter the original question which caused the name fractional question was: Can a derivative of integer order dy dxn be extended to have the meaning when n is a fractional number? This question has a golden value for the development of fractional calculus. After that time fractional calculus has been a valuable tool for expressing the nonlinear phenomenons in the nature. Thus scientists presented various applications of fractional derivatives and integrals for the physical interpretation the real world circumstances. For this interpretation most of scientists use popular and known derivative formulas such as Riemann-Liouville, Caputo [1, 2, 3]. As time progressed researchers found some insufficiencies at these definitions. It is understood that these definitions do not satisfy basic properties of the Newtonian concept, for instance the product, quotient, the Chain rules and etc. Researchers who study in the field of fractional calculus has been worried about these problems. To overcome this issue, R. Khalil et al. [4] introduced a new derivative and its anti-derivative called ’conformable derivative and integral’. Definition. Let f : [0,∞) → R be a function. The α order ”conformable derivative” of f is defined by, Tα(f)(t) = lim ε→0 f(t+ εt1−α)− f(t) ε , for all t > 0, α ∈ (0, 1). If f is α-differentiable in some (0, a), a > 0 and lim t→0+ f (t) exists then define f (0) = lim t→0+ f (t) and the ”conformable integral” of a function f starting from a ≥ 0 is defined as: I α(f)(t) = t ˆ a f(x) x1−α dx where the integral is the usual Riemann improper integral, and α ∈ (0, 1]. To test the efficiency and the accuracy of the conformable derivative, researchers made huge number of scientific articles on it. For example T. Abdeljawad [5] has presented conformable versions of the chain rule, exponential functions, Gronwalls inequality, integration by parts, Taylor power series expansions and Laplace transform. Iyiola et al. [6] expressed the analytical solution of space-time fractional Fornberg-Whitham equation, Hesameddini et al. [7] demonstrated the numerical solution Bulletin of Mathematical Sciences and Applications Submitted: 2016-08-11 ISSN: 2278-9634, Vol. 17, pp 17-23 Revised: 2016-10-21 doi:10.18052/www.scipress.com/BMSA.17.17 Accepted: 2016-10-24 2016 SciPress Ltd, Switzerland Online: 2016-11-01 SciPress applies the CC-BY 4.0 license to works we publish: https://creativecommons.org/licenses/by/4.0/ of multi-order fractional differential equations via the sinc collocation method and K. R. Prasad et al. [8] discussed prove the existence of multiple positive solutions for a coupled system of iterative type boundary value problems involving new conformable fractional order derivative. Therewithal to these studies O. Acan et al. [9] introduced conformable fractional reduced differential transform method using Conformable Calculus and reduced differential transform method. A. Kurt et al. [10] used conformable derivative while obtaining the analytical and approximate solution of Burgers’ equation. Finally A. Atangana et al. [11] introduced the new properties of conformable derivative. Hence, there are many physical applications can be considered in this new area. In the following theorem some properties of this new definition is given [4]. Theorem 1 Let α ∈ (0, 1] and f, g α-differentiable at point t > 0. Then 1. Tα(cf + dg) = cTα(f) + cTα(g) for all a, b ∈ R. 2. Tα(t) = ptp−α for all p ∈ R. 3. Tα(λ) = 0 for all constant functions f(t) = λ. 4. Tα(fg) = fTα(g) + gTα(f). 5. Tα ( f g ) = gTα(f)−fTα(g) g2 . 6. If in addition to f is differentiable, then Tα(f)(t) = t1−α df dt . Homotopy Analysis Method In this section, the authors applied HAM to the processed problem. To implement the fundamentals of the method we assume the following differential equation, N [u(x, t)] = 0 whereN is a nonlinear operator, x and t show independent variables and u(x, t) is an unknown function. By using the generalization of HAM, a zero-order deformation equation has been builded by Liao [12, 13] (1− p)L [φ(x, t; p)− u0(x, t)] = p~N [φ(x, t; p)] (1) where p ∈ [0, 1] is the embedding parameter, ~ ̸= 0 is an auxiliary parameter, L is an auxiliary linear operator, u0(x, t) is an initial guess of u(x, t), φ(x, t; p) is an unknown function, respectively. In this way, it would be thinkable to choose auxiliary parameters and operators in HAM. When p is selected as p = 0 and p = 1 then φ(x, t; 0) = u0(x, t), φ(x, t; 1) = u(x, t) are gained successively. Hence, as long as the embedding parameter p increases from 0 to 1, the solutions φ(x, t; p) show a change from the initial value u0(x, t) to the solution u(x, t). If φ(x, t; p) is expanded in Taylor series with respect to the embedding parameter p, we have: φ(x, t; p) = u0(x, t) + ∞ ∑